Belle II Software development
BkgNeutronWeight.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9#include <simulation/background/BkgNeutronWeight.h>
10#include <cmath>
11
12namespace Belle2 {
18 BkgNeutronWeight* BkgNeutronWeight::m_Instance = 0;
19
20 BkgNeutronWeight& BkgNeutronWeight::getInstance()
21 {
22 if (m_Instance == 0) m_Instance = new BkgNeutronWeight();
23 return *m_Instance;
24 }
25
26 double BkgNeutronWeight::getWeight(double ke)
27 {
28 if (m_Instance == 0) m_Instance = new BkgNeutronWeight();
29 if (ke <= 0) return 0;
30
31 double x = log(ke);
32 int i = 0;
33 int j = m_n;
34 while (j > i + 1) {
35 int k = (i + j) / 2;
36 if (x < m_x[k]) j = k;
37 if (x >= m_x[k]) i = k;
38 }
39 double dx = x - m_x[i];
40 return (m_y[i] + dx * (m_b[i] + dx * (m_c[i] + dx * m_d[i])));
41 }
42
43 BkgNeutronWeight::BkgNeutronWeight(void)
44 {
45 addPoint(log(1.025E-10), 1.575E-02);
46 addPoint(log(1.075E-10), 1.537E-02);
47 addPoint(log(1.125E-10), 1.503E-02);
48 addPoint(log(1.175E-10), 1.470E-02);
49 addPoint(log(1.238E-10), 1.432E-02);
50 addPoint(log(1.313E-10), 1.391E-02);
51 addPoint(log(1.388E-10), 1.353E-02);
52 addPoint(log(1.463E-10), 1.317E-02);
53 addPoint(log(1.550E-10), 1.280E-02);
54 addPoint(log(1.650E-10), 1.242E-02);
55 addPoint(log(1.750E-10), 1.206E-02);
56 addPoint(log(1.850E-10), 1.172E-02);
57 addPoint(log(1.950E-10), 1.142E-02);
58 addPoint(log(2.050E-10), 1.113E-02);
59 addPoint(log(2.150E-10), 1.087E-02);
60 addPoint(log(2.250E-10), 1.063E-02);
61 addPoint(log(2.350E-10), 1.039E-02);
62 addPoint(log(2.475E-10), 1.013E-02);
63 addPoint(log(2.625E-10), 9.834E-03);
64 addPoint(log(2.750E-10), 9.609E-03);
65 addPoint(log(2.900E-10), 9.356E-03);
66 addPoint(log(3.100E-10), 9.048E-03);
67 addPoint(log(3.300E-10), 8.770E-03);
68 addPoint(log(3.500E-10), 8.516E-03);
69 addPoint(log(3.700E-10), 8.290E-03);
70 addPoint(log(3.900E-10), 8.075E-03);
71 addPoint(log(4.125E-10), 7.847E-03);
72 addPoint(log(4.375E-10), 7.622E-03);
73 addPoint(log(4.625E-10), 7.410E-03);
74 addPoint(log(4.875E-10), 7.219E-03);
75 addPoint(log(5.125E-10), 7.039E-03);
76 addPoint(log(5.375E-10), 6.875E-03);
77 addPoint(log(5.625E-10), 6.718E-03);
78 addPoint(log(5.875E-10), 6.575E-03);
79 addPoint(log(6.150E-10), 6.425E-03);
80 addPoint(log(6.450E-10), 6.274E-03);
81 addPoint(log(6.750E-10), 6.134E-03);
82 addPoint(log(7.050E-10), 6.002E-03);
83 addPoint(log(7.400E-10), 5.864E-03);
84 addPoint(log(7.800E-10), 5.709E-03);
85 addPoint(log(8.200E-10), 5.567E-03);
86 addPoint(log(8.600E-10), 5.438E-03);
87 addPoint(log(9.000E-10), 5.312E-03);
88 addPoint(log(9.400E-10), 5.200E-03);
89 addPoint(log(9.800E-10), 5.091E-03);
90 addPoint(log(1.025E-09), 4.978E-03);
91 addPoint(log(1.075E-09), 4.861E-03);
92 addPoint(log(1.125E-09), 4.751E-03);
93 addPoint(log(1.175E-09), 4.649E-03);
94 addPoint(log(1.238E-09), 4.530E-03);
95 addPoint(log(1.313E-09), 4.398E-03);
96 addPoint(log(1.388E-09), 4.278E-03);
97 addPoint(log(1.463E-09), 4.171E-03);
98 addPoint(log(1.550E-09), 4.049E-03);
99 addPoint(log(1.650E-09), 3.925E-03);
100 addPoint(log(1.750E-09), 3.810E-03);
101 addPoint(log(1.850E-09), 3.706E-03);
102 addPoint(log(1.950E-09), 3.609E-03);
103 addPoint(log(2.050E-09), 3.520E-03);
104 addPoint(log(2.150E-09), 3.437E-03);
105 addPoint(log(2.250E-09), 3.360E-03);
106 addPoint(log(2.350E-09), 3.287E-03);
107 addPoint(log(2.475E-09), 3.203E-03);
108 addPoint(log(2.625E-09), 3.110E-03);
109 addPoint(log(2.750E-09), 3.039E-03);
110 addPoint(log(2.900E-09), 2.962E-03);
111 addPoint(log(3.100E-09), 2.864E-03);
112 addPoint(log(3.300E-09), 2.776E-03);
113 addPoint(log(3.500E-09), 2.695E-03);
114 addPoint(log(3.700E-09), 2.621E-03);
115 addPoint(log(3.900E-09), 2.552E-03);
116 addPoint(log(4.125E-09), 2.482E-03);
117 addPoint(log(4.375E-09), 2.409E-03);
118 addPoint(log(4.625E-09), 2.344E-03);
119 addPoint(log(4.875E-09), 2.283E-03);
120 addPoint(log(5.125E-09), 2.226E-03);
121 addPoint(log(5.375E-09), 2.174E-03);
122 addPoint(log(5.625E-09), 2.125E-03);
123 addPoint(log(5.875E-09), 2.082E-03);
124 addPoint(log(6.150E-09), 2.034E-03);
125 addPoint(log(6.450E-09), 1.985E-03);
126 addPoint(log(6.750E-09), 1.941E-03);
127 addPoint(log(7.050E-09), 1.898E-03);
128 addPoint(log(7.400E-09), 1.853E-03);
129 addPoint(log(7.800E-09), 1.805E-03);
130 addPoint(log(8.200E-09), 1.760E-03);
131 addPoint(log(8.600E-09), 1.718E-03);
132 addPoint(log(9.000E-09), 1.680E-03);
133 addPoint(log(9.400E-09), 1.644E-03);
134 addPoint(log(9.800E-09), 1.610E-03);
135 addPoint(log(1.025E-08), 1.574E-03);
136 addPoint(log(1.075E-08), 1.537E-03);
137 addPoint(log(1.125E-08), 1.503E-03);
138 addPoint(log(1.175E-08), 1.472E-03);
139 addPoint(log(1.238E-08), 1.433E-03);
140 addPoint(log(1.313E-08), 1.392E-03);
141 addPoint(log(1.388E-08), 1.353E-03);
142 addPoint(log(1.463E-08), 1.318E-03);
143 addPoint(log(1.550E-08), 1.280E-03);
144 addPoint(log(1.650E-08), 1.241E-03);
145 addPoint(log(1.750E-08), 1.204E-03);
146 addPoint(log(1.850E-08), 1.172E-03);
147 addPoint(log(1.950E-08), 1.141E-03);
148 addPoint(log(2.050E-08), 1.113E-03);
149 addPoint(log(2.150E-08), 1.087E-03);
150 addPoint(log(2.250E-08), 1.062E-03);
151 addPoint(log(2.350E-08), 1.040E-03);
152 addPoint(log(2.475E-08), 1.013E-03);
153 addPoint(log(2.625E-08), 9.841E-04);
154 addPoint(log(2.750E-08), 9.619E-04);
155 addPoint(log(2.900E-08), 9.359E-04);
156 addPoint(log(3.100E-08), 9.058E-04);
157 addPoint(log(3.300E-08), 8.773E-04);
158 addPoint(log(3.500E-08), 8.523E-04);
159 addPoint(log(3.700E-08), 8.285E-04);
160 addPoint(log(3.900E-08), 8.073E-04);
161 addPoint(log(4.125E-08), 7.847E-04);
162 addPoint(log(4.375E-08), 7.619E-04);
163 addPoint(log(4.625E-08), 7.410E-04);
164 addPoint(log(4.875E-08), 7.215E-04);
165 addPoint(log(5.125E-08), 7.036E-04);
166 addPoint(log(5.375E-08), 6.871E-04);
167 addPoint(log(5.625E-08), 6.716E-04);
168 addPoint(log(5.875E-08), 6.578E-04);
169 addPoint(log(6.150E-08), 6.428E-04);
170 addPoint(log(6.450E-08), 6.273E-04);
171 addPoint(log(6.750E-08), 6.136E-04);
172 addPoint(log(7.050E-08), 6.001E-04);
173 addPoint(log(7.400E-08), 5.861E-04);
174 addPoint(log(7.800E-08), 5.707E-04);
175 addPoint(log(8.200E-08), 5.567E-04);
176 addPoint(log(8.600E-08), 5.434E-04);
177 addPoint(log(9.000E-08), 5.313E-04);
178 addPoint(log(9.400E-08), 5.197E-04);
179 addPoint(log(9.800E-08), 5.090E-04);
180 addPoint(log(1.025E-07), 4.976E-04);
181 addPoint(log(1.075E-07), 4.864E-04);
182 addPoint(log(1.125E-07), 4.753E-04);
183 addPoint(log(1.175E-07), 4.646E-04);
184 addPoint(log(1.238E-07), 4.531E-04);
185 addPoint(log(1.313E-07), 4.397E-04);
186 addPoint(log(1.388E-07), 4.279E-04);
187 addPoint(log(1.463E-07), 4.166E-04);
188 addPoint(log(1.550E-07), 4.049E-04);
189 addPoint(log(1.650E-07), 3.924E-04);
190 addPoint(log(1.750E-07), 3.810E-04);
191 addPoint(log(1.850E-07), 3.705E-04);
192 addPoint(log(1.950E-07), 3.607E-04);
193 addPoint(log(2.050E-07), 3.518E-04);
194 addPoint(log(2.150E-07), 3.437E-04);
195 addPoint(log(2.250E-07), 3.362E-04);
196 addPoint(log(2.350E-07), 3.286E-04);
197 addPoint(log(2.475E-07), 3.204E-04);
198 addPoint(log(2.625E-07), 3.111E-04);
199 addPoint(log(2.750E-07), 3.039E-04);
200 addPoint(log(2.900E-07), 2.962E-04);
201 addPoint(log(3.100E-07), 2.865E-04);
202 addPoint(log(3.300E-07), 2.775E-04);
203 addPoint(log(3.500E-07), 2.695E-04);
204 addPoint(log(3.700E-07), 2.620E-04);
205 addPoint(log(3.900E-07), 2.551E-04);
206 addPoint(log(4.125E-07), 2.483E-04);
207 addPoint(log(4.375E-07), 2.410E-04);
208 addPoint(log(4.625E-07), 2.343E-04);
209 addPoint(log(4.875E-07), 2.284E-04);
210 addPoint(log(5.125E-07), 2.226E-04);
211 addPoint(log(5.375E-07), 2.176E-04);
212 addPoint(log(5.625E-07), 2.128E-04);
213 addPoint(log(5.875E-07), 2.082E-04);
214 addPoint(log(6.150E-07), 2.034E-04);
215 addPoint(log(6.450E-07), 1.985E-04);
216 addPoint(log(6.750E-07), 1.941E-04);
217 addPoint(log(7.050E-07), 1.898E-04);
218 addPoint(log(7.400E-07), 1.852E-04);
219 addPoint(log(7.800E-07), 1.805E-04);
220 addPoint(log(8.200E-07), 1.762E-04);
221 addPoint(log(8.600E-07), 1.718E-04);
222 addPoint(log(9.000E-07), 1.679E-04);
223 addPoint(log(9.400E-07), 1.645E-04);
224 addPoint(log(9.800E-07), 1.610E-04);
225 addPoint(log(1.025E-06), 1.575E-04);
226 addPoint(log(1.075E-06), 1.542E-04);
227 addPoint(log(1.125E-06), 1.508E-04);
228 addPoint(log(1.175E-06), 1.474E-04);
229 addPoint(log(1.238E-06), 1.434E-04);
230 addPoint(log(1.313E-06), 1.393E-04);
231 addPoint(log(1.388E-06), 1.356E-04);
232 addPoint(log(1.463E-06), 1.318E-04);
233 addPoint(log(1.550E-06), 1.281E-04);
234 addPoint(log(1.650E-06), 1.242E-04);
235 addPoint(log(1.750E-06), 1.204E-04);
236 addPoint(log(1.850E-06), 1.172E-04);
237 addPoint(log(1.950E-06), 1.141E-04);
238 addPoint(log(2.050E-06), 1.114E-04);
239 addPoint(log(2.150E-06), 1.090E-04);
240 addPoint(log(2.250E-06), 1.066E-04);
241 addPoint(log(2.350E-06), 1.042E-04);
242 addPoint(log(2.475E-06), 1.013E-04);
243 addPoint(log(2.625E-06), 9.853E-05);
244 addPoint(log(2.750E-06), 9.627E-05);
245 addPoint(log(2.900E-06), 9.360E-05);
246 addPoint(log(3.100E-06), 9.056E-05);
247 addPoint(log(3.300E-06), 8.776E-05);
248 addPoint(log(3.500E-06), 8.513E-05);
249 addPoint(log(3.700E-06), 8.287E-05);
250 addPoint(log(3.900E-06), 8.065E-05);
251 addPoint(log(4.125E-06), 7.860E-05);
252 addPoint(log(4.375E-06), 7.644E-05);
253 addPoint(log(4.625E-06), 7.428E-05);
254 addPoint(log(4.875E-06), 7.217E-05);
255 addPoint(log(5.125E-06), 7.045E-05);
256 addPoint(log(5.375E-06), 6.885E-05);
257 addPoint(log(5.625E-06), 6.725E-05);
258 addPoint(log(5.875E-06), 6.569E-05);
259 addPoint(log(6.150E-06), 6.427E-05);
260 addPoint(log(6.450E-06), 6.278E-05);
261 addPoint(log(6.750E-06), 6.129E-05);
262 addPoint(log(7.050E-06), 5.998E-05);
263 addPoint(log(7.400E-06), 5.858E-05);
264 addPoint(log(7.800E-06), 5.701E-05);
265 addPoint(log(8.200E-06), 5.572E-05);
266 addPoint(log(8.600E-06), 5.450E-05);
267 addPoint(log(9.000E-06), 5.327E-05);
268 addPoint(log(9.400E-06), 5.205E-05);
269 addPoint(log(9.800E-06), 5.087E-05);
270 addPoint(log(1.025E-05), 4.981E-05);
271 addPoint(log(1.075E-05), 4.868E-05);
272 addPoint(log(1.125E-05), 4.754E-05);
273 addPoint(log(1.175E-05), 4.644E-05);
274 addPoint(log(1.238E-05), 4.531E-05);
275 addPoint(log(1.313E-05), 4.399E-05);
276 addPoint(log(1.388E-05), 4.274E-05);
277 addPoint(log(1.463E-05), 4.167E-05);
278 addPoint(log(1.550E-05), 4.045E-05);
279 addPoint(log(1.650E-05), 3.929E-05);
280 addPoint(log(1.750E-05), 3.821E-05);
281 addPoint(log(1.850E-05), 3.712E-05);
282 addPoint(log(1.950E-05), 3.607E-05);
283 addPoint(log(2.050E-05), 3.522E-05);
284 addPoint(log(2.150E-05), 3.442E-05);
285 addPoint(log(2.250E-05), 3.361E-05);
286 addPoint(log(2.350E-05), 3.284E-05);
287 addPoint(log(2.475E-05), 3.204E-05);
288 addPoint(log(2.625E-05), 3.110E-05);
289 addPoint(log(2.750E-05), 3.035E-05);
290 addPoint(log(2.900E-05), 2.959E-05);
291 addPoint(log(3.100E-05), 2.860E-05);
292 addPoint(log(3.300E-05), 2.778E-05);
293 addPoint(log(3.500E-05), 2.702E-05);
294 addPoint(log(3.700E-05), 2.625E-05);
295 addPoint(log(3.900E-05), 2.550E-05);
296 addPoint(log(4.125E-05), 2.483E-05);
297 addPoint(log(4.375E-05), 2.412E-05);
298 addPoint(log(4.625E-05), 2.342E-05);
299 addPoint(log(4.875E-05), 2.282E-05);
300 addPoint(log(5.125E-05), 2.227E-05);
301 addPoint(log(5.375E-05), 2.172E-05);
302 addPoint(log(5.625E-05), 2.123E-05);
303 addPoint(log(5.875E-05), 2.079E-05);
304 addPoint(log(6.150E-05), 2.030E-05);
305 addPoint(log(6.450E-05), 1.985E-05);
306 addPoint(log(6.750E-05), 1.944E-05);
307 addPoint(log(7.050E-05), 1.904E-05);
308 addPoint(log(7.400E-05), 1.856E-05);
309 addPoint(log(7.800E-05), 1.803E-05);
310 addPoint(log(8.200E-05), 1.761E-05);
311 addPoint(log(8.600E-05), 1.721E-05);
312 addPoint(log(9.000E-05), 1.681E-05);
313 addPoint(log(9.400E-05), 1.642E-05);
314 addPoint(log(9.800E-05), 1.610E-05);
315 addPoint(log(1.025E-04), 1.575E-05);
316 addPoint(log(1.075E-04), 1.536E-05);
317 addPoint(log(1.125E-04), 1.502E-05);
318 addPoint(log(1.175E-04), 1.470E-05);
319 addPoint(log(1.238E-04), 1.432E-05);
320 addPoint(log(1.313E-04), 1.393E-05);
321 addPoint(log(1.388E-04), 1.357E-05);
322 addPoint(log(1.463E-04), 1.321E-05);
323 addPoint(log(1.550E-04), 1.714E-05);
324 addPoint(log(1.650E-04), 6.792E-05);
325 addPoint(log(1.750E-04), 1.315E-04);
326 addPoint(log(1.850E-04), 1.938E-04);
327 addPoint(log(1.950E-04), 2.250E-04);
328 addPoint(log(2.050E-04), 2.433E-04);
329 addPoint(log(2.150E-04), 2.614E-04);
330 addPoint(log(2.250E-04), 2.754E-04);
331 addPoint(log(2.350E-04), 2.867E-04);
332 addPoint(log(2.475E-04), 3.007E-04);
333 addPoint(log(2.625E-04), 3.177E-04);
334 addPoint(log(2.750E-04), 3.318E-04);
335 addPoint(log(2.900E-04), 3.486E-04);
336 addPoint(log(3.100E-04), 3.712E-04);
337 addPoint(log(3.300E-04), 3.937E-04);
338 addPoint(log(3.500E-04), 4.163E-04);
339 addPoint(log(3.700E-04), 4.388E-04);
340 addPoint(log(3.900E-04), 4.613E-04);
341 addPoint(log(4.125E-04), 4.865E-04);
342 addPoint(log(4.375E-04), 5.146E-04);
343 addPoint(log(4.625E-04), 5.426E-04);
344 addPoint(log(4.875E-04), 5.706E-04);
345 addPoint(log(5.125E-04), 5.985E-04);
346 addPoint(log(5.375E-04), 6.264E-04);
347 addPoint(log(5.625E-04), 6.543E-04);
348 addPoint(log(5.875E-04), 6.822E-04);
349 addPoint(log(6.150E-04), 7.128E-04);
350 addPoint(log(6.450E-04), 7.462E-04);
351 addPoint(log(6.750E-04), 7.795E-04);
352 addPoint(log(7.050E-04), 8.128E-04);
353 addPoint(log(7.400E-04), 8.515E-04);
354 addPoint(log(7.800E-04), 8.957E-04);
355 addPoint(log(8.200E-04), 9.399E-04);
356 addPoint(log(8.600E-04), 9.840E-04);
357 addPoint(log(9.000E-04), 1.028E-03);
358 addPoint(log(9.400E-04), 1.072E-03);
359 addPoint(log(9.800E-04), 1.116E-03);
360 addPoint(log(1.025E-03), 1.165E-03);
361 addPoint(log(1.075E-03), 1.220E-03);
362 addPoint(log(1.125E-03), 1.274E-03);
363 addPoint(log(1.175E-03), 1.329E-03);
364 addPoint(log(1.238E-03), 1.397E-03);
365 addPoint(log(1.313E-03), 1.478E-03);
366 addPoint(log(1.388E-03), 1.559E-03);
367 addPoint(log(1.463E-03), 1.640E-03);
368 addPoint(log(1.550E-03), 1.734E-03);
369 addPoint(log(1.650E-03), 1.842E-03);
370 addPoint(log(1.750E-03), 1.949E-03);
371 addPoint(log(1.850E-03), 2.056E-03);
372 addPoint(log(1.950E-03), 2.162E-03);
373 addPoint(log(2.050E-03), 2.269E-03);
374 addPoint(log(2.150E-03), 2.382E-03);
375 addPoint(log(2.250E-03), 4.405E-03);
376 addPoint(log(2.350E-03), 2.589E-03);
377 addPoint(log(2.475E-03), 2.717E-03);
378 addPoint(log(2.625E-03), 2.874E-03);
379 addPoint(log(2.750E-03), 3.004E-03);
380 addPoint(log(2.900E-03), 3.159E-03);
381 addPoint(log(3.100E-03), 3.366E-03);
382 addPoint(log(3.300E-03), 3.573E-03);
383 addPoint(log(3.500E-03), 3.778E-03);
384 addPoint(log(3.700E-03), 3.982E-03);
385 addPoint(log(3.900E-03), 4.186E-03);
386 addPoint(log(4.125E-03), 4.413E-03);
387 addPoint(log(4.375E-03), 4.666E-03);
388 addPoint(log(4.625E-03), 4.917E-03);
389 addPoint(log(4.875E-03), 5.492E-03);
390 addPoint(log(5.125E-03), 5.419E-03);
391 addPoint(log(5.375E-03), 5.662E-03);
392 addPoint(log(5.625E-03), 5.907E-03);
393 addPoint(log(5.875E-03), 6.151E-03);
394 addPoint(log(6.150E-03), 6.418E-03);
395 addPoint(log(6.450E-03), 6.710E-03);
396 addPoint(log(6.750E-03), 6.999E-03);
397 addPoint(log(7.050E-03), 7.285E-03);
398 addPoint(log(7.400E-03), 7.618E-03);
399 addPoint(log(7.800E-03), 7.998E-03);
400 addPoint(log(8.200E-03), 8.372E-03);
401 addPoint(log(8.600E-03), 8.745E-03);
402 addPoint(log(9.000E-03), 9.115E-03);
403 addPoint(log(9.400E-03), 9.477E-03);
404 addPoint(log(9.800E-03), 9.838E-03);
405 addPoint(log(1.025E-02), 1.024E-02);
406 addPoint(log(1.075E-02), 1.069E-02);
407 addPoint(log(1.125E-02), 1.115E-02);
408 addPoint(log(1.175E-02), 1.159E-02);
409 addPoint(log(1.238E-02), 1.213E-02);
410 addPoint(log(1.313E-02), 1.277E-02);
411 addPoint(log(1.388E-02), 1.340E-02);
412 addPoint(log(1.463E-02), 1.403E-02);
413 addPoint(log(1.550E-02), 1.557E-02);
414 addPoint(log(1.650E-02), 1.554E-02);
415 addPoint(log(1.750E-02), 1.634E-02);
416 addPoint(log(1.850E-02), 1.710E-02);
417 addPoint(log(1.950E-02), 1.782E-02);
418 addPoint(log(2.050E-02), 1.854E-02);
419 addPoint(log(2.150E-02), 1.926E-02);
420 addPoint(log(2.250E-02), 1.998E-02);
421 addPoint(log(2.350E-02), 2.069E-02);
422 addPoint(log(2.475E-02), 2.150E-02);
423 addPoint(log(2.625E-02), 2.245E-02);
424 addPoint(log(2.750E-02), 2.317E-02);
425 addPoint(log(2.900E-02), 2.400E-02);
426 addPoint(log(3.100E-02), 2.509E-02);
427 addPoint(log(3.300E-02), 2.599E-02);
428 addPoint(log(3.500E-02), 2.675E-02);
429 addPoint(log(3.700E-02), 2.743E-02);
430 addPoint(log(3.900E-02), 3.247E-02);
431 addPoint(log(4.125E-02), 2.792E-02);
432 addPoint(log(4.375E-02), 2.739E-02);
433 addPoint(log(4.625E-02), 2.598E-02);
434 addPoint(log(4.875E-02), 2.288E-02);
435 addPoint(log(5.125E-02), 1.626E-02);
436 addPoint(log(5.375E-02), 1.485E-02);
437 addPoint(log(5.625E-02), 5.078E-01);
438 addPoint(log(5.875E-02), 1.145E-01);
439 addPoint(log(6.150E-02), 7.793E-02);
440 addPoint(log(6.450E-02), 6.706E-02);
441 addPoint(log(6.750E-02), 6.455E-02);
442 addPoint(log(7.050E-02), 5.988E-02);
443 addPoint(log(7.400E-02), 5.797E-02);
444 addPoint(log(7.800E-02), 5.636E-02);
445 addPoint(log(8.200E-02), 5.494E-02);
446 addPoint(log(8.600E-02), 5.374E-02);
447 addPoint(log(9.000E-02), 5.188E-02);
448 addPoint(log(9.400E-02), 4.992E-02);
449 addPoint(log(9.800E-02), 4.790E-02);
450 addPoint(log(1.025E-01), 4.526E-02);
451 addPoint(log(1.075E-01), 4.162E-02);
452 addPoint(log(1.125E-01), 3.764E-02);
453 addPoint(log(1.175E-01), 3.297E-02);
454 addPoint(log(1.238E-01), 2.645E-02);
455 addPoint(log(1.313E-01), 1.826E-02);
456 addPoint(log(1.388E-01), 1.121E-02);
457 addPoint(log(1.463E-01), 1.090E-02);
458 addPoint(log(1.550E-01), 4.848E-02);
459 addPoint(log(1.650E-01), 2.173E-01);
460 addPoint(log(1.750E-01), 6.987E-01);
461 addPoint(log(1.850E-01), 1.196E+00);
462 addPoint(log(1.950E-01), 1.162E+00);
463 addPoint(log(2.050E-01), 9.710E-01);
464 addPoint(log(2.150E-01), 8.288E-01);
465 addPoint(log(2.250E-01), 7.334E-01);
466 addPoint(log(2.350E-01), 6.709E-01);
467 addPoint(log(2.475E-01), 6.172E-01);
468 addPoint(log(2.625E-01), 5.768E-01);
469 addPoint(log(2.750E-01), 5.547E-01);
470 addPoint(log(2.900E-01), 5.393E-01);
471 addPoint(log(3.100E-01), 5.237E-01);
472 addPoint(log(3.300E-01), 5.198E-01);
473 addPoint(log(3.500E-01), 5.097E-01);
474 addPoint(log(3.700E-01), 5.036E-01);
475 addPoint(log(3.900E-01), 5.320E-01);
476 addPoint(log(4.125E-01), 5.403E-01);
477 addPoint(log(4.375E-01), 5.451E-01);
478 addPoint(log(4.625E-01), 5.581E-01);
479 addPoint(log(4.875E-01), 5.740E-01);
480 addPoint(log(5.125E-01), 5.987E-01);
481 addPoint(log(5.375E-01), 7.554E-01);
482 addPoint(log(5.625E-01), 1.280E+00);
483 addPoint(log(5.875E-01), 5.997E-01);
484 addPoint(log(6.150E-01), 5.425E-01);
485 addPoint(log(6.450E-01), 5.536E-01);
486 addPoint(log(6.750E-01), 5.786E-01);
487 addPoint(log(7.050E-01), 6.019E-01);
488 addPoint(log(7.400E-01), 6.724E-01);
489 addPoint(log(7.800E-01), 9.209E-01);
490 addPoint(log(8.200E-01), 1.459E+00);
491 addPoint(log(8.600E-01), 8.318E-01);
492 addPoint(log(9.000E-01), 9.434E-01);
493 addPoint(log(9.400E-01), 1.172E+00);
494 addPoint(log(9.800E-01), 1.174E+00);
495 addPoint(log(1.050E+00), 8.020E-01);
496 addPoint(log(1.150E+00), 6.578E-01);
497 addPoint(log(1.250E+00), 9.680E-01);
498 addPoint(log(1.350E+00), 9.410E-01);
499 addPoint(log(1.450E+00), 1.079E+00);
500 addPoint(log(1.550E+00), 1.128E+00);
501 addPoint(log(1.650E+00), 1.766E+00);
502 addPoint(log(1.750E+00), 8.366E-01);
503 addPoint(log(1.850E+00), 1.411E+00);
504 addPoint(log(1.950E+00), 1.393E+00);
505 addPoint(log(2.050E+00), 1.022E+00);
506 addPoint(log(2.150E+00), 1.159E+00);
507 addPoint(log(2.250E+00), 1.126E+00);
508 addPoint(log(2.350E+00), 1.106E+00);
509 addPoint(log(2.450E+00), 1.267E+00);
510 addPoint(log(2.550E+00), 1.384E+00);
511 addPoint(log(2.650E+00), 1.238E+00);
512 addPoint(log(2.750E+00), 1.153E+00);
513 addPoint(log(2.850E+00), 1.437E+00);
514 addPoint(log(2.950E+00), 1.061E+00);
515 addPoint(log(3.050E+00), 1.278E+00);
516 addPoint(log(3.150E+00), 1.427E+00);
517 addPoint(log(3.250E+00), 1.281E+00);
518 addPoint(log(3.350E+00), 1.217E+00);
519 addPoint(log(3.450E+00), 1.259E+00);
520 addPoint(log(3.550E+00), 1.213E+00);
521 addPoint(log(3.650E+00), 7.506E-01);
522 addPoint(log(3.750E+00), 1.235E+00);
523 addPoint(log(3.850E+00), 1.187E+00);
524 addPoint(log(3.950E+00), 1.452E+00);
525 addPoint(log(4.050E+00), 1.431E+00);
526 addPoint(log(4.150E+00), 1.147E+00);
527 addPoint(log(4.250E+00), 1.782E+00);
528 addPoint(log(4.350E+00), 1.442E+00);
529 addPoint(log(4.450E+00), 1.509E+00);
530 addPoint(log(4.550E+00), 1.498E+00);
531 addPoint(log(4.650E+00), 1.673E+00);
532 addPoint(log(4.750E+00), 2.003E+00);
533 addPoint(log(4.850E+00), 1.717E+00);
534 addPoint(log(4.950E+00), 1.579E+00);
535 addPoint(log(5.050E+00), 1.612E+00);
536 addPoint(log(5.150E+00), 1.831E+00);
537 addPoint(log(5.250E+00), 1.605E+00);
538 addPoint(log(5.350E+00), 1.336E+00);
539 addPoint(log(5.450E+00), 1.277E+00);
540 addPoint(log(5.550E+00), 1.567E+00);
541 addPoint(log(5.650E+00), 1.616E+00);
542 addPoint(log(5.750E+00), 1.929E+00);
543 addPoint(log(5.850E+00), 1.794E+00);
544 addPoint(log(5.950E+00), 1.468E+00);
545 addPoint(log(6.050E+00), 1.666E+00);
546 addPoint(log(6.150E+00), 1.364E+00);
547 addPoint(log(6.250E+00), 1.900E+00);
548 addPoint(log(6.350E+00), 1.641E+00);
549 addPoint(log(6.450E+00), 1.514E+00);
550 addPoint(log(6.550E+00), 1.314E+00);
551 addPoint(log(6.650E+00), 1.611E+00);
552 addPoint(log(6.750E+00), 1.778E+00);
553 addPoint(log(6.850E+00), 1.577E+00);
554 addPoint(log(6.950E+00), 1.521E+00);
555 addPoint(log(7.050E+00), 1.770E+00);
556 addPoint(log(7.150E+00), 1.470E+00);
557 addPoint(log(7.250E+00), 1.803E+00);
558 addPoint(log(7.350E+00), 1.818E+00);
559 addPoint(log(7.450E+00), 1.793E+00);
560 addPoint(log(7.550E+00), 1.757E+00);
561 addPoint(log(7.650E+00), 1.794E+00);
562 addPoint(log(7.750E+00), 1.866E+00);
563 addPoint(log(7.850E+00), 1.861E+00);
564 addPoint(log(7.950E+00), 1.992E+00);
565 addPoint(log(8.050E+00), 1.806E+00);
566 addPoint(log(8.150E+00), 1.663E+00);
567 addPoint(log(8.250E+00), 1.743E+00);
568 addPoint(log(8.350E+00), 1.768E+00);
569 addPoint(log(8.450E+00), 1.782E+00);
570 addPoint(log(8.550E+00), 1.771E+00);
571 addPoint(log(8.650E+00), 1.769E+00);
572 addPoint(log(8.750E+00), 1.589E+00);
573 addPoint(log(8.850E+00), 1.685E+00);
574 addPoint(log(8.950E+00), 1.871E+00);
575 addPoint(log(9.050E+00), 1.876E+00);
576 addPoint(log(9.150E+00), 1.686E+00);
577 addPoint(log(9.250E+00), 1.619E+00);
578 addPoint(log(9.350E+00), 1.751E+00);
579 addPoint(log(9.450E+00), 1.824E+00);
580 addPoint(log(9.550E+00), 1.767E+00);
581 addPoint(log(9.650E+00), 1.669E+00);
582 addPoint(log(9.750E+00), 1.707E+00);
583 addPoint(log(9.850E+00), 1.754E+00);
584 addPoint(log(9.950E+00), 1.780E+00);
585 addPoint(log(1.005E+01), 1.760E+00);
586 addPoint(log(1.015E+01), 1.782E+00);
587 addPoint(log(1.025E+01), 1.764E+00);
588 addPoint(log(1.035E+01), 1.772E+00);
589 addPoint(log(1.045E+01), 1.709E+00);
590 addPoint(log(1.055E+01), 1.696E+00);
591 addPoint(log(1.065E+01), 1.674E+00);
592 addPoint(log(1.075E+01), 1.726E+00);
593 addPoint(log(1.085E+01), 1.765E+00);
594 addPoint(log(1.095E+01), 1.765E+00);
595 addPoint(log(1.105E+01), 1.736E+00);
596 addPoint(log(1.115E+01), 1.680E+00);
597 addPoint(log(1.125E+01), 1.692E+00);
598 addPoint(log(1.135E+01), 1.735E+00);
599 addPoint(log(1.145E+01), 1.765E+00);
600 addPoint(log(1.155E+01), 1.778E+00);
601 addPoint(log(1.165E+01), 1.767E+00);
602 addPoint(log(1.175E+01), 1.751E+00);
603 addPoint(log(1.185E+01), 1.754E+00);
604 addPoint(log(1.195E+01), 1.767E+00);
605 addPoint(log(1.205E+01), 1.784E+00);
606 addPoint(log(1.215E+01), 1.812E+00);
607 addPoint(log(1.225E+01), 1.816E+00);
608 addPoint(log(1.235E+01), 1.822E+00);
609 addPoint(log(1.245E+01), 1.791E+00);
610 addPoint(log(1.255E+01), 1.757E+00);
611 addPoint(log(1.265E+01), 1.837E+00);
612 addPoint(log(1.275E+01), 1.856E+00);
613 addPoint(log(1.285E+01), 1.850E+00);
614 addPoint(log(1.295E+01), 1.853E+00);
615 addPoint(log(1.305E+01), 1.862E+00);
616 addPoint(log(1.315E+01), 1.828E+00);
617 addPoint(log(1.325E+01), 1.793E+00);
618 addPoint(log(1.335E+01), 1.806E+00);
619 addPoint(log(1.345E+01), 1.839E+00);
620 addPoint(log(1.355E+01), 1.802E+00);
621 addPoint(log(1.365E+01), 1.780E+00);
622 addPoint(log(1.375E+01), 1.804E+00);
623 addPoint(log(1.385E+01), 1.870E+00);
624 addPoint(log(1.395E+01), 1.808E+00);
625 addPoint(log(1.405E+01), 1.787E+00);
626 addPoint(log(1.415E+01), 1.817E+00);
627 addPoint(log(1.425E+01), 1.830E+00);
628 addPoint(log(1.435E+01), 1.843E+00);
629 addPoint(log(1.445E+01), 1.845E+00);
630 addPoint(log(1.455E+01), 1.791E+00);
631 addPoint(log(1.465E+01), 1.765E+00);
632 addPoint(log(1.475E+01), 1.774E+00);
633 addPoint(log(1.485E+01), 1.808E+00);
634 addPoint(log(1.495E+01), 1.813E+00);
635 addPoint(log(1.505E+01), 1.773E+00);
636 addPoint(log(1.515E+01), 1.801E+00);
637 addPoint(log(1.525E+01), 1.842E+00);
638 addPoint(log(1.535E+01), 1.841E+00);
639 addPoint(log(1.545E+01), 1.828E+00);
640 addPoint(log(1.555E+01), 1.808E+00);
641 addPoint(log(1.565E+01), 1.806E+00);
642 addPoint(log(1.575E+01), 1.817E+00);
643 addPoint(log(1.585E+01), 1.858E+00);
644 addPoint(log(1.595E+01), 1.891E+00);
645 addPoint(log(1.605E+01), 1.899E+00);
646 addPoint(log(1.615E+01), 1.899E+00);
647 addPoint(log(1.625E+01), 1.893E+00);
648 addPoint(log(1.635E+01), 1.880E+00);
649 addPoint(log(1.645E+01), 1.869E+00);
650 addPoint(log(1.655E+01), 1.918E+00);
651 addPoint(log(1.665E+01), 1.970E+00);
652 addPoint(log(1.675E+01), 1.952E+00);
653 addPoint(log(1.685E+01), 1.929E+00);
654 addPoint(log(1.695E+01), 1.935E+00);
655 addPoint(log(1.705E+01), 1.940E+00);
656 addPoint(log(1.715E+01), 1.937E+00);
657 addPoint(log(1.725E+01), 1.934E+00);
658 addPoint(log(1.735E+01), 1.935E+00);
659 addPoint(log(1.745E+01), 1.932E+00);
660 addPoint(log(1.755E+01), 1.917E+00);
661 addPoint(log(1.765E+01), 1.906E+00);
662 addPoint(log(1.775E+01), 1.914E+00);
663 addPoint(log(1.785E+01), 1.926E+00);
664 addPoint(log(1.795E+01), 1.951E+00);
665 addPoint(log(1.805E+01), 1.974E+00);
666 addPoint(log(1.815E+01), 1.965E+00);
667 addPoint(log(1.825E+01), 1.950E+00);
668 addPoint(log(1.835E+01), 1.938E+00);
669 addPoint(log(1.845E+01), 1.926E+00);
670 addPoint(log(1.855E+01), 1.945E+00);
671 addPoint(log(1.865E+01), 1.986E+00);
672 addPoint(log(1.875E+01), 2.008E+00);
673 addPoint(log(1.885E+01), 1.999E+00);
674 addPoint(log(1.895E+01), 1.989E+00);
675 addPoint(log(1.905E+01), 1.974E+00);
676 addPoint(log(1.915E+01), 1.960E+00);
677 addPoint(log(1.925E+01), 1.955E+00);
678 addPoint(log(1.935E+01), 1.953E+00);
679 addPoint(log(1.945E+01), 1.953E+00);
680 addPoint(log(1.955E+01), 1.955E+00);
681 addPoint(log(1.965E+01), 1.958E+00);
682 addPoint(log(1.975E+01), 1.963E+00);
683 addPoint(log(1.985E+01), 1.969E+00);
684 addPoint(log(1.995E+01), 1.975E+00);
685 addPoint(log(2.000E+01), 2.071E+00);
686 addPoint(log(2.500E+01), 2.049E+00);
687 addPoint(log(3.000E+01), 2.041E+00);
688 addPoint(log(4.000E+01), 2.012E+00);
689 addPoint(log(5.000E+01), 1.805E+00);
690 addPoint(log(6.000E+01), 1.644E+00);
691 addPoint(log(7.000E+01), 1.499E+00);
692 addPoint(log(8.000E+01), 1.378E+00);
693 addPoint(log(9.000E+01), 1.264E+00);
694 addPoint(log(1.000E+02), 1.168E+00);
695 addPoint(log(1.300E+02), 9.740E-01);
696 addPoint(log(1.600E+02), 8.650E-01);
697 addPoint(log(2.000E+02), 7.910E-01);
698 addPoint(log(2.500E+02), 7.330E-01);
699 addPoint(log(3.000E+02), 6.960E-01);
700 addPoint(log(3.500E+02), 6.930E-01);
701 addPoint(log(4.000E+02), 6.850E-01);
702 addPoint(log(4.500E+02), 6.850E-01);
703 addPoint(log(5.000E+02), 6.810E-01);
704 addPoint(log(6.000E+02), 6.790E-01);
705 addPoint(log(7.000E+02), 6.770E-01);
706 addPoint(log(8.000E+02), 6.720E-01);
707 addPoint(log(8.050E+02), 6.004E-01);
708 addPoint(log(8.150E+02), 5.980E-01);
709 addPoint(log(8.250E+02), 5.959E-01);
710 addPoint(log(8.350E+02), 5.942E-01);
711 addPoint(log(8.450E+02), 5.932E-01);
712 addPoint(log(8.550E+02), 5.922E-01);
713 addPoint(log(8.650E+02), 5.912E-01);
714 addPoint(log(8.750E+02), 5.902E-01);
715 addPoint(log(8.850E+02), 5.892E-01);
716 addPoint(log(8.950E+02), 5.883E-01);
717 addPoint(log(9.050E+02), 5.873E-01);
718 addPoint(log(9.150E+02), 5.863E-01);
719 addPoint(log(9.250E+02), 5.854E-01);
720 addPoint(log(9.350E+02), 5.845E-01);
721 addPoint(log(9.450E+02), 5.836E-01);
722 addPoint(log(9.550E+02), 5.828E-01);
723 addPoint(log(9.650E+02), 5.819E-01);
724 addPoint(log(9.750E+02), 5.810E-01);
725 addPoint(log(9.850E+02), 5.801E-01);
726 addPoint(log(9.950E+02), 5.793E-01);
727 addPoint(log(1.005E+03), 5.784E-01);
728 addPoint(log(1.015E+03), 5.773E-01);
729 addPoint(log(1.025E+03), 5.763E-01);
730 addPoint(log(1.035E+03), 5.752E-01);
731 addPoint(log(1.045E+03), 5.741E-01);
732 addPoint(log(1.055E+03), 5.734E-01);
733 addPoint(log(1.065E+03), 5.728E-01);
734 addPoint(log(1.075E+03), 5.722E-01);
735 addPoint(log(1.085E+03), 5.716E-01);
736 addPoint(log(1.095E+03), 5.711E-01);
737 addPoint(log(1.105E+03), 5.708E-01);
738 addPoint(log(1.115E+03), 5.704E-01);
739 addPoint(log(1.125E+03), 5.701E-01);
740 addPoint(log(1.135E+03), 5.698E-01);
741 addPoint(log(1.145E+03), 5.694E-01);
742 addPoint(log(1.155E+03), 5.691E-01);
743 addPoint(log(1.165E+03), 5.688E-01);
744 addPoint(log(1.175E+03), 5.684E-01);
745 addPoint(log(1.185E+03), 5.681E-01);
746 addPoint(log(1.195E+03), 5.678E-01);
747 addPoint(log(1.205E+03), 5.674E-01);
748 addPoint(log(1.215E+03), 5.671E-01);
749 addPoint(log(1.225E+03), 5.668E-01);
750 addPoint(log(1.235E+03), 5.664E-01);
751 addPoint(log(1.245E+03), 5.661E-01);
752 addPoint(log(1.255E+03), 5.658E-01);
753 addPoint(log(1.265E+03), 5.655E-01);
754 addPoint(log(1.275E+03), 5.650E-01);
755 addPoint(log(1.285E+03), 5.643E-01);
756 addPoint(log(1.295E+03), 5.636E-01);
757 addPoint(log(1.305E+03), 5.628E-01);
758 addPoint(log(1.315E+03), 5.621E-01);
759 addPoint(log(1.325E+03), 5.614E-01);
760 addPoint(log(1.335E+03), 5.607E-01);
761 addPoint(log(1.345E+03), 5.600E-01);
762 addPoint(log(1.355E+03), 5.593E-01);
763 addPoint(log(1.365E+03), 5.587E-01);
764 addPoint(log(1.375E+03), 5.581E-01);
765 addPoint(log(1.385E+03), 5.575E-01);
766 addPoint(log(1.395E+03), 5.569E-01);
767 addPoint(log(1.405E+03), 5.563E-01);
768 addPoint(log(1.415E+03), 5.557E-01);
769 addPoint(log(1.425E+03), 5.551E-01);
770 addPoint(log(1.435E+03), 5.545E-01);
771 addPoint(log(1.445E+03), 5.539E-01);
772 addPoint(log(1.455E+03), 5.533E-01);
773 addPoint(log(1.465E+03), 5.527E-01);
774 addPoint(log(1.485E+03), 5.515E-01);
775 addPoint(log(1.495E+03), 5.510E-01);
776 addPoint(log(1.505E+03), 5.504E-01);
777 addPoint(log(1.515E+03), 5.498E-01);
778 addPoint(log(1.535E+03), 5.486E-01);
779 addPoint(log(1.545E+03), 5.480E-01);
780 addPoint(log(1.555E+03), 5.474E-01);
781 addPoint(log(1.565E+03), 5.468E-01);
782 addPoint(log(1.585E+03), 5.457E-01);
783 addPoint(log(1.595E+03), 5.451E-01);
784 addPoint(log(1.605E+03), 5.445E-01);
785 addPoint(log(1.615E+03), 5.439E-01);
786 addPoint(log(1.635E+03), 5.427E-01);
787 addPoint(log(1.645E+03), 5.421E-01);
788 addPoint(log(1.655E+03), 5.416E-01);
789 addPoint(log(1.665E+03), 5.410E-01);
790 addPoint(log(1.695E+03), 5.392E-01);
791 addPoint(log(1.705E+03), 5.387E-01);
792 addPoint(log(1.715E+03), 5.381E-01);
793 addPoint(log(1.735E+03), 5.369E-01);
794 addPoint(log(1.755E+03), 5.358E-01);
795 addPoint(log(1.765E+03), 5.354E-01);
796 addPoint(log(1.785E+03), 5.345E-01);
797 addPoint(log(1.795E+03), 5.341E-01);
798 addPoint(log(1.815E+03), 5.332E-01);
799 addPoint(log(1.835E+03), 5.323E-01);
800 addPoint(log(1.845E+03), 5.319E-01);
801 addPoint(log(1.855E+03), 5.315E-01);
802 addPoint(log(1.885E+03), 5.302E-01);
803 addPoint(log(1.895E+03), 5.297E-01);
804 addPoint(log(1.905E+03), 5.293E-01);
805 addPoint(log(1.915E+03), 5.289E-01);
806 addPoint(log(1.945E+03), 5.276E-01);
807 addPoint(log(1.955E+03), 5.272E-01);
808 addPoint(log(1.965E+03), 5.267E-01);
809 addPoint(log(1.985E+03), 5.259E-01);
810 addPoint(log(2.005E+03), 5.250E-01);
811 addPoint(log(2.015E+03), 5.246E-01);
812 addPoint(log(2.035E+03), 5.237E-01);
813 addPoint(log(2.045E+03), 5.233E-01);
814 addPoint(log(2.065E+03), 5.225E-01);
815 addPoint(log(2.085E+03), 5.216E-01);
816 addPoint(log(2.095E+03), 5.212E-01);
817 addPoint(log(2.105E+03), 5.208E-01);
818 addPoint(log(2.135E+03), 5.195E-01);
819 addPoint(log(2.145E+03), 5.191E-01);
820 addPoint(log(2.155E+03), 5.187E-01);
821 addPoint(log(2.165E+03), 5.182E-01);
822 addPoint(log(2.195E+03), 5.170E-01);
823 addPoint(log(2.205E+03), 5.166E-01);
824 addPoint(log(2.215E+03), 5.161E-01);
825 addPoint(log(2.235E+03), 5.153E-01);
826 addPoint(log(2.255E+03), 5.145E-01);
827 addPoint(log(2.265E+03), 5.140E-01);
828 addPoint(log(2.285E+03), 5.137E-01);
829 addPoint(log(2.295E+03), 5.136E-01);
830 addPoint(log(2.315E+03), 5.134E-01);
831 addPoint(log(2.335E+03), 5.132E-01);
832 addPoint(log(2.345E+03), 5.131E-01);
833 addPoint(log(2.355E+03), 5.130E-01);
834 addPoint(log(2.385E+03), 5.126E-01);
835 addPoint(log(2.395E+03), 5.125E-01);
836 addPoint(log(2.405E+03), 5.124E-01);
837 addPoint(log(2.415E+03), 5.123E-01);
838 addPoint(log(2.435E+03), 5.121E-01);
839 addPoint(log(2.445E+03), 5.120E-01);
840 addPoint(log(2.455E+03), 5.119E-01);
841 addPoint(log(2.465E+03), 5.118E-01);
842 addPoint(log(2.485E+03), 5.115E-01);
843 addPoint(log(2.495E+03), 5.114E-01);
844 addPoint(log(2.505E+03), 5.113E-01);
845 addPoint(log(2.515E+03), 5.112E-01);
846 addPoint(log(2.535E+03), 5.110E-01);
847 addPoint(log(2.545E+03), 5.109E-01);
848 addPoint(log(2.555E+03), 5.108E-01);
849 addPoint(log(2.565E+03), 5.107E-01);
850 addPoint(log(2.585E+03), 5.105E-01);
851 addPoint(log(2.595E+03), 5.104E-01);
852 addPoint(log(2.605E+03), 5.102E-01);
853 addPoint(log(2.615E+03), 5.101E-01);
854 addPoint(log(2.635E+03), 5.099E-01);
855 addPoint(log(2.645E+03), 5.098E-01);
856 addPoint(log(2.655E+03), 5.097E-01);
857 addPoint(log(2.665E+03), 5.096E-01);
858 addPoint(log(2.685E+03), 5.094E-01);
859 addPoint(log(2.695E+03), 5.093E-01);
860 addPoint(log(2.705E+03), 5.092E-01);
861 addPoint(log(2.715E+03), 5.091E-01);
862 addPoint(log(2.735E+03), 5.088E-01);
863 addPoint(log(2.745E+03), 5.087E-01);
864 addPoint(log(2.755E+03), 5.086E-01);
865 addPoint(log(2.765E+03), 5.085E-01);
866 addPoint(log(2.785E+03), 5.083E-01);
867 addPoint(log(2.795E+03), 5.082E-01);
868 addPoint(log(2.805E+03), 5.081E-01);
869 addPoint(log(2.815E+03), 5.080E-01);
870 addPoint(log(2.835E+03), 5.078E-01);
871 addPoint(log(2.845E+03), 5.077E-01);
872 addPoint(log(2.855E+03), 5.076E-01);
873 addPoint(log(2.865E+03), 5.074E-01);
874 addPoint(log(2.885E+03), 5.072E-01);
875 addPoint(log(2.895E+03), 5.071E-01);
876 addPoint(log(2.905E+03), 5.070E-01);
877 addPoint(log(2.915E+03), 5.069E-01);
878 addPoint(log(2.935E+03), 5.067E-01);
879 addPoint(log(2.945E+03), 5.066E-01);
880 addPoint(log(2.955E+03), 5.065E-01);
881 addPoint(log(2.965E+03), 5.064E-01);
882 addPoint(log(2.985E+03), 5.062E-01);
883 addPoint(log(2.995E+03), 5.060E-01);
884 addPoint(log(3.005E+03), 5.059E-01);
885 addPoint(log(3.015E+03), 5.058E-01);
886 addPoint(log(3.035E+03), 5.056E-01);
887 addPoint(log(3.045E+03), 5.055E-01);
888 addPoint(log(3.055E+03), 5.054E-01);
889 addPoint(log(3.065E+03), 5.053E-01);
890 addPoint(log(3.085E+03), 5.051E-01);
891 addPoint(log(3.095E+03), 5.050E-01);
892 addPoint(log(3.105E+03), 5.049E-01);
893 addPoint(log(3.115E+03), 5.048E-01);
894 addPoint(log(3.135E+03), 5.046E-01);
895 addPoint(log(3.145E+03), 5.044E-01);
896 addPoint(log(3.155E+03), 5.043E-01);
897 addPoint(log(3.165E+03), 5.042E-01);
898 addPoint(log(3.185E+03), 5.040E-01);
899 addPoint(log(3.195E+03), 5.039E-01);
900 addPoint(log(3.205E+03), 5.038E-01);
901 addPoint(log(3.215E+03), 5.037E-01);
902 addPoint(log(3.235E+03), 5.035E-01);
903 addPoint(log(3.245E+03), 5.034E-01);
904 addPoint(log(3.255E+03), 5.033E-01);
905 addPoint(log(3.265E+03), 5.032E-01);
906 addPoint(log(3.285E+03), 5.030E-01);
907 addPoint(log(3.295E+03), 5.028E-01);
908 addPoint(log(3.305E+03), 5.027E-01);
909 addPoint(log(3.315E+03), 5.026E-01);
910 addPoint(log(3.335E+03), 5.024E-01);
911 addPoint(log(3.345E+03), 5.023E-01);
912 addPoint(log(3.355E+03), 5.022E-01);
913 addPoint(log(3.365E+03), 5.021E-01);
914 addPoint(log(3.385E+03), 5.019E-01);
915 addPoint(log(3.395E+03), 5.018E-01);
916 addPoint(log(3.405E+03), 5.017E-01);
917 addPoint(log(3.415E+03), 5.016E-01);
918 addPoint(log(3.435E+03), 5.014E-01);
919 addPoint(log(3.445E+03), 5.013E-01);
920 addPoint(log(3.455E+03), 5.011E-01);
921 addPoint(log(3.465E+03), 5.010E-01);
922 addPoint(log(3.485E+03), 5.008E-01);
923 addPoint(log(3.495E+03), 5.007E-01);
924 addPoint(log(3.505E+03), 5.006E-01);
925 addPoint(log(3.515E+03), 5.005E-01);
926 addPoint(log(3.535E+03), 5.003E-01);
927 addPoint(log(3.545E+03), 5.002E-01);
928 addPoint(log(3.555E+03), 5.001E-01);
929 addPoint(log(3.565E+03), 5.000E-01);
930 addPoint(log(3.585E+03), 4.998E-01);
931 addPoint(log(3.595E+03), 4.997E-01);
932 addPoint(log(3.605E+03), 4.996E-01);
933 addPoint(log(3.615E+03), 4.995E-01);
934 addPoint(log(3.635E+03), 4.992E-01);
935 addPoint(log(3.645E+03), 4.991E-01);
936 addPoint(log(3.655E+03), 4.990E-01);
937 addPoint(log(3.665E+03), 4.989E-01);
938 addPoint(log(3.685E+03), 4.987E-01);
939 addPoint(log(3.695E+03), 4.986E-01);
940 addPoint(log(3.705E+03), 4.985E-01);
941 addPoint(log(3.715E+03), 4.984E-01);
942 addPoint(log(3.735E+03), 4.982E-01);
943 addPoint(log(3.745E+03), 4.981E-01);
944 addPoint(log(3.755E+03), 4.980E-01);
945 addPoint(log(3.765E+03), 4.979E-01);
946 addPoint(log(3.785E+03), 4.977E-01);
947 addPoint(log(3.795E+03), 4.976E-01);
948 addPoint(log(3.805E+03), 4.975E-01);
949 addPoint(log(3.815E+03), 4.974E-01);
950 addPoint(log(3.835E+03), 4.971E-01);
951 addPoint(log(3.845E+03), 4.970E-01);
952 addPoint(log(3.855E+03), 4.969E-01);
953 addPoint(log(3.865E+03), 4.968E-01);
954 addPoint(log(3.885E+03), 4.966E-01);
955 addPoint(log(3.895E+03), 4.965E-01);
956 addPoint(log(3.905E+03), 4.964E-01);
957 addPoint(log(3.915E+03), 4.963E-01);
958 addPoint(log(3.935E+03), 4.961E-01);
959 addPoint(log(3.945E+03), 4.960E-01);
960 addPoint(log(3.955E+03), 4.959E-01);
961 addPoint(log(3.965E+03), 4.958E-01);
962 addPoint(log(3.985E+03), 4.956E-01);
963 addPoint(log(3.995E+03), 4.955E-01);
964 addPoint(log(4.005E+03), 4.954E-01);
965 addPoint(log(4.015E+03), 4.953E-01);
966 addPoint(log(4.025E+03), 4.952E-01);
967 addPoint(log(4.035E+03), 4.950E-01);
968 addPoint(log(4.045E+03), 4.949E-01);
969 addPoint(log(4.055E+03), 4.948E-01);
970 addPoint(log(4.065E+03), 4.947E-01);
971 addPoint(log(4.075E+03), 4.946E-01);
972 addPoint(log(4.085E+03), 4.945E-01);
973 addPoint(log(4.095E+03), 4.944E-01);
974 addPoint(log(4.105E+03), 4.943E-01);
975 addPoint(log(4.115E+03), 4.942E-01);
976 addPoint(log(4.125E+03), 4.941E-01);
977 addPoint(log(4.135E+03), 4.940E-01);
978 addPoint(log(4.145E+03), 4.939E-01);
979 addPoint(log(4.155E+03), 4.938E-01);
980 addPoint(log(4.165E+03), 4.937E-01);
981 addPoint(log(4.175E+03), 4.936E-01);
982 addPoint(log(4.185E+03), 4.935E-01);
983 addPoint(log(4.195E+03), 4.934E-01);
984 addPoint(log(4.205E+03), 4.933E-01);
985 addPoint(log(4.215E+03), 4.932E-01);
986 addPoint(log(4.225E+03), 4.931E-01);
987 addPoint(log(4.235E+03), 4.930E-01);
988 addPoint(log(4.245E+03), 4.929E-01);
989 addPoint(log(4.255E+03), 4.928E-01);
990 addPoint(log(4.265E+03), 4.927E-01);
991 addPoint(log(4.275E+03), 4.925E-01);
992 addPoint(log(4.285E+03), 4.924E-01);
993 addPoint(log(4.295E+03), 4.923E-01);
994 addPoint(log(4.305E+03), 4.922E-01);
995 addPoint(log(4.315E+03), 4.921E-01);
996 addPoint(log(4.325E+03), 4.920E-01);
997 addPoint(log(4.335E+03), 4.919E-01);
998 addPoint(log(4.345E+03), 4.918E-01);
999 addPoint(log(4.355E+03), 4.917E-01);
1000 addPoint(log(4.365E+03), 4.916E-01);
1001 addPoint(log(4.375E+03), 4.915E-01);
1002 addPoint(log(4.385E+03), 4.914E-01);
1003 addPoint(log(4.395E+03), 4.913E-01);
1004 addPoint(log(4.405E+03), 4.912E-01);
1005 addPoint(log(4.415E+03), 4.911E-01);
1006 addPoint(log(4.425E+03), 4.910E-01);
1007 addPoint(log(4.435E+03), 4.909E-01);
1008 addPoint(log(4.445E+03), 4.908E-01);
1009 addPoint(log(4.455E+03), 4.907E-01);
1010 addPoint(log(4.465E+03), 4.906E-01);
1011 addPoint(log(4.475E+03), 4.905E-01);
1012 addPoint(log(4.485E+03), 4.904E-01);
1013 addPoint(log(4.495E+03), 4.903E-01);
1014 addPoint(log(4.505E+03), 4.902E-01);
1015 addPoint(log(4.515E+03), 4.901E-01);
1016 addPoint(log(4.525E+03), 4.900E-01);
1017 addPoint(log(4.535E+03), 4.899E-01);
1018 addPoint(log(4.545E+03), 4.898E-01);
1019 addPoint(log(4.555E+03), 4.897E-01);
1020 addPoint(log(4.565E+03), 4.895E-01);
1021 addPoint(log(4.575E+03), 4.894E-01);
1022 addPoint(log(4.585E+03), 4.893E-01);
1023 addPoint(log(4.595E+03), 4.892E-01);
1024 addPoint(log(4.605E+03), 4.891E-01);
1025 addPoint(log(4.615E+03), 4.890E-01);
1026 addPoint(log(4.625E+03), 4.889E-01);
1027 addPoint(log(4.635E+03), 4.888E-01);
1028 addPoint(log(4.645E+03), 4.887E-01);
1029 addPoint(log(4.655E+03), 4.886E-01);
1030 addPoint(log(4.665E+03), 4.885E-01);
1031 addPoint(log(4.675E+03), 4.884E-01);
1032 addPoint(log(4.685E+03), 4.883E-01);
1033 addPoint(log(4.695E+03), 4.882E-01);
1034 addPoint(log(4.705E+03), 4.881E-01);
1035 addPoint(log(4.715E+03), 4.880E-01);
1036 addPoint(log(4.725E+03), 4.879E-01);
1037 addPoint(log(4.735E+03), 4.878E-01);
1038 addPoint(log(4.745E+03), 4.877E-01);
1039 addPoint(log(4.755E+03), 4.876E-01);
1040 addPoint(log(4.765E+03), 4.875E-01);
1041 addPoint(log(4.775E+03), 4.874E-01);
1042 addPoint(log(4.785E+03), 4.873E-01);
1043 addPoint(log(4.795E+03), 4.872E-01);
1044 addPoint(log(4.805E+03), 4.871E-01);
1045 addPoint(log(4.815E+03), 4.870E-01);
1046 addPoint(log(4.825E+03), 4.869E-01);
1047 addPoint(log(4.835E+03), 4.868E-01);
1048 addPoint(log(4.845E+03), 4.867E-01);
1049 addPoint(log(4.855E+03), 4.866E-01);
1050 addPoint(log(4.865E+03), 4.865E-01);
1051 addPoint(log(4.875E+03), 4.864E-01);
1052 addPoint(log(4.885E+03), 4.863E-01);
1053 addPoint(log(4.895E+03), 4.862E-01);
1054 addPoint(log(4.905E+03), 4.861E-01);
1055 addPoint(log(4.915E+03), 4.860E-01);
1056 addPoint(log(4.925E+03), 4.859E-01);
1057 addPoint(log(4.935E+03), 4.858E-01);
1058 addPoint(log(4.945E+03), 4.857E-01);
1059 addPoint(log(4.955E+03), 4.856E-01);
1060 addPoint(log(4.965E+03), 4.854E-01);
1061 addPoint(log(4.975E+03), 4.853E-01);
1062 addPoint(log(4.985E+03), 4.852E-01);
1063 addPoint(log(4.995E+03), 4.851E-01);
1064 addPoint(log(5.005E+03), 4.850E-01);
1065 addPoint(log(5.015E+03), 4.849E-01);
1066 addPoint(log(5.025E+03), 4.848E-01);
1067 addPoint(log(5.035E+03), 4.847E-01);
1068 addPoint(log(5.045E+03), 4.846E-01);
1069 addPoint(log(5.055E+03), 4.845E-01);
1070 addPoint(log(5.065E+03), 4.844E-01);
1071 addPoint(log(5.075E+03), 4.843E-01);
1072 addPoint(log(5.085E+03), 4.842E-01);
1073 addPoint(log(5.095E+03), 4.841E-01);
1074 addPoint(log(5.105E+03), 4.840E-01);
1075 addPoint(log(5.115E+03), 4.839E-01);
1076 addPoint(log(5.125E+03), 4.838E-01);
1077 addPoint(log(5.135E+03), 4.837E-01);
1078 addPoint(log(5.145E+03), 4.836E-01);
1079 addPoint(log(5.155E+03), 4.835E-01);
1080 addPoint(log(5.165E+03), 4.834E-01);
1081 addPoint(log(5.175E+03), 4.833E-01);
1082 addPoint(log(5.185E+03), 4.832E-01);
1083 addPoint(log(5.195E+03), 4.831E-01);
1084 addPoint(log(5.205E+03), 4.830E-01);
1085 addPoint(log(5.215E+03), 4.829E-01);
1086 addPoint(log(5.225E+03), 4.828E-01);
1087 addPoint(log(5.235E+03), 4.827E-01);
1088 addPoint(log(5.245E+03), 4.826E-01);
1089 addPoint(log(5.255E+03), 4.825E-01);
1090 addPoint(log(5.265E+03), 4.824E-01);
1091 addPoint(log(5.275E+03), 4.823E-01);
1092 addPoint(log(5.285E+03), 4.822E-01);
1093 addPoint(log(5.295E+03), 4.821E-01);
1094 addPoint(log(5.305E+03), 4.820E-01);
1095 addPoint(log(5.315E+03), 4.819E-01);
1096 addPoint(log(5.325E+03), 4.818E-01);
1097 addPoint(log(5.335E+03), 4.817E-01);
1098 addPoint(log(5.345E+03), 4.816E-01);
1099 addPoint(log(5.355E+03), 4.815E-01);
1100 addPoint(log(5.365E+03), 4.814E-01);
1101 addPoint(log(5.375E+03), 4.813E-01);
1102 addPoint(log(5.385E+03), 4.812E-01);
1103 addPoint(log(5.395E+03), 4.811E-01);
1104 addPoint(log(5.405E+03), 4.810E-01);
1105 addPoint(log(5.415E+03), 4.809E-01);
1106 addPoint(log(5.425E+03), 4.808E-01);
1107 addPoint(log(5.435E+03), 4.807E-01);
1108 addPoint(log(5.445E+03), 4.806E-01);
1109 addPoint(log(5.455E+03), 4.805E-01);
1110 addPoint(log(5.465E+03), 4.804E-01);
1111 addPoint(log(5.475E+03), 4.803E-01);
1112 addPoint(log(5.485E+03), 4.802E-01);
1113 addPoint(log(5.495E+03), 4.801E-01);
1114 addPoint(log(5.505E+03), 4.800E-01);
1115 addPoint(log(5.515E+03), 4.799E-01);
1116 addPoint(log(5.525E+03), 4.798E-01);
1117 addPoint(log(5.535E+03), 4.797E-01);
1118 addPoint(log(5.545E+03), 4.796E-01);
1119 addPoint(log(5.555E+03), 4.795E-01);
1120 addPoint(log(5.565E+03), 4.794E-01);
1121 addPoint(log(5.575E+03), 4.793E-01);
1122 addPoint(log(5.585E+03), 4.792E-01);
1123 addPoint(log(5.595E+03), 4.791E-01);
1124 addPoint(log(5.605E+03), 4.790E-01);
1125 addPoint(log(5.615E+03), 4.789E-01);
1126 addPoint(log(5.625E+03), 4.788E-01);
1127 addPoint(log(5.635E+03), 4.787E-01);
1128 addPoint(log(5.645E+03), 4.786E-01);
1129 addPoint(log(5.655E+03), 4.785E-01);
1130 addPoint(log(5.665E+03), 4.784E-01);
1131 addPoint(log(5.675E+03), 4.783E-01);
1132 addPoint(log(5.685E+03), 4.782E-01);
1133 addPoint(log(5.695E+03), 4.781E-01);
1134 addPoint(log(5.705E+03), 4.780E-01);
1135 addPoint(log(5.715E+03), 4.779E-01);
1136 addPoint(log(5.725E+03), 4.778E-01);
1137 addPoint(log(5.735E+03), 4.777E-01);
1138 addPoint(log(5.745E+03), 4.776E-01);
1139 addPoint(log(5.755E+03), 4.775E-01);
1140 addPoint(log(5.765E+03), 4.774E-01);
1141 addPoint(log(5.775E+03), 4.773E-01);
1142 addPoint(log(5.785E+03), 4.772E-01);
1143 addPoint(log(5.795E+03), 4.771E-01);
1144 addPoint(log(5.805E+03), 4.770E-01);
1145 addPoint(log(5.815E+03), 4.769E-01);
1146 addPoint(log(5.825E+03), 4.768E-01);
1147 addPoint(log(5.835E+03), 4.767E-01);
1148 addPoint(log(5.845E+03), 4.766E-01);
1149 addPoint(log(5.855E+03), 4.765E-01);
1150 addPoint(log(5.865E+03), 4.764E-01);
1151 addPoint(log(5.875E+03), 4.763E-01);
1152 addPoint(log(5.885E+03), 4.762E-01);
1153 addPoint(log(5.895E+03), 4.761E-01);
1154 addPoint(log(5.905E+03), 4.760E-01);
1155 addPoint(log(5.915E+03), 4.759E-01);
1156 addPoint(log(5.925E+03), 4.758E-01);
1157 addPoint(log(5.935E+03), 4.757E-01);
1158 addPoint(log(5.945E+03), 4.756E-01);
1159 addPoint(log(5.955E+03), 4.755E-01);
1160 addPoint(log(5.965E+03), 4.754E-01);
1161 addPoint(log(5.975E+03), 4.753E-01);
1162 addPoint(log(5.985E+03), 4.752E-01);
1163 addPoint(log(5.995E+03), 4.751E-01);
1164 addPoint(log(6.005E+03), 4.750E-01);
1165 addPoint(log(6.015E+03), 4.749E-01);
1166 addPoint(log(6.025E+03), 4.748E-01);
1167 addPoint(log(6.035E+03), 4.747E-01);
1168 addPoint(log(6.045E+03), 4.746E-01);
1169 addPoint(log(6.055E+03), 4.745E-01);
1170 addPoint(log(6.065E+03), 4.744E-01);
1171 addPoint(log(6.075E+03), 4.743E-01);
1172 addPoint(log(6.085E+03), 4.742E-01);
1173 addPoint(log(6.095E+03), 4.741E-01);
1174 addPoint(log(6.105E+03), 4.740E-01);
1175 addPoint(log(6.115E+03), 4.739E-01);
1176 addPoint(log(6.125E+03), 4.738E-01);
1177 addPoint(log(6.135E+03), 4.737E-01);
1178 addPoint(log(6.145E+03), 4.736E-01);
1179 addPoint(log(6.155E+03), 4.735E-01);
1180 addPoint(log(6.165E+03), 4.734E-01);
1181 addPoint(log(6.175E+03), 4.733E-01);
1182 addPoint(log(6.185E+03), 4.732E-01);
1183 addPoint(log(6.195E+03), 4.731E-01);
1184 addPoint(log(6.205E+03), 4.730E-01);
1185 addPoint(log(6.215E+03), 4.729E-01);
1186 addPoint(log(6.225E+03), 4.728E-01);
1187 addPoint(log(6.235E+03), 4.727E-01);
1188 addPoint(log(6.245E+03), 4.726E-01);
1189 addPoint(log(6.255E+03), 4.725E-01);
1190 addPoint(log(6.265E+03), 4.724E-01);
1191 addPoint(log(6.275E+03), 4.723E-01);
1192 addPoint(log(6.285E+03), 4.722E-01);
1193 addPoint(log(6.295E+03), 4.721E-01);
1194 addPoint(log(6.305E+03), 4.720E-01);
1195 addPoint(log(6.315E+03), 4.719E-01);
1196 addPoint(log(6.325E+03), 4.718E-01);
1197 addPoint(log(6.335E+03), 4.717E-01);
1198 addPoint(log(6.345E+03), 4.716E-01);
1199 addPoint(log(6.355E+03), 4.715E-01);
1200 addPoint(log(6.365E+03), 4.714E-01);
1201 addPoint(log(6.375E+03), 4.713E-01);
1202 addPoint(log(6.385E+03), 4.712E-01);
1203 addPoint(log(6.395E+03), 4.711E-01);
1204 addPoint(log(6.405E+03), 4.710E-01);
1205 addPoint(log(6.415E+03), 4.709E-01);
1206 addPoint(log(6.425E+03), 4.708E-01);
1207 addPoint(log(6.435E+03), 4.707E-01);
1208 addPoint(log(6.445E+03), 4.706E-01);
1209 addPoint(log(6.455E+03), 4.705E-01);
1210 addPoint(log(6.465E+03), 4.704E-01);
1211 addPoint(log(6.475E+03), 4.703E-01);
1212 addPoint(log(6.485E+03), 4.702E-01);
1213 addPoint(log(6.495E+03), 4.701E-01);
1214 addPoint(log(6.505E+03), 4.700E-01);
1215 addPoint(log(6.515E+03), 4.699E-01);
1216 addPoint(log(6.525E+03), 4.698E-01);
1217 addPoint(log(6.535E+03), 4.697E-01);
1218 addPoint(log(6.545E+03), 4.696E-01);
1219 addPoint(log(6.555E+03), 4.695E-01);
1220 addPoint(log(6.565E+03), 4.695E-01);
1221 addPoint(log(6.575E+03), 4.694E-01);
1222 addPoint(log(6.585E+03), 4.693E-01);
1223 addPoint(log(6.595E+03), 4.692E-01);
1224 addPoint(log(6.605E+03), 4.691E-01);
1225 addPoint(log(6.615E+03), 4.690E-01);
1226 addPoint(log(6.625E+03), 4.689E-01);
1227 addPoint(log(6.635E+03), 4.688E-01);
1228 addPoint(log(6.645E+03), 4.687E-01);
1229 addPoint(log(6.655E+03), 4.686E-01);
1230 addPoint(log(6.665E+03), 4.685E-01);
1231 addPoint(log(6.675E+03), 4.684E-01);
1232 addPoint(log(6.685E+03), 4.683E-01);
1233 addPoint(log(6.695E+03), 4.682E-01);
1234 addPoint(log(6.705E+03), 4.681E-01);
1235 addPoint(log(6.715E+03), 4.680E-01);
1236 addPoint(log(6.725E+03), 4.679E-01);
1237 addPoint(log(6.735E+03), 4.678E-01);
1238 addPoint(log(6.745E+03), 4.677E-01);
1239 addPoint(log(6.755E+03), 4.676E-01);
1240 addPoint(log(6.765E+03), 4.675E-01);
1241 addPoint(log(6.775E+03), 4.674E-01);
1242 addPoint(log(6.785E+03), 4.673E-01);
1243 addPoint(log(6.795E+03), 4.672E-01);
1244 addPoint(log(6.805E+03), 4.671E-01);
1245 addPoint(log(6.815E+03), 4.670E-01);
1246 addPoint(log(6.825E+03), 4.669E-01);
1247 addPoint(log(6.835E+03), 4.668E-01);
1248 addPoint(log(6.845E+03), 4.667E-01);
1249 addPoint(log(6.855E+03), 4.666E-01);
1250 addPoint(log(6.865E+03), 4.665E-01);
1251 addPoint(log(6.875E+03), 4.664E-01);
1252 addPoint(log(6.885E+03), 4.663E-01);
1253 addPoint(log(6.895E+03), 4.662E-01);
1254 addPoint(log(6.905E+03), 4.661E-01);
1255 addPoint(log(6.915E+03), 4.660E-01);
1256 addPoint(log(6.925E+03), 4.659E-01);
1257 addPoint(log(6.935E+03), 4.658E-01);
1258 addPoint(log(6.945E+03), 4.657E-01);
1259 addPoint(log(6.955E+03), 4.656E-01);
1260 addPoint(log(6.965E+03), 4.656E-01);
1261 addPoint(log(6.975E+03), 4.655E-01);
1262 addPoint(log(6.985E+03), 4.654E-01);
1263 addPoint(log(6.995E+03), 4.653E-01);
1264 addPoint(log(7.005E+03), 4.652E-01);
1265 addPoint(log(7.015E+03), 4.651E-01);
1266 addPoint(log(7.025E+03), 4.650E-01);
1267 addPoint(log(7.035E+03), 4.649E-01);
1268 addPoint(log(7.045E+03), 4.648E-01);
1269 addPoint(log(7.055E+03), 4.647E-01);
1270 addPoint(log(7.065E+03), 4.646E-01);
1271 addPoint(log(7.075E+03), 4.645E-01);
1272 addPoint(log(7.085E+03), 4.644E-01);
1273 addPoint(log(7.095E+03), 4.643E-01);
1274 addPoint(log(7.105E+03), 4.642E-01);
1275 addPoint(log(7.115E+03), 4.641E-01);
1276 addPoint(log(7.125E+03), 4.640E-01);
1277 addPoint(log(7.135E+03), 4.639E-01);
1278 addPoint(log(7.145E+03), 4.638E-01);
1279 addPoint(log(7.155E+03), 4.637E-01);
1280 addPoint(log(7.165E+03), 4.636E-01);
1281 addPoint(log(7.185E+03), 4.634E-01);
1282 addPoint(log(7.195E+03), 4.633E-01);
1283 addPoint(log(7.205E+03), 4.632E-01);
1284 addPoint(log(7.225E+03), 4.630E-01);
1285 addPoint(log(7.235E+03), 4.629E-01);
1286 addPoint(log(7.245E+03), 4.628E-01);
1287 addPoint(log(7.255E+03), 4.627E-01);
1288 addPoint(log(7.275E+03), 4.626E-01);
1289 addPoint(log(7.285E+03), 4.625E-01);
1290 addPoint(log(7.295E+03), 4.624E-01);
1291 addPoint(log(7.315E+03), 4.622E-01);
1292 addPoint(log(7.325E+03), 4.621E-01);
1293 addPoint(log(7.335E+03), 4.620E-01);
1294 addPoint(log(7.345E+03), 4.619E-01);
1295 addPoint(log(7.365E+03), 4.617E-01);
1296 addPoint(log(7.375E+03), 4.616E-01);
1297 addPoint(log(7.385E+03), 4.615E-01);
1298 addPoint(log(7.395E+03), 4.614E-01);
1299 addPoint(log(7.415E+03), 4.612E-01);
1300 addPoint(log(7.425E+03), 4.611E-01);
1301 addPoint(log(7.435E+03), 4.610E-01);
1302 addPoint(log(7.445E+03), 4.609E-01);
1303 addPoint(log(7.465E+03), 4.607E-01);
1304 addPoint(log(7.475E+03), 4.606E-01);
1305 addPoint(log(7.485E+03), 4.605E-01);
1306 addPoint(log(7.495E+03), 4.604E-01);
1307 addPoint(log(7.515E+03), 4.603E-01);
1308 addPoint(log(7.525E+03), 4.602E-01);
1309 addPoint(log(7.535E+03), 4.601E-01);
1310 addPoint(log(7.545E+03), 4.600E-01);
1311 addPoint(log(7.565E+03), 4.598E-01);
1312 addPoint(log(7.575E+03), 4.597E-01);
1313 addPoint(log(7.585E+03), 4.596E-01);
1314 addPoint(log(7.595E+03), 4.595E-01);
1315 addPoint(log(7.615E+03), 4.593E-01);
1316 addPoint(log(7.625E+03), 4.592E-01);
1317 addPoint(log(7.635E+03), 4.591E-01);
1318 addPoint(log(7.645E+03), 4.590E-01);
1319 addPoint(log(7.665E+03), 4.588E-01);
1320 addPoint(log(7.675E+03), 4.587E-01);
1321 addPoint(log(7.685E+03), 4.586E-01);
1322 addPoint(log(7.695E+03), 4.585E-01);
1323 addPoint(log(7.715E+03), 4.583E-01);
1324 addPoint(log(7.725E+03), 4.582E-01);
1325 addPoint(log(7.735E+03), 4.582E-01);
1326 addPoint(log(7.745E+03), 4.581E-01);
1327 addPoint(log(7.765E+03), 4.579E-01);
1328 addPoint(log(7.775E+03), 4.578E-01);
1329 addPoint(log(7.785E+03), 4.577E-01);
1330 addPoint(log(7.795E+03), 4.576E-01);
1331 addPoint(log(7.815E+03), 4.574E-01);
1332 addPoint(log(7.825E+03), 4.573E-01);
1333 addPoint(log(7.835E+03), 4.572E-01);
1334 addPoint(log(7.845E+03), 4.571E-01);
1335 addPoint(log(7.865E+03), 4.569E-01);
1336 addPoint(log(7.875E+03), 4.568E-01);
1337 addPoint(log(7.885E+03), 4.567E-01);
1338 addPoint(log(7.895E+03), 4.566E-01);
1339 addPoint(log(7.915E+03), 4.564E-01);
1340 addPoint(log(7.925E+03), 4.564E-01);
1341 addPoint(log(7.935E+03), 4.563E-01);
1342 addPoint(log(7.945E+03), 4.562E-01);
1343 addPoint(log(7.965E+03), 4.560E-01);
1344 addPoint(log(7.975E+03), 4.559E-01);
1345 addPoint(log(7.985E+03), 4.558E-01);
1346 addPoint(log(7.995E+03), 4.557E-01);
1347 addPoint(log(8.015E+03), 4.555E-01);
1348 addPoint(log(8.025E+03), 4.554E-01);
1349 addPoint(log(8.035E+03), 4.553E-01);
1350 addPoint(log(8.045E+03), 4.552E-01);
1351 addPoint(log(8.065E+03), 4.550E-01);
1352 addPoint(log(8.075E+03), 4.549E-01);
1353 addPoint(log(8.085E+03), 4.548E-01);
1354 addPoint(log(8.095E+03), 4.547E-01);
1355 addPoint(log(8.115E+03), 4.546E-01);
1356 addPoint(log(8.125E+03), 4.545E-01);
1357 addPoint(log(8.135E+03), 4.544E-01);
1358 addPoint(log(8.145E+03), 4.543E-01);
1359 addPoint(log(8.165E+03), 4.541E-01);
1360 addPoint(log(8.175E+03), 4.540E-01);
1361 addPoint(log(8.185E+03), 4.539E-01);
1362 addPoint(log(8.195E+03), 4.538E-01);
1363 addPoint(log(8.215E+03), 4.536E-01);
1364 addPoint(log(8.225E+03), 4.535E-01);
1365 addPoint(log(8.235E+03), 4.534E-01);
1366 addPoint(log(8.245E+03), 4.533E-01);
1367 addPoint(log(8.265E+03), 4.531E-01);
1368 addPoint(log(8.275E+03), 4.531E-01);
1369 addPoint(log(8.285E+03), 4.530E-01);
1370 addPoint(log(8.295E+03), 4.529E-01);
1371 addPoint(log(8.315E+03), 4.527E-01);
1372 addPoint(log(8.325E+03), 4.526E-01);
1373 addPoint(log(8.335E+03), 4.525E-01);
1374 addPoint(log(8.345E+03), 4.524E-01);
1375 addPoint(log(8.365E+03), 4.522E-01);
1376 addPoint(log(8.375E+03), 4.521E-01);
1377 addPoint(log(8.385E+03), 4.520E-01);
1378 addPoint(log(8.395E+03), 4.519E-01);
1379 addPoint(log(8.415E+03), 4.517E-01);
1380 addPoint(log(8.425E+03), 4.516E-01);
1381 addPoint(log(8.435E+03), 4.516E-01);
1382 addPoint(log(8.445E+03), 4.515E-01);
1383 addPoint(log(8.465E+03), 4.513E-01);
1384 addPoint(log(8.475E+03), 4.512E-01);
1385 addPoint(log(8.485E+03), 4.511E-01);
1386 addPoint(log(8.495E+03), 4.510E-01);
1387 addPoint(log(8.515E+03), 4.508E-01);
1388 addPoint(log(8.525E+03), 4.507E-01);
1389 addPoint(log(8.535E+03), 4.506E-01);
1390 addPoint(log(8.545E+03), 4.505E-01);
1391 addPoint(log(8.565E+03), 4.503E-01);
1392 addPoint(log(8.575E+03), 4.502E-01);
1393 addPoint(log(8.585E+03), 4.502E-01);
1394 addPoint(log(8.595E+03), 4.501E-01);
1395 addPoint(log(8.615E+03), 4.499E-01);
1396 addPoint(log(8.625E+03), 4.498E-01);
1397 addPoint(log(8.635E+03), 4.497E-01);
1398 addPoint(log(8.645E+03), 4.496E-01);
1399 addPoint(log(8.665E+03), 4.494E-01);
1400 addPoint(log(8.675E+03), 4.493E-01);
1401 addPoint(log(8.685E+03), 4.492E-01);
1402 addPoint(log(8.695E+03), 4.491E-01);
1403 addPoint(log(8.715E+03), 4.489E-01);
1404 addPoint(log(8.725E+03), 4.489E-01);
1405 addPoint(log(8.735E+03), 4.488E-01);
1406 addPoint(log(8.745E+03), 4.487E-01);
1407 addPoint(log(8.765E+03), 4.485E-01);
1408 addPoint(log(8.775E+03), 4.484E-01);
1409 addPoint(log(8.785E+03), 4.483E-01);
1410 addPoint(log(8.795E+03), 4.482E-01);
1411 addPoint(log(8.815E+03), 4.480E-01);
1412 addPoint(log(8.825E+03), 4.479E-01);
1413 addPoint(log(8.835E+03), 4.478E-01);
1414 addPoint(log(8.845E+03), 4.477E-01);
1415 addPoint(log(8.865E+03), 4.476E-01);
1416 addPoint(log(8.875E+03), 4.475E-01);
1417 addPoint(log(8.885E+03), 4.474E-01);
1418 addPoint(log(8.895E+03), 4.473E-01);
1419 addPoint(log(8.915E+03), 4.471E-01);
1420 addPoint(log(8.925E+03), 4.470E-01);
1421 addPoint(log(8.935E+03), 4.469E-01);
1422 addPoint(log(8.945E+03), 4.468E-01);
1423 addPoint(log(8.965E+03), 4.466E-01);
1424 addPoint(log(8.975E+03), 4.465E-01);
1425 addPoint(log(8.995E+03), 4.464E-01);
1426 doSpline();
1427 };
1428
1429 void BkgNeutronWeight::addPoint(double x, double y)
1430 {
1431 m_x.push_back(x);
1432 m_y.push_back(y);
1433 m_b.push_back(0.0);
1434 m_c.push_back(0.0);
1435 m_d.push_back(0.0);
1436 }
1437
1438 void BkgNeutronWeight::doSpline()
1439 {
1440 // converted from fortran in Bevington's Data Analysis ... ugh!
1441 m_n = m_x.size();
1442 int n = m_n;
1443 int i = 0;
1444
1445 m_d[0] = m_x[1] - m_x[0];
1446 m_c[1] = (m_y[1] - m_y[0]) / m_d[0];
1447
1448 for (i = 1; i < n - 1; i++) {
1449 m_d[i] = m_x[i + 1] - m_x[i];
1450 m_b[i] = 2.0 * (m_d[i - 1] + m_d[i]);
1451 m_c[i + 1] = (m_y[i + 1] - m_y[i]) / m_d[i];
1452 m_c[i] = m_c[i + 1] - m_c[i];
1453 }
1454
1455 m_b[0] = -m_d[0];
1456 m_b[n - 1] = -m_d[n - 2];
1457
1458 m_c[0] = m_c[2] / (m_x[3] - m_x[1]) - m_c[1] / (m_x[2] - m_x[0]);
1459 m_c[0] = m_c[0] * m_d[0] * m_d[0] / (m_x[3] - m_x[0]);
1460 m_c[n - 1] = m_c[n - 2] / (m_x[n - 1] - m_x[n - 3]) - m_c[n - 3] / (m_x[n - 2] - m_x[n - 4]);
1461 m_c[n - 1] = -m_c[n - 1] * m_d[n - 2] * m_d[n - 2] / (m_x[n - 1] - m_x[n - 4]);
1462
1463 for (i = 1; i < n; i++) {
1464 double t = m_d[i - 1] / m_b[i - 1];
1465 m_b[i] = m_b[i] - t * m_d[i - 1];
1466 m_c[i] = m_c[i] - t * m_c[i - 1];
1467 }
1468
1469 m_c[n - 1] = m_c[n - 1] / m_b[n - 1];
1470
1471 for (i = n - 2; i >= 0; i--) {
1472 m_c[i] = (m_c[i] - m_d[i] * m_c[i + 1]) / m_b[i];
1473 }
1474
1475 m_b[n - 1] = (m_y[n - 1] - m_y[n - 2]) / m_d[n - 2] + m_d[n - 2] * (m_c[n - 2] + 2.0 * m_c[n - 1]);
1476
1477 for (i = 0; i < n - 1; i++) {
1478 m_b[i] = (m_y[i + 1] - m_y[i]) / m_d[i] - m_d[i] * (m_c[i + 1] + 2.0 * m_c[i]);
1479 m_d[i] = (m_c[i + 1] - m_c[i]) / m_d[i];
1480 m_c[i] = 3.0 * m_c[i];
1481 }
1482
1483 m_c[n - 1] = 3.0 * m_c[n - 1];
1484 m_d[n - 1] = m_d[n - 2];
1485
1486 }
1487
1489} // end of namespace Belle2
Belle2::BkgNeutronWeight
The class to get the weighting factor for a 1-MeV-equivalent neutron flux on silicon.
Definition: BkgNeutronWeight.h:23
Belle2::BkgNeutronWeight::m_b
std::vector< double > m_b
Vector of first derivatives in spline fit (calculated)
Definition: BkgNeutronWeight.h:54
Belle2::BkgNeutronWeight::m_c
std::vector< double > m_c
Vector of second derivatives in spline fit (calculated)
Definition: BkgNeutronWeight.h:57
Belle2::BkgNeutronWeight::m_n
int m_n
Number of tabulated points.
Definition: BkgNeutronWeight.h:45
Belle2::BkgNeutronWeight::m_y
std::vector< double > m_y
Vector of tabulated y coordinates (weighting factors)
Definition: BkgNeutronWeight.h:51
Belle2::BkgNeutronWeight::m_d
std::vector< double > m_d
Vector of third derivatives in spline fit (calculated)
Definition: BkgNeutronWeight.h:60
Belle2::BkgNeutronWeight::m_x
std::vector< double > m_x
Vector of tabulated x coordinates (logarithm of kinetic energy)
Definition: BkgNeutronWeight.h:48
Belle2::BkgNeutronWeight::getWeight
double getWeight(double ke)
Get weighting factor to convert a neutron to its 1-MeV equivalent.
Definition: BkgNeutronWeight.cc:26
Belle2::BkgNeutronWeight::addPoint
void addPoint(double, double)
Function to add one (x,y) data point to the spline table.
Definition: BkgNeutronWeight.cc:1429
Belle2::BkgNeutronWeight::m_Instance
static BkgNeutronWeight * m_Instance
static pointer to the singleton instance of this class
Definition: BkgNeutronWeight.h:63
Belle2::BkgNeutronWeight::getInstance
static BkgNeutronWeight & getInstance()
Return a reference to the singleton BkgNeutronWeight instance.
Definition: BkgNeutronWeight.cc:20
Belle2::BkgNeutronWeight::doSpline
void doSpline()
Function to calculate the cubic spline derivatives.
Definition: BkgNeutronWeight.cc:1438
Belle2::BkgNeutronWeight::BkgNeutronWeight
BkgNeutronWeight(void)
Constructor is private for this singleton.
Definition: BkgNeutronWeight.cc:43
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17