G4VCrossSectionDataSet Class Reference

#include <G4VCrossSectionDataSet.hh>

Inherited by CexmcFakeCrossSectionData, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4BGGPionInelasticXS, G4ChipsAntiBaryonElasticXS, G4ChipsAntiBaryonInelasticXS, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonMinusElasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsKaonPlusInelasticXS, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4ChipsNeutronElasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4ChipsPionPlusElasticXS, G4ChipsPionPlusInelasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4CrossSectionElastic, G4CrossSectionInelastic, G4CrossSectionPairGG, G4ElectroNuclearCrossSection, G4EMDissociationCrossSection, G4GeneralSpaceNNCrossSection, G4HadronCaptureDataSet, G4HadronElasticDataSet, G4HadronFissionDataSet, G4HadronInelasticDataSet, G4IonProtonCrossSection, G4IonsKoxCrossSection, G4IonsShenCrossSection, G4IonsSihverCrossSection, G4KokoulinMuonNuclearXS, G4LENDCrossSection, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticCrossSection, G4NeutronInelasticXS, G4NucleonNuclearCrossSection, G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, G4ParticleHPInelasticData, G4ParticleHPJENDLHEData, G4ParticleHPThermalScatteringData, G4PhotoNuclearCrossSection, G4PiNuclearCrossSection, G4ProtonInelasticCrossSection, G4TripathiCrossSection, G4TripathiLightCrossSection, G4UPiNuclearCrossSection, and G4ZeroXS.

Collaboration diagram for G4VCrossSectionDataSet:

Public Member Functions
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)
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)
const G4String &	GetName () const

Protected Member Functions
void	SetName (const G4String &)

Protected Attributes
G4int	verboseLevel

Private Member Functions
G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)
	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)

Private Attributes
G4CrossSectionDataSetRegistry *	registry
G4double	minKinEnergy
G4double	maxKinEnergy
G4String	name

Detailed Description

Definition at line 71 of file G4VCrossSectionDataSet.hh.

Constructor & Destructor Documentation

G4VCrossSectionDataSet() [1/2]

G4VCrossSectionDataSet::G4VCrossSectionDataSet

(

const G4String &

nam = ""

)

Definition at line 54 of file G4VCrossSectionDataSet.cc.

                                                                  :
  verboseLevel(0),minKinEnergy(0.0),maxKinEnergy(100*TeV),name(nam) 
{
  registry = G4CrossSectionDataSetRegistry::Instance();
  registry->Register(this);
}

Here is the call graph for this function:

~G4VCrossSectionDataSet()

G4VCrossSectionDataSet::~G4VCrossSectionDataSet ( )

virtual

Definition at line 61 of file G4VCrossSectionDataSet.cc.

{
  registry->DeRegister(this);
}

Here is the call graph for this function:

G4VCrossSectionDataSet() [2/2]

G4VCrossSectionDataSet::G4VCrossSectionDataSet ( const G4VCrossSectionDataSet &  )

private

Member Function Documentation

BuildPhysicsTable()

void G4VCrossSectionDataSet::BuildPhysicsTable ( const G4ParticleDefinition &  )

virtual

Reimplemented in G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4NeutronInelasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronCaptureXS, G4ParticleHPThermalScatteringData, G4LENDCrossSection, G4CrossSectionElastic, G4CrossSectionInelastic, G4ParticleHPInelasticData, G4NeutronElasticXS, G4ParticleHPElasticData, G4ParticleHPCaptureData, G4CrossSectionPairGG, G4ParticleHPFissionData, G4KokoulinMuonNuclearXS, CexmcFakeCrossSectionData, G4UPiNuclearCrossSection, G4ParticleHPJENDLHEData, and G4PiNuclearCrossSection.

Definition at line 193 of file G4VCrossSectionDataSet.cc.

{}

Here is the caller graph for this function:

ComputeCrossSection()

G4double G4VCrossSectionDataSet::ComputeCrossSection	(	const G4DynamicParticle *	part,
		const G4Element *	elm,
		const G4Material *	mat = 0
	)

Definition at line 84 of file G4VCrossSectionDataSet.cc.

{
  G4int Z = G4lrint(elm->GetZ());

  if (IsElementApplicable(part, Z, mat)) { 
    return GetElementCrossSection(part, Z, mat);
  }

  // isotope-wise cross section making sum over available
  // isotope cross sections, which may be incomplete, so
  // the result is corrected 
  G4int nIso = elm->GetNumberOfIsotopes();    
  G4double fact = 0.0;
  G4double xsec = 0.0;
  G4Isotope* iso = 0;

  if (0 < nIso) { 

    // user-defined isotope abundances        
    G4IsotopeVector* isoVector = elm->GetIsotopeVector();
    G4double* abundVector = elm->GetRelativeAbundanceVector();

    for (G4int j = 0; j<nIso; ++j) {
      iso = (*isoVector)[j];
      G4int A = iso->GetN();
      if(abundVector[j] > 0.0 && IsIsoApplicable(part, Z, A, elm, mat)) {
        fact += abundVector[j];
    xsec += abundVector[j]*GetIsoCrossSection(part, Z, A, iso, elm, mat);
      }
    }

  } else {

    // natural isotope abundances
    G4NistManager* nist = G4NistManager::Instance();
    G4int n0 = nist->GetNistFirstIsotopeN(Z);
    G4int nn = nist->GetNumberOfNistIsotopes(Z);
    for (G4int A = n0; A < n0+nn; ++A) {
      G4double abund = nist->GetIsotopeAbundance(Z, A);
      if(abund > 0.0 && IsIsoApplicable(part, Z, A, elm, mat)) {
        fact += abund;
    xsec += abund*GetIsoCrossSection(part, Z, A, iso, elm, mat);
      }
    }
  }
  if(fact > 0.0) { xsec /= fact; }
  return xsec;
}

Here is the call graph for this function:

Here is the caller graph for this function:

CrossSectionDescription()

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

virtual

Reimplemented in G4GeneralSpaceNNCrossSection, G4ParticleHPThermalScatteringData, G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4NeutronInelasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronCaptureXS, G4ParticleHPInelasticData, G4ParticleHPCaptureData, G4ParticleHPElasticData, G4CrossSectionElastic, G4CrossSectionInelastic, G4NeutronElasticXS, G4ParticleHPFissionData, G4IonProtonCrossSection, G4IonsShenCrossSection, G4ZeroXS, G4UPiNuclearCrossSection, G4ElectroNuclearCrossSection, G4CrossSectionPairGG, G4NucleonNuclearCrossSection, G4IonsKoxCrossSection, G4IonsSihverCrossSection, G4KokoulinMuonNuclearXS, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4NeutronInelasticCrossSection, G4ChipsPionPlusInelasticXS, G4ChipsKaonPlusInelasticXS, G4ChipsAntiBaryonElasticXS, G4ChipsKaonMinusElasticXS, G4ChipsNeutronElasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4PiNuclearCrossSection, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsPionPlusElasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4HadronCaptureDataSet, G4ChipsAntiBaryonInelasticXS, G4HadronElasticDataSet, G4HadronFissionDataSet, G4HadronInelasticDataSet, and G4PhotoNuclearCrossSection.

Definition at line 199 of file G4VCrossSectionDataSet.cc.

{
  outFile << "The description for this cross section data set has not been written yet.\n";
}

Here is the caller graph for this function:

DumpPhysicsTable()

void G4VCrossSectionDataSet::DumpPhysicsTable ( const G4ParticleDefinition &  )

virtual

Reimplemented in G4ParticleHPThermalScatteringData, G4LENDCrossSection, G4CrossSectionElastic, G4CrossSectionInelastic, G4ParticleHPInelasticData, G4ParticleHPElasticData, G4ParticleHPCaptureData, G4CrossSectionPairGG, G4ParticleHPFissionData, CexmcFakeCrossSectionData, G4UPiNuclearCrossSection, and G4ParticleHPJENDLHEData.

Definition at line 196 of file G4VCrossSectionDataSet.cc.

{}

GetCrossSection()

G4double G4VCrossSectionDataSet::GetCrossSection ( const G4DynamicParticle *  dp, const G4Element *  elm, const G4Material *  mat = 0  )

inline

Definition at line 183 of file G4VCrossSectionDataSet.hh.

{
  return ComputeCrossSection(dp, elm, mat);
}

Here is the call graph for this function:

GetElementCrossSection()

G4double G4VCrossSectionDataSet::GetElementCrossSection ( const G4DynamicParticle *  dynPart, G4int  Z, const G4Material *  mat = 0  )

virtual

Reimplemented in G4GeneralSpaceNNCrossSection, G4TripathiLightCrossSection, G4EMDissociationCrossSection, G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4CrossSectionElastic, G4CrossSectionInelastic, G4NeutronElasticXS, G4CrossSectionPairGG, G4IonProtonCrossSection, G4ElectroNuclearCrossSection, G4ZeroXS, G4KokoulinMuonNuclearXS, G4ProtonInelasticCrossSection, G4IonsShenCrossSection, G4HadronCaptureDataSet, G4NucleonNuclearCrossSection, G4HadronElasticDataSet, G4IonsKoxCrossSection, G4IonsSihverCrossSection, G4HadronFissionDataSet, G4HadronInelasticDataSet, G4PhotoNuclearCrossSection, G4NeutronInelasticCrossSection, G4TripathiCrossSection, and G4PiNuclearCrossSection.

Definition at line 136 of file G4VCrossSectionDataSet.cc.

{
  G4cout << "G4VCrossSectionDataSet::GetCrossSection per element ERROR: "
     << " there is no cross section for "
     << dynPart->GetDefinition()->GetParticleName()
     << "  E(MeV)= "  << dynPart->GetKineticEnergy()/MeV;
  if(mat) { G4cout << "  inside " << mat->GetName(); }
  G4cout << " for Z= " << Z << G4endl;
  throw G4HadronicException(__FILE__, __LINE__,
        "G4VCrossSectionDataSet::GetElementCrossSection is absent");
  return 0.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetIsoCrossSection()

G4double G4VCrossSectionDataSet::GetIsoCrossSection ( const G4DynamicParticle *  dynPart, G4int  Z, G4int  A, const G4Isotope *  iso = 0, const G4Element *  elm = 0, const G4Material *  mat = 0  )

virtual

Reimplemented in G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4LENDCrossSection, G4ParticleHPThermalScatteringData, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4IonsShenCrossSection, G4ChipsNeutronInelasticXS, G4ChipsPionPlusInelasticXS, G4ParticleHPInelasticData, G4ChipsKaonPlusInelasticXS, G4ParticleHPElasticData, G4ChipsAntiBaryonElasticXS, G4ChipsKaonMinusElasticXS, G4ChipsNeutronElasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4ParticleHPCaptureData, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsPionPlusElasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4ChipsAntiBaryonInelasticXS, and G4ParticleHPFissionData.

Definition at line 152 of file G4VCrossSectionDataSet.cc.

{
  G4cout << "G4VCrossSectionDataSet::GetCrossSection per isotope ERROR: "
     << " there is no cross section for "
     << dynPart->GetDefinition()->GetParticleName()
     << "  E(MeV)= "  << dynPart->GetKineticEnergy()/MeV;
  if(mat) { G4cout << "  inside " << mat->GetName(); }
  if(elm) { G4cout << " for " << elm->GetName(); }
  G4cout << "  Z= " << Z << " A= " << A << G4endl;
  throw G4HadronicException(__FILE__, __LINE__,
        "G4VCrossSectionDataSet::GetIsoCrossSection is absent");
  return 0.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetMaxKinEnergy()

G4double G4VCrossSectionDataSet::GetMaxKinEnergy ( ) const

inline

Definition at line 216 of file G4VCrossSectionDataSet.hh.

{
  return maxKinEnergy;
}

Here is the caller graph for this function:

GetMinKinEnergy()

G4double G4VCrossSectionDataSet::GetMinKinEnergy ( ) const

inline

Definition at line 206 of file G4VCrossSectionDataSet.hh.

{
  return minKinEnergy;
}

Here is the caller graph for this function:

GetName()

const G4String & G4VCrossSectionDataSet::GetName ( void  ) const

inline

Definition at line 221 of file G4VCrossSectionDataSet.hh.

{
  return name;
}

Here is the caller graph for this function:

GetVerboseLevel()

G4int G4VCrossSectionDataSet::GetVerboseLevel ( ) const

inlinevirtual

Reimplemented in G4ParticleHPInelasticData, G4ParticleHPElasticData, G4ParticleHPCaptureData, and G4ParticleHPFissionData.

Definition at line 191 of file G4VCrossSectionDataSet.hh.

{
  return verboseLevel;
}

IsElementApplicable()

G4bool G4VCrossSectionDataSet::IsElementApplicable ( const G4DynamicParticle *  , G4int  Z, const G4Material *  mat = 0  )

virtual

Reimplemented in G4GeneralSpaceNNCrossSection, G4TripathiLightCrossSection, G4EMDissociationCrossSection, G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4ParticleHPBGGNucleonInelasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4CrossSectionElastic, G4CrossSectionInelastic, G4NeutronCaptureXS, G4NeutronElasticXS, G4CrossSectionPairGG, G4ElectroNuclearCrossSection, G4IonProtonCrossSection, G4KokoulinMuonNuclearXS, G4ZeroXS, G4ProtonInelasticCrossSection, G4IonsShenCrossSection, G4HadronCaptureDataSet, G4NucleonNuclearCrossSection, G4UPiNuclearCrossSection, G4HadronElasticDataSet, G4IonsKoxCrossSection, G4IonsSihverCrossSection, G4HadronFissionDataSet, G4HadronInelasticDataSet, G4PhotoNuclearCrossSection, G4NeutronInelasticCrossSection, G4TripathiCrossSection, and G4PiNuclearCrossSection.

Definition at line 67 of file G4VCrossSectionDataSet.cc.

{
  return false;
}

Here is the caller graph for this function:

IsIsoApplicable()

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

virtual

Reimplemented in G4BGGPionInelasticXS, G4LENDCrossSection, G4BGGNucleonElasticXS, G4ParticleHPBGGNucleonInelasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4ParticleHPThermalScatteringData, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionPlusInelasticXS, G4ParticleHPInelasticData, G4ChipsKaonPlusInelasticXS, G4ParticleHPElasticData, G4ChipsAntiBaryonElasticXS, G4ChipsKaonMinusElasticXS, G4ChipsNeutronElasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4ParticleHPCaptureData, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsPionPlusElasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4ChipsAntiBaryonInelasticXS, and G4ParticleHPFissionData.

Definition at line 75 of file G4VCrossSectionDataSet.cc.

{
  return false;
}

Here is the caller graph for this function:

operator=()

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

private

SelectIsotope()

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

virtual

Reimplemented in G4NeutronInelasticXS, and G4NeutronCaptureXS.

Definition at line 171 of file G4VCrossSectionDataSet.cc.

{
  G4int nIso = anElement->GetNumberOfIsotopes();
  G4IsotopeVector* isoVector = anElement->GetIsotopeVector();
  G4Isotope* iso = (*isoVector)[0];

  // more than 1 isotope
  if(1 < nIso) {
    G4double* abundVector = anElement->GetRelativeAbundanceVector();
    G4double sum = 0.0;
    G4double q = G4UniformRand();
    for (G4int j = 0; j<nIso; ++j) {
      sum += abundVector[j];
      if(q <= sum) {
    iso = (*isoVector)[j];
    break;
      }
    }
  }
  return iso;
}

Here is the call graph for this function:

SetMaxKinEnergy()

void G4VCrossSectionDataSet::SetMaxKinEnergy ( G4double  value )

inline

Definition at line 211 of file G4VCrossSectionDataSet.hh.

{
  maxKinEnergy = value;
}

Here is the caller graph for this function:

SetMinKinEnergy()

void G4VCrossSectionDataSet::SetMinKinEnergy ( G4double  value )

inline

Definition at line 201 of file G4VCrossSectionDataSet.hh.

{
  minKinEnergy = value;
}

Here is the caller graph for this function:

SetName()

void G4VCrossSectionDataSet::SetName ( const G4String &  nam )

inlineprotected

Definition at line 226 of file G4VCrossSectionDataSet.hh.

{
  name = nam;
}

Here is the caller graph for this function:

SetVerboseLevel()

void G4VCrossSectionDataSet::SetVerboseLevel ( G4int  value )

inlinevirtual

Reimplemented in G4ParticleHPInelasticData, G4ParticleHPCaptureData, G4ParticleHPElasticData, and G4ParticleHPFissionData.

Definition at line 196 of file G4VCrossSectionDataSet.hh.

{
  verboseLevel = value;
}

Member Data Documentation

maxKinEnergy

G4double G4VCrossSectionDataSet::maxKinEnergy

private

Definition at line 177 of file G4VCrossSectionDataSet.hh.

minKinEnergy

G4double G4VCrossSectionDataSet::minKinEnergy

private

Definition at line 176 of file G4VCrossSectionDataSet.hh.

name

G4String G4VCrossSectionDataSet::name

private

Definition at line 179 of file G4VCrossSectionDataSet.hh.

registry

G4CrossSectionDataSetRegistry* G4VCrossSectionDataSet::registry

private

Definition at line 174 of file G4VCrossSectionDataSet.hh.

verboseLevel

G4int G4VCrossSectionDataSet::verboseLevel

protected

Definition at line 167 of file G4VCrossSectionDataSet.hh.

---

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

• G4VCrossSectionDataSet.hh
• G4VCrossSectionDataSet.cc