#include <G4RToEConvForPositron.hh>

Inheritance diagram for G4RToEConvForPositron:

Public Member Functions
	G4RToEConvForPositron ()

virtual	~G4RToEConvForPositron ()

Public Member Functions inherited from G4VRangeToEnergyConverter
	G4VRangeToEnergyConverter ()

	G4VRangeToEnergyConverter (const G4VRangeToEnergyConverter &right)

G4VRangeToEnergyConverter &	operator= (const G4VRangeToEnergyConverter &right)

virtual	~G4VRangeToEnergyConverter ()

G4int	operator== (const G4VRangeToEnergyConverter &right) const

G4int	operator!= (const G4VRangeToEnergyConverter &right) const

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

const G4ParticleDefinition *	GetParticleType () const

const G4PhysicsTable *	GetLossTable () const

virtual void	Reset ()

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

Protected Member Functions
virtual G4double	ComputeLoss (G4double AtomicNumber, G4double KineticEnergy) const

Protected Member Functions inherited from G4VRangeToEnergyConverter
virtual void	BuildLossTable ()

virtual G4double	ComputeLoss (G4double AtomicNumber, G4double KineticEnergy) const =0

virtual void	BuildRangeVector (const G4Material aMaterial, G4RangeVector rangeVector)

G4double	ConvertCutToKineticEnergy (G4RangeVector *theRangeVector, G4double theCutInLength, size_t materialIndex) const

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

static G4double	GetLowEdgeEnergy ()

static G4double	GetHighEdgeEnergy ()

static G4double	GetMaxEnergyCut ()

static void	SetMaxEnergyCut (G4double value)

Protected Types inherited from G4VRangeToEnergyConverter
typedef G4PhysicsTable	G4LossTable

typedef G4PhysicsLogVector	G4LossVector

typedef G4PhysicsLogVector	G4RangeVector

Protected Attributes inherited from G4VRangeToEnergyConverter
G4double	fMaxEnergyCut

const G4ParticleDefinition *	theParticle

G4LossTable *	theLossTable

G4int	NumberOfElements

const G4int	TotBin

std::vector< G4RangeVector * >	fRangeVectorStore

Static Protected Attributes inherited from G4VRangeToEnergyConverter
static G4double	LowestEnergy = 0.99e-3*MeV

static G4double	HighestEnergy = 100.0e6*MeV

static G4double	MaxEnergyCut = 10.0*GeV

Detailed Description

Definition at line 51 of file G4RToEConvForPositron.hh.

Constructor & Destructor Documentation

◆ G4RToEConvForPositron()

G4RToEConvForPositron::G4RToEConvForPositron ( )

Definition at line 45 of file G4RToEConvForPositron.cc.

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

◆ ~G4RToEConvForPositron()

G4RToEConvForPositron::~G4RToEConvForPositron ( )

virtual

Definition at line 58 of file G4RToEConvForPositron.cc.

59{

60}

Member Function Documentation

◆ ComputeLoss()

G4double G4RToEConvForPositron::ComputeLoss	(	G4double	AtomicNumber,
		G4double	KineticEnergy
	)		const

protectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 68 of file G4RToEConvForPositron.cc.

{
  static G4double Z;  
  static G4double taul, ionpot, ionpotlog;
  const  G4double cbr1=0.02, cbr2=-5.7e-5, cbr3=1., cbr4=0.072;
  const  G4double Tlow=10.*keV, Thigh=1.*GeV;
  static G4double bremfactor = 0.1 ;
 
  G4double Mass = theParticle->GetPDGMass();       
  //  calculate dE/dx for electrons
  if( std::fabs(AtomicNumber-Z)>0.1 ) {
    Z = AtomicNumber;
    taul = Tlow/Mass;
    ionpot = 1.6e-5*MeV*std::exp(0.9*std::log(Z))/Mass;
    ionpotlog = std::log(ionpot);
  } 
 
  G4double tau = KineticEnergy/Mass;
  G4double dEdx;
 
 if(tau<taul)
  {
    G4double t1 = taul+1.;
    G4double t2 = taul+2.;
    G4double tsq = taul*taul;
    G4double beta2 = taul*t2/(t1*t1);
    G4double     f = 2.*std::log(taul)
                     -(6.*taul+1.5*tsq-taul*(1.-tsq/3.)/t2-tsq*(0.5-tsq/12.)/
                       (t2*t2))/(t1*t1);
    dEdx = (std::log(2.*taul+4.)-2.*ionpotlog+f)/beta2;
    dEdx = twopi_mc2_rcl2*Z*dEdx;
    G4double clow = dEdx*std::sqrt(taul);
    dEdx = clow/std::sqrt(KineticEnergy/Mass);
 
  } else {
   G4double t1 = tau+1.;
    G4double t2 = tau+2.;
    G4double tsq = tau*tau;
    G4double beta2 = tau*t2/(t1*t1);
    G4double f = 2.*std::log(tau)
                 - (6.*tau+1.5*tsq-tau*(1.-tsq/3.)/t2-tsq*(0.5-tsq/12.)/
                     (t2*t2))/(t1*t1);
    dEdx = (std::log(2.*tau+4.)-2.*ionpotlog+f)/beta2;
    dEdx = twopi_mc2_rcl2*Z*dEdx;
 
    // loss from bremsstrahlung follows
    G4double cbrem = (cbr1+cbr2*Z)
                       *(cbr3+cbr4*std::log(KineticEnergy/Thigh));
    cbrem = Z*(Z+1.)*cbrem*tau/beta2;
    cbrem *= bremfactor ;
    dEdx += twopi_mc2_rcl2*cbrem;
  }
  return dEdx;
}

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

Public Member Functions