Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
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
13namespace JDETECTOR {}
14namespace JPP { using namespace JDETECTOR; }
15
16namespace 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 * (see Calibration e-log entry <a href="https://elog.km3net.de/Analysis/467">467</a>).
2055 */
2056 struct JTransitionTimeProbability67 :
2057 public JTransitionTimeProbability_t
2058 {
2059 /**
2060 * Default constructor.
2061 */
2062 JTransitionTimeProbability67()
2063 {
2064 // produced by JLegolas.cc
2065 (*this)[ -20.00] = 0.000000;
2066 (*this)[ -19.75] = 0.000000;
2067 (*this)[ -19.50] = 0.000000;
2068 (*this)[ -19.25] = 0.000000;
2069 (*this)[ -19.00] = 0.000001;
2070 (*this)[ -18.75] = 0.000002;
2071 (*this)[ -18.50] = 0.000004;
2072 (*this)[ -18.25] = 0.000007;
2073 (*this)[ -18.00] = 0.000009;
2074 (*this)[ -17.75] = 0.000011;
2075 (*this)[ -17.50] = 0.000014;
2076 (*this)[ -17.25] = 0.000017;
2077 (*this)[ -17.00] = 0.000019;
2078 (*this)[ -16.75] = 0.000022;
2079 (*this)[ -16.50] = 0.000025;
2080 (*this)[ -16.25] = 0.000028;
2081 (*this)[ -16.00] = 0.000029;
2082 (*this)[ -15.75] = 0.000029;
2083 (*this)[ -15.50] = 0.000028;
2084 (*this)[ -15.25] = 0.000027;
2085 (*this)[ -15.00] = 0.000026;
2086 (*this)[ -14.75] = 0.000024;
2087 (*this)[ -14.50] = 0.000023;
2088 (*this)[ -14.25] = 0.000022;
2089 (*this)[ -14.00] = 0.000018;
2090 (*this)[ -13.75] = 0.000016;
2091 (*this)[ -13.50] = 0.000014;
2092 (*this)[ -13.25] = 0.000012;
2093 (*this)[ -13.00] = 0.000012;
2094 (*this)[ -12.75] = 0.000011;
2095 (*this)[ -12.50] = 0.000010;
2096 (*this)[ -12.25] = 0.000009;
2097 (*this)[ -12.00] = 0.000008;
2098 (*this)[ -11.75] = 0.000007;
2099 (*this)[ -11.50] = 0.000007;
2100 (*this)[ -11.25] = 0.000006;
2101 (*this)[ -11.00] = 0.000006;
2102 (*this)[ -10.75] = 0.000005;
2103 (*this)[ -10.50] = 0.000005;
2104 (*this)[ -10.25] = 0.000005;
2105 (*this)[ -10.00] = 0.000004;
2106 (*this)[ -9.75] = 0.000004;
2107 (*this)[ -9.50] = 0.000004;
2108 (*this)[ -9.25] = 0.000004;
2109 (*this)[ -9.00] = 0.000004;
2110 (*this)[ -8.75] = 0.000004;
2111 (*this)[ -8.50] = 0.000003;
2112 (*this)[ -8.25] = 0.000003;
2113 (*this)[ -8.00] = 0.000003;
2114 (*this)[ -7.75] = 0.000004;
2115 (*this)[ -7.50] = 0.000006;
2116 (*this)[ -7.25] = 0.000007;
2117 (*this)[ -7.00] = 0.000008;
2118 (*this)[ -6.75] = 0.000010;
2119 (*this)[ -6.50] = 0.000015;
2120 (*this)[ -6.25] = 0.000020;
2121 (*this)[ -6.00] = 0.000034;
2122 (*this)[ -5.75] = 0.000048;
2123 (*this)[ -5.50] = 0.000073;
2124 (*this)[ -5.25] = 0.000099;
2125 (*this)[ -5.00] = 0.000136;
2126 (*this)[ -4.75] = 0.000176;
2127 (*this)[ -4.50] = 0.000237;
2128 (*this)[ -4.25] = 0.000311;
2129 (*this)[ -4.00] = 0.000453;
2130 (*this)[ -3.75] = 0.000653;
2131 (*this)[ -3.50] = 0.001143;
2132 (*this)[ -3.25] = 0.001841;
2133 (*this)[ -3.00] = 0.003591;
2134 (*this)[ -2.75] = 0.005936;
2135 (*this)[ -2.50] = 0.011277;
2136 (*this)[ -2.25] = 0.017862;
2137 (*this)[ -2.00] = 0.030713;
2138 (*this)[ -1.75] = 0.045033;
2139 (*this)[ -1.50] = 0.066754;
2140 (*this)[ -1.25] = 0.088953;
2141 (*this)[ -1.00] = 0.113561;
2142 (*this)[ -0.75] = 0.136391;
2143 (*this)[ -0.50] = 0.150268;
2144 (*this)[ -0.25] = 0.161153;
2145 (*this)[ +0.00] = 0.156972;
2146 (*this)[ +0.25] = 0.150445;
2147 (*this)[ +0.50] = 0.132097;
2148 (*this)[ +0.75] = 0.113745;
2149 (*this)[ +1.00] = 0.095372;
2150 (*this)[ +1.25] = 0.077917;
2151 (*this)[ +1.50] = 0.065082;
2152 (*this)[ +1.75] = 0.053069;
2153 (*this)[ +2.00] = 0.045194;
2154 (*this)[ +2.25] = 0.037467;
2155 (*this)[ +2.50] = 0.030483;
2156 (*this)[ +2.75] = 0.023742;
2157 (*this)[ +3.00] = 0.018229;
2158 (*this)[ +3.25] = 0.013161;
2159 (*this)[ +3.50] = 0.010332;
2160 (*this)[ +3.75] = 0.007812;
2161 (*this)[ +4.00] = 0.006844;
2162 (*this)[ +4.25] = 0.005969;
2163 (*this)[ +4.50] = 0.005559;
2164 (*this)[ +4.75] = 0.005147;
2165 (*this)[ +5.00] = 0.004726;
2166 (*this)[ +5.25] = 0.004306;
2167 (*this)[ +5.50] = 0.003888;
2168 (*this)[ +5.75] = 0.003486;
2169 (*this)[ +6.00] = 0.003165;
2170 (*this)[ +6.25] = 0.002854;
2171 (*this)[ +6.50] = 0.002594;
2172 (*this)[ +6.75] = 0.002344;
2173 (*this)[ +7.00] = 0.002145;
2174 (*this)[ +7.25] = 0.001955;
2175 (*this)[ +7.50] = 0.001803;
2176 (*this)[ +7.75] = 0.001658;
2177 (*this)[ +8.00] = 0.001541;
2178 (*this)[ +8.25] = 0.001428;
2179 (*this)[ +8.50] = 0.001333;
2180 (*this)[ +8.75] = 0.001241;
2181 (*this)[ +9.00] = 0.001170;
2182 (*this)[ +9.25] = 0.001099;
2183 (*this)[ +9.50] = 0.001030;
2184 (*this)[ +9.75] = 0.000965;
2185 (*this)[ +10.00] = 0.000919;
2186 (*this)[ +10.25] = 0.000873;
2187 (*this)[ +10.50] = 0.000825;
2188 (*this)[ +10.75] = 0.000779;
2189 (*this)[ +11.00] = 0.000744;
2190 (*this)[ +11.25] = 0.000708;
2191 (*this)[ +11.50] = 0.000674;
2192 (*this)[ +11.75] = 0.000640;
2193 (*this)[ +12.00] = 0.000615;
2194 (*this)[ +12.25] = 0.000589;
2195 (*this)[ +12.50] = 0.000562;
2196 (*this)[ +12.75] = 0.000537;
2197 (*this)[ +13.00] = 0.000519;
2198 (*this)[ +13.25] = 0.000501;
2199 (*this)[ +13.50] = 0.000484;
2200 (*this)[ +13.75] = 0.000468;
2201 (*this)[ +14.00] = 0.000458;
2202 (*this)[ +14.25] = 0.000449;
2203 (*this)[ +14.50] = 0.000441;
2204 (*this)[ +14.75] = 0.000435;
2205 (*this)[ +15.00] = 0.000438;
2206 (*this)[ +15.25] = 0.000442;
2207 (*this)[ +15.50] = 0.000447;
2208 (*this)[ +15.75] = 0.000452;
2209 (*this)[ +16.00] = 0.000461;
2210 (*this)[ +16.25] = 0.000469;
2211 (*this)[ +16.50] = 0.000475;
2212 (*this)[ +16.75] = 0.000482;
2213 (*this)[ +17.00] = 0.000499;
2214 (*this)[ +17.25] = 0.000513;
2215 (*this)[ +17.50] = 0.000516;
2216 (*this)[ +17.75] = 0.000522;
2217 (*this)[ +18.00] = 0.000533;
2218 (*this)[ +18.25] = 0.000544;
2219 (*this)[ +18.50] = 0.000551;
2220 (*this)[ +18.75] = 0.000559;
2221 (*this)[ +19.00] = 0.000574;
2222 (*this)[ +19.25] = 0.000587;
2223 (*this)[ +19.50] = 0.000591;
2224 (*this)[ +19.75] = 0.000595;
2225 (*this)[ +20.00] = 0.000607;
2226 (*this)[ +20.25] = 0.000618;
2227 (*this)[ +20.50] = 0.000624;
2228 (*this)[ +20.75] = 0.000630;
2229 (*this)[ +21.00] = 0.000637;
2230 (*this)[ +21.25] = 0.000643;
2231 (*this)[ +21.50] = 0.000646;
2232 (*this)[ +21.75] = 0.000649;
2233 (*this)[ +22.00] = 0.000657;
2234 (*this)[ +22.25] = 0.000664;
2235 (*this)[ +22.50] = 0.000666;
2236 (*this)[ +22.75] = 0.000670;
2237 (*this)[ +23.00] = 0.000679;
2238 (*this)[ +23.25] = 0.000687;
2239 (*this)[ +23.50] = 0.000694;
2240 (*this)[ +23.75] = 0.000702;
2241 (*this)[ +24.00] = 0.000716;
2242 (*this)[ +24.25] = 0.000730;
2243 (*this)[ +24.50] = 0.000738;
2244 (*this)[ +24.75] = 0.000748;
2245 (*this)[ +25.00] = 0.000769;
2246 (*this)[ +25.25] = 0.000789;
2247 (*this)[ +25.50] = 0.000803;
2248 (*this)[ +25.75] = 0.000819;
2249 (*this)[ +26.00] = 0.000843;
2250 (*this)[ +26.25] = 0.000864;
2251 (*this)[ +26.50] = 0.000877;
2252 (*this)[ +26.75] = 0.000891;
2253 (*this)[ +27.00] = 0.000907;
2254 (*this)[ +27.25] = 0.000921;
2255 (*this)[ +27.50] = 0.000921;
2256 (*this)[ +27.75] = 0.000921;
2257 (*this)[ +28.00] = 0.000918;
2258 (*this)[ +28.25] = 0.000911;
2259 (*this)[ +28.50] = 0.000886;
2260 (*this)[ +28.75] = 0.000860;
2261 (*this)[ +29.00] = 0.000832;
2262 (*this)[ +29.25] = 0.000804;
2263 (*this)[ +29.50] = 0.000770;
2264 (*this)[ +29.75] = 0.000736;
2265 (*this)[ +30.00] = 0.000699;
2266 (*this)[ +30.25] = 0.000660;
2267 (*this)[ +30.50] = 0.000615;
2268 (*this)[ +30.75] = 0.000572;
2269 (*this)[ +31.00] = 0.000537;
2270 (*this)[ +31.25] = 0.000501;
2271 (*this)[ +31.50] = 0.000458;
2272 (*this)[ +31.75] = 0.000417;
2273 (*this)[ +32.00] = 0.000381;
2274 (*this)[ +32.25] = 0.000345;
2275 (*this)[ +32.50] = 0.000313;
2276 (*this)[ +32.75] = 0.000282;
2277 (*this)[ +33.00] = 0.000257;
2278 (*this)[ +33.25] = 0.000231;
2279 (*this)[ +33.50] = 0.000207;
2280 (*this)[ +33.75] = 0.000185;
2281 (*this)[ +34.00] = 0.000172;
2282 (*this)[ +34.25] = 0.000158;
2283 (*this)[ +34.50] = 0.000143;
2284 (*this)[ +34.75] = 0.000128;
2285 (*this)[ +35.00] = 0.000117;
2286 (*this)[ +35.25] = 0.000106;
2287 (*this)[ +35.50] = 0.000096;
2288 (*this)[ +35.75] = 0.000087;
2289 (*this)[ +36.00] = 0.000079;
2290 (*this)[ +36.25] = 0.000072;
2291 (*this)[ +36.50] = 0.000067;
2292 (*this)[ +36.75] = 0.000063;
2293 (*this)[ +37.00] = 0.000059;
2294 (*this)[ +37.25] = 0.000055;
2295 (*this)[ +37.50] = 0.000052;
2296 (*this)[ +37.75] = 0.000049;
2297 (*this)[ +38.00] = 0.000047;
2298 (*this)[ +38.25] = 0.000046;
2299 (*this)[ +38.50] = 0.000046;
2300 (*this)[ +38.75] = 0.000046;
2301 (*this)[ +39.00] = 0.000044;
2302 (*this)[ +39.25] = 0.000043;
2303 (*this)[ +39.50] = 0.000042;
2304 (*this)[ +39.75] = 0.000041;
2305 (*this)[ +40.00] = 0.000039;
2306 (*this)[ +40.25] = 0.000038;
2307 (*this)[ +40.50] = 0.000037;
2308 (*this)[ +40.75] = 0.000035;
2309 (*this)[ +41.00] = 0.000033;
2310 (*this)[ +41.25] = 0.000031;
2311 (*this)[ +41.50] = 0.000030;
2312 (*this)[ +41.75] = 0.000029;
2313 (*this)[ +42.00] = 0.000029;
2314 (*this)[ +42.25] = 0.000028;
2315 (*this)[ +42.50] = 0.000029;
2316 (*this)[ +42.75] = 0.000029;
2317 (*this)[ +43.00] = 0.000029;
2318 (*this)[ +43.25] = 0.000029;
2319 (*this)[ +43.50] = 0.000029;
2320 (*this)[ +43.75] = 0.000029;
2321 (*this)[ +44.00] = 0.000028;
2322 (*this)[ +44.25] = 0.000026;
2323 (*this)[ +44.50] = 0.000026;
2324 (*this)[ +44.75] = 0.000025;
2325 (*this)[ +45.00] = 0.000023;
2326 (*this)[ +45.25] = 0.000022;
2327 (*this)[ +45.50] = 0.000022;
2328 (*this)[ +45.75] = 0.000022;
2329 (*this)[ +46.00] = 0.000022;
2330 (*this)[ +46.25] = 0.000022;
2331 (*this)[ +46.50] = 0.000023;
2332 (*this)[ +46.75] = 0.000023;
2333 (*this)[ +47.00] = 0.000022;
2334 (*this)[ +47.25] = 0.000022;
2335 (*this)[ +47.50] = 0.000022;
2336 (*this)[ +47.75] = 0.000022;
2337 (*this)[ +48.00] = 0.000022;
2338 (*this)[ +48.25] = 0.000022;
2339 (*this)[ +48.50] = 0.000022;
2340 (*this)[ +48.75] = 0.000021;
2341 (*this)[ +49.00] = 0.000021;
2342 (*this)[ +49.25] = 0.000021;
2343 (*this)[ +49.50] = 0.000021;
2344 (*this)[ +49.75] = 0.000021;
2345 (*this)[ +50.00] = 0.000020;
2346 (*this)[ +50.25] = 0.000020;
2347 (*this)[ +50.50] = 0.000020;
2348 (*this)[ +50.75] = 0.000019;
2349 (*this)[ +51.00] = 0.000019;
2350 (*this)[ +51.25] = 0.000018;
2351 (*this)[ +51.50] = 0.000018;
2352 (*this)[ +51.75] = 0.000017;
2353 (*this)[ +52.00] = 0.000018;
2354 (*this)[ +52.25] = 0.000018;
2355 (*this)[ +52.50] = 0.000018;
2356 (*this)[ +52.75] = 0.000018;
2357 (*this)[ +53.00] = 0.000018;
2358 (*this)[ +53.25] = 0.000019;
2359 (*this)[ +53.50] = 0.000018;
2360 (*this)[ +53.75] = 0.000018;
2361 (*this)[ +54.00] = 0.000019;
2362 (*this)[ +54.25] = 0.000019;
2363 (*this)[ +54.50] = 0.000019;
2364 (*this)[ +54.75] = 0.000018;
2365 (*this)[ +55.00] = 0.000018;
2366 (*this)[ +55.25] = 0.000017;
2367 (*this)[ +55.50] = 0.000016;
2368 (*this)[ +55.75] = 0.000016;
2369 (*this)[ +56.00] = 0.000016;
2370 (*this)[ +56.25] = 0.000016;
2371 (*this)[ +56.50] = 0.000015;
2372 (*this)[ +56.75] = 0.000015;
2373 (*this)[ +57.00] = 0.000015;
2374 (*this)[ +57.25] = 0.000015;
2375 (*this)[ +57.50] = 0.000014;
2376 (*this)[ +57.75] = 0.000013;
2377 (*this)[ +58.00] = 0.000012;
2378 (*this)[ +58.25] = 0.000011;
2379 (*this)[ +58.50] = 0.000011;
2380 (*this)[ +58.75] = 0.000010;
2381 (*this)[ +59.00] = 0.000011;
2382 (*this)[ +59.25] = 0.000011;
2383 (*this)[ +59.50] = 0.000009;
2384 (*this)[ +59.75] = 0.000008;
2385 (*this)[ +60.00] = 0.000008;
2386 (*this)[ +60.25] = 0.000008;
2387 (*this)[ +60.50] = 0.000008;
2388 (*this)[ +60.75] = 0.000008;
2389 (*this)[ +61.00] = 0.000007;
2390 (*this)[ +61.25] = 0.000006;
2391 (*this)[ +61.50] = 0.000005;
2392 (*this)[ +61.75] = 0.000004;
2393 (*this)[ +62.00] = 0.000005;
2394 (*this)[ +62.25] = 0.000005;
2395 (*this)[ +62.50] = 0.000004;
2396 (*this)[ +62.75] = 0.000004;
2397 (*this)[ +63.00] = 0.000004;
2398 (*this)[ +63.25] = 0.000004;
2399 (*this)[ +63.50] = 0.000004;
2400 (*this)[ +63.75] = 0.000004;
2401 (*this)[ +64.00] = 0.000004;
2402 (*this)[ +64.25] = 0.000003;
2403 (*this)[ +64.50] = 0.000002;
2404 (*this)[ +64.75] = 0.000002;
2405 (*this)[ +65.00] = 0.000001;
2406 (*this)[ +65.25] = 0.000001;
2407 (*this)[ +65.50] = 0.000001;
2408 (*this)[ +65.75] = 0.000001;
2409 (*this)[ +66.00] = 0.000001;
2410 (*this)[ +66.25] = 0.000001;
2411 (*this)[ +66.50] = 0.000001;
2412 (*this)[ +66.75] = 0.000002;
2413 (*this)[ +67.00] = 0.000002;
2414 (*this)[ +67.25] = 0.000002;
2415 (*this)[ +67.50] = 0.000003;
2416 (*this)[ +67.75] = 0.000004;
2417 (*this)[ +68.00] = 0.000003;
2418 (*this)[ +68.25] = 0.000003;
2419 (*this)[ +68.50] = 0.000002;
2420 (*this)[ +68.75] = 0.000001;
2421 (*this)[ +69.00] = 0.000001;
2422 (*this)[ +69.25] = 0.000000;
2423 (*this)[ +69.50] = 0.000000;
2424 (*this)[ +69.75] = 0.000000;
2425 (*this)[ +70.00] = 0.000000;
2426 (*this)[ +70.25] = 0.000000;
2427 (*this)[ +70.50] = 0.000000;
2428 (*this)[ +70.75] = 0.000000;
2429 (*this)[ +71.00] = 0.000000;
2430 (*this)[ +71.25] = 0.000001;
2431 (*this)[ +71.50] = 0.000001;
2432 (*this)[ +71.75] = 0.000002;
2433 (*this)[ +72.00] = 0.000002;
2434 (*this)[ +72.25] = 0.000002;
2435 (*this)[ +72.50] = 0.000002;
2436 (*this)[ +72.75] = 0.000002;
2437 (*this)[ +73.00] = 0.000002;
2438 (*this)[ +73.25] = 0.000002;
2439 (*this)[ +73.50] = 0.000002;
2440 (*this)[ +73.75] = 0.000002;
2441 (*this)[ +74.00] = 0.000001;
2442 (*this)[ +74.25] = 0.000001;
2443 (*this)[ +74.50] = 0.000001;
2444 (*this)[ +74.75] = 0.000001;
2445 (*this)[ +75.00] = 0.000001;
2446 (*this)[ +75.25] = 0.000000;
2447 (*this)[ +75.50] = 0.000000;
2448 (*this)[ +75.75] = 0.000000;
2449 (*this)[ +76.00] = 0.000000;
2450 (*this)[ +76.25] = 0.000000;
2451 (*this)[ +76.50] = 0.000000;
2452 (*this)[ +76.75] = 0.000000;
2453 (*this)[ +77.00] = 0.000001;
2454 (*this)[ +77.25] = 0.000001;
2455 (*this)[ +77.50] = 0.000001;
2456 (*this)[ +77.75] = 0.000001;
2457 (*this)[ +78.00] = 0.000001;
2458 (*this)[ +78.25] = 0.000001;
2459 (*this)[ +78.50] = 0.000001;
2460 (*this)[ +78.75] = 0.000001;
2461 (*this)[ +79.00] = 0.000001;
2462 (*this)[ +79.25] = 0.000001;
2463 (*this)[ +79.50] = 0.000001;
2464 (*this)[ +79.75] = 0.000001;
2465 (*this)[ +80.00] = 0.000001;
2466 (*this)[ +80.25] = 0.000001;
2467 (*this)[ +80.50] = 0.000000;
2468 (*this)[ +80.75] = 0.000000;
2469 (*this)[ +81.00] = 0.000000;
2470 (*this)[ +81.25] = 0.000000;
2471 (*this)[ +81.50] = 0.000000;
2472 (*this)[ +81.75] = 0.000000;
2473 (*this)[ +82.00] = 0.000000;
2474 (*this)[ +82.25] = 0.000000;
2475 (*this)[ +82.50] = 0.000000;
2476 (*this)[ +82.75] = 0.000000;
2477 (*this)[ +83.00] = 0.000000;
2478 (*this)[ +83.25] = 0.000000;
2479 (*this)[ +83.50] = 0.000000;
2480 (*this)[ +83.75] = 0.000000;
2481 (*this)[ +84.00] = 0.000000;
2482 (*this)[ +84.25] = 0.000000;
2483 (*this)[ +84.50] = 0.000000;
2484 (*this)[ +84.75] = 0.000000;
2485 (*this)[ +85.00] = 0.000000;
2486 (*this)[ +85.25] = 0.000000;
2487 (*this)[ +85.50] = 0.000001;
2488 (*this)[ +85.75] = 0.000001;
2489 (*this)[ +86.00] = 0.000002;
2490 (*this)[ +86.25] = 0.000002;
2491 (*this)[ +86.50] = 0.000003;
2492 (*this)[ +86.75] = 0.000003;
2493 (*this)[ +87.00] = 0.000002;
2494 (*this)[ +87.25] = 0.000002;
2495 (*this)[ +87.50] = 0.000003;
2496 (*this)[ +87.75] = 0.000003;
2497 (*this)[ +88.00] = 0.000003;
2498 (*this)[ +88.25] = 0.000003;
2499 (*this)[ +88.50] = 0.000002;
2500 (*this)[ +88.75] = 0.000002;
2501 (*this)[ +89.00] = 0.000001;
2502 (*this)[ +89.25] = 0.000001;
2503 (*this)[ +89.50] = 0.000001;
2504 (*this)[ +89.75] = 0.000002;
2505 (*this)[ +90.00] = 0.000001;
2506 (*this)[ +90.25] = 0.000001;
2507 (*this)[ +90.50] = 0.000001;
2508 (*this)[ +90.75] = 0.000001;
2509 (*this)[ +91.00] = 0.000001;
2510 (*this)[ +91.25] = 0.000002;
2511 (*this)[ +91.50] = 0.000002;
2512 (*this)[ +91.75] = 0.000002;
2513 (*this)[ +92.00] = 0.000002;
2514 (*this)[ +92.25] = 0.000002;
2515 (*this)[ +92.50] = 0.000002;
2516 (*this)[ +92.75] = 0.000002;
2517 (*this)[ +93.00] = 0.000002;
2518 (*this)[ +93.25] = 0.000002;
2519 (*this)[ +93.50] = 0.000002;
2520 (*this)[ +93.75] = 0.000001;
2521 (*this)[ +94.00] = 0.000001;
2522 (*this)[ +94.25] = 0.000001;
2523 (*this)[ +94.50] = 0.000001;
2524 (*this)[ +94.75] = 0.000000;
2525 (*this)[ +95.00] = 0.000001;
2526 (*this)[ +95.25] = 0.000001;
2527 (*this)[ +95.50] = 0.000001;
2528 (*this)[ +95.75] = 0.000001;
2529 (*this)[ +96.00] = 0.000000;
2530 (*this)[ +96.25] = 0.000000;
2531 (*this)[ +96.50] = 0.000000;
2532 (*this)[ +96.75] = 0.000000;
2533 (*this)[ +97.00] = 0.000000;
2534 (*this)[ +97.25] = 0.000000;
2535 (*this)[ +97.50] = 0.000000;
2536 (*this)[ +97.75] = 0.000000;
2537 (*this)[ +98.00] = 0.000000;
2538 (*this)[ +98.25] = 0.000000;
2539 (*this)[ +98.50] = 0.000000;
2540 (*this)[ +98.75] = 0.000000;
2541 (*this)[ +99.00] = 0.000000;
2542 (*this)[ +99.25] = 0.000000;
2543 (*this)[ +99.50] = 0.000000;
2544 (*this)[ +99.75] = 0.000000;
2545 (*this)[+100.00] = 0.000000;
2546
2547 compile();
2548 };
2549 };
2550
2551
2552 /**
2553 * Implementation of probability distibution of the PMT transition times based on measurements.
2554 */
2555 struct JTransitionTimeProbabilityXX :
2556 public JTransitionTimeProbability_t
2557 {
2558 /**
2559 * Default constructor.
2560 */
2561 JTransitionTimeProbabilityXX()
2562 {
2563
2564 compile();
2565 };
2566 };
2567
2568
2569 /**
2570 */
2571 static const struct getTransitionTimeProbability :
2572 public JTransitionTimeProbability
2573 {
2574 using JTransitionTimeProbability::operator();
2575
2576 /**
2577 * Get probability for given transit time according given option.
2578 *
2579 * The option corresponds to the detector identifier of the measurement.
2580 *
2581 * \param x transit time [ns]
2582 * \param option option
2583 * \return probability
2584 */
2585 double operator()(const double x, const int option) const
2586 {
2587 switch (option) {
2588 case 38:
2589 return p38(x);
2590
2591 case 41:
2592 return p41(x);
2593
2594 case 67:
2595 return p67(x);
2596
2597 case 1:
2598 return (*this)(x);
2599
2600 default:
2601 return (*this)(x);
2602 };
2603 }
2604
2605 const JTransitionTimeProbability38 p38;
2606 const JTransitionTimeProbability41 p41;
2607 const JTransitionTimeProbability67 p67;
2608
2609 } getTransitionTimeProbability; //!< Function object of transition time prabability.
2610
2611
2612
2613 /**
2614 * Auxiliary class to generate twofold PMT transition times.
2615 */
2616 class JTransitionTimeProbability2D :
2617 public JTOOLS::JGridSplineFunction1D_t
2618 {
2619 public:
2620 /**
2621 * Default constructor.
2622 */
2623 JTransitionTimeProbability2D() :
2624 JTOOLS::JGridSplineFunction1D_t()
2625 {
2626 using namespace std;
2627
2628 const double xmin = min( getTransitionTimeProbability.getXmin(), -getTransitionTimeProbability.getXmax());
2629 const double xmax = max(-getTransitionTimeProbability.getXmin(), getTransitionTimeProbability.getXmax());
2630 const double dx = (getTransitionTimeProbability.getXmax() - getTransitionTimeProbability.getXmin()) / (getTransitionTimeProbability.getSize() - 1);
2631
2632 double Y = 0.0;
2633
2634 for (double x1 = xmin; x1 < xmax + 0.5*dx; x1 += dx) {
2635
2636 double y = 0.0;
2637
2638 for (double x2 = xmin; x2 < xmax + 0.5*dx; x2 += dx) {
2639 y += getTransitionTimeProbability(x2) * getTransitionTimeProbability( x1 - x2);
2640 y += getTransitionTimeProbability(x2) * getTransitionTimeProbability(-x1 - x2);
2641 }
2642
2643 (*this)[x1] = y;
2644
2645 Y += y;
2646 }
2647
2648 this->div(dx*Y);
2649
2650 compile();
2651
2652 setExceptionHandler(new JDefaultResult(0.0));
2653 }
2654 };
2655
2656
2657 /**
2658 * Function object of twofold transition time prabability.
2659 */
2660 static const JTransitionTimeProbability2D getTransitionTimeProbability2D;
2661}
2662
2663#endif
JFunction1D_t.hh
JDETECTOR::JTransitionTimeProbability2D
Auxiliary class to generate twofold PMT transition times.
Definition JPMTTransitTimeProbability.hh:2618
JDETECTOR
file Auxiliary data structures and methods for detector calibration.
Definition JAnchor.hh:12
JDETECTOR::getTransitionTimeProbability2D
static const JTransitionTimeProbability2D getTransitionTimeProbability2D
Function object of twofold transition time prabability.
Definition JPMTTransitTimeProbability.hh:2660
JDETECTOR::getTransitionTimeProbability
JDETECTOR::getTransitionTimeProbability getTransitionTimeProbability
Function object of transition time prabability.
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:43
JTOOLS
Auxiliary classes and methods for multi-dimensional interpolations and histograms.
Definition JAbstractCollection.hh:11
JDETECTOR::JTransitionTimeProbability38
Implementation of probability distibution of the PMT transition times based on measurements.
Definition JPMTTransitTimeProbability.hh:1059
JDETECTOR::JTransitionTimeProbability41
Implementation of probability distibution of the PMT transition times based on measurements.
Definition JPMTTransitTimeProbability.hh:1558
JDETECTOR::JTransitionTimeProbability67
Implementation of probability distibution of the PMT transition times based on measurements (see Cali...
Definition JPMTTransitTimeProbability.hh:2058
JDETECTOR::JTransitionTimeProbabilityXX
Implementation of probability distibution of the PMT transition times based on measurements.
Definition JPMTTransitTimeProbability.hh:2557
JDETECTOR::JTransitionTimeProbability_t
Auxiliary base class for probability distibution of the PMT transition times based on measurements.
Definition JPMTTransitTimeProbability.hh:24
JDETECTOR::JTransitionTimeProbability_t::JTransitionTimeProbability_t
JTransitionTimeProbability_t()
Default constructor.
Definition JPMTTransitTimeProbability.hh:28
JDETECTOR::JTransitionTimeProbability
Implementation of probability distibution of the PMT transition times based on measurements by Oleg K...
Definition JPMTTransitTimeProbability.hh:41
JDETECTOR::getTransitionTimeProbability::operator()
double operator()(const double x, const int option) const
Get probability for given transit time according given option.
Definition JPMTTransitTimeProbability.hh:2585
JTOOLS::JGridSplineFunction1D_t
Type definition of a spline interpolation based on a JGridCollection with result type double.
Definition JFunction1D_t.hh:61
Generated by doxygen 1.12.0