34 #define INCLXX_IN_GEANT4_MODE 1
50 G4double Fmax = std::sqrt(std::pow(xN/(xN+1.),N)/(xN+1.));
53 unsigned long loopCounter = 0;
54 const unsigned long maxLoopCounter = 10000000;
57 F = std::sqrt(std::pow(chi,N)*(1.-chi));
59 }
while (loopCounter<maxLoopCounter && Fmax*
Random::shoot() > F);
69 const size_t N = particles.size();
82 Particle *restParticle = particles.front();
87 for (ParticleList::reverse_iterator
p=particles.rbegin(); k>0; ++
p, --
k) {
95 (*p)->setMomentum(PFragCM);
96 (*p)->adjustEnergyFromMomentum();
97 restParticle->
setMass(restMass);