| 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 //-------------------------------------------- 27 //---------------------------------------------------------------------------
28 // 28 //
29 // GEANT4 Class file 29 // GEANT4 Class file
30 // 30 //
31 // Description: Data on stopping power 31 // Description: Data on stopping power
32 // 32 //
33 // Description: Data on stopping power 33 // Description: Data on stopping power
34 // 34 //
35 // Author: Alexander Bagulya & Vladimir 35 // Author: Alexander Bagulya & Vladimir Ivanchenko
36 // 36 //
37 // Creation date: 23.04.2018 37 // Creation date: 23.04.2018
38 // 38 //
39 //-------------------------------------------- 39 //----------------------------------------------------------------------------
40 // 40 //
41 41
42 //....oooOO0OOooo........oooOO0OOooo........oo 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 43
44 #include "G4LindhardSorensenData.hh" 44 #include "G4LindhardSorensenData.hh"
45 #include "G4PhysicsVector.hh" 45 #include "G4PhysicsVector.hh"
46 #include "G4PhysicsLinearVector.hh" 46 #include "G4PhysicsLinearVector.hh"
47 #include "G4Log.hh" 47 #include "G4Log.hh"
48 #include "G4Pow.hh" 48 #include "G4Pow.hh"
49 49
50 const G4int zlist[9] = {1, 10, 18, 36, 54, 66, 50 const G4int zlist[9] = {1, 10, 18, 36, 54, 66, 79, 92, 109};
51 const G4int LVECT = 8; 51 const G4int LVECT = 8;
52 const G4int NPOINT = 41; 52 const G4int NPOINT = 41;
53 53
54 //....oooOO0OOooo........oooOO0OOooo........oo 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55 55
56 G4LindhardSorensenData::G4LindhardSorensenData 56 G4LindhardSorensenData::G4LindhardSorensenData()
57 { 57 {
58 g4calc = G4Pow::GetInstance(); 58 g4calc = G4Pow::GetInstance();
59 InitialiseData(); 59 InitialiseData();
60 } 60 }
61 61
62 //....oooOO0OOooo........oooOO0OOooo........oo 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
63 63
64 G4LindhardSorensenData::~G4LindhardSorensenDat 64 G4LindhardSorensenData::~G4LindhardSorensenData()
65 { 65 {
66 for(G4int i=0; i<=LVECT; ++i) { delete data[ 66 for(G4int i=0; i<=LVECT; ++i) { delete data[i]; }
67 } 67 }
68 68
69 //....oooOO0OOooo........oooOO0OOooo........oo 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
70 70
71 G4double G4LindhardSorensenData::GetDeltaL(G4i 71 G4double G4LindhardSorensenData::GetDeltaL(G4int Z, G4double gamma) const
72 { 72 {
73 G4int idx = 0; 73 G4int idx = 0;
74 for(; idx<LVECT; ++idx) { 74 for(; idx<LVECT; ++idx) {
75 if(Z < zlist[idx+1]) { break; } 75 if(Z < zlist[idx+1]) { break; }
76 } 76 }
77 idx = std::min(idx, LVECT); 77 idx = std::min(idx, LVECT);
78 78
79 G4double x = G4Log(gamma - 1.0); 79 G4double x = G4Log(gamma - 1.0);
80 G4double y = ComputeDeltaL(idx, x); 80 G4double y = ComputeDeltaL(idx, x);
81 81
82 // interpolation over Z if needed 82 // interpolation over Z if needed
83 if(idx < LVECT && Z > zlist[idx]) { 83 if(idx < LVECT && Z > zlist[idx]) {
84 G4double y1 = ComputeDeltaL(idx+1, x); 84 G4double y1 = ComputeDeltaL(idx+1, x);
85 //G4cout << "idx= " << idx << " x= " << x 85 //G4cout << "idx= " << idx << " x= " << x << " y= " << y << " y1= " << y1 << G4endl;
86 y += (y1 - y)*(Z - zlist[idx])/(G4double)( 86 y += (y1 - y)*(Z - zlist[idx])/(G4double)(zlist[idx+1] - zlist[idx]);
87 } 87 }
88 y *= g4calc->Z23(Z); 88 y *= g4calc->Z23(Z);
89 return y; 89 return y;
90 } 90 }
91 91
92 //....oooOO0OOooo........oooOO0OOooo........oo 92 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
93 93
94 G4double G4LindhardSorensenData::ComputeDeltaL 94 G4double G4LindhardSorensenData::ComputeDeltaL(G4int idx, G4double x) const
95 { 95 {
96 G4double y(0.0); 96 G4double y(0.0);
97 if(x < xmin) { 97 if(x < xmin) {
98 G4double x1 = (data[idx])->Energy(1); 98 G4double x1 = (data[idx])->Energy(1);
99 G4double ymin = (*(data[idx]))[0]; 99 G4double ymin = (*(data[idx]))[0];
100 G4double y1 = (*(data[idx]))[1]; 100 G4double y1 = (*(data[idx]))[1];
101 y = ymin + (y1 - ymin)*(x - xmin)/(x1 - xm 101 y = ymin + (y1 - ymin)*(x - xmin)/(x1 - xmin);
102 } else if(x >= xmax) { 102 } else if(x >= xmax) {
103 G4double x1 = (data[idx])->Energy(NPOINT-2 103 G4double x1 = (data[idx])->Energy(NPOINT-2);
104 G4double ymax = (*(data[idx]))[NPOINT-1]; 104 G4double ymax = (*(data[idx]))[NPOINT-1];
105 G4double y1 = (*(data[idx]))[NPOINT-2]; 105 G4double y1 = (*(data[idx]))[NPOINT-2];
106 y = y1 + (ymax - y1)*(x - x1)/(xmax - x1); 106 y = y1 + (ymax - y1)*(x - x1)/(xmax - x1);
107 } else { 107 } else {
108 y = (data[idx])->Value(x); 108 y = (data[idx])->Value(x);
109 } 109 }
110 return y; 110 return y;
111 } 111 }
112 112
113 //....oooOO0OOooo........oooOO0OOooo........oo 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
114 114
115 void G4LindhardSorensenData::InitialiseData() 115 void G4LindhardSorensenData::InitialiseData()
116 { 116 {
117 xmin = G4Log(0.02); 117 xmin = G4Log(0.02);
118 xmax = G4Log(316.22777); 118 xmax = G4Log(316.22777);
119 119
120 const G4double lsdata[9][NPOINT] = { 120 const G4double lsdata[9][NPOINT] = {
121 {0.0036181621, 121 {0.0036181621,
122 0.0042618872, 0.0040786701, 0.0039097273, 122 0.0042618872, 0.0040786701, 0.0039097273, 0.0041240731, 0.0044311195, // 0-5
123 0.0059858073, 0.0079655897, 0.0089669217, 123 0.0059858073, 0.0079655897, 0.0089669217, 0.010091248, 0.0096319233, // 5-10
124 0.0085344428, 0.0087326058, 0.009780094, 124 0.0085344428, 0.0087326058, 0.009780094, 0.010617094, 0.011277997, // 10-15
125 0.01185287, 0.012242278, 0.012507042, 125 0.01185287, 0.012242278, 0.012507042, 0.012369698, 0.012202436, // 15-20
126 0.012016446, 0.011793182, 0.011612196, 126 0.012016446, 0.011793182, 0.011612196, 0.011375017, 0.010624408, // 20-25
127 0.0096685612, 0.0093765113, 0.0091152632, 127 0.0096685612, 0.0093765113, 0.0091152632, 0.0087458522, 0.0082645153, // 25-30
128 0.0068625676, 0.0046558921, 0.0021660968, 128 0.0068625676, 0.0046558921, 0.0021660968, -0.00072760644, -0.0062653709,// 30-35
129 -0.015834368, -0.03224412, -0.062296044, 129 -0.015834368, -0.03224412, -0.062296044, -0.10327705, -0.16545368},
130 {-0.024837796, 130 {-0.024837796,
131 -0.018943357, -0.013251703, -0.007646936, 131 -0.018943357, -0.013251703, -0.007646936, -0.003751624, 7.4932758e-05, // 0-5
132 0.0027996278, 0.0057534909, 0.0078441157, 132 0.0027996278, 0.0057534909, 0.0078441157, 0.010341958, 0.012419648, // 5-10
133 0.01491866, 0.016403769, 0.018106955, 133 0.01491866, 0.016403769, 0.018106955, 0.019249728, 0.020704966, // 10-15
134 0.021636244, 0.022526456, 0.023429501, 134 0.021636244, 0.022526456, 0.023429501, 0.024251631, 0.024727856, // 15-20
135 0.024997337, 0.025225657, 0.025335138, 135 0.024997337, 0.025225657, 0.025335138, 0.025474553, 0.025306356, // 20-25
136 0.025048151, 0.02467553, 0.02416496, 136 0.025048151, 0.02467553, 0.02416496, 0.023519917, 0.022699646, // 25-30
137 0.020902526, 0.017645158, 0.013219058, 137 0.020902526, 0.017645158, 0.013219058, 0.0059989118, -0.0042559395, // 30-35
138 -0.020223315, -0.044480728, -0.074942757, 138 -0.020223315, -0.044480728, -0.074942757, -0.1108863, -0.15439805},
139 {-0.04948514, 139 {-0.04948514,
140 -0.039259418, -0.029641846, -0.022112859, 140 -0.039259418, -0.029641846, -0.022112859, -0.0152972, -0.0077213168, // 0-5
141 -0.0031032447, 0.0023506153, 0.0058573025 141 -0.0031032447, 0.0023506153, 0.0058573025, 0.0097554723, 0.012692892, // 5-10
142 0.016404575, 0.018642388, 0.021492082, 142 0.016404575, 0.018642388, 0.021492082, 0.023328047, 0.025358859, // 10-15
143 0.026839018, 0.028169751, 0.029240668, 143 0.026839018, 0.028169751, 0.029240668, 0.030035325, 0.030689491, // 15-20
144 0.031113343, 0.031518392, 0.031657748, 144 0.031113343, 0.031518392, 0.031657748, 0.031794607, 0.031672646, // 20-25
145 0.031507453, 0.031024952, 0.030410522, 145 0.031507453, 0.031024952, 0.030410522, 0.029266343, 0.027729316, // 25-30
146 0.025132583, 0.02141173, 0.016127626, 146 0.025132583, 0.02141173, 0.016127626, 0.008039769, -0.0046793613, // 30-35
147 -0.02196456, -0.043877567, -0.070852217, 147 -0.02196456, -0.043877567, -0.070852217, -0.09979099, -0.13187648},
148 {-0.07502957, 148 {-0.07502957,
149 -0.063927817, -0.053895181, -0.043621797, 149 -0.063927817, -0.053895181, -0.043621797, -0.033329203, -0.024765223, // 0-5
150 -0.015996122, -0.0084995741, -0.000174992 150 -0.015996122, -0.0084995741, -0.00017499271, 0.0064195209, 0.01117532,// 5-10
151 0.016956954, 0.02155737, 0.025213602, 151 0.016956954, 0.02155737, 0.025213602, 0.029869549, 0.031980316, // 10-15
152 0.034664781, 0.036995048, 0.038244843, 152 0.034664781, 0.036995048, 0.038244843, 0.039836367, 0.040930356, // 15-20
153 0.041677336, 0.042211864, 0.042472572, 153 0.041677336, 0.042211864, 0.042472572, 0.042617755, 0.042393662, // 20-25
154 0.042049893, 0.041341375, 0.040408112, 154 0.042049893, 0.041341375, 0.040408112, 0.03867031, 0.036457345, // 25-30
155 0.032742289, 0.027897265, 0.020305954, 155 0.032742289, 0.027897265, 0.020305954, 0.0099602019, -0.0036658833, // 30-35
156 -0.020344557, -0.039700638, -0.059781744, 156 -0.020344557, -0.039700638, -0.059781744, -0.080649004, -0.10263347},
157 {-0.081885964, 157 {-0.081885964,
158 -0.072049323, -0.06299877, -0.053400445, 158 -0.072049323, -0.06299877, -0.053400445, -0.043830227, -0.03415406, // 0-5
159 -0.024723298, -0.015706543, -0.0070305988 159 -0.024723298, -0.015706543, -0.0070305988, 0.0016621818, 0.0093560784, // 5-10
160 0.017338568, 0.023286722, 0.028930541, 160 0.017338568, 0.023286722, 0.028930541, 0.034048359, 0.037750635, // 10-15
161 0.041512165, 0.04362858, 0.046323682, 161 0.041512165, 0.04362858, 0.046323682, 0.048044285, 0.049379817, // 15-20
162 0.050501289, 0.050963981, 0.051463496, 162 0.050501289, 0.050963981, 0.051463496, 0.051292532, 0.051074821, // 20-25
163 0.050350498, 0.049370817, 0.047660517, 163 0.050350498, 0.049370817, 0.047660517, 0.045247985, 0.04188952, // 25-30
164 0.037009647, 0.030171627, 0.021186013, 164 0.037009647, 0.030171627, 0.021186013, 0.0092064517, -0.0051358689, // 30-35
165 -0.01977885, -0.036354277, -0.053187271, 165 -0.01977885, -0.036354277, -0.053187271, -0.069851489, -0.087776477},
166 {-0.083068958, 166 {-0.083068958,
167 -0.074302116, -0.065040071, -0.056500699, 167 -0.074302116, -0.065040071, -0.056500699, -0.047189921, -0.03781489, // 0-5
168 -0.028666422, -0.0192974, -0.0097949279 168 -0.028666422, -0.0192974, -0.0097949279, -0.0011649946, 0.0071301285, // 5-10
169 0.01548354, 0.023483416, 0.030515545, 169 0.01548354, 0.023483416, 0.030515545, 0.036162102, 0.040602642, // 10-15
170 0.045011738, 0.048259036, 0.051136495, 170 0.045011738, 0.048259036, 0.051136495, 0.053524791, 0.054971064, // 15-20
171 0.056247898, 0.056729496, 0.057135059, 171 0.056247898, 0.056729496, 0.057135059, 0.056923775, 0.05657171, // 20-25
172 0.055504504, 0.054134086, 0.051812977, 172 0.055504504, 0.054134086, 0.051812977, 0.048857213, 0.044525663, // 25-30
173 0.038593442, 0.03027527, 0.020452704, 173 0.038593442, 0.03027527, 0.020452704, 0.0082978722, -0.0045556908, // 30-35
174 -0.018750622, -0.033642805, -0.048788529, 174 -0.018750622, -0.033642805, -0.048788529, -0.064509221, -0.080215679},
175 {-0.083061344, 175 {-0.083061344,
176 -0.075180377, -0.066630961, -0.058114277, 176 -0.075180377, -0.066630961, -0.058114277, -0.049303999, -0.040618218, // 0-5
177 -0.031608369, -0.022553882, -0.012699373, 177 -0.031608369, -0.022553882, -0.012699373, -0.0033434259, 0.0051129584, // 5-10
178 0.013817276, 0.023132109, 0.030886148, 178 0.013817276, 0.023132109, 0.030886148, 0.037688311, 0.043478106, // 10-15
179 0.048602814, 0.052502168, 0.055969817, 179 0.048602814, 0.052502168, 0.055969817, 0.058142903, 0.060860174, // 15-20
180 0.061732795, 0.062844014, 0.062911637, 180 0.061732795, 0.062844014, 0.062911637, 0.06278651, 0.061898673, // 20-25
181 0.060500452, 0.058323088, 0.055297604, 181 0.060500452, 0.058323088, 0.055297604, 0.051242701, 0.04596475, // 25-30
182 0.038863441, 0.030302781, 0.019055839, 182 0.038863441, 0.030302781, 0.019055839, 0.0072837125, -0.0047516889, // 30-35
183 -0.017733076, -0.031547911, -0.04523347, 183 -0.017733076, -0.031547911, -0.04523347, -0.059727542, -0.074202114},
184 {-0.082449782, 184 {-0.082449782,
185 -0.07492693, -0.067041495, -0.059614338, 185 -0.07492693, -0.067041495, -0.059614338, -0.051061455, -0.042587329, // 0-5
186 -0.034242658, -0.025296495, -0.016157818, 186 -0.034242658, -0.025296495, -0.016157818, -0.0064198954, 0.0033148315, // 5-10
187 0.012927383, 0.021914897, 0.03053429, 187 0.012927383, 0.021914897, 0.03053429, 0.03854699, 0.045293871, // 10-15
188 0.05141546, 0.056473799, 0.060050334, 188 0.05141546, 0.056473799, 0.060050334, 0.063292547, 0.065215264, // 15-20
189 0.066961692, 0.067369593, 0.067623539, 189 0.066961692, 0.067369593, 0.067623539, 0.066825274, 0.065740573, // 20-25
190 0.063673411, 0.061033708, 0.057399807, 190 0.063673411, 0.061033708, 0.057399807, 0.052544441, 0.04583072, // 25-30
191 0.037931379, 0.028569029, 0.017544895, 191 0.037931379, 0.028569029, 0.017544895, 0.0068751554, -0.0038105519, // 30-35
192 -0.015233268, -0.026956839, -0.038975048, 192 -0.015233268, -0.026956839, -0.038975048, -0.051387898, -0.064344384},
193 {-0.081232852, 193 {-0.081232852,
194 -0.074976912, -0.068421001, -0.060951614, 194 -0.074976912, -0.068421001, -0.060951614, -0.053624488, -0.045755024, // 0-5
195 -0.037708151, -0.029246587, -0.02018635, 195 -0.037708151, -0.029246587, -0.02018635, -0.01073862, -0.00085387172, // 5-10
196 0.0090831897, 0.018970381, 0.02761289, 196 0.0090831897, 0.018970381, 0.02761289, 0.036215831, 0.044355404, // 10-15
197 0.051067631, 0.057143777, 0.061348214, 197 0.051067631, 0.057143777, 0.061348214, 0.065462308, 0.06715594, // 15-20
198 0.069155687, 0.069845259, 0.069828188, 198 0.069155687, 0.069845259, 0.069828188, 0.069170195, 0.067585183, // 20-25
199 0.065486015, 0.062437523, 0.058111347, 199 0.065486015, 0.062437523, 0.058111347, 0.052410482, 0.045597673, // 25-30
200 0.037571853, 0.028019555, 0.016564385, 200 0.037571853, 0.028019555, 0.016564385, 0.0060254421, -0.0043254115, // 30-35
201 -0.015675051, -0.027353574, -0.03931666, 201 -0.015675051, -0.027353574, -0.03931666, -0.050325729, -0.06174447}};
202 202
203 for(G4int i=0; i<=LVECT; ++i) { 203 for(G4int i=0; i<=LVECT; ++i) {
204 data[i] = new G4PhysicsLinearVector(xmin, << 204 data[i] = new G4PhysicsLinearVector(xmin, xmax, NPOINT-1);
>> 205 data[i]->SetSpline(true);
205 for(std::size_t j=0; j<NPOINT; ++j) { 206 for(std::size_t j=0; j<NPOINT; ++j) {
206 data[i]->PutValue(j, lsdata[i][j]); 207 data[i]->PutValue(j, lsdata[i][j]);
207 } 208 }
208 data[i]->FillSecondDerivatives(); <<
209 } 209 }
210 } 210 }
211 211
212 //....oooOO0OOooo........oooOO0OOooo........oo 212 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
213 213