G4INCL::DeuteronDensity Namespace Reference

Static class for the deuteron density. More...

Functions
G4double	densityR (const G4double r)
	PDF for a nucleon in r space.
G4double	derivDensityR (const G4double r)
	First derivative of the r-space density function.
G4double	densityP (const G4double p)
	PDF for a nucleon in p space.
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

densityP()

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);
    }

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

densityR()

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);
    }

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

derivDensityR()

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;
    }

derivWavefunctionR()

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;
    }

Referenced by derivDensityR().

wavefunctionP()

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;
    }

Referenced by densityP().

wavefunctionR()

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;
    }

Referenced by densityR(), and derivDensityR().