141 std::map < G4double , std::vector < E_isoAng* >* >::iterator itt;
142 for ( itt = it->second->begin() ; itt != it->second->end() ; itt++ )
144 std::vector< E_isoAng* >::iterator ittt;
145 for ( ittt = itt->second->begin(); ittt != itt->second->end() ; ittt++ )
156 for ( std::map <
G4int , std::map <
G4double , std::vector < std::pair< G4double , G4double >* >* >* >::iterator it =
coherentFSs->begin() ; it !=
coherentFSs->end() ; it++ )
158 std::map < G4double , std::vector < std::pair< G4double , G4double >* >* >::iterator itt;
159 for ( itt = it->second->begin() ; itt != it->second->end() ; itt++ )
161 std::vector < std::pair< G4double , G4double >* >::iterator ittt;
162 for ( ittt = itt->second->begin(); ittt != itt->second->end() ; ittt++ )
175 std::map < G4double , std::vector < E_P_E_isoAng* >* >::iterator itt;
176 for ( itt = it->second->begin() ; itt != it->second->end() ; itt++ )
178 std::vector < E_P_E_isoAng* >::iterator ittt;
179 for ( ittt = itt->second->begin(); ittt != itt->second->end() ; ittt++ )
181 std::vector < E_isoAng* >::iterator it4;
182 for ( it4 = (*ittt)->vE_isoAngle.begin() ; it4 != (*ittt)->vE_isoAngle.end() ; it4++ )
212 std::map < G4double , std::vector < std::pair< G4double , G4double >* >* >* aCoherentFSDATA =
new std::map < G4double , std::vector < std::pair< G4double , G4double >* >* >;
218 std::vector< G4double > vBraggE;
221 while ( theChannel >> dummy )
226 std::vector < std::pair< G4double , G4double >* >* anBragE_P =
new std::vector < std::pair< G4double , G4double >* >;
230 for (
G4int i = 0 ; i <
n ; i++ )
234 if ( aCoherentFSDATA->size() == 0 )
237 vBraggE.push_back( Ei );
244 anBragE_P->push_back (
new std::pair < G4double , G4double > ( Ei , Pi ) );
247 aCoherentFSDATA->insert ( std::pair <
G4double , std::vector < std::pair< G4double , G4double >* >* > ( temp , anBragE_P ) );
250 return aCoherentFSDATA;
257 std::map < G4double , std::vector < E_P_E_isoAng* >* >* anT_E_P_E_isoAng =
new std::map < G4double , std::vector < E_P_E_isoAng* >* >;
264 while ( theChannel >> dummy )
269 std::vector < E_P_E_isoAng* >* vE_P_E_isoAng =
new std::vector < E_P_E_isoAng* >;
272 for (
G4int i = 0 ; i <
n ; i++ )
276 anT_E_P_E_isoAng->insert ( std::pair <
G4double , std::vector < E_P_E_isoAng* >* > ( temp , vE_P_E_isoAng ) );
280 return anT_E_P_E_isoAng;
300 for (
G4int i = 0 ; i < aData->
n ; i++ )
307 anE_isoAng->
n = nl - 2;
308 anE_isoAng->
isoAngle.resize( anE_isoAng->
n );
310 aData->
prob.push_back( prob );
312 for (
G4int j = 0 ; j < anE_isoAng->
n ; j++ )
323 for (
G4int i = 0 ; i < aData->
n - 1 ; i++ )
328 total += ( ( aData->
prob[i] ) * dE );
339 std::map < G4double , std::vector < E_isoAng* >* >* T_E =
new std::map < G4double , std::vector < E_isoAng* >* >;
346 while ( theChannel >> dummy )
351 std::vector < E_isoAng* >* vE_isoAng =
new std::vector < E_isoAng* >;
354 for (
G4int i = 0 ; i <
n ; i++ )
356 T_E->insert ( std::pair <
G4double , std::vector < E_isoAng* >* > ( temp , vE_isoAng ) );
384 for (
G4int i = 0 ; i < aData->
n ; i++ )
410 G4bool findThermalElement =
false;
413 for (
G4int i = 0; i <
n ; i++ )
420 if (
getTS_ID( NULL , theElement ) != -1 )
422 ielement =
getTS_ID( NULL , theElement );
423 findThermalElement =
true;
426 else if (
getTS_ID( theMaterial , theElement ) != -1 )
428 ielement =
getTS_ID( theMaterial , theElement );
429 findThermalElement =
true;
435 if ( findThermalElement ==
true )
447 if ( random <= inelastic/total )
452 std::map < G4double , std::vector< E_P_E_isoAng* >* >::iterator it;
453 std::vector<G4double> v_temp;
455 for ( it =
inelasticFSs->find( ielement )->second->begin() ; it !=
inelasticFSs->find( ielement )->second->end() ; it++ )
457 v_temp.push_back( it->first );
461 std::pair < G4double , G4double > tempLH =
find_LH ( aTemp , &v_temp );
465 std::vector< E_P_E_isoAng* >* vNEP_EPM_TL = 0;
466 std::vector< E_P_E_isoAng* >* vNEP_EPM_TH = 0;
468 if ( tempLH.first != 0.0 && tempLH.second != 0.0 )
470 vNEP_EPM_TL =
inelasticFSs->find( ielement )->second->find ( tempLH.first/
kelvin )->second;
471 vNEP_EPM_TH =
inelasticFSs->find( ielement )->second->find ( tempLH.second/
kelvin )->second;
473 else if ( tempLH.first == 0.0 )
475 std::map < G4double , std::vector< E_P_E_isoAng* >* >::iterator itm;
477 vNEP_EPM_TL = itm->second;
479 vNEP_EPM_TH = itm->second;
480 tempLH.first = tempLH.second;
481 tempLH.second = itm->first;
483 else if ( tempLH.second == 0.0 )
485 std::map < G4double , std::vector< E_P_E_isoAng* >* >::iterator itm;
488 vNEP_EPM_TH = itm->second;
490 vNEP_EPM_TL = itm->second;
491 tempLH.second = tempLH.first;
492 tempLH.first = itm->first;
501 sE =
get_linear_interpolated ( aTemp , std::pair < G4double , G4double > ( tempLH.first , TL.first ) , std::pair < G4double , G4double > ( tempLH.second , TH.first ) );
505 if ( TL.second.n == TH.second.n )
508 anE_isoAng.
n = TL.second.n;
509 for (
G4int i=0 ; i < anE_isoAng.
n ; i++ )
512 angle =
get_linear_interpolated ( aTemp , std::pair< G4double , G4double > ( tempLH.first , TL.second.isoAngle[ i ] ) , std::pair< G4double , G4double > ( tempLH.second , TH.second.isoAngle[ i ] ) );
513 anE_isoAng.
isoAngle.push_back( angle );
515 mu =
getMu( &anE_isoAng );
519 throw G4HadronicException(__FILE__, __LINE__,
"A problem is found in Thermal Scattering Data! Do not yet supported");
535 std::map < G4double , std::vector< std::pair< G4double , G4double >* >* >::iterator it;
536 std::vector<G4double> v_temp;
538 for ( it =
coherentFSs->find( ielement )->second->begin() ; it !=
coherentFSs->find( ielement )->second->end() ; it++ )
540 v_temp.push_back( it->first );
544 std::pair < G4double , G4double > tempLH =
find_LH ( aTemp , &v_temp );
549 std::vector< std::pair< G4double , G4double >* >* pvE_p_TL = NULL;
550 std::vector< std::pair< G4double , G4double >* >* pvE_p_TH = NULL;
552 if ( tempLH.first != 0.0 && tempLH.second != 0.0 )
554 pvE_p_TL =
coherentFSs->find( ielement )->second->find ( tempLH.first/
kelvin )->second;
555 pvE_p_TH =
coherentFSs->find( ielement )->second->find ( tempLH.first/
kelvin )->second;
557 else if ( tempLH.first == 0.0 )
559 pvE_p_TL =
coherentFSs->find( ielement )->second->find ( v_temp[ 0 ] )->second;
560 pvE_p_TH =
coherentFSs->find( ielement )->second->find ( v_temp[ 1 ] )->second;
561 tempLH.first = tempLH.second;
562 tempLH.second = v_temp[ 1 ];
564 else if ( tempLH.second == 0.0 )
566 pvE_p_TH =
coherentFSs->find( ielement )->second->find ( v_temp.back() )->
second;
567 std::vector< G4double >::iterator itv;
571 pvE_p_TL =
coherentFSs->find( ielement )->second->find ( *itv )->second;
572 tempLH.second = tempLH.first;
578 throw G4HadronicException(__FILE__, __LINE__,
"A problem is found in Thermal Scattering Data! Unexpected temperature values in data");
581 std::vector< G4double > vE_T;
582 std::vector< G4double > vp_T;
584 G4int n1 = pvE_p_TL->size();
588 for (
G4int i=0 ; i < n1 ; i++ )
590 if ( (*pvE_p_TL)[i]->
first != (*pvE_p_TH)[i]->
first )
throw G4HadronicException(__FILE__, __LINE__,
"A problem is found in Thermal Scattering Data!");
591 vE_T.push_back ( (*pvE_p_TL)[i]->
first );
592 vp_T.push_back (
get_linear_interpolated ( aTemp , std::pair< G4double , G4double > ( tempLH.first , (*pvE_p_TL)[i]->second ) , std::pair< G4double , G4double > ( tempLH.second , (*pvE_p_TL)[i]->second ) ) );
596 for (
G4int i = 1 ; i < n1 ; i++ )
598 if ( E/
eV < vE_T[ i ] )
608 for (
G4int i = 0 ; i <= j ; i++ )
610 G4double Pi = vp_T[ i ] / vp_T[ j ];
611 if ( rand_for_mu < Pi )
622 if ( mu < -1.0 ) mu = -1.0;
634 std::map < G4double , std::vector < E_isoAng* >* >::iterator it;
635 std::vector<G4double> v_temp;
639 v_temp.push_back( it->first );
643 std::pair < G4double , G4double > tempLH =
find_LH ( aTemp , &v_temp );
652 if ( tempLH.first != 0.0 && tempLH.second != 0.0 ) {
656 }
else if ( tempLH.first == 0.0 ) {
660 tempLH.first = tempLH.second;
661 tempLH.second = v_temp[ 1 ];
662 }
else if ( tempLH.second == 0.0 ) {
665 std::vector< G4double >::iterator itv;
670 tempLH.second = tempLH.first;
678 if ( anEPM_TL_E.
n == anEPM_TH_E.
n )
680 anEPM_T_E.
n = anEPM_TL_E.
n;
681 for (
G4int i=0 ; i < anEPM_TL_E.
n ; i++ )
684 angle =
get_linear_interpolated ( aTemp , std::pair< G4double , G4double > ( tempLH.first , anEPM_TL_E.
isoAngle[ i ] ) , std::pair< G4double , G4double > ( tempLH.second , anEPM_TH_E.
isoAngle[ i ] ) );
685 anEPM_T_E.
isoAngle.push_back( angle );
687 mu =
getMu ( &anEPM_T_E );
691 throw G4HadronicException(__FILE__, __LINE__,
"A problem is found in Thermal Scattering Data! Do not yet supported");
725 G4double mu_l = (*anEPM).isoAngle[ in-1 ];
727 result = ( mu_h - mu_l ) * ( random * ( (*anEPM).n ) -
in ) + mu_l;
738 G4double mu_h = (*anEPM).isoAngle[ 0 ];
739 result = ( mu_h - mu_l ) * xx + mu_l;
743 G4double mu_l = (*anEPM).isoAngle[ (*anEPM).n - 1 ];
745 result = ( mu_h - mu_l ) * xx + mu_l;
759 if ( aVector->size() == 1 ) {
760 LL = aVector->front();
761 H = aVector->front();
766 for ( std::vector< G4double >::iterator
767 it = aVector->begin() ; it != aVector->end() ; it++ ) {
770 if ( it != aVector->begin() ) {
782 if ( H == 0.0 ) LL = aVector->back();
785 return std::pair < G4double , G4double > (
LL , H );
793 if ( High.first - Low.first != 0 ) {
794 y = ( High.second - Low.second ) / ( High.first - Low.first ) * ( x - Low.first ) + Low.second;
796 if ( High.second == Low.second ) {
799 G4cout <<
"G4ParticleHPThermalScattering liner interpolation err!!" <<
G4endl;
812 std::vector< E_isoAng* >::iterator iv;
814 std::vector< G4double > v_e;
816 for ( iv = vEPM->begin() ; iv != vEPM->end() ; iv++ )
817 v_e.push_back ( (*iv)->energy );
819 std::pair < G4double , G4double > energyLH =
find_LH ( energy , &v_e );
825 if ( energyLH.first != 0.0 && energyLH.second != 0.0 )
827 for ( iv = vEPM->begin() ; iv != vEPM->end() ; iv++ )
829 if ( energyLH.first == (*iv)->energy ) {
837 else if ( energyLH.first == 0.0 )
839 panEPM_T_EL = (*vEPM)[0];
840 panEPM_T_EH = (*vEPM)[1];
842 else if ( energyLH.second == 0.0 )
844 panEPM_T_EH = (*vEPM).back();
853 if ( ! (
check_E_isoAng (panEPM_T_EL) ) ) panEPM_T_EL= panEPM_T_EH;
854 if ( ! (
check_E_isoAng (panEPM_T_EH) ) ) panEPM_T_EH= panEPM_T_EL;
856 if ( panEPM_T_EL->
n == panEPM_T_EH->
n )
859 anEPM_T_E.
n = panEPM_T_EL->
n;
861 for (
G4int i=0 ; i < panEPM_T_EL->
n ; i++ )
864 angle =
get_linear_interpolated ( energy , std::pair< G4double , G4double > ( energyLH.first , panEPM_T_EL->
isoAngle[ i ] ) , std::pair< G4double , G4double > ( energyLH.second , panEPM_T_EH->
isoAngle[ i ] ) );
865 anEPM_T_E.
isoAngle.push_back( angle );
870 G4cout <<
"G4ParticleHPThermalScattering Do not Suuport yet." <<
G4endl;
904 for (
G4int i = 0 ; i < n-1 ; i++ )
909 sum_p += ( ( anE_P_E_isoAng->
prob[i] ) * dE );
911 if ( random <= sum_p/total )
913 secondary_energy =
get_linear_interpolated ( random , std::pair < G4double , G4double > ( sum_p_L/total , E_L ) , std::pair < G4double , G4double > ( sum_p/total , E_H ) );
914 secondary_energy = secondary_energy*
eV;
920 return secondary_energy;
928 std::map< G4double , G4int > map_energy;
930 std::vector< G4double > v_energy;
932 std::vector< E_P_E_isoAng* >::iterator itv;
934 for ( itv = vNEP_EPM->begin(); itv != vNEP_EPM->end(); itv++ )
936 v_energy.push_back( (*itv)->energy );
937 map_energy.insert( std::pair < G4double , G4int > ( (*itv)->energy , i ) );
941 std::pair < G4double , G4double > energyLH =
find_LH ( pE , &v_energy );
946 if ( energyLH.first != 0.0 && energyLH.second != 0.0 )
948 pE_P_E_isoAng_EL = (*vNEP_EPM)[ map_energy.find ( energyLH.first )->second ];
949 pE_P_E_isoAng_EH = (*vNEP_EPM)[ map_energy.find ( energyLH.second )->second ];
951 else if ( energyLH.first == 0.0 )
953 pE_P_E_isoAng_EL = (*vNEP_EPM)[ 0 ];
954 pE_P_E_isoAng_EH = (*vNEP_EPM)[ 1 ];
956 if ( energyLH.second == 0.0 )
958 pE_P_E_isoAng_EH = (*vNEP_EPM).back();
959 itv = vNEP_EPM->end();
962 pE_P_E_isoAng_EL = *itv;
974 sE =
get_linear_interpolated ( pE , std::pair < G4double , G4double > ( energyLH.first , sE_L ) , std::pair < G4double , G4double > ( energyLH.second , sE_H ) );
984 if ( E_isoAng_L.
n == E_isoAng_H.
n )
986 anE_isoAng.
n = E_isoAng_L.
n;
987 for (
G4int j=0 ; j < anE_isoAng.
n ; j++ )
991 anE_isoAng.
isoAngle.push_back( angle );
1002 return std::pair< G4double , E_isoAng >( sE , anE_isoAng);
1017 std::map < G4String , G4int > co_dic;
1022 for (
size_t i = 0 ; i < numberOfMaterials ; i++ )
1024 G4Material* material = (*theMaterialTable)[i];
1026 for (
size_t j = 0 ; j < numberOfElements ; j++ )
1031 G4int ts_ID_of_this_geometry;
1033 if ( co_dic.find ( ts_ndl_name ) != co_dic.end() )
1035 ts_ID_of_this_geometry = co_dic.find ( ts_ndl_name ) ->
second;
1039 ts_ID_of_this_geometry = co_dic.size();
1040 co_dic.insert ( std::pair< G4String , G4int >( ts_ndl_name , ts_ID_of_this_geometry ) );
1047 dic.insert( std::pair < std::pair < G4Material* , const G4Element* > ,
G4int > ( std::pair < G4Material* , const G4Element* > ( material , element ) , ts_ID_of_this_geometry ) );
1056 for (
size_t i = 0 ; i < numberOfElements ; i++ )
1058 const G4Element* element = (*theElementTable)[i];
1063 G4int ts_ID_of_this_geometry;
1065 if ( co_dic.find ( ts_ndl_name ) != co_dic.end() )
1067 ts_ID_of_this_geometry = co_dic.find ( ts_ndl_name ) ->
second;
1071 ts_ID_of_this_geometry = co_dic.size();
1072 co_dic.insert ( std::pair< G4String , G4int >( ts_ndl_name , ts_ID_of_this_geometry ) );
1079 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 ) );
1085 G4cout <<
"Neutron HP Thermal Scattering: Following material-element pairs or elements are registered." <<
G4endl;
1086 for ( std::map < std::pair < const G4Material* , const G4Element* > ,
G4int >::iterator it =
dic.begin() ; it !=
dic.end() ; it++ )
1088 if ( it->first.first != NULL )
1090 G4cout <<
"Material " << it->first.first->GetName() <<
" - Element " << it->first.second->GetName() <<
", internal thermal scattering id " << it->second <<
G4endl;
1094 G4cout <<
"Element " << it->first.second->GetName() <<
", internal thermal scattering id " << it->second <<
G4endl;
1110 if (
coherentFSs == NULL )
coherentFSs =
new std::map < G4int , std::map < G4double , std::vector < std::pair< G4double , G4double >* >* >* >;
1112 if (
inelasticFSs == NULL )
inelasticFSs =
new std::map < G4int , std::map < G4double , std::vector < E_P_E_isoAng* >* >* >;
1115 if ( !getenv(
"G4NEUTRONHPDATA" ) )
1116 throw G4HadronicException(__FILE__, __LINE__,
"Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
1117 dirName = getenv(
"G4NEUTRONHPDATA" );
1121 for ( std::map < G4String , G4int >::iterator it = co_dic.begin() ; it != co_dic.end() ; it++ )
1124 G4int ts_ID = it->second;
1127 G4String fsName =
"/ThermalScattering/Coherent/FS/";
1128 G4String fileName = dirName + fsName + tsndlName;
1132 fsName =
"/ThermalScattering/Incoherent/FS/";
1133 fileName = dirName + fsName + tsndlName;
1137 fsName =
"/ThermalScattering/Inelastic/FS/";
1138 fileName = dirName + fsName + tsndlName;
1154 if (
dic.find( std::pair < const G4Material* , const G4Element* > ( material , element ) ) !=
dic.end() )
1155 result =
dic.find( std::pair < const G4Material* , const G4Element* > ( material , element ) )->second;
1179 for (
G4int i = 0 ; i <
n ; i++ ) {
1182 if ( sum != 0.0 ) result =
true;
1189 outFile <<
"High Precision model based on thermal scattering data in evaluated nuclear data libraries for neutrons below 5eV on specific materials\n";
static constexpr double kelvin
G4double GetCrossSection(const G4DynamicParticle *, const G4Element *, const G4Material *)
G4double total(Particle const *const p1, Particle const *const p2)
G4double get_secondary_energy_from_E_P_E_isoAng(G4double, E_P_E_isoAng *)
G4double GetCoherentCrossSection(const G4DynamicParticle *, const G4Element *, const G4Material *)
static G4ParticleHPManager * GetInstance()
void SetMinEnergy(G4double anEnergy)
void BuildPhysicsTable(const G4ParticleDefinition &)
void SetMaxEnergy(const G4double anEnergy)
static size_t GetNumberOfMaterials()
std::map< G4double, std::vector< E_isoAng * > * > * readAnIncoherentFSDATA(G4String)
static const G4int LL[nN]
static G4MaterialTable * GetMaterialTable()
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels() const
const G4String & GetName() const
void SetMomentumChange(const G4ThreeVector &aV)
G4double GetInelasticCrossSection(const G4DynamicParticle *, const G4Element *, const G4Material *)
G4double getMu(E_isoAng *)
std::map< G4double, std::vector< E_P_E_isoAng * > * > * readAnInelasticFSDATA(G4String)
void AddThermalElement(G4String, G4String)
static constexpr double perCent
void clearCurrentFSData()
std::map< G4int, std::map< G4double, std::vector< E_isoAng * > * > * > * incoherentFSs
std::map< G4int, std::map< G4double, std::vector< E_P_E_isoAng * > * > * > * GetThermalScatteringInelasticFinalStates()
std::vector< G4double > prob
static constexpr double second
std::vector< G4double > isoAngle
void RegisterThermalScatteringInelasticFinalStates(std::map< G4int, std::map< G4double, std::vector< E_P_E_isoAng * > * > * > *val)
void AddUserThermalScatteringFile(G4String, G4String)
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
std::pair< G4double, E_isoAng > create_sE_and_EPM_from_pE_and_vE_P_E_isoAng(G4double, G4double, std::vector< E_P_E_isoAng * > *)
const G4String & GetName() const
void SetEnergyChange(G4double anEnergy)
void GetDataStream(G4String, std::istringstream &iss)
const G4Material * GetMaterial() const
G4bool check_E_isoAng(E_isoAng *)
std::vector< G4Element * > G4ElementTable
G4bool IsThisThermalElement(G4String)
G4ParticleHPThermalScatteringNames names
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
void BuildPhysicsTable(const G4ParticleDefinition &)
std::map< std::pair< const G4Material *, const G4Element * >, G4int > dic
static constexpr double eV
void BuildPhysicsTable(const G4ParticleDefinition &)
void RegisterThermalScatteringCoherentFinalStates(std::map< G4int, std::map< G4double, std::vector< std::pair< G4double, G4double > * > * > * > *val)
std::map< G4int, std::map< G4double, std::vector< E_P_E_isoAng * > * > * > * inelasticFSs
E_P_E_isoAng * readAnE_P_E_isoAng(std::istream *)
virtual void ModelDescription(std::ostream &outFile) const
std::vector< E_isoAng * > vE_isoAngle
const G4LorentzVector & Get4Momentum() const
G4double G4ParticleHPJENDLHEData::G4double result
std::vector< G4Material * > G4MaterialTable
G4double GetKineticEnergy() const
G4ParticleHPElastic * theHPElastic
const G4ParticleDefinition * GetDefinition() const
std::map< G4double, std::vector< std::pair< G4double, G4double > * > * > * readACoherentFSDATA(G4String)
E_isoAng * readAnE_isoAng(std::istream *)
static G4Neutron * Neutron()
std::pair< G4double, G4double > find_LH(G4double, std::vector< G4double > *)
void RegisterThermalScatteringIncoherentFinalStates(std::map< G4int, std::map< G4double, std::vector< E_isoAng * > * > * > *val)
E_isoAng create_E_isoAng_from_energy(G4double, std::vector< E_isoAng * > *)
G4GLOB_DLL std::ostream G4cout
std::map< G4int, std::map< G4double, std::vector< std::pair< G4double, G4double > * > * > * > * coherentFSs
G4ParticleHPThermalScatteringData * theXSection
~G4ParticleHPThermalScattering()
G4int getTS_ID(const G4Material *, const G4Element *)
G4String GetTS_NDL_Name(G4String nameG4Element)
static size_t GetNumberOfElements()
G4ParticleHPThermalScattering()
const G4Element * GetElement(G4int iel) const
std::map< G4int, std::map< G4double, std::vector< E_isoAng * > * > * > * GetThermalScatteringIncoherentFinalStates()
void AddUserThermalScatteringFile(G4String, G4String)
G4HadFinalState theParticleChange
std::map< G4int, std::map< G4double, std::vector< std::pair< G4double, G4double > * > * > * > * GetThermalScatteringCoherentFinalStates()
size_t GetNumberOfElements() const
G4double GetTemperature() const
static G4ElementTable * GetElementTable()
G4double get_linear_interpolated(G4double, std::pair< G4double, G4double >, std::pair< G4double, G4double >)