33 #define INCLXX_IN_GEANT4_MODE 1
46 cellSize(
std::pow(2.38*4.5*Math::
pi,1./6.)*
std::sqrt(PhysicalConstants::
hc))
55 if( !(*p)->isNucleon() )
continue;
66 const G4double rbl = pbl * r0/pFermi;
75 if(rs - maxVolR > rdeq) {
79 if(rs + maxVolR > rdeq) {
80 vol = vol * 0.5 * (rdeq - rs + maxVolR) / maxVolR;
88 for(
ParticleIter it=particles.begin(), e=particles.end(); it!=e; ++it) {
90 if( (*it)->getID() == particle->
getID() )
continue;
92 if((*it)->getType() == particle->
getType()) {
95 if(dx2 > maxVolR * maxVolR)
continue;
99 if(dp2 > maxVolP * maxVolP)
continue;
106 if(blockingProbability > 1.0)
return 1.0;
107 else if(blockingProbability < 0.0)
return 0.0;
108 else return blockingProbability;
G4int getA() const
Returns the baryon number.
ParticleList const & getParticles() const
Return the list of "active" particles (i.e.
G4double getFermiMomentum(const Particle *const p) const
Return the Fermi momentum for a particle.
const G4double hc
[MeV*fm]
const G4INCL::ThreeVector & getMomentum() const
Get the momentum vector.
const G4double cellSize
Phase-space blocking cell.
G4bool isBlocked(ParticleList const &, Nucleus const *const)
G4double mag2() const
Get the square of the length.
G4double getNuclearRadius(const ParticleType t, const G4int A, const G4int Z)
G4int getZ() const
Returns the charge number.
G4double getBlockingProbability(Particle const *const, Nucleus const *const) const
NuclearPotential::INuclearPotential const * getPotential() const
Getter for thePotential.
G4INCL::ParticleType getType() const
Get the particle type.
const G4INCL::ThreeVector & getPosition() const
Set the position vector.
G4double shoot()
Generate flat distribution of random numbers.
G4double getUniverseRadius() const
Getter for theUniverseRadius.
G4double mag() const
Get the length of the vector.
ParticleList::const_iterator ParticleIter