93 G4double eMomentum = std::sqrt(eKE*(eKE + 2.*eMass) );
98 - eMomentum*eMomentum)/(parentMass - eTE + eMomentum*cosThetaENu)/2.;
102 G4double sinTheta = std::sqrt(1.0 - cosTheta*cosTheta);
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);
164 const G4int& daughterA,
179 for (
G4int ptn = 0; ptn < npti; ptn++) {
182 p = std::sqrt(e*e - 1.);
183 f = p*e*(e0 - e + 1.)*(e0 - e + 1.);
189 f *= corrections.
ShapeFactor(betaType, p, e0-e+1.);