67 rj =
static_cast<G4double>(nParticles - nCharged)/
68 static_cast<G4double>(nParticles);
115 if (K > 50*
MeV) { Kc = 50*
MeV; }
128 G4double xs = landa*Kc + mu + nu/Kc;
160 if (nu < 0.) { nu = -nu; }
172 G4double b = p*ecsq + mu + 2.*nu/ec;
175 if (cut > 0.) { ecut = std::sqrt(cut); }
176 ecut = (ecut-
a) / (2*p);
177 if (cut < 0.) { ecut -= 2.; }
183 if (elab > ecut) { sig =
std::max(0.0,(p*elab*elab+a*elab+b) * signor); }
186 sig = (landa*elab+mu+nu/elab) * signor;
188 if (xnulam >= flow && elab >= etest) {
189 geom = std::sqrt(
theA*K);
190 geom = 1.23*ResidualAthrd + ra + 4.573/geom;
191 geom = 31.416 * geom * geom;
G4double GetOpt12(G4double K)
G4double ResidualA13() const
G4double GetOpt34(G4double K)
G4double GetOpt0(G4double ekin)
virtual G4double CrossSection(G4double ekin)
G4double Z13(G4int Z) const
virtual ~G4PreCompoundNeutron()
virtual G4double GetAlpha()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual G4double GetBeta()
virtual G4double GetRj(G4int NumberParticles, G4int NumberCharged)