G4EMDissociationSpectrum Class Reference

#include <G4EMDissociationSpectrum.hh>

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

Detailed Description

Definition at line 67 of file G4EMDissociationSpectrum.hh.

Constructor & Destructor Documentation

G4EMDissociationSpectrum::G4EMDissociationSpectrum

(

)

Definition at line 69 of file G4EMDissociationSpectrum.cc.

{
  bessel = new G4Bessel();
}

G4EMDissociationSpectrum::~G4EMDissociationSpectrum

(

)

Definition at line 75 of file G4EMDissociationSpectrum.cc.

{
  delete bessel;
}

Member Function Documentation

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:

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 caller graph for this function:

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:

---

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

• source/geant4.10.03.p03/source/processes/hadronic/cross_sections/include/G4EMDissociationSpectrum.hh
• source/geant4.10.03.p03/source/processes/hadronic/cross_sections/src/G4EMDissociationSpectrum.cc