89 default_evaluation =
"ENDF.B-VII.0";
104 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
105 it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
121 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use LEND data for particles other than neutrons!!!");
125 G4cout <<
"(Pointwise cross-section at 300 Kelvin.)" <<
G4endl;
130 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
131 it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
134 <<
"Wanted " << it->second->GetWantedEvaluation()
135 <<
", Z= " << it->second->GetWantedZ()
136 <<
", A= " << it->second->GetWantedA()
137 <<
"; Actual " << it->second->GetActualEvaluation()
138 <<
", Z= " << it->second->GetActualZ()
139 <<
", A= " << it->second->GetActualA()
140 <<
", " << it->second->GetTarget()
147 for ( ie = 0 ; ie < 130 ; ie++ )
149 G4double ke = 1.0e-5 * std::pow ( 10.0 , ie/10.0 ) *
eV;
260 void G4LENDCrossSection::recreate_used_target_map()
262 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
263 it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
267 usedTarget_map.clear();
282 for (
size_t i = 0 ; i < numberOfElements ; ++i )
285 const G4Element* anElement = (*theElementTable)[i];
288 if ( numberOfIsotope > 0 )
291 for (
G4int i_iso = 0 ; i_iso < numberOfIsotope ; i_iso++ )
299 if ( allow_nat ==
true ) aTarget->
AllowNat();
300 if ( allow_any ==
true ) aTarget->
AllowAny();
301 usedTarget_map.insert( std::pair< G4int , G4LENDUsedTarget* > ( lend_manager->
GetNucleusEncoding( iZ , iA , iIsomer ) , aTarget ) );
312 for (
G4int ii = 0 ; ii < numberOfNistIso ; ii++ )
322 if ( allow_nat ==
true ) aTarget->
AllowNat();
323 if ( allow_any ==
true ) aTarget->
AllowAny();
324 usedTarget_map.insert( std::pair< G4int , G4LENDUsedTarget* > ( lend_manager->
GetNucleusEncoding( iZ , iMass , iIsomer ) , aTarget ) );
333 G4cout <<
"Requested Evaluation, Z , A -> Actual Evaluation, Z , A(0=Nat) , Pointer of Target" <<
G4endl;
334 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
335 it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
338 <<
" " << it->second->GetWantedEvaluation()
339 <<
", " << it->second->GetWantedZ()
340 <<
", " << it->second->GetWantedA()
341 <<
" -> " << it->second->GetActualEvaluation()
342 <<
", " << it->second->GetActualZ()
343 <<
", " << it->second->GetActualA()
344 <<
", " << it->second->GetTarget()
355 G4double a = ( y2 -
y1 ) / ( 1/std::sqrt(x2) - 1/std::sqrt(x1) );
size_t GetNumberOfIsotopes() const
virtual G4double getLENDCrossSection(G4GIDI_target *, G4double, G4double)
void DumpPhysicsTable(const G4ParticleDefinition &)
G4double GetUltraLowEnergyExtrapolatedXS(G4double, G4double, G4double, G4double, G4double)
G4double GetKineticEnergy() const
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
G4NistElementBuilder * GetNistElementBuilder()
G4double G4NeutronHPJENDLHEData::G4double result
G4double GetIsotopeAbundance(G4int Z, G4int N) const
G4int GetNumberOfNistIsotopes(G4int Z) const
G4ParticleDefinition * GetDefinition() const
G4ParticleDefinition * proj
const G4String & GetName() const
G4LENDCrossSection(const G4String name="")
G4bool IsIsoApplicable(const G4DynamicParticle *, G4int, G4int, const G4Element *, const G4Material *)
void create_used_target_map()
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
G4GLOB_DLL std::ostream G4cout
static size_t GetNumberOfElements()
void SetMinKinEnergy(G4double value)
G4int GetNistFirstIsotopeN(G4int Z) const
static G4LENDManager * GetInstance()
void SetMaxKinEnergy(G4double value)
const G4Isotope * GetIsotope(G4int iso) const
G4double GetTemperature() const
std::vector< G4Element * > G4ElementTable
G4bool RequestChangeOfVerboseLevel(G4int)
static G4ElementTable * GetElementTable()
G4double GetMinKinEnergy() const
G4double GetMaxKinEnergy() const
G4double GetIsoCrossSection(const G4DynamicParticle *, G4int, G4int, const G4Isotope *, const G4Element *, const G4Material *)
void BuildPhysicsTable(const G4ParticleDefinition &)