94 G4double sintheta = std::sqrt(1.0 - costheta*costheta);
96 G4ThreeVector direction(sintheta*std::cos(phi),sintheta*std::sin(phi),
99 G4double KE = std::sqrt(cmMomentum*cmMomentum + alphaMass*alphaMass)
105 KE = std::sqrt(cmMomentum*cmMomentum + nucleusMass*nucleusMass) - nucleusMass;
134 <<
" with branching ratio " <<
GetBR() <<
"% and Q value "