33 #define INCLXX_IN_GEANT4_MODE 1
47 theMaximumRadius(
std::
min((*rpCorrelationTableProton)(1.), (*rpCorrelationTableNeutron)(1.))),
66 INCL_DEBUG(
"Interpolation table for proton local energy (A=" <<
theA <<
", Z=" <<
theZ <<
") initialised:"
70 <<
"Interpolation table for neutron local energy (A=" <<
theA <<
", Z=" <<
theZ <<
") initialised:"
74 <<
"Interpolation table for delta++ local energy (A=" <<
theA <<
", Z=" <<
theZ <<
") initialised:"
78 <<
"Interpolation table for delta+ local energy (A=" <<
theA <<
", Z=" <<
theZ <<
") initialised:"
82 <<
"Interpolation table for delta0 local energy (A=" <<
theA <<
", Z=" <<
theZ <<
") initialised:"
86 <<
"Interpolation table for delta- local energy (A=" <<
theA <<
", Z=" <<
theZ <<
") initialised:"
107 theMaximumRadius(rhs.theMaximumRadius),
108 theProtonNuclearRadius(rhs.theProtonNuclearRadius)
131 swap(temporaryDensity);
149 std::swap(
pFromR[DeltaPlusPlus], rhs.
pFromR[DeltaPlusPlus]);
156 const G4double theProtonRadius = 0.88;
G4double transmissionRadius[UnknownParticle]
G4double getMinPFromR(const ParticleType t, const G4double r) const
InverseInterpolationTable const * rFromP[UnknownParticle]
InverseInterpolationTable const * pFromR[UnknownParticle]
NuclearDensity(const G4int A, const G4int Z, InverseInterpolationTable const *const rpCorrelationTableProton, InverseInterpolationTable const *const rpCorrelationTableNeutron)
NuclearDensity & operator=(const NuclearDensity &rhs)
Assignment operator.
G4double getNuclearRadius(const ParticleType t, const G4int A, const G4int Z)
void copy(std::vector< T > &main, const std::vector< T > &data)
static const G4double A[nN]
Class for interpolating the inverse of a 1-dimensional function.
void swap(NuclearDensity &rhs)
Helper method for the assignment operator.
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double theMaximumRadius
void initializeTransmissionRadii()
Initialize the transmission radius.
void print(const std::vector< T > &data)
G4double theProtonNuclearRadius
Represents INCL4.5's R0 variable.
G4double getMaxRFromP(const ParticleType t, const G4double p) const
Get the maximum allowed radius for a given momentum.