source: trunk/source/processes/electromagnetic/muons/include/G4EnergyLossForExtrapolator.hh@ 1201

Last change on this file since 1201 was 1196, checked in by garnier, 16 years ago

update CVS release candidate geant4.9.3.01
File size: 10.6 KB
Line 
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// $Id: G4EnergyLossForExtrapolator.hh,v 1.12 2008/11/13 14:14:07 vnivanch Exp $
27// GEANT4 tag $Name: geant4-09-03-cand-01 $
28//
29//---------------------------------------------------------------------------
30//
31// ClassName: G4EnergyLossForExtrapolator
32//
33// Description: This class provide calculation of energy loss, fluctuation,
34// and msc angle
35//
36// Author: 09.12.04 V.Ivanchenko
37//
38// Modification:
39// 08-04-05 Rename Propogator -> Extrapolator
40// 16-03-06 Add muon tables
41// 21-03-06 Add verbosity defined in the constructor and Initialisation
42// start only when first public method is called (V.Ivanchenko)
43// 03-05-06 Remove unused pointer G4Material* from number of methods (VI)
44// 28-07-07 Add maxEnergyTransfer for computation of energy loss (VI)
45//
46//----------------------------------------------------------------------------
47//
48
49//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
50
51#ifndef G4EnergyLossForExtrapolator_h
52#define G4EnergyLossForExtrapolator_h 1
53
54
55#include "globals.hh"
56#include "G4PhysicsTable.hh"
57#include <vector>
58
59class G4ParticleDefinition;
60class G4Material;
61class G4MaterialCutsCouple;
62class G4ProductionCuts;
63
64//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
65
66class G4EnergyLossForExtrapolator
67{
68public:
69 G4EnergyLossForExtrapolator(G4int verb = 1);
70
71 ~G4EnergyLossForExtrapolator();
72
73 G4double ComputeDEDX(G4double kinEnergy, const G4ParticleDefinition*);
74
75 G4double ComputeRange(G4double kinEnergy, const G4ParticleDefinition*);
76
77 G4double ComputeEnergy(G4double range, const G4ParticleDefinition*);
78
79 G4double EnergyAfterStep(G4double kinEnergy, G4double step,
80 const G4Material*, const G4ParticleDefinition*);
81
82 G4double EnergyBeforeStep(G4double kinEnergy, G4double step,
83 const G4Material*, const G4ParticleDefinition*);
84
85 G4double TrueStepLength(G4double kinEnergy, G4double step,
86 const G4Material*, const G4ParticleDefinition* part);
87
88 inline G4double EnergyAfterStep(G4double kinEnergy, G4double step,
89 const G4Material*, const G4String& particleName);
90
91 inline G4double EnergyBeforeStep(G4double kinEnergy, G4double step,
92 const G4Material*, const G4String& particleName);
93
94 inline G4double AverageScatteringAngle(G4double kinEnergy, G4double step,
95 const G4Material*,
96 const G4ParticleDefinition* part);
97
98 inline G4double AverageScatteringAngle(G4double kinEnergy, G4double step,
99 const G4Material*,
100 const G4String& particleName);
101
102 inline G4double ComputeTrueStep(const G4Material*, const G4ParticleDefinition* part,
103 G4double kinEnergy, G4double stepLength);
104
105 inline G4double EnergyDispersion(G4double kinEnergy, G4double step,
106 const G4Material*, const G4ParticleDefinition*);
107
108 inline G4double EnergyDispersion(G4double kinEnergy, G4double step,
109 const G4Material*, const G4String& particleName);
110
111 inline void SetVerbose(G4int val);
112
113 inline void SetMinKinEnergy(G4double);
114
115 inline void SetMaxKinEnergy(G4double);
116
117 inline void SetMaxEnergyTransfer(G4double);
118
119private:
120
121 void Initialisation();
122
123 G4bool SetupKinematics(const G4ParticleDefinition*, const G4Material*,
124 G4double kinEnergy);
125
126 G4PhysicsTable* PrepareTable();
127
128 const G4ParticleDefinition* FindParticle(const G4String& name);
129
130 void ComputeElectronDEDX(const G4ParticleDefinition* part, G4PhysicsTable* table);
131
132 void ComputeMuonDEDX(const G4ParticleDefinition* part, G4PhysicsTable* table);
133
134 void ComputeProtonDEDX(const G4ParticleDefinition* part, G4PhysicsTable* table);
135
136 void ComputeTrasportXS(const G4ParticleDefinition* part, G4PhysicsTable* table);
137
138 inline G4double ComputeValue(G4double x, const G4PhysicsTable* table);
139
140 // hide assignment operator
141 G4EnergyLossForExtrapolator & operator=(const G4EnergyLossForExtrapolator &right);
142 G4EnergyLossForExtrapolator(const G4EnergyLossForExtrapolator&);
143
144 const G4ParticleDefinition* currentParticle;
145 const G4ParticleDefinition* electron;
146 const G4ParticleDefinition* positron;
147 const G4ParticleDefinition* muonPlus;
148 const G4ParticleDefinition* muonMinus;
149 const G4ParticleDefinition* proton;
150
151 G4ProductionCuts* cuts;
152 std::vector<const G4MaterialCutsCouple*> couples;
153
154 G4String currentParticleName;
155
156 G4PhysicsTable* dedxElectron;
157 G4PhysicsTable* dedxPositron;
158 G4PhysicsTable* dedxMuon;
159 G4PhysicsTable* dedxProton;
160 G4PhysicsTable* rangeElectron;
161 G4PhysicsTable* rangePositron;
162 G4PhysicsTable* rangeMuon;
163 G4PhysicsTable* rangeProton;
164 G4PhysicsTable* invRangeElectron;
165 G4PhysicsTable* invRangePositron;
166 G4PhysicsTable* invRangeMuon;
167 G4PhysicsTable* invRangeProton;
168 G4PhysicsTable* mscElectron;
169
170 const G4Material* currentMaterial;
171 G4int index;
172 G4double electronDensity;
173 G4double radLength;
174 G4double mass;
175 G4double charge2;
176 G4double kineticEnergy;
177 G4double gam;
178 G4double bg2;
179 G4double beta2;
180 G4double tmax;
181
182 G4double linLossLimit;
183 G4double emin;
184 G4double emax;
185 G4double maxEnergyTransfer;
186
187 G4int nbins;
188 G4int nmat;
189 G4int verbose;
190 G4bool isInitialised;
191};
192
193//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
194
195inline G4double G4EnergyLossForExtrapolator::EnergyAfterStep(G4double kinEnergy,
196 G4double step,
197 const G4Material* mat,
198 const G4String& name)
199{
200 return EnergyAfterStep(kinEnergy,step,mat,FindParticle(name));
201}
202
203//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
204
205inline G4double G4EnergyLossForExtrapolator::EnergyBeforeStep(G4double kinEnergy,
206 G4double step,
207 const G4Material* mat,
208 const G4String& name)
209{
210 return EnergyBeforeStep(kinEnergy,step,mat,FindParticle(name));
211}
212
213//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
214
215inline G4double G4EnergyLossForExtrapolator::AverageScatteringAngle(G4double kinEnergy,
216 G4double step,
217 const G4Material* mat,
218 const G4String& name)
219{
220 return AverageScatteringAngle(kinEnergy,step,mat,FindParticle(name));
221}
222
223//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
224
225inline G4double
226G4EnergyLossForExtrapolator::EnergyDispersion(G4double kinEnergy,
227 G4double step,
228 const G4Material* mat,
229 const G4String& name)
230{
231 return EnergyDispersion(kinEnergy,step,mat,FindParticle(name));
232}
233
234//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
235
236inline G4double
237G4EnergyLossForExtrapolator::AverageScatteringAngle(G4double kinEnergy,
238 G4double stepLength,
239 const G4Material* mat,
240 const G4ParticleDefinition* part)
241{
242 G4double theta = 0.0;
243 if(SetupKinematics(part, mat, kinEnergy)) {
244 G4double t = stepLength/radLength;
245 G4double y = std::max(0.001, t);
246 theta = 19.23*MeV*std::sqrt(charge2*t)*(1.0 + 0.038*std::log(y))/(beta2*gam*mass);
247 }
248 return theta;
249}
250
251//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
252
253inline G4double
254G4EnergyLossForExtrapolator::ComputeTrueStep(const G4Material* mat,
255 const G4ParticleDefinition* part,
256 G4double kinEnergy,
257 G4double stepLength)
258{
259 G4double theta = AverageScatteringAngle(kinEnergy,stepLength,mat,part);
260 return stepLength*std::sqrt(1.0 + 0.625*theta*theta);
261}
262
263//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
264
265inline
266G4double G4EnergyLossForExtrapolator::EnergyDispersion(G4double kinEnergy,
267 G4double stepLength,
268 const G4Material* mat,
269 const G4ParticleDefinition* part)
270{
271 G4double sig2 = 0.0;
272 if(SetupKinematics(part, mat, kinEnergy)) {
273 G4double step = ComputeTrueStep(mat,part,kinEnergy,stepLength);
274 sig2 = (1.0/beta2 - 0.5)*twopi_mc2_rcl2*tmax*step*electronDensity*charge2;
275 }
276 return sig2;
277}
278
279//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
280
281inline G4double G4EnergyLossForExtrapolator::ComputeValue(G4double x,
282 const G4PhysicsTable* table)
283{
284 G4double res = 0.0;
285 G4bool b;
286 if(table) res = ((*table)[index])->GetValue(x, b);
287 return res;
288}
289
290//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
291
292inline void G4EnergyLossForExtrapolator::SetVerbose(G4int val)
293{
294 verbose = val;
295}
296
297//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
298
299inline void G4EnergyLossForExtrapolator::SetMinKinEnergy(G4double val)
300{
301 emin = val;
302}
303
304//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
305
306inline void G4EnergyLossForExtrapolator::SetMaxKinEnergy(G4double val)
307{
308 emax = val;
309}
310
311//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
312
313inline void G4EnergyLossForExtrapolator::SetMaxEnergyTransfer(G4double val)
314{
315 maxEnergyTransfer = val;
316}
317
318//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
319
320#endif
321
Note: See TracBrowser for help on using the repository browser.