34 #define INCLXX_IN_GEANT4_MODE 1
52 namespace DeuteronDensity {
56 const G4int coeffTableSize = 13;
59 const G4double coeff1[coeffTableSize] = {
76 const G4double coeff2[coeffTableSize] = {
106 return r*r*(sWave*sWave + dWave*dWave);
114 return (sWave*sWaveDeriv + dWave*dWaveDeriv) /
Math::twoPi;
120 return p*p*(sWave*sWave + dWave*dWave);
130 for(
G4int i=0; i<coeffTableSize; ++i) {
133 result += coeff1[i] * std::exp(-fmr);
135 result += coeff2[i] * std::exp(-fmr) * (1.+3./fmr+3./(fmr*fmr));
139 result *= normalisationR/
r;
150 for(
G4int i=0; i<coeffTableSize; ++i) {
153 result += coeff1[i] * std::exp(-fmr) * (fmr + 1.);
155 result += coeff2[i] * std::exp(-fmr) * (fmr + 4. + 9./fmr + 9./(fmr*fmr));
159 result *= -normalisationR/(r*
r);
169 for(
G4int i=0; i<coeffTableSize; ++i) {
171 fmq = q2 + alPlusI*alPlusI;
173 result += coeff1[i] / fmq;
175 result += coeff2[i] / fmq;
179 result *= normalisationP;