94 G4double eMomentum = std::sqrt(eKE*(eKE + 2.*eMass) );
97 G4double sinTheta = std::sqrt(1.0 - cosTheta*cosTheta);
112 - eMomentum*eMomentum)/(parentMass - eTE + eMomentum*cosThetaENu)/2.;
114 G4double sinThetaENu = std::sqrt(1.0 - cosThetaENu*cosThetaENu);
120 nuDirection.
setX(sinThetaENu*cosPhiNu*cosTheta*cosPhi -
121 sinThetaENu*sinPhiNu*sinPhi + cosThetaENu*sinTheta*cosPhi);
122 nuDirection.
setY(sinThetaENu*cosPhiNu*cosTheta*sinPhi +
123 sinThetaENu*sinPhiNu*cosPhi + cosThetaENu*sinTheta*sinPhi);
124 nuDirection.
setZ(-sinThetaENu*cosPhiNu*sinTheta + cosThetaENu*cosTheta);
134 -eDirection*eMomentum - nuDirection*nuEnergy);
162 const G4int& daughterA,
178 for (
G4int ptn = 0; ptn < npti; ptn++) {
180 e = 1. + e0*(ptn + 0.5)/
G4double(npti);
181 p = std::sqrt(e*e - 1.);
182 f = p*e*(e0 - e + 1.)*(e0 - e + 1.);
188 f *= corrections.
ShapeFactor(betaType, p, e0-e+1.);