Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JPMTTransitTimeProbability.hh
Go to the documentation of this file.
1 #ifndef __JDETECTOR__JPMTTRANSITTIMEPROBABILITY__
2 #define __JDETECTOR__JPMTTRANSITTIMEPROBABILITY__
3 
4 #include <algorithm>
5 
6 #include "JTools/JFunction1D_t.hh"
7 
8 
9 /**
10  * \author mdejong
11  */
12 
13 namespace JDETECTOR {}
14 namespace JPP { using namespace JDETECTOR; }
15 
16 namespace JDETECTOR {
17 
18  /**
19  * Probability distibution of the PMT transition times based on measurements
20  * by Oleg Kalekin et al. (see also examples/JDetector/JTTS.cc).
21  *
22  * For a given transition time [ns], this function returns the correspoding probability.
23  */
24  class JTransitionTimeProbability :
25  public JTOOLS::JGridSplineFunction1D_t
26  {
27  public:
28  /**
29  * Default constructor.
30  */
31  JTransitionTimeProbability() :
32  JTOOLS::JGridSplineFunction1D_t()
33  {
34  (*this)[-11.4549] = +0.00018;
35  (*this)[-11.3869] = +0.00018;
36  (*this)[-11.3189] = +0.00018;
37  (*this)[-11.2509] = +0.00018;
38  (*this)[-11.1829] = +0.00018;
39  (*this)[-11.1149] = +0.00018;
40  (*this)[-11.0469] = +0.00018;
41  (*this)[-10.9789] = +0.00018;
42  (*this)[-10.9109] = +0.00017;
43  (*this)[-10.8429] = +0.00017;
44  (*this)[-10.7749] = +0.00017;
45  (*this)[-10.7069] = +0.00017;
46  (*this)[-10.6389] = +0.00017;
47  (*this)[-10.5709] = +0.00017;
48  (*this)[-10.5029] = +0.00017;
49  (*this)[-10.4349] = +0.00017;
50  (*this)[-10.3669] = +0.00017;
51  (*this)[-10.2989] = +0.00017;
52  (*this)[-10.2309] = +0.00017;
53  (*this)[-10.1629] = +0.00017;
54  (*this)[-10.0949] = +0.00017;
55  (*this)[-10.0269] = +0.00017;
56  (*this)[ -9.9589] = +0.00017;
57  (*this)[ -9.8909] = +0.00017;
58  (*this)[ -9.8229] = +0.00017;
59  (*this)[ -9.7549] = +0.00017;
60  (*this)[ -9.6869] = +0.00017;
61  (*this)[ -9.6189] = +0.00017;
62  (*this)[ -9.5509] = +0.00017;
63  (*this)[ -9.4829] = +0.00017;
64  (*this)[ -9.4149] = +0.00017;
65  (*this)[ -9.3469] = +0.00017;
66  (*this)[ -9.2789] = +0.00018;
67  (*this)[ -9.2109] = +0.00018;
68  (*this)[ -9.1429] = +0.00018;
69  (*this)[ -9.0749] = +0.00018;
70  (*this)[ -9.0069] = +0.00018;
71  (*this)[ -8.9389] = +0.00018;
72  (*this)[ -8.8709] = +0.00018;
73  (*this)[ -8.8029] = +0.00018;
74  (*this)[ -8.7349] = +0.00018;
75  (*this)[ -8.6669] = +0.00018;
76  (*this)[ -8.5989] = +0.00018;
77  (*this)[ -8.5309] = +0.00018;
78  (*this)[ -8.4629] = +0.00018;
79  (*this)[ -8.3949] = +0.00019;
80  (*this)[ -8.3269] = +0.00019;
81  (*this)[ -8.2589] = +0.00020;
82  (*this)[ -8.1909] = +0.00020;
83  (*this)[ -8.1229] = +0.00020;
84  (*this)[ -8.0549] = +0.00021;
85  (*this)[ -7.9869] = +0.00021;
86  (*this)[ -7.9189] = +0.00022;
87  (*this)[ -7.8509] = +0.00022;
88  (*this)[ -7.7829] = +0.00022;
89  (*this)[ -7.7149] = +0.00023;
90  (*this)[ -7.6469] = +0.00023;
91  (*this)[ -7.5789] = +0.00023;
92  (*this)[ -7.5109] = +0.00024;
93  (*this)[ -7.4429] = +0.00025;
94  (*this)[ -7.3749] = +0.00026;
95  (*this)[ -7.3069] = +0.00028;
96  (*this)[ -7.2389] = +0.00030;
97  (*this)[ -7.1709] = +0.00031;
98  (*this)[ -7.1029] = +0.00033;
99  (*this)[ -7.0349] = +0.00035;
100  (*this)[ -6.9669] = +0.00037;
101  (*this)[ -6.8989] = +0.00038;
102  (*this)[ -6.8309] = +0.00040;
103  (*this)[ -6.7629] = +0.00042;
104  (*this)[ -6.6949] = +0.00044;
105  (*this)[ -6.6269] = +0.00045;
106  (*this)[ -6.5589] = +0.00047;
107  (*this)[ -6.4909] = +0.00049;
108  (*this)[ -6.4229] = +0.00054;
109  (*this)[ -6.3549] = +0.00063;
110  (*this)[ -6.2869] = +0.00073;
111  (*this)[ -6.2189] = +0.00082;
112  (*this)[ -6.1509] = +0.00091;
113  (*this)[ -6.0829] = +0.00101;
114  (*this)[ -6.0149] = +0.00110;
115  (*this)[ -5.9469] = +0.00119;
116  (*this)[ -5.8789] = +0.00128;
117  (*this)[ -5.8109] = +0.00138;
118  (*this)[ -5.7429] = +0.00147;
119  (*this)[ -5.6749] = +0.00156;
120  (*this)[ -5.6069] = +0.00166;
121  (*this)[ -5.5389] = +0.00175;
122  (*this)[ -5.4709] = +0.00184;
123  (*this)[ -5.4029] = +0.00219;
124  (*this)[ -5.3349] = +0.00261;
125  (*this)[ -5.2669] = +0.00304;
126  (*this)[ -5.1989] = +0.00346;
127  (*this)[ -5.1309] = +0.00388;
128  (*this)[ -5.0629] = +0.00431;
129  (*this)[ -4.9949] = +0.00473;
130  (*this)[ -4.9269] = +0.00516;
131  (*this)[ -4.8589] = +0.00558;
132  (*this)[ -4.7909] = +0.00600;
133  (*this)[ -4.7229] = +0.00643;
134  (*this)[ -4.6549] = +0.00685;
135  (*this)[ -4.5869] = +0.00728;
136  (*this)[ -4.5189] = +0.00770;
137  (*this)[ -4.4509] = +0.00818;
138  (*this)[ -4.3829] = +0.00959;
139  (*this)[ -4.3149] = +0.01100;
140  (*this)[ -4.2469] = +0.01241;
141  (*this)[ -4.1789] = +0.01382;
142  (*this)[ -4.1109] = +0.01523;
143  (*this)[ -4.0429] = +0.01664;
144  (*this)[ -3.9749] = +0.01805;
145  (*this)[ -3.9069] = +0.01946;
146  (*this)[ -3.8389] = +0.02087;
147  (*this)[ -3.7709] = +0.02228;
148  (*this)[ -3.7029] = +0.02369;
149  (*this)[ -3.6349] = +0.02510;
150  (*this)[ -3.5669] = +0.02651;
151  (*this)[ -3.4989] = +0.02792;
152  (*this)[ -3.4309] = +0.02993;
153  (*this)[ -3.3629] = +0.03302;
154  (*this)[ -3.2949] = +0.03611;
155  (*this)[ -3.2269] = +0.03920;
156  (*this)[ -3.1589] = +0.04230;
157  (*this)[ -3.0909] = +0.04539;
158  (*this)[ -3.0229] = +0.04848;
159  (*this)[ -2.9549] = +0.05157;
160  (*this)[ -2.8869] = +0.05467;
161  (*this)[ -2.8189] = +0.05776;
162  (*this)[ -2.7509] = +0.06085;
163  (*this)[ -2.6829] = +0.06394;
164  (*this)[ -2.6149] = +0.06704;
165  (*this)[ -2.5469] = +0.07013;
166  (*this)[ -2.4789] = +0.07322;
167  (*this)[ -2.4109] = +0.07709;
168  (*this)[ -2.3429] = +0.08138;
169  (*this)[ -2.2749] = +0.08567;
170  (*this)[ -2.2069] = +0.08997;
171  (*this)[ -2.1389] = +0.09426;
172  (*this)[ -2.0709] = +0.09855;
173  (*this)[ -2.0029] = +0.10284;
174  (*this)[ -1.9349] = +0.10713;
175  (*this)[ -1.8669] = +0.11143;
176  (*this)[ -1.7989] = +0.11572;
177  (*this)[ -1.7309] = +0.12001;
178  (*this)[ -1.6629] = +0.12430;
179  (*this)[ -1.5949] = +0.12860;
180  (*this)[ -1.5269] = +0.13289;
181  (*this)[ -1.4589] = +0.13718;
182  (*this)[ -1.3909] = +0.14044;
183  (*this)[ -1.3229] = +0.14364;
184  (*this)[ -1.2549] = +0.14684;
185  (*this)[ -1.1869] = +0.15003;
186  (*this)[ -1.1189] = +0.15323;
187  (*this)[ -1.0509] = +0.15643;
188  (*this)[ -0.9829] = +0.15962;
189  (*this)[ -0.9149] = +0.16282;
190  (*this)[ -0.8469] = +0.16602;
191  (*this)[ -0.7789] = +0.16921;
192  (*this)[ -0.7109] = +0.17241;
193  (*this)[ -0.6429] = +0.17561;
194  (*this)[ -0.5749] = +0.17880;
195  (*this)[ -0.5069] = +0.18200;
196  (*this)[ -0.4389] = +0.18449;
197  (*this)[ -0.3709] = +0.18467;
198  (*this)[ -0.3029] = +0.18485;
199  (*this)[ -0.2349] = +0.18503;
200  (*this)[ -0.1669] = +0.18521;
201  (*this)[ -0.0989] = +0.18539;
202  (*this)[ -0.0309] = +0.18557;
203  (*this)[ +0.0371] = +0.18575;
204  (*this)[ +0.1051] = +0.18593;
205  (*this)[ +0.1731] = +0.18611;
206  (*this)[ +0.2411] = +0.18629;
207  (*this)[ +0.3091] = +0.18647;
208  (*this)[ +0.3771] = +0.18665;
209  (*this)[ +0.4451] = +0.18683;
210  (*this)[ +0.5131] = +0.18701;
211  (*this)[ +0.5811] = +0.18581;
212  (*this)[ +0.6491] = +0.18339;
213  (*this)[ +0.7171] = +0.18097;
214  (*this)[ +0.7851] = +0.17854;
215  (*this)[ +0.8531] = +0.17612;
216  (*this)[ +0.9211] = +0.17370;
217  (*this)[ +0.9891] = +0.17127;
218  (*this)[ +1.0571] = +0.16885;
219  (*this)[ +1.1251] = +0.16643;
220  (*this)[ +1.1931] = +0.16400;
221  (*this)[ +1.2611] = +0.16158;
222  (*this)[ +1.3291] = +0.15916;
223  (*this)[ +1.3971] = +0.15674;
224  (*this)[ +1.4651] = +0.15431;
225  (*this)[ +1.5331] = +0.15189;
226  (*this)[ +1.6011] = +0.14861;
227  (*this)[ +1.6691] = +0.14514;
228  (*this)[ +1.7371] = +0.14168;
229  (*this)[ +1.8051] = +0.13822;
230  (*this)[ +1.8731] = +0.13475;
231  (*this)[ +1.9411] = +0.13129;
232  (*this)[ +2.0091] = +0.12783;
233  (*this)[ +2.0771] = +0.12436;
234  (*this)[ +2.1451] = +0.12090;
235  (*this)[ +2.2131] = +0.11743;
236  (*this)[ +2.2811] = +0.11397;
237  (*this)[ +2.3491] = +0.11051;
238  (*this)[ +2.4171] = +0.10704;
239  (*this)[ +2.4851] = +0.10358;
240  (*this)[ +2.5531] = +0.10015;
241  (*this)[ +2.6211] = +0.09696;
242  (*this)[ +2.6891] = +0.09377;
243  (*this)[ +2.7571] = +0.09059;
244  (*this)[ +2.8251] = +0.08740;
245  (*this)[ +2.8931] = +0.08421;
246  (*this)[ +2.9611] = +0.08103;
247  (*this)[ +3.0291] = +0.07784;
248  (*this)[ +3.0971] = +0.07465;
249  (*this)[ +3.1651] = +0.07146;
250  (*this)[ +3.2331] = +0.06828;
251  (*this)[ +3.3011] = +0.06509;
252  (*this)[ +3.3691] = +0.06190;
253  (*this)[ +3.4371] = +0.05871;
254  (*this)[ +3.5051] = +0.05553;
255  (*this)[ +3.5731] = +0.05280;
256  (*this)[ +3.6411] = +0.05072;
257  (*this)[ +3.7091] = +0.04864;
258  (*this)[ +3.7771] = +0.04656;
259  (*this)[ +3.8451] = +0.04448;
260  (*this)[ +3.9131] = +0.04240;
261  (*this)[ +3.9811] = +0.04032;
262  (*this)[ +4.0491] = +0.03823;
263  (*this)[ +4.1171] = +0.03615;
264  (*this)[ +4.1851] = +0.03407;
265  (*this)[ +4.2531] = +0.03199;
266  (*this)[ +4.3211] = +0.02991;
267  (*this)[ +4.3891] = +0.02783;
268  (*this)[ +4.4571] = +0.02575;
269  (*this)[ +4.5251] = +0.02367;
270  (*this)[ +4.5931] = +0.02235;
271  (*this)[ +4.6611] = +0.02133;
272  (*this)[ +4.7291] = +0.02032;
273  (*this)[ +4.7971] = +0.01930;
274  (*this)[ +4.8651] = +0.01829;
275  (*this)[ +4.9331] = +0.01727;
276  (*this)[ +5.0011] = +0.01626;
277  (*this)[ +5.0691] = +0.01524;
278  (*this)[ +5.1371] = +0.01423;
279  (*this)[ +5.2051] = +0.01321;
280  (*this)[ +5.2731] = +0.01220;
281  (*this)[ +5.3411] = +0.01118;
282  (*this)[ +5.4091] = +0.01017;
283  (*this)[ +5.4771] = +0.00915;
284  (*this)[ +5.5451] = +0.00814;
285  (*this)[ +5.6131] = +0.00780;
286  (*this)[ +5.6811] = +0.00746;
287  (*this)[ +5.7491] = +0.00712;
288  (*this)[ +5.8171] = +0.00679;
289  (*this)[ +5.8851] = +0.00645;
290  (*this)[ +5.9531] = +0.00611;
291  (*this)[ +6.0211] = +0.00577;
292  (*this)[ +6.0891] = +0.00543;
293  (*this)[ +6.1571] = +0.00509;
294  (*this)[ +6.2251] = +0.00476;
295  (*this)[ +6.2931] = +0.00442;
296  (*this)[ +6.3611] = +0.00408;
297  (*this)[ +6.4291] = +0.00374;
298  (*this)[ +6.4971] = +0.00340;
299  (*this)[ +6.5651] = +0.00314;
300  (*this)[ +6.6331] = +0.00306;
301  (*this)[ +6.7011] = +0.00297;
302  (*this)[ +6.7691] = +0.00289;
303  (*this)[ +6.8371] = +0.00280;
304  (*this)[ +6.9051] = +0.00272;
305  (*this)[ +6.9731] = +0.00263;
306  (*this)[ +7.0411] = +0.00255;
307  (*this)[ +7.1091] = +0.00247;
308  (*this)[ +7.1771] = +0.00238;
309  (*this)[ +7.2451] = +0.00230;
310  (*this)[ +7.3131] = +0.00221;
311  (*this)[ +7.3811] = +0.00213;
312  (*this)[ +7.4491] = +0.00205;
313  (*this)[ +7.5171] = +0.00196;
314  (*this)[ +7.5851] = +0.00191;
315  (*this)[ +7.6531] = +0.00188;
316  (*this)[ +7.7211] = +0.00185;
317  (*this)[ +7.7891] = +0.00182;
318  (*this)[ +7.8571] = +0.00179;
319  (*this)[ +7.9251] = +0.00177;
320  (*this)[ +7.9931] = +0.00174;
321  (*this)[ +8.0611] = +0.00171;
322  (*this)[ +8.1291] = +0.00168;
323  (*this)[ +8.1971] = +0.00165;
324  (*this)[ +8.2651] = +0.00162;
325  (*this)[ +8.3331] = +0.00159;
326  (*this)[ +8.4011] = +0.00156;
327  (*this)[ +8.4691] = +0.00154;
328  (*this)[ +8.5371] = +0.00151;
329  (*this)[ +8.6051] = +0.00149;
330  (*this)[ +8.6731] = +0.00147;
331  (*this)[ +8.7411] = +0.00145;
332  (*this)[ +8.8091] = +0.00143;
333  (*this)[ +8.8771] = +0.00141;
334  (*this)[ +8.9451] = +0.00139;
335  (*this)[ +9.0131] = +0.00138;
336  (*this)[ +9.0811] = +0.00136;
337  (*this)[ +9.1491] = +0.00134;
338  (*this)[ +9.2171] = +0.00132;
339  (*this)[ +9.2851] = +0.00130;
340  (*this)[ +9.3531] = +0.00128;
341  (*this)[ +9.4211] = +0.00126;
342  (*this)[ +9.4891] = +0.00124;
343  (*this)[ +9.5571] = +0.00123;
344  (*this)[ +9.6251] = +0.00121;
345  (*this)[ +9.6931] = +0.00120;
346  (*this)[ +9.7611] = +0.00118;
347  (*this)[ +9.8291] = +0.00117;
348  (*this)[ +9.8971] = +0.00115;
349  (*this)[ +9.9651] = +0.00114;
350  (*this)[+10.0331] = +0.00113;
351  (*this)[+10.1011] = +0.00111;
352  (*this)[+10.1691] = +0.00110;
353  (*this)[+10.2371] = +0.00108;
354  (*this)[+10.3051] = +0.00107;
355  (*this)[+10.3731] = +0.00105;
356  (*this)[+10.4411] = +0.00104;
357  (*this)[+10.5091] = +0.00102;
358  (*this)[+10.5771] = +0.00101;
359  (*this)[+10.6451] = +0.00101;
360  (*this)[+10.7131] = +0.00100;
361  (*this)[+10.7811] = +0.00099;
362  (*this)[+10.8491] = +0.00098;
363  (*this)[+10.9171] = +0.00097;
364  (*this)[+10.9851] = +0.00096;
365  (*this)[+11.0531] = +0.00095;
366  (*this)[+11.1211] = +0.00095;
367  (*this)[+11.1891] = +0.00094;
368  (*this)[+11.2571] = +0.00093;
369  (*this)[+11.3251] = +0.00092;
370  (*this)[+11.3931] = +0.00091;
371  (*this)[+11.4611] = +0.00090;
372  (*this)[+11.5291] = +0.00090;
373  (*this)[+11.5971] = +0.00089;
374  (*this)[+11.6651] = +0.00089;
375  (*this)[+11.7331] = +0.00089;
376  (*this)[+11.8011] = +0.00089;
377  (*this)[+11.8691] = +0.00089;
378  (*this)[+11.9371] = +0.00088;
379  (*this)[+12.0051] = +0.00088;
380  (*this)[+12.0731] = +0.00088;
381  (*this)[+12.1411] = +0.00088;
382  (*this)[+12.2091] = +0.00088;
383  (*this)[+12.2771] = +0.00087;
384  (*this)[+12.3451] = +0.00087;
385  (*this)[+12.4131] = +0.00087;
386  (*this)[+12.4811] = +0.00087;
387  (*this)[+12.5491] = +0.00087;
388  (*this)[+12.6171] = +0.00087;
389  (*this)[+12.6851] = +0.00087;
390  (*this)[+12.7531] = +0.00087;
391  (*this)[+12.8211] = +0.00088;
392  (*this)[+12.8891] = +0.00088;
393  (*this)[+12.9571] = +0.00088;
394  (*this)[+13.0251] = +0.00089;
395  (*this)[+13.0931] = +0.00089;
396  (*this)[+13.1611] = +0.00089;
397  (*this)[+13.2291] = +0.00090;
398  (*this)[+13.2971] = +0.00090;
399  (*this)[+13.3651] = +0.00090;
400  (*this)[+13.4331] = +0.00090;
401  (*this)[+13.5011] = +0.00091;
402  (*this)[+13.5691] = +0.00091;
403  (*this)[+13.6371] = +0.00092;
404  (*this)[+13.7051] = +0.00092;
405  (*this)[+13.7731] = +0.00093;
406  (*this)[+13.8411] = +0.00093;
407  (*this)[+13.9091] = +0.00094;
408  (*this)[+13.9771] = +0.00094;
409  (*this)[+14.0451] = +0.00095;
410  (*this)[+14.1131] = +0.00095;
411  (*this)[+14.1811] = +0.00096;
412  (*this)[+14.2491] = +0.00096;
413  (*this)[+14.3171] = +0.00097;
414  (*this)[+14.3851] = +0.00097;
415  (*this)[+14.4531] = +0.00098;
416  (*this)[+14.5211] = +0.00098;
417  (*this)[+14.5891] = +0.00099;
418  (*this)[+14.6571] = +0.00099;
419  (*this)[+14.7251] = +0.00100;
420  (*this)[+14.7931] = +0.00100;
421  (*this)[+14.8611] = +0.00101;
422  (*this)[+14.9291] = +0.00102;
423  (*this)[+14.9971] = +0.00102;
424  (*this)[+15.0651] = +0.00103;
425  (*this)[+15.1331] = +0.00103;
426  (*this)[+15.2011] = +0.00104;
427  (*this)[+15.2691] = +0.00105;
428  (*this)[+15.3371] = +0.00105;
429  (*this)[+15.4051] = +0.00106;
430  (*this)[+15.4731] = +0.00106;
431  (*this)[+15.5411] = +0.00107;
432  (*this)[+15.6091] = +0.00108;
433  (*this)[+15.6771] = +0.00108;
434  (*this)[+15.7451] = +0.00109;
435  (*this)[+15.8131] = +0.00110;
436  (*this)[+15.8811] = +0.00110;
437  (*this)[+15.9491] = +0.00111;
438  (*this)[+16.0171] = +0.00112;
439  (*this)[+16.0851] = +0.00113;
440  (*this)[+16.1531] = +0.00113;
441  (*this)[+16.2211] = +0.00114;
442  (*this)[+16.2891] = +0.00115;
443  (*this)[+16.3571] = +0.00115;
444  (*this)[+16.4251] = +0.00116;
445  (*this)[+16.4931] = +0.00117;
446  (*this)[+16.5611] = +0.00118;
447  (*this)[+16.6291] = +0.00118;
448  (*this)[+16.6971] = +0.00119;
449  (*this)[+16.7651] = +0.00119;
450  (*this)[+16.8331] = +0.00120;
451  (*this)[+16.9011] = +0.00121;
452  (*this)[+16.9691] = +0.00121;
453  (*this)[+17.0371] = +0.00122;
454  (*this)[+17.1051] = +0.00123;
455  (*this)[+17.1731] = +0.00123;
456  (*this)[+17.2411] = +0.00124;
457  (*this)[+17.3091] = +0.00125;
458  (*this)[+17.3771] = +0.00125;
459  (*this)[+17.4451] = +0.00126;
460  (*this)[+17.5131] = +0.00127;
461  (*this)[+17.5811] = +0.00127;
462  (*this)[+17.6491] = +0.00128;
463  (*this)[+17.7171] = +0.00128;
464  (*this)[+17.7851] = +0.00128;
465  (*this)[+17.8531] = +0.00129;
466  (*this)[+17.9211] = +0.00129;
467  (*this)[+17.9891] = +0.00130;
468  (*this)[+18.0571] = +0.00130;
469  (*this)[+18.1251] = +0.00131;
470  (*this)[+18.1931] = +0.00131;
471  (*this)[+18.2611] = +0.00132;
472  (*this)[+18.3291] = +0.00132;
473  (*this)[+18.3971] = +0.00133;
474  (*this)[+18.4651] = +0.00133;
475  (*this)[+18.5331] = +0.00133;
476  (*this)[+18.6011] = +0.00134;
477  (*this)[+18.6691] = +0.00134;
478  (*this)[+18.7371] = +0.00134;
479  (*this)[+18.8051] = +0.00135;
480  (*this)[+18.8731] = +0.00135;
481  (*this)[+18.9411] = +0.00135;
482  (*this)[+19.0091] = +0.00136;
483  (*this)[+19.0771] = +0.00136;
484  (*this)[+19.1451] = +0.00136;
485  (*this)[+19.2131] = +0.00137;
486  (*this)[+19.2811] = +0.00137;
487  (*this)[+19.3491] = +0.00137;
488  (*this)[+19.4171] = +0.00137;
489  (*this)[+19.4851] = +0.00138;
490  (*this)[+19.5531] = +0.00138;
491  (*this)[+19.6211] = +0.00138;
492  (*this)[+19.6891] = +0.00139;
493  (*this)[+19.7571] = +0.00139;
494  (*this)[+19.8251] = +0.00139;
495  (*this)[+19.8931] = +0.00139;
496  (*this)[+19.9611] = +0.00140;
497  (*this)[+20.0291] = +0.00140;
498  (*this)[+20.0971] = +0.00140;
499  (*this)[+20.1651] = +0.00140;
500  (*this)[+20.2331] = +0.00141;
501  (*this)[+20.3011] = +0.00141;
502  (*this)[+20.3691] = +0.00141;
503  (*this)[+20.4371] = +0.00142;
504  (*this)[+20.5051] = +0.00142;
505  (*this)[+20.5731] = +0.00142;
506  (*this)[+20.6411] = +0.00142;
507  (*this)[+20.7091] = +0.00142;
508  (*this)[+20.7771] = +0.00143;
509  (*this)[+20.8451] = +0.00143;
510  (*this)[+20.9131] = +0.00143;
511  (*this)[+20.9811] = +0.00143;
512  (*this)[+21.0491] = +0.00143;
513  (*this)[+21.1171] = +0.00143;
514  (*this)[+21.1851] = +0.00143;
515  (*this)[+21.2531] = +0.00144;
516  (*this)[+21.3211] = +0.00144;
517  (*this)[+21.3891] = +0.00144;
518  (*this)[+21.4571] = +0.00144;
519  (*this)[+21.5251] = +0.00144;
520  (*this)[+21.5931] = +0.00144;
521  (*this)[+21.6611] = +0.00145;
522  (*this)[+21.7291] = +0.00145;
523  (*this)[+21.7971] = +0.00145;
524  (*this)[+21.8651] = +0.00145;
525  (*this)[+21.9331] = +0.00145;
526  (*this)[+22.0011] = +0.00145;
527  (*this)[+22.0691] = +0.00145;
528  (*this)[+22.1371] = +0.00146;
529  (*this)[+22.2051] = +0.00146;
530  (*this)[+22.2731] = +0.00146;
531  (*this)[+22.3411] = +0.00146;
532  (*this)[+22.4091] = +0.00146;
533  (*this)[+22.4771] = +0.00146;
534  (*this)[+22.5451] = +0.00146;
535  (*this)[+22.6131] = +0.00147;
536  (*this)[+22.6811] = +0.00147;
537  (*this)[+22.7491] = +0.00147;
538  (*this)[+22.8171] = +0.00148;
539  (*this)[+22.8851] = +0.00148;
540  (*this)[+22.9531] = +0.00148;
541  (*this)[+23.0211] = +0.00148;
542  (*this)[+23.0891] = +0.00149;
543  (*this)[+23.1571] = +0.00149;
544  (*this)[+23.2251] = +0.00149;
545  (*this)[+23.2931] = +0.00150;
546  (*this)[+23.3611] = +0.00150;
547  (*this)[+23.4291] = +0.00150;
548  (*this)[+23.4971] = +0.00151;
549  (*this)[+23.5651] = +0.00151;
550  (*this)[+23.6331] = +0.00151;
551  (*this)[+23.7011] = +0.00152;
552  (*this)[+23.7691] = +0.00152;
553  (*this)[+23.8371] = +0.00153;
554  (*this)[+23.9051] = +0.00153;
555  (*this)[+23.9731] = +0.00153;
556  (*this)[+24.0411] = +0.00154;
557  (*this)[+24.1091] = +0.00154;
558  (*this)[+24.1771] = +0.00155;
559  (*this)[+24.2451] = +0.00155;
560  (*this)[+24.3131] = +0.00156;
561  (*this)[+24.3811] = +0.00156;
562  (*this)[+24.4491] = +0.00156;
563  (*this)[+24.5171] = +0.00157;
564  (*this)[+24.5851] = +0.00157;
565  (*this)[+24.6531] = +0.00158;
566  (*this)[+24.7211] = +0.00158;
567  (*this)[+24.7891] = +0.00159;
568  (*this)[+24.8571] = +0.00159;
569  (*this)[+24.9251] = +0.00160;
570  (*this)[+24.9931] = +0.00160;
571  (*this)[+25.0611] = +0.00160;
572  (*this)[+25.1291] = +0.00161;
573  (*this)[+25.1971] = +0.00161;
574  (*this)[+25.2651] = +0.00162;
575  (*this)[+25.3331] = +0.00162;
576  (*this)[+25.4011] = +0.00163;
577  (*this)[+25.4691] = +0.00163;
578  (*this)[+25.5371] = +0.00164;
579  (*this)[+25.6051] = +0.00164;
580  (*this)[+25.6731] = +0.00164;
581  (*this)[+25.7411] = +0.00164;
582  (*this)[+25.8091] = +0.00164;
583  (*this)[+25.8771] = +0.00165;
584  (*this)[+25.9451] = +0.00165;
585  (*this)[+26.0131] = +0.00165;
586  (*this)[+26.0811] = +0.00165;
587  (*this)[+26.1491] = +0.00165;
588  (*this)[+26.2171] = +0.00165;
589  (*this)[+26.2851] = +0.00166;
590  (*this)[+26.3531] = +0.00166;
591  (*this)[+26.4211] = +0.00166;
592  (*this)[+26.4891] = +0.00166;
593  (*this)[+26.5571] = +0.00166;
594  (*this)[+26.6251] = +0.00166;
595  (*this)[+26.6931] = +0.00166;
596  (*this)[+26.7611] = +0.00166;
597  (*this)[+26.8291] = +0.00166;
598  (*this)[+26.8971] = +0.00166;
599  (*this)[+26.9651] = +0.00165;
600  (*this)[+27.0331] = +0.00165;
601  (*this)[+27.1011] = +0.00165;
602  (*this)[+27.1691] = +0.00165;
603  (*this)[+27.2371] = +0.00165;
604  (*this)[+27.3051] = +0.00165;
605  (*this)[+27.3731] = +0.00165;
606  (*this)[+27.4411] = +0.00164;
607  (*this)[+27.5091] = +0.00164;
608  (*this)[+27.5771] = +0.00164;
609  (*this)[+27.6451] = +0.00163;
610  (*this)[+27.7131] = +0.00163;
611  (*this)[+27.7811] = +0.00162;
612  (*this)[+27.8491] = +0.00162;
613  (*this)[+27.9171] = +0.00161;
614  (*this)[+27.9851] = +0.00160;
615  (*this)[+28.0531] = +0.00160;
616  (*this)[+28.1211] = +0.00159;
617  (*this)[+28.1891] = +0.00159;
618  (*this)[+28.2571] = +0.00158;
619  (*this)[+28.3251] = +0.00157;
620  (*this)[+28.3931] = +0.00157;
621  (*this)[+28.4611] = +0.00156;
622  (*this)[+28.5291] = +0.00156;
623  (*this)[+28.5971] = +0.00155;
624  (*this)[+28.6651] = +0.00153;
625  (*this)[+28.7331] = +0.00152;
626  (*this)[+28.8011] = +0.00151;
627  (*this)[+28.8691] = +0.00150;
628  (*this)[+28.9371] = +0.00149;
629  (*this)[+29.0051] = +0.00148;
630  (*this)[+29.0731] = +0.00147;
631  (*this)[+29.1411] = +0.00146;
632  (*this)[+29.2091] = +0.00144;
633  (*this)[+29.2771] = +0.00143;
634  (*this)[+29.3451] = +0.00142;
635  (*this)[+29.4131] = +0.00141;
636  (*this)[+29.4811] = +0.00140;
637  (*this)[+29.5491] = +0.00139;
638  (*this)[+29.6171] = +0.00137;
639  (*this)[+29.6851] = +0.00136;
640  (*this)[+29.7531] = +0.00134;
641  (*this)[+29.8211] = +0.00133;
642  (*this)[+29.8891] = +0.00131;
643  (*this)[+29.9571] = +0.00130;
644  (*this)[+30.0251] = +0.00128;
645  (*this)[+30.0931] = +0.00127;
646  (*this)[+30.1611] = +0.00125;
647  (*this)[+30.2291] = +0.00124;
648  (*this)[+30.2971] = +0.00122;
649  (*this)[+30.3651] = +0.00121;
650  (*this)[+30.4331] = +0.00119;
651  (*this)[+30.5011] = +0.00118;
652  (*this)[+30.5691] = +0.00116;
653  (*this)[+30.6371] = +0.00115;
654  (*this)[+30.7051] = +0.00113;
655  (*this)[+30.7731] = +0.00112;
656  (*this)[+30.8411] = +0.00110;
657  (*this)[+30.9091] = +0.00109;
658  (*this)[+30.9771] = +0.00107;
659  (*this)[+31.0451] = +0.00106;
660  (*this)[+31.1131] = +0.00105;
661  (*this)[+31.1811] = +0.00103;
662  (*this)[+31.2491] = +0.00102;
663  (*this)[+31.3171] = +0.00100;
664  (*this)[+31.3851] = +0.00099;
665  (*this)[+31.4531] = +0.00097;
666  (*this)[+31.5211] = +0.00096;
667  (*this)[+31.5891] = +0.00094;
668  (*this)[+31.6571] = +0.00093;
669  (*this)[+31.7251] = +0.00091;
670  (*this)[+31.7931] = +0.00090;
671  (*this)[+31.8611] = +0.00088;
672  (*this)[+31.9291] = +0.00087;
673  (*this)[+31.9971] = +0.00086;
674  (*this)[+32.0651] = +0.00084;
675  (*this)[+32.1331] = +0.00083;
676  (*this)[+32.2011] = +0.00081;
677  (*this)[+32.2691] = +0.00080;
678  (*this)[+32.3371] = +0.00079;
679  (*this)[+32.4051] = +0.00077;
680  (*this)[+32.4731] = +0.00076;
681  (*this)[+32.5411] = +0.00074;
682  (*this)[+32.6091] = +0.00073;
683  (*this)[+32.6771] = +0.00072;
684  (*this)[+32.7451] = +0.00071;
685  (*this)[+32.8131] = +0.00069;
686  (*this)[+32.8811] = +0.00068;
687  (*this)[+32.9491] = +0.00067;
688  (*this)[+33.0171] = +0.00066;
689  (*this)[+33.0851] = +0.00064;
690  (*this)[+33.1531] = +0.00063;
691  (*this)[+33.2211] = +0.00062;
692  (*this)[+33.2891] = +0.00061;
693  (*this)[+33.3571] = +0.00060;
694  (*this)[+33.4251] = +0.00058;
695  (*this)[+33.4931] = +0.00057;
696  (*this)[+33.5611] = +0.00056;
697  (*this)[+33.6291] = +0.00055;
698  (*this)[+33.6971] = +0.00054;
699  (*this)[+33.7651] = +0.00053;
700  (*this)[+33.8331] = +0.00052;
701  (*this)[+33.9011] = +0.00051;
702  (*this)[+33.9691] = +0.00050;
703  (*this)[+34.0371] = +0.00049;
704  (*this)[+34.1051] = +0.00048;
705  (*this)[+34.1731] = +0.00047;
706  (*this)[+34.2411] = +0.00046;
707  (*this)[+34.3091] = +0.00045;
708  (*this)[+34.3771] = +0.00044;
709  (*this)[+34.4451] = +0.00043;
710  (*this)[+34.5131] = +0.00042;
711  (*this)[+34.5811] = +0.00042;
712  (*this)[+34.6491] = +0.00041;
713  (*this)[+34.7171] = +0.00040;
714  (*this)[+34.7851] = +0.00039;
715  (*this)[+34.8531] = +0.00039;
716  (*this)[+34.9211] = +0.00038;
717  (*this)[+34.9891] = +0.00037;
718  (*this)[+35.0571] = +0.00036;
719  (*this)[+35.1251] = +0.00036;
720  (*this)[+35.1931] = +0.00035;
721  (*this)[+35.2611] = +0.00034;
722  (*this)[+35.3291] = +0.00034;
723  (*this)[+35.3971] = +0.00033;
724  (*this)[+35.4651] = +0.00032;
725  (*this)[+35.5331] = +0.00031;
726  (*this)[+35.6011] = +0.00031;
727  (*this)[+35.6691] = +0.00030;
728  (*this)[+35.7371] = +0.00030;
729  (*this)[+35.8051] = +0.00029;
730  (*this)[+35.8731] = +0.00029;
731  (*this)[+35.9411] = +0.00028;
732  (*this)[+36.0091] = +0.00028;
733  (*this)[+36.0771] = +0.00027;
734  (*this)[+36.1451] = +0.00026;
735  (*this)[+36.2131] = +0.00026;
736  (*this)[+36.2811] = +0.00025;
737  (*this)[+36.3491] = +0.00025;
738  (*this)[+36.4171] = +0.00024;
739  (*this)[+36.4851] = +0.00024;
740  (*this)[+36.5531] = +0.00023;
741  (*this)[+36.6211] = +0.00023;
742  (*this)[+36.6891] = +0.00022;
743  (*this)[+36.7571] = +0.00022;
744  (*this)[+36.8251] = +0.00021;
745  (*this)[+36.8931] = +0.00021;
746  (*this)[+36.9611] = +0.00021;
747  (*this)[+37.0291] = +0.00020;
748  (*this)[+37.0971] = +0.00020;
749  (*this)[+37.1651] = +0.00019;
750  (*this)[+37.2331] = +0.00019;
751  (*this)[+37.3011] = +0.00019;
752  (*this)[+37.3691] = +0.00018;
753  (*this)[+37.4371] = +0.00018;
754  (*this)[+37.5051] = +0.00017;
755  (*this)[+37.5731] = +0.00017;
756  (*this)[+37.6411] = +0.00017;
757  (*this)[+37.7091] = +0.00016;
758  (*this)[+37.7771] = +0.00016;
759  (*this)[+37.8451] = +0.00016;
760  (*this)[+37.9131] = +0.00015;
761  (*this)[+37.9811] = +0.00015;
762  (*this)[+38.0491] = +0.00015;
763  (*this)[+38.1171] = +0.00014;
764  (*this)[+38.1851] = +0.00014;
765  (*this)[+38.2531] = +0.00014;
766  (*this)[+38.3211] = +0.00013;
767  (*this)[+38.3891] = +0.00013;
768  (*this)[+38.4571] = +0.00013;
769  (*this)[+38.5251] = +0.00013;
770  (*this)[+38.5931] = +0.00012;
771  (*this)[+38.6611] = +0.00012;
772  (*this)[+38.7291] = +0.00012;
773  (*this)[+38.7971] = +0.00012;
774  (*this)[+38.8651] = +0.00012;
775  (*this)[+38.9331] = +0.00011;
776  (*this)[+39.0011] = +0.00011;
777  (*this)[+39.0691] = +0.00011;
778  (*this)[+39.1371] = +0.00011;
779  (*this)[+39.2051] = +0.00011;
780  (*this)[+39.2731] = +0.00010;
781  (*this)[+39.3411] = +0.00010;
782  (*this)[+39.4091] = +0.00010;
783  (*this)[+39.4771] = +0.00010;
784  (*this)[+39.5451] = +0.00010;
785  (*this)[+39.6131] = +0.00010;
786  (*this)[+39.6811] = +0.00009;
787  (*this)[+39.7491] = +0.00009;
788  (*this)[+39.8171] = +0.00009;
789  (*this)[+39.8851] = +0.00009;
790  (*this)[+39.9531] = +0.00009;
791  (*this)[+40.0211] = +0.00009;
792  (*this)[+40.0891] = +0.00009;
793  (*this)[+40.1571] = +0.00009;
794  (*this)[+40.2251] = +0.00008;
795  (*this)[+40.2931] = +0.00008;
796  (*this)[+40.3611] = +0.00008;
797  (*this)[+40.4291] = +0.00008;
798  (*this)[+40.4971] = +0.00008;
799  (*this)[+40.5651] = +0.00008;
800  (*this)[+40.6331] = +0.00008;
801  (*this)[+40.7011] = +0.00008;
802  (*this)[+40.7691] = +0.00008;
803  (*this)[+40.8371] = +0.00007;
804  (*this)[+40.9051] = +0.00007;
805  (*this)[+40.9731] = +0.00007;
806  (*this)[+41.0411] = +0.00007;
807  (*this)[+41.1091] = +0.00007;
808  (*this)[+41.1771] = +0.00007;
809  (*this)[+41.2451] = +0.00007;
810  (*this)[+41.3131] = +0.00007;
811  (*this)[+41.3811] = +0.00007;
812  (*this)[+41.4491] = +0.00007;
813  (*this)[+41.5171] = +0.00007;
814  (*this)[+41.5851] = +0.00006;
815  (*this)[+41.6531] = +0.00006;
816  (*this)[+41.7211] = +0.00006;
817  (*this)[+41.7891] = +0.00006;
818  (*this)[+41.8571] = +0.00006;
819  (*this)[+41.9251] = +0.00006;
820  (*this)[+41.9931] = +0.00006;
821  (*this)[+42.0611] = +0.00006;
822  (*this)[+42.1291] = +0.00006;
823  (*this)[+42.1971] = +0.00006;
824  (*this)[+42.2651] = +0.00006;
825  (*this)[+42.3331] = +0.00006;
826  (*this)[+42.4011] = +0.00005;
827  (*this)[+42.4691] = +0.00005;
828  (*this)[+42.5371] = +0.00005;
829  (*this)[+42.6051] = +0.00005;
830  (*this)[+42.6731] = +0.00005;
831  (*this)[+42.7411] = +0.00005;
832  (*this)[+42.8091] = +0.00005;
833  (*this)[+42.8771] = +0.00005;
834  (*this)[+42.9451] = +0.00005;
835  (*this)[+43.0131] = +0.00005;
836  (*this)[+43.0811] = +0.00005;
837  (*this)[+43.1491] = +0.00005;
838  (*this)[+43.2171] = +0.00005;
839  (*this)[+43.2851] = +0.00005;
840  (*this)[+43.3531] = +0.00005;
841  (*this)[+43.4211] = +0.00005;
842  (*this)[+43.4891] = +0.00005;
843  (*this)[+43.5571] = +0.00005;
844  (*this)[+43.6251] = +0.00005;
845  (*this)[+43.6931] = +0.00005;
846  (*this)[+43.7611] = +0.00005;
847  (*this)[+43.8291] = +0.00005;
848  (*this)[+43.8971] = +0.00005;
849  (*this)[+43.9651] = +0.00005;
850  (*this)[+44.0331] = +0.00005;
851  (*this)[+44.1011] = +0.00005;
852  (*this)[+44.1691] = +0.00005;
853  (*this)[+44.2371] = +0.00005;
854  (*this)[+44.3051] = +0.00005;
855  (*this)[+44.3731] = +0.00005;
856  (*this)[+44.4411] = +0.00005;
857  (*this)[+44.5091] = +0.00005;
858  (*this)[+44.5771] = +0.00005;
859  (*this)[+44.6451] = +0.00005;
860  (*this)[+44.7131] = +0.00005;
861  (*this)[+44.7811] = +0.00005;
862  (*this)[+44.8491] = +0.00005;
863  (*this)[+44.9171] = +0.00005;
864  (*this)[+44.9851] = +0.00005;
865  (*this)[+45.0531] = +0.00005;
866  (*this)[+45.1211] = +0.00005;
867  (*this)[+45.1891] = +0.00005;
868  (*this)[+45.2571] = +0.00005;
869  (*this)[+45.3251] = +0.00005;
870  (*this)[+45.3931] = +0.00005;
871  (*this)[+45.4611] = +0.00005;
872  (*this)[+45.5291] = +0.00005;
873  (*this)[+45.5971] = +0.00005;
874  (*this)[+45.6651] = +0.00005;
875  (*this)[+45.7331] = +0.00005;
876  (*this)[+45.8011] = +0.00005;
877  (*this)[+45.8691] = +0.00005;
878  (*this)[+45.9371] = +0.00005;
879  (*this)[+46.0051] = +0.00005;
880  (*this)[+46.0731] = +0.00005;
881  (*this)[+46.1411] = +0.00005;
882  (*this)[+46.2091] = +0.00005;
883  (*this)[+46.2771] = +0.00005;
884  (*this)[+46.3451] = +0.00005;
885  (*this)[+46.4131] = +0.00005;
886  (*this)[+46.4811] = +0.00005;
887  (*this)[+46.5491] = +0.00005;
888  (*this)[+46.6171] = +0.00005;
889  (*this)[+46.6851] = +0.00005;
890  (*this)[+46.7531] = +0.00005;
891  (*this)[+46.8211] = +0.00005;
892  (*this)[+46.8891] = +0.00005;
893  (*this)[+46.9571] = +0.00005;
894  (*this)[+47.0251] = +0.00005;
895  (*this)[+47.0931] = +0.00005;
896  (*this)[+47.1611] = +0.00005;
897  (*this)[+47.2291] = +0.00005;
898  (*this)[+47.2971] = +0.00005;
899  (*this)[+47.3651] = +0.00005;
900  (*this)[+47.4331] = +0.00005;
901  (*this)[+47.5011] = +0.00005;
902  (*this)[+47.5691] = +0.00005;
903  (*this)[+47.6371] = +0.00005;
904  (*this)[+47.7051] = +0.00005;
905  (*this)[+47.7731] = +0.00005;
906  (*this)[+47.8411] = +0.00005;
907  (*this)[+47.9091] = +0.00005;
908  (*this)[+47.9771] = +0.00005;
909  (*this)[+48.0451] = +0.00005;
910  (*this)[+48.1131] = +0.00005;
911  (*this)[+48.1811] = +0.00005;
912  (*this)[+48.2491] = +0.00005;
913  (*this)[+48.3171] = +0.00005;
914  (*this)[+48.3851] = +0.00005;
915  (*this)[+48.4531] = +0.00005;
916  (*this)[+48.5211] = +0.00005;
917  (*this)[+48.5891] = +0.00005;
918  (*this)[+48.6571] = +0.00005;
919  (*this)[+48.7251] = +0.00005;
920  (*this)[+48.7931] = +0.00005;
921  (*this)[+48.8611] = +0.00005;
922  (*this)[+48.9291] = +0.00005;
923  (*this)[+48.9971] = +0.00005;
924  (*this)[+49.0651] = +0.00005;
925  (*this)[+49.1331] = +0.00005;
926  (*this)[+49.2011] = +0.00005;
927  (*this)[+49.2691] = +0.00005;
928  (*this)[+49.3371] = +0.00005;
929  (*this)[+49.4051] = +0.00005;
930  (*this)[+49.4731] = +0.00005;
931  (*this)[+49.5411] = +0.00005;
932  (*this)[+49.6091] = +0.00005;
933  (*this)[+49.6771] = +0.00005;
934  (*this)[+49.7451] = +0.00005;
935  (*this)[+49.8131] = +0.00005;
936  (*this)[+49.8811] = +0.00005;
937  (*this)[+49.9491] = +0.00005;
938  (*this)[+50.0171] = +0.00005;
939  (*this)[+50.0851] = +0.00005;
940  (*this)[+50.1531] = +0.00005;
941  (*this)[+50.2211] = +0.00005;
942  (*this)[+50.2891] = +0.00005;
943  (*this)[+50.3571] = +0.00005;
944  (*this)[+50.4251] = +0.00005;
945  (*this)[+50.4931] = +0.00005;
946  (*this)[+50.5611] = +0.00005;
947  (*this)[+50.6291] = +0.00005;
948  (*this)[+50.6971] = +0.00005;
949  (*this)[+50.7651] = +0.00005;
950  (*this)[+50.8331] = +0.00005;
951  (*this)[+50.9011] = +0.00005;
952  (*this)[+50.9691] = +0.00005;
953  (*this)[+51.0371] = +0.00005;
954  (*this)[+51.1051] = +0.00005;
955  (*this)[+51.1731] = +0.00005;
956  (*this)[+51.2411] = +0.00005;
957  (*this)[+51.3091] = +0.00005;
958  (*this)[+51.3771] = +0.00005;
959  (*this)[+51.4451] = +0.00005;
960  (*this)[+51.5131] = +0.00004;
961  (*this)[+51.5811] = +0.00004;
962  (*this)[+51.6491] = +0.00004;
963  (*this)[+51.7171] = +0.00004;
964  (*this)[+51.7851] = +0.00004;
965  (*this)[+51.8531] = +0.00004;
966  (*this)[+51.9211] = +0.00004;
967  (*this)[+51.9891] = +0.00004;
968  (*this)[+52.0571] = +0.00004;
969  (*this)[+52.1251] = +0.00004;
970  (*this)[+52.1931] = +0.00004;
971  (*this)[+52.2611] = +0.00004;
972  (*this)[+52.3291] = +0.00004;
973  (*this)[+52.3971] = +0.00004;
974  (*this)[+52.4651] = +0.00004;
975  (*this)[+52.5331] = +0.00004;
976  (*this)[+52.6011] = +0.00004;
977  (*this)[+52.6691] = +0.00004;
978  (*this)[+52.7371] = +0.00004;
979  (*this)[+52.8051] = +0.00004;
980  (*this)[+52.8731] = +0.00004;
981  (*this)[+52.9411] = +0.00004;
982  (*this)[+53.0091] = +0.00004;
983  (*this)[+53.0771] = +0.00004;
984  (*this)[+53.1451] = +0.00004;
985  (*this)[+53.2131] = +0.00004;
986  (*this)[+53.2811] = +0.00005;
987  (*this)[+53.3491] = +0.00005;
988  (*this)[+53.4171] = +0.00005;
989  (*this)[+53.4851] = +0.00005;
990  (*this)[+53.5531] = +0.00005;
991  (*this)[+53.6211] = +0.00005;
992  (*this)[+53.6891] = +0.00005;
993  (*this)[+53.7571] = +0.00004;
994  (*this)[+53.8251] = +0.00004;
995  (*this)[+53.8931] = +0.00004;
996  (*this)[+53.9611] = +0.00004;
997  (*this)[+54.0291] = +0.00004;
998  (*this)[+54.0971] = +0.00004;
999  (*this)[+54.1651] = +0.00004;
1000  (*this)[+54.2331] = +0.00004;
1001  (*this)[+54.3011] = +0.00004;
1002  (*this)[+54.3691] = +0.00004;
1003  (*this)[+54.4371] = +0.00004;
1004  (*this)[+54.5051] = +0.00004;
1005  (*this)[+54.5731] = +0.00004;
1006  (*this)[+54.6411] = +0.00004;
1007  (*this)[+54.7091] = +0.00004;
1008  (*this)[+54.7771] = +0.00004;
1009  (*this)[+54.8451] = +0.00004;
1010  (*this)[+54.9131] = +0.00004;
1011  (*this)[+54.9811] = +0.00004;
1012  (*this)[+55.0491] = +0.00004;
1013  (*this)[+55.1171] = +0.00004;
1014  (*this)[+55.1851] = +0.00004;
1015  (*this)[+55.2531] = +0.00004;
1016  (*this)[+55.3211] = +0.00004;
1017  (*this)[+55.3891] = +0.00004;
1018  (*this)[+55.4571] = +0.00004;
1019  (*this)[+55.5251] = +0.00004;
1020  (*this)[+55.5931] = +0.00004;
1021  (*this)[+55.6611] = +0.00003;
1022  (*this)[+55.7291] = +0.00003;
1023  (*this)[+55.7971] = +0.00003;
1024  (*this)[+55.8651] = +0.00003;
1025  (*this)[+55.9331] = +0.00002;
1026  (*this)[+56.0011] = +0.00002;
1027  (*this)[+56.0691] = +0.00002;
1028  (*this)[+56.1371] = +0.00002;
1029  (*this)[+56.2051] = +0.00001;
1030  (*this)[+56.2731] = +0.00001;
1031  (*this)[+56.3411] = +0.00001;
1032  (*this)[+56.4091] = +0.00001;
1033  (*this)[+56.4771] = +0.00000;
1034  (*this)[+56.5451] = +0.00000;
1035 
1036  compile();
1037 
1038  setExceptionHandler(new JDefaultResult(0.0));
1039  };
1040  };
1041 
1042 
1043  /**
1044  * Function object of transition time prabability.
1045  */
1046  static const JTransitionTimeProbability getTransitionTimeProbability;
1047 
1048 
1049  /**
1050  * Auxiliary class to generate twofold PMT transition times.
1051  */
1052  class JTransitionTimeProbability2D :
1053  public JTOOLS::JGridSplineFunction1D_t
1054  {
1055  public:
1056  /**
1057  * Default constructor.
1058  */
1059  JTransitionTimeProbability2D() :
1060  JTOOLS::JGridSplineFunction1D_t()
1061  {
1062  using namespace std;
1063 
1064  const double xmin = min( getTransitionTimeProbability.getXmin(), -getTransitionTimeProbability.getXmax());
1065  const double xmax = max(-getTransitionTimeProbability.getXmin(), getTransitionTimeProbability.getXmax());
1066  const double dx = (getTransitionTimeProbability.getXmax() - getTransitionTimeProbability.getXmin()) / (getTransitionTimeProbability.getSize() - 1);
1067 
1068  double Y = 0.0;
1069 
1070  for (double x1 = xmin; x1 < xmax + 0.5*dx; x1 += dx) {
1071 
1072  double y = 0.0;
1073 
1074  for (double x2 = xmin; x2 < xmax + 0.5*dx; x2 += dx) {
1075  y += getTransitionTimeProbability(x2) * getTransitionTimeProbability( x1 - x2);
1076  y += getTransitionTimeProbability(x2) * getTransitionTimeProbability(-x1 - x2);
1077  }
1078 
1079  (*this)[x1] = y;
1080 
1081  Y += y;
1082  }
1083 
1084  this->div(dx*Y);
1085 
1086  compile();
1087 
1088  setExceptionHandler(new JDefaultResult(0.0));
1089  }
1090  };
1091 
1092 
1093  /**
1094  * Function object of twofold transition time prabability.
1095  */
1096  static const JTransitionTimeProbability2D getTransitionTimeProbability2D;
1097 }
1098 
1099 #endif
JDETECTOR::JTransitionTimeProbability2D
Auxiliary class to generate twofold PMT transition times.
Definition: JPMTTransitTimeProbability.hh:1052
JDETECTOR::JTransitionTimeProbability
Probability distibution of the PMT transition times based on measurements by Oleg Kalekin et al...
Definition: JPMTTransitTimeProbability.hh:24
JTOOLS::JGridSplineFunction1D_t
Type definition of a spline interpolation based on a JGridCollection with result type double...
Definition: JFunction1D_t.hh:59
JDETECTOR::getTransitionTimeProbability
static const JTransitionTimeProbability getTransitionTimeProbability
Function object of transition time prabability.
Definition: JPMTTransitTimeProbability.hh:1046
JDETECTOR::getTransitionTimeProbability2D
static const JTransitionTimeProbability2D getTransitionTimeProbability2D
Function object of twofold transition time prabability.
Definition: JPMTTransitTimeProbability.hh:1096
JFunction1D_t.hh
JTOOLS::JFunctional< JSplineElement2D< double, double >::abscissa_type, double >::setExceptionHandler
void setExceptionHandler(const JSupervisor &supervisor)
Set the supervisor for handling of exceptions.
Definition: JFunctional.hh:286
Generated by   doxygen 1.8.5