79 throw G4HadronicException(__FILE__, __LINE__,
"G4NeutronHPNBodyPhaseSpace: Unknown ion case 2");
92 else if(theTotalCount==4)
96 else if(theTotalCount==5)
102 throw G4HadronicException(__FILE__, __LINE__,
"NeutronHP Phase-space distribution cannot cope with this number of particles");
105 G4double rand0 = Prob(max, maxE, theTotalCount);
112 testit = Prob(energy, maxE, theTotalCount);
114 while(rand > testit);
123 G4ThreeVector tempVector(mtot*sinth*std::cos(phi), mtot*sinth*std::sin(phi), mtot*std::cos(theta) );
126 result->
Lorentz(*result, -1.*aCMS);