#include <G4ParticleHPInelasticData.hh>

Inheritance diagram for G4ParticleHPInelasticData:

Public Member Functions
void	BuildPhysicsTable (const G4ParticleDefinition &)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

virtual void	CrossSectionDescription (std::ostream &) const

void	DumpPhysicsTable (const G4ParticleDefinition &)

bool	ForAllAtomsAndEnergies () const

	G4ParticleHPInelasticData (G4ParticleDefinition *projectile=G4Neutron::Neutron())

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , G4double aT)

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat=nullptr)

G4double	GetIsoCrossSection (const G4DynamicParticle , G4int, G4int, const G4Isotope , const G4Element , const G4Material )

G4double	GetMaxKinEnergy () const

G4double	GetMinKinEnergy () const

const G4String &	GetName () const

G4ParticleDefinition *	GetProjectile ()

G4int	GetVerboseLevel () const

virtual G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material mat=nullptr)

G4bool	IsIsoApplicable (const G4DynamicParticle , G4int, G4int, const G4Element , const G4Material *)

virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)

void	SetForAllAtomsAndEnergies (G4bool val)

void	SetMaxKinEnergy (G4double value)

void	SetMinKinEnergy (G4double value)

void	SetName (const G4String &nam)

void	SetVerboseLevel (G4int)

	~G4ParticleHPInelasticData ()

Protected Attributes
G4String	name

G4int	verboseLevel

Private Attributes
const G4Element *	element_cache

G4bool	instanceOfWorker

G4bool	isForAllAtomsAndEnergies

G4double	ke_cache

const G4Material *	material_cache

G4double	maxKinEnergy

G4double	minKinEnergy

G4CrossSectionDataSetRegistry *	registry

G4PhysicsTable *	theCrossSections

G4ParticleHPData *	theHPData

G4ParticleDefinition *	theProjectile

G4double	xs_cache

Detailed Description

Definition at line 53 of file G4ParticleHPInelasticData.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPInelasticData()

G4ParticleHPInelasticData::G4ParticleHPInelasticData ( G4ParticleDefinition * projectile = G4Neutron::Neutron() )

Definition at line 47 of file G4ParticleHPInelasticData.cc.

  : G4VCrossSectionDataSet("")
{
  const char* dataDirVariable;
  G4String particleName;
  if( projectile == G4Neutron::Neutron() ) {
    dataDirVariable = "G4NEUTRONHPDATA";
  }else if( projectile == G4Proton::Proton() ) {
    dataDirVariable = "G4PROTONHPDATA";
    particleName = "Proton";
  }else if( projectile == G4Deuteron::Deuteron() ) {
    dataDirVariable = "G4DEUTERONHPDATA";
    particleName = "Deuteron";
  }else if( projectile == G4Triton::Triton() ) {
    dataDirVariable = "G4TRITONHPDATA";
    particleName = "Triton";
  }else if( projectile == G4He3::He3() ) {
    dataDirVariable = "G4HE3HPDATA";
    particleName = "He3";
  }else if( projectile == G4Alpha::Alpha() ) {
    dataDirVariable = "G4ALPHAHPDATA";
    particleName = "Alpha";
  } else {
    G4String message("G4ParticleHPInelasticData may only be called for neutron, proton, deuteron, triton, He3 or alpha, while it is called for " + projectile->GetParticleName());
    throw G4HadronicException(__FILE__, __LINE__,message.c_str());
  }
  //  G4cout << this << " G4ParticleHPInelasticData::G4ParticleHPInelasticData " << projectile->GetParticleName() << " DATADIR " << dataDirVariable << G4endl;//GDEB
  G4String dataName = projectile->GetParticleName()+"HPInelasticXS";
  dataName.at(0) = toupper(dataName.at(0)) ;
  SetName( dataName );
 
  if ( !std::getenv(dataDirVariable) && !std::getenv( "G4PARTICLEHPDATA" ) ){
    G4String message("Please setenv G4PARTICLEHPDATA (recommended) or, at least setenv " +
                     G4String(dataDirVariable) + " to point to the " + projectile->GetParticleName() + " cross-section files.");
    throw G4HadronicException(__FILE__, __LINE__,message.c_str());
  }
 
  G4String dirName;
  if ( std::getenv(dataDirVariable) ) {
    dirName = std::getenv(dataDirVariable);
  } else {
    G4String baseName = std::getenv( "G4PARTICLEHPDATA" );
    dirName = baseName + "/" + particleName;
  }
  #ifdef G4VERBOSE
  if ( G4HadronicParameters::Instance()->GetVerboseLevel() > 0 ) {
    G4cout << "@@@ G4ParticleHPInelasticData instantiated for particle " << projectile->GetParticleName() << " data directory variable is " << dataDirVariable << " pointing to " << dirName << G4endl;
  }
  #endif
 
  SetMinKinEnergy( 0*CLHEP::MeV );                                   
  SetMaxKinEnergy( 20*CLHEP::MeV );                                   
 
  theCrossSections = 0;
  theProjectile=projectile;
 
  theHPData = NULL;
  instanceOfWorker = false;
  if ( G4Threading::IsMasterThread() ) {
    theHPData = new G4ParticleHPData( theProjectile ); 
  } else {
    instanceOfWorker = true;
  }
  element_cache = NULL;
  material_cache = NULL;
  ke_cache = 0.0; 
  xs_cache = 0.0; 
}

References G4Alpha::Alpha(), G4Deuteron::Deuteron(), element_cache, G4cout, G4endl, G4ParticleDefinition::GetParticleName(), GetVerboseLevel(), G4He3::He3(), G4HadronicParameters::Instance(), instanceOfWorker, G4Threading::IsMasterThread(), ke_cache, material_cache, CLHEP::MeV, G4Neutron::Neutron(), G4Proton::Proton(), G4VCrossSectionDataSet::SetMaxKinEnergy(), G4VCrossSectionDataSet::SetMinKinEnergy(), G4VCrossSectionDataSet::SetName(), theCrossSections, theHPData, theProjectile, G4Triton::Triton(), and xs_cache.

◆ ~G4ParticleHPInelasticData()

G4ParticleHPInelasticData::~G4ParticleHPInelasticData ( )

Definition at line 116 of file G4ParticleHPInelasticData.cc.

{
   if ( theCrossSections != NULL && instanceOfWorker != true ) {
     theCrossSections->clearAndDestroy();
     delete theCrossSections;
     theCrossSections = NULL;
   }
   if ( theHPData != NULL && instanceOfWorker != true ) {
     delete theHPData;
     theHPData = NULL;
   }
}

References G4PhysicsTable::clearAndDestroy(), instanceOfWorker, theCrossSections, and theHPData.

Member Function Documentation

◆ BuildPhysicsTable()

void G4ParticleHPInelasticData::BuildPhysicsTable ( const G4ParticleDefinition & projectile )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 169 of file G4ParticleHPInelasticData.cc.

{
  //  if(&projectile!=G4Neutron::Neutron()) 
  //     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
 
   if ( G4Threading::IsWorkerThread() ) {
      theCrossSections = G4ParticleHPManager::GetInstance()->GetInelasticCrossSections( &projectile );
      return;
   } else {
      if ( theHPData == NULL ) theHPData = G4ParticleHPData::Instance( const_cast<G4ParticleDefinition*> ( &projectile ) ); 
   }
 
  size_t numberOfElements = G4Element::GetNumberOfElements();
//  theCrossSections = new G4PhysicsTable( numberOfElements );
// TKDB
  //if ( theCrossSections == 0 )
  //{ theCrossSections = new G4PhysicsTable( numberOfElements ); }
  if ( theCrossSections == NULL ) 
    theCrossSections = new G4PhysicsTable( numberOfElements );
  else
    theCrossSections->clearAndDestroy();
 
  // make a PhysicsVector for each element
 
  //G4ParticleHPData* hpData = new G4ParticleHPData(projectile); //NEW
  static G4ThreadLocal G4ElementTable *theElementTable  = 0 ;
  if (!theElementTable) theElementTable= G4Element::GetElementTable();
  for( size_t i=0; i<numberOfElements; ++i )
  {
    //NEW    G4PhysicsVector* physVec = G4ParticleHPData::
    //NEW      Instance(projectile, dataDirVariable)->MakePhysicsVector((*theElementTable)[i], this);
    //G4PhysicsVector* physVec = hpData->MakePhysicsVector((*theElementTable)[i], this);
    G4PhysicsVector* physVec = theHPData->MakePhysicsVector((*theElementTable)[i], this);
    theCrossSections->push_back(physVec);
  }
 
   G4ParticleHPManager::GetInstance()->RegisterInelasticCrossSections( &projectile , theCrossSections );
}

References G4PhysicsTable::clearAndDestroy(), G4ThreadLocal, G4Element::GetElementTable(), G4ParticleHPManager::GetInelasticCrossSections(), G4ParticleHPManager::GetInstance(), G4Element::GetNumberOfElements(), G4ParticleHPData::Instance(), G4Threading::IsWorkerThread(), G4ParticleHPData::MakePhysicsVector(), G4PhysicsTable::push_back(), G4ParticleHPManager::RegisterInelasticCrossSections(), theCrossSections, and theHPData.

◆ ComputeCrossSection()

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

inherited

Definition at line 81 of file G4VCrossSectionDataSet.cc.

{
  G4int Z = elm->GetZasInt();
 
  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 
  size_t nIso = elm->GetNumberOfIsotopes();    
  G4double fact = 0.0;
  G4double xsec = 0.0;
 
  // user-defined isotope abundances        
  const G4IsotopeVector* isoVector = elm->GetIsotopeVector();
  const G4double* abundVector = elm->GetRelativeAbundanceVector();
 
  for (size_t j=0; j<nIso; ++j) {
    const G4Isotope* 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);
    }
  }
  return (fact > 0.0) ? xsec/fact : 0.0; 
}

References A, G4VCrossSectionDataSet::GetElementCrossSection(), G4VCrossSectionDataSet::GetIsoCrossSection(), G4Element::GetIsotopeVector(), G4Isotope::GetN(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZasInt(), G4VCrossSectionDataSet::IsElementApplicable(), G4VCrossSectionDataSet::IsIsoApplicable(), and Z.

Referenced by G4VCrossSectionDataSet::GetCrossSection().

◆ CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 370 of file G4ParticleHPInelasticData.cc.

{
   outFile << "Extension of High Precision cross section for inelastic reaction of proton, deuteron, triton, He3 and alpha below 20MeV\n";
}

◆ DumpPhysicsTable()

void G4ParticleHPInelasticData::DumpPhysicsTable ( const G4ParticleDefinition & projectile )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 208 of file G4ParticleHPInelasticData.cc.

{
  if(&projectile!=theProjectile) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use ParticleHP data for a wrong projectile!!!");  
 
  #ifdef G4VERBOSE
  if ( G4HadronicParameters::Instance()->GetVerboseLevel() == 0 ) return;
//
// Dump element based cross section
// range 10e-5 eV to 20 MeV
// 10 point per decade
// in barn
//
 
   G4cout << G4endl;
   G4cout << G4endl;
   G4cout << "Inelastic Cross Section of Neutron HP"<< G4endl;
   G4cout << "(Pointwise cross-section at 0 Kelvin.)" << G4endl;
   G4cout << G4endl;
   G4cout << "Name of Element" << G4endl;
   G4cout << "Energy[eV]  XS[barn]" << G4endl;
   G4cout << G4endl;
 
   size_t numberOfElements = G4Element::GetNumberOfElements();
   static G4ThreadLocal G4ElementTable *theElementTable  = 0 ;
   if (!theElementTable) theElementTable= G4Element::GetElementTable();
 
   for ( size_t i = 0 ; i < numberOfElements ; ++i )
   {
 
      G4cout << (*theElementTable)[i]->GetName() << G4endl;
 
      G4int ie = 0;
 
      for ( ie = 0 ; ie < 130 ; ie++ )
      {
    G4double eKinetic = 1.0e-5 * G4Pow::GetInstance()->powA ( 10.0 , ie/10.0 ) *CLHEP::eV;
         G4bool outOfRange = false;
 
         if ( eKinetic < 20*CLHEP::MeV )
         {
       G4cout << eKinetic/CLHEP::eV << " " << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/CLHEP::barn << G4endl;
         }
 
      }
 
      G4cout << G4endl;
   }
 
  //G4cout << "G4ParticleHPInelasticData::DumpPhysicsTable still to be implemented"<<G4endl;
  #endif
}

References CLHEP::barn, CLHEP::eV, G4cout, G4endl, G4ThreadLocal, G4Element::GetElementTable(), G4Pow::GetInstance(), G4Element::GetNumberOfElements(), GetVerboseLevel(), G4HadronicParameters::Instance(), CLHEP::MeV, G4Pow::powA(), and theProjectile.

◆ ForAllAtomsAndEnergies()

bool G4VCrossSectionDataSet::ForAllAtomsAndEnergies ( ) const

inlineinherited

Definition at line 230 of file G4VCrossSectionDataSet.hh.

{
  return isForAllAtomsAndEnergies;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4CrossSectionDataStore::AddDataSet().

◆ GetCrossSection() [1/2]

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

inlineinherited

Definition at line 187 of file G4VCrossSectionDataSet.hh.

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

References G4VCrossSectionDataSet::ComputeCrossSection().

◆ GetCrossSection() [2/2]

G4double G4ParticleHPInelasticData::GetCrossSection	(	const G4DynamicParticle *	projectile,
		const G4Element *	anE,
		G4double	aT
	)

Definition at line 263 of file G4ParticleHPInelasticData.cc.

{
  G4double result = 0;
  G4bool outOfRange;
  G4int index = anE->GetIndex();
 
  // prepare neutron
  G4double eKinetic = projectile->GetKineticEnergy();
 
   if ( G4ParticleHPManager::GetInstance()->GetNeglectDoppler() )
  {
     //NEGLECT_DOPPLER
     G4double factor = 1.0;
     if ( eKinetic < aT * CLHEP::k_Boltzmann ) 
     {
        // below 0.1 eV neutrons 
        // Have to do some, but now just igonre.   
        // Will take care after performance check.  
        // factor = factor * targetV;
     }
     return ( (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) )* factor;
  }   
 
  G4ReactionProduct theNeutron( projectile->GetDefinition() );
  theNeutron.SetMomentum( projectile->GetMomentum() );
  theNeutron.SetKineticEnergy( eKinetic );
 
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eps = 0.0001;
  G4double theA = anE->GetN();
  G4double theZ = anE->GetZ();
  G4double eleMass; 
  eleMass = G4NucleiProperties::GetNuclearMass(static_cast<G4int>(theA+eps), static_cast<G4int>(theZ+eps) );
  
  G4ReactionProduct boosted;
  G4double aXsection;
  
  // MC integration loop
  G4int counter = 0;
  G4int failCount = 0;
  G4double buffer = 0; G4int size = G4int(std::max(10., aT/60*CLHEP::kelvin));
  G4ThreeVector neutronVelocity = 1./theProjectile->GetPDGMass()*theNeutron.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
 
  //  G4cout << " G4ParticleHPInelasticData 2 " << size << G4endl;//GDEB
#ifndef G4PHP_DOPPLER_LOOP_ONCE
  while(counter == 0 || std::abs(buffer-result/std::max(1,counter)) > 0.01*buffer) // Loop checking, 11.05.2015, T. Koi
  {
    if(counter) buffer = result/counter;
    while (counter<size) // Loop checking, 11.05.2015, T. Koi
    {
      counter ++;
#endif
      //G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass/theProjectile->GetPDGMass(), aT );
      //G4Nucleus::GetThermalNucleus requests normalization of mass in neutron mass
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass/G4Neutron::Neutron()->GetPDGMass(), aT );
      boosted.Lorentz(theNeutron, aThermalNuc);
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
      //       G4cout << " G4ParticleHPInelasticData aXsection " << aXsection << " index " << index << " theEkin " << theEkin << " outOfRange " << outOfRange <<G4endl;//GDEB
     if(aXsection <0) 
      {
        if(failCount<1000)
    {
      failCount++;
#ifndef G4PHP_DOPPLER_LOOP_ONCE
      counter--;
      continue;
#endif
    }
    else
    {
      aXsection = 0;
    }
      }
      // velocity correction.
      G4ThreeVector targetVelocity = 1./aThermalNuc.GetMass()*aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity-neutronVelocity).mag()/neutronVMag;
      result += aXsection;
    }
#ifndef G4PHP_DOPPLER_LOOP_ONCE
    size += size;
  }
  result /= counter;
#endif
 
/*
  // Checking impact of  G4PHP_NEGLECT_DOPPLER
  G4cout << " result " << result << " " 
         << (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) << " " 
         << (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) /result << G4endl;
*/
//   G4cout << this << " G4ParticleHPInelasticData result " << result << G4endl; //GDEB
 
  return result;
}

References buffer, eps, G4DynamicParticle::GetDefinition(), G4Element::GetIndex(), G4ParticleHPManager::GetInstance(), G4DynamicParticle::GetKineticEnergy(), G4ReactionProduct::GetKineticEnergy(), G4ReactionProduct::GetMass(), G4DynamicParticle::GetMomentum(), G4ReactionProduct::GetMomentum(), G4Element::GetN(), G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGMass(), G4Nucleus::GetThermalNucleus(), G4Element::GetZ(), CLHEP::k_Boltzmann, CLHEP::kelvin, G4ReactionProduct::Lorentz(), CLHEP::Hep3Vector::mag(), G4INCL::Math::max(), G4Neutron::Neutron(), G4ReactionProduct::SetKineticEnergy(), G4ReactionProduct::SetMomentum(), theCrossSections, and theProjectile.

Referenced by GetIsoCrossSection().

◆ GetElementCrossSection()

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

virtualinherited

Reimplemented in G4EMDissociationCrossSection, G4IonsShenCrossSection, G4NeutrinoElectronCcXsc, G4NeutrinoElectronNcXsc, G4NeutrinoElectronTotXsc, G4NeutronElectronElXsc, G4PhotoNuclearCrossSection, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticXS, G4ElectroNuclearCrossSection, G4BGGNucleonElasticXS, G4BGGPionElasticXS, G4BGGPionInelasticXS, G4BGGNucleonInelasticXS, G4CrossSectionElastic, G4CrossSectionInelastic, G4GammaNuclearXS, G4ParticleInelasticXS, G4ZeroXS, G4NucleonNuclearCrossSection, G4MuNeutrinoNucleusTotXsc, and G4KokoulinMuonNuclearXS.

Definition at line 114 of file G4VCrossSectionDataSet.cc.

{
  G4ExceptionDescription ed;
  ed << "GetElementCrossSection is not implemented in <" << name << ">\n"
     << "Particle: " << dynPart->GetDefinition()->GetParticleName()
     << "  Ekin(MeV)= "  << dynPart->GetKineticEnergy()/MeV;
  if(mat) { ed << "  material: " << mat->GetName(); }
  ed << " target Z= " << Z << G4endl;
  G4Exception("G4VCrossSectionDataSet::GetElementCrossSection", "had001", 
              FatalException, ed);
  return 0.0;
}

References FatalException, G4endl, G4Exception(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4Material::GetName(), G4ParticleDefinition::GetParticleName(), MeV, G4VCrossSectionDataSet::name, and Z.

Referenced by G4QMDReaction::ApplyYourself(), G4VCrossSectionDataSet::ComputeCrossSection(), G4GammaNuclearXS::GetElementCrossSection(), G4GammaNuclearXS::GetIsoCrossSection(), and G4GammaNuclearXS::Initialise().

◆ GetIsoCrossSection()

G4double G4ParticleHPInelasticData::GetIsoCrossSection	(	const G4DynamicParticle *	dp,
		G4int	,
		G4int	,
		const G4Isotope *	,
		const G4Element *	element,
		const G4Material *	material
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 142 of file G4ParticleHPInelasticData.cc.

{
   if ( dp->GetKineticEnergy() == ke_cache && element == element_cache &&  material == material_cache ) return xs_cache;
 
   ke_cache = dp->GetKineticEnergy();
   element_cache = element;
   material_cache = material;
   G4double xs = GetCrossSection( dp , element , material->GetTemperature() );
   xs_cache = xs;
   return xs;
}

References element_cache, GetCrossSection(), G4DynamicParticle::GetKineticEnergy(), ke_cache, eplot::material, material_cache, and xs_cache.

◆ GetMaxKinEnergy()

G4double G4VCrossSectionDataSet::GetMaxKinEnergy ( ) const

inlineinherited

Definition at line 220 of file G4VCrossSectionDataSet.hh.

{
  return maxKinEnergy;
}

References G4VCrossSectionDataSet::maxKinEnergy.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4LENDCrossSection::IsIsoApplicable(), G4ParticleHPCaptureData::IsIsoApplicable(), G4ParticleHPElasticData::IsIsoApplicable(), G4ParticleHPFissionData::IsIsoApplicable(), and IsIsoApplicable().

◆ GetMinKinEnergy()

G4double G4VCrossSectionDataSet::GetMinKinEnergy ( ) const

inlineinherited

Definition at line 210 of file G4VCrossSectionDataSet.hh.

{
  return minKinEnergy;
}

References G4VCrossSectionDataSet::minKinEnergy.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4LENDCrossSection::IsIsoApplicable(), G4ParticleHPCaptureData::IsIsoApplicable(), G4ParticleHPElasticData::IsIsoApplicable(), G4ParticleHPFissionData::IsIsoApplicable(), and IsIsoApplicable().

◆ GetName()

const G4String & G4VCrossSectionDataSet::GetName ( ) const

inlineinherited

Definition at line 225 of file G4VCrossSectionDataSet.hh.

{
  return name;
}

References G4VCrossSectionDataSet::name.

Referenced by G4LENDCrossSection::DumpLENDTargetInfo(), G4LENDCrossSection::DumpPhysicsTable(), G4LENDCrossSection::GetIsoCrossSection(), and G4CrossSectionDataStore::PrintCrossSectionHtml().

◆ GetProjectile()

G4ParticleDefinition * G4ParticleHPInelasticData::GetProjectile ( )

inline

Definition at line 85 of file G4ParticleHPInelasticData.hh.

85{return theProjectile;}

References theProjectile.

◆ GetVerboseLevel()

G4int G4ParticleHPInelasticData::GetVerboseLevel ( ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 362 of file G4ParticleHPInelasticData.cc.

{
   return G4ParticleHPManager::GetInstance()->GetVerboseLevel();
}

References G4ParticleHPManager::GetInstance(), and G4ParticleHPManager::GetVerboseLevel().

Referenced by DumpPhysicsTable(), and G4ParticleHPInelasticData().

◆ IsElementApplicable()

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

virtualinherited

◆ IsIsoApplicable()

G4bool G4ParticleHPInelasticData::IsIsoApplicable	(	const G4DynamicParticle *	dp,
		G4int	,
		G4int	,
		const G4Element *	,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 129 of file G4ParticleHPInelasticData.cc.

{
   G4double eKin = dp->GetKineticEnergy();
   if ( eKin > GetMaxKinEnergy() 
     || eKin < GetMinKinEnergy() 
     || dp->GetDefinition() != theProjectile ) return false;                                   
 
   return true;
}

References G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4VCrossSectionDataSet::GetMaxKinEnergy(), G4VCrossSectionDataSet::GetMinKinEnergy(), and theProjectile.

◆ SelectIsotope()

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

virtualinherited

Reimplemented in G4GammaNuclearXS, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticXS, and G4ParticleInelasticXS.

Definition at line 149 of file G4VCrossSectionDataSet.cc.

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

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

◆ SetForAllAtomsAndEnergies()

void G4VCrossSectionDataSet::SetForAllAtomsAndEnergies ( G4bool val )

inlineinherited

Definition at line 235 of file G4VCrossSectionDataSet.hh.

{
  isForAllAtomsAndEnergies = val;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4BGGNucleonElasticXS::G4BGGNucleonElasticXS(), G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4ElectroNuclearCrossSection::G4ElectroNuclearCrossSection(), G4GammaNuclearXS::G4GammaNuclearXS(), G4NeutronElasticXS::G4NeutronElasticXS(), G4NeutronInelasticXS::G4NeutronInelasticXS(), G4ParticleInelasticXS::G4ParticleInelasticXS(), and G4PhotoNuclearCrossSection::G4PhotoNuclearCrossSection().

◆ SetMaxKinEnergy()

void G4VCrossSectionDataSet::SetMaxKinEnergy ( G4double value )

inlineinherited

Definition at line 215 of file G4VCrossSectionDataSet.hh.

{
  maxKinEnergy = value;
}

References G4VCrossSectionDataSet::maxKinEnergy.

Referenced by G4AlphaPHPBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4He3PHPBuilder::Build(), G4ProtonPHPBuilder::Build(), G4TritonPHPBuilder::Build(), G4CrossSectionElastic::BuildPhysicsTable(), G4CrossSectionInelastic::BuildPhysicsTable(), G4ParticleInelasticXS::BuildPhysicsTable(), G4IonPhysicsPHP::ConstructProcess(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionElastic::G4CrossSectionElastic(), G4CrossSectionInelastic::G4CrossSectionInelastic(), G4LENDCrossSection::G4LENDCrossSection(), G4ParticleHPCaptureData::G4ParticleHPCaptureData(), G4ParticleHPElasticData::G4ParticleHPElasticData(), G4ParticleHPFissionData::G4ParticleHPFissionData(), G4ParticleHPInelasticData(), and G4ParticleHPThermalScatteringData::G4ParticleHPThermalScatteringData().

◆ SetMinKinEnergy()

void G4VCrossSectionDataSet::SetMinKinEnergy ( G4double value )

inlineinherited

Definition at line 205 of file G4VCrossSectionDataSet.hh.

{
  minKinEnergy = value;
}

References G4VCrossSectionDataSet::minKinEnergy.

Referenced by G4AlphaPHPBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4He3PHPBuilder::Build(), G4ProtonPHPBuilder::Build(), G4TritonPHPBuilder::Build(), G4IonElasticPhysics::ConstructProcess(), G4IonPhysicsPHP::ConstructProcess(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionElastic::G4CrossSectionElastic(), G4CrossSectionInelastic::G4CrossSectionInelastic(), G4LENDCrossSection::G4LENDCrossSection(), G4ParticleHPCaptureData::G4ParticleHPCaptureData(), G4ParticleHPElasticData::G4ParticleHPElasticData(), G4ParticleHPFissionData::G4ParticleHPFissionData(), G4ParticleHPInelasticData(), and G4ParticleHPThermalScatteringData::G4ParticleHPThermalScatteringData().

◆ SetName()

void G4VCrossSectionDataSet::SetName ( const G4String & nam )

inlineinherited

Definition at line 240 of file G4VCrossSectionDataSet.hh.

{
  name = nam;
}

References G4VCrossSectionDataSet::name.

Referenced by G4ParticleHPInelasticData().

◆ SetVerboseLevel()

void G4ParticleHPInelasticData::SetVerboseLevel ( G4int newValue )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 366 of file G4ParticleHPInelasticData.cc.

{
   G4ParticleHPManager::GetInstance()->SetVerboseLevel(newValue);
}

References G4ParticleHPManager::GetInstance(), and G4ParticleHPManager::SetVerboseLevel().

Field Documentation

◆ element_cache

const G4Element* G4ParticleHPInelasticData::element_cache

private

Definition at line 104 of file G4ParticleHPInelasticData.hh.

Referenced by G4ParticleHPInelasticData(), and GetIsoCrossSection().

◆ instanceOfWorker

G4bool G4ParticleHPInelasticData::instanceOfWorker

private

Definition at line 100 of file G4ParticleHPInelasticData.hh.

Referenced by G4ParticleHPInelasticData(), and ~G4ParticleHPInelasticData().

◆ isForAllAtomsAndEnergies

G4bool G4VCrossSectionDataSet::isForAllAtomsAndEnergies

privateinherited

Definition at line 182 of file G4VCrossSectionDataSet.hh.

Referenced by G4VCrossSectionDataSet::ForAllAtomsAndEnergies(), and G4VCrossSectionDataSet::SetForAllAtomsAndEnergies().

◆ ke_cache

G4double G4ParticleHPInelasticData::ke_cache

private

Definition at line 102 of file G4ParticleHPInelasticData.hh.

Referenced by G4ParticleHPInelasticData(), and GetIsoCrossSection().

◆ material_cache

const G4Material* G4ParticleHPInelasticData::material_cache

private

Definition at line 105 of file G4ParticleHPInelasticData.hh.

Referenced by G4ParticleHPInelasticData(), and GetIsoCrossSection().

◆ maxKinEnergy

G4double G4VCrossSectionDataSet::maxKinEnergy

privateinherited

Definition at line 180 of file G4VCrossSectionDataSet.hh.

Referenced by G4VCrossSectionDataSet::GetMaxKinEnergy(), and G4VCrossSectionDataSet::SetMaxKinEnergy().

◆ minKinEnergy

G4double G4VCrossSectionDataSet::minKinEnergy

privateinherited

Definition at line 179 of file G4VCrossSectionDataSet.hh.

Referenced by G4VCrossSectionDataSet::GetMinKinEnergy(), and G4VCrossSectionDataSet::SetMinKinEnergy().

◆ name

G4String G4VCrossSectionDataSet::name

protectedinherited

Definition at line 170 of file G4VCrossSectionDataSet.hh.

Referenced by source.g4viscp.G4Scene::create_scene(), G4VCrossSectionDataSet::GetElementCrossSection(), G4VCrossSectionDataSet::GetIsoCrossSection(), G4VCrossSectionDataSet::GetName(), mcscore.MCParticle::printout(), G4VCrossSectionDataSet::SetName(), and source.g4viscp.G4Scene::update_scene().

◆ registry

G4CrossSectionDataSetRegistry* G4VCrossSectionDataSet::registry

privateinherited

Definition at line 177 of file G4VCrossSectionDataSet.hh.

Referenced by G4VCrossSectionDataSet::G4VCrossSectionDataSet(), and G4VCrossSectionDataSet::~G4VCrossSectionDataSet().

◆ theCrossSections

G4PhysicsTable* G4ParticleHPInelasticData::theCrossSections

private

Definition at line 94 of file G4ParticleHPInelasticData.hh.

Referenced by BuildPhysicsTable(), G4ParticleHPInelasticData(), GetCrossSection(), and ~G4ParticleHPInelasticData().

◆ theHPData

G4ParticleHPData* G4ParticleHPInelasticData::theHPData

private

Definition at line 98 of file G4ParticleHPInelasticData.hh.

Referenced by BuildPhysicsTable(), G4ParticleHPInelasticData(), and ~G4ParticleHPInelasticData().

◆ theProjectile

G4ParticleDefinition* G4ParticleHPInelasticData::theProjectile

private

Definition at line 96 of file G4ParticleHPInelasticData.hh.

Referenced by DumpPhysicsTable(), G4ParticleHPInelasticData(), GetCrossSection(), GetProjectile(), and IsIsoApplicable().

◆ verboseLevel

G4int G4VCrossSectionDataSet::verboseLevel

protectedinherited

Definition at line 168 of file G4VCrossSectionDataSet.hh.

Referenced by G4BGGNucleonElasticXS::BuildPhysicsTable(), G4BGGPionElasticXS::BuildPhysicsTable(), G4BGGPionInelasticXS::BuildPhysicsTable(), G4GammaNuclearXS::BuildPhysicsTable(), G4NeutronCaptureXS::BuildPhysicsTable(), G4NeutronElasticXS::BuildPhysicsTable(), G4NeutronInelasticXS::BuildPhysicsTable(), G4ParticleInelasticXS::BuildPhysicsTable(), G4BGGNucleonInelasticXS::BuildPhysicsTable(), G4LENDCrossSection::create_used_target_map(), G4BGGNucleonElasticXS::G4BGGNucleonElasticXS(), G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4GammaNuclearXS::G4GammaNuclearXS(), G4NeutronCaptureXS::G4NeutronCaptureXS(), G4NeutronElasticXS::G4NeutronElasticXS(), G4NeutronInelasticXS::G4NeutronInelasticXS(), G4ParticleInelasticXS::G4ParticleInelasticXS(), G4NeutronCaptureXS::GetElementCrossSection(), G4NeutronElasticXS::GetElementCrossSection(), G4NeutronInelasticXS::GetElementCrossSection(), G4BGGNucleonElasticXS::GetElementCrossSection(), G4BGGPionElasticXS::GetElementCrossSection(), G4BGGPionInelasticXS::GetElementCrossSection(), G4BGGNucleonInelasticXS::GetElementCrossSection(), G4GammaNuclearXS::GetElementCrossSection(), G4ParticleInelasticXS::GetElementCrossSection(), G4BGGNucleonElasticXS::GetIsoCrossSection(), G4BGGPionElasticXS::GetIsoCrossSection(), G4BGGPionInelasticXS::GetIsoCrossSection(), G4GammaNuclearXS::GetIsoCrossSection(), G4BGGNucleonInelasticXS::GetIsoCrossSection(), G4VCrossSectionDataSet::GetVerboseLevel(), G4NeutronElasticXS::Initialise(), G4ParticleInelasticXS::IsoCrossSection(), G4NeutronCaptureXS::IsoCrossSection(), G4NeutronInelasticXS::IsoCrossSection(), G4GammaNuclearXS::RetrieveVector(), G4NeutronCaptureXS::RetrieveVector(), G4NeutronInelasticXS::RetrieveVector(), G4ParticleInelasticXS::RetrieveVector(), and G4VCrossSectionDataSet::SetVerboseLevel().

◆ xs_cache

G4double G4ParticleHPInelasticData::xs_cache

private

Definition at line 103 of file G4ParticleHPInelasticData.hh.

Referenced by G4ParticleHPInelasticData(), and GetIsoCrossSection().

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

source/processes/hadronic/models/particle_hp/include/G4ParticleHPInelasticData.hh
source/processes/hadronic/models/particle_hp/src/G4ParticleHPInelasticData.cc

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ParticleHPInelasticData()

◆ ~G4ParticleHPInelasticData()

Member Function Documentation

◆ BuildPhysicsTable()

◆ ComputeCrossSection()

◆ CrossSectionDescription()

◆ DumpPhysicsTable()

◆ ForAllAtomsAndEnergies()

◆ GetCrossSection() [1/2]

◆ GetCrossSection() [2/2]

◆ GetElementCrossSection()

◆ GetIsoCrossSection()

◆ GetMaxKinEnergy()

◆ GetMinKinEnergy()

◆ GetName()

◆ GetProjectile()

◆ GetVerboseLevel()

◆ IsElementApplicable()

◆ IsIsoApplicable()

◆ SelectIsotope()

◆ SetForAllAtomsAndEnergies()

◆ SetMaxKinEnergy()

◆ SetMinKinEnergy()

◆ SetName()

◆ SetVerboseLevel()

Field Documentation

◆ element_cache

◆ instanceOfWorker

◆ isForAllAtomsAndEnergies

◆ ke_cache

◆ material_cache

◆ maxKinEnergy

◆ minKinEnergy

◆ name

◆ registry

◆ theCrossSections

◆ theHPData

◆ theProjectile

◆ verboseLevel

◆ xs_cache