#include <G4RToEConvForElectron.hh>

Inheritance diagram for G4RToEConvForElectron:

Public Member Functions
	G4RToEConvForElectron ()

virtual	~G4RToEConvForElectron ()

Public Member Functions inherited from G4VRangeToEnergyConverter
	G4VRangeToEnergyConverter ()

virtual	~G4VRangeToEnergyConverter ()

	G4VRangeToEnergyConverter (const G4VRangeToEnergyConverter &r)=delete

G4VRangeToEnergyConverter &	operator= (const G4VRangeToEnergyConverter &r)=delete

G4bool	operator== (const G4VRangeToEnergyConverter &r) const =delete

G4bool	operator!= (const G4VRangeToEnergyConverter &r) const =delete

virtual G4double	Convert (const G4double rangeCut, const G4Material *material)

const G4ParticleDefinition *	GetParticleType () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

Protected Member Functions
G4double	ComputeValue (const G4int Z, const G4double kinEnergy) final

virtual G4double	ComputeValue (const G4int Z, const G4double kinEnergy)=0

Additional Inherited Members
Static Public Member Functions inherited from G4VRangeToEnergyConverter
static void	SetEnergyRange (const G4double lowedge, const G4double highedge)

static G4double	GetLowEdgeEnergy ()

static G4double	GetHighEdgeEnergy ()

static G4double	GetMaxEnergyCut ()

static void	SetMaxEnergyCut (const G4double value)

Protected Attributes inherited from G4VRangeToEnergyConverter
const G4ParticleDefinition *	theParticle = nullptr

G4int	fPDG = 0

Detailed Description

Definition at line 43 of file G4RToEConvForElectron.hh.

Constructor & Destructor Documentation

◆ G4RToEConvForElectron()

G4RToEConvForElectron::G4RToEConvForElectron ( )

explicit

Definition at line 42 of file G4RToEConvForElectron.cc.

  : G4VRangeToEnergyConverter()
{    
  theParticle = G4ParticleTable::GetParticleTable()->FindParticle("e-");
  if (theParticle == nullptr )
  {
#ifdef G4VERBOSE
    if (GetVerboseLevel()>0)
    {
      G4cout << "G4RToEConvForElectron::G4RToEConvForElectron() - ";
      G4cout << "Electron is not defined !!" << G4endl;
    }
#endif
  }
  else
  {
    fPDG = theParticle->GetPDGEncoding();
  }
}

◆ ~G4RToEConvForElectron()

G4RToEConvForElectron::~G4RToEConvForElectron ( )

virtual

Definition at line 63 of file G4RToEConvForElectron.cc.

64{}

Member Function Documentation

◆ ComputeValue()

G4double G4RToEConvForElectron::ComputeValue	(	const G4int	Z,
		const G4double	kinEnergy
	)

finalprotectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 67 of file G4RToEConvForElectron.cc.

{
  const G4double cbr1=0.02, cbr2=-5.7e-5, cbr3=1., cbr4=0.072;
  const G4double Tlow=10.*CLHEP::keV, Thigh=1.*CLHEP::GeV;
  const G4double taul = Tlow/CLHEP::electron_mass_c2;
  const G4double log05 = G4Log(0.5);
  const G4double taul12 = std::sqrt(taul);
  const G4double bremfactor = 0.1;
 
  const G4double Zlog = G4Pow::GetInstance()->logZ(Z);
  const G4double ionpot = 
    1.6e-5*CLHEP::MeV*G4Exp(0.9*Zlog)/CLHEP::electron_mass_c2;
  const G4double ionpotlog = G4Log(ionpot);
 
  const G4double tau = kinEnergy/CLHEP::electron_mass_c2;
  G4double dEdx = 0.0;
 
  if(tau<taul)
  {
    G4double t1 = taul+1.;
    G4double t2 = taul+2.;
    G4double tsq = taul*taul;
    G4double beta2 = taul*t2/(t1*t1);
    G4double f = 1.-beta2+G4Log(tsq/2.)
               +(0.5+0.25*tsq+(1.+2.*taul)*log05)/(t1*t1);
    dEdx = Z*(G4Log(2.*taul+4.)-2.*ionpotlog+f)/beta2;
    dEdx *= taul12/std::sqrt(tau);
  }
  else
  {
    G4double t1 = tau+1.;
    G4double t2 = tau+2.;
    G4double tsq = tau*tau;
    G4double beta2 = tau*t2/(t1*t1);
    G4double f = 1.-beta2+G4Log(tsq/2.)
                   +(0.5+0.25*tsq+(1.+2.*tau)*log05)/(t1*t1);
    dEdx = Z*(G4Log(2.*tau+4.)-2.*ionpotlog+f)/beta2;
 
    // loss from bremsstrahlung follows
    G4double cbrem = (cbr1+cbr2*Z)*(cbr3+cbr4*G4Log(kinEnergy/Thigh));
    dEdx += Z*(Z+1)*cbrem*bremfactor*tau/beta2;
  }
 
  return dEdx*CLHEP::twopi_mc2_rcl2;
}

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

Public Member Functions