G4RToEConvForGamma Class Reference

#include <G4RToEConvForGamma.hh>

Inheritance diagram for G4RToEConvForGamma:

Public Member Functions
	G4RToEConvForGamma ()
virtual	~G4RToEConvForGamma ()
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 Types
typedef G4LossTable	G4CrossSectionTable
Protected Types inherited from G4VRangeToEnergyConverter
typedef G4PhysicsTable	G4LossTable
typedef G4PhysicsLogVector	G4LossVector
typedef G4PhysicsLogVector	G4RangeVector

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 ()
G4double	ConvertCutToKineticEnergy (G4RangeVector *theRangeVector, G4double theCutInLength, size_t materialIndex) const

Protected Attributes
G4double	Z
G4double	s200keV
G4double	s1keV
G4double	tmin
G4double	tlow
G4double	smin
G4double	slow
G4double	cmin
G4double	clow
G4double	chigh
Protected Attributes inherited from G4VRangeToEnergyConverter
G4double	fMaxEnergyCut
const G4ParticleDefinition *	theParticle
G4LossTable *	theLossTable
G4int	NumberOfElements
const G4int	TotBin
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 51 of file G4RToEConvForGamma.hh.

Member Typedef Documentation

typedef G4LossTable G4RToEConvForGamma::G4CrossSectionTable

protected

Definition at line 74 of file G4RToEConvForGamma.hh.

Constructor & Destructor Documentation

G4RToEConvForGamma::G4RToEConvForGamma

(

)

Definition at line 44 of file G4RToEConvForGamma.cc.

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

  : G4VRangeToEnergyConverter(),
    Z(-1),  
    s200keV(0.), s1keV(0.),
    tmin(0.),    tlow(0.), 
    smin(0.),    slow(0.),
    cmin(0.),    clow(0.), chigh(0.)
{    
  theParticle =  G4ParticleTable::GetParticleTable()->FindParticle("gamma");
  if (theParticle ==0) {
#ifdef G4VERBOSE
    if (GetVerboseLevel()>0) {
      G4cout << " G4RToEConvForGamma::G4RToEConvForGamma() ";
      G4cout << " Gamma is not defined !!" << G4endl;
    }
#endif
  } 
}

G4RToEConvForGamma::~G4RToEConvForGamma ( )

virtual

Definition at line 63 of file G4RToEConvForGamma.cc.

{ 
}

Member Function Documentation

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

protected

Definition at line 71 of file G4RToEConvForGamma.cc.

References G4Material::GetAtomicNumDensityVector(), G4Material::GetElementVector(), G4Material::GetNumberOfElements(), G4PhysicsVector::PutValue(), G4VRangeToEnergyConverter::theLossTable, and G4VRangeToEnergyConverter::TotBin.

Referenced by BuildRangeVector().

{
  // 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 (size_t ibin=0; ibin<size_t(TotBin); ibin++) {
    G4double SIGMA = 0. ;
    for (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;
  }
}

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

inlineprotectedvirtual

Reimplemented from G4VRangeToEnergyConverter.

Definition at line 98 of file G4RToEConvForGamma.hh.

References BuildAbsorptionLengthVector().

{
  BuildAbsorptionLengthVector(aMaterial, rangeVector);
}

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

protected

Definition at line 105 of file G4RToEConvForGamma.cc.

References python.hepunit::barn, chigh, clow, cmin, python.hepunit::keV, python.hepunit::MeV, s1keV, s200keV, slow, smin, tlow, tmin, and Z.

Referenced by ComputeLoss().

{
  //  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;
}

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

inlineprotectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 91 of file G4RToEConvForGamma.hh.

References ComputeCrossSection().

{
  return ComputeCrossSection(AtomicNumber,KineticEnergy);
}

Field Documentation

G4double G4RToEConvForGamma::chigh

protected

Definition at line 86 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::clow

protected

Definition at line 86 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::cmin

protected

Definition at line 86 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::s1keV

protected

Definition at line 83 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::s200keV

protected

Definition at line 83 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::slow

protected

Definition at line 85 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::smin

protected

Definition at line 85 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::tlow

protected

Definition at line 84 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::tmin

protected

Definition at line 84 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

G4double G4RToEConvForGamma::Z

protected

Definition at line 82 of file G4RToEConvForGamma.hh.

Referenced by ComputeCrossSection().

---

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

• G4RToEConvForGamma.hh
• G4RToEConvForGamma.cc