48 theCrossSectionHandler(0)
82 G4cout <<
"Entering in method G4PenelopeIonisationCrossSection::CrossSection()" <<
G4endl;
91 ed <<
"The method has been called with a null G4Material pointer" <<
G4endl;
92 G4Exception(
"G4PenelopeIonisationCrossSection::CrossSection()",
"em2042",
104 if(
G4int(shell) < nmax &&
128 ed <<
"There is something wrong here: it looks like the index is wrong" <<
G4endl;
129 ed <<
"Requested: shell " <<
G4int(shell) <<
" and Z = " << Z <<
G4endl;
131 G4Exception(
"G4PenelopeIonisationCrossSection::CrossSection()",
"em2043",
143 cross = crossPerMolecule/atomsPerMolec;
148 G4cout <<
"Cross section of shell " <<
G4int(shell) <<
" and Z= " <<
Z;
149 G4cout <<
" of material: " << material->
GetName() <<
" and energy = " << incidentEnergy/
keV <<
" keV" <<
G4endl;
155 G4cout <<
"Cross section per molecule: " << crossPerMolecule/
barn <<
" barn" <<
G4endl;
156 G4cout <<
"Atoms " << Z <<
" per molecule: " << atomsPerMolec <<
G4endl;
165 std::vector<G4double>
172 std::vector<G4double> vec(nmax,0.0);
181 std::vector<G4double>
189 size_t n = vec.size();
192 for(i=0; i<
n; ++i) { sum += vec[i]; }
195 for(i=0; i<
n; ++i) { vec[i] = vec[i]*
sum; }
207 G4cout <<
"Entering in method G4PenelopeIonisationCrossSection::FindShellIDIndex()" <<
G4endl;
212 std::pair<const G4Material*,G4int> theKey = std::make_pair(mat,Z);
220 G4DataVector* theVec =
shellIDTable->find(theKey)->second;
222 if (ishell>=0 && ishell < (
G4int) theVec->size())
223 result = (
G4int) (*theVec)[ishell];
227 ed <<
"Shell ID: " << ishell <<
" not available for material " << mat->
GetName() <<
" and Z = " <<
240 size_t numberOfOscillators = theTable->size();
244 G4DataVector* dat =
new G4DataVector(
nMaxLevels,-1);
245 for (
size_t iosc=0;iosc<numberOfOscillators;iosc++)
256 if ((shFlag-1) == ishell)
257 result = (
G4int) iosc;
264 G4cout <<
"Leaving method G4PenelopeIonisationCrossSection::FindShellIDIndex() with index = " << result <<
G4endl;
G4double fHighEnergyLimit
std::vector< G4PenelopeOscillator * > G4PenelopeOscillatorTable
G4PenelopeOscillator * GetOscillatorIonisation(const G4Material *, G4int)
std::ostringstream G4ExceptionDescription
G4int NumberOfShells(G4int Z) const
G4PenelopeIonisationXSHandler * theCrossSectionHandler
static constexpr double keV
G4double GetIonisationEnergy()
std::vector< G4double > Probabilities(G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
G4PenelopeOscillatorTable * GetOscillatorTableIonisation(const G4Material *)
G4double GetResonanceEnergy() const
std::vector< G4double > GetCrossSection(G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
G4double GetShellCrossSection(size_t shellID, G4double energy) const
Returns the hard cross section for the given shell (per molecule)
G4int FindShellIDIndex(const G4Material *mat, G4int Z, G4AtomicShellEnumerator shell)
const G4String & GetName() const
void BuildXSTable(const G4Material *, G4double cut, const G4ParticleDefinition *, G4bool isMaster=true)
This can be inkoved only by the master.
static G4PenelopeOscillatorManager * GetOscillatorManager()
G4PenelopeIonisationCrossSection()
Constructor.
static constexpr double eV
static G4Electron * Electron()
std::map< std::pair< const G4Material *, G4int >, G4DataVector * > * shellIDTable
G4double G4ParticleHPJENDLHEData::G4double result
G4double CrossSection(G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
const G4AtomicTransitionManager * transitionManager
static constexpr double barn
G4GLOB_DLL std::ostream G4cout
const G4PenelopeCrossSection * GetCrossSectionTableForCouple(const G4ParticleDefinition *, const G4Material *, const G4double cut) const
~G4PenelopeIonisationCrossSection()
Destructor. Clean all tables.
static G4AtomicTransitionManager * Instance()
G4PenelopeOscillatorManager * oscManager
G4double GetNumberOfZAtomsPerMolecule(const G4Material *, G4int Z)
static constexpr double GeV
T min(const T t1, const T t2)
brief Return the smallest of the two arguments