33 #define INCLXX_IN_GEANT4_MODE 1
48 if(!(*p)->isNucleon())
continue;
53 const G4double pFermiSquared = pFermi*pFermi;
54 if((*p)->getMomentum().mag2() > pFermiSquared)
continue;
59 for(
ParticleIter i=particles.begin(), end=particles.end(); i!=end; ++i) {
60 if((*i)->getType() != t)
continue;
61 const G4double pmod2 = (*i)->getMomentum().mag2();
62 if(pmod2<pFermiSquared) nSea++;
G4int getInitialZ() const
ParticleList const & getParticles() const
G4double getFermiMomentum(const Particle *const p) const
Return the Fermi momentum for a particle.
G4int getInitialA() const
G4bool isBlocked(ParticleList const &, Nucleus const *const)
NuclearPotential::INuclearPotential const * getPotential() const
Getter for thePotential.
ParticleList::const_iterator ParticleIter