34#define INCLXX_IN_GEANT4_MODE 1
52 namespace DeuteronDensity {
106 return r*r*(sWave*sWave + dWave*dWave);
114 return (sWave*sWaveDeriv + dWave*dWaveDeriv) /
Math::twoPi;
120 return p*p*(sWave*sWave + dWave*dWave);
133 result +=
coeff1[i] * std::exp(-fmr);
135 result +=
coeff2[i] * std::exp(-fmr) * (1.+3./fmr+3./(fmr*fmr));
153 result +=
coeff1[i] * std::exp(-fmr) * (fmr + 1.);
155 result +=
coeff2[i] * std::exp(-fmr) * (fmr + 4. + 9./fmr + 9./(fmr*fmr));
171 fmq = q2 + alPlusI*alPlusI;
173 result +=
coeff1[i] / fmq;
175 result +=
coeff2[i] / fmq;
Deuteron density in r and p according to the Paris potential.
const G4int coeffTableSize
const G4double coeff1[coeffTableSize]
Coefficients for the deuteron wave function.
const G4double coeff2[coeffTableSize]
Coefficients for the deuteron wave function.
const G4double normalisationR
Normalisation coefficient for the r-space deuteron wave function.
const G4double al
Mysterious coefficient that appears in the wavefunctions.
const G4double normalisationP
Normalisation coefficient for the p-space deuteron wave function.
G4double wavefunctionR(const G4int l, const G4double r)
G4double derivDensityR(const G4double r)
First derivative of the r-space density function.
G4double derivWavefunctionR(const G4int l, const G4double r)
G4double wavefunctionP(const G4int l, const G4double p)
G4double densityR(const G4double r)
PDF for a nucleon in r space.
G4double densityP(const G4double p)
PDF for a nucleon in p space.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
const G4double hc
[MeV*fm]