34 #define INCLXX_IN_GEANT4_MODE 1
47 cellSize(std::pow(2.38*4.5*Math::
pi,1./6.)*std::sqrt(PhysicalConstants::
hc))
56 if( !(*p)->isNucleon() )
continue;
67 const G4double rbl = pbl * r0/pFermi;
76 if(rs - maxVolR > rdeq) {
80 if(rs + maxVolR > rdeq) {
81 vol = vol * 0.5 * (rdeq - rs + maxVolR) / maxVolR;
89 for(
ParticleIter it=particles.begin(),
e=particles.end(); it!=
e; ++it) {
91 if( (*it)->getID() == particle->
getID() )
continue;
93 if((*it)->getType() == particle->
getType()) {
96 if(dx2 > maxVolR * maxVolR)
continue;
100 if(dp2 > maxVolP * maxVolP)
continue;
107 if(blockingProbability > 1.0)
return 1.0;
108 else if(blockingProbability < 0.0)
return 0.0;
109 else return blockingProbability;
ParticleList::const_iterator ParticleIter
NuclearPotential::INuclearPotential const * getPotential() const
Getter for thePotential.
const G4INCL::ThreeVector & getPosition() const
ParticleList const & getParticles() const
G4bool isBlocked(ParticleList const &, Nucleus const *const)
const G4INCL::ThreeVector & getMomentum() const
G4double getBlockingProbability(Particle const *const, Nucleus const *const) const
G4double getFermiMomentum(const Particle *const p) const
Return the Fermi momentum for a particle.
const G4double cellSize
Phase-space blocking cell.
G4int getZ() const
Returns the charge number.
G4double getUniverseRadius() const
Getter for theUniverseRadius.
G4int getA() const
Returns the baryon number.
G4INCL::ParticleType getType() const
static constexpr double pi
G4double getNuclearRadius(const ParticleType t, const G4int A, const G4int Z)
const G4double hc
[MeV*fm]