G4PhotoNuclearCrossSection Class Reference

#include <G4PhotoNuclearCrossSection.hh>

Inheritance diagram for G4PhotoNuclearCrossSection:

Public Member Functions
	G4PhotoNuclearCrossSection ()
virtual	~G4PhotoNuclearCrossSection ()
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *particle, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
virtual G4bool	IsElementApplicable (const G4DynamicParticle *particle, G4int Z, const G4Material *)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *dynPart, G4int Z, G4int A, const G4Isotope *, const G4Element *elm, const G4Material *mat)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *)
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, 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	CrossSectionDescription (std::ostream &) const
virtual G4int	GetVerboseLevel () const
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

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 42 of file G4PhotoNuclearCrossSection.hh.

Constructor & Destructor Documentation

G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::G4PhotoNuclearCrossSection

(

)

Definition at line 1494 of file G4PhotoNuclearCrossSection.cc.

 : G4VCrossSectionDataSet(Default_Name()), lastZ(0), lastSig(0), lastGDR(0),
   lastHEN(0), lastE(0), lastTH(0), lastSP(0),
   deuteron_GDR(0), deuteron_HR(0), deuteron_TH(0), deuteron_SP(0),
   triton_GDR(0), triton_HR(0), triton_TH(0), triton_SP(0),
   he3_GDR(0), he3_HR(0), he3_TH(0), he3_SP(0),
   mNeut(G4NucleiProperties::GetNuclearMass(1,0)),
   mProt(G4NucleiProperties::GetNuclearMass(1,1))
{
  SetForAllAtomsAndEnergies(true);
  nistmngr = G4NistManager::Instance();
    
  for (G4int i = 0; i < 120; i++) {
    GDR.push_back(0);
    HEN.push_back(0);
    spA.push_back(0);
    eTH.push_back(0);
  }
}

~G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::~G4PhotoNuclearCrossSection ( )

virtual

Definition at line 1515 of file G4PhotoNuclearCrossSection.cc.

{
  std::vector<G4double*>::iterator posi;
  for (posi = GDR.begin(); posi < GDR.end(); posi++)
    { delete [] *posi; }
  GDR.clear();
  for (posi = HEN.begin(); posi < HEN.end(); posi++)
    { delete [] *posi; }
  HEN.clear();
}

Member Function Documentation

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1528 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().

GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1649 of file G4PhotoNuclearCrossSection.cc.

{
    const G4double Energy = aPart->GetKineticEnergy()/MeV;
    
    if (Energy<THmin) return 0.;
    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 = (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
    
    if(sigma<0.) return 0.;
    return sigma*millibarn;
}

Referenced by GetIsoCrossSection().

GetIsoCrossSection()

G4double G4PhotoNuclearCrossSection::GetIsoCrossSection ( const G4DynamicParticle *  dynPart, G4int  Z, G4int  A, const G4Isotope *  , const G4Element *  elm, const G4Material *  mat  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1562 of file G4PhotoNuclearCrossSection.cc.

{
  const G4double Energy = aPart->GetKineticEnergy()/MeV;
  if (Energy < THmin) return 0.;
 
  G4double sigma;
  G4double lE;
  if (Z == 1 && A == 2) { 
    // init the XS table if need be
    if (deuteron_GDR == NULL) {
      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 == NULL) {
      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 == NULL) {
      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 GetElementCrossSection(aPart, Z, mat);
  }
 
  if(sigma < 0.) sigma = 0.;
  return sigma*millibarn;
}

IsElementApplicable()

G4bool G4PhotoNuclearCrossSection::IsElementApplicable ( const G4DynamicParticle *  particle, G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1553 of file G4PhotoNuclearCrossSection.cc.

{
  return true;
}

IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1539 of file G4PhotoNuclearCrossSection.cc.

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

---

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

• G4PhotoNuclearCrossSection.hh
• G4PhotoNuclearCrossSection.cc