69 if(kineticEnergy <= 0.0) {
return 0.0; }
78 G4double alog10E = std::log10(kineticEnergy);
80 static const G4double nuleonRadius=1.36e-15;
87 G4double crossSection = 1.0E31*fac*fac2*(1. + 1./a13 - fac1);
90 crossSection *= (1 - 0.15*
G4Exp(-kineticEnergy))/(1.0 - 0.0007*a);
97 G4double ff4= 1.0 - (1.0/(1+
G4Exp(-8*ff1*(alog10E + 1.37*ff2))));
99 crossSection *= (1 + ff3*ff4);
103 ff1=1. - 1./a - 0.001*
a;
104 ff2=1.17 - 2.7/a - 0.0014*
a;
106 ff4=-8.*ff1*(alog10E + 2.0*ff2);