G4RToEConvForPositron Class Reference

#include <G4RToEConvForPositron.hh>

Inheritance diagram for G4RToEConvForPositron:

Public Member Functions
	G4RToEConvForPositron ()
virtual	~G4RToEConvForPositron ()
Public Member Functions inherited from G4VRangeToEnergyConverter
	G4VRangeToEnergyConverter ()
	G4VRangeToEnergyConverter (const G4VRangeToEnergyConverter &r)
G4VRangeToEnergyConverter &	operator= (const G4VRangeToEnergyConverter &r)
virtual	~G4VRangeToEnergyConverter ()
G4bool	operator== (const G4VRangeToEnergyConverter &r) const
G4bool	operator!= (const G4VRangeToEnergyConverter &r) 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)
Protected Member Functions inherited from G4VRangeToEnergyConverter
virtual void	BuildLossTable ()
virtual G4double	ComputeLoss (G4double AtomicNumber, G4double KineticEnergy)=0
virtual void	BuildRangeVector (const G4Material *aMaterial, G4RangeVector *rangeVector)
G4double	ConvertCutToKineticEnergy (G4RangeVector *theRangeVector, G4double theCutInLength, std::size_t materialIndex) const

Protected Attributes
G4double	Mass = 0.0
G4double	Z = -1.0
G4double	taul = 0.0
G4double	ionpot = 0.0
G4double	ionpotlog = -1.0e-10
G4double	bremfactor = 0.1
Protected Attributes inherited from G4VRangeToEnergyConverter
G4double	fMaxEnergyCut = 0.0
const G4ParticleDefinition *	theParticle = nullptr
G4LossTable *	theLossTable = nullptr
G4int	NumberOfElements = 0
const G4int	TotBin = 300
std::vector< G4RangeVector * >	fRangeVectorStore

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
using	G4LossVector = G4PhysicsLogVector
using	G4RangeVector = G4PhysicsLogVector
using	G4LossTable = G4PhysicsTable
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 43 of file G4RToEConvForPositron.hh.

Constructor & Destructor Documentation

G4RToEConvForPositron()

G4RToEConvForPositron::G4RToEConvForPositron

(

)

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
  {
    Mass = theParticle->GetPDGMass();
  }
}

~G4RToEConvForPositron()

G4RToEConvForPositron::~G4RToEConvForPositron ( )

virtual

Definition at line 63 of file G4RToEConvForPositron.cc.

{ 
}

Member Function Documentation

ComputeLoss()

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

protectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 70 of file G4RToEConvForPositron.cc.

{
  const G4double cbr1=0.02, cbr2=-5.7e-5, cbr3=1., cbr4=0.072;
  const G4double Tlow=10.*keV, Thigh=1.*GeV;
 
  // 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;
}

Member Data Documentation

bremfactor

G4double G4RToEConvForPositron::bremfactor = 0.1

protected

Definition at line 63 of file G4RToEConvForPositron.hh.

Referenced by ComputeLoss().

ionpot

G4double G4RToEConvForPositron::ionpot = 0.0

protected

Definition at line 61 of file G4RToEConvForPositron.hh.

Referenced by ComputeLoss().

ionpotlog

G4double G4RToEConvForPositron::ionpotlog = -1.0e-10

protected

Definition at line 62 of file G4RToEConvForPositron.hh.

Referenced by ComputeLoss().

Mass

G4double G4RToEConvForPositron::Mass = 0.0

protected

Definition at line 58 of file G4RToEConvForPositron.hh.

Referenced by ComputeLoss(), and G4RToEConvForPositron().

taul

G4double G4RToEConvForPositron::taul = 0.0

protected

Definition at line 60 of file G4RToEConvForPositron.hh.

Referenced by ComputeLoss().

Z

G4double G4RToEConvForPositron::Z = -1.0

protected

Definition at line 59 of file G4RToEConvForPositron.hh.

Referenced by ComputeLoss().

---

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

• G4RToEConvForPositron.hh
• G4RToEConvForPositron.cc