34 #define INCLXX_IN_GEANT4_MODE 1
184 maxGeneratedWeight(0.)
195 for(
size_t i=1; i<
wMaxNP; ++i) {
216 const G4int maxIter = 500;
224 }
while(++iter<maxIter && r*weightMax>weight);
226 #ifndef INCLXX_IN_GEANT4_MODE
228 INCL_WARN(
"Number of tries exceeded in PhaseSpaceRauboldLynch::generate()\n"
249 if(availableEnergy<0.)
250 availableEnergy = 0.;
263 eMMin += masses[i-1];
271 #ifndef INCLXX_IN_GEANT4_MODE
273 INCL_WARN(
"The requested number of particles (" <<
nParticles <<
") requires extrapolation the tables in PhaseSpaceRauboldLynch. YMMV." <<
'\n');
277 INCL_WARN(
"The requested available energy (" <<
availableEnergy <<
" MeV) requires extrapolation the tables in PhaseSpaceRauboldLynch. YMMV." <<
'\n');
284 const G4double correction = (*wMaxCorrection)(reducedSqrtS);
297 rnd[nParticles-1] = 1.;
298 std::sort(
rnd.begin()+1,
rnd.begin()+nParticles-1);
307 for(
size_t i=1; i<nParticles-1; ++i) {
313 weight *= momentumCM;
338 const G4double recoilE = std::sqrt(momentum.
mag2() + iM*iM);
339 boostV = -momentum/recoilE;
340 for(
size_t j=0; j<=i; ++j)
341 particles[j]->boost(boostV);