33 #define INCLXX_IN_GEANT4_MODE 1
53 namespace NuclearPotential {
65 void NuclearPotentialIsospin::initialize() {
79 vProton = theProtonFermiEnergy + theProtonSeparationEnergy;
87 vNeutron = theNeutronFermiEnergy + theNeutronSeparationEnergy;
90 vDeltaZero = vNeutron;
91 vDeltaPlusPlus = 2.*vDeltaPlus - vDeltaZero;
92 vDeltaMinus = 2.*vDeltaZero - vDeltaPlus;
108 INCL_DEBUG(
"Table of separation energies [MeV] for A=" << theA <<
", Z=" <<
theZ <<
":" << std::endl
109 <<
" proton: " << separationEnergy[
Proton] << std::endl
110 <<
" neutron: " << separationEnergy[
Neutron] << std::endl
111 <<
" delta++: " << separationEnergy[
DeltaPlusPlus] << std::endl
112 <<
" delta+: " << separationEnergy[
DeltaPlus] << std::endl
113 <<
" delta0: " << separationEnergy[
DeltaZero] << std::endl
114 <<
" delta-: " << separationEnergy[
DeltaMinus] << std::endl
115 <<
" pi+: " << separationEnergy[
PiPlus] << std::endl
116 <<
" pi0: " << separationEnergy[
PiZero] << std::endl
117 <<
" pi-: " << separationEnergy[
PiMinus] << std::endl
120 INCL_DEBUG(
"Table of Fermi energies [MeV] for A=" << theA <<
", Z=" <<
theZ <<
":" << std::endl
129 INCL_DEBUG(
"Table of Fermi momenta [MeV/c] for A=" << theA <<
", Z=" <<
theZ <<
":" << std::endl
153 return vDeltaPlusPlus;
165 INCL_ERROR(
"No potential computed for particle of type Cluster.");
169 INCL_ERROR(
"Trying to compute potential energy for an unknown particle.");
174 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()
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)
G4INCL::ParticleType getType() const
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