G4CrossSectionHP Class Reference

#include <G4CrossSectionHP.hh>

Inheritance diagram for G4CrossSectionHP:

Public Member Functions
	G4CrossSectionHP (const G4ParticleDefinition *, const G4String &nameData, const G4String &nameDir, G4double emaxHP, G4int zmin, G4int zmax)
	~G4CrossSectionHP () override=default
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *, const G4Material *) override
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *, const G4Element *, const G4Material *) override
G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, G4int Z, G4int A, const G4Isotope *iso, const G4Element *elm, const G4Material *mat) override
const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE) override
void	BuildPhysicsTable (const G4ParticleDefinition &) override
void	DumpPhysicsTable (const G4ParticleDefinition &) override
G4CrossSectionHP &	operator= (const G4CrossSectionHP &right)=delete
	G4CrossSectionHP (const G4CrossSectionHP &)=delete
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, 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	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual void	CrossSectionDescription (std::ostream &) 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
void	SetName (const G4String &nam)
G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete
	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Protected Member Functions
G4double	GetMaxHPEnergy () const
void	SetBinSearch (G4int n)

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

Detailed Description

Definition at line 50 of file G4CrossSectionHP.hh.

Constructor & Destructor Documentation

G4CrossSectionHP() [1/2]

G4CrossSectionHP::G4CrossSectionHP ( const G4ParticleDefinition * p, const G4String & nameData, const G4String & nameDir, G4double emaxHP, G4int zmin, G4int zmax )

explicit

Definition at line 61 of file G4CrossSectionHP.cc.

  : G4VCrossSectionDataSet(nameData),
    fParticle(p),
    fNeutron(G4Neutron::Neutron()),
    fManagerHP(G4ParticleHPManager::GetInstance()),
    emax(emaxHP),
    emaxT(fManagerHP->GetMaxEnergyDoppler()),
    elimit(1.0e-11*CLHEP::eV),
    logElimit(G4Log(elimit)),
    minZ(zmin),
    maxZ(zmax),
    fDataName(nameData),
    fDataDirectory(nameDir)
{
  //  verboseLevel = 1;
  if (verboseLevel > 1) {
    G4cout << "G4CrossSectionHP::G4CrossSectionHP: Initialise for "
       << fDataName << "  " << minZ << " < Z < " << maxZ
       << "  EmaxT(MeV)=" << emaxT << G4endl;
    G4cout << "Data directory: " << fDataDirectory << G4endl; 
  }
  auto ptr = G4ElementDataRegistry::Instance();
  auto data = ptr->GetElementDataByName(fDataName);
  if (nullptr == data) { 
    data = new G4ElementData(maxZ - minZ + 1);
    data->SetName(fDataName);
  }
  fData = data;
}

~G4CrossSectionHP()

G4CrossSectionHP::~G4CrossSectionHP ( )

overridedefault

G4CrossSectionHP() [2/2]

G4CrossSectionHP::G4CrossSectionHP ( const G4CrossSectionHP & )

delete

Member Function Documentation

BuildPhysicsTable()

void G4CrossSectionHP::BuildPhysicsTable ( const G4ParticleDefinition & p )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 265 of file G4CrossSectionHP.cc.

{
  if (verboseLevel > 1){
    G4cout << "G4CrossSectionHP::BuildPhysicsTable for "
       << p.GetParticleName() << " and " << fDataName << G4endl;
  }
 
  // it is possible re-initialisation for the second run
  const G4ElementTable* table = G4Element::GetElementTable();
 
  // Access to elements
  for ( auto const & elm : *table ) {
    G4int Z = elm->GetZasInt();
    if (Z >= minZ && Z <= maxZ &&
        nullptr == fData->GetElementData(Z - minZ)) { 
      InitialiseOnFly(Z);
    }
  }
 
  // prepare isotope selection
  std::size_t nmax = 0;
  std::size_t imax = 0;
  for ( auto const & mat : *G4Material::GetMaterialTable() ) {
    std::size_t n = 0; 
    for ( auto const & elm : *mat->GetElementVector() ) {
      std::size_t niso = elm->GetNumberOfIsotopes();
      n += niso;
      if(niso > imax) { imax = niso; }
    }
    if (n > nmax) { nmax = n; }
  }
  fTemp.resize(imax, 0.0);
  fZA.clear();
  fZA.reserve(nmax);
  fIsoXS.resize(nmax, 0.0);
}

ComputeIsoCrossSection()

G4double G4CrossSectionHP::ComputeIsoCrossSection ( G4double kinEnergy, G4double loge, const G4ParticleDefinition * , G4int Z, G4int A, const G4Isotope * iso, const G4Element * elm, const G4Material * mat )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 120 of file G4CrossSectionHP.cc.

{
  G4double ekin = e;
  G4double loge = le;
  if (ekin < elimit) {
    ekin = elimit;
    loge = logElimit;
  }
  if (mat != fCurrentMat) { PrepareCache(mat); }
 
  return IsoCrossSection(ekin, loge, Z, A, mat->GetTemperature()); 
}

DumpPhysicsTable()

void G4CrossSectionHP::DumpPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 302 of file G4CrossSectionHP.cc.

{
  if (fManagerHP->GetVerboseLevel() == 0 || fPrinted)
    return;
 
  //
  // Dump element based cross section
  // range 10e-5 eV to 20 MeV
  // 10 point per decade
  // in barn
  //
  fPrinted = true;
  G4cout << G4endl;
  G4cout << "HP Cross Section " << fDataName << " for "
     << fParticle->GetParticleName() << G4endl;
  G4cout << "(Pointwise cross-section at 0 Kelvin.)" << G4endl;
  G4cout << G4endl;
  G4cout << "Name of Element" << G4endl;
  G4cout << "Energy[eV]  XS[barn]" << G4endl;
  G4cout << G4endl;
 
  const G4ElementTable* table = G4Element::GetElementTable();
  for ( auto const & elm : *table ) {
    G4int Z = elm->GetZasInt();
    if (Z >= minZ && Z <= maxZ && nullptr != fData->GetElementData(Z - minZ)) {
      G4cout << "---------------------------------------------------" << G4endl;
      G4cout << elm->GetName() << G4endl;
      std::size_t n = fData->GetNumberOfComponents(Z);
      for (size_t i=0; i<n; ++i) {
    G4cout << "##  Z=" << Z << "  A=" << fData->GetComponentID(Z - minZ, i);
    G4cout << *(fData->GetComponentDataByIndex(Z - minZ, i)) << G4endl;
      }
    }
  }
}

GetIsoCrossSection()

G4double G4CrossSectionHP::GetIsoCrossSection ( const G4DynamicParticle * dp, G4int Z, G4int A, const G4Isotope * , const G4Element * , const G4Material * mat )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 102 of file G4CrossSectionHP.cc.

{
  G4double ekin = dp->GetKineticEnergy();
  G4double loge = dp->GetLogKineticEnergy();
  if (ekin < elimit) {
    ekin = elimit;
    loge = logElimit;
  }
  if (mat != fCurrentMat) { PrepareCache(mat); }
 
  return IsoCrossSection(ekin, loge, Z, A, mat->GetTemperature());
}

GetMaxHPEnergy()

G4double G4CrossSectionHP::GetMaxHPEnergy ( ) const

inlineprotected

Definition at line 160 of file G4CrossSectionHP.hh.

{ 
  return emax;
}

Referenced by G4NeutronHPCaptureXS::IsElementApplicable(), and G4NeutronHPFissionXS::IsElementApplicable().

IsIsoApplicable()

G4bool G4CrossSectionHP::IsIsoApplicable ( const G4DynamicParticle * dp, G4int Z, G4int A, const G4Element * , const G4Material *  )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 94 of file G4CrossSectionHP.cc.

{
  return (dp->GetKineticEnergy() <= emax);
}

operator=()

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

delete

SelectIsotope()

const G4Isotope * G4CrossSectionHP::SelectIsotope ( const G4Element * elm, G4double kinEnergy, G4double logE )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 214 of file G4CrossSectionHP.cc.

{
  std::size_t nIso = elm->GetNumberOfIsotopes();
  const G4Isotope* iso = elm->GetIsotope(0);
 
  //G4cout << "SelectIsotope NIso= " << nIso << G4endl;
  if(1 == nIso) { return iso; }
 
  // more than 1 isotope
  G4int Z = elm->GetZasInt();
  if (Z >= minZ && Z <= maxZ && nullptr == fData->GetElementData(Z - minZ)) {
    InitialiseOnFly(Z);
  }
  
  const G4double* abundVector = elm->GetRelativeAbundanceVector();
  G4double q = G4UniformRand();
  G4double sum = 0.0;
 
  // is there cached isotope wise cross section?
  std::size_t j;
  if (Z < minZ || Z > maxZ || !CheckCache(Z) ||
      0 == fData->GetNumberOfComponents(Z - minZ)) {
    for (j = 0; j<nIso; ++j) {
      sum += abundVector[j];
      if(q <= sum) {
    iso = elm->GetIsotope((G4int)j);
    break;
      }
    }
    return iso;
  } 
  std::size_t nn = fTemp.size();
  if (nn < nIso) { fTemp.resize(nIso, 0.); }
 
  // reuse cache
  for (j=0; j<nIso; ++j) {
    sum += abundVector[j]*
      GetCrossSection(Z - minZ, elm->GetIsotope((G4int)j)->GetN());
    fTemp[j] = sum;
  }
  sum *= q;
  for (j = 0; j<nIso; ++j) {
    if (fTemp[j] >= sum) {
      iso = elm->GetIsotope((G4int)j);
      break;
    }
  }
  return iso;
}

SetBinSearch()

void G4CrossSectionHP::SetBinSearch ( G4int n )

inlineprotected

Definition at line 165 of file G4CrossSectionHP.hh.

{
  if (n > 0) { binSearch = n; }
}

---

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

• G4CrossSectionHP.hh
• G4CrossSectionHP.cc