34 #define INCLXX_IN_GEANT4_MODE 1
54 namespace NuclearPotential {
81 vProton = theProtonFermiEnergy + theProtonSeparationEnergy;
89 vNeutron = theNeutronFermiEnergy + theNeutronSeparationEnergy;
91 const G4double separationEnergyDeltaPlusPlus = 2.*theProtonSeparationEnergy - theNeutronSeparationEnergy;
95 const G4double separationEnergyDeltaMinus = 2.*theNeutronSeparationEnergy - theProtonSeparationEnergy;
110 if(asy>0.11)
vLambda = 56.549-678.73*asy+4905.35*std::pow(asy,2.)-9789.1*std::pow(asy,3.);
131 separationEnergy[
KMinus] = 2.*theNeutronSeparationEnergy - theProtonSeparationEnergy-theLambdaSeparationEnergy;
151 INCL_DEBUG(
"Table of separation energies [MeV] for A=" << theA <<
", Z=" <<
theZ <<
":" <<
'\n'
152 <<
" proton: " << separationEnergy[
Proton] <<
'\n'
153 <<
" neutron: " << separationEnergy[
Neutron] <<
'\n'
155 <<
" delta+: " << separationEnergy[
DeltaPlus] <<
'\n'
156 <<
" delta0: " << separationEnergy[
DeltaZero] <<
'\n'
157 <<
" delta-: " << separationEnergy[
DeltaMinus] <<
'\n'
158 <<
" pi+: " << separationEnergy[
PiPlus] <<
'\n'
159 <<
" pi0: " << separationEnergy[
PiZero] <<
'\n'
160 <<
" pi-: " << separationEnergy[
PiMinus] <<
'\n'
161 <<
" eta: " << separationEnergy[
Eta] <<
'\n'
162 <<
" omega: " << separationEnergy[
Omega] <<
'\n'
163 <<
" etaprime:" << separationEnergy[
EtaPrime] <<
'\n'
164 <<
" photon: " << separationEnergy[
Photon] <<
'\n'
165 <<
" lambda: " << separationEnergy[
Lambda] <<
'\n'
166 <<
" sigmaplus: " << separationEnergy[
SigmaPlus] <<
'\n'
167 <<
" sigmazero: " << separationEnergy[
SigmaZero] <<
'\n'
168 <<
" sigmaminus: " << separationEnergy[
SigmaMinus] <<
'\n'
169 <<
" kplus: " << separationEnergy[
KPlus] <<
'\n'
170 <<
" kzero: " << separationEnergy[
KZero] <<
'\n'
171 <<
" kzerobar: " << separationEnergy[
KZeroBar] <<
'\n'
172 <<
" kminus: " << separationEnergy[
KMinus] <<
'\n'
173 <<
" kshort: " << separationEnergy[
KShort] <<
'\n'
174 <<
" klong: " << separationEnergy[
KLong] <<
'\n'
177 INCL_DEBUG(
"Table of Fermi energies [MeV] for A=" << theA <<
", Z=" <<
theZ <<
":" <<
'\n'
190 INCL_DEBUG(
"Table of Fermi momenta [MeV/c] for A=" << theA <<
", Z=" <<
theZ <<
":" <<
'\n'
258 INCL_ERROR(
"No potential computed for particle of type Cluster.");
262 INCL_ERROR(
"Trying to compute potential energy for an unknown particle.");
267 INCL_ERROR(
"There is no potential for this type of particle.");