36 alphaZ = fine_structure_const*
Z;
39 Rnuc = 0.5*fine_structure_const*std::pow(
A, 0.33333);
42 V0 = 1.13*fine_structure_const*fine_structure_const
43 *std::pow(std::abs(Z), 1.33333);
68 if (Wprime <= 1.00001) Wprime = 1.00001;
71 G4double p_e = std::sqrt(Wprime*Wprime - 1.);
82 G4double factor3 = (Wprime/W)*std::sqrt( (Wprime*Wprime - 1.)/(W*W - 1.) );
84 return factor1*factor2*factor3;
96 G4double factor1 = std::pow( (1+re)*(1+re) + im*im, re+0.5);
97 G4double factor2 = std::exp(2*im * std::atan(im/(1+re)));
98 G4double factor3 = std::exp(2*(1+re));
100 G4double factor5 = std::exp( (1+re)/( (1+re)*(1+re) + im*im)/6 );
102 return factor1*factor4*factor5/factor2/factor3/factor6;
119 for (
G4int i = 1; i < 6; i++) sum = sum*x +
gc[i];
145 G4double w = std::sqrt(1. + p_e*p_e);
146 factor = 1. + c1*w + c2/w + c3*w*w;
156 G4double term2 = 12.*(2. + gamma1)*p_e*p_e
157 *std::pow(twoPR, 2.*(gamma1-
gamma0-1) )
160 factor = term1 + term2;
171 G4double gamma2 = std::sqrt(9. - alphaZ*alphaZ);
177 G4double term2 = 4.*(2. + gamma1)*e_nu*e_nu*p_e*p_e
178 *std::pow(twoPR, 2.*(gamma1-
gamma0-1.) )
182 G4double term3 = 180.*(3.+gamma2)*p_e*p_e*p_e*p_e
183 *std::pow(twoPR, 2.*(gamma2-
gamma0-2) )
187 factor = term1 + term2 + term3;
198 G4double gamma2 = std::sqrt(9. - alphaZ*alphaZ);
199 G4double gamma3 = std::sqrt(16. - alphaZ*alphaZ);
205 G4double term1 = e_nu*e_nu*e_nu*e_nu*e_nu*e_nu*(1. +
gamma0)/1260.;
207 G4double term2 = 2.*(2. + gamma1)*e_nu*e_nu*e_nu*e_nu*p_e*p_e
208 *std::pow(twoPR, 2.*(gamma1-
gamma0-1.) )
212 G4double term3 = 60.*(3.+gamma2)*p_e*p_e*p_e*p_e*e_nu*e_nu
213 *std::pow(twoPR, 2.*(gamma2-
gamma0-2.) )
217 G4double term4 = 2240.*p_e*p_e*p_e*p_e*p_e*p_e*(4. + gamma3)
218 *std::pow(twoPR, 2.*(gamma3-
gamma0-3.) )
222 factor = term1 + term2 + term3 + term4;
227 G4Exception(
"G4BetaDecayCorrections::ShapeFactor()",
"HAD_RDM_010",
229 "Transition not yet implemented - using allowed shape");
static c2_factory< G4double > c2
static const G4double fac
static const G4double A[nN]
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4double ModSquared(const G4double &x, const G4double &y)
G4double ShapeFactor(const G4BetaDecayType &, const G4double &p_e, const G4double &e_nu)
G4double Gamma(const G4double &arg)
G4BetaDecayCorrections(G4int Z, G4int A)
G4double FermiFunction(const G4double &W)