Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PAIySection.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // $Id: G4PAIySection.hh 96934 2016-05-18 09:10:41Z gcosmo $
28 //
29 //
30 // G4PAIySection.hh -- header file
31 //
32 //
33 // Preparation of ionizing collision cross section according to Photo Absorption
34 // Ionization (PAI) model for simulation of ionization energy losses in very thin
35 // absorbers. Author: Vladimir.Grichine@cern.ch
36 //
37 // History:
38 //
39 // 01.10.07, V.Ivanchenko create using V.Grichine G4PAIxSection class
40 // 21.11.10, V.Grichine fVerbose and SetVerbose added
41 // 28.10.11, V.Ivanchenko Migration of exceptions to the new design
42 
43 #ifndef G4PAIYSECTION_HH
44 #define G4PAIYSECTION_HH
45 
46 #include "G4ios.hh"
47 #include "globals.hh"
48 #include "Randomize.hh"
49 
50 #include "G4SandiaTable.hh"
51 
53 {
54 public:
55 
56  explicit G4PAIySection();
57 
59 
60  void Initialize(const G4Material* material, G4double maxEnergyTransfer,
61  G4double betaGammaSq, G4SandiaTable*);
62 
63  void ComputeLowEnergyCof(const G4Material* material);
64 
65  void InitPAI();
66 
67  void NormShift( G4double betaGammaSq );
68 
69  void SplainPAI( G4double betaGammaSq );
70 
71  // Physical methods
72  G4double RutherfordIntegral( G4int intervalNumber,
73  G4double limitLow,
74  G4double limitHigh );
75 
76  G4double ImPartDielectricConst( G4int intervalNumber,
77  G4double energy );
78 
80 
81  G4double DifPAIySection( G4int intervalNumber,
82  G4double betaGammaSq );
83 
84  G4double PAIdNdxCerenkov( G4int intervalNumber,
85  G4double betaGammaSq );
86 
87  G4double PAIdNdxPlasmon( G4int intervalNumber,
88  G4double betaGammaSq );
89 
90  void IntegralPAIySection();
91  void IntegralCerenkov();
92  void IntegralPlasmon();
93 
94  G4double SumOverInterval(G4int intervalNumber);
95  G4double SumOverIntervaldEdx(G4int intervalNumber);
96  G4double SumOverInterCerenkov(G4int intervalNumber);
97  G4double SumOverInterPlasmon(G4int intervalNumber);
98 
99  G4double SumOverBorder( G4int intervalNumber,
100  G4double energy );
101  G4double SumOverBorderdEdx( G4int intervalNumber,
102  G4double energy );
103  G4double SumOverBordCerenkov( G4int intervalNumber,
104  G4double energy );
105  G4double SumOverBordPlasmon( G4int intervalNumber,
106  G4double energy );
107 
111 
112  // Inline access functions
113 
114  inline G4int GetNumberOfGammas() const { return fNumberOfGammas; }
115 
116  inline G4int GetSplineSize() const { return fSplineNumber; }
117 
118  inline G4int GetIntervalNumber() const { return fIntervalNumber; }
119 
120  inline G4double GetEnergyInterval(G4int i){ return fEnergyInterval[i]; }
121 
122  inline G4double GetDifPAIySection(G4int i){ return fDifPAIySection[i]; }
123  inline G4double GetPAIdNdxCrenkov(G4int i){ return fdNdxCerenkov[i]; }
124  inline G4double GetPAIdNdxPlasmon(G4int i){ return fdNdxPlasmon[i]; }
125 
126  inline G4double GetMeanEnergyLoss() const {return fIntegralPAIySection[0]; }
127  inline G4double GetMeanCerenkovLoss() const {return fIntegralCerenkov[0]; }
128  inline G4double GetMeanPlasmonLoss() const {return fIntegralPlasmon[0]; }
129 
130  inline G4double GetNormalizationCof() const { return fNormalizationCof; }
131 
132  inline G4double GetPAItable(G4int i,G4int j) const;
133 
134  inline G4double GetLorentzFactor(G4int i) const;
135 
136  inline G4double GetSplineEnergy(G4int i) const;
137 
138  inline G4double GetIntegralPAIySection(G4int i) const;
139  inline G4double GetIntegralPAIdEdx(G4int i) const;
140  inline G4double GetIntegralCerenkov(G4int i) const;
141  inline G4double GetIntegralPlasmon(G4int i) const;
142 
143  inline void SetVerbose(G4int v) { fVerbose = v; };
144 
145 private :
146 
147  void CallError(G4int i, const G4String& methodName) const;
148 
149  // Local class constants
150 
151  static const G4double fDelta; // energy shift from interval border = 0.001
152  static const G4double fError; // error in lin-log approximation = 0.005
153 
154  static G4int fNumberOfGammas; // = 111;
155  static const G4double fLorentzFactor[112]; // static gamma array
156 
157  static
158  const G4int fRefGammaNumber; // The number of gamma for creation of spline (15)
159 
160  G4int fIntervalNumber ; // The number of energy intervals
161  G4double fNormalizationCof; // Normalization cof for PhotoAbsorptionXsection
162 
163  G4double betaBohr;
164  G4double betaBohr4;
165 
166  G4double fDensity; // Current density
167  G4double fElectronDensity; // Current electron (number) density
168  G4double fLowEnergyCof; // Correction cof for low energy region
169  G4int fSplineNumber; // Current size of spline
170  G4int fVerbose; // verbose flag
171 
172  G4SandiaTable* fSandia;
173 
174  G4DataVector fEnergyInterval;
175  G4DataVector fA1;
176  G4DataVector fA2;
177  G4DataVector fA3;
178  G4DataVector fA4;
179 
180  static
181  const G4int fMaxSplineSize; // Max size of output splain arrays = 500
182 
183  G4DataVector fSplineEnergy; // energy points of splain
184  G4DataVector fRePartDielectricConst; // Real part of dielectric const
185  G4DataVector fImPartDielectricConst; // Imaginary part of dielectric const
186  G4DataVector fIntegralTerm; // Integral term in PAI cross section
187  G4DataVector fDifPAIySection; // Differential PAI cross section
188  G4DataVector fdNdxCerenkov; // dNdx of Cerenkov collisions
189  G4DataVector fdNdxPlasmon; // dNdx of Plasmon collisions
190 
191  G4DataVector fIntegralPAIySection; // Integral PAI cross section ?
192  G4DataVector fIntegralPAIdEdx; // Integral PAI dEdx ?
193  G4DataVector fIntegralCerenkov; // Integral Cerenkov N>omega ?
194  G4DataVector fIntegralPlasmon; // Integral Plasmon N>omega ?
195 
196  G4double fPAItable[500][112]; // Output array
197 };
198 
200 {
201  return fPAItable[i][j];
202 }
203 
205 {
206  return fLorentzFactor[j];
207 }
208 
210 {
211  if(i < 1 || i > fSplineNumber) { CallError(i, "GetSplineEnergy"); }
212  return fSplineEnergy[i];
213 }
214 
216 {
217  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralPAIySection"); }
218  return fIntegralPAIySection[i];
219 }
220 
222 {
223  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralPAIdEdx"); }
224  return fIntegralPAIdEdx[i];
225 }
226 
228 {
229  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralCerenkov"); }
230  return fIntegralCerenkov[i];
231 }
232 
234 {
235  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralPlasmon"); }
236  return fIntegralPlasmon[i];
237 }
238 
239 #endif
240 
241 // ----------------- end of G4PAIySection header file -------------------
G4double SumOverIntervaldEdx(G4int intervalNumber)
G4double RutherfordIntegral(G4int intervalNumber, G4double limitLow, G4double limitHigh)
G4double GetDifPAIySection(G4int i)
G4double GetPAIdNdxCrenkov(G4int i)
G4double GetNormalizationCof() const
void NormShift(G4double betaGammaSq)
G4double GetMeanEnergyLoss() const
G4double ImPartDielectricConst(G4int intervalNumber, G4double energy)
G4double SumOverInterCerenkov(G4int intervalNumber)
G4double GetMeanPlasmonLoss() const
G4double PAIdNdxCerenkov(G4int intervalNumber, G4double betaGammaSq)
void ComputeLowEnergyCof(const G4Material *material)
void IntegralCerenkov()
G4int GetNumberOfGammas() const
G4double SumOverBorderdEdx(G4int intervalNumber, G4double energy)
int G4int
Definition: G4Types.hh:78
void IntegralPlasmon()
G4double GetIntegralPAIdEdx(G4int i) const
G4double GetIntegralPlasmon(G4int i) const
string material
Definition: eplot.py:19
G4double SumOverBordPlasmon(G4int intervalNumber, G4double energy)
G4double GetStepEnergyLoss(G4double step)
void Initialize(const G4Material *material, G4double maxEnergyTransfer, G4double betaGammaSq, G4SandiaTable *)
G4double GetSplineEnergy(G4int i) const
G4double GetIntegralPAIySection(G4int i) const
tuple v
Definition: test.py:18
G4double GetIntegralCerenkov(G4int i) const
G4double GetStepCerenkovLoss(G4double step)
G4int GetIntervalNumber() const
void SetVerbose(G4int v)
G4double SumOverBorder(G4int intervalNumber, G4double energy)
G4double energy(const ThreeVector &p, const G4double m)
G4double GetPAItable(G4int i, G4int j) const
void SplainPAI(G4double betaGammaSq)
G4double SumOverBordCerenkov(G4int intervalNumber, G4double energy)
G4int GetSplineSize() const
double G4double
Definition: G4Types.hh:76
G4double RePartDielectricConst(G4double energy)
G4double GetStepPlasmonLoss(G4double step)
G4double GetEnergyInterval(G4int i)
G4double DifPAIySection(G4int intervalNumber, G4double betaGammaSq)
G4double SumOverInterval(G4int intervalNumber)
G4double PAIdNdxPlasmon(G4int intervalNumber, G4double betaGammaSq)
G4double GetMeanCerenkovLoss() const
void IntegralPAIySection()
G4double GetPAIdNdxPlasmon(G4int i)
G4double SumOverInterPlasmon(G4int intervalNumber)
G4double GetLorentzFactor(G4int i) const