GENIEGenerator
Loading...
Searching...
No Matches
BetheBlochMaterialParams.h
Go to the documentation of this file.
1//____________________________________________________________________________
2/*!
3
4\class genie::mueloss::BetheBlochMaterialParams
5
6\brief Bethe Bloch parameters for various materials
7
8\ref W.Lohmann, R.Kopp and R.Voss,
9 Energy Loss of Muons in the Energy Range 1-10000 GeV, CERN 85-03
10
11\author Costas Andreopoulos <c.andreopoulos \at cern.ch>
12 University of Liverpool
13
14\created December 10, 2003
15
16\cpright Copyright (c) 2003-2025, The GENIE Collaboration
17 For the full text of the license visit http://copyright.genie-mc.org
18*/
19//____________________________________________________________________________
20
21#ifndef _BETHE_BLOCH_MATERIAL_PARAMS_H_
22#define _BETHE_BLOCH_MATERIAL_PARAMS_H_
23
24namespace genie {
25namespace mueloss {
26
27#include "Physics/MuonEnergyLoss/MuELMaterial.h"
28
29class BetheBlochMaterialParams
30{
31public:
32 //____________________________________________________________________
33 static double IonizationPotential(MuELMaterial_t material)
34 {
35 // returns the ionization potential for the input material (in eV)
36
37 switch(material) {
38
39 //- compound materials
40 case eMuBariumFluoride: return 375.9; break;
41 case eMuBismuthGermanate: return 534.1; break;
42 case eMuPyrex: return 134.0; break;
43 case eMuCalciumCarbonate: return 136.4; break;
44 case eMuConcrete: return 135.2; break;
45 case eMuFreon12: return 143.0; break;
46 case eMuFreon13B1: return 210.5; break;
47 case eMuLeadOxide: return 766.7; break;
48 case eMuLithiumFluoride: return 94.0; break;
49 case eMuLucite: return 74.0; break;
50 case eMuPolyethylene: return 57.4; break;
51 case eMuPolystyrene: return 68.7; break;
52 case eMuLiquidPropane: return 52.0; break;
53 case eMuSiliconDioxide: return 139.2; break;
54 case eMuSodiumIodide: return 452.0; break;
55 case eMuStandardRock: return 136.4; break;
56 case eMuUraniumOxide: return 720.6; break;
57 case eMuWater: return 75.0; break;
58
59 //- elements
60 case eMuHydrogen: return 21.8; break;
61 case eMuDeuterium: return 21.8; break;
62 case eMuHelium: return 41.8; break;
63 case eMuLithium: return 40.0; break;
64 case eMuBeryllium: return 63.7; break;
65 case eMuBoron: return 76.0; break;
66 case eMuCarbon: return 78.0; break;
67 case eMuNitrogen: return 82.0; break;
68 case eMuOxygen: return 95.0; break;
69 case eMuFluorine: return 115.0; break;
70 case eMuNeon: return 137.0; break;
71 case eMuSodium: return 149.0; break;
72 case eMuMagnesium: return 156.0; break;
73 case eMuAluminium: return 166.0; break;
74 case eMuSilicon: return 173.0; break;
75 case eMuSulphur: return 180.0; break;
76 case eMuChlorine: return 174.0; break;
77 case eMuArgon: return 188.0; break;
78 case eMuPotassium: return 190.0; break;
79 case eMuCalcium: return 191.0; break;
80 case eMuChromium: return 257.0; break;
81 case eMuManganese: return 272.0; break;
82 case eMuIron: return 286.0; break;
83 case eMuNickel: return 311.0; break;
84 case eMuCopper: return 322.0; break;
85 case eMuZinc: return 330.0; break;
86 case eMuGermanium: return 350.0; break;
87 case eMuBromine: return 343.0; break;
88 case eMuTin: return 488.0; break;
89 case eMuIodine: return 491.0; break;
90 case eMuBarium: return 491.0; break;
91 case eMuTungsten: return 727.0; break;
92 case eMuLead: return 823.0; break;
93 case eMuBismuth: return 823.0; break;
94 case eMuUranium: return 890.0; break;
95
96 case eMuUndefined:
97 default:
98 return 0;
99 }
100 return 0;
101 }
102 //____________________________________________________________________
103 static double DensityCorrection_C(MuELMaterial_t material)
104 {
105 // returns the density correction factor C
106
107 switch(material) {
108
109 //- compound materials
110 case eMuBariumFluoride: return -5.412; break;
111 case eMuBismuthGermanate: return -5.741; break;
112 case eMuPyrex: return -3.971; break;
113 case eMuCalciumCarbonate: return -3.774; break;
114 case eMuConcrete: return -3.946; break;
115 case eMuFreon12: return -4.825; break;
116 case eMuFreon13B1: return -5.356; break;
117 case eMuLeadOxide: return -6.216; break;
118 case eMuLithiumFluoride: return -3.167; break;
119 case eMuLucite: return -3.330; break;
120 case eMuPolyethylene: return -3.002; break;
121 case eMuPolystyrene: return -3.300; break;
122 case eMuLiquidPropane: return -3.553; break;
123 case eMuSiliconDioxide: return -4.003; break;
124 case eMuSodiumIodide: return -6.057; break;
125 case eMuStandardRock: return -3.774; break;
126 case eMuUraniumOxide: return -5.961; break;
127 case eMuWater: return -3.502; break;
128
129 //- elements
130 case eMuHydrogen: return -3.263; break;
131 case eMuDeuterium: return -2.942; break;
132 case eMuHelium: return -4.517; break;
133 case eMuLithium: return -3.122; break;
134 case eMuBeryllium: return -2.785; break;
135 case eMuBoron: return -2.848; break;
136 case eMuCarbon: return -2.868; break;
137 case eMuNitrogen: return -3.998; break;
138 case eMuOxygen: return -3.948; break;
139 case eMuFluorine: return -4.413; break;
140 case eMuNeon: return -4.632; break;
141 case eMuSodium: return -5.053; break;
142 case eMuMagnesium: return -4.530; break;
143 case eMuAluminium: return -4.240; break;
144 case eMuSilicon: return -4.435; break;
145 case eMuSulphur: return -4.666; break;
146 case eMuChlorine: return -4.887; break;
147 case eMuArgon: return -5.217; break;
148 case eMuPotassium: return -5.642; break;
149 case eMuCalcium: return -5.040; break;
150 case eMuChromium: return -4.178; break;
151 case eMuManganese: return -4.270; break;
152 case eMuIron: return -4.291; break;
153 case eMuNickel: return -4.312; break;
154 case eMuCopper: return -4.419; break;
155 case eMuZinc: return -4.691; break;
156 case eMuGermanium: return -5.141; break;
157 case eMuBromine: return -5.641; break;
158 case eMuTin: return -5.534; break;
159 case eMuIodine: return -5.949; break;
160 case eMuBarium: return -6.315; break;
161 case eMuTungsten: return -5.406; break;
162 case eMuLead: return -6.202; break;
163 case eMuBismuth: return -6.351; break;
164 case eMuUranium: return -5.869; break;
165
166 case eMuUndefined:
167 default:
168 return 0;
169 }
170 return 0;
171 }
172 //____________________________________________________________________
173 static double DensityCorrection_X0(MuELMaterial_t material)
174 {
175 // returns the density correction constant X0
176
177 switch(material) {
178
179 //- compound materials
180 case eMuBariumFluoride: return -0.010; break;
181 case eMuBismuthGermanate: return 0.046; break;
182 case eMuPyrex: return 0.148; break;
183 case eMuCalciumCarbonate: return 0.049; break;
184 case eMuConcrete: return 0.130; break;
185 case eMuFreon12: return 0.304; break;
186 case eMuFreon13B1: return 0.352; break;
187 case eMuLeadOxide: return 0.036; break;
188 case eMuLithiumFluoride: return 0.017; break;
189 case eMuLucite: return 0.182; break;
190 case eMuPolyethylene: return 0.137; break;
191 case eMuPolystyrene: return 0.165; break;
192 case eMuLiquidPropane: return 0.286; break;
193 case eMuSiliconDioxide: return 0.139; break;
194 case eMuSodiumIodide: return 0.120; break;
195 case eMuStandardRock: return 0.049; break;
196 case eMuUraniumOxide: return -0.194; break;
197 case eMuWater: return 0.240; break;
198
199 //- elements
200 case eMuHydrogen: return 0.476; break;
201 case eMuDeuterium: return 0.200; break;
202 case eMuHelium: return 0.473; break;
203 case eMuLithium: return 0.130; break;
204 case eMuBeryllium: return 0.059; break;
205 case eMuBoron: return 0.031; break;
206 case eMuCarbon: return -0.018; break;
207 case eMuNitrogen: return 0.304; break;
208 case eMuOxygen: return 0.287; break;
209 case eMuFluorine: return 0.200; break;
210 case eMuNeon: return 0.200; break;
211 case eMuSodium: return 0.288; break;
212 case eMuMagnesium: return 0.150; break;
213 case eMuAluminium: return 0.171; break;
214 case eMuSilicon: return 0.201; break;
215 case eMuSulphur: return 0.158; break;
216 case eMuChlorine: return 0.200; break;
217 case eMuArgon: return 0.201; break;
218 case eMuPotassium: return 0.385; break;
219 case eMuCalcium: return 0.323; break;
220 case eMuChromium: return 0.034; break;
221 case eMuManganese: return 0.045; break;
222 case eMuIron: return -0.001; break;
223 case eMuNickel: return -0.057; break;
224 case eMuCopper: return -0.025; break;
225 case eMuZinc: return 0.005; break;
226 case eMuGermanium: return 0.338; break;
227 case eMuBromine: return 0.339; break;
228 case eMuTin: return 0.288; break;
229 case eMuIodine: return 0.055; break;
230 case eMuBarium: return 0.419; break;
231 case eMuTungsten: return 0.217; break;
232 case eMuLead: return 0.378; break;
233 case eMuBismuth: return 0.415; break;
234 case eMuUranium: return 0.226; break;
235
236 case eMuUndefined:
237 default:
238 return 0;
239 }
240 return 0;
241 }
242 //____________________________________________________________________
243 static double DensityCorrection_X1(MuELMaterial_t material)
244 {
245 // returns the density correction constant X1
246
247 switch(material) {
248
249 //- compound materials
250 case eMuBariumFluoride: return 3.387; break;
251 case eMuBismuthGermanate: return 3.782; break;
252 case eMuPyrex: return 2.993; break;
253 case eMuCalciumCarbonate: return 3.055; break;
254 case eMuConcrete: return 3.047; break;
255 case eMuFreon12: return 3.266; break;
256 case eMuFreon13B1: return 3.755; break;
257 case eMuLeadOxide: return 3.546; break;
258 case eMuLithiumFluoride: return 2.705; break;
259 case eMuLucite: return 2.668; break;
260 case eMuPolyethylene: return 2.518; break;
261 case eMuPolystyrene: return 2.503; break;
262 case eMuLiquidPropane: return 2.657; break;
263 case eMuSiliconDioxide: return 3.003; break;
264 case eMuSodiumIodide: return 3.592; break;
265 case eMuStandardRock: return 3.055; break;
266 case eMuUraniumOxide: return 3.529; break;
267 case eMuWater: return 2.800; break;
268
269 //- elements
270 case eMuHydrogen: return 1.922; break;
271 case eMuDeuterium: return 2.000; break;
272 case eMuHelium: return 2.000; break;
273 case eMuLithium: return 1.640; break;
274 case eMuBeryllium: return 1.692; break;
275 case eMuBoron: return 1.969; break;
276 case eMuCarbon: return 2.342; break;
277 case eMuNitrogen: return 2.000; break;
278 case eMuOxygen: return 2.000; break;
279 case eMuFluorine: return 3.000; break;
280 case eMuNeon: return 3.000; break;
281 case eMuSodium: return 3.196; break;
282 case eMuMagnesium: return 3.067; break;
283 case eMuAluminium: return 3.013; break;
284 case eMuSilicon: return 2.872; break;
285 case eMuSulphur: return 2.716; break;
286 case eMuChlorine: return 3.000; break;
287 case eMuArgon: return 3.000; break;
288 case eMuPotassium: return 3.172; break;
289 case eMuCalcium: return 3.119; break;
290 case eMuChromium: return 3.045; break;
291 case eMuManganese: return 3.107; break;
292 case eMuIron: return 3.153; break;
293 case eMuNickel: return 3.185; break;
294 case eMuCopper: return 3.279; break;
295 case eMuZinc: return 3.367; break;
296 case eMuGermanium: return 3.610; break;
297 case eMuBromine: return 3.000; break;
298 case eMuTin: return 3.296; break;
299 case eMuIodine: return 3.260; break;
300 case eMuBarium: return 3.455; break;
301 case eMuTungsten: return 3.496; break;
302 case eMuLead: return 3.807; break;
303 case eMuBismuth: return 3.825; break;
304 case eMuUranium: return 3.372; break;
305
306 case eMuUndefined:
307 default:
308 return 0;
309 }
310 return 0;
311 }
312 //____________________________________________________________________
313 static double DensityCorrection_a(MuELMaterial_t material)
314 {
315 // returns the density correction constant a
316
317 switch(material) {
318
319 //- compound materials
320 case eMuBariumFluoride: return 0.160; break;
321 case eMuBismuthGermanate: return 0.096; break;
322 case eMuPyrex: return 0.083; break;
323 case eMuCalciumCarbonate: return 0.083; break;
324 case eMuConcrete: return 0.075; break;
325 case eMuFreon12: return 0.080; break;
326 case eMuFreon13B1: return 0.039; break;
327 case eMuLeadOxide: return 0.196; break;
328 case eMuLithiumFluoride: return 0.076; break;
329 case eMuLucite: return 0.114; break;
330 case eMuPolyethylene: return 0.121; break;
331 case eMuPolystyrene: return 0.165; break;
332 case eMuLiquidPropane: return 0.103; break;
333 case eMuSiliconDioxide: return 0.084; break;
334 case eMuSodiumIodide: return 0.125; break;
335 case eMuStandardRock: return 0.083; break;
336 case eMuUraniumOxide: return 0.205; break;
337 case eMuWater: return 0.091; break;
338
339 //- elements
340 case eMuHydrogen: return 0.135; break;
341 case eMuDeuterium: return 0.347; break;
342 case eMuHelium: return 0.657; break;
343 case eMuLithium: return 0.951; break;
344 case eMuBeryllium: return 0.804; break;
345 case eMuBoron: return 0.562; break;
346 case eMuCarbon: return 0.261; break;
347 case eMuNitrogen: return 0.533; break;
348 case eMuOxygen: return 0.523; break;
349 case eMuFluorine: return 0.159; break;
350 case eMuNeon: return 0.169; break;
351 case eMuSodium: return 0.078; break;
352 case eMuMagnesium: return 0.082; break;
353 case eMuAluminium: return 0.080; break;
354 case eMuSilicon: return 0.146; break;
355 case eMuSulphur: return 0.340; break;
356 case eMuChlorine: return 0.181; break;
357 case eMuArgon: return 0.196; break;
358 case eMuPotassium: return 0.198; break;
359 case eMuCalcium: return 0.156; break;
360 case eMuChromium: return 0.154; break;
361 case eMuManganese: return 0.150; break;
362 case eMuIron: return 0.147; break;
363 case eMuNickel: return 0.165; break;
364 case eMuCopper: return 0.143; break;
365 case eMuZinc: return 0.147; break;
366 case eMuGermanium: return 0.072; break;
367 case eMuBromine: return 0.217; break;
368 case eMuTin: return 0.187; break;
369 case eMuIodine: return 0.238; break;
370 case eMuBarium: return 0.183; break;
371 case eMuTungsten: return 0.155; break;
372 case eMuLead: return 0.094; break;
373 case eMuBismuth: return 0.094; break;
374 case eMuUranium: return 0.197; break;
375
376 case eMuUndefined:
377 default:
378 return 0;
379 }
380 return 0;
381 }
382 //____________________________________________________________________
383 static double DensityCorrection_m(MuELMaterial_t material)
384 {
385 // returns the density correction constant m
386
387 switch(material) {
388
389 //- compound materials
390 case eMuBariumFluoride: return 2.887; break;
391 case eMuBismuthGermanate: return 3.078; break;
392 case eMuPyrex: return 3.522; break;
393 case eMuCalciumCarbonate: return 3.412; break;
394 case eMuConcrete: return 3.547; break;
395 case eMuFreon12: return 3.463; break;
396 case eMuFreon13B1: return 3.719; break;
397 case eMuLeadOxide: return 2.730; break;
398 case eMuLithiumFluoride: return 3.748; break;
399 case eMuLucite: return 3.384; break;
400 case eMuPolyethylene: return 3.429; break;
401 case eMuPolystyrene: return 3.222; break;
402 case eMuLiquidPropane: return 3.562; break;
403 case eMuSiliconDioxide: return 3.506; break;
404 case eMuSodiumIodide: return 3.040; break;
405 case eMuStandardRock: return 3.412; break;
406 case eMuUraniumOxide: return 2.671; break;
407 case eMuWater: return 3.477; break;
408
409 //- elements
410 case eMuHydrogen: return 5.625; break;
411 case eMuDeuterium: return 3.000; break;
412 case eMuHelium: return 3.000; break;
413 case eMuLithium: return 2.499; break;
414 case eMuBeryllium: return 2.434; break;
415 case eMuBoron: return 2.451; break;
416 case eMuCarbon: return 2.870; break;
417 case eMuNitrogen: return 3.000; break;
418 case eMuOxygen: return 3.000; break;
419 case eMuFluorine: return 3.000; break;
420 case eMuNeon: return 3.000; break;
421 case eMuSodium: return 3.645; break;
422 case eMuMagnesium: return 3.617; break;
423 case eMuAluminium: return 3.635; break;
424 case eMuSilicon: return 3.255; break;
425 case eMuSulphur: return 2.646; break;
426 case eMuChlorine: return 3.000; break;
427 case eMuArgon: return 3.000; break;
428 case eMuPotassium: return 2.923; break;
429 case eMuCalcium: return 3.075; break;
430 case eMuChromium: return 2.990; break;
431 case eMuManganese: return 2.980; break;
432 case eMuIron: return 2.963; break;
433 case eMuNickel: return 2.843; break;
434 case eMuCopper: return 2.904; break;
435 case eMuZinc: return 2.865; break;
436 case eMuGermanium: return 3.331; break;
437 case eMuBromine: return 3.000; break;
438 case eMuTin: return 2.858; break;
439 case eMuIodine: return 2.728; break;
440 case eMuBarium: return 2.891; break;
441 case eMuTungsten: return 2.845; break;
442 case eMuLead: return 3.161; break;
443 case eMuBismuth: return 3.167; break;
444 case eMuUranium: return 2.817; break;
445
446 case eMuUndefined:
447 default:
448 return 0;
449 }
450 return 0;
451 }
452 //____________________________________________________________________
453};
454
455} // mueloss namespace
456} // genie namespace
457
458#endif // _BETHE_BLOCH_MATERIAL_PARAMS_H_
genie::mueloss::BetheBlochMaterialParams
Bethe Bloch parameters for various materials.
Definition BetheBlochMaterialParams.h:30
genie::mueloss::BetheBlochMaterialParams::DensityCorrection_a
static double DensityCorrection_a(MuELMaterial_t material)
Definition BetheBlochMaterialParams.h:313
genie::mueloss::BetheBlochMaterialParams::DensityCorrection_m
static double DensityCorrection_m(MuELMaterial_t material)
Definition BetheBlochMaterialParams.h:383
genie::mueloss::BetheBlochMaterialParams::DensityCorrection_X0
static double DensityCorrection_X0(MuELMaterial_t material)
Definition BetheBlochMaterialParams.h:173
genie::mueloss::BetheBlochMaterialParams::DensityCorrection_C
static double DensityCorrection_C(MuELMaterial_t material)
Definition BetheBlochMaterialParams.h:103
genie::mueloss::BetheBlochMaterialParams::IonizationPotential
static double IonizationPotential(MuELMaterial_t material)
Definition BetheBlochMaterialParams.h:33
genie::mueloss::BetheBlochMaterialParams::DensityCorrection_X1
static double DensityCorrection_X1(MuELMaterial_t material)
Definition BetheBlochMaterialParams.h:243
genie::mueloss::MuELMaterial_t
enum genie::mueloss::EMuELMaterial MuELMaterial_t
genie
THE MAIN GENIE PROJECT NAMESPACE
Definition AlgCmp.h:25
Generated on for GENIEGenerator by doxygen 1.14.0