Jpp  pmt_effective_area_update
the software that should make you happy
 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  * Auxiliary base class for probability distibution of the PMT transition times based on measurements.\n
20  * For a given transition time [ns], this function returns the correspoding probability.
21  */
22  struct JTransitionTimeProbability_t :
23  public JTOOLS::JGridSplineFunction1D_t
24  {
25  /**
26  * Default constructor.
27  */
28  JTransitionTimeProbability_t()
29  {
30  setExceptionHandler(new JDefaultResult(0.0));
31  }
32  };
33 
34 
35  /**
36  * Implementation of probability distibution of the PMT transition times based on measurements
37  * by Oleg Kalekin et al. (see also examples/JDetector/JTTS.cc).
38  */
39  struct JTransitionTimeProbability :
40  public JTransitionTimeProbability_t
41  {
42  /**
43  * Default constructor.
44  */
45  JTransitionTimeProbability()
46  {
47  (*this)[-11.4549] = +0.00018;
48  (*this)[-11.3869] = +0.00018;
49  (*this)[-11.3189] = +0.00018;
50  (*this)[-11.2509] = +0.00018;
51  (*this)[-11.1829] = +0.00018;
52  (*this)[-11.1149] = +0.00018;
53  (*this)[-11.0469] = +0.00018;
54  (*this)[-10.9789] = +0.00018;
55  (*this)[-10.9109] = +0.00017;
56  (*this)[-10.8429] = +0.00017;
57  (*this)[-10.7749] = +0.00017;
58  (*this)[-10.7069] = +0.00017;
59  (*this)[-10.6389] = +0.00017;
60  (*this)[-10.5709] = +0.00017;
61  (*this)[-10.5029] = +0.00017;
62  (*this)[-10.4349] = +0.00017;
63  (*this)[-10.3669] = +0.00017;
64  (*this)[-10.2989] = +0.00017;
65  (*this)[-10.2309] = +0.00017;
66  (*this)[-10.1629] = +0.00017;
67  (*this)[-10.0949] = +0.00017;
68  (*this)[-10.0269] = +0.00017;
69  (*this)[ -9.9589] = +0.00017;
70  (*this)[ -9.8909] = +0.00017;
71  (*this)[ -9.8229] = +0.00017;
72  (*this)[ -9.7549] = +0.00017;
73  (*this)[ -9.6869] = +0.00017;
74  (*this)[ -9.6189] = +0.00017;
75  (*this)[ -9.5509] = +0.00017;
76  (*this)[ -9.4829] = +0.00017;
77  (*this)[ -9.4149] = +0.00017;
78  (*this)[ -9.3469] = +0.00017;
79  (*this)[ -9.2789] = +0.00018;
80  (*this)[ -9.2109] = +0.00018;
81  (*this)[ -9.1429] = +0.00018;
82  (*this)[ -9.0749] = +0.00018;
83  (*this)[ -9.0069] = +0.00018;
84  (*this)[ -8.9389] = +0.00018;
85  (*this)[ -8.8709] = +0.00018;
86  (*this)[ -8.8029] = +0.00018;
87  (*this)[ -8.7349] = +0.00018;
88  (*this)[ -8.6669] = +0.00018;
89  (*this)[ -8.5989] = +0.00018;
90  (*this)[ -8.5309] = +0.00018;
91  (*this)[ -8.4629] = +0.00018;
92  (*this)[ -8.3949] = +0.00019;
93  (*this)[ -8.3269] = +0.00019;
94  (*this)[ -8.2589] = +0.00020;
95  (*this)[ -8.1909] = +0.00020;
96  (*this)[ -8.1229] = +0.00020;
97  (*this)[ -8.0549] = +0.00021;
98  (*this)[ -7.9869] = +0.00021;
99  (*this)[ -7.9189] = +0.00022;
100  (*this)[ -7.8509] = +0.00022;
101  (*this)[ -7.7829] = +0.00022;
102  (*this)[ -7.7149] = +0.00023;
103  (*this)[ -7.6469] = +0.00023;
104  (*this)[ -7.5789] = +0.00023;
105  (*this)[ -7.5109] = +0.00024;
106  (*this)[ -7.4429] = +0.00025;
107  (*this)[ -7.3749] = +0.00026;
108  (*this)[ -7.3069] = +0.00028;
109  (*this)[ -7.2389] = +0.00030;
110  (*this)[ -7.1709] = +0.00031;
111  (*this)[ -7.1029] = +0.00033;
112  (*this)[ -7.0349] = +0.00035;
113  (*this)[ -6.9669] = +0.00037;
114  (*this)[ -6.8989] = +0.00038;
115  (*this)[ -6.8309] = +0.00040;
116  (*this)[ -6.7629] = +0.00042;
117  (*this)[ -6.6949] = +0.00044;
118  (*this)[ -6.6269] = +0.00045;
119  (*this)[ -6.5589] = +0.00047;
120  (*this)[ -6.4909] = +0.00049;
121  (*this)[ -6.4229] = +0.00054;
122  (*this)[ -6.3549] = +0.00063;
123  (*this)[ -6.2869] = +0.00073;
124  (*this)[ -6.2189] = +0.00082;
125  (*this)[ -6.1509] = +0.00091;
126  (*this)[ -6.0829] = +0.00101;
127  (*this)[ -6.0149] = +0.00110;
128  (*this)[ -5.9469] = +0.00119;
129  (*this)[ -5.8789] = +0.00128;
130  (*this)[ -5.8109] = +0.00138;
131  (*this)[ -5.7429] = +0.00147;
132  (*this)[ -5.6749] = +0.00156;
133  (*this)[ -5.6069] = +0.00166;
134  (*this)[ -5.5389] = +0.00175;
135  (*this)[ -5.4709] = +0.00184;
136  (*this)[ -5.4029] = +0.00219;
137  (*this)[ -5.3349] = +0.00261;
138  (*this)[ -5.2669] = +0.00304;
139  (*this)[ -5.1989] = +0.00346;
140  (*this)[ -5.1309] = +0.00388;
141  (*this)[ -5.0629] = +0.00431;
142  (*this)[ -4.9949] = +0.00473;
143  (*this)[ -4.9269] = +0.00516;
144  (*this)[ -4.8589] = +0.00558;
145  (*this)[ -4.7909] = +0.00600;
146  (*this)[ -4.7229] = +0.00643;
147  (*this)[ -4.6549] = +0.00685;
148  (*this)[ -4.5869] = +0.00728;
149  (*this)[ -4.5189] = +0.00770;
150  (*this)[ -4.4509] = +0.00818;
151  (*this)[ -4.3829] = +0.00959;
152  (*this)[ -4.3149] = +0.01100;
153  (*this)[ -4.2469] = +0.01241;
154  (*this)[ -4.1789] = +0.01382;
155  (*this)[ -4.1109] = +0.01523;
156  (*this)[ -4.0429] = +0.01664;
157  (*this)[ -3.9749] = +0.01805;
158  (*this)[ -3.9069] = +0.01946;
159  (*this)[ -3.8389] = +0.02087;
160  (*this)[ -3.7709] = +0.02228;
161  (*this)[ -3.7029] = +0.02369;
162  (*this)[ -3.6349] = +0.02510;
163  (*this)[ -3.5669] = +0.02651;
164  (*this)[ -3.4989] = +0.02792;
165  (*this)[ -3.4309] = +0.02993;
166  (*this)[ -3.3629] = +0.03302;
167  (*this)[ -3.2949] = +0.03611;
168  (*this)[ -3.2269] = +0.03920;
169  (*this)[ -3.1589] = +0.04230;
170  (*this)[ -3.0909] = +0.04539;
171  (*this)[ -3.0229] = +0.04848;
172  (*this)[ -2.9549] = +0.05157;
173  (*this)[ -2.8869] = +0.05467;
174  (*this)[ -2.8189] = +0.05776;
175  (*this)[ -2.7509] = +0.06085;
176  (*this)[ -2.6829] = +0.06394;
177  (*this)[ -2.6149] = +0.06704;
178  (*this)[ -2.5469] = +0.07013;
179  (*this)[ -2.4789] = +0.07322;
180  (*this)[ -2.4109] = +0.07709;
181  (*this)[ -2.3429] = +0.08138;
182  (*this)[ -2.2749] = +0.08567;
183  (*this)[ -2.2069] = +0.08997;
184  (*this)[ -2.1389] = +0.09426;
185  (*this)[ -2.0709] = +0.09855;
186  (*this)[ -2.0029] = +0.10284;
187  (*this)[ -1.9349] = +0.10713;
188  (*this)[ -1.8669] = +0.11143;
189  (*this)[ -1.7989] = +0.11572;
190  (*this)[ -1.7309] = +0.12001;
191  (*this)[ -1.6629] = +0.12430;
192  (*this)[ -1.5949] = +0.12860;
193  (*this)[ -1.5269] = +0.13289;
194  (*this)[ -1.4589] = +0.13718;
195  (*this)[ -1.3909] = +0.14044;
196  (*this)[ -1.3229] = +0.14364;
197  (*this)[ -1.2549] = +0.14684;
198  (*this)[ -1.1869] = +0.15003;
199  (*this)[ -1.1189] = +0.15323;
200  (*this)[ -1.0509] = +0.15643;
201  (*this)[ -0.9829] = +0.15962;
202  (*this)[ -0.9149] = +0.16282;
203  (*this)[ -0.8469] = +0.16602;
204  (*this)[ -0.7789] = +0.16921;
205  (*this)[ -0.7109] = +0.17241;
206  (*this)[ -0.6429] = +0.17561;
207  (*this)[ -0.5749] = +0.17880;
208  (*this)[ -0.5069] = +0.18200;
209  (*this)[ -0.4389] = +0.18449;
210  (*this)[ -0.3709] = +0.18467;
211  (*this)[ -0.3029] = +0.18485;
212  (*this)[ -0.2349] = +0.18503;
213  (*this)[ -0.1669] = +0.18521;
214  (*this)[ -0.0989] = +0.18539;
215  (*this)[ -0.0309] = +0.18557;
216  (*this)[ +0.0371] = +0.18575;
217  (*this)[ +0.1051] = +0.18593;
218  (*this)[ +0.1731] = +0.18611;
219  (*this)[ +0.2411] = +0.18629;
220  (*this)[ +0.3091] = +0.18647;
221  (*this)[ +0.3771] = +0.18665;
222  (*this)[ +0.4451] = +0.18683;
223  (*this)[ +0.5131] = +0.18701;
224  (*this)[ +0.5811] = +0.18581;
225  (*this)[ +0.6491] = +0.18339;
226  (*this)[ +0.7171] = +0.18097;
227  (*this)[ +0.7851] = +0.17854;
228  (*this)[ +0.8531] = +0.17612;
229  (*this)[ +0.9211] = +0.17370;
230  (*this)[ +0.9891] = +0.17127;
231  (*this)[ +1.0571] = +0.16885;
232  (*this)[ +1.1251] = +0.16643;
233  (*this)[ +1.1931] = +0.16400;
234  (*this)[ +1.2611] = +0.16158;
235  (*this)[ +1.3291] = +0.15916;
236  (*this)[ +1.3971] = +0.15674;
237  (*this)[ +1.4651] = +0.15431;
238  (*this)[ +1.5331] = +0.15189;
239  (*this)[ +1.6011] = +0.14861;
240  (*this)[ +1.6691] = +0.14514;
241  (*this)[ +1.7371] = +0.14168;
242  (*this)[ +1.8051] = +0.13822;
243  (*this)[ +1.8731] = +0.13475;
244  (*this)[ +1.9411] = +0.13129;
245  (*this)[ +2.0091] = +0.12783;
246  (*this)[ +2.0771] = +0.12436;
247  (*this)[ +2.1451] = +0.12090;
248  (*this)[ +2.2131] = +0.11743;
249  (*this)[ +2.2811] = +0.11397;
250  (*this)[ +2.3491] = +0.11051;
251  (*this)[ +2.4171] = +0.10704;
252  (*this)[ +2.4851] = +0.10358;
253  (*this)[ +2.5531] = +0.10015;
254  (*this)[ +2.6211] = +0.09696;
255  (*this)[ +2.6891] = +0.09377;
256  (*this)[ +2.7571] = +0.09059;
257  (*this)[ +2.8251] = +0.08740;
258  (*this)[ +2.8931] = +0.08421;
259  (*this)[ +2.9611] = +0.08103;
260  (*this)[ +3.0291] = +0.07784;
261  (*this)[ +3.0971] = +0.07465;
262  (*this)[ +3.1651] = +0.07146;
263  (*this)[ +3.2331] = +0.06828;
264  (*this)[ +3.3011] = +0.06509;
265  (*this)[ +3.3691] = +0.06190;
266  (*this)[ +3.4371] = +0.05871;
267  (*this)[ +3.5051] = +0.05553;
268  (*this)[ +3.5731] = +0.05280;
269  (*this)[ +3.6411] = +0.05072;
270  (*this)[ +3.7091] = +0.04864;
271  (*this)[ +3.7771] = +0.04656;
272  (*this)[ +3.8451] = +0.04448;
273  (*this)[ +3.9131] = +0.04240;
274  (*this)[ +3.9811] = +0.04032;
275  (*this)[ +4.0491] = +0.03823;
276  (*this)[ +4.1171] = +0.03615;
277  (*this)[ +4.1851] = +0.03407;
278  (*this)[ +4.2531] = +0.03199;
279  (*this)[ +4.3211] = +0.02991;
280  (*this)[ +4.3891] = +0.02783;
281  (*this)[ +4.4571] = +0.02575;
282  (*this)[ +4.5251] = +0.02367;
283  (*this)[ +4.5931] = +0.02235;
284  (*this)[ +4.6611] = +0.02133;
285  (*this)[ +4.7291] = +0.02032;
286  (*this)[ +4.7971] = +0.01930;
287  (*this)[ +4.8651] = +0.01829;
288  (*this)[ +4.9331] = +0.01727;
289  (*this)[ +5.0011] = +0.01626;
290  (*this)[ +5.0691] = +0.01524;
291  (*this)[ +5.1371] = +0.01423;
292  (*this)[ +5.2051] = +0.01321;
293  (*this)[ +5.2731] = +0.01220;
294  (*this)[ +5.3411] = +0.01118;
295  (*this)[ +5.4091] = +0.01017;
296  (*this)[ +5.4771] = +0.00915;
297  (*this)[ +5.5451] = +0.00814;
298  (*this)[ +5.6131] = +0.00780;
299  (*this)[ +5.6811] = +0.00746;
300  (*this)[ +5.7491] = +0.00712;
301  (*this)[ +5.8171] = +0.00679;
302  (*this)[ +5.8851] = +0.00645;
303  (*this)[ +5.9531] = +0.00611;
304  (*this)[ +6.0211] = +0.00577;
305  (*this)[ +6.0891] = +0.00543;
306  (*this)[ +6.1571] = +0.00509;
307  (*this)[ +6.2251] = +0.00476;
308  (*this)[ +6.2931] = +0.00442;
309  (*this)[ +6.3611] = +0.00408;
310  (*this)[ +6.4291] = +0.00374;
311  (*this)[ +6.4971] = +0.00340;
312  (*this)[ +6.5651] = +0.00314;
313  (*this)[ +6.6331] = +0.00306;
314  (*this)[ +6.7011] = +0.00297;
315  (*this)[ +6.7691] = +0.00289;
316  (*this)[ +6.8371] = +0.00280;
317  (*this)[ +6.9051] = +0.00272;
318  (*this)[ +6.9731] = +0.00263;
319  (*this)[ +7.0411] = +0.00255;
320  (*this)[ +7.1091] = +0.00247;
321  (*this)[ +7.1771] = +0.00238;
322  (*this)[ +7.2451] = +0.00230;
323  (*this)[ +7.3131] = +0.00221;
324  (*this)[ +7.3811] = +0.00213;
325  (*this)[ +7.4491] = +0.00205;
326  (*this)[ +7.5171] = +0.00196;
327  (*this)[ +7.5851] = +0.00191;
328  (*this)[ +7.6531] = +0.00188;
329  (*this)[ +7.7211] = +0.00185;
330  (*this)[ +7.7891] = +0.00182;
331  (*this)[ +7.8571] = +0.00179;
332  (*this)[ +7.9251] = +0.00177;
333  (*this)[ +7.9931] = +0.00174;
334  (*this)[ +8.0611] = +0.00171;
335  (*this)[ +8.1291] = +0.00168;
336  (*this)[ +8.1971] = +0.00165;
337  (*this)[ +8.2651] = +0.00162;
338  (*this)[ +8.3331] = +0.00159;
339  (*this)[ +8.4011] = +0.00156;
340  (*this)[ +8.4691] = +0.00154;
341  (*this)[ +8.5371] = +0.00151;
342  (*this)[ +8.6051] = +0.00149;
343  (*this)[ +8.6731] = +0.00147;
344  (*this)[ +8.7411] = +0.00145;
345  (*this)[ +8.8091] = +0.00143;
346  (*this)[ +8.8771] = +0.00141;
347  (*this)[ +8.9451] = +0.00139;
348  (*this)[ +9.0131] = +0.00138;
349  (*this)[ +9.0811] = +0.00136;
350  (*this)[ +9.1491] = +0.00134;
351  (*this)[ +9.2171] = +0.00132;
352  (*this)[ +9.2851] = +0.00130;
353  (*this)[ +9.3531] = +0.00128;
354  (*this)[ +9.4211] = +0.00126;
355  (*this)[ +9.4891] = +0.00124;
356  (*this)[ +9.5571] = +0.00123;
357  (*this)[ +9.6251] = +0.00121;
358  (*this)[ +9.6931] = +0.00120;
359  (*this)[ +9.7611] = +0.00118;
360  (*this)[ +9.8291] = +0.00117;
361  (*this)[ +9.8971] = +0.00115;
362  (*this)[ +9.9651] = +0.00114;
363  (*this)[+10.0331] = +0.00113;
364  (*this)[+10.1011] = +0.00111;
365  (*this)[+10.1691] = +0.00110;
366  (*this)[+10.2371] = +0.00108;
367  (*this)[+10.3051] = +0.00107;
368  (*this)[+10.3731] = +0.00105;
369  (*this)[+10.4411] = +0.00104;
370  (*this)[+10.5091] = +0.00102;
371  (*this)[+10.5771] = +0.00101;
372  (*this)[+10.6451] = +0.00101;
373  (*this)[+10.7131] = +0.00100;
374  (*this)[+10.7811] = +0.00099;
375  (*this)[+10.8491] = +0.00098;
376  (*this)[+10.9171] = +0.00097;
377  (*this)[+10.9851] = +0.00096;
378  (*this)[+11.0531] = +0.00095;
379  (*this)[+11.1211] = +0.00095;
380  (*this)[+11.1891] = +0.00094;
381  (*this)[+11.2571] = +0.00093;
382  (*this)[+11.3251] = +0.00092;
383  (*this)[+11.3931] = +0.00091;
384  (*this)[+11.4611] = +0.00090;
385  (*this)[+11.5291] = +0.00090;
386  (*this)[+11.5971] = +0.00089;
387  (*this)[+11.6651] = +0.00089;
388  (*this)[+11.7331] = +0.00089;
389  (*this)[+11.8011] = +0.00089;
390  (*this)[+11.8691] = +0.00089;
391  (*this)[+11.9371] = +0.00088;
392  (*this)[+12.0051] = +0.00088;
393  (*this)[+12.0731] = +0.00088;
394  (*this)[+12.1411] = +0.00088;
395  (*this)[+12.2091] = +0.00088;
396  (*this)[+12.2771] = +0.00087;
397  (*this)[+12.3451] = +0.00087;
398  (*this)[+12.4131] = +0.00087;
399  (*this)[+12.4811] = +0.00087;
400  (*this)[+12.5491] = +0.00087;
401  (*this)[+12.6171] = +0.00087;
402  (*this)[+12.6851] = +0.00087;
403  (*this)[+12.7531] = +0.00087;
404  (*this)[+12.8211] = +0.00088;
405  (*this)[+12.8891] = +0.00088;
406  (*this)[+12.9571] = +0.00088;
407  (*this)[+13.0251] = +0.00089;
408  (*this)[+13.0931] = +0.00089;
409  (*this)[+13.1611] = +0.00089;
410  (*this)[+13.2291] = +0.00090;
411  (*this)[+13.2971] = +0.00090;
412  (*this)[+13.3651] = +0.00090;
413  (*this)[+13.4331] = +0.00090;
414  (*this)[+13.5011] = +0.00091;
415  (*this)[+13.5691] = +0.00091;
416  (*this)[+13.6371] = +0.00092;
417  (*this)[+13.7051] = +0.00092;
418  (*this)[+13.7731] = +0.00093;
419  (*this)[+13.8411] = +0.00093;
420  (*this)[+13.9091] = +0.00094;
421  (*this)[+13.9771] = +0.00094;
422  (*this)[+14.0451] = +0.00095;
423  (*this)[+14.1131] = +0.00095;
424  (*this)[+14.1811] = +0.00096;
425  (*this)[+14.2491] = +0.00096;
426  (*this)[+14.3171] = +0.00097;
427  (*this)[+14.3851] = +0.00097;
428  (*this)[+14.4531] = +0.00098;
429  (*this)[+14.5211] = +0.00098;
430  (*this)[+14.5891] = +0.00099;
431  (*this)[+14.6571] = +0.00099;
432  (*this)[+14.7251] = +0.00100;
433  (*this)[+14.7931] = +0.00100;
434  (*this)[+14.8611] = +0.00101;
435  (*this)[+14.9291] = +0.00102;
436  (*this)[+14.9971] = +0.00102;
437  (*this)[+15.0651] = +0.00103;
438  (*this)[+15.1331] = +0.00103;
439  (*this)[+15.2011] = +0.00104;
440  (*this)[+15.2691] = +0.00105;
441  (*this)[+15.3371] = +0.00105;
442  (*this)[+15.4051] = +0.00106;
443  (*this)[+15.4731] = +0.00106;
444  (*this)[+15.5411] = +0.00107;
445  (*this)[+15.6091] = +0.00108;
446  (*this)[+15.6771] = +0.00108;
447  (*this)[+15.7451] = +0.00109;
448  (*this)[+15.8131] = +0.00110;
449  (*this)[+15.8811] = +0.00110;
450  (*this)[+15.9491] = +0.00111;
451  (*this)[+16.0171] = +0.00112;
452  (*this)[+16.0851] = +0.00113;
453  (*this)[+16.1531] = +0.00113;
454  (*this)[+16.2211] = +0.00114;
455  (*this)[+16.2891] = +0.00115;
456  (*this)[+16.3571] = +0.00115;
457  (*this)[+16.4251] = +0.00116;
458  (*this)[+16.4931] = +0.00117;
459  (*this)[+16.5611] = +0.00118;
460  (*this)[+16.6291] = +0.00118;
461  (*this)[+16.6971] = +0.00119;
462  (*this)[+16.7651] = +0.00119;
463  (*this)[+16.8331] = +0.00120;
464  (*this)[+16.9011] = +0.00121;
465  (*this)[+16.9691] = +0.00121;
466  (*this)[+17.0371] = +0.00122;
467  (*this)[+17.1051] = +0.00123;
468  (*this)[+17.1731] = +0.00123;
469  (*this)[+17.2411] = +0.00124;
470  (*this)[+17.3091] = +0.00125;
471  (*this)[+17.3771] = +0.00125;
472  (*this)[+17.4451] = +0.00126;
473  (*this)[+17.5131] = +0.00127;
474  (*this)[+17.5811] = +0.00127;
475  (*this)[+17.6491] = +0.00128;
476  (*this)[+17.7171] = +0.00128;
477  (*this)[+17.7851] = +0.00128;
478  (*this)[+17.8531] = +0.00129;
479  (*this)[+17.9211] = +0.00129;
480  (*this)[+17.9891] = +0.00130;
481  (*this)[+18.0571] = +0.00130;
482  (*this)[+18.1251] = +0.00131;
483  (*this)[+18.1931] = +0.00131;
484  (*this)[+18.2611] = +0.00132;
485  (*this)[+18.3291] = +0.00132;
486  (*this)[+18.3971] = +0.00133;
487  (*this)[+18.4651] = +0.00133;
488  (*this)[+18.5331] = +0.00133;
489  (*this)[+18.6011] = +0.00134;
490  (*this)[+18.6691] = +0.00134;
491  (*this)[+18.7371] = +0.00134;
492  (*this)[+18.8051] = +0.00135;
493  (*this)[+18.8731] = +0.00135;
494  (*this)[+18.9411] = +0.00135;
495  (*this)[+19.0091] = +0.00136;
496  (*this)[+19.0771] = +0.00136;
497  (*this)[+19.1451] = +0.00136;
498  (*this)[+19.2131] = +0.00137;
499  (*this)[+19.2811] = +0.00137;
500  (*this)[+19.3491] = +0.00137;
501  (*this)[+19.4171] = +0.00137;
502  (*this)[+19.4851] = +0.00138;
503  (*this)[+19.5531] = +0.00138;
504  (*this)[+19.6211] = +0.00138;
505  (*this)[+19.6891] = +0.00139;
506  (*this)[+19.7571] = +0.00139;
507  (*this)[+19.8251] = +0.00139;
508  (*this)[+19.8931] = +0.00139;
509  (*this)[+19.9611] = +0.00140;
510  (*this)[+20.0291] = +0.00140;
511  (*this)[+20.0971] = +0.00140;
512  (*this)[+20.1651] = +0.00140;
513  (*this)[+20.2331] = +0.00141;
514  (*this)[+20.3011] = +0.00141;
515  (*this)[+20.3691] = +0.00141;
516  (*this)[+20.4371] = +0.00142;
517  (*this)[+20.5051] = +0.00142;
518  (*this)[+20.5731] = +0.00142;
519  (*this)[+20.6411] = +0.00142;
520  (*this)[+20.7091] = +0.00142;
521  (*this)[+20.7771] = +0.00143;
522  (*this)[+20.8451] = +0.00143;
523  (*this)[+20.9131] = +0.00143;
524  (*this)[+20.9811] = +0.00143;
525  (*this)[+21.0491] = +0.00143;
526  (*this)[+21.1171] = +0.00143;
527  (*this)[+21.1851] = +0.00143;
528  (*this)[+21.2531] = +0.00144;
529  (*this)[+21.3211] = +0.00144;
530  (*this)[+21.3891] = +0.00144;
531  (*this)[+21.4571] = +0.00144;
532  (*this)[+21.5251] = +0.00144;
533  (*this)[+21.5931] = +0.00144;
534  (*this)[+21.6611] = +0.00145;
535  (*this)[+21.7291] = +0.00145;
536  (*this)[+21.7971] = +0.00145;
537  (*this)[+21.8651] = +0.00145;
538  (*this)[+21.9331] = +0.00145;
539  (*this)[+22.0011] = +0.00145;
540  (*this)[+22.0691] = +0.00145;
541  (*this)[+22.1371] = +0.00146;
542  (*this)[+22.2051] = +0.00146;
543  (*this)[+22.2731] = +0.00146;
544  (*this)[+22.3411] = +0.00146;
545  (*this)[+22.4091] = +0.00146;
546  (*this)[+22.4771] = +0.00146;
547  (*this)[+22.5451] = +0.00146;
548  (*this)[+22.6131] = +0.00147;
549  (*this)[+22.6811] = +0.00147;
550  (*this)[+22.7491] = +0.00147;
551  (*this)[+22.8171] = +0.00148;
552  (*this)[+22.8851] = +0.00148;
553  (*this)[+22.9531] = +0.00148;
554  (*this)[+23.0211] = +0.00148;
555  (*this)[+23.0891] = +0.00149;
556  (*this)[+23.1571] = +0.00149;
557  (*this)[+23.2251] = +0.00149;
558  (*this)[+23.2931] = +0.00150;
559  (*this)[+23.3611] = +0.00150;
560  (*this)[+23.4291] = +0.00150;
561  (*this)[+23.4971] = +0.00151;
562  (*this)[+23.5651] = +0.00151;
563  (*this)[+23.6331] = +0.00151;
564  (*this)[+23.7011] = +0.00152;
565  (*this)[+23.7691] = +0.00152;
566  (*this)[+23.8371] = +0.00153;
567  (*this)[+23.9051] = +0.00153;
568  (*this)[+23.9731] = +0.00153;
569  (*this)[+24.0411] = +0.00154;
570  (*this)[+24.1091] = +0.00154;
571  (*this)[+24.1771] = +0.00155;
572  (*this)[+24.2451] = +0.00155;
573  (*this)[+24.3131] = +0.00156;
574  (*this)[+24.3811] = +0.00156;
575  (*this)[+24.4491] = +0.00156;
576  (*this)[+24.5171] = +0.00157;
577  (*this)[+24.5851] = +0.00157;
578  (*this)[+24.6531] = +0.00158;
579  (*this)[+24.7211] = +0.00158;
580  (*this)[+24.7891] = +0.00159;
581  (*this)[+24.8571] = +0.00159;
582  (*this)[+24.9251] = +0.00160;
583  (*this)[+24.9931] = +0.00160;
584  (*this)[+25.0611] = +0.00160;
585  (*this)[+25.1291] = +0.00161;
586  (*this)[+25.1971] = +0.00161;
587  (*this)[+25.2651] = +0.00162;
588  (*this)[+25.3331] = +0.00162;
589  (*this)[+25.4011] = +0.00163;
590  (*this)[+25.4691] = +0.00163;
591  (*this)[+25.5371] = +0.00164;
592  (*this)[+25.6051] = +0.00164;
593  (*this)[+25.6731] = +0.00164;
594  (*this)[+25.7411] = +0.00164;
595  (*this)[+25.8091] = +0.00164;
596  (*this)[+25.8771] = +0.00165;
597  (*this)[+25.9451] = +0.00165;
598  (*this)[+26.0131] = +0.00165;
599  (*this)[+26.0811] = +0.00165;
600  (*this)[+26.1491] = +0.00165;
601  (*this)[+26.2171] = +0.00165;
602  (*this)[+26.2851] = +0.00166;
603  (*this)[+26.3531] = +0.00166;
604  (*this)[+26.4211] = +0.00166;
605  (*this)[+26.4891] = +0.00166;
606  (*this)[+26.5571] = +0.00166;
607  (*this)[+26.6251] = +0.00166;
608  (*this)[+26.6931] = +0.00166;
609  (*this)[+26.7611] = +0.00166;
610  (*this)[+26.8291] = +0.00166;
611  (*this)[+26.8971] = +0.00166;
612  (*this)[+26.9651] = +0.00165;
613  (*this)[+27.0331] = +0.00165;
614  (*this)[+27.1011] = +0.00165;
615  (*this)[+27.1691] = +0.00165;
616  (*this)[+27.2371] = +0.00165;
617  (*this)[+27.3051] = +0.00165;
618  (*this)[+27.3731] = +0.00165;
619  (*this)[+27.4411] = +0.00164;
620  (*this)[+27.5091] = +0.00164;
621  (*this)[+27.5771] = +0.00164;
622  (*this)[+27.6451] = +0.00163;
623  (*this)[+27.7131] = +0.00163;
624  (*this)[+27.7811] = +0.00162;
625  (*this)[+27.8491] = +0.00162;
626  (*this)[+27.9171] = +0.00161;
627  (*this)[+27.9851] = +0.00160;
628  (*this)[+28.0531] = +0.00160;
629  (*this)[+28.1211] = +0.00159;
630  (*this)[+28.1891] = +0.00159;
631  (*this)[+28.2571] = +0.00158;
632  (*this)[+28.3251] = +0.00157;
633  (*this)[+28.3931] = +0.00157;
634  (*this)[+28.4611] = +0.00156;
635  (*this)[+28.5291] = +0.00156;
636  (*this)[+28.5971] = +0.00155;
637  (*this)[+28.6651] = +0.00153;
638  (*this)[+28.7331] = +0.00152;
639  (*this)[+28.8011] = +0.00151;
640  (*this)[+28.8691] = +0.00150;
641  (*this)[+28.9371] = +0.00149;
642  (*this)[+29.0051] = +0.00148;
643  (*this)[+29.0731] = +0.00147;
644  (*this)[+29.1411] = +0.00146;
645  (*this)[+29.2091] = +0.00144;
646  (*this)[+29.2771] = +0.00143;
647  (*this)[+29.3451] = +0.00142;
648  (*this)[+29.4131] = +0.00141;
649  (*this)[+29.4811] = +0.00140;
650  (*this)[+29.5491] = +0.00139;
651  (*this)[+29.6171] = +0.00137;
652  (*this)[+29.6851] = +0.00136;
653  (*this)[+29.7531] = +0.00134;
654  (*this)[+29.8211] = +0.00133;
655  (*this)[+29.8891] = +0.00131;
656  (*this)[+29.9571] = +0.00130;
657  (*this)[+30.0251] = +0.00128;
658  (*this)[+30.0931] = +0.00127;
659  (*this)[+30.1611] = +0.00125;
660  (*this)[+30.2291] = +0.00124;
661  (*this)[+30.2971] = +0.00122;
662  (*this)[+30.3651] = +0.00121;
663  (*this)[+30.4331] = +0.00119;
664  (*this)[+30.5011] = +0.00118;
665  (*this)[+30.5691] = +0.00116;
666  (*this)[+30.6371] = +0.00115;
667  (*this)[+30.7051] = +0.00113;
668  (*this)[+30.7731] = +0.00112;
669  (*this)[+30.8411] = +0.00110;
670  (*this)[+30.9091] = +0.00109;
671  (*this)[+30.9771] = +0.00107;
672  (*this)[+31.0451] = +0.00106;
673  (*this)[+31.1131] = +0.00105;
674  (*this)[+31.1811] = +0.00103;
675  (*this)[+31.2491] = +0.00102;
676  (*this)[+31.3171] = +0.00100;
677  (*this)[+31.3851] = +0.00099;
678  (*this)[+31.4531] = +0.00097;
679  (*this)[+31.5211] = +0.00096;
680  (*this)[+31.5891] = +0.00094;
681  (*this)[+31.6571] = +0.00093;
682  (*this)[+31.7251] = +0.00091;
683  (*this)[+31.7931] = +0.00090;
684  (*this)[+31.8611] = +0.00088;
685  (*this)[+31.9291] = +0.00087;
686  (*this)[+31.9971] = +0.00086;
687  (*this)[+32.0651] = +0.00084;
688  (*this)[+32.1331] = +0.00083;
689  (*this)[+32.2011] = +0.00081;
690  (*this)[+32.2691] = +0.00080;
691  (*this)[+32.3371] = +0.00079;
692  (*this)[+32.4051] = +0.00077;
693  (*this)[+32.4731] = +0.00076;
694  (*this)[+32.5411] = +0.00074;
695  (*this)[+32.6091] = +0.00073;
696  (*this)[+32.6771] = +0.00072;
697  (*this)[+32.7451] = +0.00071;
698  (*this)[+32.8131] = +0.00069;
699  (*this)[+32.8811] = +0.00068;
700  (*this)[+32.9491] = +0.00067;
701  (*this)[+33.0171] = +0.00066;
702  (*this)[+33.0851] = +0.00064;
703  (*this)[+33.1531] = +0.00063;
704  (*this)[+33.2211] = +0.00062;
705  (*this)[+33.2891] = +0.00061;
706  (*this)[+33.3571] = +0.00060;
707  (*this)[+33.4251] = +0.00058;
708  (*this)[+33.4931] = +0.00057;
709  (*this)[+33.5611] = +0.00056;
710  (*this)[+33.6291] = +0.00055;
711  (*this)[+33.6971] = +0.00054;
712  (*this)[+33.7651] = +0.00053;
713  (*this)[+33.8331] = +0.00052;
714  (*this)[+33.9011] = +0.00051;
715  (*this)[+33.9691] = +0.00050;
716  (*this)[+34.0371] = +0.00049;
717  (*this)[+34.1051] = +0.00048;
718  (*this)[+34.1731] = +0.00047;
719  (*this)[+34.2411] = +0.00046;
720  (*this)[+34.3091] = +0.00045;
721  (*this)[+34.3771] = +0.00044;
722  (*this)[+34.4451] = +0.00043;
723  (*this)[+34.5131] = +0.00042;
724  (*this)[+34.5811] = +0.00042;
725  (*this)[+34.6491] = +0.00041;
726  (*this)[+34.7171] = +0.00040;
727  (*this)[+34.7851] = +0.00039;
728  (*this)[+34.8531] = +0.00039;
729  (*this)[+34.9211] = +0.00038;
730  (*this)[+34.9891] = +0.00037;
731  (*this)[+35.0571] = +0.00036;
732  (*this)[+35.1251] = +0.00036;
733  (*this)[+35.1931] = +0.00035;
734  (*this)[+35.2611] = +0.00034;
735  (*this)[+35.3291] = +0.00034;
736  (*this)[+35.3971] = +0.00033;
737  (*this)[+35.4651] = +0.00032;
738  (*this)[+35.5331] = +0.00031;
739  (*this)[+35.6011] = +0.00031;
740  (*this)[+35.6691] = +0.00030;
741  (*this)[+35.7371] = +0.00030;
742  (*this)[+35.8051] = +0.00029;
743  (*this)[+35.8731] = +0.00029;
744  (*this)[+35.9411] = +0.00028;
745  (*this)[+36.0091] = +0.00028;
746  (*this)[+36.0771] = +0.00027;
747  (*this)[+36.1451] = +0.00026;
748  (*this)[+36.2131] = +0.00026;
749  (*this)[+36.2811] = +0.00025;
750  (*this)[+36.3491] = +0.00025;
751  (*this)[+36.4171] = +0.00024;
752  (*this)[+36.4851] = +0.00024;
753  (*this)[+36.5531] = +0.00023;
754  (*this)[+36.6211] = +0.00023;
755  (*this)[+36.6891] = +0.00022;
756  (*this)[+36.7571] = +0.00022;
757  (*this)[+36.8251] = +0.00021;
758  (*this)[+36.8931] = +0.00021;
759  (*this)[+36.9611] = +0.00021;
760  (*this)[+37.0291] = +0.00020;
761  (*this)[+37.0971] = +0.00020;
762  (*this)[+37.1651] = +0.00019;
763  (*this)[+37.2331] = +0.00019;
764  (*this)[+37.3011] = +0.00019;
765  (*this)[+37.3691] = +0.00018;
766  (*this)[+37.4371] = +0.00018;
767  (*this)[+37.5051] = +0.00017;
768  (*this)[+37.5731] = +0.00017;
769  (*this)[+37.6411] = +0.00017;
770  (*this)[+37.7091] = +0.00016;
771  (*this)[+37.7771] = +0.00016;
772  (*this)[+37.8451] = +0.00016;
773  (*this)[+37.9131] = +0.00015;
774  (*this)[+37.9811] = +0.00015;
775  (*this)[+38.0491] = +0.00015;
776  (*this)[+38.1171] = +0.00014;
777  (*this)[+38.1851] = +0.00014;
778  (*this)[+38.2531] = +0.00014;
779  (*this)[+38.3211] = +0.00013;
780  (*this)[+38.3891] = +0.00013;
781  (*this)[+38.4571] = +0.00013;
782  (*this)[+38.5251] = +0.00013;
783  (*this)[+38.5931] = +0.00012;
784  (*this)[+38.6611] = +0.00012;
785  (*this)[+38.7291] = +0.00012;
786  (*this)[+38.7971] = +0.00012;
787  (*this)[+38.8651] = +0.00012;
788  (*this)[+38.9331] = +0.00011;
789  (*this)[+39.0011] = +0.00011;
790  (*this)[+39.0691] = +0.00011;
791  (*this)[+39.1371] = +0.00011;
792  (*this)[+39.2051] = +0.00011;
793  (*this)[+39.2731] = +0.00010;
794  (*this)[+39.3411] = +0.00010;
795  (*this)[+39.4091] = +0.00010;
796  (*this)[+39.4771] = +0.00010;
797  (*this)[+39.5451] = +0.00010;
798  (*this)[+39.6131] = +0.00010;
799  (*this)[+39.6811] = +0.00009;
800  (*this)[+39.7491] = +0.00009;
801  (*this)[+39.8171] = +0.00009;
802  (*this)[+39.8851] = +0.00009;
803  (*this)[+39.9531] = +0.00009;
804  (*this)[+40.0211] = +0.00009;
805  (*this)[+40.0891] = +0.00009;
806  (*this)[+40.1571] = +0.00009;
807  (*this)[+40.2251] = +0.00008;
808  (*this)[+40.2931] = +0.00008;
809  (*this)[+40.3611] = +0.00008;
810  (*this)[+40.4291] = +0.00008;
811  (*this)[+40.4971] = +0.00008;
812  (*this)[+40.5651] = +0.00008;
813  (*this)[+40.6331] = +0.00008;
814  (*this)[+40.7011] = +0.00008;
815  (*this)[+40.7691] = +0.00008;
816  (*this)[+40.8371] = +0.00007;
817  (*this)[+40.9051] = +0.00007;
818  (*this)[+40.9731] = +0.00007;
819  (*this)[+41.0411] = +0.00007;
820  (*this)[+41.1091] = +0.00007;
821  (*this)[+41.1771] = +0.00007;
822  (*this)[+41.2451] = +0.00007;
823  (*this)[+41.3131] = +0.00007;
824  (*this)[+41.3811] = +0.00007;
825  (*this)[+41.4491] = +0.00007;
826  (*this)[+41.5171] = +0.00007;
827  (*this)[+41.5851] = +0.00006;
828  (*this)[+41.6531] = +0.00006;
829  (*this)[+41.7211] = +0.00006;
830  (*this)[+41.7891] = +0.00006;
831  (*this)[+41.8571] = +0.00006;
832  (*this)[+41.9251] = +0.00006;
833  (*this)[+41.9931] = +0.00006;
834  (*this)[+42.0611] = +0.00006;
835  (*this)[+42.1291] = +0.00006;
836  (*this)[+42.1971] = +0.00006;
837  (*this)[+42.2651] = +0.00006;
838  (*this)[+42.3331] = +0.00006;
839  (*this)[+42.4011] = +0.00005;
840  (*this)[+42.4691] = +0.00005;
841  (*this)[+42.5371] = +0.00005;
842  (*this)[+42.6051] = +0.00005;
843  (*this)[+42.6731] = +0.00005;
844  (*this)[+42.7411] = +0.00005;
845  (*this)[+42.8091] = +0.00005;
846  (*this)[+42.8771] = +0.00005;
847  (*this)[+42.9451] = +0.00005;
848  (*this)[+43.0131] = +0.00005;
849  (*this)[+43.0811] = +0.00005;
850  (*this)[+43.1491] = +0.00005;
851  (*this)[+43.2171] = +0.00005;
852  (*this)[+43.2851] = +0.00005;
853  (*this)[+43.3531] = +0.00005;
854  (*this)[+43.4211] = +0.00005;
855  (*this)[+43.4891] = +0.00005;
856  (*this)[+43.5571] = +0.00005;
857  (*this)[+43.6251] = +0.00005;
858  (*this)[+43.6931] = +0.00005;
859  (*this)[+43.7611] = +0.00005;
860  (*this)[+43.8291] = +0.00005;
861  (*this)[+43.8971] = +0.00005;
862  (*this)[+43.9651] = +0.00005;
863  (*this)[+44.0331] = +0.00005;
864  (*this)[+44.1011] = +0.00005;
865  (*this)[+44.1691] = +0.00005;
866  (*this)[+44.2371] = +0.00005;
867  (*this)[+44.3051] = +0.00005;
868  (*this)[+44.3731] = +0.00005;
869  (*this)[+44.4411] = +0.00005;
870  (*this)[+44.5091] = +0.00005;
871  (*this)[+44.5771] = +0.00005;
872  (*this)[+44.6451] = +0.00005;
873  (*this)[+44.7131] = +0.00005;
874  (*this)[+44.7811] = +0.00005;
875  (*this)[+44.8491] = +0.00005;
876  (*this)[+44.9171] = +0.00005;
877  (*this)[+44.9851] = +0.00005;
878  (*this)[+45.0531] = +0.00005;
879  (*this)[+45.1211] = +0.00005;
880  (*this)[+45.1891] = +0.00005;
881  (*this)[+45.2571] = +0.00005;
882  (*this)[+45.3251] = +0.00005;
883  (*this)[+45.3931] = +0.00005;
884  (*this)[+45.4611] = +0.00005;
885  (*this)[+45.5291] = +0.00005;
886  (*this)[+45.5971] = +0.00005;
887  (*this)[+45.6651] = +0.00005;
888  (*this)[+45.7331] = +0.00005;
889  (*this)[+45.8011] = +0.00005;
890  (*this)[+45.8691] = +0.00005;
891  (*this)[+45.9371] = +0.00005;
892  (*this)[+46.0051] = +0.00005;
893  (*this)[+46.0731] = +0.00005;
894  (*this)[+46.1411] = +0.00005;
895  (*this)[+46.2091] = +0.00005;
896  (*this)[+46.2771] = +0.00005;
897  (*this)[+46.3451] = +0.00005;
898  (*this)[+46.4131] = +0.00005;
899  (*this)[+46.4811] = +0.00005;
900  (*this)[+46.5491] = +0.00005;
901  (*this)[+46.6171] = +0.00005;
902  (*this)[+46.6851] = +0.00005;
903  (*this)[+46.7531] = +0.00005;
904  (*this)[+46.8211] = +0.00005;
905  (*this)[+46.8891] = +0.00005;
906  (*this)[+46.9571] = +0.00005;
907  (*this)[+47.0251] = +0.00005;
908  (*this)[+47.0931] = +0.00005;
909  (*this)[+47.1611] = +0.00005;
910  (*this)[+47.2291] = +0.00005;
911  (*this)[+47.2971] = +0.00005;
912  (*this)[+47.3651] = +0.00005;
913  (*this)[+47.4331] = +0.00005;
914  (*this)[+47.5011] = +0.00005;
915  (*this)[+47.5691] = +0.00005;
916  (*this)[+47.6371] = +0.00005;
917  (*this)[+47.7051] = +0.00005;
918  (*this)[+47.7731] = +0.00005;
919  (*this)[+47.8411] = +0.00005;
920  (*this)[+47.9091] = +0.00005;
921  (*this)[+47.9771] = +0.00005;
922  (*this)[+48.0451] = +0.00005;
923  (*this)[+48.1131] = +0.00005;
924  (*this)[+48.1811] = +0.00005;
925  (*this)[+48.2491] = +0.00005;
926  (*this)[+48.3171] = +0.00005;
927  (*this)[+48.3851] = +0.00005;
928  (*this)[+48.4531] = +0.00005;
929  (*this)[+48.5211] = +0.00005;
930  (*this)[+48.5891] = +0.00005;
931  (*this)[+48.6571] = +0.00005;
932  (*this)[+48.7251] = +0.00005;
933  (*this)[+48.7931] = +0.00005;
934  (*this)[+48.8611] = +0.00005;
935  (*this)[+48.9291] = +0.00005;
936  (*this)[+48.9971] = +0.00005;
937  (*this)[+49.0651] = +0.00005;
938  (*this)[+49.1331] = +0.00005;
939  (*this)[+49.2011] = +0.00005;
940  (*this)[+49.2691] = +0.00005;
941  (*this)[+49.3371] = +0.00005;
942  (*this)[+49.4051] = +0.00005;
943  (*this)[+49.4731] = +0.00005;
944  (*this)[+49.5411] = +0.00005;
945  (*this)[+49.6091] = +0.00005;
946  (*this)[+49.6771] = +0.00005;
947  (*this)[+49.7451] = +0.00005;
948  (*this)[+49.8131] = +0.00005;
949  (*this)[+49.8811] = +0.00005;
950  (*this)[+49.9491] = +0.00005;
951  (*this)[+50.0171] = +0.00005;
952  (*this)[+50.0851] = +0.00005;
953  (*this)[+50.1531] = +0.00005;
954  (*this)[+50.2211] = +0.00005;
955  (*this)[+50.2891] = +0.00005;
956  (*this)[+50.3571] = +0.00005;
957  (*this)[+50.4251] = +0.00005;
958  (*this)[+50.4931] = +0.00005;
959  (*this)[+50.5611] = +0.00005;
960  (*this)[+50.6291] = +0.00005;
961  (*this)[+50.6971] = +0.00005;
962  (*this)[+50.7651] = +0.00005;
963  (*this)[+50.8331] = +0.00005;
964  (*this)[+50.9011] = +0.00005;
965  (*this)[+50.9691] = +0.00005;
966  (*this)[+51.0371] = +0.00005;
967  (*this)[+51.1051] = +0.00005;
968  (*this)[+51.1731] = +0.00005;
969  (*this)[+51.2411] = +0.00005;
970  (*this)[+51.3091] = +0.00005;
971  (*this)[+51.3771] = +0.00005;
972  (*this)[+51.4451] = +0.00005;
973  (*this)[+51.5131] = +0.00004;
974  (*this)[+51.5811] = +0.00004;
975  (*this)[+51.6491] = +0.00004;
976  (*this)[+51.7171] = +0.00004;
977  (*this)[+51.7851] = +0.00004;
978  (*this)[+51.8531] = +0.00004;
979  (*this)[+51.9211] = +0.00004;
980  (*this)[+51.9891] = +0.00004;
981  (*this)[+52.0571] = +0.00004;
982  (*this)[+52.1251] = +0.00004;
983  (*this)[+52.1931] = +0.00004;
984  (*this)[+52.2611] = +0.00004;
985  (*this)[+52.3291] = +0.00004;
986  (*this)[+52.3971] = +0.00004;
987  (*this)[+52.4651] = +0.00004;
988  (*this)[+52.5331] = +0.00004;
989  (*this)[+52.6011] = +0.00004;
990  (*this)[+52.6691] = +0.00004;
991  (*this)[+52.7371] = +0.00004;
992  (*this)[+52.8051] = +0.00004;
993  (*this)[+52.8731] = +0.00004;
994  (*this)[+52.9411] = +0.00004;
995  (*this)[+53.0091] = +0.00004;
996  (*this)[+53.0771] = +0.00004;
997  (*this)[+53.1451] = +0.00004;
998  (*this)[+53.2131] = +0.00004;
999  (*this)[+53.2811] = +0.00005;
1000  (*this)[+53.3491] = +0.00005;
1001  (*this)[+53.4171] = +0.00005;
1002  (*this)[+53.4851] = +0.00005;
1003  (*this)[+53.5531] = +0.00005;
1004  (*this)[+53.6211] = +0.00005;
1005  (*this)[+53.6891] = +0.00005;
1006  (*this)[+53.7571] = +0.00004;
1007  (*this)[+53.8251] = +0.00004;
1008  (*this)[+53.8931] = +0.00004;
1009  (*this)[+53.9611] = +0.00004;
1010  (*this)[+54.0291] = +0.00004;
1011  (*this)[+54.0971] = +0.00004;
1012  (*this)[+54.1651] = +0.00004;
1013  (*this)[+54.2331] = +0.00004;
1014  (*this)[+54.3011] = +0.00004;
1015  (*this)[+54.3691] = +0.00004;
1016  (*this)[+54.4371] = +0.00004;
1017  (*this)[+54.5051] = +0.00004;
1018  (*this)[+54.5731] = +0.00004;
1019  (*this)[+54.6411] = +0.00004;
1020  (*this)[+54.7091] = +0.00004;
1021  (*this)[+54.7771] = +0.00004;
1022  (*this)[+54.8451] = +0.00004;
1023  (*this)[+54.9131] = +0.00004;
1024  (*this)[+54.9811] = +0.00004;
1025  (*this)[+55.0491] = +0.00004;
1026  (*this)[+55.1171] = +0.00004;
1027  (*this)[+55.1851] = +0.00004;
1028  (*this)[+55.2531] = +0.00004;
1029  (*this)[+55.3211] = +0.00004;
1030  (*this)[+55.3891] = +0.00004;
1031  (*this)[+55.4571] = +0.00004;
1032  (*this)[+55.5251] = +0.00004;
1033  (*this)[+55.5931] = +0.00004;
1034  (*this)[+55.6611] = +0.00003;
1035  (*this)[+55.7291] = +0.00003;
1036  (*this)[+55.7971] = +0.00003;
1037  (*this)[+55.8651] = +0.00003;
1038  (*this)[+55.9331] = +0.00002;
1039  (*this)[+56.0011] = +0.00002;
1040  (*this)[+56.0691] = +0.00002;
1041  (*this)[+56.1371] = +0.00002;
1042  (*this)[+56.2051] = +0.00001;
1043  (*this)[+56.2731] = +0.00001;
1044  (*this)[+56.3411] = +0.00001;
1045  (*this)[+56.4091] = +0.00001;
1046  (*this)[+56.4771] = +0.00000;
1047  (*this)[+56.5451] = +0.00000;
1048 
1049  compile();
1050  };
1051  };
1052 
1053 
1054  /**
1055  * Implementation of probability distibution of the PMT transition times based on measurements.
1056  */
1057  struct JTransitionTimeProbability38 :
1058  public JTransitionTimeProbability_t
1059  {
1060  /**
1061  * Default constructor.
1062  */
1063  JTransitionTimeProbability38()
1064  {
1065  // produced by JLegolas.cc
1066  (*this)[ -20.00] = 0.000000;
1067  (*this)[ -19.75] = 0.000000;
1068  (*this)[ -19.50] = 0.000000;
1069  (*this)[ -19.25] = 0.000002;
1070  (*this)[ -19.00] = 0.000004;
1071  (*this)[ -18.75] = 0.000006;
1072  (*this)[ -18.50] = 0.000009;
1073  (*this)[ -18.25] = 0.000011;
1074  (*this)[ -18.00] = 0.000015;
1075  (*this)[ -17.75] = 0.000018;
1076  (*this)[ -17.50] = 0.000023;
1077  (*this)[ -17.25] = 0.000028;
1078  (*this)[ -17.00] = 0.000032;
1079  (*this)[ -16.75] = 0.000035;
1080  (*this)[ -16.50] = 0.000039;
1081  (*this)[ -16.25] = 0.000043;
1082  (*this)[ -16.00] = 0.000045;
1083  (*this)[ -15.75] = 0.000047;
1084  (*this)[ -15.50] = 0.000047;
1085  (*this)[ -15.25] = 0.000048;
1086  (*this)[ -15.00] = 0.000047;
1087  (*this)[ -14.75] = 0.000045;
1088  (*this)[ -14.50] = 0.000044;
1089  (*this)[ -14.25] = 0.000044;
1090  (*this)[ -14.00] = 0.000041;
1091  (*this)[ -13.75] = 0.000039;
1092  (*this)[ -13.50] = 0.000038;
1093  (*this)[ -13.25] = 0.000037;
1094  (*this)[ -13.00] = 0.000034;
1095  (*this)[ -12.75] = 0.000032;
1096  (*this)[ -12.50] = 0.000030;
1097  (*this)[ -12.25] = 0.000028;
1098  (*this)[ -12.00] = 0.000025;
1099  (*this)[ -11.75] = 0.000023;
1100  (*this)[ -11.50] = 0.000022;
1101  (*this)[ -11.25] = 0.000021;
1102  (*this)[ -11.00] = 0.000020;
1103  (*this)[ -10.75] = 0.000019;
1104  (*this)[ -10.50] = 0.000020;
1105  (*this)[ -10.25] = 0.000020;
1106  (*this)[ -10.00] = 0.000022;
1107  (*this)[ -9.75] = 0.000024;
1108  (*this)[ -9.50] = 0.000027;
1109  (*this)[ -9.25] = 0.000030;
1110  (*this)[ -9.00] = 0.000035;
1111  (*this)[ -8.75] = 0.000039;
1112  (*this)[ -8.50] = 0.000045;
1113  (*this)[ -8.25] = 0.000050;
1114  (*this)[ -8.00] = 0.000058;
1115  (*this)[ -7.75] = 0.000065;
1116  (*this)[ -7.50] = 0.000074;
1117  (*this)[ -7.25] = 0.000084;
1118  (*this)[ -7.00] = 0.000096;
1119  (*this)[ -6.75] = 0.000109;
1120  (*this)[ -6.50] = 0.000125;
1121  (*this)[ -6.25] = 0.000141;
1122  (*this)[ -6.00] = 0.000167;
1123  (*this)[ -5.75] = 0.000194;
1124  (*this)[ -5.50] = 0.000243;
1125  (*this)[ -5.25] = 0.000292;
1126  (*this)[ -5.00] = 0.000394;
1127  (*this)[ -4.75] = 0.000500;
1128  (*this)[ -4.50] = 0.000778;
1129  (*this)[ -4.25] = 0.001069;
1130  (*this)[ -4.00] = 0.001918;
1131  (*this)[ -3.75] = 0.002797;
1132  (*this)[ -3.50] = 0.005052;
1133  (*this)[ -3.25] = 0.007357;
1134  (*this)[ -3.00] = 0.012074;
1135  (*this)[ -2.75] = 0.016857;
1136  (*this)[ -2.50] = 0.024827;
1137  (*this)[ -2.25] = 0.032883;
1138  (*this)[ -2.00] = 0.045017;
1139  (*this)[ -1.75] = 0.057241;
1140  (*this)[ -1.50] = 0.073679;
1141  (*this)[ -1.25] = 0.090146;
1142  (*this)[ -1.00] = 0.107977;
1143  (*this)[ -0.75] = 0.125687;
1144  (*this)[ -0.50] = 0.137677;
1145  (*this)[ -0.25] = 0.149420;
1146  (*this)[ +0.00] = 0.149505;
1147  (*this)[ +0.25] = 0.149327;
1148  (*this)[ +0.50] = 0.136745;
1149  (*this)[ +0.75] = 0.124047;
1150  (*this)[ +1.00] = 0.105908;
1151  (*this)[ +1.25] = 0.087797;
1152  (*this)[ +1.50] = 0.071017;
1153  (*this)[ +1.75] = 0.054352;
1154  (*this)[ +2.00] = 0.043056;
1155  (*this)[ +2.25] = 0.031857;
1156  (*this)[ +2.50] = 0.025196;
1157  (*this)[ +2.75] = 0.018592;
1158  (*this)[ +3.00] = 0.014719;
1159  (*this)[ +3.25] = 0.010873;
1160  (*this)[ +3.50] = 0.008306;
1161  (*this)[ +3.75] = 0.005762;
1162  (*this)[ +4.00] = 0.004358;
1163  (*this)[ +4.25] = 0.002972;
1164  (*this)[ +4.50] = 0.002403;
1165  (*this)[ +4.75] = 0.001842;
1166  (*this)[ +5.00] = 0.001643;
1167  (*this)[ +5.25] = 0.001446;
1168  (*this)[ +5.50] = 0.001319;
1169  (*this)[ +5.75] = 0.001192;
1170  (*this)[ +6.00] = 0.001064;
1171  (*this)[ +6.25] = 0.000937;
1172  (*this)[ +6.50] = 0.000826;
1173  (*this)[ +6.75] = 0.000716;
1174  (*this)[ +7.00] = 0.000636;
1175  (*this)[ +7.25] = 0.000557;
1176  (*this)[ +7.50] = 0.000503;
1177  (*this)[ +7.75] = 0.000448;
1178  (*this)[ +8.00] = 0.000412;
1179  (*this)[ +8.25] = 0.000376;
1180  (*this)[ +8.50] = 0.000351;
1181  (*this)[ +8.75] = 0.000327;
1182  (*this)[ +9.00] = 0.000309;
1183  (*this)[ +9.25] = 0.000291;
1184  (*this)[ +9.50] = 0.000280;
1185  (*this)[ +9.75] = 0.000268;
1186  (*this)[ +10.00] = 0.000260;
1187  (*this)[ +10.25] = 0.000252;
1188  (*this)[ +10.50] = 0.000248;
1189  (*this)[ +10.75] = 0.000243;
1190  (*this)[ +11.00] = 0.000240;
1191  (*this)[ +11.25] = 0.000237;
1192  (*this)[ +11.50] = 0.000236;
1193  (*this)[ +11.75] = 0.000236;
1194  (*this)[ +12.00] = 0.000237;
1195  (*this)[ +12.25] = 0.000237;
1196  (*this)[ +12.50] = 0.000241;
1197  (*this)[ +12.75] = 0.000244;
1198  (*this)[ +13.00] = 0.000249;
1199  (*this)[ +13.25] = 0.000253;
1200  (*this)[ +13.50] = 0.000259;
1201  (*this)[ +13.75] = 0.000265;
1202  (*this)[ +14.00] = 0.000270;
1203  (*this)[ +14.25] = 0.000275;
1204  (*this)[ +14.50] = 0.000281;
1205  (*this)[ +14.75] = 0.000287;
1206  (*this)[ +15.00] = 0.000290;
1207  (*this)[ +15.25] = 0.000294;
1208  (*this)[ +15.50] = 0.000297;
1209  (*this)[ +15.75] = 0.000300;
1210  (*this)[ +16.00] = 0.000302;
1211  (*this)[ +16.25] = 0.000304;
1212  (*this)[ +16.50] = 0.000307;
1213  (*this)[ +16.75] = 0.000309;
1214  (*this)[ +17.00] = 0.000311;
1215  (*this)[ +17.25] = 0.000313;
1216  (*this)[ +17.50] = 0.000315;
1217  (*this)[ +17.75] = 0.000317;
1218  (*this)[ +18.00] = 0.000317;
1219  (*this)[ +18.25] = 0.000317;
1220  (*this)[ +18.50] = 0.000318;
1221  (*this)[ +18.75] = 0.000318;
1222  (*this)[ +19.00] = 0.000316;
1223  (*this)[ +19.25] = 0.000314;
1224  (*this)[ +19.50] = 0.000312;
1225  (*this)[ +19.75] = 0.000311;
1226  (*this)[ +20.00] = 0.000310;
1227  (*this)[ +20.25] = 0.000309;
1228  (*this)[ +20.50] = 0.000311;
1229  (*this)[ +20.75] = 0.000312;
1230  (*this)[ +21.00] = 0.000314;
1231  (*this)[ +21.25] = 0.000317;
1232  (*this)[ +21.50] = 0.000321;
1233  (*this)[ +21.75] = 0.000326;
1234  (*this)[ +22.00] = 0.000330;
1235  (*this)[ +22.25] = 0.000335;
1236  (*this)[ +22.50] = 0.000344;
1237  (*this)[ +22.75] = 0.000352;
1238  (*this)[ +23.00] = 0.000360;
1239  (*this)[ +23.25] = 0.000368;
1240  (*this)[ +23.50] = 0.000379;
1241  (*this)[ +23.75] = 0.000390;
1242  (*this)[ +24.00] = 0.000403;
1243  (*this)[ +24.25] = 0.000416;
1244  (*this)[ +24.50] = 0.000432;
1245  (*this)[ +24.75] = 0.000448;
1246  (*this)[ +25.00] = 0.000463;
1247  (*this)[ +25.25] = 0.000479;
1248  (*this)[ +25.50] = 0.000495;
1249  (*this)[ +25.75] = 0.000511;
1250  (*this)[ +26.00] = 0.000525;
1251  (*this)[ +26.25] = 0.000538;
1252  (*this)[ +26.50] = 0.000549;
1253  (*this)[ +26.75] = 0.000559;
1254  (*this)[ +27.00] = 0.000563;
1255  (*this)[ +27.25] = 0.000567;
1256  (*this)[ +27.50] = 0.000566;
1257  (*this)[ +27.75] = 0.000565;
1258  (*this)[ +28.00] = 0.000560;
1259  (*this)[ +28.25] = 0.000555;
1260  (*this)[ +28.50] = 0.000548;
1261  (*this)[ +28.75] = 0.000541;
1262  (*this)[ +29.00] = 0.000530;
1263  (*this)[ +29.25] = 0.000519;
1264  (*this)[ +29.50] = 0.000508;
1265  (*this)[ +29.75] = 0.000496;
1266  (*this)[ +30.00] = 0.000482;
1267  (*this)[ +30.25] = 0.000467;
1268  (*this)[ +30.50] = 0.000452;
1269  (*this)[ +30.75] = 0.000437;
1270  (*this)[ +31.00] = 0.000417;
1271  (*this)[ +31.25] = 0.000397;
1272  (*this)[ +31.50] = 0.000379;
1273  (*this)[ +31.75] = 0.000360;
1274  (*this)[ +32.00] = 0.000341;
1275  (*this)[ +32.25] = 0.000322;
1276  (*this)[ +32.50] = 0.000307;
1277  (*this)[ +32.75] = 0.000292;
1278  (*this)[ +33.00] = 0.000276;
1279  (*this)[ +33.25] = 0.000260;
1280  (*this)[ +33.50] = 0.000247;
1281  (*this)[ +33.75] = 0.000235;
1282  (*this)[ +34.00] = 0.000222;
1283  (*this)[ +34.25] = 0.000209;
1284  (*this)[ +34.50] = 0.000197;
1285  (*this)[ +34.75] = 0.000186;
1286  (*this)[ +35.00] = 0.000173;
1287  (*this)[ +35.25] = 0.000161;
1288  (*this)[ +35.50] = 0.000150;
1289  (*this)[ +35.75] = 0.000140;
1290  (*this)[ +36.00] = 0.000130;
1291  (*this)[ +36.25] = 0.000121;
1292  (*this)[ +36.50] = 0.000114;
1293  (*this)[ +36.75] = 0.000108;
1294  (*this)[ +37.00] = 0.000102;
1295  (*this)[ +37.25] = 0.000096;
1296  (*this)[ +37.50] = 0.000092;
1297  (*this)[ +37.75] = 0.000087;
1298  (*this)[ +38.00] = 0.000084;
1299  (*this)[ +38.25] = 0.000080;
1300  (*this)[ +38.50] = 0.000077;
1301  (*this)[ +38.75] = 0.000073;
1302  (*this)[ +39.00] = 0.000069;
1303  (*this)[ +39.25] = 0.000066;
1304  (*this)[ +39.50] = 0.000063;
1305  (*this)[ +39.75] = 0.000060;
1306  (*this)[ +40.00] = 0.000057;
1307  (*this)[ +40.25] = 0.000055;
1308  (*this)[ +40.50] = 0.000054;
1309  (*this)[ +40.75] = 0.000053;
1310  (*this)[ +41.00] = 0.000053;
1311  (*this)[ +41.25] = 0.000052;
1312  (*this)[ +41.50] = 0.000052;
1313  (*this)[ +41.75] = 0.000052;
1314  (*this)[ +42.00] = 0.000052;
1315  (*this)[ +42.25] = 0.000052;
1316  (*this)[ +42.50] = 0.000051;
1317  (*this)[ +42.75] = 0.000051;
1318  (*this)[ +43.00] = 0.000050;
1319  (*this)[ +43.25] = 0.000049;
1320  (*this)[ +43.50] = 0.000048;
1321  (*this)[ +43.75] = 0.000047;
1322  (*this)[ +44.00] = 0.000046;
1323  (*this)[ +44.25] = 0.000045;
1324  (*this)[ +44.50] = 0.000045;
1325  (*this)[ +44.75] = 0.000045;
1326  (*this)[ +45.00] = 0.000045;
1327  (*this)[ +45.25] = 0.000045;
1328  (*this)[ +45.50] = 0.000046;
1329  (*this)[ +45.75] = 0.000046;
1330  (*this)[ +46.00] = 0.000046;
1331  (*this)[ +46.25] = 0.000046;
1332  (*this)[ +46.50] = 0.000046;
1333  (*this)[ +46.75] = 0.000046;
1334  (*this)[ +47.00] = 0.000045;
1335  (*this)[ +47.25] = 0.000044;
1336  (*this)[ +47.50] = 0.000043;
1337  (*this)[ +47.75] = 0.000042;
1338  (*this)[ +48.00] = 0.000041;
1339  (*this)[ +48.25] = 0.000041;
1340  (*this)[ +48.50] = 0.000040;
1341  (*this)[ +48.75] = 0.000040;
1342  (*this)[ +49.00] = 0.000040;
1343  (*this)[ +49.25] = 0.000040;
1344  (*this)[ +49.50] = 0.000040;
1345  (*this)[ +49.75] = 0.000040;
1346  (*this)[ +50.00] = 0.000040;
1347  (*this)[ +50.25] = 0.000039;
1348  (*this)[ +50.50] = 0.000039;
1349  (*this)[ +50.75] = 0.000039;
1350  (*this)[ +51.00] = 0.000038;
1351  (*this)[ +51.25] = 0.000038;
1352  (*this)[ +51.50] = 0.000037;
1353  (*this)[ +51.75] = 0.000037;
1354  (*this)[ +52.00] = 0.000036;
1355  (*this)[ +52.25] = 0.000036;
1356  (*this)[ +52.50] = 0.000036;
1357  (*this)[ +52.75] = 0.000036;
1358  (*this)[ +53.00] = 0.000036;
1359  (*this)[ +53.25] = 0.000036;
1360  (*this)[ +53.50] = 0.000036;
1361  (*this)[ +53.75] = 0.000036;
1362  (*this)[ +54.00] = 0.000035;
1363  (*this)[ +54.25] = 0.000035;
1364  (*this)[ +54.50] = 0.000035;
1365  (*this)[ +54.75] = 0.000034;
1366  (*this)[ +55.00] = 0.000033;
1367  (*this)[ +55.25] = 0.000033;
1368  (*this)[ +55.50] = 0.000032;
1369  (*this)[ +55.75] = 0.000031;
1370  (*this)[ +56.00] = 0.000031;
1371  (*this)[ +56.25] = 0.000030;
1372  (*this)[ +56.50] = 0.000029;
1373  (*this)[ +56.75] = 0.000029;
1374  (*this)[ +57.00] = 0.000028;
1375  (*this)[ +57.25] = 0.000028;
1376  (*this)[ +57.50] = 0.000027;
1377  (*this)[ +57.75] = 0.000027;
1378  (*this)[ +58.00] = 0.000026;
1379  (*this)[ +58.25] = 0.000026;
1380  (*this)[ +58.50] = 0.000025;
1381  (*this)[ +58.75] = 0.000024;
1382  (*this)[ +59.00] = 0.000023;
1383  (*this)[ +59.25] = 0.000022;
1384  (*this)[ +59.50] = 0.000021;
1385  (*this)[ +59.75] = 0.000021;
1386  (*this)[ +60.00] = 0.000020;
1387  (*this)[ +60.25] = 0.000019;
1388  (*this)[ +60.50] = 0.000018;
1389  (*this)[ +60.75] = 0.000018;
1390  (*this)[ +61.00] = 0.000017;
1391  (*this)[ +61.25] = 0.000016;
1392  (*this)[ +61.50] = 0.000016;
1393  (*this)[ +61.75] = 0.000015;
1394  (*this)[ +62.00] = 0.000015;
1395  (*this)[ +62.25] = 0.000014;
1396  (*this)[ +62.50] = 0.000013;
1397  (*this)[ +62.75] = 0.000013;
1398  (*this)[ +63.00] = 0.000012;
1399  (*this)[ +63.25] = 0.000012;
1400  (*this)[ +63.50] = 0.000011;
1401  (*this)[ +63.75] = 0.000010;
1402  (*this)[ +64.00] = 0.000010;
1403  (*this)[ +64.25] = 0.000009;
1404  (*this)[ +64.50] = 0.000009;
1405  (*this)[ +64.75] = 0.000009;
1406  (*this)[ +65.00] = 0.000008;
1407  (*this)[ +65.25] = 0.000008;
1408  (*this)[ +65.50] = 0.000007;
1409  (*this)[ +65.75] = 0.000007;
1410  (*this)[ +66.00] = 0.000007;
1411  (*this)[ +66.25] = 0.000006;
1412  (*this)[ +66.50] = 0.000006;
1413  (*this)[ +66.75] = 0.000006;
1414  (*this)[ +67.00] = 0.000005;
1415  (*this)[ +67.25] = 0.000005;
1416  (*this)[ +67.50] = 0.000005;
1417  (*this)[ +67.75] = 0.000005;
1418  (*this)[ +68.00] = 0.000005;
1419  (*this)[ +68.25] = 0.000004;
1420  (*this)[ +68.50] = 0.000004;
1421  (*this)[ +68.75] = 0.000004;
1422  (*this)[ +69.00] = 0.000004;
1423  (*this)[ +69.25] = 0.000004;
1424  (*this)[ +69.50] = 0.000004;
1425  (*this)[ +69.75] = 0.000004;
1426  (*this)[ +70.00] = 0.000004;
1427  (*this)[ +70.25] = 0.000003;
1428  (*this)[ +70.50] = 0.000003;
1429  (*this)[ +70.75] = 0.000003;
1430  (*this)[ +71.00] = 0.000003;
1431  (*this)[ +71.25] = 0.000003;
1432  (*this)[ +71.50] = 0.000003;
1433  (*this)[ +71.75] = 0.000003;
1434  (*this)[ +72.00] = 0.000002;
1435  (*this)[ +72.25] = 0.000002;
1436  (*this)[ +72.50] = 0.000002;
1437  (*this)[ +72.75] = 0.000002;
1438  (*this)[ +73.00] = 0.000002;
1439  (*this)[ +73.25] = 0.000002;
1440  (*this)[ +73.50] = 0.000002;
1441  (*this)[ +73.75] = 0.000002;
1442  (*this)[ +74.00] = 0.000002;
1443  (*this)[ +74.25] = 0.000002;
1444  (*this)[ +74.50] = 0.000002;
1445  (*this)[ +74.75] = 0.000002;
1446  (*this)[ +75.00] = 0.000002;
1447  (*this)[ +75.25] = 0.000002;
1448  (*this)[ +75.50] = 0.000002;
1449  (*this)[ +75.75] = 0.000002;
1450  (*this)[ +76.00] = 0.000002;
1451  (*this)[ +76.25] = 0.000002;
1452  (*this)[ +76.50] = 0.000002;
1453  (*this)[ +76.75] = 0.000001;
1454  (*this)[ +77.00] = 0.000001;
1455  (*this)[ +77.25] = 0.000001;
1456  (*this)[ +77.50] = 0.000001;
1457  (*this)[ +77.75] = 0.000001;
1458  (*this)[ +78.00] = 0.000001;
1459  (*this)[ +78.25] = 0.000001;
1460  (*this)[ +78.50] = 0.000001;
1461  (*this)[ +78.75] = 0.000001;
1462  (*this)[ +79.00] = 0.000001;
1463  (*this)[ +79.25] = 0.000001;
1464  (*this)[ +79.50] = 0.000001;
1465  (*this)[ +79.75] = 0.000001;
1466  (*this)[ +80.00] = 0.000001;
1467  (*this)[ +80.25] = 0.000001;
1468  (*this)[ +80.50] = 0.000001;
1469  (*this)[ +80.75] = 0.000001;
1470  (*this)[ +81.00] = 0.000001;
1471  (*this)[ +81.25] = 0.000001;
1472  (*this)[ +81.50] = 0.000001;
1473  (*this)[ +81.75] = 0.000001;
1474  (*this)[ +82.00] = 0.000001;
1475  (*this)[ +82.25] = 0.000001;
1476  (*this)[ +82.50] = 0.000001;
1477  (*this)[ +82.75] = 0.000001;
1478  (*this)[ +83.00] = 0.000001;
1479  (*this)[ +83.25] = 0.000001;
1480  (*this)[ +83.50] = 0.000001;
1481  (*this)[ +83.75] = 0.000001;
1482  (*this)[ +84.00] = 0.000001;
1483  (*this)[ +84.25] = 0.000001;
1484  (*this)[ +84.50] = 0.000001;
1485  (*this)[ +84.75] = 0.000001;
1486  (*this)[ +85.00] = 0.000001;
1487  (*this)[ +85.25] = 0.000001;
1488  (*this)[ +85.50] = 0.000001;
1489  (*this)[ +85.75] = 0.000001;
1490  (*this)[ +86.00] = 0.000001;
1491  (*this)[ +86.25] = 0.000001;
1492  (*this)[ +86.50] = 0.000001;
1493  (*this)[ +86.75] = 0.000001;
1494  (*this)[ +87.00] = 0.000001;
1495  (*this)[ +87.25] = 0.000001;
1496  (*this)[ +87.50] = 0.000001;
1497  (*this)[ +87.75] = 0.000001;
1498  (*this)[ +88.00] = 0.000000;
1499  (*this)[ +88.25] = 0.000000;
1500  (*this)[ +88.50] = 0.000000;
1501  (*this)[ +88.75] = 0.000000;
1502  (*this)[ +89.00] = 0.000000;
1503  (*this)[ +89.25] = 0.000000;
1504  (*this)[ +89.50] = 0.000000;
1505  (*this)[ +89.75] = 0.000001;
1506  (*this)[ +90.00] = 0.000000;
1507  (*this)[ +90.25] = 0.000000;
1508  (*this)[ +90.50] = 0.000000;
1509  (*this)[ +90.75] = 0.000000;
1510  (*this)[ +91.00] = 0.000000;
1511  (*this)[ +91.25] = 0.000000;
1512  (*this)[ +91.50] = 0.000000;
1513  (*this)[ +91.75] = 0.000000;
1514  (*this)[ +92.00] = 0.000000;
1515  (*this)[ +92.25] = 0.000000;
1516  (*this)[ +92.50] = 0.000000;
1517  (*this)[ +92.75] = 0.000000;
1518  (*this)[ +93.00] = 0.000000;
1519  (*this)[ +93.25] = 0.000000;
1520  (*this)[ +93.50] = 0.000000;
1521  (*this)[ +93.75] = 0.000000;
1522  (*this)[ +94.00] = 0.000000;
1523  (*this)[ +94.25] = 0.000000;
1524  (*this)[ +94.50] = 0.000000;
1525  (*this)[ +94.75] = 0.000000;
1526  (*this)[ +95.00] = 0.000000;
1527  (*this)[ +95.25] = 0.000000;
1528  (*this)[ +95.50] = 0.000000;
1529  (*this)[ +95.75] = 0.000000;
1530  (*this)[ +96.00] = 0.000000;
1531  (*this)[ +96.25] = 0.000000;
1532  (*this)[ +96.50] = 0.000000;
1533  (*this)[ +96.75] = 0.000000;
1534  (*this)[ +97.00] = 0.000000;
1535  (*this)[ +97.25] = 0.000000;
1536  (*this)[ +97.50] = 0.000000;
1537  (*this)[ +97.75] = 0.000000;
1538  (*this)[ +98.00] = 0.000000;
1539  (*this)[ +98.25] = 0.000000;
1540  (*this)[ +98.50] = 0.000000;
1541  (*this)[ +98.75] = 0.000000;
1542  (*this)[ +99.00] = 0.000000;
1543  (*this)[ +99.25] = 0.000000;
1544  (*this)[ +99.50] = 0.000000;
1545  (*this)[ +99.75] = 0.000000;
1546  (*this)[+100.00] = 0.000000;
1547 
1548  compile();
1549  };
1550  };
1551 
1552 
1553  /**
1554  * Implementation of probability distibution of the PMT transition times based on measurements.
1555  */
1556  struct JTransitionTimeProbability41 :
1557  public JTransitionTimeProbability_t
1558  {
1559  /**
1560  * Default constructor.
1561  */
1562  JTransitionTimeProbability41()
1563  {
1564  // produced by JLegolas.cc
1565  (*this)[ -20.00] = 0.000000;
1566  (*this)[ -19.75] = 0.000000;
1567  (*this)[ -19.50] = 0.000000;
1568  (*this)[ -19.25] = 0.000000;
1569  (*this)[ -19.00] = 0.000000;
1570  (*this)[ -18.75] = 0.000001;
1571  (*this)[ -18.50] = 0.000001;
1572  (*this)[ -18.25] = 0.000002;
1573  (*this)[ -18.00] = 0.000002;
1574  (*this)[ -17.75] = 0.000003;
1575  (*this)[ -17.50] = 0.000003;
1576  (*this)[ -17.25] = 0.000004;
1577  (*this)[ -17.00] = 0.000005;
1578  (*this)[ -16.75] = 0.000006;
1579  (*this)[ -16.50] = 0.000007;
1580  (*this)[ -16.25] = 0.000008;
1581  (*this)[ -16.00] = 0.000008;
1582  (*this)[ -15.75] = 0.000008;
1583  (*this)[ -15.50] = 0.000009;
1584  (*this)[ -15.25] = 0.000010;
1585  (*this)[ -15.00] = 0.000010;
1586  (*this)[ -14.75] = 0.000010;
1587  (*this)[ -14.50] = 0.000010;
1588  (*this)[ -14.25] = 0.000011;
1589  (*this)[ -14.00] = 0.000011;
1590  (*this)[ -13.75] = 0.000011;
1591  (*this)[ -13.50] = 0.000011;
1592  (*this)[ -13.25] = 0.000011;
1593  (*this)[ -13.00] = 0.000011;
1594  (*this)[ -12.75] = 0.000011;
1595  (*this)[ -12.50] = 0.000011;
1596  (*this)[ -12.25] = 0.000010;
1597  (*this)[ -12.00] = 0.000010;
1598  (*this)[ -11.75] = 0.000009;
1599  (*this)[ -11.50] = 0.000008;
1600  (*this)[ -11.25] = 0.000008;
1601  (*this)[ -11.00] = 0.000008;
1602  (*this)[ -10.75] = 0.000007;
1603  (*this)[ -10.50] = 0.000007;
1604  (*this)[ -10.25] = 0.000007;
1605  (*this)[ -10.00] = 0.000008;
1606  (*this)[ -9.75] = 0.000008;
1607  (*this)[ -9.50] = 0.000009;
1608  (*this)[ -9.25] = 0.000010;
1609  (*this)[ -9.00] = 0.000011;
1610  (*this)[ -8.75] = 0.000012;
1611  (*this)[ -8.50] = 0.000014;
1612  (*this)[ -8.25] = 0.000016;
1613  (*this)[ -8.00] = 0.000019;
1614  (*this)[ -7.75] = 0.000022;
1615  (*this)[ -7.50] = 0.000027;
1616  (*this)[ -7.25] = 0.000033;
1617  (*this)[ -7.00] = 0.000041;
1618  (*this)[ -6.75] = 0.000050;
1619  (*this)[ -6.50] = 0.000066;
1620  (*this)[ -6.25] = 0.000085;
1621  (*this)[ -6.00] = 0.000121;
1622  (*this)[ -5.75] = 0.000168;
1623  (*this)[ -5.50] = 0.000286;
1624  (*this)[ -5.25] = 0.000431;
1625  (*this)[ -5.00] = 0.000760;
1626  (*this)[ -4.75] = 0.001161;
1627  (*this)[ -4.50] = 0.002050;
1628  (*this)[ -4.25] = 0.003081;
1629  (*this)[ -4.00] = 0.005082;
1630  (*this)[ -3.75] = 0.007327;
1631  (*this)[ -3.50] = 0.011240;
1632  (*this)[ -3.25] = 0.015453;
1633  (*this)[ -3.00] = 0.021710;
1634  (*this)[ -2.75] = 0.028254;
1635  (*this)[ -2.50] = 0.036764;
1636  (*this)[ -2.25] = 0.045460;
1637  (*this)[ -2.00] = 0.055429;
1638  (*this)[ -1.75] = 0.065516;
1639  (*this)[ -1.50] = 0.076414;
1640  (*this)[ -1.25] = 0.087142;
1641  (*this)[ -1.00] = 0.096712;
1642  (*this)[ -0.75] = 0.105711;
1643  (*this)[ -0.50] = 0.110816;
1644  (*this)[ -0.25] = 0.115116;
1645  (*this)[ +0.00] = 0.113934;
1646  (*this)[ +0.25] = 0.112007;
1647  (*this)[ +0.50] = 0.105015;
1648  (*this)[ +0.75] = 0.097604;
1649  (*this)[ +1.00] = 0.087346;
1650  (*this)[ +1.25] = 0.077089;
1651  (*this)[ +1.50] = 0.066838;
1652  (*this)[ +1.75] = 0.056813;
1653  (*this)[ +2.00] = 0.048319;
1654  (*this)[ +2.25] = 0.040125;
1655  (*this)[ +2.50] = 0.033974;
1656  (*this)[ +2.75] = 0.028059;
1657  (*this)[ +3.00] = 0.023746;
1658  (*this)[ +3.25] = 0.019606;
1659  (*this)[ +3.50] = 0.016643;
1660  (*this)[ +3.75] = 0.013777;
1661  (*this)[ +4.00] = 0.011570;
1662  (*this)[ +4.25] = 0.009444;
1663  (*this)[ +4.50] = 0.007876;
1664  (*this)[ +4.75] = 0.006383;
1665  (*this)[ +5.00] = 0.005394;
1666  (*this)[ +5.25] = 0.004466;
1667  (*this)[ +5.50] = 0.003949;
1668  (*this)[ +5.75] = 0.003463;
1669  (*this)[ +6.00] = 0.003183;
1670  (*this)[ +6.25] = 0.002914;
1671  (*this)[ +6.50] = 0.002716;
1672  (*this)[ +6.75] = 0.002520;
1673  (*this)[ +7.00] = 0.002335;
1674  (*this)[ +7.25] = 0.002154;
1675  (*this)[ +7.50] = 0.001999;
1676  (*this)[ +7.75] = 0.001848;
1677  (*this)[ +8.00] = 0.001719;
1678  (*this)[ +8.25] = 0.001594;
1679  (*this)[ +8.50] = 0.001495;
1680  (*this)[ +8.75] = 0.001398;
1681  (*this)[ +9.00] = 0.001316;
1682  (*this)[ +9.25] = 0.001238;
1683  (*this)[ +9.50] = 0.001179;
1684  (*this)[ +9.75] = 0.001122;
1685  (*this)[ +10.00] = 0.001073;
1686  (*this)[ +10.25] = 0.001025;
1687  (*this)[ +10.50] = 0.000991;
1688  (*this)[ +10.75] = 0.000957;
1689  (*this)[ +11.00] = 0.000927;
1690  (*this)[ +11.25] = 0.000898;
1691  (*this)[ +11.50] = 0.000879;
1692  (*this)[ +11.75] = 0.000860;
1693  (*this)[ +12.00] = 0.000844;
1694  (*this)[ +12.25] = 0.000829;
1695  (*this)[ +12.50] = 0.000822;
1696  (*this)[ +12.75] = 0.000814;
1697  (*this)[ +13.00] = 0.000808;
1698  (*this)[ +13.25] = 0.000802;
1699  (*this)[ +13.50] = 0.000803;
1700  (*this)[ +13.75] = 0.000803;
1701  (*this)[ +14.00] = 0.000802;
1702  (*this)[ +14.25] = 0.000801;
1703  (*this)[ +14.50] = 0.000805;
1704  (*this)[ +14.75] = 0.000808;
1705  (*this)[ +15.00] = 0.000809;
1706  (*this)[ +15.25] = 0.000810;
1707  (*this)[ +15.50] = 0.000814;
1708  (*this)[ +15.75] = 0.000818;
1709  (*this)[ +16.00] = 0.000817;
1710  (*this)[ +16.25] = 0.000817;
1711  (*this)[ +16.50] = 0.000818;
1712  (*this)[ +16.75] = 0.000818;
1713  (*this)[ +17.00] = 0.000813;
1714  (*this)[ +17.25] = 0.000809;
1715  (*this)[ +17.50] = 0.000806;
1716  (*this)[ +17.75] = 0.000802;
1717  (*this)[ +18.00] = 0.000795;
1718  (*this)[ +18.25] = 0.000788;
1719  (*this)[ +18.50] = 0.000784;
1720  (*this)[ +18.75] = 0.000779;
1721  (*this)[ +19.00] = 0.000771;
1722  (*this)[ +19.25] = 0.000763;
1723  (*this)[ +19.50] = 0.000757;
1724  (*this)[ +19.75] = 0.000751;
1725  (*this)[ +20.00] = 0.000742;
1726  (*this)[ +20.25] = 0.000733;
1727  (*this)[ +20.50] = 0.000726;
1728  (*this)[ +20.75] = 0.000719;
1729  (*this)[ +21.00] = 0.000709;
1730  (*this)[ +21.25] = 0.000700;
1731  (*this)[ +21.50] = 0.000694;
1732  (*this)[ +21.75] = 0.000688;
1733  (*this)[ +22.00] = 0.000682;
1734  (*this)[ +22.25] = 0.000678;
1735  (*this)[ +22.50] = 0.000679;
1736  (*this)[ +22.75] = 0.000679;
1737  (*this)[ +23.00] = 0.000680;
1738  (*this)[ +23.25] = 0.000681;
1739  (*this)[ +23.50] = 0.000689;
1740  (*this)[ +23.75] = 0.000698;
1741  (*this)[ +24.00] = 0.000706;
1742  (*this)[ +24.25] = 0.000716;
1743  (*this)[ +24.50] = 0.000732;
1744  (*this)[ +24.75] = 0.000747;
1745  (*this)[ +25.00] = 0.000762;
1746  (*this)[ +25.25] = 0.000778;
1747  (*this)[ +25.50] = 0.000800;
1748  (*this)[ +25.75] = 0.000821;
1749  (*this)[ +26.00] = 0.000838;
1750  (*this)[ +26.25] = 0.000856;
1751  (*this)[ +26.50] = 0.000877;
1752  (*this)[ +26.75] = 0.000897;
1753  (*this)[ +27.00] = 0.000910;
1754  (*this)[ +27.25] = 0.000923;
1755  (*this)[ +27.50] = 0.000935;
1756  (*this)[ +27.75] = 0.000945;
1757  (*this)[ +28.00] = 0.000946;
1758  (*this)[ +28.25] = 0.000947;
1759  (*this)[ +28.50] = 0.000944;
1760  (*this)[ +28.75] = 0.000939;
1761  (*this)[ +29.00] = 0.000924;
1762  (*this)[ +29.25] = 0.000910;
1763  (*this)[ +29.50] = 0.000892;
1764  (*this)[ +29.75] = 0.000873;
1765  (*this)[ +30.00] = 0.000846;
1766  (*this)[ +30.25] = 0.000819;
1767  (*this)[ +30.50] = 0.000791;
1768  (*this)[ +30.75] = 0.000762;
1769  (*this)[ +31.00] = 0.000731;
1770  (*this)[ +31.25] = 0.000700;
1771  (*this)[ +31.50] = 0.000671;
1772  (*this)[ +31.75] = 0.000641;
1773  (*this)[ +32.00] = 0.000610;
1774  (*this)[ +32.25] = 0.000578;
1775  (*this)[ +32.50] = 0.000550;
1776  (*this)[ +32.75] = 0.000522;
1777  (*this)[ +33.00] = 0.000493;
1778  (*this)[ +33.25] = 0.000465;
1779  (*this)[ +33.50] = 0.000439;
1780  (*this)[ +33.75] = 0.000413;
1781  (*this)[ +34.00] = 0.000386;
1782  (*this)[ +34.25] = 0.000360;
1783  (*this)[ +34.50] = 0.000338;
1784  (*this)[ +34.75] = 0.000316;
1785  (*this)[ +35.00] = 0.000294;
1786  (*this)[ +35.25] = 0.000273;
1787  (*this)[ +35.50] = 0.000255;
1788  (*this)[ +35.75] = 0.000237;
1789  (*this)[ +36.00] = 0.000220;
1790  (*this)[ +36.25] = 0.000203;
1791  (*this)[ +36.50] = 0.000190;
1792  (*this)[ +36.75] = 0.000177;
1793  (*this)[ +37.00] = 0.000164;
1794  (*this)[ +37.25] = 0.000152;
1795  (*this)[ +37.50] = 0.000142;
1796  (*this)[ +37.75] = 0.000132;
1797  (*this)[ +38.00] = 0.000124;
1798  (*this)[ +38.25] = 0.000115;
1799  (*this)[ +38.50] = 0.000108;
1800  (*this)[ +38.75] = 0.000101;
1801  (*this)[ +39.00] = 0.000095;
1802  (*this)[ +39.25] = 0.000089;
1803  (*this)[ +39.50] = 0.000085;
1804  (*this)[ +39.75] = 0.000081;
1805  (*this)[ +40.00] = 0.000077;
1806  (*this)[ +40.25] = 0.000073;
1807  (*this)[ +40.50] = 0.000070;
1808  (*this)[ +40.75] = 0.000067;
1809  (*this)[ +41.00] = 0.000065;
1810  (*this)[ +41.25] = 0.000062;
1811  (*this)[ +41.50] = 0.000060;
1812  (*this)[ +41.75] = 0.000058;
1813  (*this)[ +42.00] = 0.000055;
1814  (*this)[ +42.25] = 0.000053;
1815  (*this)[ +42.50] = 0.000052;
1816  (*this)[ +42.75] = 0.000050;
1817  (*this)[ +43.00] = 0.000049;
1818  (*this)[ +43.25] = 0.000048;
1819  (*this)[ +43.50] = 0.000047;
1820  (*this)[ +43.75] = 0.000046;
1821  (*this)[ +44.00] = 0.000045;
1822  (*this)[ +44.25] = 0.000044;
1823  (*this)[ +44.50] = 0.000044;
1824  (*this)[ +44.75] = 0.000043;
1825  (*this)[ +45.00] = 0.000042;
1826  (*this)[ +45.25] = 0.000041;
1827  (*this)[ +45.50] = 0.000041;
1828  (*this)[ +45.75] = 0.000040;
1829  (*this)[ +46.00] = 0.000039;
1830  (*this)[ +46.25] = 0.000038;
1831  (*this)[ +46.50] = 0.000037;
1832  (*this)[ +46.75] = 0.000037;
1833  (*this)[ +47.00] = 0.000036;
1834  (*this)[ +47.25] = 0.000035;
1835  (*this)[ +47.50] = 0.000035;
1836  (*this)[ +47.75] = 0.000034;
1837  (*this)[ +48.00] = 0.000033;
1838  (*this)[ +48.25] = 0.000033;
1839  (*this)[ +48.50] = 0.000032;
1840  (*this)[ +48.75] = 0.000032;
1841  (*this)[ +49.00] = 0.000031;
1842  (*this)[ +49.25] = 0.000031;
1843  (*this)[ +49.50] = 0.000030;
1844  (*this)[ +49.75] = 0.000029;
1845  (*this)[ +50.00] = 0.000029;
1846  (*this)[ +50.25] = 0.000029;
1847  (*this)[ +50.50] = 0.000028;
1848  (*this)[ +50.75] = 0.000028;
1849  (*this)[ +51.00] = 0.000027;
1850  (*this)[ +51.25] = 0.000026;
1851  (*this)[ +51.50] = 0.000026;
1852  (*this)[ +51.75] = 0.000026;
1853  (*this)[ +52.00] = 0.000026;
1854  (*this)[ +52.25] = 0.000025;
1855  (*this)[ +52.50] = 0.000025;
1856  (*this)[ +52.75] = 0.000025;
1857  (*this)[ +53.00] = 0.000024;
1858  (*this)[ +53.25] = 0.000024;
1859  (*this)[ +53.50] = 0.000024;
1860  (*this)[ +53.75] = 0.000023;
1861  (*this)[ +54.00] = 0.000023;
1862  (*this)[ +54.25] = 0.000023;
1863  (*this)[ +54.50] = 0.000023;
1864  (*this)[ +54.75] = 0.000022;
1865  (*this)[ +55.00] = 0.000022;
1866  (*this)[ +55.25] = 0.000022;
1867  (*this)[ +55.50] = 0.000022;
1868  (*this)[ +55.75] = 0.000022;
1869  (*this)[ +56.00] = 0.000021;
1870  (*this)[ +56.25] = 0.000021;
1871  (*this)[ +56.50] = 0.000021;
1872  (*this)[ +56.75] = 0.000020;
1873  (*this)[ +57.00] = 0.000020;
1874  (*this)[ +57.25] = 0.000020;
1875  (*this)[ +57.50] = 0.000019;
1876  (*this)[ +57.75] = 0.000019;
1877  (*this)[ +58.00] = 0.000019;
1878  (*this)[ +58.25] = 0.000018;
1879  (*this)[ +58.50] = 0.000018;
1880  (*this)[ +58.75] = 0.000018;
1881  (*this)[ +59.00] = 0.000017;
1882  (*this)[ +59.25] = 0.000017;
1883  (*this)[ +59.50] = 0.000017;
1884  (*this)[ +59.75] = 0.000016;
1885  (*this)[ +60.00] = 0.000016;
1886  (*this)[ +60.25] = 0.000015;
1887  (*this)[ +60.50] = 0.000015;
1888  (*this)[ +60.75] = 0.000015;
1889  (*this)[ +61.00] = 0.000014;
1890  (*this)[ +61.25] = 0.000014;
1891  (*this)[ +61.50] = 0.000013;
1892  (*this)[ +61.75] = 0.000013;
1893  (*this)[ +62.00] = 0.000012;
1894  (*this)[ +62.25] = 0.000012;
1895  (*this)[ +62.50] = 0.000011;
1896  (*this)[ +62.75] = 0.000011;
1897  (*this)[ +63.00] = 0.000010;
1898  (*this)[ +63.25] = 0.000010;
1899  (*this)[ +63.50] = 0.000010;
1900  (*this)[ +63.75] = 0.000009;
1901  (*this)[ +64.00] = 0.000009;
1902  (*this)[ +64.25] = 0.000009;
1903  (*this)[ +64.50] = 0.000009;
1904  (*this)[ +64.75] = 0.000008;
1905  (*this)[ +65.00] = 0.000008;
1906  (*this)[ +65.25] = 0.000008;
1907  (*this)[ +65.50] = 0.000007;
1908  (*this)[ +65.75] = 0.000007;
1909  (*this)[ +66.00] = 0.000007;
1910  (*this)[ +66.25] = 0.000006;
1911  (*this)[ +66.50] = 0.000006;
1912  (*this)[ +66.75] = 0.000006;
1913  (*this)[ +67.00] = 0.000005;
1914  (*this)[ +67.25] = 0.000005;
1915  (*this)[ +67.50] = 0.000005;
1916  (*this)[ +67.75] = 0.000004;
1917  (*this)[ +68.00] = 0.000004;
1918  (*this)[ +68.25] = 0.000004;
1919  (*this)[ +68.50] = 0.000004;
1920  (*this)[ +68.75] = 0.000004;
1921  (*this)[ +69.00] = 0.000004;
1922  (*this)[ +69.25] = 0.000004;
1923  (*this)[ +69.50] = 0.000004;
1924  (*this)[ +69.75] = 0.000003;
1925  (*this)[ +70.00] = 0.000003;
1926  (*this)[ +70.25] = 0.000003;
1927  (*this)[ +70.50] = 0.000003;
1928  (*this)[ +70.75] = 0.000003;
1929  (*this)[ +71.00] = 0.000003;
1930  (*this)[ +71.25] = 0.000003;
1931  (*this)[ +71.50] = 0.000003;
1932  (*this)[ +71.75] = 0.000003;
1933  (*this)[ +72.00] = 0.000003;
1934  (*this)[ +72.25] = 0.000002;
1935  (*this)[ +72.50] = 0.000002;
1936  (*this)[ +72.75] = 0.000002;
1937  (*this)[ +73.00] = 0.000002;
1938  (*this)[ +73.25] = 0.000002;
1939  (*this)[ +73.50] = 0.000002;
1940  (*this)[ +73.75] = 0.000002;
1941  (*this)[ +74.00] = 0.000002;
1942  (*this)[ +74.25] = 0.000002;
1943  (*this)[ +74.50] = 0.000002;
1944  (*this)[ +74.75] = 0.000002;
1945  (*this)[ +75.00] = 0.000002;
1946  (*this)[ +75.25] = 0.000002;
1947  (*this)[ +75.50] = 0.000002;
1948  (*this)[ +75.75] = 0.000002;
1949  (*this)[ +76.00] = 0.000002;
1950  (*this)[ +76.25] = 0.000001;
1951  (*this)[ +76.50] = 0.000001;
1952  (*this)[ +76.75] = 0.000001;
1953  (*this)[ +77.00] = 0.000001;
1954  (*this)[ +77.25] = 0.000001;
1955  (*this)[ +77.50] = 0.000001;
1956  (*this)[ +77.75] = 0.000001;
1957  (*this)[ +78.00] = 0.000001;
1958  (*this)[ +78.25] = 0.000001;
1959  (*this)[ +78.50] = 0.000001;
1960  (*this)[ +78.75] = 0.000001;
1961  (*this)[ +79.00] = 0.000001;
1962  (*this)[ +79.25] = 0.000001;
1963  (*this)[ +79.50] = 0.000001;
1964  (*this)[ +79.75] = 0.000001;
1965  (*this)[ +80.00] = 0.000001;
1966  (*this)[ +80.25] = 0.000001;
1967  (*this)[ +80.50] = 0.000001;
1968  (*this)[ +80.75] = 0.000001;
1969  (*this)[ +81.00] = 0.000001;
1970  (*this)[ +81.25] = 0.000001;
1971  (*this)[ +81.50] = 0.000001;
1972  (*this)[ +81.75] = 0.000001;
1973  (*this)[ +82.00] = 0.000001;
1974  (*this)[ +82.25] = 0.000001;
1975  (*this)[ +82.50] = 0.000001;
1976  (*this)[ +82.75] = 0.000001;
1977  (*this)[ +83.00] = 0.000001;
1978  (*this)[ +83.25] = 0.000001;
1979  (*this)[ +83.50] = 0.000001;
1980  (*this)[ +83.75] = 0.000001;
1981  (*this)[ +84.00] = 0.000001;
1982  (*this)[ +84.25] = 0.000001;
1983  (*this)[ +84.50] = 0.000001;
1984  (*this)[ +84.75] = 0.000001;
1985  (*this)[ +85.00] = 0.000001;
1986  (*this)[ +85.25] = 0.000001;
1987  (*this)[ +85.50] = 0.000001;
1988  (*this)[ +85.75] = 0.000001;
1989  (*this)[ +86.00] = 0.000001;
1990  (*this)[ +86.25] = 0.000001;
1991  (*this)[ +86.50] = 0.000001;
1992  (*this)[ +86.75] = 0.000001;
1993  (*this)[ +87.00] = 0.000001;
1994  (*this)[ +87.25] = 0.000001;
1995  (*this)[ +87.50] = 0.000001;
1996  (*this)[ +87.75] = 0.000001;
1997  (*this)[ +88.00] = 0.000001;
1998  (*this)[ +88.25] = 0.000001;
1999  (*this)[ +88.50] = 0.000001;
2000  (*this)[ +88.75] = 0.000001;
2001  (*this)[ +89.00] = 0.000000;
2002  (*this)[ +89.25] = 0.000000;
2003  (*this)[ +89.50] = 0.000000;
2004  (*this)[ +89.75] = 0.000000;
2005  (*this)[ +90.00] = 0.000000;
2006  (*this)[ +90.25] = 0.000000;
2007  (*this)[ +90.50] = 0.000000;
2008  (*this)[ +90.75] = 0.000000;
2009  (*this)[ +91.00] = 0.000000;
2010  (*this)[ +91.25] = 0.000000;
2011  (*this)[ +91.50] = 0.000000;
2012  (*this)[ +91.75] = 0.000000;
2013  (*this)[ +92.00] = 0.000000;
2014  (*this)[ +92.25] = 0.000001;
2015  (*this)[ +92.50] = 0.000000;
2016  (*this)[ +92.75] = 0.000000;
2017  (*this)[ +93.00] = 0.000000;
2018  (*this)[ +93.25] = 0.000000;
2019  (*this)[ +93.50] = 0.000000;
2020  (*this)[ +93.75] = 0.000000;
2021  (*this)[ +94.00] = 0.000000;
2022  (*this)[ +94.25] = 0.000001;
2023  (*this)[ +94.50] = 0.000000;
2024  (*this)[ +94.75] = 0.000000;
2025  (*this)[ +95.00] = 0.000000;
2026  (*this)[ +95.25] = 0.000000;
2027  (*this)[ +95.50] = 0.000000;
2028  (*this)[ +95.75] = 0.000000;
2029  (*this)[ +96.00] = 0.000000;
2030  (*this)[ +96.25] = 0.000000;
2031  (*this)[ +96.50] = 0.000000;
2032  (*this)[ +96.75] = 0.000000;
2033  (*this)[ +97.00] = 0.000000;
2034  (*this)[ +97.25] = 0.000000;
2035  (*this)[ +97.50] = 0.000000;
2036  (*this)[ +97.75] = 0.000000;
2037  (*this)[ +98.00] = 0.000000;
2038  (*this)[ +98.25] = 0.000000;
2039  (*this)[ +98.50] = 0.000000;
2040  (*this)[ +98.75] = 0.000000;
2041  (*this)[ +99.00] = 0.000000;
2042  (*this)[ +99.25] = 0.000000;
2043  (*this)[ +99.50] = 0.000000;
2044  (*this)[ +99.75] = 0.000000;
2045  (*this)[+100.00] = 0.000000;
2046 
2047  compile();
2048  };
2049  };
2050 
2051 
2052  /**
2053  * Implementation of probability distibution of the PMT transition times based on measurements.
2054  */
2055  struct JTransitionTimeProbabilityXX :
2056  public JTransitionTimeProbability_t
2057  {
2058  /**
2059  * Default constructor.
2060  */
2061  JTransitionTimeProbabilityXX()
2062  {
2063 
2064  compile();
2065  };
2066  };
2067 
2068 
2069  /**
2070  */
2071  static const struct :
2072  public JTransitionTimeProbability
2073  {
2074  using JTransitionTimeProbability::operator();
2075 
2076  /**
2077  * Get probability for given transit time according given option.
2078  *
2079  * The option corresponds to the detector identifier of the measurement.
2080  *
2081  * \param x transit time [ns]
2082  * \param option option
2083  * \return probability
2084  */
2085  double operator()(const double x, const int option) const
2086  {
2087  switch (option) {
2088  case 38:
2089  return p38(x);
2090 
2091  case 41:
2092  return p41(x);
2093 
2094  default:
2095  return (*this)(x);
2096  };
2097  }
2098 
2099  const JTransitionTimeProbability38 p38;
2100  const JTransitionTimeProbability41 p41;
2101 
2102  } getTransitionTimeProbability; //!< Function object of transition time prabability.
2103 
2104 
2105 
2106  /**
2107  * Auxiliary class to generate twofold PMT transition times.
2108  */
2109  class JTransitionTimeProbability2D :
2110  public JTOOLS::JGridSplineFunction1D_t
2111  {
2112  public:
2113  /**
2114  * Default constructor.
2115  */
2116  JTransitionTimeProbability2D() :
2117  JTOOLS::JGridSplineFunction1D_t()
2118  {
2119  using namespace std;
2120 
2121  const double xmin = min( getTransitionTimeProbability.getXmin(), -getTransitionTimeProbability.getXmax());
2122  const double xmax = max(-getTransitionTimeProbability.getXmin(), getTransitionTimeProbability.getXmax());
2123  const double dx = (getTransitionTimeProbability.getXmax() - getTransitionTimeProbability.getXmin()) / (getTransitionTimeProbability.getSize() - 1);
2124 
2125  double Y = 0.0;
2126 
2127  for (double x1 = xmin; x1 < xmax + 0.5*dx; x1 += dx) {
2128 
2129  double y = 0.0;
2130 
2131  for (double x2 = xmin; x2 < xmax + 0.5*dx; x2 += dx) {
2132  y += getTransitionTimeProbability(x2) * getTransitionTimeProbability( x1 - x2);
2133  y += getTransitionTimeProbability(x2) * getTransitionTimeProbability(-x1 - x2);
2134  }
2135 
2136  (*this)[x1] = y;
2137 
2138  Y += y;
2139  }
2140 
2141  this->div(dx*Y);
2142 
2143  compile();
2144 
2145  setExceptionHandler(new JDefaultResult(0.0));
2146  }
2147  };
2148 
2149 
2150  /**
2151  * Function object of twofold transition time prabability.
2152  */
2153  static const JTransitionTimeProbability2D getTransitionTimeProbability2D;
2154 }
2155 
2156 #endif
JDETECTOR::JTransitionTimeProbability2D
Auxiliary class to generate twofold PMT transition times.
Definition: JPMTTransitTimeProbability.hh:2109
JDETECTOR::p38
const JTransitionTimeGenerator38 p38
Definition: JPMTTransitTimeGenerator.hh:3156
JDETECTOR::JTransitionTimeProbability38
Implementation of probability distibution of the PMT transition times based on measurements.
Definition: JPMTTransitTimeProbability.hh:1057
JTOOLS::JGridSplineFunction1D_t
Type definition of a spline interpolation based on a JGridCollection with result type double...
Definition: JFunction1D_t.hh:59
Y
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
Definition: detector-XY:fit1d.sh:33
JDETECTOR::getTransitionTimeProbability2D
static const JTransitionTimeProbability2D getTransitionTimeProbability2D
Function object of twofold transition time prabability.
Definition: JPMTTransitTimeProbability.hh:2153
JFunction1D_t.hh
JDETECTOR::JTransitionTimeProbability_t
Auxiliary base class for probability distibution of the PMT transition times based on measurements...
Definition: JPMTTransitTimeProbability.hh:22
JDETECTOR::getTransitionTimeProbability
getTransitionTimeProbability
Function object of transition time prabability.
JTOOLS::JFunctional< JSplineElement2D< double, double >::abscissa_type, double >::setExceptionHandler
void setExceptionHandler(const JSupervisor &supervisor)
Set the supervisor for handling of exceptions.
Definition: JFunctional.hh:288
JDETECTOR::JTransitionTimeProbability
Implementation of probability distibution of the PMT transition times based on measurements by Oleg K...
Definition: JPMTTransitTimeProbability.hh:39
JDETECTOR::p41
const JTransitionTimeGenerator41 p41
Definition: JPMTTransitTimeGenerator.hh:3157
JDETECTOR::JTransitionTimeProbability41
Implementation of probability distibution of the PMT transition times based on measurements.
Definition: JPMTTransitTimeProbability.hh:1556
JDETECTOR::JTransitionTimeProbability_t::JTransitionTimeProbability_t
JTransitionTimeProbability_t()
Default constructor.
Definition: JPMTTransitTimeProbability.hh:28
JDETECTOR::JTransitionTimeProbabilityXX
Implementation of probability distibution of the PMT transition times based on measurements.
Definition: JPMTTransitTimeProbability.hh:2055
Generated by   doxygen 1.8.5