G4EMDissociationSpectrum Class Reference

#include <G4EMDissociationSpectrum.hh>

Collaboration diagram for G4EMDissociationSpectrum:

Public Member Functions
	G4EMDissociationSpectrum ()
	~G4EMDissociationSpectrum ()
G4double	GetGeneralE1Spectrum (G4double, G4double, G4double)
G4double	GetGeneralE2Spectrum (G4double, G4double, G4double)
G4double	GetClosestApproach (const G4double, const G4double, G4double, G4double, G4double)

Private Attributes
G4Bessel *	bessel

Detailed Description

Definition at line 67 of file G4EMDissociationSpectrum.hh.

Constructor & Destructor Documentation

G4EMDissociationSpectrum()

G4EMDissociationSpectrum::G4EMDissociationSpectrum

(

)

Definition at line 69 of file G4EMDissociationSpectrum.cc.

{
  bessel = new G4Bessel();
}

~G4EMDissociationSpectrum()

G4EMDissociationSpectrum::~G4EMDissociationSpectrum

(

)

Definition at line 75 of file G4EMDissociationSpectrum.cc.

{
  delete bessel;
}

Here is the call graph for this function:

Member Function Documentation

GetClosestApproach()

G4double G4EMDissociationSpectrum::GetClosestApproach	(	const G4double	AP,
		const G4double	ZP,
		G4double	AT,
		G4double	ZT,
		G4double	b
	)

Definition at line 112 of file G4EMDissociationSpectrum.cc.

{
  G4double bsq     = b * b;
  G4double gg       = 1.0/std::sqrt(1-bsq);
  G4double AProot3 = G4Pow::GetInstance()->powA(AP,1.0/3.0);
  G4double ATroot3 = G4Pow::GetInstance()->powA(AT,1.0/3.0);
  G4double bc      = 1.34 * fermi * (AProot3+ATroot3 - 0.75 *(1.0/AProot3+1.0/ATroot3));
//  G4double a0      = ZP * ZT * classic_electr_radius/bsq;
  G4double a0      = ZP * ZT * elm_coupling / (AT*AP*amu_c2/(AT+AP)) / bsq;
  G4double bmin    = 1.25 * bc + halfpi*a0/gg;
  return bmin;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetGeneralE1Spectrum()

G4double G4EMDissociationSpectrum::GetGeneralE1Spectrum	(	G4double	Eg,
		G4double	b,
		G4double	bmin
	)

Definition at line 82 of file G4EMDissociationSpectrum.cc.

{
  G4double b2 = b*b;
  G4double gg  = 1.0/std::sqrt(1.0-b2);
  G4double xi = Eg * bmin / gg / b / hbarc;
  G4double K0 = bessel->K0(xi);
  G4double K1 = bessel->K1(xi);
  G4double n  = 2.0 * fine_structure_const / pi / b2 / Eg *
    (xi*K0*K1 - xi*xi*b2/2.0*(K1*K1-K0*K0));
  return n;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetGeneralE2Spectrum()

G4double G4EMDissociationSpectrum::GetGeneralE2Spectrum	(	G4double	Eg,
		G4double	b,
		G4double	bmin
	)

Definition at line 96 of file G4EMDissociationSpectrum.cc.

{
  G4double b2 = b * b;
  G4double b4 = b2 * b2;
  G4double gg  = 1.0/std::sqrt(1.0-b2);
  G4double xi = Eg * bmin / gg / b / hbarc;
  G4double K0 = bessel->K0(xi);
  G4double K1 = bessel->K1(xi);
  G4double n  = 2.0 * fine_structure_const / pi / b4 / Eg *
    (2.0*(1.0-b2)*K1*K1 + xi*G4Pow::GetInstance()->powA((2.0-b2),2.0)*K0*K1 -
     xi*xi*b4/2.0*(K1*K1-K0*K0));
  return n;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

bessel

G4Bessel* G4EMDissociationSpectrum::bessel

private

Definition at line 78 of file G4EMDissociationSpectrum.hh.

---

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

• G4EMDissociationSpectrum.hh
• G4EMDissociationSpectrum.cc