Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JPMTTransitTimeGenerator.hh
Go to the documentation of this file.
1 #ifndef __JDETECTOR__JPMTTRANSITTIMEGENERATOR__
2 #define __JDETECTOR__JPMTTRANSITTIMEGENERATOR__
3 
4 #include "TRandom3.h"
5 
7 #include "JLang/JComparator.hh"
8 #include "JLang/JComparison.hh"
9 
10 
11 /**
12  * \author mdejong
13  */
14 
15 namespace JDETECTOR {}
16 namespace JPP { using namespace JDETECTOR; }
17 
18 namespace JDETECTOR {
19 
20  /**
21  * Auxiliary class to generate PMT transition times based on measurements
22  * by Oleg Kalekin et al. (see also examples/JDetector/JTTS.cc).
23  *
24  * For a given random number, this function returns a transition time [ns].
25  */
28  {
29  public:
30  /**
31  * Default constructor.
32  */
34  JTOOLS::JGridSplineFunction1D_t()
35  {
36  (*this)[ 0.000] = -10.955 ;
37  (*this)[+0.001] = -6.5027 ;
38  (*this)[+0.002] = -5.5647 ;
39  (*this)[+0.003] = -5.0281 ;
40  (*this)[+0.004] = -4.8482 ;
41  (*this)[+0.005] = -4.7248 ;
42  (*this)[+0.006] = -4.6013 ;
43  (*this)[+0.007] = -4.4778 ;
44  (*this)[+0.008] = -4.3544 ;
45  (*this)[+0.009] = -4.2309 ;
46  (*this)[+0.01 ] = -4.1074 ;
47  (*this)[+0.011] = -3.984 ;
48  (*this)[+0.012] = -3.9284 ;
49  (*this)[+0.013] = -3.8937 ;
50  (*this)[+0.014] = -3.859 ;
51  (*this)[+0.015] = -3.8243 ;
52  (*this)[+0.016] = -3.7896 ;
53  (*this)[+0.017] = -3.755 ;
54  (*this)[+0.018] = -3.7203 ;
55  (*this)[+0.019] = -3.6856 ;
56  (*this)[+0.02 ] = -3.6509 ;
57  (*this)[+0.021] = -3.6163 ;
58  (*this)[+0.022] = -3.5816 ;
59  (*this)[+0.023] = -3.5469 ;
60  (*this)[+0.024] = -3.5122 ;
61  (*this)[+0.025] = -3.4775 ;
62  (*this)[+0.026] = -3.4429 ;
63  (*this)[+0.027] = -3.4082 ;
64  (*this)[+0.028] = -3.3735 ;
65  (*this)[+0.029] = -3.3388 ;
66  (*this)[+0.03 ] = -3.3041 ;
67  (*this)[+0.031] = -3.2695 ;
68  (*this)[+0.032] = -3.2348 ;
69  (*this)[+0.033] = -3.2001 ;
70  (*this)[+0.034] = -3.1654 ;
71  (*this)[+0.035] = -3.1307 ;
72  (*this)[+0.036] = -3.0961 ;
73  (*this)[+0.037] = -3.0614 ;
74  (*this)[+0.038] = -3.0267 ;
75  (*this)[+0.039] = -2.992 ;
76  (*this)[+0.04 ] = -2.9573 ;
77  (*this)[+0.041] = -2.9424 ;
78  (*this)[+0.042] = -2.9289 ;
79  (*this)[+0.043] = -2.9155 ;
80  (*this)[+0.044] = -2.902 ;
81  (*this)[+0.045] = -2.8886 ;
82  (*this)[+0.046] = -2.8751 ;
83  (*this)[+0.047] = -2.8616 ;
84  (*this)[+0.048] = -2.8482 ;
85  (*this)[+0.049] = -2.8347 ;
86  (*this)[+0.05 ] = -2.8213 ;
87  (*this)[+0.051] = -2.8078 ;
88  (*this)[+0.052] = -2.7944 ;
89  (*this)[+0.053] = -2.7809 ;
90  (*this)[+0.054] = -2.7674 ;
91  (*this)[+0.055] = -2.754 ;
92  (*this)[+0.056] = -2.7405 ;
93  (*this)[+0.057] = -2.7271 ;
94  (*this)[+0.058] = -2.7136 ;
95  (*this)[+0.059] = -2.7002 ;
96  (*this)[+0.06 ] = -2.6867 ;
97  (*this)[+0.061] = -2.6732 ;
98  (*this)[+0.062] = -2.6598 ;
99  (*this)[+0.063] = -2.6463 ;
100  (*this)[+0.064] = -2.6329 ;
101  (*this)[+0.065] = -2.6194 ;
102  (*this)[+0.066] = -2.606 ;
103  (*this)[+0.067] = -2.5925 ;
104  (*this)[+0.068] = -2.579 ;
105  (*this)[+0.069] = -2.5656 ;
106  (*this)[+0.07 ] = -2.5521 ;
107  (*this)[+0.071] = -2.5387 ;
108  (*this)[+0.072] = -2.5252 ;
109  (*this)[+0.073] = -2.5118 ;
110  (*this)[+0.074] = -2.4983 ;
111  (*this)[+0.075] = -2.4848 ;
112  (*this)[+0.076] = -2.4714 ;
113  (*this)[+0.077] = -2.4579 ;
114  (*this)[+0.078] = -2.4445 ;
115  (*this)[+0.079] = -2.431 ;
116  (*this)[+0.08 ] = -2.4176 ;
117  (*this)[+0.081] = -2.4041 ;
118  (*this)[+0.082] = -2.3906 ;
119  (*this)[+0.083] = -2.3772 ;
120  (*this)[+0.084] = -2.3637 ;
121  (*this)[+0.085] = -2.3503 ;
122  (*this)[+0.086] = -2.3368 ;
123  (*this)[+0.087] = -2.3234 ;
124  (*this)[+0.088] = -2.3099 ;
125  (*this)[+0.089] = -2.2965 ;
126  (*this)[+0.09 ] = -2.283 ;
127  (*this)[+0.091] = -2.2695 ;
128  (*this)[+0.092] = -2.2561 ;
129  (*this)[+0.093] = -2.2426 ;
130  (*this)[+0.094] = -2.2292 ;
131  (*this)[+0.095] = -2.2157 ;
132  (*this)[+0.096] = -2.2023 ;
133  (*this)[+0.097] = -2.1888 ;
134  (*this)[+0.098] = -2.1753 ;
135  (*this)[+0.099] = -2.1619 ;
136  (*this)[+0.1 ] = -2.1484 ;
137  (*this)[+0.101] = -2.135 ;
138  (*this)[+0.102] = -2.1215 ;
139  (*this)[+0.103] = -2.1081 ;
140  (*this)[+0.104] = -2.0946 ;
141  (*this)[+0.105] = -2.0811 ;
142  (*this)[+0.106] = -2.0677 ;
143  (*this)[+0.107] = -2.0542 ;
144  (*this)[+0.108] = -2.0408 ;
145  (*this)[+0.109] = -2.0273 ;
146  (*this)[+0.11 ] = -2.0139 ;
147  (*this)[+0.111] = -2.0004 ;
148  (*this)[+0.112] = -1.9869 ;
149  (*this)[+0.113] = -1.9735 ;
150  (*this)[+0.114] = -1.96 ;
151  (*this)[+0.115] = -1.9504 ;
152  (*this)[+0.116] = -1.9431 ;
153  (*this)[+0.117] = -1.9358 ;
154  (*this)[+0.118] = -1.9286 ;
155  (*this)[+0.119] = -1.9213 ;
156  (*this)[+0.12 ] = -1.914 ;
157  (*this)[+0.121] = -1.9067 ;
158  (*this)[+0.122] = -1.8995 ;
159  (*this)[+0.123] = -1.8922 ;
160  (*this)[+0.124] = -1.8849 ;
161  (*this)[+0.125] = -1.8776 ;
162  (*this)[+0.126] = -1.8703 ;
163  (*this)[+0.127] = -1.8631 ;
164  (*this)[+0.128] = -1.8558 ;
165  (*this)[+0.129] = -1.8485 ;
166  (*this)[+0.13 ] = -1.8412 ;
167  (*this)[+0.131] = -1.834 ;
168  (*this)[+0.132] = -1.8267 ;
169  (*this)[+0.133] = -1.8194 ;
170  (*this)[+0.134] = -1.8121 ;
171  (*this)[+0.135] = -1.8049 ;
172  (*this)[+0.136] = -1.7976 ;
173  (*this)[+0.137] = -1.7903 ;
174  (*this)[+0.138] = -1.783 ;
175  (*this)[+0.139] = -1.7758 ;
176  (*this)[+0.14 ] = -1.7685 ;
177  (*this)[+0.141] = -1.7612 ;
178  (*this)[+0.142] = -1.7539 ;
179  (*this)[+0.143] = -1.7467 ;
180  (*this)[+0.144] = -1.7394 ;
181  (*this)[+0.145] = -1.7321 ;
182  (*this)[+0.146] = -1.7248 ;
183  (*this)[+0.147] = -1.7175 ;
184  (*this)[+0.148] = -1.7103 ;
185  (*this)[+0.149] = -1.703 ;
186  (*this)[+0.15 ] = -1.6957 ;
187  (*this)[+0.151] = -1.6884 ;
188  (*this)[+0.152] = -1.6812 ;
189  (*this)[+0.153] = -1.6739 ;
190  (*this)[+0.154] = -1.6666 ;
191  (*this)[+0.155] = -1.6593 ;
192  (*this)[+0.156] = -1.6521 ;
193  (*this)[+0.157] = -1.6448 ;
194  (*this)[+0.158] = -1.6375 ;
195  (*this)[+0.159] = -1.6302 ;
196  (*this)[+0.16 ] = -1.623 ;
197  (*this)[+0.161] = -1.6157 ;
198  (*this)[+0.162] = -1.6084 ;
199  (*this)[+0.163] = -1.6011 ;
200  (*this)[+0.164] = -1.5939 ;
201  (*this)[+0.165] = -1.5866 ;
202  (*this)[+0.166] = -1.5793 ;
203  (*this)[+0.167] = -1.572 ;
204  (*this)[+0.168] = -1.5647 ;
205  (*this)[+0.169] = -1.5575 ;
206  (*this)[+0.17 ] = -1.5502 ;
207  (*this)[+0.171] = -1.5429 ;
208  (*this)[+0.172] = -1.5356 ;
209  (*this)[+0.173] = -1.5284 ;
210  (*this)[+0.174] = -1.5211 ;
211  (*this)[+0.175] = -1.5138 ;
212  (*this)[+0.176] = -1.5065 ;
213  (*this)[+0.177] = -1.4993 ;
214  (*this)[+0.178] = -1.492 ;
215  (*this)[+0.179] = -1.4847 ;
216  (*this)[+0.18 ] = -1.4774 ;
217  (*this)[+0.181] = -1.4702 ;
218  (*this)[+0.182] = -1.4629 ;
219  (*this)[+0.183] = -1.4556 ;
220  (*this)[+0.184] = -1.4483 ;
221  (*this)[+0.185] = -1.441 ;
222  (*this)[+0.186] = -1.4338 ;
223  (*this)[+0.187] = -1.4265 ;
224  (*this)[+0.188] = -1.4192 ;
225  (*this)[+0.189] = -1.4119 ;
226  (*this)[+0.19 ] = -1.4047 ;
227  (*this)[+0.191] = -1.3974 ;
228  (*this)[+0.192] = -1.3901 ;
229  (*this)[+0.193] = -1.3828 ;
230  (*this)[+0.194] = -1.3756 ;
231  (*this)[+0.195] = -1.3683 ;
232  (*this)[+0.196] = -1.361 ;
233  (*this)[+0.197] = -1.3537 ;
234  (*this)[+0.198] = -1.3465 ;
235  (*this)[+0.199] = -1.3392 ;
236  (*this)[+0.2 ] = -1.3319 ;
237  (*this)[+0.201] = -1.3246 ;
238  (*this)[+0.202] = -1.3174 ;
239  (*this)[+0.203] = -1.3101 ;
240  (*this)[+0.204] = -1.3028 ;
241  (*this)[+0.205] = -1.2955 ;
242  (*this)[+0.206] = -1.2882 ;
243  (*this)[+0.207] = -1.281 ;
244  (*this)[+0.208] = -1.2737 ;
245  (*this)[+0.209] = -1.2664 ;
246  (*this)[+0.21 ] = -1.2591 ;
247  (*this)[+0.211] = -1.2519 ;
248  (*this)[+0.212] = -1.2446 ;
249  (*this)[+0.213] = -1.2373 ;
250  (*this)[+0.214] = -1.23 ;
251  (*this)[+0.215] = -1.2228 ;
252  (*this)[+0.216] = -1.2155 ;
253  (*this)[+0.217] = -1.2082 ;
254  (*this)[+0.218] = -1.2009 ;
255  (*this)[+0.219] = -1.1937 ;
256  (*this)[+0.22 ] = -1.1864 ;
257  (*this)[+0.221] = -1.1791 ;
258  (*this)[+0.222] = -1.1718 ;
259  (*this)[+0.223] = -1.1645 ;
260  (*this)[+0.224] = -1.1573 ;
261  (*this)[+0.225] = -1.15 ;
262  (*this)[+0.226] = -1.1427 ;
263  (*this)[+0.227] = -1.1354 ;
264  (*this)[+0.228] = -1.1282 ;
265  (*this)[+0.229] = -1.1209 ;
266  (*this)[+0.23 ] = -1.1136 ;
267  (*this)[+0.231] = -1.1063 ;
268  (*this)[+0.232] = -1.0991 ;
269  (*this)[+0.233] = -1.0918 ;
270  (*this)[+0.234] = -1.0845 ;
271  (*this)[+0.235] = -1.0772 ;
272  (*this)[+0.236] = -1.07 ;
273  (*this)[+0.237] = -1.0627 ;
274  (*this)[+0.238] = -1.0554 ;
275  (*this)[+0.239] = -1.0481 ;
276  (*this)[+0.24 ] = -1.0409 ;
277  (*this)[+0.241] = -1.0336 ;
278  (*this)[+0.242] = -1.0263 ;
279  (*this)[+0.243] = -1.019 ;
280  (*this)[+0.244] = -1.0117 ;
281  (*this)[+0.245] = -1.0045 ;
282  (*this)[+0.246] = -0.99719 ;
283  (*this)[+0.247] = -0.98992 ;
284  (*this)[+0.248] = -0.98264 ;
285  (*this)[+0.249] = -0.97537 ;
286  (*this)[+0.25 ] = -0.96809 ;
287  (*this)[+0.251] = -0.96081 ;
288  (*this)[+0.252] = -0.95388 ;
289  (*this)[+0.253] = -0.94846 ;
290  (*this)[+0.254] = -0.94304 ;
291  (*this)[+0.255] = -0.93761 ;
292  (*this)[+0.256] = -0.93219 ;
293  (*this)[+0.257] = -0.92677 ;
294  (*this)[+0.258] = -0.92135 ;
295  (*this)[+0.259] = -0.91593 ;
296  (*this)[+0.26 ] = -0.91051 ;
297  (*this)[+0.261] = -0.90508 ;
298  (*this)[+0.262] = -0.89966 ;
299  (*this)[+0.263] = -0.89424 ;
300  (*this)[+0.264] = -0.88882 ;
301  (*this)[+0.265] = -0.8834 ;
302  (*this)[+0.266] = -0.87798 ;
303  (*this)[+0.267] = -0.87255 ;
304  (*this)[+0.268] = -0.86713 ;
305  (*this)[+0.269] = -0.86171 ;
306  (*this)[+0.27 ] = -0.85629 ;
307  (*this)[+0.271] = -0.85087 ;
308  (*this)[+0.272] = -0.84545 ;
309  (*this)[+0.273] = -0.84002 ;
310  (*this)[+0.274] = -0.8346 ;
311  (*this)[+0.275] = -0.82918 ;
312  (*this)[+0.276] = -0.82376 ;
313  (*this)[+0.277] = -0.81834 ;
314  (*this)[+0.278] = -0.81292 ;
315  (*this)[+0.279] = -0.80749 ;
316  (*this)[+0.28 ] = -0.80207 ;
317  (*this)[+0.281] = -0.79665 ;
318  (*this)[+0.282] = -0.79123 ;
319  (*this)[+0.283] = -0.78581 ;
320  (*this)[+0.284] = -0.78038 ;
321  (*this)[+0.285] = -0.77496 ;
322  (*this)[+0.286] = -0.76954 ;
323  (*this)[+0.287] = -0.76412 ;
324  (*this)[+0.288] = -0.7587 ;
325  (*this)[+0.289] = -0.75328 ;
326  (*this)[+0.29 ] = -0.74785 ;
327  (*this)[+0.291] = -0.74243 ;
328  (*this)[+0.292] = -0.73701 ;
329  (*this)[+0.293] = -0.73159 ;
330  (*this)[+0.294] = -0.72617 ;
331  (*this)[+0.295] = -0.72075 ;
332  (*this)[+0.296] = -0.71532 ;
333  (*this)[+0.297] = -0.7099 ;
334  (*this)[+0.298] = -0.70448 ;
335  (*this)[+0.299] = -0.69906 ;
336  (*this)[+0.3 ] = -0.69364 ;
337  (*this)[+0.301] = -0.68822 ;
338  (*this)[+0.302] = -0.68279 ;
339  (*this)[+0.303] = -0.67737 ;
340  (*this)[+0.304] = -0.67195 ;
341  (*this)[+0.305] = -0.66653 ;
342  (*this)[+0.306] = -0.66111 ;
343  (*this)[+0.307] = -0.65569 ;
344  (*this)[+0.308] = -0.65026 ;
345  (*this)[+0.309] = -0.64484 ;
346  (*this)[+0.31 ] = -0.63942 ;
347  (*this)[+0.311] = -0.634 ;
348  (*this)[+0.312] = -0.62858 ;
349  (*this)[+0.313] = -0.62316 ;
350  (*this)[+0.314] = -0.61773 ;
351  (*this)[+0.315] = -0.61231 ;
352  (*this)[+0.316] = -0.60689 ;
353  (*this)[+0.317] = -0.60147 ;
354  (*this)[+0.318] = -0.59605 ;
355  (*this)[+0.319] = -0.59063 ;
356  (*this)[+0.32 ] = -0.5852 ;
357  (*this)[+0.321] = -0.57978 ;
358  (*this)[+0.322] = -0.57436 ;
359  (*this)[+0.323] = -0.56894 ;
360  (*this)[+0.324] = -0.56352 ;
361  (*this)[+0.325] = -0.5581 ;
362  (*this)[+0.326] = -0.55267 ;
363  (*this)[+0.327] = -0.54725 ;
364  (*this)[+0.328] = -0.54183 ;
365  (*this)[+0.329] = -0.53641 ;
366  (*this)[+0.33 ] = -0.53099 ;
367  (*this)[+0.331] = -0.52557 ;
368  (*this)[+0.332] = -0.52014 ;
369  (*this)[+0.333] = -0.51472 ;
370  (*this)[+0.334] = -0.5093 ;
371  (*this)[+0.335] = -0.50388 ;
372  (*this)[+0.336] = -0.49846 ;
373  (*this)[+0.337] = -0.49304 ;
374  (*this)[+0.338] = -0.48761 ;
375  (*this)[+0.339] = -0.48219 ;
376  (*this)[+0.34 ] = -0.47677 ;
377  (*this)[+0.341] = -0.47135 ;
378  (*this)[+0.342] = -0.46593 ;
379  (*this)[+0.343] = -0.46051 ;
380  (*this)[+0.344] = -0.45508 ;
381  (*this)[+0.345] = -0.44966 ;
382  (*this)[+0.346] = -0.44424 ;
383  (*this)[+0.347] = -0.43882 ;
384  (*this)[+0.348] = -0.4334 ;
385  (*this)[+0.349] = -0.42798 ;
386  (*this)[+0.35 ] = -0.42255 ;
387  (*this)[+0.351] = -0.41713 ;
388  (*this)[+0.352] = -0.41171 ;
389  (*this)[+0.353] = -0.40629 ;
390  (*this)[+0.354] = -0.40087 ;
391  (*this)[+0.355] = -0.39545 ;
392  (*this)[+0.356] = -0.39002 ;
393  (*this)[+0.357] = -0.3846 ;
394  (*this)[+0.358] = -0.37918 ;
395  (*this)[+0.359] = -0.37376 ;
396  (*this)[+0.36 ] = -0.36834 ;
397  (*this)[+0.361] = -0.36292 ;
398  (*this)[+0.362] = -0.35749 ;
399  (*this)[+0.363] = -0.35207 ;
400  (*this)[+0.364] = -0.34665 ;
401  (*this)[+0.365] = -0.34123 ;
402  (*this)[+0.366] = -0.33581 ;
403  (*this)[+0.367] = -0.33039 ;
404  (*this)[+0.368] = -0.32496 ;
405  (*this)[+0.369] = -0.31954 ;
406  (*this)[+0.37 ] = -0.31412 ;
407  (*this)[+0.371] = -0.3087 ;
408  (*this)[+0.372] = -0.30328 ;
409  (*this)[+0.373] = -0.29786 ;
410  (*this)[+0.374] = -0.29243 ;
411  (*this)[+0.375] = -0.28701 ;
412  (*this)[+0.376] = -0.28159 ;
413  (*this)[+0.377] = -0.27617 ;
414  (*this)[+0.378] = -0.27075 ;
415  (*this)[+0.379] = -0.26533 ;
416  (*this)[+0.38 ] = -0.2599 ;
417  (*this)[+0.381] = -0.25448 ;
418  (*this)[+0.382] = -0.24906 ;
419  (*this)[+0.383] = -0.24364 ;
420  (*this)[+0.384] = -0.23822 ;
421  (*this)[+0.385] = -0.2328 ;
422  (*this)[+0.386] = -0.22737 ;
423  (*this)[+0.387] = -0.22195 ;
424  (*this)[+0.388] = -0.21653 ;
425  (*this)[+0.389] = -0.21111 ;
426  (*this)[+0.39 ] = -0.20569 ;
427  (*this)[+0.391] = -0.20027 ;
428  (*this)[+0.392] = -0.19484 ;
429  (*this)[+0.393] = -0.18942 ;
430  (*this)[+0.394] = -0.184 ;
431  (*this)[+0.395] = -0.17858 ;
432  (*this)[+0.396] = -0.17316 ;
433  (*this)[+0.397] = -0.16774 ;
434  (*this)[+0.398] = -0.16231 ;
435  (*this)[+0.399] = -0.15689 ;
436  (*this)[+0.4 ] = -0.15147 ;
437  (*this)[+0.401] = -0.14605 ;
438  (*this)[+0.402] = -0.14063 ;
439  (*this)[+0.403] = -0.13521 ;
440  (*this)[+0.404] = -0.12978 ;
441  (*this)[+0.405] = -0.12436 ;
442  (*this)[+0.406] = -0.11894 ;
443  (*this)[+0.407] = -0.11352 ;
444  (*this)[+0.408] = -0.1081 ;
445  (*this)[+0.409] = -0.10268 ;
446  (*this)[+0.41 ] = -0.097253;
447  (*this)[+0.411] = -0.091832;
448  (*this)[+0.412] = -0.08641 ;
449  (*this)[+0.413] = -0.080988;
450  (*this)[+0.414] = -0.075567;
451  (*this)[+0.415] = -0.070145;
452  (*this)[+0.416] = -0.064723;
453  (*this)[+0.417] = -0.059302;
454  (*this)[+0.418] = -0.05388 ;
455  (*this)[+0.419] = -0.048458;
456  (*this)[+0.42 ] = -0.043037;
457  (*this)[+0.421] = -0.037615;
458  (*this)[+0.422] = -0.032193;
459  (*this)[+0.423] = -0.026772;
460  (*this)[+0.424] = -0.02135 ;
461  (*this)[+0.425] = -0.015928;
462  (*this)[+0.426] = -0.010506;
463  (*this)[+0.427] = -0.0050848;
464  (*this)[+0.428] = +0.00033689;
465  (*this)[+0.429] = +0.0057586;
466  (*this)[+0.43 ] = +0.01118 ;
467  (*this)[+0.431] = +0.016602;
468  (*this)[+0.432] = +0.022024;
469  (*this)[+0.433] = +0.027445;
470  (*this)[+0.434] = +0.032867;
471  (*this)[+0.435] = +0.038289;
472  (*this)[+0.436] = +0.04371 ;
473  (*this)[+0.437] = +0.049075;
474  (*this)[+0.438] = +0.05442 ;
475  (*this)[+0.439] = +0.059765;
476  (*this)[+0.44 ] = +0.06511 ;
477  (*this)[+0.441] = +0.070454;
478  (*this)[+0.442] = +0.075799;
479  (*this)[+0.443] = +0.081144;
480  (*this)[+0.444] = +0.086489;
481  (*this)[+0.445] = +0.091834;
482  (*this)[+0.446] = +0.097179;
483  (*this)[+0.447] = +0.10252 ;
484  (*this)[+0.448] = +0.10787 ;
485  (*this)[+0.449] = +0.11321 ;
486  (*this)[+0.45 ] = +0.11856 ;
487  (*this)[+0.451] = +0.1239 ;
488  (*this)[+0.452] = +0.12925 ;
489  (*this)[+0.453] = +0.13459 ;
490  (*this)[+0.454] = +0.13994 ;
491  (*this)[+0.455] = +0.14528 ;
492  (*this)[+0.456] = +0.15063 ;
493  (*this)[+0.457] = +0.15597 ;
494  (*this)[+0.458] = +0.16132 ;
495  (*this)[+0.459] = +0.16666 ;
496  (*this)[+0.46 ] = +0.17201 ;
497  (*this)[+0.461] = +0.17735 ;
498  (*this)[+0.462] = +0.1827 ;
499  (*this)[+0.463] = +0.18804 ;
500  (*this)[+0.464] = +0.19339 ;
501  (*this)[+0.465] = +0.19873 ;
502  (*this)[+0.466] = +0.20407 ;
503  (*this)[+0.467] = +0.20942 ;
504  (*this)[+0.468] = +0.21476 ;
505  (*this)[+0.469] = +0.22011 ;
506  (*this)[+0.47 ] = +0.22545 ;
507  (*this)[+0.471] = +0.2308 ;
508  (*this)[+0.472] = +0.23614 ;
509  (*this)[+0.473] = +0.24149 ;
510  (*this)[+0.474] = +0.24683 ;
511  (*this)[+0.475] = +0.25218 ;
512  (*this)[+0.476] = +0.25752 ;
513  (*this)[+0.477] = +0.26287 ;
514  (*this)[+0.478] = +0.26821 ;
515  (*this)[+0.479] = +0.27356 ;
516  (*this)[+0.48 ] = +0.2789 ;
517  (*this)[+0.481] = +0.28425 ;
518  (*this)[+0.482] = +0.28959 ;
519  (*this)[+0.483] = +0.29494 ;
520  (*this)[+0.484] = +0.30028 ;
521  (*this)[+0.485] = +0.30563 ;
522  (*this)[+0.486] = +0.31097 ;
523  (*this)[+0.487] = +0.31632 ;
524  (*this)[+0.488] = +0.32166 ;
525  (*this)[+0.489] = +0.32701 ;
526  (*this)[+0.49 ] = +0.33235 ;
527  (*this)[+0.491] = +0.3377 ;
528  (*this)[+0.492] = +0.34304 ;
529  (*this)[+0.493] = +0.34838 ;
530  (*this)[+0.494] = +0.35373 ;
531  (*this)[+0.495] = +0.35907 ;
532  (*this)[+0.496] = +0.36442 ;
533  (*this)[+0.497] = +0.36976 ;
534  (*this)[+0.498] = +0.37511 ;
535  (*this)[+0.499] = +0.38045 ;
536  (*this)[+0.5 ] = +0.3858 ;
537  (*this)[+0.501] = +0.39114 ;
538  (*this)[+0.502] = +0.39649 ;
539  (*this)[+0.503] = +0.40183 ;
540  (*this)[+0.504] = +0.40718 ;
541  (*this)[+0.505] = +0.41252 ;
542  (*this)[+0.506] = +0.41787 ;
543  (*this)[+0.507] = +0.42321 ;
544  (*this)[+0.508] = +0.42856 ;
545  (*this)[+0.509] = +0.4339 ;
546  (*this)[+0.51 ] = +0.43925 ;
547  (*this)[+0.511] = +0.44459 ;
548  (*this)[+0.512] = +0.44994 ;
549  (*this)[+0.513] = +0.45528 ;
550  (*this)[+0.514] = +0.46063 ;
551  (*this)[+0.515] = +0.46597 ;
552  (*this)[+0.516] = +0.47132 ;
553  (*this)[+0.517] = +0.47666 ;
554  (*this)[+0.518] = +0.48201 ;
555  (*this)[+0.519] = +0.48735 ;
556  (*this)[+0.52 ] = +0.4927 ;
557  (*this)[+0.521] = +0.49804 ;
558  (*this)[+0.522] = +0.50338 ;
559  (*this)[+0.523] = +0.50873 ;
560  (*this)[+0.524] = +0.51407 ;
561  (*this)[+0.525] = +0.51942 ;
562  (*this)[+0.526] = +0.52476 ;
563  (*this)[+0.527] = +0.53011 ;
564  (*this)[+0.528] = +0.53545 ;
565  (*this)[+0.529] = +0.5408 ;
566  (*this)[+0.53 ] = +0.54614 ;
567  (*this)[+0.531] = +0.55149 ;
568  (*this)[+0.532] = +0.55683 ;
569  (*this)[+0.533] = +0.56218 ;
570  (*this)[+0.534] = +0.56752 ;
571  (*this)[+0.535] = +0.57287 ;
572  (*this)[+0.536] = +0.57821 ;
573  (*this)[+0.537] = +0.58356 ;
574  (*this)[+0.538] = +0.5889 ;
575  (*this)[+0.539] = +0.59425 ;
576  (*this)[+0.54 ] = +0.59959 ;
577  (*this)[+0.541] = +0.60494 ;
578  (*this)[+0.542] = +0.61028 ;
579  (*this)[+0.543] = +0.61563 ;
580  (*this)[+0.544] = +0.62097 ;
581  (*this)[+0.545] = +0.62632 ;
582  (*this)[+0.546] = +0.63166 ;
583  (*this)[+0.547] = +0.63701 ;
584  (*this)[+0.548] = +0.64235 ;
585  (*this)[+0.549] = +0.64769 ;
586  (*this)[+0.55 ] = +0.65304 ;
587  (*this)[+0.551] = +0.65838 ;
588  (*this)[+0.552] = +0.66373 ;
589  (*this)[+0.553] = +0.66907 ;
590  (*this)[+0.554] = +0.67442 ;
591  (*this)[+0.555] = +0.67976 ;
592  (*this)[+0.556] = +0.68511 ;
593  (*this)[+0.557] = +0.69045 ;
594  (*this)[+0.558] = +0.6958 ;
595  (*this)[+0.559] = +0.70114 ;
596  (*this)[+0.56 ] = +0.70649 ;
597  (*this)[+0.561] = +0.71183 ;
598  (*this)[+0.562] = +0.71718 ;
599  (*this)[+0.563] = +0.72252 ;
600  (*this)[+0.564] = +0.72787 ;
601  (*this)[+0.565] = +0.73321 ;
602  (*this)[+0.566] = +0.73856 ;
603  (*this)[+0.567] = +0.7439 ;
604  (*this)[+0.568] = +0.74925 ;
605  (*this)[+0.569] = +0.75459 ;
606  (*this)[+0.57 ] = +0.75994 ;
607  (*this)[+0.571] = +0.76528 ;
608  (*this)[+0.572] = +0.77063 ;
609  (*this)[+0.573] = +0.77597 ;
610  (*this)[+0.574] = +0.78132 ;
611  (*this)[+0.575] = +0.78666 ;
612  (*this)[+0.576] = +0.792 ;
613  (*this)[+0.577] = +0.79735 ;
614  (*this)[+0.578] = +0.80269 ;
615  (*this)[+0.579] = +0.80804 ;
616  (*this)[+0.58 ] = +0.81338 ;
617  (*this)[+0.581] = +0.81873 ;
618  (*this)[+0.582] = +0.82407 ;
619  (*this)[+0.583] = +0.82942 ;
620  (*this)[+0.584] = +0.83476 ;
621  (*this)[+0.585] = +0.84011 ;
622  (*this)[+0.586] = +0.84545 ;
623  (*this)[+0.587] = +0.8508 ;
624  (*this)[+0.588] = +0.85614 ;
625  (*this)[+0.589] = +0.86149 ;
626  (*this)[+0.59 ] = +0.86683 ;
627  (*this)[+0.591] = +0.87218 ;
628  (*this)[+0.592] = +0.87752 ;
629  (*this)[+0.593] = +0.88287 ;
630  (*this)[+0.594] = +0.88821 ;
631  (*this)[+0.595] = +0.89356 ;
632  (*this)[+0.596] = +0.8989 ;
633  (*this)[+0.597] = +0.90425 ;
634  (*this)[+0.598] = +0.90959 ;
635  (*this)[+0.599] = +0.91494 ;
636  (*this)[+0.6 ] = +0.92028 ;
637  (*this)[+0.601] = +0.92563 ;
638  (*this)[+0.602] = +0.93097 ;
639  (*this)[+0.603] = +0.93631 ;
640  (*this)[+0.604] = +0.94166 ;
641  (*this)[+0.605] = +0.947 ;
642  (*this)[+0.606] = +0.95235 ;
643  (*this)[+0.607] = +0.95769 ;
644  (*this)[+0.608] = +0.96304 ;
645  (*this)[+0.609] = +0.96838 ;
646  (*this)[+0.61 ] = +0.97373 ;
647  (*this)[+0.611] = +0.97907 ;
648  (*this)[+0.612] = +0.98442 ;
649  (*this)[+0.613] = +0.98976 ;
650  (*this)[+0.614] = +0.99511 ;
651  (*this)[+0.615] = +1.0005 ;
652  (*this)[+0.616] = +1.0058 ;
653  (*this)[+0.617] = +1.0111 ;
654  (*this)[+0.618] = +1.0165 ;
655  (*this)[+0.619] = +1.0218 ;
656  (*this)[+0.62 ] = +1.0272 ;
657  (*this)[+0.621] = +1.0325 ;
658  (*this)[+0.622] = +1.0379 ;
659  (*this)[+0.623] = +1.0432 ;
660  (*this)[+0.624] = +1.0494 ;
661  (*this)[+0.625] = +1.056 ;
662  (*this)[+0.626] = +1.0626 ;
663  (*this)[+0.627] = +1.0692 ;
664  (*this)[+0.628] = +1.0758 ;
665  (*this)[+0.629] = +1.0824 ;
666  (*this)[+0.63 ] = +1.089 ;
667  (*this)[+0.631] = +1.0956 ;
668  (*this)[+0.632] = +1.1022 ;
669  (*this)[+0.633] = +1.1088 ;
670  (*this)[+0.634] = +1.1154 ;
671  (*this)[+0.635] = +1.122 ;
672  (*this)[+0.636] = +1.1286 ;
673  (*this)[+0.637] = +1.1352 ;
674  (*this)[+0.638] = +1.1418 ;
675  (*this)[+0.639] = +1.1484 ;
676  (*this)[+0.64 ] = +1.155 ;
677  (*this)[+0.641] = +1.1616 ;
678  (*this)[+0.642] = +1.1682 ;
679  (*this)[+0.643] = +1.1748 ;
680  (*this)[+0.644] = +1.1814 ;
681  (*this)[+0.645] = +1.188 ;
682  (*this)[+0.646] = +1.1946 ;
683  (*this)[+0.647] = +1.2012 ;
684  (*this)[+0.648] = +1.2078 ;
685  (*this)[+0.649] = +1.2144 ;
686  (*this)[+0.65 ] = +1.221 ;
687  (*this)[+0.651] = +1.2276 ;
688  (*this)[+0.652] = +1.2342 ;
689  (*this)[+0.653] = +1.2408 ;
690  (*this)[+0.654] = +1.2474 ;
691  (*this)[+0.655] = +1.254 ;
692  (*this)[+0.656] = +1.2606 ;
693  (*this)[+0.657] = +1.2672 ;
694  (*this)[+0.658] = +1.2738 ;
695  (*this)[+0.659] = +1.2804 ;
696  (*this)[+0.66 ] = +1.287 ;
697  (*this)[+0.661] = +1.2937 ;
698  (*this)[+0.662] = +1.3003 ;
699  (*this)[+0.663] = +1.3069 ;
700  (*this)[+0.664] = +1.3135 ;
701  (*this)[+0.665] = +1.3201 ;
702  (*this)[+0.666] = +1.3267 ;
703  (*this)[+0.667] = +1.3333 ;
704  (*this)[+0.668] = +1.3399 ;
705  (*this)[+0.669] = +1.3465 ;
706  (*this)[+0.67 ] = +1.3531 ;
707  (*this)[+0.671] = +1.3597 ;
708  (*this)[+0.672] = +1.3663 ;
709  (*this)[+0.673] = +1.3729 ;
710  (*this)[+0.674] = +1.3795 ;
711  (*this)[+0.675] = +1.3861 ;
712  (*this)[+0.676] = +1.3927 ;
713  (*this)[+0.677] = +1.3993 ;
714  (*this)[+0.678] = +1.4059 ;
715  (*this)[+0.679] = +1.4125 ;
716  (*this)[+0.68 ] = +1.4191 ;
717  (*this)[+0.681] = +1.4257 ;
718  (*this)[+0.682] = +1.4323 ;
719  (*this)[+0.683] = +1.4389 ;
720  (*this)[+0.684] = +1.4455 ;
721  (*this)[+0.685] = +1.4521 ;
722  (*this)[+0.686] = +1.4587 ;
723  (*this)[+0.687] = +1.4653 ;
724  (*this)[+0.688] = +1.4719 ;
725  (*this)[+0.689] = +1.4785 ;
726  (*this)[+0.69 ] = +1.4851 ;
727  (*this)[+0.691] = +1.4917 ;
728  (*this)[+0.692] = +1.4983 ;
729  (*this)[+0.693] = +1.5049 ;
730  (*this)[+0.694] = +1.5115 ;
731  (*this)[+0.695] = +1.5181 ;
732  (*this)[+0.696] = +1.5247 ;
733  (*this)[+0.697] = +1.5313 ;
734  (*this)[+0.698] = +1.5379 ;
735  (*this)[+0.699] = +1.5445 ;
736  (*this)[+0.7 ] = +1.5511 ;
737  (*this)[+0.701] = +1.5577 ;
738  (*this)[+0.702] = +1.5643 ;
739  (*this)[+0.703] = +1.571 ;
740  (*this)[+0.704] = +1.5776 ;
741  (*this)[+0.705] = +1.5842 ;
742  (*this)[+0.706] = +1.5908 ;
743  (*this)[+0.707] = +1.5974 ;
744  (*this)[+0.708] = +1.604 ;
745  (*this)[+0.709] = +1.6106 ;
746  (*this)[+0.71 ] = +1.6172 ;
747  (*this)[+0.711] = +1.6238 ;
748  (*this)[+0.712] = +1.6304 ;
749  (*this)[+0.713] = +1.637 ;
750  (*this)[+0.714] = +1.6436 ;
751  (*this)[+0.715] = +1.6502 ;
752  (*this)[+0.716] = +1.6568 ;
753  (*this)[+0.717] = +1.6634 ;
754  (*this)[+0.718] = +1.67 ;
755  (*this)[+0.719] = +1.6766 ;
756  (*this)[+0.72 ] = +1.6832 ;
757  (*this)[+0.721] = +1.6898 ;
758  (*this)[+0.722] = +1.6964 ;
759  (*this)[+0.723] = +1.703 ;
760  (*this)[+0.724] = +1.7096 ;
761  (*this)[+0.725] = +1.7162 ;
762  (*this)[+0.726] = +1.7228 ;
763  (*this)[+0.727] = +1.7294 ;
764  (*this)[+0.728] = +1.736 ;
765  (*this)[+0.729] = +1.7426 ;
766  (*this)[+0.73 ] = +1.7492 ;
767  (*this)[+0.731] = +1.7558 ;
768  (*this)[+0.732] = +1.7624 ;
769  (*this)[+0.733] = +1.769 ;
770  (*this)[+0.734] = +1.7756 ;
771  (*this)[+0.735] = +1.7822 ;
772  (*this)[+0.736] = +1.7888 ;
773  (*this)[+0.737] = +1.7954 ;
774  (*this)[+0.738] = +1.802 ;
775  (*this)[+0.739] = +1.8086 ;
776  (*this)[+0.74 ] = +1.8152 ;
777  (*this)[+0.741] = +1.8218 ;
778  (*this)[+0.742] = +1.8284 ;
779  (*this)[+0.743] = +1.835 ;
780  (*this)[+0.744] = +1.8416 ;
781  (*this)[+0.745] = +1.8483 ;
782  (*this)[+0.746] = +1.8549 ;
783  (*this)[+0.747] = +1.8615 ;
784  (*this)[+0.748] = +1.8681 ;
785  (*this)[+0.749] = +1.8747 ;
786  (*this)[+0.75 ] = +1.8813 ;
787  (*this)[+0.751] = +1.8879 ;
788  (*this)[+0.752] = +1.8945 ;
789  (*this)[+0.753] = +1.9011 ;
790  (*this)[+0.754] = +1.9077 ;
791  (*this)[+0.755] = +1.9143 ;
792  (*this)[+0.756] = +1.9209 ;
793  (*this)[+0.757] = +1.9275 ;
794  (*this)[+0.758] = +1.9341 ;
795  (*this)[+0.759] = +1.9407 ;
796  (*this)[+0.76 ] = +1.9473 ;
797  (*this)[+0.761] = +1.9539 ;
798  (*this)[+0.762] = +1.9605 ;
799  (*this)[+0.763] = +1.9671 ;
800  (*this)[+0.764] = +1.9737 ;
801  (*this)[+0.765] = +1.9803 ;
802  (*this)[+0.766] = +1.9869 ;
803  (*this)[+0.767] = +1.9935 ;
804  (*this)[+0.768] = +2.0001 ;
805  (*this)[+0.769] = +2.0067 ;
806  (*this)[+0.77 ] = +2.0133 ;
807  (*this)[+0.771] = +2.0199 ;
808  (*this)[+0.772] = +2.0265 ;
809  (*this)[+0.773] = +2.0331 ;
810  (*this)[+0.774] = +2.0397 ;
811  (*this)[+0.775] = +2.0469 ;
812  (*this)[+0.776] = +2.0569 ;
813  (*this)[+0.777] = +2.0668 ;
814  (*this)[+0.778] = +2.0768 ;
815  (*this)[+0.779] = +2.0867 ;
816  (*this)[+0.78 ] = +2.0967 ;
817  (*this)[+0.781] = +2.1066 ;
818  (*this)[+0.782] = +2.1166 ;
819  (*this)[+0.783] = +2.1265 ;
820  (*this)[+0.784] = +2.1365 ;
821  (*this)[+0.785] = +2.1464 ;
822  (*this)[+0.786] = +2.1564 ;
823  (*this)[+0.787] = +2.1663 ;
824  (*this)[+0.788] = +2.1763 ;
825  (*this)[+0.789] = +2.1862 ;
826  (*this)[+0.79 ] = +2.1961 ;
827  (*this)[+0.791] = +2.2061 ;
828  (*this)[+0.792] = +2.216 ;
829  (*this)[+0.793] = +2.226 ;
830  (*this)[+0.794] = +2.2359 ;
831  (*this)[+0.795] = +2.2459 ;
832  (*this)[+0.796] = +2.2558 ;
833  (*this)[+0.797] = +2.2658 ;
834  (*this)[+0.798] = +2.2757 ;
835  (*this)[+0.799] = +2.2857 ;
836  (*this)[+0.8 ] = +2.2956 ;
837  (*this)[+0.801] = +2.3056 ;
838  (*this)[+0.802] = +2.3155 ;
839  (*this)[+0.803] = +2.3255 ;
840  (*this)[+0.804] = +2.3354 ;
841  (*this)[+0.805] = +2.3454 ;
842  (*this)[+0.806] = +2.3553 ;
843  (*this)[+0.807] = +2.3653 ;
844  (*this)[+0.808] = +2.3752 ;
845  (*this)[+0.809] = +2.3852 ;
846  (*this)[+0.81 ] = +2.3951 ;
847  (*this)[+0.811] = +2.4051 ;
848  (*this)[+0.812] = +2.415 ;
849  (*this)[+0.813] = +2.425 ;
850  (*this)[+0.814] = +2.4349 ;
851  (*this)[+0.815] = +2.4448 ;
852  (*this)[+0.816] = +2.4548 ;
853  (*this)[+0.817] = +2.4647 ;
854  (*this)[+0.818] = +2.4747 ;
855  (*this)[+0.819] = +2.4846 ;
856  (*this)[+0.82 ] = +2.4946 ;
857  (*this)[+0.821] = +2.5045 ;
858  (*this)[+0.822] = +2.5145 ;
859  (*this)[+0.823] = +2.5244 ;
860  (*this)[+0.824] = +2.5344 ;
861  (*this)[+0.825] = +2.5443 ;
862  (*this)[+0.826] = +2.5543 ;
863  (*this)[+0.827] = +2.5642 ;
864  (*this)[+0.828] = +2.5742 ;
865  (*this)[+0.829] = +2.5841 ;
866  (*this)[+0.83 ] = +2.5941 ;
867  (*this)[+0.831] = +2.604 ;
868  (*this)[+0.832] = +2.614 ;
869  (*this)[+0.833] = +2.6239 ;
870  (*this)[+0.834] = +2.6339 ;
871  (*this)[+0.835] = +2.6438 ;
872  (*this)[+0.836] = +2.6538 ;
873  (*this)[+0.837] = +2.6637 ;
874  (*this)[+0.838] = +2.6736 ;
875  (*this)[+0.839] = +2.6836 ;
876  (*this)[+0.84 ] = +2.6935 ;
877  (*this)[+0.841] = +2.7035 ;
878  (*this)[+0.842] = +2.7134 ;
879  (*this)[+0.843] = +2.7234 ;
880  (*this)[+0.844] = +2.7333 ;
881  (*this)[+0.845] = +2.7433 ;
882  (*this)[+0.846] = +2.7532 ;
883  (*this)[+0.847] = +2.7632 ;
884  (*this)[+0.848] = +2.7731 ;
885  (*this)[+0.849] = +2.7831 ;
886  (*this)[+0.85 ] = +2.793 ;
887  (*this)[+0.851] = +2.803 ;
888  (*this)[+0.852] = +2.8129 ;
889  (*this)[+0.853] = +2.8229 ;
890  (*this)[+0.854] = +2.8328 ;
891  (*this)[+0.855] = +2.8428 ;
892  (*this)[+0.856] = +2.8527 ;
893  (*this)[+0.857] = +2.8627 ;
894  (*this)[+0.858] = +2.8726 ;
895  (*this)[+0.859] = +2.8826 ;
896  (*this)[+0.86 ] = +2.8925 ;
897  (*this)[+0.861] = +2.9025 ;
898  (*this)[+0.862] = +2.9124 ;
899  (*this)[+0.863] = +2.9223 ;
900  (*this)[+0.864] = +2.9323 ;
901  (*this)[+0.865] = +2.9422 ;
902  (*this)[+0.866] = +2.9522 ;
903  (*this)[+0.867] = +2.9621 ;
904  (*this)[+0.868] = +2.9721 ;
905  (*this)[+0.869] = +2.982 ;
906  (*this)[+0.87 ] = +2.992 ;
907  (*this)[+0.871] = +3.0019 ;
908  (*this)[+0.872] = +3.0119 ;
909  (*this)[+0.873] = +3.0218 ;
910  (*this)[+0.874] = +3.0318 ;
911  (*this)[+0.875] = +3.0417 ;
912  (*this)[+0.876] = +3.0574 ;
913  (*this)[+0.877] = +3.076 ;
914  (*this)[+0.878] = +3.0947 ;
915  (*this)[+0.879] = +3.1133 ;
916  (*this)[+0.88 ] = +3.1319 ;
917  (*this)[+0.881] = +3.1506 ;
918  (*this)[+0.882] = +3.1692 ;
919  (*this)[+0.883] = +3.1879 ;
920  (*this)[+0.884] = +3.2065 ;
921  (*this)[+0.885] = +3.2251 ;
922  (*this)[+0.886] = +3.2438 ;
923  (*this)[+0.887] = +3.2624 ;
924  (*this)[+0.888] = +3.2811 ;
925  (*this)[+0.889] = +3.2997 ;
926  (*this)[+0.89 ] = +3.3183 ;
927  (*this)[+0.891] = +3.337 ;
928  (*this)[+0.892] = +3.3556 ;
929  (*this)[+0.893] = +3.3742 ;
930  (*this)[+0.894] = +3.3929 ;
931  (*this)[+0.895] = +3.4115 ;
932  (*this)[+0.896] = +3.4302 ;
933  (*this)[+0.897] = +3.4488 ;
934  (*this)[+0.898] = +3.4674 ;
935  (*this)[+0.899] = +3.4861 ;
936  (*this)[+0.9 ] = +3.5047 ;
937  (*this)[+0.901] = +3.5234 ;
938  (*this)[+0.902] = +3.542 ;
939  (*this)[+0.903] = +3.5606 ;
940  (*this)[+0.904] = +3.5793 ;
941  (*this)[+0.905] = +3.5979 ;
942  (*this)[+0.906] = +3.6165 ;
943  (*this)[+0.907] = +3.6352 ;
944  (*this)[+0.908] = +3.6538 ;
945  (*this)[+0.909] = +3.6725 ;
946  (*this)[+0.91 ] = +3.6911 ;
947  (*this)[+0.911] = +3.7097 ;
948  (*this)[+0.912] = +3.7284 ;
949  (*this)[+0.913] = +3.747 ;
950  (*this)[+0.914] = +3.7657 ;
951  (*this)[+0.915] = +3.7843 ;
952  (*this)[+0.916] = +3.8029 ;
953  (*this)[+0.917] = +3.8216 ;
954  (*this)[+0.918] = +3.8402 ;
955  (*this)[+0.919] = +3.8588 ;
956  (*this)[+0.92 ] = +3.8775 ;
957  (*this)[+0.921] = +3.8961 ;
958  (*this)[+0.922] = +3.9148 ;
959  (*this)[+0.923] = +3.9334 ;
960  (*this)[+0.924] = +3.952 ;
961  (*this)[+0.925] = +3.9707 ;
962  (*this)[+0.926] = +3.9893 ;
963  (*this)[+0.927] = +4.008 ;
964  (*this)[+0.928] = +4.0266 ;
965  (*this)[+0.929] = +4.0454 ;
966  (*this)[+0.93 ] = +4.0887 ;
967  (*this)[+0.931] = +4.1321 ;
968  (*this)[+0.932] = +4.1755 ;
969  (*this)[+0.933] = +4.2188 ;
970  (*this)[+0.934] = +4.2622 ;
971  (*this)[+0.935] = +4.3055 ;
972  (*this)[+0.936] = +4.3489 ;
973  (*this)[+0.937] = +4.3923 ;
974  (*this)[+0.938] = +4.4356 ;
975  (*this)[+0.939] = +4.479 ;
976  (*this)[+0.94 ] = +4.5224 ;
977  (*this)[+0.941] = +4.5657 ;
978  (*this)[+0.942] = +4.6091 ;
979  (*this)[+0.943] = +4.6524 ;
980  (*this)[+0.944] = +4.6958 ;
981  (*this)[+0.945] = +4.7392 ;
982  (*this)[+0.946] = +4.7825 ;
983  (*this)[+0.947] = +4.8259 ;
984  (*this)[+0.948] = +4.8692 ;
985  (*this)[+0.949] = +4.9126 ;
986  (*this)[+0.95 ] = +4.956 ;
987  (*this)[+0.951] = +4.9993 ;
988  (*this)[+0.952] = +5.0427 ;
989  (*this)[+0.953] = +5.1611 ;
990  (*this)[+0.954] = +5.284 ;
991  (*this)[+0.955] = +5.4068 ;
992  (*this)[+0.956] = +5.5297 ;
993  (*this)[+0.957] = +5.6525 ;
994  (*this)[+0.958] = +5.7754 ;
995  (*this)[+0.959] = +5.8982 ;
996  (*this)[+0.96 ] = +6.0211 ;
997  (*this)[+0.961] = +6.2994 ;
998  (*this)[+0.962] = +6.6154 ;
999  (*this)[+0.963] = +6.9315 ;
1000  (*this)[+0.964] = +7.3775 ;
1001  (*this)[+0.965] = +7.8967 ;
1002  (*this)[+0.966] = +8.5201 ;
1003  (*this)[+0.967] = +9.2166 ;
1004  (*this)[+0.968] = +10.03 ;
1005  (*this)[+0.969] = +11.01 ;
1006  (*this)[+0.97 ] = +12.125 ;
1007  (*this)[+0.971] = +13.269 ;
1008  (*this)[+0.972] = +14.346 ;
1009  (*this)[+0.973] = +15.338 ;
1010  (*this)[+0.974] = +16.253 ;
1011  (*this)[+0.975] = +17.1 ;
1012  (*this)[+0.976] = +17.888 ;
1013  (*this)[+0.977] = +18.645 ;
1014  (*this)[+0.978] = +19.383 ;
1015  (*this)[+0.979] = +20.105 ;
1016  (*this)[+0.98 ] = +20.81 ;
1017  (*this)[+0.981] = +21.506 ;
1018  (*this)[+0.982] = +22.197 ;
1019  (*this)[+0.983] = +22.88 ;
1020  (*this)[+0.984] = +23.548 ;
1021  (*this)[+0.985] = +24.205 ;
1022  (*this)[+0.986] = +24.841 ;
1023  (*this)[+0.987] = +25.461 ;
1024  (*this)[+0.988] = +26.072 ;
1025  (*this)[+0.989] = +26.673 ;
1026  (*this)[+0.99 ] = +27.277 ;
1027  (*this)[+0.991] = +27.886 ;
1028  (*this)[+0.992] = +28.521 ;
1029  (*this)[+0.993] = +29.178 ;
1030  (*this)[+0.994] = +29.899 ;
1031  (*this)[+0.995] = +30.726 ;
1032  (*this)[+0.996] = +31.705 ;
1033  (*this)[+0.997] = +32.956 ;
1034  (*this)[+0.998] = +34.933 ;
1035  (*this)[+0.999] = +40.25 ;
1036  (*this)[+1.000] = +62.276 ;
1037 
1038  compile();
1039  };
1040 
1041 
1042  /**
1043  * Get integrated probablity from \f$-\infty\f$ to given time.
1044  *
1045  * \param t_ns time [ns]
1046  * \return probability
1047  */
1048  inline double getProbability(const double t_ns) const
1049  {
1050  using namespace std;
1051  using namespace JPP;
1052 
1053  const_iterator p = std::lower_bound(this->begin(),
1054  this->end(),
1055  t_ns,
1056  make_comparator(&value_type::getY, JComparison::lt()));
1057 
1058  if (p == this->end())
1059  return 1.0;
1060  else if (p == this->begin())
1061  return 0.0;
1062  else
1063  return p->getX();
1064  }
1065  };
1066 
1067 
1068  /**
1069  * Function object to generate transition time.
1070  */
1072 }
1073 
1074 #endif
JComparator< JResult_t T::*, JComparison::lt > make_comparator(JResult_t T::*member)
Helper method to create comparator between values of data member.
Definition: JComparator.hh:143
Type definition of a spline interpolation based on a JGridCollection with result type double...
static const JTransitionTimeGenerator getTransitionTime
Function object to generate transition time.
double getProbability(const double t_ns) const
Get integrated probablity from to given time.
Auxiliary class to generate PMT transition times based on measurements by Oleg Kalekin et al...