| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 // 23 //
27 // ------------------------------------------- 24 // -------------------------------------------------------------------
28 // 25 //
29 // GEANT4 Class file 26 // GEANT4 Class file
30 // 27 //
31 // 28 //
32 // File name: G4hZiegler1985p 29 // File name: G4hZiegler1985p
33 // 30 //
34 // Author: V.Ivanchenko (Vladimir.Ivanc 31 // Author: V.Ivanchenko (Vladimir.Ivanchenko@cern.ch)
35 // 32 //
36 // Creation date: 20 July 2000 33 // Creation date: 20 July 2000
37 // 34 //
38 // Modifications: 35 // Modifications:
39 // 20/07/2000 V.Ivanchenko First implementati 36 // 20/07/2000 V.Ivanchenko First implementation
40 // 18/09/2000 V.Ivanchenko clean up - all var 37 // 18/09/2000 V.Ivanchenko clean up - all variable are the same as in ICRU
41 // 10/05/2001 V.Ivanchenko Clean up againist 38 // 10/05/2001 V.Ivanchenko Clean up againist Linux compilation with -Wall
42 // 39 //
43 // Class Description: 40 // Class Description:
44 // 41 //
45 // Electronic stopping power parametrised acco 42 // Electronic stopping power parametrised according to
46 // J.F.Ziegler, J.P.Biersack, U.Littmark The S 43 // J.F.Ziegler, J.P.Biersack, U.Littmark The Stoping and
47 // Range of Ions in Solids, Vol.1, Pergamon Pr 44 // Range of Ions in Solids, Vol.1, Pergamon Press, 1985
48 // 45 //
49 // Class Description: End 46 // Class Description: End
50 // 47 //
51 // ------------------------------------------- 48 // -------------------------------------------------------------------
52 // 49 //
53 //....oooOO0OOooo........oooOO0OOooo........oo 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
54 51
55 #include "G4hZiegler1985p.hh" 52 #include "G4hZiegler1985p.hh"
56 <<
57 #include "globals.hh" <<
58 #include "G4SystemOfUnits.hh" <<
59 #include "G4UnitsTable.hh" 53 #include "G4UnitsTable.hh"
>> 54 #include "globals.hh"
60 #include "G4Material.hh" 55 #include "G4Material.hh"
61 56
62 const G4double G4hZiegler1985p::a[92][8] = { << 57 G4double G4hZiegler1985p::a[92][8] = {
63 {0.0091827, 0.0053496, 0.69741, 0.48493, 31 58 {0.0091827, 0.0053496, 0.69741, 0.48493, 316.07 , 1.0143, 9329.3, 0.053989},
64 {0.11393, 0.0051984, 1.0822, 0.39252, 10 59 {0.11393, 0.0051984, 1.0822, 0.39252, 1081.0 , 1.0645, 4068.5, 0.017699},
65 {0.85837, 0.0050147, 1.6044, 0.38844, 13 60 {0.85837, 0.0050147, 1.6044, 0.38844, 1337.3 , 1.047, 2659.2, 0.01898},
66 {0.8781, 0.0051049, 5.4232, 0.2032 , 12 61 {0.8781, 0.0051049, 5.4232, 0.2032 , 1200.6 , 1.0211, 1401.8, 0.038529},
67 {1.4608, 0.0048836, 2.338, 0.44249, 18 62 {1.4608, 0.0048836, 2.338, 0.44249, 1801.3 , 1.0352, 1784.1, 0.02024},
68 {3.2579, 0.0049148, 2.7156, 0.36473, 20 63 {3.2579, 0.0049148, 2.7156, 0.36473, 2092.2 , 1.0291, 2643.6, 0.018237},
69 {0.59674, 0.0050837, 4.2073, 0.30612, 23 64 {0.59674, 0.0050837, 4.2073, 0.30612, 2394.2 , 1.0255, 4892.1, 0.016006},
70 {0.75253, 0.0050314, 4.0824, 0.30067, 24 65 {0.75253, 0.0050314, 4.0824, 0.30067, 2455.8 , 1.0181, 5069.7, 0.017426},
71 {1.226, 0.0051385, 3.2246, 0.32703, 25 66 {1.226, 0.0051385, 3.2246, 0.32703, 2525.0 , 1.0142, 7563.6, 0.019469},
72 {1.0332, 0.0051645, 3.004, 0.33889, 23 67 {1.0332, 0.0051645, 3.004, 0.33889, 2338.6 ,0.99997, 6991.2, 0.021799},
73 68
74 {6.0972, 0.0044292, 3.1929, 0.45763, 13 69 {6.0972, 0.0044292, 3.1929, 0.45763, 1363.3 , 0.95182, 2380.6, 0.081835},
75 {14.013, 0.0043646, 2.2641, 0.36326, 21 70 {14.013, 0.0043646, 2.2641, 0.36326, 2187.4 , 0.99098, 6264.8, 0.0462},
76 {0.039001, 0.0045415, 5.5463, 0.39562, 15 71 {0.039001, 0.0045415, 5.5463, 0.39562, 1589.2 , 0.95316, 816.16, 0.047484},
77 {2.072, 0.0044516, 3.5585, 0.53933, 15 72 {2.072, 0.0044516, 3.5585, 0.53933, 1515.2 , 0.93161, 1790.3, 0.035189},
78 {17.575, 0.0038346, 0.078694,1.2388, 28 73 {17.575, 0.0038346, 0.078694,1.2388, 2806.0 , 0.97284, 1037.6, 0.012879},
79 {16.126, 0.0038315, 0.054164,1.3104, 28 74 {16.126, 0.0038315, 0.054164,1.3104, 2813.3 , 0.96587, 1251.4, 0.011847},
80 {3.217, 0.0044579, 3.6696, 0.5091, 27 75 {3.217, 0.0044579, 3.6696, 0.5091, 2734.6 , 0.96253, 2187.5, 0.016907},
81 {2.0379, 0.0044775, 3.0743, 0.54773, 35 76 {2.0379, 0.0044775, 3.0743, 0.54773, 3505.0 , 0.96575, 1714.0, 0.011701},
82 {0.74171, 0.0043051, 1.1515, 0.95083, 91 77 {0.74171, 0.0043051, 1.1515, 0.95083, 917.21 , 0.8782, 389.93, 0.18926},
83 {9.1316, 0.0043809, 5.4611, 0.31327, 38 78 {9.1316, 0.0043809, 5.4611, 0.31327, 3891.8 , 0.97933, 6267.9, 0.015196},
84 79
85 {7.2247, 0.0043718, 6.1017, 0.37511, 28 80 {7.2247, 0.0043718, 6.1017, 0.37511, 2829.2 , 0.95218, 6376.1, 0.020398},
86 {0.147, 0.0048456, 6.3485, 0.41057, 21 81 {0.147, 0.0048456, 6.3485, 0.41057, 2164.1 , 0.94028, 5292.6, 0.050263},
87 {5.0611, 0.0039867, 2.6174, 0.57957, 22 82 {5.0611, 0.0039867, 2.6174, 0.57957, 2218.9 , 0.92361, 6323.0, 0.025669},
88 {0.53267, 0.0042968, 0.39005, 1.2725, 18 83 {0.53267, 0.0042968, 0.39005, 1.2725, 1872.7 , 0.90776, 64.166, 0.030107},
89 {0.47697, 0.0043038, 0.31452, 1.3289, 19 84 {0.47697, 0.0043038, 0.31452, 1.3289, 1920.5 , 0.90649, 45.576, 0.027469},
90 {0.027426, 0.0035443, 0.031563,2.1755, 19 85 {0.027426, 0.0035443, 0.031563,2.1755, 1919.5 , 0.90099, 23.902, 0.025363},
91 {0.16383, 0.0043042, 0.073454,1.8592, 19 86 {0.16383, 0.0043042, 0.073454,1.8592, 1918.4 , 0.89678, 27.61, 0.023184},
92 {4.2562, 0.0043737, 1.5606, 0.72067, 15 87 {4.2562, 0.0043737, 1.5606, 0.72067, 1546.8 , 0.87958, 302.02, 0.040944},
93 {2.3508, 0.0043237, 2.882, 0.50113, 18 88 {2.3508, 0.0043237, 2.882, 0.50113, 1837.7 , 0.89992, 2377.0, 0.04965},
94 {3.1095, 0.0038455, 0.11477, 1.5037, 21 89 {3.1095, 0.0038455, 0.11477, 1.5037, 2184.7 , 0.89309, 67.306, 0.016588},
95 90
96 {15.322, 0.0040306, 0.65391, 0.67668, 30 91 {15.322, 0.0040306, 0.65391, 0.67668, 3001.7 , 0.92484, 3344.2, 0.016366},
97 {3.6932, 0.0044813, 8.608, 0.27638, 29 92 {3.6932, 0.0044813, 8.608, 0.27638, 2982.7 , 0.9276, 3166.6, 0.030874},
98 {7.1373, 0.0043134, 9.4247, 0.27937, 27 93 {7.1373, 0.0043134, 9.4247, 0.27937, 2725.8 , 0.91597, 3166.1, 0.025008},
99 {4.8979, 0.0042937, 3.7793, 0.50004, 28 94 {4.8979, 0.0042937, 3.7793, 0.50004, 2824.5 , 0.91028, 1282.4, 0.017061},
100 {1.3683, 0.0043024, 2.5679, 0.60822, 69 95 {1.3683, 0.0043024, 2.5679, 0.60822, 6907.8 , 0.9817, 628.01, 0.0068055},
101 {1.8301, 0.0042983, 2.9057, 0.6038, 47 96 {1.8301, 0.0042983, 2.9057, 0.6038, 4744.6 , 0.94722, 936.64, 0.0092242},
102 {0.42056, 0.0041169, 0.01695, 2.3616, 22 97 {0.42056, 0.0041169, 0.01695, 2.3616, 2252.7 , 0.89192, 39.752, 0.027757},
103 {30.78, 0.0037736, 0.55813, 0.76816, 71 98 {30.78, 0.0037736, 0.55813, 0.76816, 7113.2 , 0.97697, 1604.4, 0.0065268},
104 {11.576, 0.0042119, 7.0244, 0.37764, 47 99 {11.576, 0.0042119, 7.0244, 0.37764, 4713.5 , 0.94264, 2493.2, 0.01127},
105 {6.2406, 0.0041916, 5.2701, 0.49453, 42 100 {6.2406, 0.0041916, 5.2701, 0.49453, 4234.6 , 0.93232, 2063.9, 0.011844},
106 101
107 {0.33073, 0.0041243, 1.7246, 1.1062, 19 102 {0.33073, 0.0041243, 1.7246, 1.1062, 1930.2 , 0.86907, 27.416, 0.038208},
108 {0.017747, 0.0041715, 0.14586, 1.7305, 18 103 {0.017747, 0.0041715, 0.14586, 1.7305, 1803.6 , 0.86315, 29.669, 0.032123},
109 {3.7229, 0.0041768, 4.6286, 0.56769, 16 104 {3.7229, 0.0041768, 4.6286, 0.56769, 1678.0 , 0.86202, 3094.0, 0.06244},
110 {0.13998, 0.0041329, 0.25573, 1.4241, 19 105 {0.13998, 0.0041329, 0.25573, 1.4241, 1919.3 , 0.86326, 72.797, 0.032235},
111 {0.2859, 0.0041386, 0.31301, 1.3424, 19 106 {0.2859, 0.0041386, 0.31301, 1.3424, 1954.8 , 0.86175, 115.18, 0.029342},
112 {0.76002, 0.0042179, 3.386, 0.76285, 18 107 {0.76002, 0.0042179, 3.386, 0.76285, 1867.4 , 0.85805, 69.994, 0.036448},
113 {6.3957, 0.0041935, 5.4689, 0.41378, 17 108 {6.3957, 0.0041935, 5.4689, 0.41378, 1712.6 , 0.85397, 18493., 0.056471},
114 {3.4717, 0.0041344, 3.2337, 0.63788, 11 109 {3.4717, 0.0041344, 3.2337, 0.63788, 1116.4 , 0.81959, 4766.0, 0.1179},
115 {2.5265, 0.0042282, 4.532, 0.53562, 10 110 {2.5265, 0.0042282, 4.532, 0.53562, 1030.8 , 0.81652, 16252., 0.19722},
116 {7.3683, 0.0041007, 4.6791, 0.51428, 11 111 {7.3683, 0.0041007, 4.6791, 0.51428, 1160.0 , 0.82454, 17956., 0.13316},
117 112
118 {7.7197, 0.004388, 3.242, 0.68434, 14 113 {7.7197, 0.004388, 3.242, 0.68434, 1428.1 , 0.83389, 1786.7, 0.066512},
119 {16.78, 0.0041918, 9.3198, 0.29569, 33 114 {16.78, 0.0041918, 9.3198, 0.29569, 3370.9 , 0.90298, 7431.7, 0.02616},
120 {4.2132, 0.0042098, 4.6753, 0.57945, 35 115 {4.2132, 0.0042098, 4.6753, 0.57945, 3503.9 , 0.89261, 1468.9, 0.014359},
121 {4.0818, 0.004214, 4.4425, 0.58393, 39 116 {4.0818, 0.004214, 4.4425, 0.58393, 3945.3 , 0.90281, 1340.5, 0.013414},
122 {0.18517, 0.0036215,0.00058788,3.5315, 29 117 {0.18517, 0.0036215,0.00058788,3.5315, 2931.3 , 0.88936, 26.18, 0.026393},
123 {4.8248, 0.0041458, 6.0934, 0.57026, 23 118 {4.8248, 0.0041458, 6.0934, 0.57026, 2300.1 , 0.86359, 2980.7, 0.038679},
124 {0.49857, 0.0041054, 1.9775, 0.95877, 78 119 {0.49857, 0.0041054, 1.9775, 0.95877, 786.55 , 0.78509, 806.6, 0.40882},
125 {3.2754, 0.0042177, 5.768, 0.54054, 66 120 {3.2754, 0.0042177, 5.768, 0.54054, 6631.3 , 0.94282, 744.07, 0.0083026},
126 {2.9978, 0.0040901, 4.5299, 0.62025, 21 121 {2.9978, 0.0040901, 4.5299, 0.62025, 2161.2 , 0.85669, 1268.6, 0.043031},
127 {2.8701, 0.004096, 4.2568, 0.6138, 21 122 {2.8701, 0.004096, 4.2568, 0.6138, 2130.4 , 0.85235, 1704.1, 0.039385},
128 123
129 {10.853, 0.0041149, 5.8906, 0.46834, 28 124 {10.853, 0.0041149, 5.8906, 0.46834, 2857.2 , 0.87550, 3654.2, 0.029955},
130 {3.6407, 0.0041782, 4.8742, 0.57861, 12 125 {3.6407, 0.0041782, 4.8742, 0.57861, 1267.7 , 0.82211, 3508.2, 0.24174},
131 {17.645, 0.0040992, 6.5855, 0.32734, 39 126 {17.645, 0.0040992, 6.5855, 0.32734, 3931.3 , 0.90754, 5156.7, 0.036278},
132 {7.5309, 0.0040814, 4.9389, 0.50679, 25 127 {7.5309, 0.0040814, 4.9389, 0.50679, 2519.7 , 0.85819, 3314.6, 0.030514},
133 {5.4742, 0.0040829, 4.897, 0.51113, 23 128 {5.4742, 0.0040829, 4.897, 0.51113, 2340.1 , 0.85296, 2342.7, 0.035662},
134 {4.2661, 0.0040667, 4.5032, 0.55257, 20 129 {4.2661, 0.0040667, 4.5032, 0.55257, 2076.4 , 0.84151, 1666.6, 0.040801},
135 {6.8313, 0.0040486, 4.3987, 0.51675, 20 130 {6.8313, 0.0040486, 4.3987, 0.51675, 2003.0 , 0.83437, 1410.4, 0.03478},
136 {1.2707, 0.0040553, 4.6295, 0.57428, 16 131 {1.2707, 0.0040553, 4.6295, 0.57428, 1626.3 , 0.81858, 995.68, 0.055319},
137 {5.7561, 0.0040491, 4.357, 0.52496, 22 132 {5.7561, 0.0040491, 4.357, 0.52496, 2207.3 , 0.83796, 1579.5, 0.027165},
138 {14.127, 0.0040596, 5.8304, 0.37755, 36 133 {14.127, 0.0040596, 5.8304, 0.37755, 3645.9 , 0.87823, 3411.8, 0.016392},
139 134
140 {6.6948, 0.0040603, 4.9361, 0.47961, 27 135 {6.6948, 0.0040603, 4.9361, 0.47961, 2719.0 , 0.85249, 1885.8, 0.019713},
141 {3.0619, 0.0040511, 3.5803, 0.59082, 23 136 {3.0619, 0.0040511, 3.5803, 0.59082, 2346.1 , 0.83713, 1222.0, 0.020072},
142 {10.811, 0.0033008, 1.3767, 0.76512, 20 137 {10.811, 0.0033008, 1.3767, 0.76512, 2003.7 , 0.82269, 1110.6, 0.024958},
143 {2.7101, 0.0040961, 1.2289, 0.98598, 12 138 {2.7101, 0.0040961, 1.2289, 0.98598, 1232.4 , 0.79066, 155.42, 0.047294},
144 {0.52345, 0.0040244, 1.4038, 0.8551, 14 139 {0.52345, 0.0040244, 1.4038, 0.8551, 1461.4 , 0.79677, 503.34, 0.036789},
145 {0.4616, 0.0040203, 1.3014, 0.87043, 14 140 {0.4616, 0.0040203, 1.3014, 0.87043, 1473.5 , 0.79687, 443.09, 0.036301},
146 {0.97814, 0.0040374, 2.0127, 0.7225, 18 141 {0.97814, 0.0040374, 2.0127, 0.7225, 1890.8 , 0.81747, 930.7, 0.02769},
147 {3.2086, 0.0040510, 3.6658, 0.53618, 30 142 {3.2086, 0.0040510, 3.6658, 0.53618, 3091.2 , 0.85602, 1508.1, 0.015401},
148 {2.0035, 0.0040431, 7.4882, 0.3561, 44 143 {2.0035, 0.0040431, 7.4882, 0.3561, 4464.3 , 0.88836, 3966.5, 0.012839},
149 {15.43, 0.0039432, 1.1237, 0.70703, 45 144 {15.43, 0.0039432, 1.1237, 0.70703, 4595.7 , 0.88437, 1576.5, 0.0088534},
150 145
151 {3.1512, 0.0040524, 4.0996, 0.5425, 32 146 {3.1512, 0.0040524, 4.0996, 0.5425, 3246.3 , 0.85772, 1691.8, 0.015058},
152 {7.1896, 0.0040588, 8.6927, 0.35842, 47 147 {7.1896, 0.0040588, 8.6927, 0.35842, 4760.6 , 0.88833, 2888.3, 0.011029},
153 {9.3209, 0.0040540, 11.543, 0.32027, 48 148 {9.3209, 0.0040540, 11.543, 0.32027, 4866.2 , 0.89124, 3213.4, 0.011935},
154 {29.242, 0.0036195, 0.16864, 1.1226, 56 149 {29.242, 0.0036195, 0.16864, 1.1226, 5688.0 , 0.89812, 1033.3, 0.0071303},
155 {1.8522, 0.0039973, 3.1556, 0.65096, 37 150 {1.8522, 0.0039973, 3.1556, 0.65096, 3755.0 , 0.86383, 1602.0, 0.012042},
156 {3.222, 0.0040041, 5.9024, 0.52678, 40 151 {3.222, 0.0040041, 5.9024, 0.52678, 4040.2 , 0.86804, 1658.4, 0.011747},
157 {9.3412, 0.0039661, 7.921, 0.42977, 51 152 {9.3412, 0.0039661, 7.921, 0.42977, 5180.9 , 0.88773, 2173.2, 0.0092007},
158 {36.183, 0.0036003, 0.58341, 0.86747, 69 153 {36.183, 0.0036003, 0.58341, 0.86747, 6990.2 , 0.91082, 1417.1, 0.0062187},
159 {5.9284, 0.0039695, 6.4082, 0.52122, 46 154 {5.9284, 0.0039695, 6.4082, 0.52122, 4619.5 , 0.88083, 2323.5, 0.011627},
160 {5.2454, 0.0039744, 6.7969, 0.48542, 45 155 {5.2454, 0.0039744, 6.7969, 0.48542, 4586.3 , 0.87794, 2481.5, 0.011282},
161 156
162 {33.702, 0.0036901, 0.47257, 0.89235, 52 157 {33.702, 0.0036901, 0.47257, 0.89235, 5295.7 , 0.8893, 2053.3, 0.0091908},
163 {2.7589, 0.0039806, 3.2092, 0.66122, 25 158 {2.7589, 0.0039806, 3.2092, 0.66122, 2505.4 , 0.82863, 2065.1, 0.022816}
164 }; 159 };
165 160
166 //....oooOO0OOooo........oooOO0OOooo........oo 161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
167 162
168 G4hZiegler1985p::G4hZiegler1985p():G4VhElectro 163 G4hZiegler1985p::G4hZiegler1985p():G4VhElectronicStoppingPower(),
169 protonMassAMU(1.007276) 164 protonMassAMU(1.007276)
170 {;} 165 {;}
171 166
172 //....oooOO0OOooo........oooOO0OOooo........oo 167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
173 168
174 G4hZiegler1985p::~G4hZiegler1985p() 169 G4hZiegler1985p::~G4hZiegler1985p()
175 {;} 170 {;}
176 171
177 //....oooOO0OOooo........oooOO0OOooo........oo 172 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
178 173
179 G4bool G4hZiegler1985p::HasMaterial(const G4Ma 174 G4bool G4hZiegler1985p::HasMaterial(const G4Material* material)
180 { 175 {
181 if(1 == (material->GetNumberOfElements())) r 176 if(1 == (material->GetNumberOfElements())) return true;
182 return false ; 177 return false ;
183 } 178 }
184 179
185 //....oooOO0OOooo........oooOO0OOooo........oo 180 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
186 181
187 G4double G4hZiegler1985p::StoppingPower(const 182 G4double G4hZiegler1985p::StoppingPower(const G4Material* material,
188 183 G4double kineticEnergy)
189 { 184 {
190 G4double ionloss = 0.0 ; 185 G4double ionloss = 0.0 ;
191 186
192 // pure material (normally not the case for 187 // pure material (normally not the case for this function)
193 if(1 == (material->GetNumberOfElements())) { 188 if(1 == (material->GetNumberOfElements())) {
194 G4double z = material->GetZ() ; 189 G4double z = material->GetZ() ;
195 ionloss = ElectronicStoppingPower( z, kine 190 ionloss = ElectronicStoppingPower( z, kineticEnergy ) ;
196 } 191 }
197 192
198 return ionloss; 193 return ionloss;
199 } 194 }
200 195
201 //....oooOO0OOooo........oooOO0OOooo........oo 196 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
202 197
203 G4double G4hZiegler1985p::ElectronicStoppingPo 198 G4double G4hZiegler1985p::ElectronicStoppingPower(G4double z,
204 199 G4double kineticEnergy) const
205 { 200 {
206 G4double ionloss ; 201 G4double ionloss ;
207 G4int i = G4int(z) - 1 ; // index of atom 202 G4int i = G4int(z) - 1 ; // index of atom
208 if(i < 0) i = 0 ; 203 if(i < 0) i = 0 ;
209 if(i > 91) i = 91 ; 204 if(i > 91) i = 91 ;
210 205
211 // The data and the fit from: 206 // The data and the fit from:
212 // J.F.Ziegler, J.P.Biersack, U.Littmark The 207 // J.F.Ziegler, J.P.Biersack, U.Littmark The Stoping and
213 // Range of Ions in Solids, Vol.1, Pergamon 208 // Range of Ions in Solids, Vol.1, Pergamon Press, 1985
214 // Proton kinetic energy for parametrisation 209 // Proton kinetic energy for parametrisation in Ziegler's units (keV/amu)
215 210
216 G4double T = kineticEnergy/(keV*protonMassAM 211 G4double T = kineticEnergy/(keV*protonMassAMU) ;
217 212
218 213
219 214
220 G4double e = T ; 215 G4double e = T ;
221 if ( T < 25.0 ) e = 25.0 ; 216 if ( T < 25.0 ) e = 25.0 ;
222 217
223 // universal approximation 218 // universal approximation
224 G4double slow = a[i][0] * std::pow(e, a[i][ << 219 G4double slow = a[i][0] * pow(e, a[i][1]) + a[i][2] * pow(e, a[i][3]) ;
225 G4double shigh = std::log( a[i][6]/e + a[i][ << 220 G4double shigh = log( a[i][6]/e + a[i][7]*e ) * a[i][4] / pow(e, a[i][5]) ;
226 ionloss = slow*shigh / (slow + shigh) ; 221 ionloss = slow*shigh / (slow + shigh) ;
227 222
228 // low energy region 223 // low energy region
229 if ( T < 25.0 ) { 224 if ( T < 25.0 ) {
230 225
231 G4double sLocal = 0.45 ; << 226 G4double s = 0.45 ;
232 // light elements 227 // light elements
233 if(6.5 > z) sLocal = 0.25 ; << 228 if(6.5 > z) s = 0.25 ;
234 // semiconductors 229 // semiconductors
235 if(5 == i || 13 == i || 31 == i) sLocal = << 230 if(5 == i || 13 == i || 31 == i) s = 0.375 ;
236 231
237 ionloss *= std::pow(T/25.0, sLocal) ; << 232 ionloss *= pow(T/25.0, s) ;
238 } 233 }
239 234
240 if ( ionloss < 0.0) ionloss = 0.0 ; 235 if ( ionloss < 0.0) ionloss = 0.0 ;
241 236
242 return ionloss; 237 return ionloss;
243 } 238 }
244 239
245 240