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

Member Function Documentation

GetClosestApproach()

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

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

Referenced by G4EMDissociation::ApplyYourself(), and G4EMDissociationCrossSection::GetElementCrossSection().

GetGeneralE1Spectrum()

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

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

Referenced by G4EMDissociationCrossSection::GetCrossSectionForProjectile().

GetGeneralE2Spectrum()

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

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

Referenced by G4EMDissociationCrossSection::GetCrossSectionForProjectile().

---

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

• G4EMDissociationSpectrum.hh
• G4EMDissociationSpectrum.cc