34 #define INCLXX_IN_GEANT4_MODE 1
49 #ifndef G4INCLINUCLEARPOTENTIAL_HH
50 #define G4INCLINUCLEARPOTENTIAL_HH 1
60 namespace NuclearPotential {
110 std::map<ParticleType, G4double>::const_iterator i =
fermiEnergy.find(t);
132 std::map<ParticleType, G4double>::const_iterator i =
separationEnergy.find(t);
145 return std::sqrt(Tf*(Tf+2.*
m));
160 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.
double A(double temperature)
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.
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.