Changeset 1315 for trunk/source/processes/hadronic/models/chiral_inv_phase_space/processes/src/G4QInelastic.cc

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/processes/hadronic/models/chiral_inv_phase_space/processes/src/G4QInelastic.cc (modified) (5 diffs)

trunk/source/processes/hadronic/models/chiral_inv_phase_space/processes/src/G4QInelastic.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4QInelastic.cc,v 1.5 2009/12/03 18:09:07 mkossov Exp $
-// GEANT4 tag $Name: geant4-09-03 $
+// $Id: G4QInelastic.cc,v 1.8 2010/06/09 12:56:01 mkossov Exp $
+// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
 //
 //      ---------------- G4QInelastic class -----------------
@@ -68 +68 @@
 
 G4bool   G4QInelastic::manualFlag=false; // If false then standard parameters are used
-G4double G4QInelastic::Temperature=210.; // Critical Temperature (sensitive at High En)
+G4double G4QInelastic::Temperature=200.; // Critical Temperature (sensitive at High En)
 G4double G4QInelastic::SSin2Gluons=0.3;  // Supression of s-quarks (in respect to u&d)
 G4double G4QInelastic::EtaEtaprime=0.3;  // Supression of eta mesons (gg->qq/3g->qq)
@@ -74 +74 @@
 G4double G4QInelastic::freeDib=.05;      // Percentage of free diBaryons on the surface
 G4double G4QInelastic::clustProb=5.;     // Nuclear clusterization parameter
-G4double G4QInelastic::mediRatio=1.;     // medium/vacuum hadronization ratio
+G4double G4QInelastic::mediRatio=10.;    // medium/vacuum hadronization ratio
 G4int    G4QInelastic::nPartCWorld=152;  // The#of particles initialized in CHIPS World
 G4double G4QInelastic::SolidAngle=0.5;   // Part of Solid Angle to capture (@@A-dep.)
@@ -757 +757 @@
   G4double sr=std::sqrt(am);
   G4double dsr=0.01*(sr+sr);
+  //G4double medRA=mediRatio*pow(am,third);
+  G4double medRA=mediRatio;
   if(dsr<dd)dsr=dd;
-  if(manualFlag) G4QNucleus::SetParameters(freeNuc,freeDib,clustProb,mediRatio); // ManualP
+  if(manualFlag) G4QNucleus::SetParameters(freeNuc,freeDib,clustProb,medRA); // ManualP
   else if(projPDG==-2212) G4QNucleus::SetParameters(1.-dsr-dsr,dd+dd,5.,10.);//aP ClustPars
-  else if(projPDG==-211)  G4QNucleus::SetParameters(.67-dsr,.32-dsr,5.,9.);//Pi- ClustPars
+  else if(projPDG==-211)  G4QNucleus::SetParameters(.67-dsr,.32-dsr,5.,9.); //Pi- ClustPars
 #ifdef debug
   G4cout<<"G4QInelastic::PostStepDoIt: N="<<N<<" for element with Z="<<Z<<G4endl;
@@ -2106 +2108 @@
     }
     if(!theDefinition)
-    {
-#ifdef debug
+    { // !! Commenting of this print costs days of searching for E/mom nonconservation !!
+      //#ifdef debug
       G4cout<<"---Warning---G4QInelastic::PostStepDoIt: drop PDG="<<PDGCode<<G4endl;
-#endif
+      //#endif
       delete hadr;
       continue;