#include <G4LENDInelasticCrossSection.hh>

Inheritance diagram for G4LENDInelasticCrossSection:

Public Member Functions
void	AllowAnyCandidateTarget ()

void	AllowNaturalAbundanceTarget ()

void	BuildPhysicsTable (const G4ParticleDefinition &)

void	ChangeDefaultEvaluation (G4String name_tmp)

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

virtual void	CrossSectionDescription (std::ostream &) const

void	DumpLENDTargetInfo (G4bool force=false)

void	DumpPhysicsTable (const G4ParticleDefinition &)

bool	ForAllAtomsAndEnergies () const

	G4LENDInelasticCrossSection ()

	G4LENDInelasticCrossSection (G4ParticleDefinition *pd)

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

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

virtual 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)

virtual void	SetVerboseLevel (G4int value)

	~G4LENDInelasticCrossSection ()

Protected Member Functions
void	create_used_target_map ()

G4GIDI_target *	get_target_from_map (G4int nuclear_code)

G4double	GetUltraLowEnergyExtrapolatedXS (G4double, G4double, G4double, G4double, G4double)

Protected Attributes
G4String	name

G4ParticleDefinition *	proj

G4int	verboseLevel

Private Member Functions
G4double	getLENDCrossSection (G4GIDI_target *, G4double, G4double)

void	recreate_used_target_map ()

Private Attributes
G4bool	allow_any

G4bool	allow_nat

G4String	default_evaluation

G4bool	isForAllAtomsAndEnergies

G4LENDManager *	lend_manager

G4double	maxKinEnergy

G4double	minKinEnergy

G4CrossSectionDataSetRegistry *	registry

std::map< G4int, G4LENDUsedTarget * >	usedTarget_map

Detailed Description

Definition at line 45 of file G4LENDInelasticCrossSection.hh.

Constructor & Destructor Documentation

◆ G4LENDInelasticCrossSection() [1/2]

G4LENDInelasticCrossSection::G4LENDInelasticCrossSection ( )

inline

Definition at line 50 of file G4LENDInelasticCrossSection.hh.

51 {;};

◆ G4LENDInelasticCrossSection() [2/2]

G4LENDInelasticCrossSection::G4LENDInelasticCrossSection ( G4ParticleDefinition * pd )

inline

Definition at line 52 of file G4LENDInelasticCrossSection.hh.

      :G4LENDCrossSection("LENDInelasitcCrossSection")
      {
         proj = pd; 
         //name = "LEND Inelastic Cross Section for ";
         //name += proj->GetParticleName();
         //create_used_target_map();
      };

References G4LENDCrossSection::proj.

◆ ~G4LENDInelasticCrossSection()

G4LENDInelasticCrossSection::~G4LENDInelasticCrossSection ( )

inline

Definition at line 61 of file G4LENDInelasticCrossSection.hh.

61{;};

Member Function Documentation

◆ AllowAnyCandidateTarget()

void G4LENDCrossSection::AllowAnyCandidateTarget ( )

inlineinherited

Definition at line 96 of file G4LENDCrossSection.hh.

96{ allow_any = true; };

G4LENDCrossSection::allow_any

G4bool allow_any

Definition: G4LENDCrossSection.hh:120

References G4LENDCrossSection::allow_any.

◆ AllowNaturalAbundanceTarget()

void G4LENDCrossSection::AllowNaturalAbundanceTarget ( )

inlineinherited

Definition at line 95 of file G4LENDCrossSection.hh.

95{ allow_nat = true; };

G4LENDCrossSection::allow_nat

G4bool allow_nat

Definition: G4LENDCrossSection.hh:119

References G4LENDCrossSection::allow_nat.

Referenced by G4NeutronLENDBuilder::Build(), and G4HadronElasticPhysicsLEND::ConstructProcess().

◆ BuildPhysicsTable()

void G4LENDCrossSection::BuildPhysicsTable ( const G4ParticleDefinition & )

virtualinherited

Reimplemented from G4VCrossSectionDataSet.

Reimplemented in G4LENDGammaCrossSection.

Definition at line 148 of file G4LENDCrossSection.cc.

{
   create_used_target_map();
}

References G4LENDCrossSection::create_used_target_map().

Referenced by G4LENDCombinedCrossSection::BuildPhysicsTable(), and G4LENDGammaCrossSection::BuildPhysicsTable().

◆ ChangeDefaultEvaluation()

void G4LENDCrossSection::ChangeDefaultEvaluation ( G4String name_tmp )

inlineinherited

Definition at line 94 of file G4LENDCrossSection.hh.

94{ default_evaluation = name_tmp; };

G4LENDCrossSection::default_evaluation

G4String default_evaluation

Definition: G4LENDCrossSection.hh:118

References G4LENDCrossSection::default_evaluation.

Referenced by G4NeutronLENDBuilder::Build(), and G4HadronElasticPhysicsLEND::ConstructProcess().

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

◆ create_used_target_map()

void G4LENDCrossSection::create_used_target_map ( )

protectedinherited

Definition at line 310 of file G4LENDCrossSection.cc.

{
 
   lend_manager->RequestChangeOfVerboseLevel( verboseLevel );
 
   size_t numberOfElements = G4Element::GetNumberOfElements();
   static const G4ElementTable* theElementTable = G4Element::GetElementTable();
 
   for ( size_t i = 0 ; i < numberOfElements ; ++i )
   {
 
      const G4Element* anElement = (*theElementTable)[i];
      G4int numberOfIsotope = anElement->GetNumberOfIsotopes(); 
 
      if ( numberOfIsotope > 0 )
      {
      // User Defined Abundances   
         for ( G4int i_iso = 0 ; i_iso < numberOfIsotope ; i_iso++ )
         {
            G4int iZ = anElement->GetIsotope( i_iso )->GetZ();
            G4int iA = anElement->GetIsotope( i_iso )->GetN();
            G4int iIsomer = anElement->GetIsotope( i_iso )->Getm();
 
            //G4LENDUsedTarget* aTarget = new G4LENDUsedTarget ( G4Neutron::Neutron() , default_evaluation , iZ , iA );  
            G4LENDUsedTarget* aTarget = new G4LENDUsedTarget ( proj , default_evaluation , iZ , iA , iIsomer );  
            if ( allow_nat == true ) aTarget->AllowNat();
            if ( allow_any == true ) aTarget->AllowAny();
            usedTarget_map.insert( std::pair< G4int , G4LENDUsedTarget* > ( lend_manager->GetNucleusEncoding( iZ , iA , iIsomer ) , aTarget ) );
         }
      }
      else
      {
      // Natural Abundances   
         G4NistElementBuilder* nistElementBuild = lend_manager->GetNistElementBuilder();
         G4int iZ = int ( anElement->GetZ() );
         //G4cout << nistElementBuild->GetNumberOfNistIsotopes( int ( anElement->GetZ() ) ) << G4endl;
         G4int numberOfNistIso = nistElementBuild->GetNumberOfNistIsotopes( int ( anElement->GetZ() ) ); 
 
         for ( G4int ii = 0 ; ii < numberOfNistIso ; ii++ )
         {
            //G4cout << nistElementBuild->GetIsotopeAbundance( iZ , nistElementBuild->GetNistFirstIsotopeN( iZ ) + i ) << G4endl;
            if ( nistElementBuild->GetIsotopeAbundance( iZ , nistElementBuild->GetNistFirstIsotopeN( iZ ) + ii ) > 0 )
            {
               G4int iMass = nistElementBuild->GetNistFirstIsotopeN( iZ ) + ii;  
               //G4cout << iZ << " " << nistElementBuild->GetNistFirstIsotopeN( iZ ) + i << " " << nistElementBuild->GetIsotopeAbundance ( iZ , iMass ) << G4endl;  
               G4int iIsomer = 0; 
 
               G4LENDUsedTarget* aTarget = new G4LENDUsedTarget ( proj , default_evaluation , iZ , iMass );  
               if ( allow_nat == true ) aTarget->AllowNat();
               if ( allow_any == true ) aTarget->AllowAny();
               usedTarget_map.insert( std::pair< G4int , G4LENDUsedTarget* > ( lend_manager->GetNucleusEncoding( iZ , iMass , iIsomer ) , aTarget ) );
 
            }
 
         }
      }
   }
   DumpLENDTargetInfo();
}

References G4LENDCrossSection::allow_any, G4LENDCrossSection::allow_nat, G4LENDUsedTarget::AllowAny(), G4LENDUsedTarget::AllowNat(), G4LENDCrossSection::default_evaluation, G4LENDCrossSection::DumpLENDTargetInfo(), G4Element::GetElementTable(), G4Element::GetIsotope(), G4NistElementBuilder::GetIsotopeAbundance(), G4Isotope::Getm(), G4Isotope::GetN(), G4LENDManager::GetNistElementBuilder(), G4NistElementBuilder::GetNistFirstIsotopeN(), G4LENDManager::GetNucleusEncoding(), G4Element::GetNumberOfElements(), G4Element::GetNumberOfIsotopes(), G4NistElementBuilder::GetNumberOfNistIsotopes(), G4Element::GetZ(), G4Isotope::GetZ(), G4LENDCrossSection::lend_manager, G4LENDCrossSection::proj, G4LENDManager::RequestChangeOfVerboseLevel(), G4LENDCrossSection::usedTarget_map, and G4VCrossSectionDataSet::verboseLevel.

Referenced by G4LENDCombinedCrossSection::BuildPhysicsTable(), G4LENDCrossSection::BuildPhysicsTable(), G4LENDGammaCrossSection::BuildPhysicsTable(), G4LENDCrossSection::DumpLENDTargetInfo(), and G4LENDCrossSection::recreate_used_target_map().

◆ CrossSectionDescription()

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

virtualinherited

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

◆ DumpLENDTargetInfo()

void G4LENDCrossSection::DumpLENDTargetInfo ( G4bool force = false )

inherited

Definition at line 389 of file G4LENDCrossSection.cc.

                                                          {
 
   if ( lend_manager->GetVerboseLevel() >= 1 || force ) {
      if ( usedTarget_map.size() == 0 ) create_used_target_map(); 
      G4cout << "Dumping UsedTarget of " << GetName() << " for " << proj->GetParticleName() << G4endl;
      G4cout << "Requested Evaluation, Z , A -> Actual Evaluation, Z , A(0=Nat) " << G4endl;
      for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ ) {
         G4cout 
         << " " << it->second->GetWantedEvaluation() 
         << ", " << it->second->GetWantedZ() 
         << ", " << it->second->GetWantedA() 
         << " -> " << it->second->GetActualEvaluation() 
         << ", " << it->second->GetActualZ() 
         << ", " << it->second->GetActualA() 
         << G4endl; 
      } 
   }
}

References G4LENDCrossSection::create_used_target_map(), G4cout, G4endl, G4VCrossSectionDataSet::GetName(), G4ParticleDefinition::GetParticleName(), G4LENDManager::GetVerboseLevel(), G4LENDCrossSection::lend_manager, G4LENDCrossSection::proj, and G4LENDCrossSection::usedTarget_map.

Referenced by G4LENDCrossSection::create_used_target_map().

◆ DumpPhysicsTable()

void G4LENDCrossSection::DumpPhysicsTable ( const G4ParticleDefinition & aP )

virtualinherited

Reimplemented from G4VCrossSectionDataSet.

Definition at line 153 of file G4LENDCrossSection.cc.

{
 
  if ( &aP != proj ) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use LEND data for particles other than neutrons!!!");  
 
   G4cout << G4endl;
   G4cout << "Dump Cross Sections of " << GetName() << G4endl;
   G4cout << "(Pointwise cross-section at 300 Kelvin.)" << G4endl;
   G4cout << G4endl;
 
   G4cout << "Target informaiton " << G4endl;
 
   for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
   {
      G4cout 
         << "Wanted " << it->second->GetWantedEvaluation() 
         << ", Z= " << it->second->GetWantedZ() 
         << ", A= " << it->second->GetWantedA() 
         << "; Actual " << it->second->GetActualEvaluation() 
         << ", Z= " << it->second->GetActualZ() 
         << ", A= " << it->second->GetActualA() 
         << ", " << it->second->GetTarget() 
         << G4endl; 
 
      G4int ie = 0;
 
      G4GIDI_target* aTarget = it->second->GetTarget();
      G4double aT = 300;
      for ( ie = 0 ; ie < 130 ; ie++ )
      {
         G4double ke = 1.0e-5 * G4Pow::GetInstance()->powA ( 10.0 , ie/10.0 ) *eV;
 
         if ( ke < 20*MeV )
         {
            G4cout << "  " << GetName() << ", cross section at " << ke/eV << " [eV] = " << getLENDCrossSection ( aTarget , ke , aT )/barn << " [barn] " << G4endl;
         }
      }
      G4cout << G4endl;
 
   }
 
}

References barn, eV, G4cout, G4endl, G4Pow::GetInstance(), G4LENDCrossSection::getLENDCrossSection(), G4VCrossSectionDataSet::GetName(), MeV, G4Pow::powA(), G4LENDCrossSection::proj, and G4LENDCrossSection::usedTarget_map.

◆ ForAllAtomsAndEnergies()

bool G4VCrossSectionDataSet::ForAllAtomsAndEnergies ( ) const

inlineinherited

Definition at line 230 of file G4VCrossSectionDataSet.hh.

{
  return isForAllAtomsAndEnergies;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4CrossSectionDataStore::AddDataSet().

◆ get_target_from_map()

G4GIDI_target * G4LENDCrossSection::get_target_from_map ( G4int nuclear_code )

protectedinherited

Definition at line 381 of file G4LENDCrossSection.cc.

                                                                           {
   G4GIDI_target* target = NULL;
   if ( usedTarget_map.find( nuclear_code ) != usedTarget_map.end() ) {
      target = usedTarget_map.find( nuclear_code )->second->GetTarget();
   }
   return target;
}

References G4LENDCrossSection::usedTarget_map.

Referenced by G4LENDCrossSection::GetIsoCrossSection(), and G4LENDCrossSection::IsIsoApplicable().

◆ GetCrossSection()

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

◆ 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 G4LENDCrossSection::GetIsoCrossSection	(	const G4DynamicParticle *	dp,
		G4int	iZ,
		G4int	iA,
		const G4Isotope *	isotope,
		const G4Element *	,
		const G4Material *	material
	)

virtualinherited

Reimplemented from G4VCrossSectionDataSet.

Reimplemented in G4LENDGammaCrossSection.

Definition at line 83 of file G4LENDCrossSection.cc.

{
 
   G4double xs = 0.0;
   G4double ke = dp->GetKineticEnergy();
   G4double temp = material->GetTemperature();
   G4int iM = 0;
   if ( isotope != NULL ) iM = isotope->Getm();
 
   G4GIDI_target* aTarget = get_target_from_map( lend_manager->GetNucleusEncoding( iZ , iA , iM ) );
   if ( aTarget == NULL ) {
      G4String message;
      message = this->GetName();
      message += " is unexpectedly called.";
      //G4Exception( "G4LEND::GetIsoCrossSection(,)" , "LENDCrossSection-01" , JustWarning ,
      G4Exception( "G4LEND::GetIsoCrossSection(,)" , "LENDCrossSection-01" , FatalException ,
                  message );
   }
   xs = getLENDCrossSection ( aTarget , ke , temp );
 
   return xs;
}

References FatalException, G4Exception(), G4LENDCrossSection::get_target_from_map(), G4DynamicParticle::GetKineticEnergy(), G4LENDCrossSection::getLENDCrossSection(), G4Isotope::Getm(), G4VCrossSectionDataSet::GetName(), G4LENDManager::GetNucleusEncoding(), G4LENDCrossSection::lend_manager, and eplot::material.

Referenced by G4LENDCombinedCrossSection::GetIsoCrossSection(), G4LENDGammaCrossSection::GetIsoCrossSection(), G4LENDCombinedCrossSection::SelectChannel(), and G4LENDGammaCrossSection::SelectChannel().

◆ getLENDCrossSection()

G4double G4LENDInelasticCrossSection::getLENDCrossSection	(	G4GIDI_target *	target,
		G4double	ke,
		G4double	temperature
	)

privatevirtual

Reimplemented from G4LENDCrossSection.

Definition at line 30 of file G4LENDInelasticCrossSection.cc.

{
   if ( target == NULL ) return 0.0;
// 090407
//   return target->getElasticCrossSection( ke/MeV , temperature )*barn;
   //return target->getOthersCrossSectionAtE( ke/MeV , temperature )*barn;
   G4double result = target->getOthersCrossSectionAtE( ke/MeV , temperature )*barn;
   if ( result == 0.0 && ke/eV < 1.0e-4) 
   {
      G4double el = 1.0e-4*eV;
      G4double eh = 2.0e-4*eV;
      G4double xs_el = target->getOthersCrossSectionAtE( el/MeV , temperature )*barn;
      G4double xs_eh = target->getOthersCrossSectionAtE( eh/MeV , temperature )*barn;
      result = GetUltraLowEnergyExtrapolatedXS( el , eh , xs_el , xs_eh , ke );
   }
   return result;
}

References barn, eV, G4GIDI_target::getOthersCrossSectionAtE(), G4LENDCrossSection::GetUltraLowEnergyExtrapolatedXS(), and MeV.

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

◆ GetUltraLowEnergyExtrapolatedXS()

G4double G4LENDCrossSection::GetUltraLowEnergyExtrapolatedXS	(	G4double	x1,
		G4double	x2,
		G4double	y1,
		G4double	y2,
		G4double	ke
	)

protectedinherited

Definition at line 371 of file G4LENDCrossSection.cc.

{
   //XS propotinal to 1/v at low energy -> 1/root(E) 
   //XS = a * 1/root(E) + b  
   G4double a = ( y2 - y1 ) / ( 1/std::sqrt(x2) - 1/std::sqrt(x1) );
   G4double b = y1 - a * 1/std::sqrt(x1);
   G4double result = a * 1/std::sqrt(ke) + b;
   return result;
}

Referenced by G4LENDCaptureCrossSection::getLENDCrossSection(), G4LENDElasticCrossSection::getLENDCrossSection(), G4LENDFissionCrossSection::getLENDCrossSection(), and getLENDCrossSection().

◆ 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.

◆ IsElementApplicable()

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

virtualinherited

◆ IsIsoApplicable()

G4bool G4LENDCrossSection::IsIsoApplicable	(	const G4DynamicParticle *	dp,
		G4int	iZ,
		G4int	iA,
		const G4Element *	element,
		const G4Material *
	)

virtualinherited

Reimplemented from G4VCrossSectionDataSet.

Definition at line 46 of file G4LENDCrossSection.cc.

{
   G4double eKin = dp->GetKineticEnergy();
   if ( dp->GetDefinition() != proj ) return false;
   if ( eKin > GetMaxKinEnergy() || eKin < GetMinKinEnergy() ) return false;
 
   //G4cout << "G4LENDCrossSection::GetIsoIsIsoApplicable this->GetName() = " << this->GetName() << ", iZ = " << iZ << ", iA = " << iA << ", allow_nat = " << allow_nat << G4endl;
   //Check existence of target data
   if ( element != NULL ) { 
      if ( element->GetNumberOfIsotopes() != 0 ) { 
         std::vector< const G4Isotope*> vIsotope;
         for ( size_t i = 0 ; i != element->GetNumberOfIsotopes() ; i++ ) {
            if ( element->GetIsotope( i )->GetN() == iA ) vIsotope.push_back( element->GetIsotope( i ) ); 
         }
         for ( size_t i = 0 ; i != vIsotope.size() ; i++ ) { 
            G4int iM = vIsotope[i]->Getm(); 
            if ( get_target_from_map( lend_manager->GetNucleusEncoding( iZ , iA , iM ) ) != NULL ) return true;
         } 
         //No isomer has data
         //Check natural aboundance data for the element
         if ( get_target_from_map( lend_manager->GetNucleusEncoding( iZ , 0 , 0 ) ) != NULL ) return true;
      } else {
         //Check for iZ and iA under assuming iM = 0
         if ( get_target_from_map( lend_manager->GetNucleusEncoding( iZ , iA , 0 ) ) != NULL ) return true;
         //Check natural aboundance data for the element
         if ( get_target_from_map( lend_manager->GetNucleusEncoding( iZ , 0 , 0 ) ) != NULL ) return true;
      }
   } else {
      //Check for iZ and iA under assuming iM = 0
      if ( get_target_from_map( lend_manager->GetNucleusEncoding( iZ , iA , 0 ) ) != NULL ) return true;
      //Check natural aboundance data for iZ
      if ( get_target_from_map( lend_manager->GetNucleusEncoding( iZ , 0 , 0 ) ) != NULL ) return true;
   }
   return false;
}

References G4LENDCrossSection::get_target_from_map(), G4DynamicParticle::GetDefinition(), G4Element::GetIsotope(), G4DynamicParticle::GetKineticEnergy(), G4VCrossSectionDataSet::GetMaxKinEnergy(), G4VCrossSectionDataSet::GetMinKinEnergy(), G4Isotope::GetN(), G4LENDManager::GetNucleusEncoding(), G4Element::GetNumberOfIsotopes(), G4LENDCrossSection::lend_manager, and G4LENDCrossSection::proj.

◆ recreate_used_target_map()

void G4LENDCrossSection::recreate_used_target_map ( )

privateinherited

Definition at line 296 of file G4LENDCrossSection.cc.

{
   for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
   { 
      delete it->second;  
   }
   usedTarget_map.clear();
 
   create_used_target_map();
}

References G4LENDCrossSection::create_used_target_map(), and G4LENDCrossSection::usedTarget_map.

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