60 for (
G4int i=6; i<93; i++)
66 for (
G4int i=6; i<93; i++)
93 if (energyIncident > 0.2*
MeV && energyIncident < 5.*
MeV && zTarget < 93 && zTarget > 5) {
96 if (sigma !=0 && energyIncident >
protonDataSetMap[zTarget]->GetEnergies(0).back()*
MeV)
return 0.;
100 else if (massIncident == aAlpha->
GetPDGMass())
102 if (energyIncident > 0.2*
MeV && energyIncident < 40.*
MeV && zTarget < 93 && zTarget > 5) {
105 if (sigma !=0 && energyIncident >
alphaDataSetMap[zTarget]->GetEnergies(0).back()*
MeV)
return 0.;
static constexpr double MeV
G4GLOB_DLL std::ostream G4cout
virtual ~G4ANSTOecpssrKxsModel()
G4VDataSetAlgorithm * interpolation
std::map< G4int, G4VEMDataSet * > alphaDataSetMap
G4double CalculateCrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)
std::map< G4int, G4VEMDataSet * > protonDataSetMap
G4double GetPDGMass() const
static G4Proton * Proton()