G4ParticleInelasticXS Class Reference

#include <G4ParticleInelasticXS.hh>

Inheritance diagram for G4ParticleInelasticXS:

Public Member Functions
	G4ParticleInelasticXS (const G4ParticleDefinition *)
	~G4ParticleInelasticXS () override=default
G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr) final
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr) final
G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr) final
G4double	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *) final
G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, G4int Z, G4int A, const G4Isotope *iso, const G4Element *elm, const G4Material *mat) 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
void	CrossSectionDescription (std::ostream &) const final
G4double	ElementCrossSection (G4double kinEnergy, G4double loge, G4int Z)
G4double	IsoCrossSection (G4double ekin, G4double logE, G4int Z, G4int A)
G4ParticleInelasticXS &	operator= (const G4ParticleInelasticXS &right)=delete
	G4ParticleInelasticXS (const G4ParticleInelasticXS &)=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 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

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

Detailed Description

Definition at line 56 of file G4ParticleInelasticXS.hh.

Constructor & Destructor Documentation

G4ParticleInelasticXS() [1/2]

G4ParticleInelasticXS::G4ParticleInelasticXS ( const G4ParticleDefinition * part )

explicit

Definition at line 69 of file G4ParticleInelasticXS.cc.

  : G4VCrossSectionDataSet("G4ParticleInelasticXS"),
    particle(part),
    elimit(20*CLHEP::MeV)
{
  auto xsr = G4CrossSectionDataSetRegistry::Instance();
  if (nullptr == part) {
    G4Exception("G4ParticleInelasticXS::G4ParticleInelasticXS(..)","had015",
        FatalException, "NO particle definition in constructor");
  } else {
    verboseLevel = 0;
    const G4String& particleName = particle->GetParticleName();
    if(verboseLevel > 1) {
      G4cout << "G4ParticleInelasticXS::G4ParticleInelasticXS for " 
         << particleName << " on atoms with Z < " << MAXZINELP << G4endl;
    }
    if (particleName == "proton") {
      highEnergyXsection = xsr->GetComponentCrossSection("Glauber-Gribov");
      if(highEnergyXsection == nullptr) {
    highEnergyXsection = new G4ComponentGGHadronNucleusXsc();
      }
    } else {
      highEnergyXsection = 
        xsr->GetComponentCrossSection("Glauber-Gribov Nucl-nucl");
      if(highEnergyXsection == nullptr) {
    highEnergyXsection = new G4ComponentGGNuclNuclXsc();
      }
      for (index=1; index<5; ++index) {
        if (particleName == pname[index]) { break; }
      }
      index = std::min(index, 4);
      if (1 < index) { SetMaxKinEnergy(25.6*CLHEP::PeV); }
    }
  }
  // this should never happens but ...
  if (nullptr == highEnergyXsection) {
    highEnergyXsection = xsr->GetComponentCrossSection("Glauber-Gribov");
  }
  SetForAllAtomsAndEnergies(true);
  if (gDataDirectory.empty()) {
    gDataDirectory = G4HadronicParameters::Instance()->GetDirPARTICLEXS();
  }
  G4String ss = pname[index] + "ParticleXS";
  SetName(ss);
  if (data[index] == nullptr) {
    data[index] = new G4ElementData(MAXZINELP);
    data[index]->SetName(pname[index] + "PartInel");
    for (G4int Z=1; Z<MAXZINELP; ++Z) { Initialise(Z); }
  }
}

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

~G4ParticleInelasticXS()

G4ParticleInelasticXS::~G4ParticleInelasticXS ( )

overridedefault

G4ParticleInelasticXS() [2/2]

G4ParticleInelasticXS::G4ParticleInelasticXS ( const G4ParticleInelasticXS & )

delete

Member Function Documentation

BuildPhysicsTable()

void G4ParticleInelasticXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 301 of file G4ParticleInelasticXS.cc.

{
  if (verboseLevel > 0){
    G4cout << "G4ParticleInelasticXS::BuildPhysicsTable for " 
       << p.GetParticleName() << G4endl;
  }
  if (&p != particle) { 
    G4ExceptionDescription ed;
    ed << p.GetParticleName() << " is a wrong particle type -"
       << particle->GetParticleName() << " is expected";
    G4Exception("G4ParticleInelasticXS::BuildPhysicsTable(..)","had012",
        FatalException, ed, "");
    return; 
  }
 
  // it is possible re-initialisation for the new run
  const G4ElementTable* table = G4Element::GetElementTable();
 
  // prepare isotope selection
  std::size_t nIso = temp.size();
 
  // Access to elements
  for ( auto const & elm : *table ) {
    std::size_t n = elm->GetNumberOfIsotopes();
    if (n > nIso) { nIso = n; }
    G4int Z = std::min( elm->GetZasInt(), MAXZINELP-1);
    if ( nullptr == (data[index])->GetElementData(Z) ) {
      InitialiseOnFly(Z);
    }
  }
  temp.resize(nIso, 0.0);
}

Referenced by G4InterfaceToXS::G4InterfaceToXS(), and ~G4ParticleInelasticXS().

ComputeCrossSectionPerElement()

G4double G4ParticleInelasticXS::ComputeCrossSectionPerElement ( G4double kinEnergy, G4double loge, const G4ParticleDefinition * , const G4Element * elm, const G4Material *  )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 153 of file G4ParticleInelasticXS.cc.

{
  return ElementCrossSection(ekin, loge, elm->GetZasInt());
}

Referenced by ~G4ParticleInelasticXS().

ComputeIsoCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 188 of file G4ParticleInelasticXS.cc.

{
  return IsoCrossSection(ekin, loge, Z, A); 
}

Referenced by ~G4ParticleInelasticXS().

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 120 of file G4ParticleInelasticXS.cc.

{
  outFile << "G4ParticleInelasticXS calculates n, p, d, t, he3, he4 inelastic\n"
          << "cross section on nuclei using data from the high precision\n"
          << "neutron database.  These data are simplified and smoothed over\n"
          << "the resonance region in order to reduce CPU time.\n"
          << "For high energy Glauber-Gribov cross section model is used.\n";
}

Referenced by ~G4ParticleInelasticXS().

ElementCrossSection()

G4double G4ParticleInelasticXS::ElementCrossSection	(	G4double	kinEnergy,
		G4double	loge,
		G4int	Z )

Definition at line 161 of file G4ParticleInelasticXS.cc.

{
  G4int Z = (ZZ >= MAXZINELP) ? MAXZINELP - 1 : ZZ; 
 
  // element data is always valid pointer by construction of XS
  auto pv = GetPhysicsVector(Z);
 
  // set to null x-section below lowest energy in the table
  G4double xs = 0.0;
  if (ekin > pv->Energy(0)) {
    xs = (ekin <= pv->GetMaxEnergy()) ? pv->LogVectorValue(ekin, loge) 
    : coeff[Z][index]*highEnergyXsection->GetInelasticElementCrossSection(particle,
                    ekin, Z, aeff[Z]);
  }
 
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout  << "ElmXS: Z= " << Z << " Ekin(MeV)= " << ekin/CLHEP::MeV 
        << " xs(bn)= " << xs/CLHEP::barn << " element data for "
        << particle->GetParticleName()
            << " idx= " << index << G4endl;
  }
#endif
  return xs;
}

Referenced by ComputeCrossSectionPerElement(), GetElementCrossSection(), and ~G4ParticleInelasticXS().

GetElementCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 145 of file G4ParticleInelasticXS.cc.

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

Referenced by ~G4ParticleInelasticXS().

GetIsoCrossSection()

G4double G4ParticleInelasticXS::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 197 of file G4ParticleInelasticXS.cc.

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

Referenced by ~G4ParticleInelasticXS().

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 130 of file G4ParticleInelasticXS.cc.

{
  return true;
}

IsIsoApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 137 of file G4ParticleInelasticXS.cc.

{
  return true;
}

Referenced by ~G4ParticleInelasticXS().

IsoCrossSection()

G4double G4ParticleInelasticXS::IsoCrossSection	(	G4double	ekin,
		G4double	logE,
		G4int	Z,
		G4int	A )

Definition at line 206 of file G4ParticleInelasticXS.cc.

{
  G4double xs = 0.0;
  G4int Z = (ZZ >= MAXZINELP) ? MAXZINELP - 1 : ZZ; 
 
  // needed here to gurantee upload data for Z
  auto pv = GetPhysicsVector(Z);
 
  // compute isotope cross section if applicable 
  if (ekin <= elimit && data[index]->GetNumberOfComponents(Z) > 0) {
    auto pviso = data[index]->GetComponentDataByID(Z, A);
    if (pviso != nullptr && ekin > pviso->Energy(0)) { 
      xs = pviso->LogVectorValue(ekin, logE);
#ifdef G4VERBOSE
      if(verboseLevel > 1) {
    G4cout << "G4ParticleInelasticXS::IsoXS: for " 
               << particle->GetParticleName() << " Ekin(MeV)= " 
               << ekin/CLHEP::MeV << "  xs(b)= " << xs/CLHEP::barn 
           << "  Z= " << Z << "  A= " << A 
               << " idx= " << index << G4endl;
      }
#endif
      return xs;
    }
  }
  // use element x-section
  if (ekin > pv->Energy(0)) {
    xs = (ekin <= pv->GetMaxEnergy()) ? pv->LogVectorValue(ekin, logE) :
      coeff[Z][index] *
      highEnergyXsection->GetInelasticElementCrossSection(particle, ekin, Z, aeff[Z])
      * A/aeff[Z];
  }
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout  << "IsoXS for " << particle->GetParticleName() 
        << " Target Z= " << Z << " A= " << A 
        << " Ekin(MeV)= " << ekin/CLHEP::MeV 
        << " xs(bn)= " << xs/CLHEP::barn
            << " idx= " << index << G4endl;
  }
#endif
  return xs;
}

Referenced by ComputeIsoCrossSection(), GetIsoCrossSection(), SelectIsotope(), and ~G4ParticleInelasticXS().

operator=()

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

delete

SelectIsotope()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 251 of file G4ParticleInelasticXS.cc.

{
  G4int nIso = (G4int)anElement->GetNumberOfIsotopes();
  const G4Isotope* iso = anElement->GetIsotope(0);
 
  if (1 == nIso) { return iso; }
 
  // more than 1 isotope
  G4int Z = anElement->GetZasInt();
 
  // initialisation for given Z
  GetPhysicsVector(Z);
 
  const G4double* abundVector = anElement->GetRelativeAbundanceVector();
  G4double q = G4UniformRand();
  G4double sum = 0.0;
  G4int j;
 
  // isotope wise cross section not available
  if (Z >= MAXZINELP || 0 == data[index]->GetNumberOfComponents(Z)) {
    for (j=0; j<nIso; ++j) {
      sum += abundVector[j];
      if(q <= sum) {
    iso = anElement->GetIsotope(j);
    break;
      }
    }
    return iso;
  }
 
  G4int nn = (G4int)temp.size();
  if (nn < nIso) { temp.resize(nIso, 0.); }
 
  for (j=0; j<nIso; ++j) {
    sum += abundVector[j]*IsoCrossSection(kinEnergy, logE, Z, 
                      anElement->GetIsotope(j)->GetN());
    temp[j] = sum;
  }
  sum *= q;
  for (j=0; j<nIso; ++j) {
    if (temp[j] >= sum) {
      iso = anElement->GetIsotope(j);
      break;
    }
  }
  return iso;
}

Referenced by ~G4ParticleInelasticXS().

---

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

• G4ParticleInelasticXS.hh
• G4ParticleInelasticXS.cc