49 : polarFlag(
false), fDirection(0.,0.,0.), fTwoJMAX(10), fVerbose(0)
82 G4cout <<
"G4GammaTransition::GenerateGamma - Etrans(MeV)= "
83 << etrans <<
" Eexnew= " << newExcEnergy
84 <<
" Ebond= " << bond_energy <<
G4endl;
87 etrans += bond_energy;
122 G4double mom = (emass > 0.0) ? std::sqrt((energy - emass)*(energy + emass))
130 energy =
std::max(ecm - energy, mass);
146 G4cout <<
"G4GammaTransition::SampleTransition : " << *result <<
G4endl;
158 G4cout <<
"G4GammaTransition::SampleDirection : 2J1= " << twoJ1
159 <<
" 2J2= " << twoJ2 <<
" ratio= " << ratio
160 <<
" mp= " << mp <<
G4endl;
177 G4int L0 = 0, Lp = 0;
189 G4double sinTheta = std::sqrt((1.-cosTheta)*(1.+cosTheta));
190 fDirection.
set(sinTheta*std::cos(phi),sinTheta*std::sin(phi),cosTheta);