57 const G4double tGamma = 87.65*fermi*fermi*fermi*fermi*fermi*
fermi;
63 G4int nBarion = Barions.size();
78 Etot += Vo*0.5/r12*
G4Exp(1/(4*Alpha*GammaY*GammaY))*
79 (
G4Exp(-r12/GammaY)*(1 -
Erf(0.5/GammaY/std::sqrt(Alpha) - std::sqrt(Alpha)*r12)) -
80 G4Exp( r12/GammaY)*(1 -
Erf(0.5/GammaY/std::sqrt(Alpha) + std::sqrt(Alpha)*r12)));
88 for(
G4int c3 =
c2 + 1; c3 < nBarion; c3++)
115 const G4double P2[5] = { 7.3738883, 6.8650185, 3.0317993, 0.56316962, 4.3187787e-5 };
116 const G4double Q2[5] = { 7.3739609, 15.184908, 12.79553, 5.3542168, 1. };
130 H = X*(P10 + P11*
Y)/(Q10+Y);
138 for(c1 = 3; c1 >= 0; c1--)
143 H = 1 -
G4Exp(-V*V)*AP/AQ;
148 H = 1 -
G4Exp(-V*V)*(C1+Y*(P30 + P31*
Y)/(Q30 + Y))/V;
212 for(
size_t cParticle = 0; cParticle < theBarions.size(); cParticle++)