#include <G4ParticleHPThermalScatteringData.hh>

Inheritance diagram for G4ParticleHPThermalScatteringData:

Public Member Functions
void	AddUserThermalScatteringFile (G4String, G4String)

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

	G4ParticleHPThermalScatteringData ()

G4double	GetCoherentCrossSection (const G4DynamicParticle , const G4Element , const G4Material *)

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *)

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

G4double	GetIncoherentCrossSection (const G4DynamicParticle , const G4Element , const G4Material *)

G4double	GetInelasticCrossSection (const G4DynamicParticle , const G4Element , const G4Material *)

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

G4double	GetMaxKinEnergy () const

G4double	GetMinKinEnergy () const

const G4String &	GetName () const

virtual G4int	GetVerboseLevel () const

G4bool	IsApplicable (const G4DynamicParticle , const G4Element )

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)

virtual void	SetVerboseLevel (G4int value)

	~G4ParticleHPThermalScatteringData ()

Protected Attributes
G4String	name

G4int	verboseLevel

Private Member Functions
void	clearCurrentXSData ()

G4int	getTS_ID (const G4Material , const G4Element )

G4double	GetX (const G4DynamicParticle , G4double aT, std::map< G4double, G4ParticleHPVector > *)

std::map< G4double, G4ParticleHPVector * > *	readData (G4String)

Private Attributes
std::map< G4int, std::map< G4double, G4ParticleHPVector * > * > *	coherent

std::map< std::pair< const G4Material , const G4Element >, G4int >	dic

const G4Element *	element_cache

G4double	emax

std::map< G4int, std::map< G4double, G4ParticleHPVector * > * > *	incoherent

std::vector< G4int >	indexOfThermalElement

std::map< G4int, std::map< G4double, G4ParticleHPVector * > * > *	inelastic

G4bool	isForAllAtomsAndEnergies

G4double	ke_cache

const G4Material *	material_cache

G4double	maxKinEnergy

G4double	minKinEnergy

G4ParticleHPThermalScatteringNames *	names

G4CrossSectionDataSetRegistry *	registry

G4double	xs_cache

Detailed Description

Definition at line 61 of file G4ParticleHPThermalScatteringData.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPThermalScatteringData()

G4ParticleHPThermalScatteringData::G4ParticleHPThermalScatteringData ( )

Definition at line 56 of file G4ParticleHPThermalScatteringData.cc.

:G4VCrossSectionDataSet("NeutronHPThermalScatteringData")
,coherent(NULL)
,incoherent(NULL)
,inelastic(NULL)
{
// Upper limit of neutron energy 
   emax = 4*eV;
   SetMinKinEnergy( 0*MeV );                                   
   SetMaxKinEnergy( emax );                                   
 
   ke_cache = 0.0;
   xs_cache = 0.0;
   element_cache = NULL;
   material_cache = NULL;
 
   indexOfThermalElement.clear(); 
 
   names = new G4ParticleHPThermalScatteringNames();
}

References element_cache, emax, eV, indexOfThermalElement, ke_cache, material_cache, MeV, names, G4VCrossSectionDataSet::SetMaxKinEnergy(), G4VCrossSectionDataSet::SetMinKinEnergy(), and xs_cache.

◆ ~G4ParticleHPThermalScatteringData()

G4ParticleHPThermalScatteringData::~G4ParticleHPThermalScatteringData ( )

Definition at line 77 of file G4ParticleHPThermalScatteringData.cc.

{
 
   clearCurrentXSData();
 
   delete names;
}

References clearCurrentXSData(), and names.

Member Function Documentation

◆ AddUserThermalScatteringFile()

void G4ParticleHPThermalScatteringData::AddUserThermalScatteringFile	(	G4String	nameG4Element,
		G4String	filename
	)

Definition at line 568 of file G4ParticleHPThermalScatteringData.cc.

{
   names->AddThermalElement( nameG4Element , filename );
}

References G4ParticleHPThermalScatteringNames::AddThermalElement(), and names.

Referenced by G4ParticleHPThermalScattering::AddUserThermalScatteringFile().

◆ BuildPhysicsTable()

void G4ParticleHPThermalScatteringData::BuildPhysicsTable ( const G4ParticleDefinition & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 215 of file G4ParticleHPThermalScatteringData.cc.

{
 
   if ( &aP != G4Neutron::Neutron() ) 
      throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
 
   //std::map < std::pair < G4Material* , const G4Element* > , G4int > dic;   
   //
   dic.clear();   
   if ( G4Threading::IsMasterThread() ) clearCurrentXSData();
 
   std::map < G4String , G4int > co_dic;   
 
   //Searching Nist Materials
   static G4ThreadLocal G4MaterialTable* theMaterialTable  = 0 ; if (!theMaterialTable) theMaterialTable= G4Material::GetMaterialTable();
   size_t numberOfMaterials = G4Material::GetNumberOfMaterials();
   for ( size_t i = 0 ; i < numberOfMaterials ; i++ )
   {
      G4Material* material = (*theMaterialTable)[i];
      size_t numberOfElements = material->GetNumberOfElements();
      for ( size_t j = 0 ; j < numberOfElements ; j++ )
      {
         const G4Element* element = material->GetElement(j);
         if ( names->IsThisThermalElement ( material->GetName() , element->GetName() ) )
         {                                    
            G4int ts_ID_of_this_geometry; 
            G4String ts_ndl_name = names->GetTS_NDL_Name( material->GetName() , element->GetName() ); 
            if ( co_dic.find ( ts_ndl_name ) != co_dic.end() )
            {
               ts_ID_of_this_geometry = co_dic.find ( ts_ndl_name ) -> second;
            }
            else
            {
               ts_ID_of_this_geometry = co_dic.size();
               co_dic.insert ( std::pair< G4String , G4int >( ts_ndl_name , ts_ID_of_this_geometry ) );
            }
 
            //G4cout << "Neutron HP Thermal Scattering Data : Registering a material-element pair of " 
            //       << material->GetName() << " " << element->GetName() 
            //       << " as internal thermal scattering id of  " <<  ts_ID_of_this_geometry << "." << G4endl;
 
            dic.insert( std::pair < std::pair < G4Material* , const G4Element* > , G4int > ( std::pair < G4Material* , const G4Element* > ( material , element ) , ts_ID_of_this_geometry ) );
         }
      }
   }
 
   //Searching TS Elements 
   static G4ThreadLocal G4ElementTable* theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
   size_t numberOfElements = G4Element::GetNumberOfElements();
   //size_t numberOfThermalElements = 0; 
   for ( size_t i = 0 ; i < numberOfElements ; i++ )
   {
      const G4Element* element = (*theElementTable)[i];
      if ( names->IsThisThermalElement ( element->GetName() ) )
      {
         if ( names->IsThisThermalElement ( element->GetName() ) )
         {                                    
            G4int ts_ID_of_this_geometry; 
            G4String ts_ndl_name = names->GetTS_NDL_Name( element->GetName() ); 
            if ( co_dic.find ( ts_ndl_name ) != co_dic.end() )
            {
               ts_ID_of_this_geometry = co_dic.find ( ts_ndl_name ) -> second;
            }
            else
            {
               ts_ID_of_this_geometry = co_dic.size();
               co_dic.insert ( std::pair< G4String , G4int >( ts_ndl_name , ts_ID_of_this_geometry ) );
            }
 
            //G4cout << "Neutron HP Thermal Scattering: Registering an element of " 
            //       << material->GetName() << " " << element->GetName() 
            //       << " as internal thermal scattering id of  " <<  ts_ID_of_this_geometry << "." << G4endl;
 
            dic.insert( std::pair < std::pair < const G4Material* , const G4Element* > , G4int > ( std::pair < const G4Material* , const G4Element* > ( (G4Material*)NULL , element ) ,  ts_ID_of_this_geometry ) );
         }
      }
   }
 
   G4cout << G4endl;
   G4cout << "Neutron HP Thermal Scattering Data: Following material-element pairs and/or elements are registered." << G4endl;
   for ( std::map < std::pair < const G4Material* , const G4Element* > , G4int >::iterator it = dic.begin() ; it != dic.end() ; it++ )   
   {
      if ( it->first.first != NULL ) 
      {
         G4cout << "Material " << it->first.first->GetName() << " - Element " << it->first.second->GetName() << ",  internal thermal scattering id " << it->second << G4endl;
      }
      else
      {
         G4cout << "Element " << it->first.second->GetName() << ",  internal thermal scattering id " << it->second << G4endl;
      }
   }
   G4cout << G4endl;
 
 
   //G4cout << "Neutron HP Thermal Scattering Data: Following NDL thermal scattering files are assigned to the internal thermal scattering id." << G4endl;
   //for ( std::map < G4String , G4int >::iterator it = co_dic.begin() ; it != co_dic.end() ; it++ )  
   //{
   //   G4cout << "NDL file name " << it->first << ", internal thermal scattering id " << it->second << G4endl;
   //}
 
   G4ParticleHPManager* hpmanager = G4ParticleHPManager::GetInstance();
 
   coherent = hpmanager->GetThermalScatteringCoherentCrossSections();
   incoherent = hpmanager->GetThermalScatteringIncoherentCrossSections();
   inelastic = hpmanager->GetThermalScatteringInelasticCrossSections();
 
   if ( G4Threading::IsMasterThread() ) {
 
      if ( coherent == NULL ) coherent = new std::map< G4int , std::map< G4double , G4ParticleHPVector* >* >;
      if ( incoherent == NULL ) incoherent = new std::map< G4int , std::map< G4double , G4ParticleHPVector* >* >;
      if ( inelastic == NULL ) inelastic = new std::map< G4int , std::map< G4double , G4ParticleHPVector* >* >;
 
 
      // Read Cross Section Data files
 
      G4String dirName;
      if ( !std::getenv( "G4NEUTRONHPDATA" ) ) 
         throw G4HadronicException(__FILE__, __LINE__, "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
      G4String baseName = std::getenv( "G4NEUTRONHPDATA" );
 
      dirName = baseName + "/ThermalScattering";
 
      G4String ndl_filename;
      G4String full_name;
 
      for ( std::map < G4String , G4int >::iterator it = co_dic.begin() ; it != co_dic.end() ; it++ )  
      {
 
         ndl_filename = it->first;
         G4int ts_ID = it->second;
 
         // Coherent
         full_name = dirName + "/Coherent/CrossSection/" + ndl_filename; 
         std::map< G4double , G4ParticleHPVector* >*  coh_amapTemp_EnergyCross = readData( full_name );
         coherent->insert ( std::pair < G4int , std::map< G4double , G4ParticleHPVector* >* > ( ts_ID , coh_amapTemp_EnergyCross ) );
 
         // Incoherent
         full_name = dirName + "/Incoherent/CrossSection/" + ndl_filename; 
         std::map< G4double , G4ParticleHPVector* >*  incoh_amapTemp_EnergyCross = readData( full_name );
         incoherent->insert ( std::pair < G4int , std::map< G4double , G4ParticleHPVector* >* > ( ts_ID , incoh_amapTemp_EnergyCross ) );
 
         // Inelastic
         full_name = dirName + "/Inelastic/CrossSection/" + ndl_filename; 
         std::map< G4double , G4ParticleHPVector* >*  inela_amapTemp_EnergyCross = readData( full_name );
         inelastic->insert ( std::pair < G4int , std::map< G4double , G4ParticleHPVector* >* > ( ts_ID , inela_amapTemp_EnergyCross ) );
 
      }
      hpmanager->RegisterThermalScatteringCoherentCrossSections( coherent );
      hpmanager->RegisterThermalScatteringIncoherentCrossSections( incoherent );
      hpmanager->RegisterThermalScatteringInelasticCrossSections( inelastic );
   } 
}

References clearCurrentXSData(), coherent, dic, G4cout, G4endl, G4ThreadLocal, G4Element::GetElementTable(), G4ParticleHPManager::GetInstance(), G4Material::GetMaterialTable(), G4Element::GetName(), G4Element::GetNumberOfElements(), G4Material::GetNumberOfMaterials(), G4ParticleHPManager::GetThermalScatteringCoherentCrossSections(), G4ParticleHPManager::GetThermalScatteringIncoherentCrossSections(), G4ParticleHPManager::GetThermalScatteringInelasticCrossSections(), G4ParticleHPThermalScatteringNames::GetTS_NDL_Name(), incoherent, inelastic, G4Threading::IsMasterThread(), G4ParticleHPThermalScatteringNames::IsThisThermalElement(), anonymous_namespace{G4QuasiElRatios.cc}::map, eplot::material, names, G4Neutron::Neutron(), readData(), G4ParticleHPManager::RegisterThermalScatteringCoherentCrossSections(), G4ParticleHPManager::RegisterThermalScatteringIncoherentCrossSections(), G4ParticleHPManager::RegisterThermalScatteringInelasticCrossSections(), and second.

Referenced by G4ParticleHPThermalScattering::buildPhysicsTable().

◆ clearCurrentXSData()

void G4ParticleHPThermalScatteringData::clearCurrentXSData ( )

private

Definition at line 128 of file G4ParticleHPThermalScatteringData.cc.

{
   std::map< G4int , std::map< G4double , G4ParticleHPVector* >* >::iterator it;
   std::map< G4double , G4ParticleHPVector* >::iterator itt;
 
   if ( coherent != NULL ) {
   for ( it = coherent->begin() ; it != coherent->end() ; it++ )
      {
         if ( it->second != NULL )
         {
            for ( itt = it->second->begin() ; itt != it->second->end() ; itt++ )
            {
               delete itt->second;
            }
         }
         delete it->second;
      }
      coherent->clear();
   }
 
   if ( incoherent != NULL ) {
      for ( it = incoherent->begin() ; it != incoherent->end() ; it++ )
      {
         if ( it->second != NULL )
         {
            for ( itt = it->second->begin() ; itt != it->second->end() ; itt++ )
            {
               delete itt->second;
            }
         }
         delete it->second;
      }
      incoherent->clear();
   }
 
   if ( inelastic != NULL ) {
      for ( it = inelastic->begin() ; it != inelastic->end() ; it++ )
      {
         if ( it->second != NULL )
         {
            for ( itt = it->second->begin() ; itt != it->second->end() ; itt++ )
            {
               delete itt->second;
            }
         }
         delete it->second;
      }
      inelastic->clear();
   }
 
}

References coherent, incoherent, and inelastic.

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

◆ 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 G4ParticleHPThermalScatteringData::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 572 of file G4ParticleHPThermalScatteringData.cc.

{
    outFile << "High Precision cross data based on thermal scattering data in evaluated nuclear data libraries for neutrons below 5eV on specific materials\n" ;
}

◆ DumpPhysicsTable()

void G4ParticleHPThermalScatteringData::DumpPhysicsTable ( const G4ParticleDefinition & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 401 of file G4ParticleHPThermalScatteringData.cc.

{
   if( &aP != G4Neutron::Neutron() ) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
//  G4cout << "G4ParticleHPThermalScatteringData::DumpPhysicsTable still to be implemented"<<G4endl;
}

References G4Neutron::Neutron().

◆ ForAllAtomsAndEnergies()

bool G4VCrossSectionDataSet::ForAllAtomsAndEnergies ( ) const

inlineinherited

Definition at line 230 of file G4VCrossSectionDataSet.hh.

{
  return isForAllAtomsAndEnergies;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4CrossSectionDataStore::AddDataSet().

◆ GetCoherentCrossSection()

G4double G4ParticleHPThermalScatteringData::GetCoherentCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		const G4Material *	aM
	)

Definition at line 481 of file G4ParticleHPThermalScatteringData.cc.

{
   G4double result = 0;
   G4int ts_id = getTS_ID( aM , anE );
   G4double aT = aM->GetTemperature();
   result = GetX ( aP , aT , coherent->find( ts_id )->second );
   return result;
}

References coherent, G4Material::GetTemperature(), getTS_ID(), and GetX().

Referenced by G4ParticleHPThermalScattering::ApplyYourself().

◆ GetCrossSection()

G4double G4ParticleHPThermalScatteringData::GetCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		const G4Material *	aM
	)

Definition at line 450 of file G4ParticleHPThermalScatteringData.cc.

{
   G4double result = 0;
   
   G4int ts_id =getTS_ID( aM , anE );
 
   if ( ts_id == -1 ) return result;
 
   G4double aT = aM->GetTemperature();
 
   G4double Xcoh = GetX ( aP , aT , coherent->find(ts_id)->second );
   G4double Xincoh = GetX ( aP , aT , incoherent->find(ts_id)->second );
   G4double Xinela = GetX ( aP , aT , inelastic->find(ts_id)->second );
 
   result = Xcoh + Xincoh + Xinela;
 
   //G4cout << "G4ParticleHPThermalScatteringData::GetCrossSection  Tot= " << result/barn << " Coherent= " << Xcoh/barn << " Incoherent= " << Xincoh/barn << " Inelastic= " << Xinela/barn << G4endl;
 
   return result;
}

References coherent, G4Material::GetTemperature(), getTS_ID(), GetX(), incoherent, and inelastic.

Referenced by G4ParticleHPThermalScattering::ApplyYourself(), and 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().

◆ GetIncoherentCrossSection()

G4double G4ParticleHPThermalScatteringData::GetIncoherentCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		const G4Material *	aM
	)

Definition at line 490 of file G4ParticleHPThermalScatteringData.cc.

{
   G4double result = 0;
   G4int ts_id = getTS_ID( aM , anE );
   G4double aT = aM->GetTemperature();
   result = GetX ( aP , aT , incoherent->find( ts_id )->second );
   return result;
}

References G4Material::GetTemperature(), getTS_ID(), GetX(), and incoherent.

◆ GetInelasticCrossSection()

G4double G4ParticleHPThermalScatteringData::GetInelasticCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		const G4Material *	aM
	)

Definition at line 472 of file G4ParticleHPThermalScatteringData.cc.

{
   G4double result = 0;
   G4int ts_id = getTS_ID( aM , anE );
   G4double aT = aM->GetTemperature();
   result = GetX ( aP , aT , inelastic->find( ts_id )->second );
   return result;
}

References G4Material::GetTemperature(), getTS_ID(), GetX(), and inelastic.

Referenced by G4ParticleHPThermalScattering::ApplyYourself().

◆ GetIsoCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 110 of file G4ParticleHPThermalScatteringData.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() );
   G4double xs = GetCrossSection( dp , element , material );
   xs_cache = xs;
   return xs;
   //return GetCrossSection( dp , element , material->GetTemperature() );
}

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 G4ParticleHPInelasticData::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 G4ParticleHPInelasticData::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().

◆ getTS_ID()

G4int G4ParticleHPThermalScatteringData::getTS_ID	(	const G4Material *	material,
		const G4Element *	element
	)

private

Definition at line 501 of file G4ParticleHPThermalScatteringData.cc.

{
   G4int result = -1;
   if ( dic.find( std::pair < const G4Material* , const G4Element* > ( (G4Material*)NULL , element ) ) != dic.end() ) 
      return dic.find( std::pair < const G4Material* , const G4Element* > ( (G4Material*)NULL , element ) )->second; 
   if ( dic.find( std::pair < const G4Material* , const G4Element* > ( material , element ) ) != dic.end() ) 
      return dic.find( std::pair < const G4Material* , const G4Element* > ( material , element ) )->second; 
   return result; 
}

References dic, and eplot::material.

Referenced by GetCoherentCrossSection(), GetCrossSection(), GetIncoherentCrossSection(), and GetInelasticCrossSection().

◆ GetVerboseLevel()

G4int G4VCrossSectionDataSet::GetVerboseLevel ( ) const

inlinevirtualinherited

Reimplemented in G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, and G4ParticleHPInelasticData.

Definition at line 195 of file G4VCrossSectionDataSet.hh.

{
  return verboseLevel;
}

References G4VCrossSectionDataSet::verboseLevel.

◆ GetX()

G4double G4ParticleHPThermalScatteringData::GetX	(	const G4DynamicParticle *	aP,
		G4double	aT,
		std::map< G4double, G4ParticleHPVector * > *	amapTemp_EnergyCross
	)

private

Definition at line 514 of file G4ParticleHPThermalScatteringData.cc.

{
 
   G4double result = 0;
   if ( amapTemp_EnergyCross->size() == 0 ) return result;
 
 
   G4double eKinetic = aP->GetKineticEnergy();
 
   if ( amapTemp_EnergyCross->size() == 1 ) { 
      if ( std::fabs ( aT - amapTemp_EnergyCross->begin()->first ) / amapTemp_EnergyCross->begin()->first > 0.1 ) {
         G4cout << "G4ParticleHPThermalScatteringData:: The temperature of material (" 
                << aT/kelvin << "K) is different more than 10% from temperature of thermal scattering file expected (" 
                << amapTemp_EnergyCross->begin()->first << "K). Result may not be reliable."
         << G4endl;
      }
      result = amapTemp_EnergyCross->begin()->second->GetXsec ( eKinetic ); 
      return result;
   }
 
   std::map< G4double , G4ParticleHPVector* >::iterator it; 
   for ( it = amapTemp_EnergyCross->begin() ; it != amapTemp_EnergyCross->end() ; it++ ) {
       if ( aT < it->first ) break;
   } 
   //if ( it == amapTemp_EnergyCross->begin() && it != amapTemp_EnergyCross->end() ) it++;  // lower than the first
   //if ( it != amapTemp_EnergyCross->begin() && it == amapTemp_EnergyCross->end() ) it--;  // upper than the last
   if ( it == amapTemp_EnergyCross->begin() ) {
      it++;  // lower than the first
   } else if ( it == amapTemp_EnergyCross->end() ) {
      it--;  // upper than the last
   }
 
   G4double TH = it->first;
   G4double XH = it->second->GetXsec ( eKinetic ); 
 
   //G4cout << "G4ParticleHPThermalScatteringData::GetX TH " << TH << " E " << eKinetic <<  " XH " << XH << G4endl;
 
   if ( it != amapTemp_EnergyCross->begin() ) it--;
   G4double TL = it->first;
   G4double XL = it->second->GetXsec ( eKinetic ); 
 
   //G4cout << "G4ParticleHPThermalScatteringData::GetX TL " << TL << " E " << eKinetic <<  " XL " << XL << G4endl;
 
   if ( TH == TL )  
      throw G4HadronicException(__FILE__, __LINE__, "Thermal Scattering Data Error!");  
 
   G4double T = aT;
   G4double X = ( XH - XL ) / ( TH - TL ) * ( T - TL ) + XL;
   result = X;
  
   return result;
}

References G4cout, G4endl, G4DynamicParticle::GetKineticEnergy(), and kelvin.

Referenced by GetCoherentCrossSection(), GetCrossSection(), GetIncoherentCrossSection(), and GetInelasticCrossSection().

◆ IsApplicable()

G4bool G4ParticleHPThermalScatteringData::IsApplicable	(	const G4DynamicParticle *	aP,
		const G4Element *	anEle
	)

Definition at line 182 of file G4ParticleHPThermalScatteringData.cc.

{
   G4bool result = false;
 
   G4double eKin = aP->GetKineticEnergy();
   // Check energy 
   if ( eKin < emax )
   {
      // Check Particle Species
      if ( aP->GetDefinition() == G4Neutron::Neutron() ) 
      {
        // anEle is one of Thermal elements 
         G4int ie = (G4int) anEle->GetIndex();
         std::vector < G4int >::iterator it; 
         for ( it = indexOfThermalElement.begin() ; it != indexOfThermalElement.end() ; it++ )
         {
             if ( ie == *it ) return true;
         }
      }
   }
 
/*
   if ( names->IsThisThermalElement ( anEle->GetName() ) )
   {
      // Check energy and projectile species 
      G4double eKin = aP->GetKineticEnergy();
      if ( eKin < emax && aP->GetDefinition() == G4Neutron::Neutron() ) result = true; 
   }
*/
   return result;
}

References emax, G4DynamicParticle::GetDefinition(), G4Element::GetIndex(), G4DynamicParticle::GetKineticEnergy(), indexOfThermalElement, and G4Neutron::Neutron().

◆ IsElementApplicable()

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

virtualinherited

◆ IsIsoApplicable()

G4bool G4ParticleHPThermalScatteringData::IsIsoApplicable	(	const G4DynamicParticle *	dp,
		G4int	,
		G4int	,
		const G4Element *	element,
		const G4Material *	material
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 85 of file G4ParticleHPThermalScatteringData.cc.

{
   G4double eKin = dp->GetKineticEnergy();
   if ( eKin > 4.0*eV //GetMaxKinEnergy() 
     || eKin < 0 //GetMinKinEnergy() 
     || dp->GetDefinition() != G4Neutron::Neutron() ) return false;                                   
 
   if ( dic.find( std::pair < const G4Material* , const G4Element* > ( (G4Material*)NULL , element ) ) != dic.end() 
     || dic.find( std::pair < const G4Material* , const G4Element* > ( material , element ) ) != dic.end() ) return true;
 
   return false;
 
//   return IsApplicable( dp , element );
/*
   G4double eKin = dp->GetKineticEnergy();
   if ( eKin > 4.0*eV //GetMaxKinEnergy() 
     || eKin < 0 //GetMinKinEnergy() 
     || dp->GetDefinition() != G4Neutron::Neutron() ) return false;                                   
   return true;
*/
}

References dic, eV, G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), eplot::material, and G4Neutron::Neutron().

◆ readData()

std::map< G4double, G4ParticleHPVector * > * G4ParticleHPThermalScatteringData::readData ( G4String full_name )

private

Definition at line 370 of file G4ParticleHPThermalScatteringData.cc.

{
 
   std::map< G4double , G4ParticleHPVector* >*  aData = new std::map< G4double , G4ParticleHPVector* >; 
   
   //std::ifstream theChannel( full_name.c_str() );
   std::istringstream theChannel;
   G4ParticleHPManager::GetInstance()->GetDataStream(full_name,theChannel);
 
   //G4cout << "G4ParticleHPThermalScatteringData " << name << G4endl;
 
   G4int dummy; 
   while ( theChannel >> dummy )   // MF // Loop checking, 11.05.2015, T. Koi
   {
      theChannel >> dummy;   // MT
      G4double temp; 
      theChannel >> temp;   
      G4ParticleHPVector* anEnergyCross = new G4ParticleHPVector;
      G4int nData;
      theChannel >> nData;
      anEnergyCross->Init ( theChannel , nData , eV , barn );
      aData->insert ( std::pair < G4double , G4ParticleHPVector* > ( temp , anEnergyCross ) );
   }
   //theChannel.close();
 
   return aData;
 
} 

References barn, eV, G4ParticleHPManager::GetDataStream(), G4ParticleHPManager::GetInstance(), and G4ParticleHPVector::Init().

Referenced by BuildPhysicsTable().

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