G4ParticleHPJENDLHEElasticData Class Reference

#include <G4ParticleHPJENDLHEElasticData.hh>

Inheritance diagram for G4ParticleHPJENDLHEElasticData:

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
	G4ParticleHPJENDLHEElasticData ()
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)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
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)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
G4bool	IsZAApplicable (const G4DynamicParticle *, G4double, G4double)
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)

Protected Attributes
G4String	name
G4int	verboseLevel

Private Member Functions
G4double	getXSfromThisIsotope (G4int, G4int, G4double)
G4bool	isThisInMap (G4int, G4int)
G4bool	isThisNewIsotope (G4int z, G4int a)
G4PhysicsVector *	readAFile (std::fstream *)
void	registAPhysicsVector (G4int, G4int, G4PhysicsVector *)

Private Attributes
G4bool	isForAllAtomsAndEnergies
G4double	maxKinEnergy
G4double	minKinEnergy
std::map< G4int, std::map< G4int, G4PhysicsVector * > * >	mIsotope
G4String	particleName
G4String	reactionName
G4CrossSectionDataSetRegistry *	registry
std::vector< G4bool >	vElement

Detailed Description

Definition at line 41 of file G4ParticleHPJENDLHEElasticData.hh.

Constructor & Destructor Documentation

G4ParticleHPJENDLHEElasticData()

G4ParticleHPJENDLHEElasticData::G4ParticleHPJENDLHEElasticData

(

)

Definition at line 38 of file G4ParticleHPJENDLHEElasticData.cc.

:G4ParticleHPJENDLHEData( "Elastic" , G4Neutron::Neutron() ) 
{
   ;
}

Member Function Documentation

BuildPhysicsTable()

void G4ParticleHPJENDLHEData::BuildPhysicsTable ( const G4ParticleDefinition &  aP )

virtualinherited

Reimplemented from G4VCrossSectionDataSet.

Definition at line 100 of file G4ParticleHPJENDLHEData.cc.

{
 
//   if ( &aP != G4Neutron::Neutron() ) 
//      throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
   particleName = aP.GetParticleName();
 
   G4String baseName = std::getenv( "G4NEUTRONHPDATA" );
   G4String dirName = baseName+"/JENDL_HE/"+particleName+"/"+reactionName ;
   G4String aFSType = "/CrossSection/";
   G4ParticleHPNames theNames; 
 
   G4String filename;
 
// Create JENDL_HE data 
// Create map element or isotope  
 
   size_t numberOfElements = G4Element::GetNumberOfElements();
   //theCrossSections = new G4PhysicsTable( numberOfElements );
 
   // make a PhysicsVector for each element
 
   static G4ThreadLocal G4ElementTable *theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
   vElement.clear();
   vElement.resize( numberOfElements );
   for ( size_t i = 0; i < numberOfElements; ++i )
   {
 
      G4Element* theElement = (*theElementTable)[i];
      vElement[i] = false;
 
      // isotope
      G4int nIso = (*theElementTable)[i]->GetNumberOfIsotopes();
      G4int Z = static_cast<G4int> ((*theElementTable)[i]->GetZ());
      if ( nIso!=0 )
      {
         G4bool found_at_least_one = false; 
         for ( G4int i1 = 0; i1 < nIso; i1++ )
         {
             G4int A = theElement->GetIsotope(i1)->GetN();
 
             if ( isThisNewIsotope( Z , A ) ) 
             { 
 
                std::stringstream ss; 
                ss << dirName << aFSType << Z << "_" << A << "_" << theNames.GetName( Z-1 );
                filename = ss.str();
                std::fstream file;
                file.open ( filename , std::fstream::in );
                G4int dummy;
                file >> dummy;
                if ( file.good() ) 
                {
 
                   //G4cout << "Found file for Z=" << Z << ", A=" << A << ", as " << filename << G4endl;
                   found_at_least_one = true;
 
                   // read the file
                   G4PhysicsVector* aPhysVec = readAFile ( &file );
 
                   //Regist 
 
                   registAPhysicsVector( Z , A , aPhysVec );
 
                }
                else 
                { 
                   //G4cout << "No file for "<< reactionType << " Z=" << Z << ", A=" << A << G4endl;
                }
 
                file.close();
 
             }
             else
             {
                found_at_least_one = TRUE;
             }
          }
 
          if ( found_at_least_one ) vElement[i] = true;
 
       }
       else
       {
          G4StableIsotopes theStableOnes;
          G4int first = theStableOnes.GetFirstIsotope( Z );
          G4bool found_at_least_one = FALSE; 
          for ( G4int i1 = 0; i1 < theStableOnes.GetNumberOfIsotopes( static_cast<G4int>(theElement->GetZ() ) ); i1++)
          {
             G4int A = theStableOnes.GetIsotopeNucleonCount( first+i1 );
 
             if ( isThisNewIsotope( Z , A ) ) 
             {
 
                std::stringstream ss; 
                ss << dirName << aFSType << Z << "_" << A << "_" << theNames.GetName( Z-1 );
                filename = ss.str();
 
                std::fstream file;
                file.open ( filename , std::fstream::in );
                G4int dummy;
                file >> dummy;
                if ( file.good() ) 
                {
                   //G4cout << "Found file for Z=" << Z << ", A=" << A << ", as " << filename << G4endl;
                   found_at_least_one = TRUE;
                   //Read the file
 
                   G4PhysicsVector* aPhysVec = readAFile ( &file );
 
                   //Regist the PhysicsVector
                   registAPhysicsVector( Z , A , aPhysVec );
 
                }
                else 
                { 
                   //G4cout << "No file for "<< reactionType << " Z=" << Z << ", A=" << A << G4endl;
                }
 
                file.close();
             }
             else
             {
                found_at_least_one = TRUE;
             }
          }
 
          if ( found_at_least_one ) vElement[i] = true;
 
       }
 
   }
 
}

References A, FALSE, geant4_check_module_cycles::file, G4ThreadLocal, G4Element::GetElementTable(), G4StableIsotopes::GetFirstIsotope(), G4Element::GetIsotope(), G4StableIsotopes::GetIsotopeNucleonCount(), G4Isotope::GetN(), G4ParticleHPNames::GetName(), G4Element::GetNumberOfElements(), G4StableIsotopes::GetNumberOfIsotopes(), G4ParticleDefinition::GetParticleName(), G4Element::GetZ(), G4ParticleHPJENDLHEData::isThisNewIsotope(), G4ParticleHPJENDLHEData::particleName, G4ParticleHPJENDLHEData::reactionName, G4ParticleHPJENDLHEData::readAFile(), G4ParticleHPJENDLHEData::registAPhysicsVector(), TRUE, G4ParticleHPJENDLHEData::vElement, and Z.

Referenced by G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData().

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

virtualinherited

Reimplemented in G4ChipsAntiBaryonElasticXS, G4ChipsAntiBaryonInelasticXS, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonMinusElasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsKaonPlusInelasticXS, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4ChipsNeutronElasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4ChipsPionPlusElasticXS, G4ChipsPionPlusInelasticXS, G4ChipsProtonElasticXS, G4ChipsProtonInelasticXS, G4CrossSectionElastic, G4CrossSectionInelastic, G4ElectroNuclearCrossSection, G4ElNucleusSFcs, G4IonsShenCrossSection, G4KokoulinMuonNuclearXS, G4PhotoNuclearCrossSection, G4ZeroXS, G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, G4ParticleHPInelasticData, G4ParticleHPThermalScatteringData, G4BGGNucleonElasticXS, G4BGGPionElasticXS, G4BGGPionInelasticXS, G4GammaNuclearXS, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticXS, G4NucleonNuclearCrossSection, G4ParticleInelasticXS, G4UPiNuclearCrossSection, and G4BGGNucleonInelasticXS.

Definition at line 177 of file G4VCrossSectionDataSet.cc.

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

Referenced by G4CrossSectionDataStore::PrintCrossSectionHtml().

DumpPhysicsTable()

void G4ParticleHPJENDLHEData::DumpPhysicsTable ( const G4ParticleDefinition &  aP )

virtualinherited

Reimplemented from G4VCrossSectionDataSet.

Definition at line 237 of file G4ParticleHPJENDLHEData.cc.

{
  if(&aP!=G4Neutron::Neutron()) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
//  G4cout << "G4ParticleHPJENDLHEData::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().

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 G4ParticleHPJENDLHEData::GetCrossSection ( const G4DynamicParticle *  aP, const G4Element *  anE, G4double  aT  )

inherited

Definition at line 246 of file G4ParticleHPJENDLHEData.cc.

{
 
   // Primary energy >20MeV
   // Thus
   // Not take account of Doppler broadening 
   // also
   // Not take account of Target thermal motions
 
   G4double result = 0;
 
   G4double ek = aP->GetKineticEnergy();
 
   G4int nIso = anE->GetNumberOfIsotopes();
   G4int Z = static_cast<G4int> ( anE->GetZ() );
   if ( nIso!=0 )
   {
      for ( G4int i1 = 0; i1 < nIso; i1++ )
      {
 
         G4int A = anE->GetIsotope(i1)->GetN();
         G4double frac = anE->GetRelativeAbundanceVector()[ i1 ];   // This case do NOT request "*perCent".    
 
         result += frac * getXSfromThisIsotope( Z , A , ek );
         //G4cout << reactionType << " XS in barn " << Z << " " << A << " " << frac << " " << getXSfromThisIsotope( Z , A , ek )/barn << G4endl; 
 
      }
   }
   else
   {
 
      G4StableIsotopes theStableOnes;
      G4int first = theStableOnes.GetFirstIsotope( Z );
      for ( G4int i1 = 0; i1 < theStableOnes.GetNumberOfIsotopes( static_cast<G4int>(anE->GetZ() ) ); i1++)
      {
         
         G4int A = theStableOnes.GetIsotopeNucleonCount( first+i1 );
         G4double frac = theStableOnes.GetAbundance( first+i1 )*perCent;  // This case request "*perCent". 
 
         result += frac * getXSfromThisIsotope( Z , A , ek );
         //G4cout << reactionType << " XS in barn " << Z << " " << A << " " << frac << " " << getXSfromThisIsotope( Z , A , ek )/barn << G4endl; 
          
      }
   }
   return result;
 
}

References A, G4StableIsotopes::GetAbundance(), G4StableIsotopes::GetFirstIsotope(), G4Element::GetIsotope(), G4StableIsotopes::GetIsotopeNucleonCount(), G4DynamicParticle::GetKineticEnergy(), G4Isotope::GetN(), G4Element::GetNumberOfIsotopes(), G4StableIsotopes::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4ParticleHPJENDLHEData::getXSfromThisIsotope(), G4Element::GetZ(), perCent, and Z.

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 G4VCrossSectionDataSet::GetIsoCrossSection ( const G4DynamicParticle *  dynPart, G4int  Z, G4int  A, const G4Isotope *  iso = nullptr, const G4Element *  elm = nullptr, const G4Material *  mat = nullptr  )

virtualinherited

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

Definition at line 130 of file G4VCrossSectionDataSet.cc.

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

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

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

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

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.

getXSfromThisIsotope()

G4double G4ParticleHPJENDLHEData::getXSfromThisIsotope ( G4int  Z, G4int  A, G4double  ek  )

privateinherited

Definition at line 364 of file G4ParticleHPJENDLHEData.cc.

{
 
   G4double aXSection = 0.0;
   G4bool outOfRange;
 
   G4PhysicsVector* aPhysVec;
   if ( mIsotope.find ( Z )->second->find ( A ) != mIsotope.find ( Z )->second->end() )
   {
  
      aPhysVec = mIsotope.find ( Z )->second->find ( A )->second;
      aXSection = aPhysVec->GetValue( ek , outOfRange );
 
   } 
   else
   { 
 
      //Select closest one in the same Z
      std::map < G4int , G4PhysicsVector* >::iterator it; 
      G4int delta0 = 99; // no mean for 99 
      for ( it = mIsotope.find ( Z )->second->begin() ; it != mIsotope.find ( Z )->second->end() ; it++ )
      {
         G4int delta = std::abs( A - it->first );
         if ( delta < delta0 ) delta0 = delta;
      }
 
      // Randomize of selection larger or smaller than A
      if ( G4UniformRand() < 0.5 ) delta0 *= -1;
      G4int A1 = A + delta0;
      if ( mIsotope.find ( Z )->second->find ( A1 ) != mIsotope.find ( Z )->second->end() )
      {
         aPhysVec = mIsotope.find ( Z )->second->find ( A1 )->second;
      }
      else
      {
         A1 = A - delta0;
         aPhysVec = mIsotope.find ( Z )->second->find ( A1 )->second;
      }
 
      aXSection = aPhysVec->GetValue( ek , outOfRange );
      // X^(2/3) factor
      //aXSection *= std::pow ( 1.0*A/ A1 , 2.0 / 3.0 );
      aXSection *= G4Pow::GetInstance()->A23( 1.0*A/ A1 );
 
   }
 
   return aXSection;
}

References A, G4Pow::A23(), G4UniformRand, G4Pow::GetInstance(), G4PhysicsVector::GetValue(), G4ParticleHPJENDLHEData::mIsotope, and Z.

Referenced by G4ParticleHPJENDLHEData::GetCrossSection().

IsApplicable()

G4bool G4ParticleHPJENDLHEData::IsApplicable ( const G4DynamicParticle *  aP, const G4Element *  anE  )

inherited

Definition at line 41 of file G4ParticleHPJENDLHEData.cc.

{
 
   G4bool result = true;
   G4double eKin = aP->GetKineticEnergy();
   //if(eKin>20*MeV||aP->GetDefinition()!=G4Neutron::Neutron()) result = false;
   if ( eKin < 20*MeV || 3*GeV < eKin || aP->GetDefinition()!=G4Neutron::Neutron() ) 
   {
      result = false;
   } 
// Element Check 
   else if ( !(vElement[ anE->GetIndex() ]) ) result = false;
 
   return result;
 
}

References G4Element::GetIndex(), G4DynamicParticle::GetKineticEnergy(), MeV, G4Neutron::Neutron(), and G4ParticleHPJENDLHEData::vElement.

IsElementApplicable()

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

virtualinherited

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

Definition at line 64 of file G4VCrossSectionDataSet.cc.

{
  return false;
}

Referenced by G4VCrossSectionDataSet::ComputeCrossSection().

IsIsoApplicable()

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

virtualinherited

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

Definition at line 72 of file G4VCrossSectionDataSet.cc.

{
  return false;
}

Referenced by G4VCrossSectionDataSet::ComputeCrossSection().

isThisInMap()

G4bool G4ParticleHPJENDLHEData::isThisInMap ( G4int  z, G4int  a  )

privateinherited

Definition at line 333 of file G4ParticleHPJENDLHEData.cc.

{
   if ( mIsotope.find ( z ) == mIsotope.end() ) return false;
   if ( mIsotope.find ( z ) -> second->find ( a ) ==  mIsotope.find ( z ) -> second->end() ) return false;
   return true; 
}

References G4ParticleHPJENDLHEData::mIsotope, and second.

Referenced by G4ParticleHPJENDLHEData::isThisNewIsotope().

isThisNewIsotope()

G4bool G4ParticleHPJENDLHEData::isThisNewIsotope ( G4int  z, G4int  a  )

inlineprivateinherited

Definition at line 76 of file G4ParticleHPJENDLHEData.hh.

{ return !( isThisInMap( z , a ) ); };

References G4ParticleHPJENDLHEData::isThisInMap().

Referenced by G4ParticleHPJENDLHEData::BuildPhysicsTable().

IsZAApplicable()

G4bool G4ParticleHPJENDLHEData::IsZAApplicable ( const G4DynamicParticle *  , G4double  , G4double    )

inlineinherited

Definition at line 60 of file G4ParticleHPJENDLHEData.hh.

   { return false; }

readAFile()

G4PhysicsVector * G4ParticleHPJENDLHEData::readAFile ( std::fstream *  file )

privateinherited

Definition at line 298 of file G4ParticleHPJENDLHEData.cc.

{
 
   G4int dummy;
   G4int len;
   *file >> dummy;
   *file >> len;
 
   std::vector< G4double > v_e; 
   std::vector< G4double > v_xs; 
 
   for ( G4int i = 0 ; i < len ; i++ )
   {
      G4double e;
      G4double xs;
 
      *file >> e; 
      *file >> xs;
      // data are written in eV and barn.    
      v_e.push_back( e*eV );
      v_xs.push_back( xs*barn );
   }
 
   G4PhysicsFreeVector* aPhysVec = new G4PhysicsFreeVector( static_cast< size_t >( len ) , v_e.front() , v_e.back() );
 
   for ( G4int i = 0 ; i < len ; i++ )
   {
      aPhysVec->PutValues( static_cast< size_t >( i ) , v_e[ i ] , v_xs[ i ] );   
   }
 
   return aPhysVec;
}

References barn, eV, geant4_check_module_cycles::file, and G4PhysicsFreeVector::PutValues().

Referenced by G4ParticleHPJENDLHEData::BuildPhysicsTable().

registAPhysicsVector()

void G4ParticleHPJENDLHEData::registAPhysicsVector ( G4int  Z, G4int  A, G4PhysicsVector *  aPhysVec  )

privateinherited

Definition at line 342 of file G4ParticleHPJENDLHEData.cc.

{
 
    std::pair< G4int , G4PhysicsVector* > aPair = std::pair < G4int , G4PhysicsVector* > ( A , aPhysVec );  
 
    std::map < G4int , std::map< G4int , G4PhysicsVector* >* >::iterator itm; 
    itm = mIsotope.find ( Z );
    if ( itm !=  mIsotope.end() ) 
    { 
       itm->second->insert ( aPair ); 
    }  
    else
    {
       std::map< G4int , G4PhysicsVector* >* aMap = new std::map< G4int , G4PhysicsVector* >;
       aMap->insert ( aPair ); 
       mIsotope.insert( std::pair< G4int , std::map< G4int , G4PhysicsVector* >* > ( Z , aMap ) );
    }
 
}

References A, G4ParticleHPJENDLHEData::mIsotope, and Z.

Referenced by G4ParticleHPJENDLHEData::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;
}

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

SetVerboseLevel()

void G4VCrossSectionDataSet::SetVerboseLevel ( G4int  value )

inlinevirtualinherited

Reimplemented in G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, and G4ParticleHPInelasticData.

Definition at line 200 of file G4VCrossSectionDataSet.hh.

{
  verboseLevel = value;
}

References G4VCrossSectionDataSet::verboseLevel.

Referenced by G4VCrossSectionRatio::G4VCrossSectionRatio().

Field Documentation

isForAllAtomsAndEnergies

G4bool G4VCrossSectionDataSet::isForAllAtomsAndEnergies

privateinherited

Definition at line 182 of file G4VCrossSectionDataSet.hh.

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

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

mIsotope

std::map< G4int , std::map< G4int , G4PhysicsVector* >* > G4ParticleHPJENDLHEData::mIsotope

privateinherited

Definition at line 73 of file G4ParticleHPJENDLHEData.hh.

Referenced by G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData(), G4ParticleHPJENDLHEData::getXSfromThisIsotope(), G4ParticleHPJENDLHEData::isThisInMap(), and G4ParticleHPJENDLHEData::registAPhysicsVector().

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

particleName

G4String G4ParticleHPJENDLHEData::particleName

privateinherited

Definition at line 83 of file G4ParticleHPJENDLHEData.hh.

Referenced by G4ParticleHPJENDLHEData::BuildPhysicsTable().

reactionName

G4String G4ParticleHPJENDLHEData::reactionName

privateinherited

Definition at line 82 of file G4ParticleHPJENDLHEData.hh.

Referenced by G4ParticleHPJENDLHEData::BuildPhysicsTable(), and G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData().

registry

G4CrossSectionDataSetRegistry* G4VCrossSectionDataSet::registry

privateinherited

Definition at line 177 of file G4VCrossSectionDataSet.hh.

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

vElement

std::vector< G4bool > G4ParticleHPJENDLHEData::vElement

privateinherited

Definition at line 71 of file G4ParticleHPJENDLHEData.hh.

Referenced by G4ParticleHPJENDLHEData::BuildPhysicsTable(), and G4ParticleHPJENDLHEData::IsApplicable().

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

---

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

• source/processes/hadronic/models/particle_hp/include/G4ParticleHPJENDLHEElasticData.hh
• source/processes/hadronic/models/particle_hp/src/G4ParticleHPJENDLHEElasticData.cc