| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
>> 27 // $Id: G4hRDEnergyLoss.hh,v 1.2 2010/11/19 17:16:09 pia Exp $
>> 28 // GEANT4 tag $Name: geant4-09-04 $
27 // 29 //
>> 30 // $Id:
28 // ------------------------------------------- 31 // ------------------------------------------------------------
29 // GEANT 4 class header file 32 // GEANT 4 class header file
30 // 33 //
31 // History: first implementation, based o 34 // History: first implementation, based on object model of
32 // 2nd December 1995, G.Cosmo 35 // 2nd December 1995, G.Cosmo
33 // ---------- G4hEnergyLoss physics proce 36 // ---------- G4hEnergyLoss physics process -----------
34 // by Laszlo Urban, 30 May 1997 37 // by Laszlo Urban, 30 May 1997
35 // 38 //
36 // ******************************************* 39 // ************************************************************
37 // It is the first implementation of the NEW U 40 // It is the first implementation of the NEW UNIFIED ENERGY LOSS PROCESS.
38 // It calculates the continuous energy loss fo 41 // It calculates the continuous energy loss for charged hadrons.
39 // Processes giving contribution to the contin 42 // Processes giving contribution to the continuous loss :
40 // ionisation (= cont.ion.loss + delta ray p 43 // ionisation (= cont.ion.loss + delta ray production)
41 // can be added more easily .......... 44 // can be added more easily ..........
42 // This class creates static proton/antiproton 45 // This class creates static proton/antiproton dE/dx and range tables ,
43 // which tables can be used by other processes 46 // which tables can be used by other processes.
44 // The energy loss for other charged hadrons i 47 // The energy loss for other charged hadrons is calculated from the p/pbar
45 // tables with scaled kinetic energy. 48 // tables with scaled kinetic energy.
46 // 49 //
47 // 7/10/98 L.Urban some bugs fixed + some cle 50 // 7/10/98 L.Urban some bugs fixed + some cleanup
48 // 22/10/98 L.Urban cleanup 51 // 22/10/98 L.Urban cleanup
49 // 02/02/99 L.Urban several bugs fixed 52 // 02/02/99 L.Urban several bugs fixed
50 // 31/03/00 V.Ivanchenko rename to lowenergy a 53 // 31/03/00 V.Ivanchenko rename to lowenergy as G4hLowEnergyLoss.hh
51 // 09/08/00 V.Ivanchenko remove GetContinuousS 54 // 09/08/00 V.Ivanchenko remove GetContinuousStepLimit and IsApplicable
52 // 23/11/01 V.Ivanchenko Move static member-fu 55 // 23/11/01 V.Ivanchenko Move static member-functions from header to source
53 // 22/01/03 V.Ivanchenko Cuts per region 56 // 22/01/03 V.Ivanchenko Cuts per region
54 // 18/04/03 V.Ivanchenko Make dRoverRange prot 57 // 18/04/03 V.Ivanchenko Make dRoverRange protected
55 // 58 //
56 // 31 Jul 2008 MGP Short term supply of en 59 // 31 Jul 2008 MGP Short term supply of energy loss of hadrons through clone of
57 // former G4hLowEnergyLoss 60 // former G4hLowEnergyLoss (with some initial cleaning)
58 // To be replaced by rewor 61 // To be replaced by reworked class to deal with condensed/discrete
59 // issues properly 62 // issues properly
60 // 63 //
61 // ------------------------------------------- 64 // --------------------------------------------------------------
62 65
63 // Class description: 66 // Class description:
64 // Short term supply of energy loss of hadrons 67 // Short term supply of energy loss of hadrons through clone of former G4hLowEnergyLoss
65 // (with some initial cleaning) 68 // (with some initial cleaning)
66 // To be replaced by reworked class to deal wi 69 // To be replaced by reworked class to deal with condensed/discrete issues properly
67 70
68 // ------------------------------------------- 71 // --------------------------------------------------------------
69 72
70 73
71 #ifndef G4HRDENERGYLOSS_HH 74 #ifndef G4HRDENERGYLOSS_HH
72 #define G4HRDENERGYLOSS_HH 1 75 #define G4HRDENERGYLOSS_HH 1
73 76
74 #include "G4ios.hh" 77 #include "G4ios.hh"
75 #include "globals.hh" 78 #include "globals.hh"
76 #include "Randomize.hh" 79 #include "Randomize.hh"
77 #include "G4VContinuousDiscreteProcess.hh" 80 #include "G4VContinuousDiscreteProcess.hh"
78 #include "G4Material.hh" 81 #include "G4Material.hh"
79 #include "G4Element.hh" 82 #include "G4Element.hh"
80 #include "G4Proton.hh" 83 #include "G4Proton.hh"
81 #include "G4AntiProton.hh" 84 #include "G4AntiProton.hh"
82 #include "G4Electron.hh" 85 #include "G4Electron.hh"
83 #include "G4VParticleChange.hh" 86 #include "G4VParticleChange.hh"
84 #include "G4Track.hh" 87 #include "G4Track.hh"
85 #include "G4Step.hh" 88 #include "G4Step.hh"
86 #include "G4PhysicsLogVector.hh" 89 #include "G4PhysicsLogVector.hh"
87 #include "G4PhysicsLinearVector.hh" 90 #include "G4PhysicsLinearVector.hh"
88 91
89 class G4EnergyLossMessenger; 92 class G4EnergyLossMessenger;
90 93
91 class G4hRDEnergyLoss : public G4VContinuousDi 94 class G4hRDEnergyLoss : public G4VContinuousDiscreteProcess
92 95
93 { 96 {
94 public: 97 public:
95 98
96 G4hRDEnergyLoss(const G4String& ); 99 G4hRDEnergyLoss(const G4String& );
97 100
98 ~G4hRDEnergyLoss(); 101 ~G4hRDEnergyLoss();
99 102
100 virtual G4double GetMeanFreePath( 103 virtual G4double GetMeanFreePath(
101 const G4Track& track, 104 const G4Track& track,
102 G4double previousStepSize, 105 G4double previousStepSize,
103 enum G4ForceCondition* condition 106 enum G4ForceCondition* condition
104 ) = 0 ; 107 ) = 0 ;
105 108
106 virtual G4VParticleChange* PostStepDoIt(cons 109 virtual G4VParticleChange* PostStepDoIt(const G4Track& track,
107 const G4Step& Step) = 0 ; 110 const G4Step& Step) = 0 ;
108 111
109 // ---- MGP ---- All this static stuff is ex 112 // ---- MGP ---- All this static stuff is expected to disappear in a future
110 // development cycle 113 // development cycle
111 114
112 // get the number of processes contributing 115 // get the number of processes contributing to the cont.energy loss
113 static G4int GetNumberOfProcesses(); 116 static G4int GetNumberOfProcesses();
114 117
115 // set the number of processes contributing 118 // set the number of processes contributing to the cont.energy loss
116 static void SetNumberOfProcesses(G4int numbe 119 static void SetNumberOfProcesses(G4int number);
117 120
118 // Increment the number of processes contri 121 // Increment the number of processes contributing to the cont.energy loss
119 static void PlusNumberOfProcesses(); 122 static void PlusNumberOfProcesses();
120 123
121 // decrement the number of processes contri 124 // decrement the number of processes contributing to the cont.energy loss
122 static void MinusNumberOfProcesses(); 125 static void MinusNumberOfProcesses();
123 126
124 static void SetdRoverRange(G4double value); 127 static void SetdRoverRange(G4double value);
125 static void SetRndmStep (G4bool value) 128 static void SetRndmStep (G4bool value);
126 static void SetEnlossFluc (G4bool value) 129 static void SetEnlossFluc (G4bool value);
127 static void SetStepFunction (G4double c1, G4 130 static void SetStepFunction (G4double c1, G4double c2);
128 131
129 protected: 132 protected:
130 133
131 G4bool CutsWhereModified(); 134 G4bool CutsWhereModified();
132 135
133 // G4Material *lastMaterial ; 136 // G4Material *lastMaterial ;
134 const G4double MaxExcitationNumber ; 137 const G4double MaxExcitationNumber ;
135 const G4double probLimFluct ; 138 const G4double probLimFluct ;
136 const long nmaxDirectFluct,nmaxCont1,nmaxCon 139 const long nmaxDirectFluct,nmaxCont1,nmaxCont2 ;
137 140
138 static void BuildDEDXTable(const G4ParticleD 141 static void BuildDEDXTable(const G4ParticleDefinition& aParticleType);
139 142
140 protected: 143 protected:
141 144
142 static G4ThreadLocal G4PhysicsTable* theDEDX << 145 static G4PhysicsTable* theDEDXpTable ;
143 static G4ThreadLocal G4PhysicsTable* theDEDX << 146 static G4PhysicsTable* theDEDXpbarTable ;
144 static G4ThreadLocal G4PhysicsTable* theRang << 147 static G4PhysicsTable* theRangepTable ;
145 static G4ThreadLocal G4PhysicsTable* theRang << 148 static G4PhysicsTable* theRangepbarTable ;
146 149
147 //inverse of the range tables 150 //inverse of the range tables
148 static G4ThreadLocal G4PhysicsTable* theInve << 151 static G4PhysicsTable* theInverseRangepTable ;
149 static G4ThreadLocal G4PhysicsTable* theInve << 152 static G4PhysicsTable* theInverseRangepbarTable ;
150 153
151 //lab and proper time tables 154 //lab and proper time tables
152 static G4ThreadLocal G4PhysicsTable* theLabT << 155 static G4PhysicsTable* theLabTimepTable ;
153 static G4ThreadLocal G4PhysicsTable* theLabT << 156 static G4PhysicsTable* theLabTimepbarTable ;
154 157
155 static G4ThreadLocal G4PhysicsTable* theProp << 158 static G4PhysicsTable* theProperTimepTable ;
156 static G4ThreadLocal G4PhysicsTable* theProp << 159 static G4PhysicsTable* theProperTimepbarTable ;
157 160
158 // processes inherited from G4hRDEnergyLoss 161 // processes inherited from G4hRDEnergyLoss
159 // register themselves in the static arra 162 // register themselves in the static array Recorder
160 static G4ThreadLocal G4PhysicsTable** Record << 163 static G4PhysicsTable** RecorderOfpProcess;
161 static G4ThreadLocal G4PhysicsTable** Record << 164 static G4PhysicsTable** RecorderOfpbarProcess;
162 static G4ThreadLocal G4int CounterOfpProcess << 165 static G4int CounterOfpProcess ;
163 static G4ThreadLocal G4int CounterOfpbarProc << 166 static G4int CounterOfpbarProcess ;
164 167
165 // particle mass 168 // particle mass
166 static G4ThreadLocal G4double ParticleMass ; << 169 static G4double ParticleMass ;
167 170
168 // cut in range 171 // cut in range
169 static G4ThreadLocal G4double ptableElectron << 172 static G4double ptableElectronCutInRange;
170 static G4ThreadLocal G4double pbartableElect << 173 static G4double pbartableElectronCutInRange;
171 174
172 static G4ThreadLocal G4double Charge ; << 175 static G4double Charge ;
173 176
174 static G4ThreadLocal G4double LowestKineticE << 177 static G4double LowestKineticEnergy;
175 static G4ThreadLocal G4double HighestKinetic << 178 static G4double HighestKineticEnergy;
176 static G4ThreadLocal G4int TotBin; // number << 179 static G4int TotBin; // number of bins in table,
177 // calculated in BuildPhysicsTable 180 // calculated in BuildPhysicsTable
178 181
179 static G4ThreadLocal G4double RTable,LOGRTab << 182 static G4double RTable,LOGRTable; // LOGRTable=std::log(HighestKineticEnergy
180 // /LowestKineticEnergy)/TotBin 183 // /LowestKineticEnergy)/TotBin
181 // RTable = std::exp(LOGRTable) 184 // RTable = std::exp(LOGRTable)
182 185
183 G4PhysicsTable* theLossTable ; 186 G4PhysicsTable* theLossTable ;
184 187
185 G4double linLossLimit ; 188 G4double linLossLimit ;
186 189
187 G4double MinKineticEnergy ; 190 G4double MinKineticEnergy ;
188 191
189 static G4ThreadLocal G4double dRoverRange ; << 192 static G4double dRoverRange ; // maximum allowed deltarange/range
190 // in one step 193 // in one step
191 static G4ThreadLocal G4double finalRange ; << 194 static G4double finalRange ; // last step before stop
192 static G4ThreadLocal G4double c1lim,c2lim,c3 << 195 static G4double c1lim,c2lim,c3lim ; // coeffs for computing steplimit
193 196
194 static G4ThreadLocal G4bool rndmStepFlag ; << 197 static G4bool rndmStepFlag ;
195 static G4ThreadLocal G4bool EnlossFlucFlag ; << 198 static G4bool EnlossFlucFlag ;
196 199
197 200
198 private: 201 private:
199 202
200 // hide assignment operator 203 // hide assignment operator
201 204
202 G4hRDEnergyLoss(G4hRDEnergyLoss &); 205 G4hRDEnergyLoss(G4hRDEnergyLoss &);
203 G4hRDEnergyLoss & operator=(const G4hRDEnerg 206 G4hRDEnergyLoss & operator=(const G4hRDEnergyLoss &right);
204 207
205 // variables for the integration routines 208 // variables for the integration routines
206 static G4ThreadLocal G4double Mass,taulow,ta << 209 static G4double Mass,taulow,tauhigh,ltaulow,ltauhigh;
207 210
208 // ========================================= 211 // ====================================================================
209 // static part of the class 212 // static part of the class
210 213
211 static void BuildRangeTable(const G4Particle 214 static void BuildRangeTable(const G4ParticleDefinition& aParticleType);
212 215
213 static void BuildInverseRangeTable(const G4P 216 static void BuildInverseRangeTable(const G4ParticleDefinition& aParticleType);
214 217
215 static void BuildTimeTables(const G4Particle 218 static void BuildTimeTables(const G4ParticleDefinition& aParticleType);
216 219
217 static void BuildLabTimeVector(G4int materia 220 static void BuildLabTimeVector(G4int materialIndex,
218 G4PhysicsLogVector* rangeVector); 221 G4PhysicsLogVector* rangeVector);
219 222
220 static void BuildProperTimeVector(G4int mate 223 static void BuildProperTimeVector(G4int materialIndex,
221 G4PhysicsLogVector* rangeVector); 224 G4PhysicsLogVector* rangeVector);
222 225
223 static void InvertRangeVector(G4int material 226 static void InvertRangeVector(G4int materialIndex,
224 G4PhysicsLogVector* rangeVector); 227 G4PhysicsLogVector* rangeVector);
225 228
226 static void BuildRangeVector(G4int materialI 229 static void BuildRangeVector(G4int materialIndex,
227 G4PhysicsLogVector* rangeVector); 230 G4PhysicsLogVector* rangeVector);
228 231
229 static G4double LabTimeIntLog(G4PhysicsVecto 232 static G4double LabTimeIntLog(G4PhysicsVector* physicsVector, G4int nbin);
230 233
231 static G4double ProperTimeIntLog(G4PhysicsVe 234 static G4double ProperTimeIntLog(G4PhysicsVector* physicsVector, G4int nbin);
232 235
233 static G4double RangeIntLin(G4PhysicsVector* 236 static G4double RangeIntLin(G4PhysicsVector* physicsVector, G4int nbin);
234 237
235 static G4double RangeIntLog(G4PhysicsVector* 238 static G4double RangeIntLog(G4PhysicsVector* physicsVector, G4int nbin);
236 239
237 static void BuildRangeCoeffATable( const G4P 240 static void BuildRangeCoeffATable( const G4ParticleDefinition& aParticleType);
238 static void BuildRangeCoeffBTable( const G4P 241 static void BuildRangeCoeffBTable( const G4ParticleDefinition& aParticleType);
239 static void BuildRangeCoeffCTable(const G4Pa 242 static void BuildRangeCoeffCTable(const G4ParticleDefinition& aParticleType);
240 243
241 // ========================================= 244 // ====================================================================
242 245
243 static G4ThreadLocal G4PhysicsTable* theDEDX << 246 static G4PhysicsTable* theDEDXTable;
244 247
245 static G4ThreadLocal G4PhysicsTable* theRang << 248 static G4PhysicsTable* theRangeTable;
246 static G4ThreadLocal G4PhysicsTable* theInve << 249 static G4PhysicsTable* theInverseRangeTable;
247 250
248 static G4ThreadLocal G4PhysicsTable* theLabT << 251 static G4PhysicsTable* theLabTimeTable;
249 static G4ThreadLocal G4PhysicsTable* theProp << 252 static G4PhysicsTable* theProperTimeTable;
250 253
251 static G4ThreadLocal G4PhysicsTable** Record << 254 static G4PhysicsTable** RecorderOfProcess;
252 static G4ThreadLocal G4int CounterOfProcess; << 255 static G4int CounterOfProcess;
253 << 256
254 static G4ThreadLocal G4PhysicsTable* thepRan << 257 static G4PhysicsTable* thepRangeCoeffATable;
255 static G4ThreadLocal G4PhysicsTable* thepRan << 258 static G4PhysicsTable* thepRangeCoeffBTable;
256 static G4ThreadLocal G4PhysicsTable* thepRan << 259 static G4PhysicsTable* thepRangeCoeffCTable;
257 static G4ThreadLocal G4PhysicsTable* thepbar << 260 static G4PhysicsTable* thepbarRangeCoeffATable;
258 static G4ThreadLocal G4PhysicsTable* thepbar << 261 static G4PhysicsTable* thepbarRangeCoeffBTable;
259 static G4ThreadLocal G4PhysicsTable* thepbar << 262 static G4PhysicsTable* thepbarRangeCoeffCTable;
260 << 263
261 static G4ThreadLocal G4PhysicsTable* theRang << 264 static G4PhysicsTable* theRangeCoeffATable;
262 static G4ThreadLocal G4PhysicsTable* theRang << 265 static G4PhysicsTable* theRangeCoeffBTable;
263 static G4ThreadLocal G4PhysicsTable* theRang << 266 static G4PhysicsTable* theRangeCoeffCTable;
264 static G4ThreadLocal G4int NumberOfProcesses << 267 static G4int NumberOfProcesses ;
265 268
266 }; 269 };
267 270
268 #endif 271 #endif
>> 272
>> 273
>> 274
269 275