Changeset 1055 for trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4He3EvaporationProbability.cc

Timestamp:

May 28, 2009, 4:26:57 PM (15 years ago)

Author:

garnier

Message:

maj sur la beta de geant 4.9.3

File:

• trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4He3EvaporationProbability.cc (modified) (4 diffs)

trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4He3EvaporationProbability.cc

@@ -104 +104 @@
 //OPT=3,4 Kalbach's parameterization 
 // 
- G4double G4He3EvaporationProbability::CrossSection(const  G4Fragment & fragment, const  G4double K)
+G4double G4He3EvaporationProbability::CrossSection(const  G4Fragment & fragment, const  G4double K)
 {
        theA=GetA();
@@ -117 +117 @@
 
        if (OPTxs==0) {std::ostringstream errOs;
-         errOs << "We should'n be here (OPT =0) at evaporation cross section calculation (He3's)!!"  <<G4endl;
+         errOs << "We should'n be here (OPT =0) at evaporation cross section calculation (He3's)!!"  
+               <<G4endl;
          throw G4HadronicException(__FILE__, __LINE__, errOs.str());
          return 0.;}
@@ -200 +201 @@
   G4double      ra=0.80;
   
-  ec = 1.44 * theZ * ResidualZ / (1.5*ResidualAthrd+ra);
+  //JMQ 13/02/09 increase of reduced radius to lower the barrier
+  // ec = 1.44 * theZ * ResidualZ / (1.5*ResidualAthrd+ra);
+  ec = 1.44 * theZ * ResidualZ / (1.7*ResidualAthrd+ra);
   ecsq = ec * ec;
   p = p0 + p1/ec + p2/ecsq;
@@ -221 +224 @@
   elab = K * FragmentA / ResidualA;
   sig = 0.;
-  
+
   if (elab <= ec) { //start for E<Ec
-    if (elab > ecut2)  sig = (p*elab*elab+a*elab+b) * signor;    
+    if (elab > ecut2)  sig = (p*elab*elab+a*elab+b) * signor;
   }           //end for E<Ec
   else {           //start for E>Ec