G4PhotoNuclearCrossSection Class Reference

#include <G4PhotoNuclearCrossSection.hh>

Inheritance diagram for G4PhotoNuclearCrossSection:

Public Member Functions
	G4PhotoNuclearCrossSection ()
	~G4PhotoNuclearCrossSection () override
void	CrossSectionDescription (std::ostream &) const override
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr) override
G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr) override
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr) override
G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *) override
G4double	ComputeElementXSection (G4double energy, G4int Z)
G4double	ComputeIsoXSection (G4double energy, G4int Z, G4int A)
G4PhotoNuclearCrossSection &	operator= (const G4PhotoNuclearCrossSection &right)=delete
	G4PhotoNuclearCrossSection (const G4PhotoNuclearCrossSection &)=delete
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
bool	ForAllAtomsAndEnergies () const
void	SetForAllAtomsAndEnergies (G4bool val)
const G4String &	GetName () const
void	SetName (const G4String &nam)
G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete
	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel {0}
G4String	name

Detailed Description

Definition at line 42 of file G4PhotoNuclearCrossSection.hh.

Constructor & Destructor Documentation

G4PhotoNuclearCrossSection() [1/2]

G4PhotoNuclearCrossSection::G4PhotoNuclearCrossSection

(

)

Definition at line 1486 of file G4PhotoNuclearCrossSection.cc.

 : G4VCrossSectionDataSet(Default_Name()), 
   mNeut(G4NucleiProperties::GetNuclearMass(1,0)),
   mProt(G4NucleiProperties::GetNuclearMass(1,1))
{
  SetForAllAtomsAndEnergies(true);
  nistmngr = G4NistManager::Instance();
    
  GDR.resize(120, nullptr);
  HEN.resize(120, nullptr);
  spA.resize(120, 0.0);
  eTH.resize(120, 0.0);
}

Referenced by G4PhotoNuclearCrossSection(), and operator=().

~G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::~G4PhotoNuclearCrossSection ( )

override

Definition at line 1501 of file G4PhotoNuclearCrossSection.cc.

{
  for (auto const & ptr : GDR) { delete [] ptr; }
  for (auto const & ptr : HEN) { delete [] ptr; }
 
  delete [] deuteron_GDR;
  delete [] deuteron_HR;
  delete [] triton_GDR;
  delete [] triton_HR;
  delete [] he3_GDR;
  delete [] he3_HR;
}

G4PhotoNuclearCrossSection() [2/2]

G4PhotoNuclearCrossSection::G4PhotoNuclearCrossSection ( const G4PhotoNuclearCrossSection & )

delete

Member Function Documentation

ComputeElementXSection()

G4double G4PhotoNuclearCrossSection::ComputeElementXSection	(	G4double	energy,
		G4int	Z )

Definition at line 1645 of file G4PhotoNuclearCrossSection.cc.

{
  if (Energy <= THmin) { return 0.0; }
  G4int ZZ = std::min(Z, nL - 1);
  
    G4double sigma=0.;
    
    if(ZZ!=lastZ)   // Otherwise the set of parameters is ready
    {
        lastZ = ZZ;             // The last Z of calculated nucleus
        if(GDR[ZZ]) // Calculated nuclei in DB
        {
            lastGDR=GDR[ZZ];   // Pointer to prepared GDR cross sections
            lastHEN=HEN[ZZ];   // Pointer to prepared High Energy cross sections
            lastTH =eTH[ZZ];   // Energy Threshold
            lastSP =spA[ZZ];   // Shadowing coefficient for UHE
        }
        else
        {
            G4double Aa = nistmngr->GetAtomicMassAmu(ZZ); // average A
            G4int N = static_cast<G4int>(Aa) - ZZ;
            
            G4double lnA=G4Log(Aa); // The nucleus is not found in DB. It is new.
            if(Aa==1.) lastSP=1.;      // The Reggeon shadowing (A=1)
            else lastSP=Aa*(1.-shc*lnA);      // The Reggeon shadowing
            lastTH=ThresholdEnergy(ZZ, N); // Energy Threshold
            lastGDR = new G4double[nL];  // Allocate memory for the new
            // GDR cross sections
            lastHEN = new G4double[nH];  // Allocate memory for the new
            // HEN cross sections
            G4int er=GetFunctions(Aa,lastGDR,lastHEN); // new ZeroPosition and
            // filling of the functions
            if(er<1) G4cerr << "***G4PhotoNucCrossSection::GetCrossSection: A="
                << Aa << " failed" << G4endl;
            
            GDR[ZZ]=lastGDR;     // added GDR, found by AH 10/7/02
            HEN[ZZ]=lastHEN;     // added HEN, found by AH 10/7/02
            eTH[ZZ]=lastTH;      // Threshold Energy
            spA[ZZ]=lastSP;      // Pomeron Shadowing
        } // End of creation of the new set of parameters
    }// End of parameters udate
    
    //
    // =================== NOW the Magic Formula ===================
    //
    if (Energy<lastTH)
    {
        lastE=Energy;
        lastSig=0.;
        return 0.;
    }
    else if (Energy<Emin)   // GDR region (approximated in E, not in lnE)
    {
        sigma=EquLinearFit(Energy,nL,THmin,dE,lastGDR);
    }
    else if (Energy<Emax)                     // High Energy region
    {
        G4double lE=G4Log(Energy);
        sigma=EquLinearFit(lE,nH,milE,dlE,lastHEN);
    }
    else               // UHE region (calculation, but not so frequent)
    {
        G4double lE=G4Log(Energy);
        sigma=lastSP*(poc*(lE-pos)+shd*G4Exp(-reg*lE));
    }
    // End of "sigma" calculation
    
    sigma = std::max(sigma, 0.);
    return sigma*millibarn;
}

Referenced by ComputeIsoXSection(), and GetElementCrossSection().

ComputeIsoXSection()

G4double G4PhotoNuclearCrossSection::ComputeIsoXSection	(	G4double	energy,
		G4int	Z,
		G4int	A )

Definition at line 1557 of file G4PhotoNuclearCrossSection.cc.

{
  if (Energy <= THmin) { return 0.0; }
  G4double sigma;
  G4double lE;
  if (Z == 1 && A == 2) { 
    // init the XS table if need be
    if (deuteron_GDR == nullptr) {
      deuteron_TH = ThresholdEnergy(1,1); // threshold calculation is correct
      deuteron_GDR = new G4double[nL];    // direct copies
      for (G4int i = 0; i < nL; ++i) deuteron_GDR[i] = SL[0][i];   // A = 2 -> SL0
      deuteron_HR = new G4double[nH];
      for (G4int i = 0; i < nH; ++i) deuteron_HR[i] = SH[1][i];    // A = 2 -> SH1
      deuteron_SP = 1;                    // as would be assigned for proton
    }
    if (Energy < deuteron_TH) {
      sigma = 0.;
    } else if (Energy < Emin) {      // GDR region (approximated in E, not in lnE)
      sigma = EquLinearFit(Energy,nL,THmin,dE,deuteron_GDR);
    } else if (Energy < Emax) {      // High Energy region
      lE = G4Log(Energy);
      sigma = EquLinearFit(lE,nH,milE,dlE,deuteron_HR);
    } else {                         // Very high energy region
      lE = G4Log(Energy);
      sigma = deuteron_SP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
    }
 
  } else if (Z == 1 && A == 3) {
    if (triton_GDR == nullptr) {
      triton_TH = ThresholdEnergy(1,2);
      triton_GDR = new G4double[nL];      // same as for deuteron since no A = 3 entry
      for (G4int i = 0; i < nL; ++i) triton_GDR[i] = SL[0][i];     // A = 3 -> SL0
      triton_HR = new G4double[nH];
      for (G4int i = 0; i < nH; ++i) triton_HR[i] = SH[2][i];      // A = 3 -> SH2
      triton_SP = 1;
    }
    if (Energy < triton_TH) {
      sigma = 0.;
    } else if (Energy < Emin) {      // GDR region
      sigma = EquLinearFit(Energy,nL,THmin,dE,triton_GDR);
    } else if (Energy < Emax) {      // High Energy region
      lE = G4Log(Energy);
      sigma = EquLinearFit(lE,nH,milE,dlE,triton_HR);
    } else {
      lE = G4Log(Energy);
      sigma = triton_SP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
    }
 
  } else if (Z == 2 && A == 3) {
    if (he3_GDR == nullptr) {
      he3_TH = ThresholdEnergy(2,1);
      he3_GDR = new G4double[nL]; // same as for deuteron since no A = 3 entry
      for (G4int i = 0 ; i < nL ; ++i) he3_GDR[i] = SL[0][i]; // A = 3 -> SL0
      he3_HR = new G4double[nH];
      for (G4int i = 0 ; i < nH ; ++i) he3_HR[i] = SH[2][i];  // A = 3 -> SH2
      he3_SP = 2;
    }
    if (Energy < he3_TH) {
      sigma = 0.;
    } else if (Energy < Emin) {      // GDR region
      sigma = EquLinearFit(Energy,nL,THmin,dE,he3_GDR);
    } else if (Energy < Emax) {      // High Energy region
      lE = G4Log(Energy);
      sigma = EquLinearFit(lE,nH,milE,dlE,he3_HR);
    } else {
      lE = G4Log(Energy);
      sigma = he3_SP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
    }
 
  } else {
    return ComputeElementXSection(Energy, Z);
  }
 
  if(sigma < 0.) sigma = 0.;
  return sigma*millibarn;
}

Referenced by GetIsoCrossSection().

CrossSectionDescription()

void G4PhotoNuclearCrossSection::CrossSectionDescription ( std::ostream & outFile ) const

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1516 of file G4PhotoNuclearCrossSection.cc.

{
  outFile << "G4PhotoNuclearCrossSection provides the total inelastic\n"
          << "cross section for photon interactions with nuclei.  The\n"
          << "cross section is a parameterization of data which covers\n"
          << "all incident gamma energies.\n";
}

Default_Name()

static const char * G4PhotoNuclearCrossSection::Default_Name ( )

inlinestatic

Definition at line 49 of file G4PhotoNuclearCrossSection.hh.

{return "PhotoNuclearXS";}

Referenced by G4ElectroVDNuclearModel::G4ElectroVDNuclearModel(), and G4PhotoNuclearCrossSection().

GetElementCrossSection()

G4double G4PhotoNuclearCrossSection::GetElementCrossSection ( const G4DynamicParticle * aPart, G4int Z, const G4Material *  )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1638 of file G4PhotoNuclearCrossSection.cc.

{
  return ComputeElementXSection(aPart->GetKineticEnergy(), Z);
}

GetIsoCrossSection()

G4double G4PhotoNuclearCrossSection::GetIsoCrossSection ( const G4DynamicParticle * aPart, G4int Z, G4int A, const G4Isotope * iso = nullptr, const G4Element * elm = nullptr, const G4Material * mat = nullptr )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1547 of file G4PhotoNuclearCrossSection.cc.

{
  return ComputeIsoXSection(aPart->GetKineticEnergy(), Z, A);
}

IsElementApplicable()

G4bool G4PhotoNuclearCrossSection::IsElementApplicable ( const G4DynamicParticle * , G4int Z, const G4Material * mat = nullptr )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1538 of file G4PhotoNuclearCrossSection.cc.

{
  return true;
}

IsIsoApplicable()

G4bool G4PhotoNuclearCrossSection::IsIsoApplicable ( const G4DynamicParticle * , G4int Z, G4int A, const G4Element * elm = nullptr, const G4Material * mat = nullptr )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1527 of file G4PhotoNuclearCrossSection.cc.

{
  // explicitly allow deuterium, tritium, and He3
  return ((Z == 1 && A == 2) || (Z == 1 && A == 3) || (Z == 2 && A == 3)); 
}

operator=()

G4PhotoNuclearCrossSection & G4PhotoNuclearCrossSection::operator= ( const G4PhotoNuclearCrossSection & right )

delete

---

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

• G4PhotoNuclearCrossSection.hh
• G4PhotoNuclearCrossSection.cc