33 #define INCLXX_IN_GEANT4_MODE 1
48 #ifndef G4INCLINUCLEARPOTENTIAL_HH
49 #define G4INCLINUCLEARPOTENTIAL_HH 1
59 namespace NuclearPotential {
109 std::map<ParticleType, G4double>::const_iterator i =
fermiEnergy.find(t);
131 std::map<ParticleType, G4double>::const_iterator i =
separationEnergy.find(t);
144 return std::sqrt(Tf*(Tf+2.*
m));
159 std::map<ParticleType, G4double>::const_iterator i =
fermiMomentum.find(t);
G4double getMass() const
Get the cached particle mass.
const G4double eSquared
Coulomb conversion factor [MeV*fm].
std::map< ParticleType, G4double > fermiEnergy
G4double getFermiMomentum(const Particle *const p) const
Return the Fermi momentum for a particle.
static const G4double vPionDefault
G4bool isDelta() const
Is it a Delta?
G4double getFermiEnergy(const ParticleType t) const
Return the Fermi energy for a particle type.
G4double getFermiEnergy(const Particle *const p) const
Return the Fermi energy for a particle.
G4double getFermiMomentum(const ParticleType t) const
Return the Fermi momentum for a particle type.
void clearCache()
Clear the INuclearPotential cache.
const G4bool pionPotential
G4double getSeparationEnergy(const Particle *const p) const
Return the separation energy for a particle.
Deuteron density in r and p according to the Paris potential.
G4double getSeparationEnergy(const ParticleType t) const
Return the separation energy for a particle type.
const G4int theA
The mass number of the nucleus.
G4bool isOutOfWell() const
Check if the particle is out of its potential well.
static const G4double A[nN]
G4bool hasPionPotential() const
Do we have a pion potential?
INuclearPotential(const G4int A, const G4int Z, const G4bool pionPot)
G4INCL::ParticleType getType() const
Get the particle type.
std::map< ParticleType, G4double > separationEnergy
const G4int theZ
The charge number of the nucleus.
G4double computePionPotentialEnergy(const Particle *const p) const
Compute the potential energy for the given pion.
virtual ~INuclearPotential()
virtual G4double computePotentialEnergy(const Particle *const p) const =0
std::map< ParticleType, G4double > fermiMomentum
G4double pow13(G4double x)
INuclearPotential const * createPotential(const PotentialType type, const G4int theA, const G4int theZ, const G4bool pionPotential)
Create an INuclearPotential object.