G4INCL::DeuteronDensity Class Reference

Static class for the deuteron density. More...

#include <G4INCLDeuteronDensity.hh>

Static Public Member Functions
static G4double	densityR (const G4double r)
	PDF for a nucleon in r space.
static G4double	derivDensityR (const G4double r)
	First derivative of the r-space density function.
static G4double	densityP (const G4double p)
	PDF for a nucleon in p space.
static G4double	wavefunctionR (const G4int l, const G4double r)
static G4double	wavefunctionP (const G4int l, const G4double p)
static G4double	derivWavefunctionR (const G4int l, const G4double r)

Detailed Description

Static class for the deuteron density.

Definition at line 55 of file G4INCLDeuteronDensity.hh.

Member Function Documentation

densityP()

G4double G4INCL::DeuteronDensity::densityP ( const G4double  p )

static

PDF for a nucleon in p space.

The distribution is normalised to 1.

Parameters

p	nucleon momentum [MeV/c]

Returns

4 * pi * p^2 * d|psi|^2/dp

Definition at line 110 of file G4INCLDeuteronDensity.cc.

                                                     {
    const G4double sWave = wavefunctionP(0, p);
    const G4double dWave = wavefunctionP(2, p);
    return p*p*(sWave*sWave + dWave*dWave);
  }

Referenced by G4INCL::NuclearDensityFunctions::ParisP::operator()().

densityR()

G4double G4INCL::DeuteronDensity::densityR ( const G4double  r )

static

PDF for a nucleon in r space.

The distribution is normalised to 1.

Parameters

r	distance from the deuteron centre [fm]

Returns

4 * pi * r^2 * |psi|^2

Definition at line 96 of file G4INCLDeuteronDensity.cc.

                                                     {
    const G4double sWave = wavefunctionR(0, r);
    const G4double dWave = wavefunctionR(2, r);
    return r*r*(sWave*sWave + dWave*dWave);
  }

Referenced by G4INCL::NuclearDensityFunctions::ParisR::operator()().

derivDensityR()

G4double G4INCL::DeuteronDensity::derivDensityR ( const G4double  r )

static

First derivative of the r-space density function.

Parameters

r	distance from the deuteron centre [fm]

Returns

d|psi|^2/dr

Definition at line 102 of file G4INCLDeuteronDensity.cc.

                                                          {
    const G4double sWave = wavefunctionR(0, r);
    const G4double dWave = wavefunctionR(2, r);
    const G4double sWaveDeriv = derivWavefunctionR(0, r);
    const G4double dWaveDeriv = derivWavefunctionR(2, r);
    return (sWave*sWaveDeriv + dWave*dWaveDeriv) / Math::twoPi;
  }

derivWavefunctionR()

G4double G4INCL::DeuteronDensity::derivWavefunctionR ( const G4int  l, const G4double  r  )

static

Definition at line 136 of file G4INCLDeuteronDensity.cc.

                                                                                 {
// assert(l==0 || l==2); // only s- and d-waves in a deuteron
    const G4double r = 2. * std::max(theR, 1.e-4);
 
    G4double result = 0.;
    G4double fmr;
 
    for(G4int i=0; i<coeffTableSize; ++i) {
      fmr = r * (al+i);
      if(l==0) { // s-wave
        result += coeff1[i] * std::exp(-fmr) * (fmr + 1.);
      } else { // d-wave
        result += coeff2[i] * std::exp(-fmr) * (fmr + 4. + 9./fmr + 9./(fmr*fmr));
      }
    }
 
    result *= -normalisationR/(r*r);
    return result;
  }

Referenced by derivDensityR().

wavefunctionP()

G4double G4INCL::DeuteronDensity::wavefunctionP ( const G4int  l, const G4double  p  )

static

Definition at line 156 of file G4INCLDeuteronDensity.cc.

                                                                            {
// assert(l==0 || l==2); // only s- and d-waves in a deuteron
    const G4double q = theQ / PhysicalConstants::hc;
    const G4double q2 = q*q;
    G4double result=0.;
    G4double fmq, alPlusI;
    for(G4int i=0; i<coeffTableSize; ++i) {
      alPlusI = al+i;
      fmq = q2 + alPlusI*alPlusI;
      if(l==0) { // s-wave
        result += coeff1[i] / fmq;
      } else { // d-wave
        result += coeff2[i] / fmq;
      }
    }
 
    result *= normalisationP;
    return result;
  }

Referenced by densityP().

wavefunctionR()

G4double G4INCL::DeuteronDensity::wavefunctionR ( const G4int  l, const G4double  r  )

static

Definition at line 116 of file G4INCLDeuteronDensity.cc.

                                                                            {
// assert(l==0 || l==2); // only s- and d-waves in a deuteron
    const G4double r = 2. * std::max(theR, 1.e-4);
 
    G4double result = 0.;
    G4double fmr;
 
    for(G4int i=0; i<coeffTableSize; ++i) {
      fmr = r * (al+i);
      if(l==0) { // s-wave
        result += coeff1[i] * std::exp(-fmr);
      } else { // d-wave
        result += coeff2[i] * std::exp(-fmr) * (1.+3./fmr+3./(fmr*fmr));
      }
    }
 
    result *= normalisationR/r;
    return result;
  }

Referenced by densityR(), and derivDensityR().

---

The documentation for this class was generated from the following files:

• G4INCLDeuteronDensity.hh
• G4INCLDeuteronDensity.cc