77 const G4double signc = (c < 0.) ? -1.0 : 1.0;
82 const G4double beta = std::sqrt(tau*(tau + 2.))/(tau + 1.);
84 return (cosTheta + beta)/(1. + cosTheta*beta);
95 const G4double sinTheta = std::sqrt((1. - cosTheta)*(1. + cosTheta));
126 const G4double sinp = std::sin(phi);
127 const G4double cosp = std::cos(phi);
130 G4double sint = std::sqrt((1. - cost)*(1. + cost));
132 dirElectron.
set(sint*cosp, sint*sinp, cost);
136 sint = std::sqrt((1. - cost)*(1. + cost));
138 dirPositron.
set(-sint*cosp, -sint*sinp, cost);
147 G4cout <<
"Angular Generator based on classical formula from" <<
G4endl;
148 G4cout <<
"J.D. Jackson, Classical Electrodynamics, Wiley, New York 1975"