#include <G4EMDissociationSpectrum.hh>

Public Member Functions
	G4EMDissociationSpectrum ()

G4double	GetClosestApproach (const G4double, const G4double, G4double, G4double, G4double)

G4double	GetGeneralE1Spectrum (G4double, G4double, G4double)

G4double	GetGeneralE2Spectrum (G4double, G4double, G4double)

	~G4EMDissociationSpectrum ()

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();
}

References bessel.

◆ ~G4EMDissociationSpectrum()

G4EMDissociationSpectrum::~G4EMDissociationSpectrum ( )

Definition at line 75 of file G4EMDissociationSpectrum.cc.

{
  delete bessel;
}

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

References a0, source.hepunit::amu_c2, anonymous_namespace{paraMaker.cc}::AP, source.hepunit::elm_coupling, fermi, G4Pow::GetInstance(), halfpi, and G4Pow::powA().

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

References bessel, source.hepunit::fine_structure_const, source.hepunit::hbarc, G4Bessel::K0(), G4Bessel::K1(), CLHEP::detail::n, and pi.

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