G4VCrossSectionDataSet Class Reference

#include <G4VCrossSectionDataSet.hh>

Inheritance 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

Detailed Description

Definition at line 70 of file G4VCrossSectionDataSet.hh.

Constructor & Destructor Documentation

G4VCrossSectionDataSet::G4VCrossSectionDataSet

(

const G4String &

nam = ""

)

Definition at line 54 of file G4VCrossSectionDataSet.cc.

References G4CrossSectionDataSetRegistry::Instance(), and G4CrossSectionDataSetRegistry::Register().

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

G4VCrossSectionDataSet::~G4VCrossSectionDataSet ( )

virtual

Definition at line 60 of file G4VCrossSectionDataSet.cc.

References G4CrossSectionDataSetRegistry::DeRegister(), and G4CrossSectionDataSetRegistry::Instance().

{
  G4CrossSectionDataSetRegistry::Instance()->DeRegister(this);
}

Member Function Documentation

void G4VCrossSectionDataSet::BuildPhysicsTable ( const G4ParticleDefinition &  )

virtual

Reimplemented in G4DPMJET2_5CrossSection, G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4NeutronHPThermalScatteringData, G4LENDCrossSection, G4CrossSectionElastic, G4CrossSectionInelastic, G4NeutronHPElasticData, G4CrossSectionPairGG, G4NeutronElasticXS, G4NeutronHPCaptureData, G4NeutronHPInelasticData, G4GGNuclNuclCrossSection, G4NeutronHPFissionData, CexmcFakeCrossSectionData, G4UPiNuclearCrossSection, G4NeutronHPJENDLHEData, and G4PiNuclearCrossSection.

Definition at line 192 of file G4VCrossSectionDataSet.cc.

Referenced by G4CrossSectionPairGG::BuildPhysicsTable(), G4BGGNucleonInelasticXS::BuildPhysicsTable(), G4BGGPionElasticXS::BuildPhysicsTable(), G4BGGNucleonElasticXS::BuildPhysicsTable(), and G4BGGPionInelasticXS::BuildPhysicsTable().

{}

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

Definition at line 83 of file G4VCrossSectionDataSet.cc.

References G4lrint(), GetElementCrossSection(), GetIsoCrossSection(), G4NistManager::GetIsotopeAbundance(), G4Element::GetIsotopeVector(), G4Isotope::GetN(), G4NistManager::GetNistFirstIsotopeN(), G4Element::GetNumberOfIsotopes(), G4NistManager::GetNumberOfNistIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), G4NistManager::Instance(), IsElementApplicable(), IsIsoApplicable(), and G4InuclParticleNames::nn.

Referenced by GetCrossSection().

{
  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;
}

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

virtual

Reimplemented in G4BGGPionInelasticXS, G4GlauberGribovCrossSection, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4CrossSectionElastic, G4CrossSectionInelastic, G4GGNuclNuclCrossSection, G4NeutronElasticXS, G4IonProtonCrossSection, G4IonsShenCrossSection, G4UPiNuclearCrossSection, G4ElectroNuclearCrossSection, G4CrossSectionPairGG, G4NucleonNuclearCrossSection, G4IonsKoxCrossSection, G4IonsSihverCrossSection, G4NeutronInelasticCrossSection, G4HadronCaptureDataSet, G4PiNuclearCrossSection, G4HadronElasticDataSet, G4HadronFissionDataSet, G4HadronInelasticDataSet, and G4PhotoNuclearCrossSection.

Definition at line 198 of file G4VCrossSectionDataSet.cc.

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

void G4VCrossSectionDataSet::DumpPhysicsTable ( const G4ParticleDefinition &  )

virtual

Reimplemented in G4DPMJET2_5CrossSection, G4NeutronHPThermalScatteringData, G4LENDCrossSection, G4CrossSectionElastic, G4CrossSectionInelastic, G4GGNuclNuclCrossSection, G4NeutronHPElasticData, G4NeutronHPCaptureData, G4NeutronHPInelasticData, G4CrossSectionPairGG, G4NeutronHPFissionData, CexmcFakeCrossSectionData, G4UPiNuclearCrossSection, and G4NeutronHPJENDLHEData.

Definition at line 195 of file G4VCrossSectionDataSet.cc.

{}

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

inline

Definition at line 180 of file G4VCrossSectionDataSet.hh.

References ComputeCrossSection().

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

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, G4CrossSectionPairGG, G4IonProtonCrossSection, G4NeutronElasticXS, G4ElectroNuclearCrossSection, G4KokoulinMuonNuclearXS, G4GGNuclNuclCrossSection, G4ProtonInelasticCrossSection, G4IonsShenCrossSection, G4HadronCaptureDataSet, G4NucleonNuclearCrossSection, G4HadronElasticDataSet, G4IonsKoxCrossSection, G4IonsSihverCrossSection, G4HadronFissionDataSet, G4HadronInelasticDataSet, G4PhotoNuclearCrossSection, G4NeutronInelasticCrossSection, G4TripathiCrossSection, and G4PiNuclearCrossSection.

Definition at line 135 of file G4VCrossSectionDataSet.cc.

References G4cout, G4endl, G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4Material::GetName(), G4ParticleDefinition::GetParticleName(), and python.hepunit::MeV.

Referenced by G4CrossSectionPairGG::BuildPhysicsTable(), ComputeCrossSection(), and G4CrossSectionPairGG::GetElementCrossSection().

{
  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;
}

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 G4DPMJET2_5CrossSection, G4BGGPionInelasticXS, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4LENDCrossSection, G4NeutronHPThermalScatteringData, G4GlauberGribovCrossSection, G4IonsShenCrossSection, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionPlusInelasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusInelasticXS, G4NeutronHPElasticData, G4ChipsAntiBaryonElasticXS, G4ChipsKaonMinusElasticXS, G4ChipsNeutronElasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4NeutronHPCaptureData, G4NeutronHPInelasticData, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsPionPlusElasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4NeutronHPFissionData, and G4ChipsAntiBaryonInelasticXS.

Definition at line 151 of file G4VCrossSectionDataSet.cc.

References G4cout, G4endl, G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4Element::GetName(), G4Material::GetName(), G4ParticleDefinition::GetParticleName(), and python.hepunit::MeV.

Referenced by G4QMDReaction::ApplyYourself(), and ComputeCrossSection().

{
  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;
}

G4double G4VCrossSectionDataSet::GetMaxKinEnergy ( ) const

inline

Definition at line 213 of file G4VCrossSectionDataSet.hh.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4NeutronHPFissionData::IsIsoApplicable(), G4NeutronHPInelasticData::IsIsoApplicable(), G4NeutronHPCaptureData::IsIsoApplicable(), G4NeutronHPElasticData::IsIsoApplicable(), and G4LENDCrossSection::IsIsoApplicable().

{
  return maxKinEnergy;
}

G4double G4VCrossSectionDataSet::GetMinKinEnergy ( ) const

inline

Definition at line 203 of file G4VCrossSectionDataSet.hh.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4NeutronHPFissionData::IsIsoApplicable(), G4NeutronHPInelasticData::IsIsoApplicable(), G4NeutronHPCaptureData::IsIsoApplicable(), G4NeutronHPElasticData::IsIsoApplicable(), and G4LENDCrossSection::IsIsoApplicable().

{
  return minKinEnergy;
}

const G4String & G4VCrossSectionDataSet::GetName ( void  ) const

inline

Definition at line 218 of file G4VCrossSectionDataSet.hh.

Referenced by G4CrossSectionPairGG::BuildPhysicsTable(), G4LENDCrossSection::create_used_target_map(), G4CrossSectionPairGG::DumpPhysicsTable(), G4LENDCrossSection::DumpPhysicsTable(), and G4CrossSectionDataSetRegistry::GetCrossSectionDataSet().

{
  return name;
}

G4int G4VCrossSectionDataSet::GetVerboseLevel ( ) const

inlinevirtual

Reimplemented in G4NeutronHPElasticData, G4NeutronHPCaptureData, G4NeutronHPInelasticData, and G4NeutronHPFissionData.

Definition at line 188 of file G4VCrossSectionDataSet.hh.

References verboseLevel.

{
  return verboseLevel;
}

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

virtual

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

Definition at line 66 of file G4VCrossSectionDataSet.cc.

Referenced by G4CrossSectionPairGG::BuildPhysicsTable(), ComputeCrossSection(), and G4CrossSectionPairGG::IsElementApplicable().

{
  return false;
}

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

virtual

Reimplemented in G4DPMJET2_5CrossSection, G4BGGPionInelasticXS, G4LENDCrossSection, G4BGGNucleonElasticXS, G4BGGNucleonInelasticXS, G4BGGPionElasticXS, G4NeutronHPBGGNucleonInelasticXS, G4NeutronInelasticXS, G4NeutronCaptureXS, G4NeutronHPThermalScatteringData, G4GlauberGribovCrossSection, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionPlusInelasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusInelasticXS, G4NeutronHPElasticData, G4ChipsAntiBaryonElasticXS, G4ChipsKaonMinusElasticXS, G4ChipsNeutronElasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4NeutronHPCaptureData, G4NeutronHPInelasticData, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsPionPlusElasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4NeutronHPFissionData, and G4ChipsAntiBaryonInelasticXS.

Definition at line 74 of file G4VCrossSectionDataSet.cc.

Referenced by ComputeCrossSection().

{
  return false;
}

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

virtual

Reimplemented in G4NeutronInelasticXS, and G4NeutronCaptureXS.

Definition at line 170 of file G4VCrossSectionDataSet.cc.

References G4UniformRand, G4Element::GetIsotopeVector(), G4Element::GetNumberOfIsotopes(), and G4Element::GetRelativeAbundanceVector().

{
  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;
}

void G4VCrossSectionDataSet::SetMaxKinEnergy ( G4double  value )

inline

Definition at line 208 of file G4VCrossSectionDataSet.hh.

Referenced by G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionElastic::G4CrossSectionElastic(), G4CrossSectionInelastic::G4CrossSectionInelastic(), G4LENDCrossSection::G4LENDCrossSection(), G4NeutronHPCaptureData::G4NeutronHPCaptureData(), G4NeutronHPElasticData::G4NeutronHPElasticData(), G4NeutronHPFissionData::G4NeutronHPFissionData(), G4NeutronHPInelasticData::G4NeutronHPInelasticData(), G4NeutronHPThermalScatteringData::G4NeutronHPThermalScatteringData(), and G4PiNuclearCrossSection::G4PiNuclearCrossSection().

{
  maxKinEnergy = value;
}

void G4VCrossSectionDataSet::SetMinKinEnergy ( G4double  value )

inline

Definition at line 198 of file G4VCrossSectionDataSet.hh.

Referenced by G4IonElasticPhysics::ConstructProcess(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionElastic::G4CrossSectionElastic(), G4CrossSectionInelastic::G4CrossSectionInelastic(), G4LENDCrossSection::G4LENDCrossSection(), G4NeutronHPCaptureData::G4NeutronHPCaptureData(), G4NeutronHPElasticData::G4NeutronHPElasticData(), G4NeutronHPFissionData::G4NeutronHPFissionData(), G4NeutronHPInelasticData::G4NeutronHPInelasticData(), G4NeutronHPThermalScatteringData::G4NeutronHPThermalScatteringData(), and G4PiNuclearCrossSection::G4PiNuclearCrossSection().

{
  minKinEnergy = value;
}

void G4VCrossSectionDataSet::SetName ( const G4String &  nam )

inlineprotected

Definition at line 223 of file G4VCrossSectionDataSet.hh.

{
  name = nam;
}

void G4VCrossSectionDataSet::SetVerboseLevel ( G4int  value )

inlinevirtual

Reimplemented in G4NeutronHPCaptureData, G4NeutronHPInelasticData, G4NeutronHPElasticData, and G4NeutronHPFissionData.

Definition at line 193 of file G4VCrossSectionDataSet.hh.

References verboseLevel.

{
  verboseLevel = value;
}

Field Documentation

G4int G4VCrossSectionDataSet::verboseLevel

protected

Definition at line 166 of file G4VCrossSectionDataSet.hh.

Referenced by G4CrossSectionPairGG::BuildPhysicsTable(), G4NeutronElasticXS::BuildPhysicsTable(), G4NeutronCaptureXS::BuildPhysicsTable(), G4NeutronInelasticXS::BuildPhysicsTable(), G4BGGNucleonInelasticXS::BuildPhysicsTable(), G4BGGPionElasticXS::BuildPhysicsTable(), G4BGGNucleonElasticXS::BuildPhysicsTable(), G4BGGPionInelasticXS::BuildPhysicsTable(), G4LENDCrossSection::create_used_target_map(), G4BGGNucleonElasticXS::G4BGGNucleonElasticXS(), G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionPairGG::G4CrossSectionPairGG(), G4NeutronCaptureXS::G4NeutronCaptureXS(), G4NeutronElasticXS::G4NeutronElasticXS(), G4NeutronInelasticXS::G4NeutronInelasticXS(), G4NeutronElasticXS::GetElementCrossSection(), G4CrossSectionPairGG::GetElementCrossSection(), G4NeutronCaptureXS::GetElementCrossSection(), G4NeutronInelasticXS::GetElementCrossSection(), G4BGGNucleonInelasticXS::GetElementCrossSection(), G4BGGPionElasticXS::GetElementCrossSection(), G4BGGNucleonElasticXS::GetElementCrossSection(), G4BGGPionInelasticXS::GetElementCrossSection(), G4BGGNucleonInelasticXS::GetIsoCrossSection(), G4BGGPionElasticXS::GetIsoCrossSection(), G4BGGNucleonElasticXS::GetIsoCrossSection(), G4BGGPionInelasticXS::GetIsoCrossSection(), GetVerboseLevel(), and SetVerboseLevel().

---

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

• G4VCrossSectionDataSet.hh
• G4VCrossSectionDataSet.cc