Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PenelopeIonisationCrossSection Class Reference

#include <G4PenelopeIonisationCrossSection.hh>

Inheritance diagram for G4PenelopeIonisationCrossSection:
Collaboration diagram for G4PenelopeIonisationCrossSection:

Public Member Functions

 G4PenelopeIonisationCrossSection ()
 Constructor. More...
 
 ~G4PenelopeIonisationCrossSection ()
 Destructor. Clean all tables. More...
 
std::vector< G4doubleGetCrossSection (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
 
G4double CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat)
 
std::vector< G4doubleProbabilities (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
 
void SetVerbosityLevel (G4int vl)
 Getter/setter for the verbosity level. More...
 
G4int GetVerbosityLevel ()
 
- Public Member Functions inherited from G4VhShellCrossSection
 G4VhShellCrossSection (const G4String &xname="")
 
virtual ~G4VhShellCrossSection ()
 
G4int SelectRandomShell (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
 
virtual void SetTotalCS (G4double)
 
const G4StringGetName () const
 

Detailed Description

Definition at line 58 of file G4PenelopeIonisationCrossSection.hh.

Constructor & Destructor Documentation

G4PenelopeIonisationCrossSection::G4PenelopeIonisationCrossSection ( )

Constructor.

NOTICE: working only for e- at the moment (no interface available for e+)

Definition at line 46 of file G4PenelopeIonisationCrossSection.cc.

46  :
47  G4VhShellCrossSection("Penelope"),shellIDTable(0),
48  theCrossSectionHandler(0)
49 {
51  nMaxLevels = 9;
52 
53  // Verbosity scale:
54  // 0 = nothing
55  // 1 = calculation of cross sections, file openings, sampling of atoms
56  // 2 = entering in methods
57  verboseLevel = 0;
58 
59  fLowEnergyLimit = 10.0*eV;
60  fHighEnergyLimit = 100.0*GeV;
61 
62  transitionManager = G4AtomicTransitionManager::Instance();
63 }
G4VhShellCrossSection(const G4String &xname="")
static constexpr double eV
Definition: G4SIunits.hh:215
static G4PenelopeOscillatorManager * GetOscillatorManager()
static constexpr double GeV
Definition: G4SIunits.hh:217
static G4AtomicTransitionManager * Instance()

Here is the call graph for this function:

G4PenelopeIonisationCrossSection::~G4PenelopeIonisationCrossSection ( )

Destructor. Clean all tables.

Definition at line 67 of file G4PenelopeIonisationCrossSection.cc.

68 {
69  if (theCrossSectionHandler)
70  delete theCrossSectionHandler;
71 }

Member Function Documentation

G4double G4PenelopeIonisationCrossSection::CrossSection ( G4int  Z,
G4AtomicShellEnumerator  shell,
G4double  incidentEnergy,
G4double  mass,
const G4Material mat 
)
virtual

Purely virtual method from the base interface. Returns the cross section for the given shell in the element Z of material mat at the specified energy

Implements G4VhShellCrossSection.

Definition at line 75 of file G4PenelopeIonisationCrossSection.cc.

80 {
81  if (verboseLevel > 1)
82  G4cout << "Entering in method G4PenelopeIonisationCrossSection::CrossSection()" << G4endl;
83 
84  G4double cross = 0.;
85 
86  //Material pointer is not available
87  if (!material)
88  {
89  //CRASH!
91  ed << "The method has been called with a null G4Material pointer" << G4endl;
92  G4Exception("G4PenelopeIonisationCrossSection::CrossSection()","em2042",
93  FatalException,ed);
94  return cross;
95  }
96 
97  if (!theCrossSectionHandler)
98  theCrossSectionHandler = new G4PenelopeIonisationXSHandler();
99 
100  theCrossSectionHandler->BuildXSTable(material,0.,G4Electron::Electron());
101 
102  G4int nmax = std::min(nMaxLevels,transitionManager->NumberOfShells(Z));
103 
104  if(G4int(shell) < nmax &&
105  incidentEnergy >= fLowEnergyLimit && incidentEnergy <= fHighEnergyLimit)
106  {
107  //The shells in Penelope are organized per *material*, rather than per
108  //element, so given a material one has to find the proper index for the
109  //given Z and shellID. An appropriate lookup table is used to avoid
110  //recalculation.
111  G4int index = FindShellIDIndex(material,Z,shell);
112 
113  //Index is not available!
114  if (index < 0)
115  return cross;
116 
117  const G4PenelopeCrossSection* theXS =
118  theCrossSectionHandler->GetCrossSectionTableForCouple(G4Electron::Electron(),
119  material,
120  0.);
121 
122  //Cross check that everything is fine:
123  G4PenelopeOscillator* theOsc = oscManager->GetOscillatorIonisation(material,index);
124  if (theOsc->GetParentZ() != Z || theOsc->GetShellFlag()-1 != G4int(shell))
125  {
126  //something went wrong!
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;
130  ed << "Retrieved: " << theOsc->GetShellFlag()-1 << " and Z = " << theOsc->GetParentZ() << G4endl;
131  G4Exception("G4PenelopeIonisationCrossSection::CrossSection()","em2043",
132  JustWarning,ed);
133  return cross;
134  }
135 
136 
137  G4double crossPerMolecule = (theXS) ? theXS->GetShellCrossSection(index,incidentEnergy) : 0.;
138 
139  //Now it must be converted in cross section per atom. I need the number of
140  //atoms of the given Z per molecule.
141  G4double atomsPerMolec = oscManager->GetNumberOfZAtomsPerMolecule(material,Z);
142  if (atomsPerMolec)
143  cross = crossPerMolecule/atomsPerMolec;
144 
145 
146  if (verboseLevel > 0)
147  {
148  G4cout << "Cross section of shell " << G4int(shell) << " and Z= " << Z;
149  G4cout << " of material: " << material->GetName() << " and energy = " << incidentEnergy/keV << " keV" << G4endl;
150  G4cout << "--> " << cross/barn << " barn" << G4endl;
151  G4cout << "Shell binding energy: " << theOsc->GetIonisationEnergy()/eV << " eV;" ;
152  G4cout << " resonance energy: " << theOsc->GetResonanceEnergy()/eV << "eV" << G4endl;
153  if (verboseLevel > 2)
154  {
155  G4cout << "Cross section per molecule: " << crossPerMolecule/barn << " barn" << G4endl;
156  G4cout << "Atoms " << Z << " per molecule: " << atomsPerMolec << G4endl;
157  }
158  }
159  }
160 
161  return cross;
162 }
std::ostringstream G4ExceptionDescription
Definition: globals.hh:76
const G4PenelopeCrossSection * GetCrossSectionTableForCouple(const G4ParticleDefinition *, const G4Material *, const G4double cut) const
G4PenelopeOscillator * GetOscillatorIonisation(const G4Material *, G4int)
int G4int
Definition: G4Types.hh:78
string material
Definition: eplot.py:19
G4GLOB_DLL std::ostream G4cout
static constexpr double eV
Definition: G4SIunits.hh:215
const G4int nmax
G4double GetResonanceEnergy() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
void BuildXSTable(const G4Material *, G4double cut, const G4ParticleDefinition *, G4bool isMaster=true)
This can be inkoved only by the master.
double G4double
Definition: G4Types.hh:76
G4double GetNumberOfZAtomsPerMolecule(const G4Material *, G4int Z)
static constexpr double barn
Definition: G4SIunits.hh:105
static constexpr double keV
Definition: G4SIunits.hh:216
G4double GetShellCrossSection(size_t shellID, G4double energy) const
Returns the hard cross section for the given shell (per molecule)

Here is the call graph for this function:

Here is the caller graph for this function:

std::vector< G4double > G4PenelopeIonisationCrossSection::GetCrossSection ( G4int  Z,
G4double  incidentEnergy,
G4double  mass,
G4double  deltaEnergy,
const G4Material mat 
)
virtual

Purely virtual method from the base interface. Returns the cross section for all levels of element Z in material mat at the given energy

Implements G4VhShellCrossSection.

Definition at line 166 of file G4PenelopeIonisationCrossSection.cc.

170 {
171  G4int nmax = std::min(nMaxLevels,transitionManager->NumberOfShells(Z));
172  std::vector<G4double> vec(nmax,0.0);
173  for(G4int i=0; i<nmax; ++i) {
174  vec[i] = CrossSection(Z, G4AtomicShellEnumerator(i), kinEnergy,0.,mat);
175  }
176  return vec;
177 }
int G4int
Definition: G4Types.hh:78
const G4int nmax
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double CrossSection(G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat)
G4AtomicShellEnumerator

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4PenelopeIonisationCrossSection::GetVerbosityLevel ( )
inline

Definition at line 95 of file G4PenelopeIonisationCrossSection.hh.

95 {return verboseLevel;};
std::vector< G4double > G4PenelopeIonisationCrossSection::Probabilities ( G4int  Z,
G4double  incidentEnergy,
G4double  mass,
G4double  deltaEnergy,
const G4Material mat 
)
virtual

Purely virtual method from the base interface. Returns the shell ionisation probabilities for the given Z in the material mat at the specified energy.

Implements G4VhShellCrossSection.

Definition at line 182 of file G4PenelopeIonisationCrossSection.cc.

187 {
188  std::vector<G4double> vec = GetCrossSection(Z, kinEnergy,0,0,mat);
189  size_t n = vec.size();
190  size_t i=0;
191  G4double sum = 0.0;
192  for(i=0; i<n; ++i) { sum += vec[i]; }
193  if(sum > 0.0) {
194  sum = 1.0/sum;
195  for(i=0; i<n; ++i) { vec[i] = vec[i]*sum; }
196  }
197  return vec;
198 }
const G4int n
std::vector< G4double > GetCrossSection(G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

void G4PenelopeIonisationCrossSection::SetVerbosityLevel ( G4int  vl)
inline

Getter/setter for the verbosity level.

Definition at line 94 of file G4PenelopeIonisationCrossSection.hh.

94 {verboseLevel = vl;};

The documentation for this class was generated from the following files: