#include <G4RToEConvForPositron.hh>

Inheritance diagram for G4RToEConvForPositron:

Public Member Functions
virtual G4double	Convert (const G4double rangeCut, const G4Material *material)

	G4RToEConvForPositron ()

const G4ParticleDefinition *	GetParticleType () const

G4int	GetVerboseLevel () const

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

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

void	SetVerboseLevel (G4int value)

virtual	~G4RToEConvForPositron ()

Static Public Member Functions
static G4double	GetHighEdgeEnergy ()

static G4double	GetLowEdgeEnergy ()

static G4double	GetMaxEnergyCut ()

static void	SetEnergyRange (const G4double lowedge, const G4double highedge)

static void	SetMaxEnergyCut (const G4double value)

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

Protected Attributes
G4int	fPDG = 0

G4bool	isFirstInstance = false

const G4ParticleDefinition *	theParticle = nullptr

G4int	verboseLevel = 1

Static Protected Attributes
static G4double	Emax = 0.0

static G4double	Emin = 0.0

static std::vector< G4double > *	Energy = nullptr

static G4int	Nbin = 350

static G4int	NbinPerDecade = 50

Private Member Functions
G4double	ConvertForElectron (const G4double rangeCut, const G4Material *material)

G4double	ConvertForGamma (const G4double rangeCut, const G4Material *material)

G4double	LiniearInterpolation (const G4double e1, const G4double e2, const G4double r1, const G4double r2, const G4double r)

Static Private Member Functions
static void	FillEnergyVector (const G4double emin, const G4double emax)

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

References G4ParticleTable::FindParticle(), G4VRangeToEnergyConverter::fPDG, G4cout, G4endl, G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGEncoding(), G4VRangeToEnergyConverter::GetVerboseLevel(), and G4VRangeToEnergyConverter::theParticle.

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

References CLHEP::electron_mass_c2, G4Exp(), G4Log(), G4Pow::GetInstance(), CLHEP::GeV, CLHEP::keV, G4Pow::logZ(), CLHEP::MeV, CLHEP::twopi_mc2_rcl2, and Z.

◆ Convert()

G4double G4VRangeToEnergyConverter::Convert	(	const G4double	rangeCut,
		const G4Material *	material
	)

virtualinherited

Reimplemented in G4RToEConvForProton.

Definition at line 87 of file G4VRangeToEnergyConverter.cc.

{
#ifdef G4VERBOSE
  if (GetVerboseLevel()>3)
  {
    G4cout << "G4VRangeToEnergyConverter::Convert() - ";
    G4cout << "Convert for " << material->GetName() 
       << " with Range Cut " << rangeCut/mm << "[mm]" << G4endl;
  }
#endif
 
  G4double cut = 0.0;
  if(fPDG == 22) 
  {  
    cut = ConvertForGamma(rangeCut, material);
  }
  else 
  {
    cut = ConvertForElectron(rangeCut, material);
 
    const G4double tune = 0.025*CLHEP::mm*CLHEP::g/CLHEP::cm3;
    const G4double lowen = 30.*CLHEP::keV; 
    if(cut < lowen)
    {
      //  corr. should be switched on smoothly   
      cut /= (1.+(1.-cut/lowen)*tune/(rangeCut*material->GetDensity())); 
    }
  }
 
  cut = std::max(Emin, std::min(cut, Emax));
  return cut;
}

References CLHEP::cm3, G4VRangeToEnergyConverter::ConvertForElectron(), G4VRangeToEnergyConverter::ConvertForGamma(), G4VRangeToEnergyConverter::Emax, G4VRangeToEnergyConverter::Emin, G4VRangeToEnergyConverter::fPDG, CLHEP::g, G4cout, G4endl, G4VRangeToEnergyConverter::GetVerboseLevel(), CLHEP::keV, eplot::material, G4INCL::Math::max(), G4INCL::Math::min(), CLHEP::mm, and mm.

Referenced by G4ProductionCutsTable::ConvertRangeToEnergy(), and G4ProductionCutsTable::UpdateCoupleTable().

◆ ConvertForElectron()

G4double G4VRangeToEnergyConverter::ConvertForElectron	(	const G4double	rangeCut,
		const G4Material *	material
	)

privateinherited

Definition at line 225 of file G4VRangeToEnergyConverter.cc.

{
  const G4ElementVector* elm = material->GetElementVector();
  const G4double* dens = material->GetAtomicNumDensityVector();
 
  // fill absorption length vector
  G4int nelm = material->GetNumberOfElements();
  G4double dedx1 = 0.0;
  G4double dedx2 = 0.0;
  G4double range1 = 0.0;
  G4double range2 = 0.0;
  G4double e1 = 0.0;
  G4double e2 = 0.0;
  G4double range = 0.;
  for (G4int i=0; i<Nbin; ++i)
  {
    e2 = (*Energy)[i];
    dedx2 = 0.0;
    for (G4int j=0; j<nelm; ++j)
    {
      dedx2 += dens[j]*ComputeValue((*elm)[j]->GetZasInt(), e2); 
    }
    range += (dedx1 + dedx2 > 0.0) ? 2*(e2 - e1)/(dedx1 + dedx2) : 0.0;
    range2 = range;
    if(range2 < rangeCut)
    {
      e1 = e2;
      dedx1 = dedx2;
      range1 = range2;
    }
    else
    {
      break;
    }
  }
  return LiniearInterpolation(e1, e2, range1, range2, rangeCut);
}

References G4VRangeToEnergyConverter::ComputeValue(), e1, e2, G4VRangeToEnergyConverter::LiniearInterpolation(), eplot::material, and G4VRangeToEnergyConverter::Nbin.

Referenced by G4VRangeToEnergyConverter::Convert().

◆ ConvertForGamma()

G4double G4VRangeToEnergyConverter::ConvertForGamma	(	const G4double	rangeCut,
		const G4Material *	material
	)

privateinherited

Definition at line 188 of file G4VRangeToEnergyConverter.cc.

{
  const G4ElementVector* elm = material->GetElementVector();
  const G4double* dens = material->GetAtomicNumDensityVector();
 
  // fill absorption length vector
  G4int nelm = material->GetNumberOfElements();
  G4double range1 = 0.0;
  G4double range2 = 0.0;
  G4double e1 = 0.0;
  G4double e2 = 0.0;
  for (G4int i=0; i<Nbin; ++i)
  {
    e2 = (*Energy)[i];
    G4double sig = 0.;
    
    for (G4int j=0; j<nelm; ++j)
    {
      sig += dens[j]*ComputeValue((*elm)[j]->GetZasInt(), e2); 
    }
    range2 = (sig > 0.0) ? 5./sig : DBL_MAX;
    if(i == 0 || range2 < rangeCut)
    {
      e1 = e2;
      range1 = range2;
    }
    else
    {
      break;
    }
  }
  return LiniearInterpolation(e1, e2, range1, range2, rangeCut);
}

References G4VRangeToEnergyConverter::ComputeValue(), DBL_MAX, e1, e2, G4VRangeToEnergyConverter::LiniearInterpolation(), eplot::material, and G4VRangeToEnergyConverter::Nbin.

Referenced by G4VRangeToEnergyConverter::Convert().

◆ FillEnergyVector()

void G4VRangeToEnergyConverter::FillEnergyVector	(	const G4double	emin,
		const G4double	emax
	)

staticprivateinherited

Definition at line 161 of file G4VRangeToEnergyConverter.cc.

{
  if(emin != Emin || emax != Emax) 
  {
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&theMutex);
    if(emin != Emin || emax != Emax) 
    { 
#endif
      Emin = emin;
      Emax = emax;
      Nbin = NbinPerDecade*static_cast<G4int>(std::log10(emax/emin));
      Energy->resize(Nbin + 1);
      (*Energy)[0] = emin;
      (*Energy)[Nbin] = emax;
      G4double fact = G4Log(emax/emin)/Nbin;
      for(G4int i=1; i<Nbin; ++i) { (*Energy)[i] = emin*G4Exp(i * fact); }
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&theMutex);
#endif
  }
}

References G4VRangeToEnergyConverter::Emax, emax, G4VRangeToEnergyConverter::Emin, G4VRangeToEnergyConverter::Energy, G4Exp(), G4Log(), G4MUTEXLOCK, G4MUTEXUNLOCK, G4VRangeToEnergyConverter::Nbin, and G4VRangeToEnergyConverter::NbinPerDecade.

Referenced by G4VRangeToEnergyConverter::G4VRangeToEnergyConverter(), G4VRangeToEnergyConverter::SetEnergyRange(), and G4VRangeToEnergyConverter::SetMaxEnergyCut().

◆ GetHighEdgeEnergy()

G4double G4VRangeToEnergyConverter::GetHighEdgeEnergy ( )

staticinherited

Definition at line 139 of file G4VRangeToEnergyConverter.cc.

{
  return Emax;
}

References G4VRangeToEnergyConverter::Emax.

Referenced by G4ProductionCutsTable::GetHighEdgeEnergy().

◆ GetLowEdgeEnergy()

G4double G4VRangeToEnergyConverter::GetLowEdgeEnergy ( )

staticinherited

Definition at line 133 of file G4VRangeToEnergyConverter.cc.

{
  return Emin;
}

References G4VRangeToEnergyConverter::Emin.

Referenced by G4ProductionCutsTable::GetLowEdgeEnergy().

◆ GetMaxEnergyCut()

G4double G4VRangeToEnergyConverter::GetMaxEnergyCut ( )

staticinherited

Definition at line 146 of file G4VRangeToEnergyConverter.cc.

{
  return Emax;
}

References G4VRangeToEnergyConverter::Emax.

Referenced by G4ProductionCutsTable::GetMaxEnergyCut().

◆ GetParticleType()

const G4ParticleDefinition * G4VRangeToEnergyConverter::GetParticleType ( ) const

inlineinherited

Definition at line 140 of file G4VRangeToEnergyConverter.hh.

{
  return theParticle;
}

References G4VRangeToEnergyConverter::theParticle.

◆ GetVerboseLevel()

G4int G4VRangeToEnergyConverter::GetVerboseLevel ( ) const

inlineinherited

Definition at line 134 of file G4VRangeToEnergyConverter.hh.

{
  return verboseLevel;
}

References G4VRangeToEnergyConverter::verboseLevel.

Referenced by G4VRangeToEnergyConverter::Convert(), G4RToEConvForProton::Convert(), G4RToEConvForElectron::G4RToEConvForElectron(), G4RToEConvForGamma::G4RToEConvForGamma(), G4RToEConvForPositron(), and G4RToEConvForProton::G4RToEConvForProton().

◆ LiniearInterpolation()

G4double G4VRangeToEnergyConverter::LiniearInterpolation	(	const G4double	e1,
		const G4double	e2,
		const G4double	r1,
		const G4double	r2,
		const G4double	r
	)

inlineprivateinherited

Definition at line 145 of file G4VRangeToEnergyConverter.hh.

{
  return (r1 == r2) ? e1 : e1 + (e2 - e1)*(r - r1)/(r2 - r1);
}

References e1, and e2.

Referenced by G4VRangeToEnergyConverter::ConvertForElectron(), and G4VRangeToEnergyConverter::ConvertForGamma().

◆ operator!=()

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

deleteinherited

◆ operator==()

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

deleteinherited

◆ SetEnergyRange()

void G4VRangeToEnergyConverter::SetEnergyRange	(	const G4double	lowedge,
		const G4double	highedge
	)

staticinherited

Definition at line 122 of file G4VRangeToEnergyConverter.cc.

{
  G4double ehigh = std::min(Emax, highedge);
  if(ehigh > lowedge)
  {
    FillEnergyVector(lowedge, ehigh);
  }  
}

References G4VRangeToEnergyConverter::Emax, G4VRangeToEnergyConverter::FillEnergyVector(), and G4INCL::Math::min().

Referenced by G4ProductionCutsTable::SetEnergyRange().

◆ SetMaxEnergyCut()

void G4VRangeToEnergyConverter::SetMaxEnergyCut ( const G4double value )

staticinherited

Definition at line 152 of file G4VRangeToEnergyConverter.cc.

{
  if(value > Emin)
  {
    FillEnergyVector(Emin, value);
  }
}

References G4VRangeToEnergyConverter::Emin, and G4VRangeToEnergyConverter::FillEnergyVector().

Referenced by G4ProductionCutsTable::SetMaxEnergyCut().

◆ SetVerboseLevel()

void G4VRangeToEnergyConverter::SetVerboseLevel ( G4int value )

inlineinherited

Definition at line 128 of file G4VRangeToEnergyConverter.hh.

{
  verboseLevel = value;
}

References G4VRangeToEnergyConverter::verboseLevel.

Referenced by G4ProductionCutsTable::SetVerboseLevel(), and G4ProductionCutsTable::UpdateCoupleTable().

Field Documentation

◆ Emax

G4double G4VRangeToEnergyConverter::Emax = 0.0

staticprotectedinherited

Definition at line 111 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::Convert(), G4VRangeToEnergyConverter::FillEnergyVector(), G4VRangeToEnergyConverter::GetHighEdgeEnergy(), G4VRangeToEnergyConverter::GetMaxEnergyCut(), G4VRangeToEnergyConverter::SetEnergyRange(), and G4VRangeToEnergyConverter::~G4VRangeToEnergyConverter().

◆ Emin

G4double G4VRangeToEnergyConverter::Emin = 0.0

staticprotectedinherited

Definition at line 110 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::Convert(), G4VRangeToEnergyConverter::FillEnergyVector(), G4VRangeToEnergyConverter::GetLowEdgeEnergy(), G4VRangeToEnergyConverter::SetMaxEnergyCut(), and G4VRangeToEnergyConverter::~G4VRangeToEnergyConverter().

◆ Energy

std::vector< G4double > * G4VRangeToEnergyConverter::Energy = nullptr

staticprotectedinherited

Definition at line 112 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::FillEnergyVector(), G4VRangeToEnergyConverter::G4VRangeToEnergyConverter(), and G4VRangeToEnergyConverter::~G4VRangeToEnergyConverter().

◆ fPDG

G4int G4VRangeToEnergyConverter::fPDG = 0

protectedinherited

Definition at line 118 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::Convert(), G4RToEConvForElectron::G4RToEConvForElectron(), G4RToEConvForGamma::G4RToEConvForGamma(), G4RToEConvForPositron(), and G4RToEConvForProton::G4RToEConvForProton().

◆ isFirstInstance

G4bool G4VRangeToEnergyConverter::isFirstInstance = false

protectedinherited

Definition at line 120 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::G4VRangeToEnergyConverter(), and G4VRangeToEnergyConverter::~G4VRangeToEnergyConverter().

◆ Nbin

G4int G4VRangeToEnergyConverter::Nbin = 350

staticprotectedinherited

Definition at line 114 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::ConvertForElectron(), G4VRangeToEnergyConverter::ConvertForGamma(), G4VRangeToEnergyConverter::FillEnergyVector(), and G4VRangeToEnergyConverter::G4VRangeToEnergyConverter().

◆ NbinPerDecade

G4int G4VRangeToEnergyConverter::NbinPerDecade = 50

staticprotectedinherited

Definition at line 113 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::FillEnergyVector().

◆ theParticle

const G4ParticleDefinition* G4VRangeToEnergyConverter::theParticle = nullptr

protectedinherited

Definition at line 116 of file G4VRangeToEnergyConverter.hh.

Referenced by G4RToEConvForElectron::G4RToEConvForElectron(), G4RToEConvForGamma::G4RToEConvForGamma(), G4RToEConvForPositron(), G4RToEConvForProton::G4RToEConvForProton(), and G4VRangeToEnergyConverter::GetParticleType().

◆ verboseLevel

G4int G4VRangeToEnergyConverter::verboseLevel = 1

protectedinherited

Definition at line 119 of file G4VRangeToEnergyConverter.hh.

Referenced by G4VRangeToEnergyConverter::GetVerboseLevel(), and G4VRangeToEnergyConverter::SetVerboseLevel().

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

source/processes/cuts/include/G4RToEConvForPositron.hh
source/processes/cuts/src/G4RToEConvForPositron.cc

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Static Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ G4RToEConvForPositron()

◆ ~G4RToEConvForPositron()

Member Function Documentation

◆ ComputeValue()

◆ Convert()

◆ ConvertForElectron()

◆ ConvertForGamma()

◆ FillEnergyVector()

◆ GetHighEdgeEnergy()

◆ GetLowEdgeEnergy()

◆ GetMaxEnergyCut()

◆ GetParticleType()

◆ GetVerboseLevel()

◆ LiniearInterpolation()

◆ operator!=()

◆ operator==()

◆ SetEnergyRange()

◆ SetMaxEnergyCut()

◆ SetVerboseLevel()

Field Documentation

◆ Emax

◆ Emin

◆ Energy

◆ fPDG

◆ isFirstInstance

◆ Nbin

◆ NbinPerDecade

◆ theParticle

◆ verboseLevel