G4GammaNuclearXS Class Reference

#include <G4GammaNuclearXS.hh>

Inheritance diagram for G4GammaNuclearXS:

Public Member Functions
	G4GammaNuclearXS ()
	~G4GammaNuclearXS () final
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) final
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso, const G4Element *elm, const G4Material *mat) final
const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE) final
void	BuildPhysicsTable (const G4ParticleDefinition &) final
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 ()
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 57 of file G4GammaNuclearXS.hh.

Constructor & Destructor Documentation

G4GammaNuclearXS() [1/2]

G4GammaNuclearXS::G4GammaNuclearXS ( )

explicit

Definition at line 67 of file G4GammaNuclearXS.cc.

 : G4VCrossSectionDataSet(Default_Name()),
   ggXsection(nullptr),
   gamma(G4Gamma::Gamma()),
   isMaster(false)
{
  //  verboseLevel = 0;
  if(verboseLevel > 0){
    G4cout  << "G4GammaNuclearXS::G4GammaNuclearXS Initialise for Z < " 
        << MAXZEL << G4endl;
  }
  ggXsection = G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet("PhotoNuclearXS");
  if(ggXsection == nullptr) ggXsection = new G4PhotoNuclearCrossSection();
  SetForAllAtomsAndEnergies(true);
}

~G4GammaNuclearXS()

G4GammaNuclearXS::~G4GammaNuclearXS ( )

final

Definition at line 83 of file G4GammaNuclearXS.cc.

{
  if(isMaster) {
    for(G4int i=0; i<MAXZEL; ++i) {
      delete data[i];
      data[i] = nullptr;
    }
  }
}

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 181 of file G4GammaNuclearXS.cc.

{
  if(verboseLevel > 0){
    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; 
  }
  if(0. == coeff[0]) { 
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&gNuclearXSMutex);
    if(0. == coeff[0]) { 
#endif
      coeff[0] = 1.0;
      isMaster = true;
      FindDirectoryPath();
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&gNuclearXSMutex);
#endif
  }
 
  // it is possible re-initialisation for the second run
  if(isMaster) {
 
    // Access to elements
    auto theCoupleTable = G4ProductionCutsTable::GetProductionCutsTable();
    size_t numOfCouples = theCoupleTable->GetTableSize();
    for(size_t j=0; j<numOfCouples; ++j) {
      auto mat = theCoupleTable->GetMaterialCutsCouple(j)->GetMaterial();
      auto elmVec = mat->GetElementVector();
      size_t numOfElem = mat->GetNumberOfElements();
      for (size_t ie = 0; ie < numOfElem; ++ie) {
    G4int Z = std::max(1,std::min(((*elmVec)[ie])->GetZasInt(), MAXZEL-1));
    if(data[Z] == nullptr) { Initialise(Z); }
      }
    }
  }
}

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 93 of file G4GammaNuclearXS.cc.

{
  outFile << "G4GammaNuclearXS calculates the gamma nuclear\n"
          << "cross section on nuclei using data from the high precision\n"
          << "LEND gamma database. The data are simplified and smoothed over\n"
          << "the resonance region in order to reduce CPU time.\n"
          << "For high energies Glauber-Gribiv cross section is used.\n";
}

Default_Name()

static const char * G4GammaNuclearXS::Default_Name ( )

inlinestatic

Definition at line 65 of file G4GammaNuclearXS.hh.

{return "G4GammaNuclearXS";}

GetElementCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 117 of file G4GammaNuclearXS.cc.

{
  G4double xs = 0.0;
  G4double ekin = aParticle->GetKineticEnergy();
 
  G4int Z = (ZZ >= MAXZEL) ? MAXZEL - 1 : ZZ; 
 
  auto pv = GetPhysicsVector(Z);
  if(pv == nullptr) { return xs; }
  //  G4cout  << "G4GammaNuclearXS::GetCrossSection e= " << ekin 
  // << " Z= " << Z << G4endl;
 
  if(ekin <= pv->GetMaxEnergy()) { 
    xs = pv->LogVectorValue(ekin, aParticle->GetLogKineticEnergy()); 
  } else {          
    xs = coeff[Z]*ggXsection->GetElementCrossSection(aParticle, Z, mat);
  }
 
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout  << "Z= " << Z << " Ekin(MeV)= " << ekin/CLHEP::MeV 
        << ",  nElmXS(b)= " << xs/CLHEP::barn 
        << G4endl;
  }
#endif
  return xs;
}

Referenced by GetIsoCrossSection().

GetIsoCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 146 of file G4GammaNuclearXS.cc.

{
  return GetElementCrossSection(aParticle, Z, mat) * A/aeff[Z];
}

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 103 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 109 of file G4GammaNuclearXS.cc.

{
  return true;
}

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 155 of file G4GammaNuclearXS.cc.

{
  size_t nIso = anElement->GetNumberOfIsotopes();
  const G4Isotope* iso = anElement->GetIsotope(0);
 
  //G4cout << "SelectIsotope NIso= " << nIso << G4endl;
  if(1 == nIso) { return iso; }
 
  const G4double* abundVector = anElement->GetRelativeAbundanceVector();
  G4double q = G4UniformRand();
  G4double sum = 0.0;
  size_t j;
 
  // isotope wise cross section not used
  for (j=0; j<nIso; ++j) {
    sum += abundVector[j];
    if(q <= sum) {
      iso = anElement->GetIsotope(j);
      break;
    }
  }
  return iso;
}

---

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

• G4GammaNuclearXS.hh
• G4GammaNuclearXS.cc