G4INCL::DeuteronDensity Namespace Reference

Static class for the deuteron density. More...

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

Detailed Description

Static class for the deuteron density.

Function Documentation

G4double G4INCL::DeuteronDensity::densityP

(

const G4double

p

)

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 117 of file G4INCLDeuteronDensity.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4INCL::DeuteronDensity::densityR

(

const G4double

r

)

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 103 of file G4INCLDeuteronDensity.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4INCL::DeuteronDensity::derivDensityR

(

const G4double

r

)

First derivative of the r-space density function.

Parameters

r	distance from the deuteron centre [fm]

Returns

d|psi|^2/dr

Definition at line 109 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;
    }

Here is the call graph for this function:

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

Definition at line 143 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;
    }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 163 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;
    }

Here is the caller graph for this function:

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

Definition at line 123 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;
    }

Here is the call graph for this function:

Here is the caller graph for this function: