43 ResidualA = ResidualZ = theA = theZ = FragmentA = 0;
44 ResidualAthrd = FragmentAthrd = 0.0;
58 CalcAlphaParam(fragment);
79 if (
OPTxs==0) {std::ostringstream errOs;
80 errOs <<
"We should'n be here (OPT =0) at evaporation cross section calculation (neutrons)!!" <<
G4endl;
83 else if(
OPTxs==1 ||
OPTxs==2)
return GetOpt12( K);
84 else if (
OPTxs==3 ||
OPTxs==4)
return GetOpt34( K);
86 std::ostringstream errOs;
87 errOs <<
"BAD NEUTRON CROSS SECTION OPTION AT EVAPORATION!!" <<
G4endl;
103 if (K > 50*
MeV) { Kc = 50*
MeV; }
105 G4double landa, landa0, landa1, mu, mum0, mu1,nu, nu0, nu1, nu2,xs;
114 landa = landa0/ResidualAthrd + landa1;
115 mu = mum0*ResidualAthrd + mu1*ResidualAthrd*ResidualAthrd;
116 nu = nu0*ResidualAthrd*ResidualA + nu1*ResidualAthrd*ResidualAthrd + nu2 ;
117 xs=landa*Kc + mu + nu/Kc;
119 std::ostringstream errOs;
120 G4cout<<
"WARNING: NEGATIVE OPT=1 neutron cross section "<<
G4endl;
121 errOs <<
"RESIDUAL: Ar=" << ResidualA <<
" Zr=" << ResidualZ <<
G4endl;
122 errOs <<
" xsec("<<Kc<<
" MeV) ="<<xs <<
G4endl;
131 G4double landa, landa0, landa1, mu, mum0, mu1,nu, nu0, nu1, nu2;
133 G4double flow,ec,ecsq,xnulam,etest(0.),ra(0.),
a,signor(1.),sig;
158 if (ResidualA < 40) { signor =0.7 + ResidualA*0.0075; }
159 if (ResidualA > 210) { signor = 1. + (ResidualA-210)/250.; }
160 landa = landa0/ResidualAthrd + landa1;
161 mu = mum0*ResidualAthrd + mu1*ResidualAthrd*ResidualAthrd;
162 nu = nu0*ResidualAthrd*ResidualA + nu1*ResidualAthrd*ResidualAthrd + nu2;
165 if (nu < 0.) { nu=-nu; }
176 a = -2.*p*ec + landa - nu/ecsq;
177 b = p*ecsq + mu + 2.*nu/ec;
180 if (cut > 0.) { ecut = std::sqrt(cut); }
181 ecut = (ecut-
a) / (p+p);
183 if (cut < 0.) { ecut2 = ecut - 2.; }
184 elab = K * FragmentA /
G4double(ResidualA);
187 if (elab > ecut2) { sig = (p*elab*elab+a*elab+
b) * signor; }
190 sig = (landa*elab+mu+nu/elab) * signor;
192 if (xnulam < flow || elab < etest) {
return sig; }
193 geom = std::sqrt(theA*K);
194 geom = 1.23*ResidualAthrd + ra + 4.573/geom;
195 geom = 31.416 * geom * geom;
G4GLOB_DLL std::ostream G4cout
G4double Z13(G4int Z) const
virtual ~G4NeutronEvaporationProbability()
G4NeutronEvaporationProbability()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double Z23(G4int Z) const