G4RToEConvForGamma Class Reference

#include <G4RToEConvForGamma.hh>

Inheritance diagram for G4RToEConvForGamma:

Public Member Functions
virtual G4double	Convert (const G4double rangeCut, const G4Material *material)
	G4RToEConvForGamma ()
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	~G4RToEConvForGamma ()

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 G4RToEConvForGamma.hh.

Constructor & Destructor Documentation

G4RToEConvForGamma()

G4RToEConvForGamma::G4RToEConvForGamma ( )

explicit

Definition at line 40 of file G4RToEConvForGamma.cc.

  : G4VRangeToEnergyConverter()
{    
  theParticle = G4ParticleTable::GetParticleTable()->FindParticle("gamma");
  if (theParticle == nullptr)
  {
#ifdef G4VERBOSE
    if (GetVerboseLevel()>0)
    {
      G4cout << " G4RToEConvForGamma::G4RToEConvForGamma() - ";
      G4cout << "Gamma 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.

~G4RToEConvForGamma()

G4RToEConvForGamma::~G4RToEConvForGamma ( )

virtual

Definition at line 61 of file G4RToEConvForGamma.cc.

{}

Member Function Documentation

ComputeValue()

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

finalprotectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 65 of file G4RToEConvForGamma.cc.

{
  // Compute the "absorption" cross-section of the photon "absorption".
  // Cross-section means here the sum of the cross-sections of the
  // pair production, Compton scattering and photoelectric processes
 
  const G4double t1keV = 1.*CLHEP::keV;
  const G4double t200keV = 200.*CLHEP::keV;
  const G4double t100MeV = 100.*CLHEP::MeV;
 
  G4double Zsquare = Z*Z;
  G4double Zlog = G4Pow::GetInstance()->logZ(Z);
  G4double Zlogsquare = Zlog*Zlog;
 
  G4double tmin = (0.552+218.5/Z+557.17/Zsquare)*CLHEP::MeV;
  G4double tlow = 0.2*G4Exp(-7.355/std::sqrt(Z))*CLHEP::MeV;
 
  G4double smin = (0.01239+0.005585*Zlog-0.000923*Zlogsquare)*G4Exp(1.5*Zlog);
  G4double s200keV = (0.2651-0.1501*Zlog+0.02283*Zlogsquare)*Zsquare;
 
  G4double cminlog = G4Log(tmin/t200keV); 
  G4double cmin = G4Log(s200keV/smin)/(cminlog*cminlog);
 
  G4double slowlog = G4Log(t200keV/tlow);
  G4double slow = s200keV * G4Exp(0.042*Z*slowlog*slowlog);
  G4double logtlow = G4Log(tlow/t1keV);
  G4double clow = G4Log(300.*Zsquare/slow)/logtlow;
  G4double chigh = (7.55e-5 - 0.0542e-5*Z)*Zsquare*Z/G4Log(t100MeV/tmin);
 
  // Calculate the cross-section (using an approximate empirical formula)
  G4double xs;
  if ( energy < tlow )
  {
    xs = (energy < t1keV) ? slow*G4Exp(clow*logtlow) :
      slow*G4Exp(clow*G4Log(tlow/energy));
  }
  else if ( energy < t200keV )
  {
    G4double x = G4Log(t200keV/energy);
    xs = s200keV * G4Exp(0.042*Z*x*x);
  }
  else if( energy<tmin )
  {
    const G4double x = G4Log(tmin/energy);
    xs = smin * G4Exp(cmin*x*x);
  }
  else
  {
    xs = smin + chigh*G4Log(energy/tmin);
  }
  return xs * CLHEP::barn;
}

References CLHEP::barn, G4INCL::KinematicsUtils::energy(), G4Exp(), G4Log(), G4Pow::GetInstance(), CLHEP::keV, G4Pow::logZ(), CLHEP::MeV, tlow, 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(), G4RToEConvForPositron::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(), G4RToEConvForPositron::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(), G4RToEConvForPositron::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/G4RToEConvForGamma.hh
• source/processes/cuts/src/G4RToEConvForGamma.cc