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

References wavefunctionP().

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

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

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

References wavefunctionR().

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

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

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

References derivWavefunctionR(), G4INCL::Math::twoPi, and wavefunctionR().

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

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

Definition at line 142 of file G4INCLDeuteronDensity.cc.

References G4INCL::Math::max().

Referenced by derivDensityR().

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

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

Definition at line 162 of file G4INCLDeuteronDensity.cc.

References G4INCL::PhysicalConstants::hc.

Referenced by densityP().

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

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

Definition at line 122 of file G4INCLDeuteronDensity.cc.

References G4INCL::Math::max().

Referenced by densityR(), and derivDensityR().

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