#include <G4RToEConvForPositron.hh>

Inheritance diagram for G4RToEConvForPositron:

Public Member Functions
	G4RToEConvForPositron ()

virtual	~G4RToEConvForPositron ()

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 40 of file G4RToEConvForPositron.hh.

Constructor & Destructor Documentation

◆ G4RToEConvForPositron()

G4RToEConvForPositron::G4RToEConvForPositron ( )

explicit

Definition at line 42 of file G4RToEConvForPositron.cc.

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

◆ ~G4RToEConvForPositron()

G4RToEConvForPositron::~G4RToEConvForPositron ( )

virtual

Definition at line 63 of file G4RToEConvForPositron.cc.

64{}

Member Function Documentation

◆ ComputeValue()

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

finalprotectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 67 of file G4RToEConvForPositron.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 logtaul = G4Log(taul);
  const G4double taul12 = std::sqrt(taul);
  const G4double bremfactor = 0.1;
 
  G4double Zlog = G4Pow::GetInstance()->logZ(Z);
  G4double ionpot = 
    1.6e-5*CLHEP::MeV*G4Exp(0.9*Zlog)/CLHEP::electron_mass_c2;
  G4double ionpotlog = G4Log(ionpot);
 
  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 = 2.*logtaul - 
      (6.*taul+1.5*tsq-taul*(1.-tsq/3.)/t2
       -tsq*(0.5-tsq/12.)/(t2*t2))/(t1*t1);
    dEdx = (G4Log(2.*taul+4.)-2.*ionpotlog+f)/beta2;
    dEdx *= Z*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 = 2.*G4Log(tau) - (6.*tau+1.5*tsq-tau*(1.-tsq/3.)/t2
             -tsq*(0.5-tsq/12.)/(t2*t2))/(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 += cbrem*Z*(Z+1.)*bremfactor*tau/beta2;
  }
  return dEdx*CLHEP::twopi_mc2_rcl2;
}

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

Public Member Functions