34 #define INCLXX_IN_GEANT4_MODE 1
49 if(!(*p)->isNucleon())
continue;
54 const G4double pFermiSquared = pFermi*pFermi;
55 if((*p)->getMomentum().mag2() > pFermiSquared)
continue;
60 for(
ParticleIter i=particles.begin(), end=particles.end(); i!=end; ++i) {
61 if((*i)->getType() !=
t)
continue;
62 const G4double pmod2 = (*i)->getMomentum().mag2();
63 if(pmod2<pFermiSquared) nSea++;