GENIEGenerator
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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-2024, 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 
24 namespace genie {
25 namespace mueloss {
26 
28 
30 {
31 public:
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_
Bethe Bloch parameters for various materials.
static double DensityCorrection_m(MuELMaterial_t material)
static double DensityCorrection_C(MuELMaterial_t material)
static double DensityCorrection_a(MuELMaterial_t material)
static double IonizationPotential(MuELMaterial_t material)
static double DensityCorrection_X0(MuELMaterial_t material)
static double DensityCorrection_X1(MuELMaterial_t material)
enum genie::mueloss::EMuELMaterial MuELMaterial_t