Belle II Software  release-08-00-10
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 
12 namespace 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