34 #define INCLXX_IN_GEANT4_MODE 1
54 namespace NuclearPotential {
80 vProton = theProtonFermiEnergy + theProtonSeparationEnergy;
88 vNeutron = theNeutronFermiEnergy + theNeutronSeparationEnergy;
90 const G4double separationEnergyDeltaPlusPlus = 2.*theProtonSeparationEnergy - theNeutronSeparationEnergy;
94 const G4double separationEnergyDeltaMinus = 2.*theNeutronSeparationEnergy - theProtonSeparationEnergy;
112 INCL_DEBUG(
"Table of separation energies [MeV] for A=" << theA <<
", Z=" <<
theZ <<
":" <<
'\n'
113 <<
" proton: " << separationEnergy[
Proton] <<
'\n'
114 <<
" neutron: " << separationEnergy[
Neutron] <<
'\n'
116 <<
" delta+: " << separationEnergy[
DeltaPlus] <<
'\n'
117 <<
" delta0: " << separationEnergy[
DeltaZero] <<
'\n'
118 <<
" delta-: " << separationEnergy[
DeltaMinus] <<
'\n'
119 <<
" pi+: " << separationEnergy[
PiPlus] <<
'\n'
120 <<
" pi0: " << separationEnergy[
PiZero] <<
'\n'
121 <<
" pi-: " << separationEnergy[
PiMinus] <<
'\n'
124 INCL_DEBUG(
"Table of Fermi energies [MeV] for A=" << theA <<
", Z=" <<
theZ <<
":" <<
'\n'
133 INCL_DEBUG(
"Table of Fermi momenta [MeV/c] for A=" << theA <<
", Z=" <<
theZ <<
":" <<
'\n'
169 INCL_ERROR(
"No potential computed for particle of type Cluster.");
173 INCL_ERROR(
"Trying to compute potential energy for an unknown particle.");
178 INCL_ERROR(
"There is no potential for this type of particle.");
virtual G4double computePotentialEnergy(const Particle *const p) const
std::map< ParticleType, G4double > fermiEnergy
virtual ~NuclearPotentialIsospin()
double A(double temperature)
const G4int theA
The mass number of the nucleus.
Isospin-dependent nuclear potential.
Isospin- and energy-independent nuclear potential.
G4ThreadLocal SeparationEnergyFn getSeparationEnergy
Static pointer to the separation-energy function.
G4double getINCLMass(const G4int A, const G4int Z)
Get INCL nuclear mass (in MeV/c^2)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4INCL::ParticleType getType() const
Get the particle type.
NuclearPotentialIsospin(const G4int A, const G4int Z, const G4bool pionPotential)
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.
std::map< ParticleType, G4double > fermiMomentum
G4double pow13(G4double x)
G4ThreadLocal FermiMomentumFn getFermiMomentum