G4GammaNuclearXS Class Reference

#include <G4GammaNuclearXS.hh>

Inheritance diagram for G4GammaNuclearXS:

Public Member Functions
	G4GammaNuclearXS ()
	~G4GammaNuclearXS () override=default
G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *) final
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *, const G4Material *mat) final
G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr) final
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr) final
const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE) final
void	BuildPhysicsTable (const G4ParticleDefinition &) final
G4double	IsoCrossSection (const G4double ekin, const G4int Z, const G4int A)
G4double	ElementCrossSection (const G4double ekin, const G4int Z)
G4double	LowEnergyCrossSection (G4double ekin, G4int Z)
void	CrossSectionDescription (std::ostream &) const final
G4GammaNuclearXS &	operator= (const G4GammaNuclearXS &right)=delete
	G4GammaNuclearXS (const G4GammaNuclearXS &)=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 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 60 of file G4GammaNuclearXS.hh.

Constructor & Destructor Documentation

G4GammaNuclearXS() [1/2]

G4GammaNuclearXS::G4GammaNuclearXS

(

)

Definition at line 78 of file G4GammaNuclearXS.cc.

  : G4VCrossSectionDataSet(Default_Name()), gamma(G4Gamma::Gamma())
{
  verboseLevel = 0;
  if (verboseLevel > 0) {
    G4cout << "G4GammaNuclearXS::G4GammaNuclearXS Initialise for Z < " 
       << MAXZGAMMAXS << G4endl;
  }
  ggXsection = dynamic_cast<G4PhotoNuclearCrossSection*>
    (G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet("PhotoNuclearXS"));
  if (ggXsection == nullptr) {
    ggXsection = new G4PhotoNuclearCrossSection();
  }
  SetForAllAtomsAndEnergies(true);
 
  // full data set is uploaded once
  if (nullptr == data) { 
    data = new G4ElementData(MAXZGAMMAXS);
    data->SetName("gNuclear");
    for (G4int Z=1; Z<MAXZGAMMAXS; ++Z) {
      Initialise(Z);
    }
  }
}

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

~G4GammaNuclearXS()

G4GammaNuclearXS::~G4GammaNuclearXS ( )

overridedefault

G4GammaNuclearXS() [2/2]

G4GammaNuclearXS::G4GammaNuclearXS ( const G4GammaNuclearXS & )

delete

Member Function Documentation

BuildPhysicsTable()

void G4GammaNuclearXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 255 of file G4GammaNuclearXS.cc.

{
  if(verboseLevel > 1) {
    G4cout << "G4GammaNuclearXS::BuildPhysicsTable for " 
       << p.GetParticleName() << G4endl;
  }
  if(p.GetParticleName() != "gamma") { 
    G4ExceptionDescription ed;
    ed << p.GetParticleName() << " is a wrong particle type -"
       << " only gamma is allowed";
    G4Exception("G4GammaNuclearXS::BuildPhysicsTable(..)","had012",
        FatalException, ed, "");
    return; 
  }
 
  // prepare isotope selection
  const G4ElementTable* table = G4Element::GetElementTable();
  std::size_t nIso = temp.size();
  for (auto const & elm : *table ) {
    std::size_t n = elm->GetNumberOfIsotopes();
    if (n > nIso) { nIso = n; }
  }
  temp.resize(nIso, 0.0);   
}

Referenced by G4InterfaceToXS::G4InterfaceToXS().

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 103 of file G4GammaNuclearXS.cc.

{
  outFile << "G4GammaNuclearXS calculates the gamma nuclear\n"
          << "cross-section for GDR energy region on nuclei using "
      << "data from the high precision\n"
          << "IAEA photonuclear database (2019). Then liniear connection\n"
      << "implemented with previous CHIPS photonuclear model." << G4endl;
}

Default_Name()

static const char * G4GammaNuclearXS::Default_Name ( )

inlinestatic

Definition at line 68 of file G4GammaNuclearXS.hh.

{ return "GammaNuclearXS"; }

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

ElementCrossSection()

G4double G4GammaNuclearXS::ElementCrossSection	(	const G4double	ekin,
		const G4int	Z )

Definition at line 134 of file G4GammaNuclearXS.cc.

{
  // check cache
  const G4int Z = (ZZ < MAXZGAMMAXS) ? ZZ : MAXZGAMMAXS - 1;
  if(Z == fZ && ekin == fEkin) { return fXS; }
  fZ = Z;
  fEkin = ekin;
 
  auto pv = data->GetElementData(Z);
  const G4double limCHIPS1 = 25*CLHEP::MeV;
  const G4double limCHIPS2 = 16*CLHEP::MeV;
  if (pv == nullptr || 1 == Z || Z == 40 || Z == 74 ||
     (Z == 24 && ekin >= limCHIPS1) ||
     (Z == 39 && ekin >= limCHIPS1) ||
     (Z == 50 && ekin >= limCHIPS2) ||
     (Z == 64 && ekin >= limCHIPS2)
     ) {
    fXS = ggXsection->ComputeElementXSection(ekin, Z);
    return fXS;
  }
  const G4double emax = pv->GetMaxEnergy();
 
  // low energy based on data
  if(ekin <= emax) {
    fXS = pv->Value(ekin);
    // high energy CHIPS parameterisation
  } else if(ekin >= eTransitionBound) {
    fXS = ggXsection->ComputeElementXSection(ekin, Z);
    // linear interpolation
  } else {
    const G4double rxs = xs150[Z];
    const G4double lxs = pv->Value(emax);
    fXS = lxs + (ekin - emax)*(rxs - lxs)/(eTransitionBound - emax);
  }
 
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout  << "Z= " << Z << " Ekin(MeV)= " << ekin/CLHEP::MeV 
        << ",  nElmXS(b)= " << fXS/CLHEP::barn 
        << G4endl;
  }
#endif
  return fXS;
}

Referenced by GetElementCrossSection(), and IsoCrossSection().

GetElementCrossSection()

G4double G4GammaNuclearXS::GetElementCrossSection ( const G4DynamicParticle * aParticle, G4int Z, const G4Material * mat = nullptr )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 127 of file G4GammaNuclearXS.cc.

{
  return ElementCrossSection(aParticle->GetKineticEnergy(), Z); 
}

GetIsoCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 186 of file G4GammaNuclearXS.cc.

{
  return IsoCrossSection(aParticle->GetKineticEnergy(), Z, A);
}

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 113 of file G4GammaNuclearXS.cc.

{
  return true;
}

IsIsoApplicable()

G4bool G4GammaNuclearXS::IsIsoApplicable ( const G4DynamicParticle * , G4int Z, G4int A, const G4Element * , const G4Material * mat )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 119 of file G4GammaNuclearXS.cc.

{
  return true;
}

IsoCrossSection()

G4double G4GammaNuclearXS::IsoCrossSection	(	const G4double	ekin,
		const G4int	Z,
		const G4int	A )

Definition at line 195 of file G4GammaNuclearXS.cc.

{
  const G4int Z = (ZZ < MAXZGAMMAXS) ? ZZ : MAXZGAMMAXS - 1;
  // cross section per element
  G4double xs = ElementCrossSection(ekin, Z);
 
  if (Z > 2) {
    xs *= A/aeff[Z];
  } else {
    G4int AA = A - amin[Z];
    if(ekin >= ehigh && AA >=0 && AA <=2) { 
      xs *= coeff[Z][AA];
    } else {
      xs = ggXsection->ComputeIsoXSection(ekin, Z, A);
    }
  }
 
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout  << "G4GammaNuclearXS::IsoXS: Z= " << Z << " A= " << A 
        << " Ekin(MeV)= " << ekin/CLHEP::MeV 
        << ", ElmXS(b)= " << xs/CLHEP::barn << G4endl;
  }
#endif
  return xs;
}

Referenced by GetIsoCrossSection(), and SelectIsotope().

LowEnergyCrossSection()

G4double G4GammaNuclearXS::LowEnergyCrossSection	(	G4double	ekin,
		G4int	Z )

Definition at line 179 of file G4GammaNuclearXS.cc.

{    
  const G4int Z = (ZZ < MAXZGAMMAXS) ? ZZ : MAXZGAMMAXS - 1;
  auto pv = data->GetElementData(Z);
  return pv->Value(ekin);
}

operator=()

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

delete

SelectIsotope()

const G4Isotope * G4GammaNuclearXS::SelectIsotope ( const G4Element * anElement, G4double kinEnergy, G4double logE )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 222 of file G4GammaNuclearXS.cc.

{
  std::size_t nIso = anElement->GetNumberOfIsotopes();
  const G4Isotope* iso = anElement->GetIsotope(0);
 
  if(1 == nIso) { return iso; }
 
  const G4double* abundVector = anElement->GetRelativeAbundanceVector();
  G4double sum = 0.0;
  G4int Z = anElement->GetZasInt();
 
  // use isotope cross sections
  std::size_t nn = temp.size();
  if(nn < nIso) { temp.resize(nIso, 0.); }
  
  for (std::size_t j=0; j<nIso; ++j) {
    //G4cout << j << "-th isotope " << (*isoVector)[j]->GetN() 
    //       <<  " abund= " << abundVector[j] << G4endl;
    sum += abundVector[j]*
      IsoCrossSection(kinEnergy, Z, anElement->GetIsotope((G4int)j)->GetN());
    temp[j] = sum;
  }
  sum *= G4UniformRand();
  for (std::size_t j = 0; j<nIso; ++j) {
    if(temp[j] >= sum) {
      iso = anElement->GetIsotope((G4int)j);
      break;
    }
  }
  return iso;
}

---

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

• G4GammaNuclearXS.hh
• G4GammaNuclearXS.cc