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 67 of file G4EMDissociationSpectrum.cc.

{
  bessel = new G4Bessel();
}

G4EMDissociationSpectrum::~G4EMDissociationSpectrum

(

)

Definition at line 73 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 110 of file G4EMDissociationSpectrum.cc.

References python.hepunit::amu_c2, test::b, python.hepunit::elm_coupling, python.hepunit::fermi, and python.hepunit::halfpi.

{
  G4double bsq     = b * b;
  G4double gg       = 1.0/std::sqrt(1-bsq);
  G4double AProot3 = std::pow(AP,1.0/3.0);
  G4double ATroot3 = std::pow(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;
}

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

Definition at line 80 of file G4EMDissociationSpectrum.cc.

References test::b, python.hepunit::fine_structure_const, python.hepunit::hbarc, n, and python.hepunit::pi.

Referenced by G4EMDissociationCrossSection::GetCrossSectionForProjectile().

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

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

Definition at line 94 of file G4EMDissociationSpectrum.cc.

References test::b, python.hepunit::fine_structure_const, python.hepunit::hbarc, n, and python.hepunit::pi.

Referenced by G4EMDissociationCrossSection::GetCrossSectionForProjectile().

{
  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*std::pow((2.0-b2),2.0)*K0*K1 -
     xi*xi*b4/2.0*(K1*K1-K0*K0));
  return n;
}

---

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

• G4EMDissociationSpectrum.hh
• G4EMDissociationSpectrum.cc