75 G4double beta = std::sqrt(tau*(tau + 2.0))/gamma;
80 G4double a1 = 0.5*beta*gamma*tau*(gamma-2.0);
94 tsam = 2.0*ac * (2.0*rand + a2*std::sqrt(rand)) / (a2*a2 - 4.0*rand);
95 gtr = (2.0 - tsam) * (a1 + 1.0/(ac+tsam));
99 costheta = 1.0 - tsam;
101 G4double sint = std::sqrt(tsam*(2.0 - tsam));
113 G4cout <<
"Non-polarized photoelectric effect angular generator." <<
G4endl;
114 G4cout <<
"The Sauter-Gavrila distribution for the K-shell is used."<<
G4endl;
115 G4cout <<
"Originally developed by M.Maire for Geant3"