#include <G4RToEConvForGamma.hh>

Inheritance diagram for G4RToEConvForGamma:

Public Member Functions
	G4RToEConvForGamma ()

virtual	~G4RToEConvForGamma ()

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 Types
using	G4CrossSectionTable = G4LossTable

Protected Types inherited from G4VRangeToEnergyConverter
using	G4LossVector = G4PhysicsLogVector

using	G4RangeVector = G4PhysicsLogVector

using	G4LossTable = G4PhysicsTable

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

virtual void	BuildRangeVector (const G4Material aMaterial, G4RangeVector rangeVector)

void	BuildAbsorptionLengthVector (const G4Material aMaterial, G4RangeVector rangeVector)

G4double	ComputeCrossSection (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	Z = -1.0

G4double	s200keV = 0.0

G4double	s1keV = 0.0

G4double	tmin = 0.0

G4double	tlow = 0.0

G4double	smin = 0.0

G4double	slow = 0.0

G4double	cmin = 0.0

G4double	clow = 0.0

G4double	chigh = 0.0

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)

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

Member Typedef Documentation

◆ G4CrossSectionTable

using G4RToEConvForGamma::G4CrossSectionTable = G4LossTable

protected

Definition at line 55 of file G4RToEConvForGamma.hh.

Constructor & Destructor Documentation

◆ G4RToEConvForGamma()

G4RToEConvForGamma::G4RToEConvForGamma ( )

Definition at line 41 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
  } 
}

◆ ~G4RToEConvForGamma()

G4RToEConvForGamma::~G4RToEConvForGamma ( )

virtual

Definition at line 57 of file G4RToEConvForGamma.cc.

58{

59}

Member Function Documentation

◆ BuildAbsorptionLengthVector()

void G4RToEConvForGamma::BuildAbsorptionLengthVector	(	const G4Material *	aMaterial,
		G4RangeVector *	rangeVector
	)

protected

Definition at line 64 of file G4RToEConvForGamma.cc.

{
  // fill the absorption length vector for this material
  // absorption length is defined here as:
  //   absorption length = 5./ macroscopic absorption cross section
  //
  const G4CrossSectionTable* aCrossSectionTable
        = (G4CrossSectionTable*)(theLossTable);
  const G4ElementVector* elementVector
        = aMaterial->GetElementVector();
  const G4double* atomicNumDensityVector
        = aMaterial->GetAtomicNumDensityVector();
 
  // fill absorption length vector
  G4int NumEl = aMaterial->GetNumberOfElements();
  G4double absorptionLengthMax = 0.0;
  for (std::size_t ibin=0; ibin<size_t(TotBin); ++ibin)
  {
    G4double SIGMA = 0.;
    for (std::size_t iel=0; iel<size_t(NumEl); ++iel)
    {
      G4int IndEl = (*elementVector)[iel]->GetIndex();
      SIGMA += atomicNumDensityVector[iel]
             * (*((*aCrossSectionTable)[IndEl]))[ibin];
    }
    //  absorption length=5./SIGMA
    absorptionLengthVector->PutValue(ibin, 5./SIGMA);
    if (absorptionLengthMax < 5./SIGMA )
      absorptionLengthMax = 5./SIGMA;
  }
}

Referenced by BuildRangeVector().

◆ BuildRangeVector()

void G4RToEConvForGamma::BuildRangeVector	(	const G4Material *	aMaterial,
		G4RangeVector *	rangeVector
	)

inlineprotectedvirtual

Reimplemented from G4VRangeToEnergyConverter.

Definition at line 89 of file G4RToEConvForGamma.hh.

{
  BuildAbsorptionLengthVector(aMaterial, rangeVector);
}

◆ ComputeCrossSection()

G4double G4RToEConvForGamma::ComputeCrossSection	(	G4double	AtomicNumber,
		G4double	KineticEnergy
	)

protected

Definition at line 101 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.*keV;
  const  G4double t200keV = 200.*keV;
  const  G4double t100MeV = 100.*MeV;
 
  // Compute Z dependent quantities in the case of a new AtomicNumber
  if(std::abs(AtomicNumber-Z)>0.1)
  {
    Z = AtomicNumber;
    G4double Zsquare = Z*Z;
    G4double Zlog = std::log(Z);
    G4double Zlogsquare = Zlog*Zlog;
 
    s200keV = (0.2651-0.1501*Zlog+0.02283*Zlogsquare)*Zsquare;
    tmin = (0.552+218.5/Z+557.17/Zsquare)*MeV;
    smin = (0.01239+0.005585*Zlog-0.000923*Zlogsquare)*std::exp(1.5*Zlog);
    cmin = std::log(s200keV/smin)
         /(std::log(tmin/t200keV)*std::log(tmin/t200keV));
    tlow = 0.2*std::exp(-7.355/std::sqrt(Z))*MeV;
    slow = s200keV
         * std::exp(0.042*Z*std::log(t200keV/tlow)*std::log(t200keV/tlow));
    s1keV = 300.*Zsquare;
    clow = std::log(s1keV/slow)/std::log(tlow/t1keV);
    chigh = (7.55e-5-0.0542e-5*Z)*Zsquare*Z/std::log(t100MeV/tmin);
  }
 
  // Calculate the cross-section (using an approximate empirical formula)
  G4double xs;
  if ( KineticEnergy<tlow )
  {
    if(KineticEnergy<t1keV) xs = slow*std::exp(clow*std::log(tlow/t1keV));
    else  xs = slow*std::exp(clow*std::log(tlow/KineticEnergy));
  }
  else if ( KineticEnergy<t200keV )
  {
    xs = s200keV
       * std::exp(0.042*Z*std::log(t200keV/KineticEnergy)
                         *std::log(t200keV/KineticEnergy));
  }
  else if( KineticEnergy<tmin )
  {
    xs = smin
       * std::exp(cmin*std::log(tmin/KineticEnergy)
                      *std::log(tmin/KineticEnergy));
  }
  else
  {
    xs = smin + chigh*std::log(KineticEnergy/tmin);
  }
  return xs * barn;
}

Referenced by ComputeLoss().

◆ ComputeLoss()

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

inlineprotectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 82 of file G4RToEConvForGamma.hh.

{
  return ComputeCrossSection(AtomicNumber,KineticEnergy);
}

Member Data Documentation

◆ chigh

G4double G4RToEConvForGamma::chigh = 0.0

protected

Definition at line 74 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ clow

G4double G4RToEConvForGamma::clow = 0.0

protected

Definition at line 74 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ cmin

G4double G4RToEConvForGamma::cmin = 0.0

protected

Definition at line 74 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ s1keV

G4double G4RToEConvForGamma::s1keV = 0.0

protected

Definition at line 71 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ s200keV

G4double G4RToEConvForGamma::s200keV = 0.0

protected

Definition at line 71 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ slow

G4double G4RToEConvForGamma::slow = 0.0

protected

Definition at line 73 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ smin

G4double G4RToEConvForGamma::smin = 0.0

protected

Definition at line 73 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ tlow

G4double G4RToEConvForGamma::tlow = 0.0

protected

Definition at line 72 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ tmin

G4double G4RToEConvForGamma::tmin = 0.0

protected

Definition at line 72 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

◆ Z

G4double G4RToEConvForGamma::Z = -1.0

protected

Definition at line 70 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

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

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ G4CrossSectionTable

Constructor & Destructor Documentation

◆ G4RToEConvForGamma()

◆ ~G4RToEConvForGamma()

Member Function Documentation

◆ BuildAbsorptionLengthVector()

◆ BuildRangeVector()

◆ ComputeCrossSection()

◆ ComputeLoss()

Member Data Documentation

◆ chigh

◆ clow

◆ cmin

◆ s1keV

◆ s200keV

◆ slow

◆ smin

◆ tlow

◆ tmin

◆ Z