Changeset 1055 for trunk/source/processes/hadronic/cross_sections/src/G4IonsShenCrossSection.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/cross_sections/src/G4IonsShenCrossSection.cc (modified) (3 diffs)

trunk/source/processes/hadronic/cross_sections/src/G4IonsShenCrossSection.cc

@@ -47 +47 @@
 
 // Apply energy check, if less than lower limit then 0 value is returned
-   if (  ke_per_N < lowerLimit )
-      return xsection;
+   // if (  ke_per_N < lowerLimit ) return xsection;
 
    G4int At = G4int(AA);
@@ -63 +62 @@
    G4double r = Rt + Rp + 3.2;   // in fm
    G4double b = 1.0;   // in MeV/fm
+   G4double targ_mass = G4ParticleTable::GetParticleTable()->GetIonTable()->GetIonMass( Zt , At );
+   G4double proj_mass = aParticle->GetMass();
+   G4double proj_momentum = aParticle->GetMomentum().mag();
+
+   G4double Ecm = calEcmValue ( proj_mass , targ_mass , proj_momentum ); 
 
    G4double B = 1.44 * Zt * Zp / r - b * Rt * Rp / ( Rt + Rp ); 
-
+   if(Ecm <= B) return xsection;
    //G4double ke_per_N = aParticle->GetKineticEnergy() / Ap; 
 
@@ -74 +78 @@
    G4double R2 = 1.0 * ( At - 2 * Zt ) * Zp / ( Ap * At );
 
-   G4double targ_mass = G4ParticleTable::GetParticleTable()->GetIonTable()->GetIonMass( Zt , At );
-   G4double proj_mass = aParticle->GetMass();
-   G4double proj_momentum = aParticle->GetMomentum().mag();
-
-   G4double Ecm = calEcmValue ( proj_mass , targ_mass , proj_momentum ); 
 
    G4double R3 = 0.176 / std::pow ( G4double(Ecm) , G4double(one_third) ) * cubicrAt * cubicrAp / ( cubicrAt + cubicrAp );