65 ratio1(1),ratio2(1),
delta(0)
93 G4double beta = std::sqrt((gamma - 1)*(gamma + 1))/gamma;
107 y = q*ymax/(1 + ymax*(1 - q));
111 if(gfun > gMax &&
nwarn >= 20) {
113 G4cout <<
"### WARNING in G4Generator2BS: Etot(MeV)= " << energy/
MeV
114 <<
" Egamma(MeV)" << (energy - final_energy)/
MeV
115 <<
" gMax= " << gMax <<
" < " << gfun
116 <<
" results are not reliable!"
119 G4cout <<
" WARNING in G4Generator2BS is closed" <<
G4endl;
127 G4double sint = std::sqrt((1 - cost)*(1 + cost));
139 G4cout <<
"Bremsstrahlung Angular Generator is 2BS Generator "
140 <<
"from 2BS Koch & Motz distribution (Rev Mod Phys 31(4), 920 (1959))" <<
G4endl;
141 G4cout <<
"Sampling algorithm adapted from PIRS-0203" <<
G4endl;