BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDtoKKpiPlot.cc
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2//
3// Environment:
4// This software is part of models developed at BES collaboration
5// based on the EvtGen framework. If you use all or part
6// of it, please give an appropriate acknowledgement.
7//
8// Copyright Information: See EvtGen/BesCopyright
9// Copyright (A) 2006 Ping Rong-Gang @IHEP
10//
11// Module: EvtDtoKKpiPlot.cc
12//
13// Modification history:
14//
15// Liaoyuan Dong August, 2022 Module created
16//
17//------------------------------------------------------------------------
18#include "EvtDtoKKpiPlot.hh"
19#include "../EvtGenBase/EvtGenKine.hh"
20#include "../EvtGenBase/EvtPDL.hh"
21#include "../EvtGenBase/EvtParticle.hh"
22#include "../EvtGenBase/EvtPatches.hh"
23#include "../EvtGenBase/EvtRandom.hh"
24#include "../EvtGenBase/EvtReport.hh"
25#include <stdlib.h>
26#include <string>
27
28#include "TApplication.h"
29#include "TAxis.h"
30#include "TFile.h"
31#include "TH1.h"
32#include "TH2.h"
33#include "TROOT.h"
34using std::endl;
35
36EvtDtoKKpiPlot::~EvtDtoKKpiPlot() {}
37
38void EvtDtoKKpiPlot::getName( std::string& model_name ) { model_name = "DtoKKpiPlot"; }
39
40EvtDecayBase* EvtDtoKKpiPlot::clone() { return new EvtDtoKKpiPlot; }
41
42void EvtDtoKKpiPlot::init() {
43
44 checkNArg( 0 );
45
46 bool idN = getDaugs()[0] == EvtPDL::getId( std::string( "K-" ) ) ||
47 getDaugs()[0] == EvtPDL::getId( std::string( "K+" ) );
48 bool idKs = getDaugs()[1] == EvtPDL::getId( std::string( "K+" ) ) ||
49 getDaugs()[1] == EvtPDL::getId( std::string( "K-" ) );
50 bool idPi = getDaugs()[2] == EvtPDL::getId( std::string( "pi+" ) ) ||
51 getDaugs()[2] == EvtPDL::getId( std::string( "pi-" ) );
52 if ( !( idN && idKs && idPi ) )
53 {
54 std::cout << "EvtDtoKKpiPlot: the daughter sequence should be K- K+ pi+" << std::endl;
55 abort();
56 }
57 EvtSpinType::spintype parenttype = EvtPDL::getSpinType( getParentId() );
58
59 Xmin = 0.35;
60 Xmax = 1.91;
61 Xwid = 0.026;
62 Ymin = 0.35;
63 Ymax = 1.91;
64 Ywid = 0.026;
65 avm1 = 4.0;
66 double HisPDFtmp[62][62] = {
67 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
68 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
69 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
70 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
71 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
72 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
73 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
74 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
75 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
76 { 0, 0, 0, 0, 0, 0, 0,
77 0, 0, 0, 0, 0, 0, 0,
78 0, 0.0139227, 0.0242858, 0.0102631, 0.019465, 0.00821469, 0.00941803,
79 0.0196361, 0.00814691, 0.0270675, 0.0188643, 0.0186632, 0.0107977, 0,
80 0, 0, 0, 0, 0, 0, 0,
81 0, 0, 0, 0, 0, 0, 0,
82 0, 0, 0, 0, 0, 0, 0,
83 0, 0, 0, 0, 0, 0, 0,
84 0, 0, 0, 0, 0, 0 },
85 { 0, 0, 0, 0, 0, 0,
86 0, 0, 0, 0, 0, 0,
87 0, 0.0140485, -0.00172283, 0, 0, 0,
88 0, 0, 0.00962505, 0, 0, -0.000873327,
89 0, 0.0188526, 0.00761928, 0.0154787, 0, 0.0088296,
90 0.00889335, 0.0150876, 0.0114478, 0.0180075, 0.0107661, 0.00843499,
91 0.0130738, 0.0455281, 0, 0, 0, 0,
92 0, 0, 0, 0, 0, 0,
93 0, 0, 0, 0, 0, 0,
94 0, 0, 0, 0, 0, 0,
95 0, 0 },
96 { 0,
97 0,
98 0,
99 0,
100 0,
101 0,
102 0,
103 0,
104 0,
105 0,
106 0.0134615,
107 0,
108 0,
109 0,
110 0,
111 0.0272361,
112 0.00927996,
113 0.00887302,
114 0.00990029,
115 0.0113715,
116 0.00920733,
117 -0.000819338,
118 0.0132542,
119 0.0195743,
120 0.00879929,
121 0.0218734,
122 0.00831937,
123 0,
124 0.0182388,
125 -0.000808949,
126 -0.000794811,
127 0.0166262,
128 0,
129 0.0184469,
130 -0.000680394,
131 0.00879748,
132 0,
133 -0.00190808,
134 0.00985711,
135 -0.00151276,
136 0.0136869,
137 -0.00110282,
138 0,
139 0,
140 0,
141 0,
142 0,
143 0,
144 0,
145 0,
146 0,
147 0,
148 0,
149 0,
150 0,
151 0,
152 0,
153 0,
154 0,
155 0,
156 0,
157 0 },
158 { 0,
159 0,
160 0,
161 0,
162 0,
163 0,
164 0,
165 0,
166 -0.00224346,
167 -0.00158143,
168 0,
169 0.0111645,
170 0.00963723,
171 0.00837308,
172 0,
173 0,
174 -0.00162554,
175 -0.000811716,
176 -0.00103891,
177 0,
178 0.00906974,
179 -0.00080413,
180 -0.000816707,
181 0,
182 0,
183 0.00980291,
184 -0.00183136,
185 0.0101183,
186 -0.000930091,
187 0,
188 0.0185092,
189 0.0157808,
190 -0.000819735,
191 0.0196731,
192 -0.000909308,
193 0,
194 0.0084228,
195 0,
196 0,
197 -0.00270583,
198 -0.000871727,
199 0,
200 -0.000823917,
201 0,
202 0,
203 -0.00290892,
204 0,
205 0,
206 0,
207 0,
208 0,
209 0,
210 0,
211 0,
212 0,
213 0,
214 0,
215 0,
216 0,
217 0,
218 0,
219 0 },
220 { 0,
221 0,
222 0,
223 0,
224 0,
225 0,
226 0,
227 0,
228 0,
229 0,
230 0,
231 0,
232 0,
233 0,
234 -0.000883683,
235 0,
236 0,
237 0,
238 -0.00099019,
239 0,
240 0,
241 -0.000890028,
242 0.0107049,
243 0.00887844,
244 0,
245 0.00868849,
246 0.0078208,
247 0,
248 -0.00102005,
249 0.0107161,
250 0,
251 0.00602897,
252 0,
253 0.028913,
254 -0.00101468,
255 0,
256 0.00864125,
257 0,
258 0,
259 -0.000888573,
260 -0.000867477,
261 -0.000792421,
262 0.00912084,
263 0.0171944,
264 -0.000897794,
265 0,
266 0.0104999,
267 0,
268 0,
269 -0.00449721,
270 0,
271 0,
272 0,
273 0,
274 0,
275 0,
276 0,
277 0,
278 0,
279 0,
280 0,
281 0 },
282 { 0,
283 0,
284 0,
285 0,
286 0,
287 0,
288 0,
289 -0.00106087,
290 0,
291 0,
292 0,
293 0,
294 0,
295 0,
296 -0.00102738,
297 0,
298 0.0200809,
299 0,
300 0,
301 -0.000936822,
302 -0.00102728,
303 0,
304 0,
305 -0.000944456,
306 0,
307 -0.000737237,
308 0.00781316,
309 0,
310 0,
311 0,
312 -0.000795664,
313 0,
314 -0.00140018,
315 0,
316 -0.000874954,
317 0.00909987,
318 -0.00215512,
319 -0.000967849,
320 -0.000805499,
321 0,
322 0,
323 0.00965563,
324 -0.0019969,
325 0,
326 0.00854241,
327 0,
328 0.0407456,
329 0.00778027,
330 0.00917508,
331 0.0078959,
332 0.0118948,
333 0,
334 0,
335 0,
336 0,
337 0,
338 0,
339 0,
340 0,
341 0,
342 0,
343 0 },
344 { 0,
345 0,
346 0,
347 0,
348 0,
349 0,
350 0,
351 -0.00101511,
352 0,
353 0,
354 0,
355 0.00886395,
356 0.00957819,
357 0,
358 0,
359 0,
360 0.0184697,
361 0,
362 -0.00101187,
363 0,
364 -0.00172169,
365 0.00763349,
366 -0.000895323,
367 0,
368 -0.00252987,
369 0.019748,
370 0.00947828,
371 0,
372 0.0192401,
373 0.00849129,
374 -0.000888861,
375 0.01189,
376 0.00896836,
377 0.00997021,
378 0.00931974,
379 0,
380 -0.00170736,
381 0,
382 0,
383 -0.000825473,
384 0,
385 0.00736716,
386 0.00858116,
387 0.0414927,
388 -0.000862929,
389 0.0108653,
390 0.00711648,
391 0,
392 0.0424946,
393 0.02056,
394 0.00888727,
395 0.0183671,
396 0.0300872,
397 0,
398 0,
399 0,
400 0,
401 0,
402 0,
403 0,
404 0,
405 0 },
406 { 0,
407 0,
408 0,
409 0,
410 0,
411 0,
412 0,
413 -0.00121692,
414 0,
415 0,
416 0,
417 0,
418 0.00980403,
419 0,
420 0,
421 0,
422 -0.000906888,
423 0.00818758,
424 -0.000931843,
425 0,
426 -0.000784871,
427 0,
428 -0.00103353,
429 0,
430 0,
431 -0.000859876,
432 0.00731771,
433 -0.00067335,
434 0,
435 -0.000905807,
436 -0.000790764,
437 0.00889752,
438 -0.00070466,
439 0,
440 0.0190332,
441 0,
442 0,
443 0.00846915,
444 0.0110697,
445 0.00897852,
446 0.00804933,
447 0.0201357,
448 0.00851472,
449 0,
450 0.0364704,
451 0.0211352,
452 0.0120208,
453 0.0163991,
454 0,
455 0.0236853,
456 0.0129102,
457 0.00920625,
458 -0.00154954,
459 0.0123733,
460 0.0351667,
461 0,
462 0,
463 0,
464 0,
465 0,
466 0,
467 0 },
468 { 0,
469 0,
470 0,
471 0,
472 0,
473 0,
474 0,
475 0,
476 -0.000847885,
477 0,
478 0,
479 0,
480 0,
481 0,
482 0,
483 0,
484 0.00934519,
485 0.00877748,
486 -0.00110685,
487 0.0100934,
488 0,
489 0.0087101,
490 0.0101893,
491 0.00770844,
492 0,
493 0,
494 0.01073,
495 0.00907416,
496 -0.00108313,
497 -0.000690736,
498 0,
499 0.0198926,
500 0,
501 -0.000960642,
502 0,
503 0.0114602,
504 -0.000867882,
505 0,
506 -0.0039687,
507 0.0108437,
508 0,
509 0,
510 -0.00214432,
511 0.0335963,
512 0.00767771,
513 0.00919692,
514 0.0221316,
515 0.0167495,
516 0.00842359,
517 0,
518 -0.000850634,
519 0.0139215,
520 0.010582,
521 0.00922736,
522 0.00934376,
523 0.0640697,
524 0,
525 0,
526 0,
527 0,
528 0,
529 0 },
530 { 0,
531 0,
532 0,
533 0,
534 -0.00139959,
535 0,
536 0,
537 0,
538 0.0312581,
539 0.0111713,
540 0.00825327,
541 -0.000705715,
542 0.00968094,
543 0.0079305,
544 0,
545 0,
546 0.0118746,
547 0,
548 0,
549 0.00892758,
550 0,
551 0.0102401,
552 0.0193089,
553 -0.000782156,
554 -0.000885337,
555 0,
556 0,
557 0.00795678,
558 0,
559 -0.000709818,
560 0,
561 0.0155405,
562 0.00799595,
563 0,
564 -0.00111399,
565 0.0105511,
566 0.0191352,
567 0,
568 0.00911081,
569 0.00939029,
570 0.0280986,
571 0.0187336,
572 -0.000793975,
573 0,
574 0,
575 0.0082419,
576 0.0074632,
577 0.009078,
578 0.0357759,
579 0.0102688,
580 0.0155454,
581 0.0498312,
582 0.017454,
583 0.0439061,
584 0.0371816,
585 0.036918,
586 -0.00125509,
587 0,
588 0,
589 0,
590 0,
591 0 },
592 { 0,
593 0,
594 0,
595 0,
596 0.0101589,
597 0,
598 0.0395604,
599 0.0114227,
600 -0.00100364,
601 0.0213129,
602 0.0228551,
603 0,
604 0,
605 0,
606 0.0104772,
607 0,
608 0,
609 0,
610 -0.00091241,
611 0,
612 -0.0026225,
613 0.0163946,
614 -0.00064703,
615 -0.000948443,
616 -0.000943019,
617 0,
618 0,
619 0.00841675,
620 0,
621 0.00975579,
622 0,
623 0,
624 0,
625 -0.00113992,
626 0,
627 0,
628 -0.000855679,
629 0.0107493,
630 0.0101579,
631 0.00928391,
632 0.00901179,
633 0,
634 0.0129731,
635 0.0198509,
636 0.0237226,
637 0.0195224,
638 0.0268942,
639 0.0243642,
640 0.0211605,
641 0.0359483,
642 0.0801728,
643 0.0314633,
644 0.0332897,
645 0.0300546,
646 0.0429153,
647 0.0455037,
648 0.0694133,
649 0.0888124,
650 0,
651 0,
652 0,
653 0 },
654 { 0,
655 0,
656 0,
657 0.0489109,
658 0.00944428,
659 0,
660 0.0173018,
661 0.0162537,
662 -0.00171801,
663 0.00829744,
664 0.0160517,
665 0,
666 0.0108203,
667 0.023278,
668 0.00844664,
669 0,
670 0,
671 -0.00106547,
672 0.00857296,
673 0,
674 0,
675 -0.000827661,
676 -0.000848852,
677 0,
678 -0.000990982,
679 0.00832698,
680 -0.00100791,
681 -0.00111151,
682 0,
683 0,
684 0.00859719,
685 0,
686 -0.000750289,
687 0.0095587,
688 0,
689 0.00928949,
690 0.0111698,
691 0.00905906,
692 0.0154512,
693 -0.000941774,
694 0.015085,
695 0.00894646,
696 0.0198978,
697 0.0191158,
698 0.0114518,
699 0.0109502,
700 0.0208813,
701 0.0151289,
702 0.0259387,
703 0.0446887,
704 0.056272,
705 0.0624139,
706 0.075904,
707 0.0216623,
708 0.092942,
709 0.0574898,
710 0.0116131,
711 0.0468348,
712 0.042885,
713 0,
714 0,
715 0 },
716 { 0,
717 0,
718 0,
719 0.0825589,
720 0.0221408,
721 0.0392638,
722 0.0372548,
723 0,
724 0.0281633,
725 0.0218577,
726 0,
727 0.0384219,
728 0.0113054,
729 0,
730 0,
731 0.0101912,
732 0.00908258,
733 0,
734 -0.00165883,
735 -0.000676169,
736 0,
737 -0.000752131,
738 0,
739 0,
740 0,
741 0,
742 -0.000812433,
743 0,
744 0,
745 -0.000921229,
746 0,
747 -0.00102465,
748 -0.000717399,
749 0,
750 -0.00123682,
751 0,
752 0.00834694,
753 -0.00118759,
754 0.00792323,
755 0.0186343,
756 0,
757 0.0301626,
758 0.0097278,
759 0,
760 0.018005,
761 0.0310391,
762 0.0104501,
763 0.0496578,
764 0.0281955,
765 0.0165578,
766 0.0455259,
767 0.035723,
768 0.0701144,
769 0.0881308,
770 0.0787734,
771 0.144527,
772 0.0778648,
773 0.0562271,
774 0.0186825,
775 0,
776 0,
777 0 },
778 { 0,
779 0,
780 0,
781 0.0698628,
782 0.16651,
783 0.105634,
784 0.0477428,
785 0.06,
786 0.0159642,
787 0.048638,
788 0.0439435,
789 0.0187432,
790 -0.000797068,
791 -0.000719579,
792 0.0207966,
793 0.0185039,
794 -0.00211158,
795 0.00953271,
796 0.0212814,
797 0,
798 0.0107911,
799 0,
800 0.0104335,
801 0.00800181,
802 0.00937909,
803 0,
804 0,
805 0,
806 0,
807 0,
808 0,
809 0.0102564,
810 0.0161653,
811 -0.000845697,
812 0,
813 0.00823757,
814 0.0289298,
815 0.0598826,
816 0.0111307,
817 0.00817876,
818 0.0318672,
819 0.0770479,
820 0.0260889,
821 0.0568727,
822 0.0695652,
823 0.0331282,
824 0.0441718,
825 0.0661653,
826 0.0929879,
827 0.0956522,
828 0.0896045,
829 0.0614806,
830 0.146617,
831 0.0557448,
832 0.0764612,
833 0.15177,
834 0.11039,
835 0.0760611,
836 0.131421,
837 0,
838 0,
839 0 },
840 { 0,
841 0,
842 0,
843 0.0844833,
844 0.0848093,
845 0.0990991,
846 0.0516922,
847 0.0821633,
848 0.0580325,
849 0.0579193,
850 0.057767,
851 0.0835098,
852 0.0265978,
853 0.00721743,
854 -0.000881274,
855 0.0530418,
856 0.0379593,
857 0.00763684,
858 0.010234,
859 -0.000857086,
860 0.0197255,
861 0,
862 0.00623986,
863 0.0290557,
864 0,
865 0.007842,
866 -0.00105588,
867 -0.000780333,
868 0.00930577,
869 0,
870 0,
871 0.00649399,
872 0.00729907,
873 0.00778127,
874 0.0308451,
875 0.0164112,
876 0.0479435,
877 0.00776708,
878 0.0518251,
879 0.0613207,
880 0.0714823,
881 0.07,
882 0.114109,
883 0.0487135,
884 0.0616141,
885 0.102306,
886 0.112254,
887 0.148481,
888 0.103235,
889 0.0919618,
890 0.130377,
891 0.205457,
892 0.116028,
893 0.166836,
894 0.1459,
895 0.128874,
896 0.128328,
897 0.106844,
898 0.111886,
899 0.0611888,
900 0,
901 0 },
902 { 0, 0, 0, 0.0835314,
903 0.0221192, 0.112568, 0.0762611, 0.058406,
904 0.0732344, 0.0692292, 0.055456, 0.0457485,
905 0.0486681, 0.0481529, 0.0340297, 0.0407047,
906 0.0280471, 0.0178287, 0.0172723, -0.000836188,
907 0.00740479, 0.00645997, 0.00971044, 0.00563776,
908 0.00868042, 0, -0.00153558, -0.000762348,
909 -0.00119712, 0.00938967, 0.00671685, -0.000859723,
910 -0.00099193, 0.00906069, 0.0252852, 0.00772808,
911 0, 0.0179303, 0.0499147, 0.0057039,
912 0.0629322, 0.028777, 0.0574818, 0.113005,
913 0.0967862, 0.0470863, 0.127491, 0.0842513,
914 0.195622, 0.0993212, 0.185804, 0.12411,
915 0.146365, 0.170778, 0.197943, 0.259737,
916 0.107905, 0.154536, 0.177368, 0.0925204,
917 0, 0 },
918 { 0,
919 0,
920 0,
921 0.111431,
922 0.0654541,
923 0.0807881,
924 0.0784929,
925 0.00795811,
926 0.0194426,
927 0.0339753,
928 0.0677406,
929 0.0450864,
930 0.0176048,
931 0.0196949,
932 -0.000871355,
933 0.0218863,
934 0.0161006,
935 0.0533049,
936 0.0188992,
937 0.0378347,
938 0.00799579,
939 0.00794461,
940 -0.0009632,
941 0.0086048,
942 0.0195132,
943 0.0211498,
944 -0.000866141,
945 0,
946 -0.00272947,
947 -0.00176606,
948 0.00767947,
949 -0.000982462,
950 0.0104655,
951 0.0193164,
952 -0.000979882,
953 0.0171381,
954 0.00888631,
955 0.00887055,
956 0.00893744,
957 0.0207668,
958 0.0102099,
959 0.0718604,
960 0.052395,
961 0.0580129,
962 0.10178,
963 0.061127,
964 0.0282911,
965 0.0639713,
966 0.0528053,
967 0.142383,
968 0.107398,
969 0.0833401,
970 0.138912,
971 0.101017,
972 0.0642299,
973 0.0739587,
974 0.0799474,
975 0.0678708,
976 0.109511,
977 0,
978 0,
979 0 },
980 { 0,
981 0,
982 0,
983 0.0345705,
984 0.0293617,
985 0.0399925,
986 0.0544398,
987 0.0200117,
988 0.00733173,
989 0.0304692,
990 0.0426771,
991 0.0171076,
992 0.0449947,
993 0.0242359,
994 0.0218323,
995 0.00666036,
996 0.00687741,
997 0.0301849,
998 0.00875201,
999 -0.000985995,
1000 0.00922699,
1001 -0.000853764,
1002 -0.00083159,
1003 0.00878674,
1004 0,
1005 0,
1006 -0.00159877,
1007 0,
1008 -0.000724172,
1009 -0.000823628,
1010 0,
1011 0.0179337,
1012 0.00984429,
1013 0.034952,
1014 0.0105243,
1015 -0.000820197,
1016 0.0176194,
1017 0.0336909,
1018 0.0344162,
1019 0.0406053,
1020 0.0239255,
1021 0.0249303,
1022 0.0207042,
1023 0.0205343,
1024 0.0578339,
1025 0.0451814,
1026 0.0176266,
1027 0.0284754,
1028 0.0276654,
1029 0.0825147,
1030 0.0652409,
1031 0.0672455,
1032 0.119474,
1033 0.062069,
1034 0.0100719,
1035 0.0201064,
1036 0.0101344,
1037 0.0135435,
1038 0.0301031,
1039 0.0376412,
1040 0,
1041 0 },
1042 { 0,
1043 0,
1044 0,
1045 0.0100569,
1046 0.0093066,
1047 0.00717103,
1048 0.0273388,
1049 0.0070312,
1050 0.0170305,
1051 0.01851,
1052 0.0316539,
1053 0.0186476,
1054 0.00801258,
1055 0.0298422,
1056 0.0162819,
1057 0,
1058 0.00823547,
1059 0,
1060 0.0246154,
1061 -0.000997358,
1062 0,
1063 0.0247961,
1064 -0.000903635,
1065 0,
1066 0,
1067 -0.000988488,
1068 0,
1069 0,
1070 0,
1071 0.00784142,
1072 0.00947002,
1073 0,
1074 0,
1075 0,
1076 -0.000766148,
1077 0.0201519,
1078 0.0107802,
1079 0.00888957,
1080 0.0252956,
1081 0.0092967,
1082 0,
1083 0.0165119,
1084 0.0311205,
1085 0.00844633,
1086 0,
1087 0.00759511,
1088 0.00904772,
1089 0.0569211,
1090 0.0219707,
1091 0.0201324,
1092 0.0250272,
1093 0.0178588,
1094 -0.00105352,
1095 0.0476843,
1096 0.0251409,
1097 0.0207664,
1098 0.0307567,
1099 0.0220836,
1100 0.0205404,
1101 0.0465057,
1102 0,
1103 0 },
1104 { 0,
1105 0,
1106 0,
1107 0.00929319,
1108 0,
1109 0.00957557,
1110 -0.000775727,
1111 0,
1112 0.0103202,
1113 0,
1114 0,
1115 0,
1116 0.00601536,
1117 0.0330814,
1118 0.01049,
1119 -0.00177138,
1120 0.0117269,
1121 0,
1122 0.00818093,
1123 0.00919637,
1124 0.00743863,
1125 0.00665334,
1126 0.00949445,
1127 0,
1128 0,
1129 0,
1130 0,
1131 -0.000980402,
1132 0.00944467,
1133 0.0103012,
1134 0,
1135 0.010084,
1136 0.0178529,
1137 0,
1138 0.00904177,
1139 -0.000788797,
1140 0,
1141 0.011637,
1142 0.0174232,
1143 0.00806195,
1144 0.0261961,
1145 0,
1146 0.0184518,
1147 0.0206251,
1148 0,
1149 0.0115893,
1150 0.0170932,
1151 0.0312629,
1152 0.0340296,
1153 0.0237025,
1154 0,
1155 0.0275816,
1156 0.0356968,
1157 0.0179104,
1158 0.0328156,
1159 0.0123295,
1160 0.0637828,
1161 0,
1162 0,
1163 0,
1164 0,
1165 0 },
1166 { 0,
1167 0,
1168 0,
1169 0,
1170 -0.00101025,
1171 -0.000858735,
1172 0.00870872,
1173 0.0202326,
1174 0,
1175 0,
1176 0,
1177 0,
1178 -0.000847364,
1179 0,
1180 0.0389124,
1181 0,
1182 0,
1183 0.0069042,
1184 0.0117152,
1185 0,
1186 0,
1187 0.00749609,
1188 0.0156619,
1189 0.0180838,
1190 0.00901454,
1191 0.00925037,
1192 0,
1193 -0.000838362,
1194 0,
1195 0.0104666,
1196 0,
1197 0.0089198,
1198 -0.00105854,
1199 0,
1200 0,
1201 0.00941842,
1202 0.0175091,
1203 0.0457886,
1204 -0.00157401,
1205 0.0180961,
1206 0.0178693,
1207 0,
1208 0.0258238,
1209 0.0236452,
1210 0.0299774,
1211 0.0208024,
1212 0,
1213 0.0278071,
1214 0.00925399,
1215 0.0437028,
1216 0,
1217 0.0163869,
1218 0.0413868,
1219 0.00910784,
1220 0.0192142,
1221 0.0366326,
1222 0.00841407,
1223 -0.0012014,
1224 0.0166667,
1225 0.0436019,
1226 0,
1227 0 },
1228 { 0,
1229 0,
1230 0,
1231 0.0203152,
1232 0.00949449,
1233 0,
1234 -0.000807642,
1235 0.0100947,
1236 0,
1237 0.00798749,
1238 0,
1239 0.0104772,
1240 -0.000759592,
1241 0.0285344,
1242 -0.000809837,
1243 0.00857392,
1244 0.00926016,
1245 0,
1246 0,
1247 0,
1248 0.0085656,
1249 0,
1250 -0.00090472,
1251 0.0101704,
1252 0.00903735,
1253 -0.00268799,
1254 0,
1255 -0.000837338,
1256 0.00714379,
1257 0.0190639,
1258 0.0199808,
1259 0.0147623,
1260 0.0104001,
1261 0.0297673,
1262 0.0588285,
1263 0,
1264 0.0238077,
1265 0.0195378,
1266 0.0230434,
1267 0.0106557,
1268 0.0297837,
1269 0.0146021,
1270 0,
1271 0.00961409,
1272 0.0100678,
1273 0.0105691,
1274 0.0106211,
1275 0.00962249,
1276 -0.00119916,
1277 0.0276335,
1278 0.0223966,
1279 0,
1280 0.0306002,
1281 0.0202501,
1282 -0.0010004,
1283 0,
1284 0.0134923,
1285 0.0132922,
1286 0,
1287 0,
1288 0,
1289 0 },
1290 { 0,
1291 0,
1292 0,
1293 0,
1294 0.0164774,
1295 0.00842688,
1296 0.00911396,
1297 0.00949977,
1298 -0.000850423,
1299 0,
1300 0,
1301 0,
1302 -0.00105277,
1303 -0.000862762,
1304 0.0105271,
1305 -0.000807488,
1306 0,
1307 0.00810359,
1308 0,
1309 0,
1310 0,
1311 0.0101164,
1312 0,
1313 0,
1314 0,
1315 0,
1316 0.00842445,
1317 0.0114901,
1318 0.0266553,
1319 0,
1320 0,
1321 0.00915325,
1322 0.0164874,
1323 0.00967876,
1324 0,
1325 0,
1326 0.023019,
1327 0.00997442,
1328 0.0181452,
1329 0.0352483,
1330 0.0235438,
1331 0.00947753,
1332 0.0210715,
1333 0.010443,
1334 0.00770524,
1335 0.0121609,
1336 0.0179744,
1337 0.0279225,
1338 0.0380326,
1339 0.0326071,
1340 0.0259083,
1341 0,
1342 0.0476909,
1343 0.0203499,
1344 0,
1345 0.0108538,
1346 0.0438583,
1347 0.0124589,
1348 0.0649872,
1349 0,
1350 0,
1351 0 },
1352 { 0,
1353 0,
1354 0,
1355 0,
1356 0,
1357 -0.000677822,
1358 0,
1359 -0.00196125,
1360 -0.00170722,
1361 0.00917402,
1362 0.0168193,
1363 0.0160035,
1364 0.0187455,
1365 0,
1366 0,
1367 0.0105069,
1368 0.00653315,
1369 0.00946888,
1370 0.00939433,
1371 -0.000951599,
1372 0.00793651,
1373 0,
1374 0,
1375 0.00800202,
1376 0,
1377 0,
1378 0,
1379 0,
1380 0.0240316,
1381 0,
1382 -0.000909594,
1383 0,
1384 -0.000833545,
1385 0.00925294,
1386 0.0407852,
1387 0.00822491,
1388 0,
1389 0,
1390 -0.000796518,
1391 0.0428913,
1392 0,
1393 0,
1394 0.00864616,
1395 0.0081278,
1396 0.0198069,
1397 -0.000900054,
1398 0.0129539,
1399 0.0291967,
1400 -0.00107986,
1401 0.0204132,
1402 0.0175669,
1403 0.00921452,
1404 0.020746,
1405 0.00976493,
1406 0,
1407 0,
1408 -0.00179696,
1409 0,
1410 0,
1411 0.0412541,
1412 0,
1413 0 },
1414 { 0,
1415 0,
1416 0,
1417 0.0203955,
1418 0.0106986,
1419 0.00881007,
1420 0.0191054,
1421 0,
1422 0,
1423 0.00945626,
1424 0.0104001,
1425 0.0222121,
1426 0.0085498,
1427 0.00947087,
1428 0.00954733,
1429 0,
1430 0.0116447,
1431 0,
1432 0,
1433 0,
1434 0.0112924,
1435 0,
1436 0.00894737,
1437 -0.000928966,
1438 -0.000824152,
1439 0,
1440 0.0169718,
1441 0.00829953,
1442 -0.0008101,
1443 0.0087907,
1444 0.0212713,
1445 0.00962365,
1446 0.00818169,
1447 0.0186238,
1448 0.0443356,
1449 0.0088327,
1450 0.00907982,
1451 0.0092243,
1452 -0.00154134,
1453 0.00745146,
1454 0.0181104,
1455 0,
1456 0.0174172,
1457 0,
1458 0.0188855,
1459 0.0111528,
1460 0.0104965,
1461 0.0187198,
1462 0.00907934,
1463 0,
1464 0.0101033,
1465 0.0132822,
1466 0,
1467 0.0255092,
1468 0.0117876,
1469 0.00884271,
1470 0.0388322,
1471 0.0348326,
1472 0.0268515,
1473 0,
1474 0,
1475 0 },
1476 { 0,
1477 0,
1478 0,
1479 0.0126706,
1480 0.0172168,
1481 0.0209784,
1482 0.0115417,
1483 0.0198827,
1484 0,
1485 -0.000613454,
1486 0,
1487 0,
1488 -0.0012482,
1489 0.00666148,
1490 0,
1491 0,
1492 -0.000914029,
1493 -0.000802899,
1494 0.00751363,
1495 0,
1496 0.0210709,
1497 0,
1498 0.00956299,
1499 0.00755021,
1500 0.0108458,
1501 0,
1502 0.00851648,
1503 0.00970513,
1504 0.00984537,
1505 0.0114522,
1506 0.00811426,
1507 0.014856,
1508 0,
1509 0,
1510 0,
1511 0.0356321,
1512 0.00915935,
1513 0.0079448,
1514 0.0112329,
1515 0,
1516 0.0195858,
1517 0.0361888,
1518 0.0197721,
1519 0.00967225,
1520 0.0190026,
1521 0,
1522 0.0116118,
1523 0.0173056,
1524 0.0256392,
1525 0,
1526 0,
1527 0.0101425,
1528 0,
1529 0.00986239,
1530 0.0204131,
1531 0.0313983,
1532 0.0328703,
1533 0.0150389,
1534 0,
1535 0,
1536 0,
1537 0 },
1538 { 0,
1539 0,
1540 0,
1541 0,
1542 0,
1543 0,
1544 0.00925089,
1545 -0.000860052,
1546 0.0113342,
1547 0.0200476,
1548 0,
1549 0,
1550 0.0109967,
1551 0,
1552 0.0092987,
1553 0.00817958,
1554 0,
1555 0,
1556 -0.000883344,
1557 0.0114325,
1558 0.0287869,
1559 0.0182725,
1560 0.00878291,
1561 0.0194701,
1562 -0.000979345,
1563 0.0103969,
1564 0,
1565 0.008714,
1566 0,
1567 0.0102515,
1568 -0.000783578,
1569 0.00882312,
1570 0.0181932,
1571 0,
1572 -0.00082272,
1573 0.0251304,
1574 0.0243868,
1575 0,
1576 0.0176539,
1577 0.0276354,
1578 0.011233,
1579 0.0169369,
1580 0,
1581 0,
1582 0,
1583 -0.000914585,
1584 0,
1585 0.0108333,
1586 0.0105352,
1587 0.0189954,
1588 0.026487,
1589 0.0219047,
1590 0.0214605,
1591 0,
1592 0.0117922,
1593 0.0112767,
1594 0.0208396,
1595 0.0206242,
1596 0,
1597 0,
1598 0,
1599 0 },
1600 { 0, 0,
1601 0, 0.0215772,
1602 0, 0.0469387,
1603 -0.00199938, 0,
1604 0.00806495, 0.0113939,
1605 0.00784219, 0.0102418,
1606 0.00848231, 0,
1607 0, -0.00238053,
1608 0, 0,
1609 0.0108451, 0.0224638,
1610 0.00895348, 0.0104752,
1611 0.00931011, 0,
1612 0.00806999, -0.000835198,
1613 0.00828157, 0,
1614 0, 0.010863,
1615 0, 0.00966184,
1616 0, 0.00880196,
1617 0, 0,
1618 0.00967279, 0,
1619 -0.00104673, 0,
1620 0.0193632, 0.0294464,
1621 0.0207633, 0.0180889,
1622 0, 0.0196807,
1623 0.0113285, 0.0207271,
1624 0.00870174, 0.0122184,
1625 0.0084899, 0,
1626 0.0191684, 0.0251521,
1627 0.0244355, 0.0117797,
1628 0.0433962, 0.107924,
1629 0, 0,
1630 0, 0 },
1631 { 0, 0,
1632 0, 0,
1633 0.0242275, 0,
1634 0, 0.0096565,
1635 -0.00155124, 0.0208565,
1636 0, 0,
1637 0.0294917, -0.00180637,
1638 0.00820202, -0.000719485,
1639 0.017647, 0.00939541,
1640 0.00995025, 0,
1641 0.031661, 0.0179949,
1642 0.0088753, 0.00956716,
1643 0, 0.00942728,
1644 0.0182113, 0.019223,
1645 0.0100184, 0,
1646 0.00826353, 0.0109375,
1647 0, 0.0095989,
1648 0.0088364, 0,
1649 0.0177186, 0.0176286,
1650 0, 0.00918478,
1651 0.00979515, 0.0718151,
1652 0.00712146, 0.0388281,
1653 0, 0.0169886,
1654 0.027248, -0.000841679,
1655 0, 0.011249,
1656 0.010187, 0,
1657 0, 0,
1658 0.064453, 0.0347804,
1659 0.14096, -0.0104849,
1660 0, 0,
1661 0, 0 },
1662 { 0,
1663 0,
1664 0,
1665 0.0304846,
1666 0.00895394,
1667 0,
1668 0,
1669 0.0095143,
1670 -0.000817347,
1671 0.0100501,
1672 0.00882872,
1673 0,
1674 0.00960198,
1675 0.00871875,
1676 0.024524,
1677 0.0173022,
1678 0,
1679 0.0126656,
1680 0.0101352,
1681 0.0331152,
1682 0.0228368,
1683 0.0179239,
1684 0,
1685 0.00906672,
1686 0.0325413,
1687 0,
1688 0.00905138,
1689 0,
1690 0.00894021,
1691 0.0172658,
1692 0,
1693 0,
1694 0.00862803,
1695 0,
1696 0.0324659,
1697 0,
1698 0.0155981,
1699 0.0206338,
1700 0.0230017,
1701 0.0108395,
1702 -0.000768815,
1703 0.0380435,
1704 0.00846202,
1705 0.0431805,
1706 0.0186941,
1707 0.0185413,
1708 0.0203616,
1709 0.0412492,
1710 0.0285326,
1711 0.0174774,
1712 0.0220788,
1713 0.00855435,
1714 0.0185348,
1715 0.0155923,
1716 0.0847581,
1717 0.17733,
1718 0.384751,
1719 0,
1720 0,
1721 0,
1722 0,
1723 0 },
1724 { 0, 0,
1725 0, 0.0549451,
1726 0, -0.000809818,
1727 0.0341507, 0.00924517,
1728 0, 0,
1729 0, 0.025915,
1730 0, 0,
1731 0, 0.0100309,
1732 0, 0.00773839,
1733 0, -0.000922284,
1734 0.00824142, 0.00998856,
1735 0.0185869, 0,
1736 0, -0.000930184,
1737 0.00774234, 0.0271866,
1738 0.0261441, 0.00859834,
1739 0.0101072, 0.0106873,
1740 0.0511486, 0.00915786,
1741 0.00992035, 0.00784197,
1742 0.0202428, 0.00991919,
1743 0.0294389, 0,
1744 0.00836111, 0.0168758,
1745 0.0119083, 0.0087859,
1746 0, 0.0168192,
1747 0.00894654, 0,
1748 0.0111521, 0,
1749 0, 0.020193,
1750 0.00745592, 0.043302,
1751 0.33012, 0.756053,
1752 4, 0,
1753 0, 0,
1754 0, 0 },
1755 { 0, 0,
1756 0, 0,
1757 0.0193347, 0,
1758 0, 0.00877543,
1759 0.00926148, -0.000978804,
1760 0, 0,
1761 0.00616838, 0.0480958,
1762 -0.00110401, 0,
1763 -0.00094806, 0,
1764 0.0266272, 0.00926523,
1765 0, 0.0470984,
1766 0.0316436, 0.0102558,
1767 0.0137878, 0.00956818,
1768 0, 0.0109435,
1769 0.00795692, 0.00996638,
1770 0, 0.0262011,
1771 0, 0,
1772 0, 0.00821991,
1773 0.00822732, 0.0178063,
1774 0, 0.0131448,
1775 0.0187651, 0.0376607,
1776 0, 0.00832224,
1777 0.035451, 0.0111331,
1778 0.00716381, 0,
1779 0.0200805, 0,
1780 0.0247405, 0.0343828,
1781 0.073858, 0.162882,
1782 0.753372, 0.877119,
1783 0, 0,
1784 0, 0,
1785 0, 0 },
1786 { 0,
1787 0,
1788 0,
1789 0,
1790 0,
1791 0.00892485,
1792 0.0087992,
1793 0.00866911,
1794 0,
1795 0,
1796 0.0183706,
1797 0.00892193,
1798 0.0295628,
1799 0,
1800 0.029044,
1801 0,
1802 0.0264391,
1803 0.0130717,
1804 -0.000957429,
1805 -0.00102977,
1806 0,
1807 0.0217497,
1808 0,
1809 0.0119743,
1810 0,
1811 0.00814037,
1812 -0.00216151,
1813 0,
1814 0.00901897,
1815 0.0800344,
1816 0.0311763,
1817 0,
1818 -0.000999707,
1819 0.0175548,
1820 0.00776979,
1821 0.00997534,
1822 0.0306324,
1823 0.0163388,
1824 0.00966387,
1825 0,
1826 0.0313085,
1827 0.0438893,
1828 0.0189919,
1829 0.0370669,
1830 0.032723,
1831 0.0253304,
1832 0.0216285,
1833 -0.000945404,
1834 0.0126512,
1835 0.0105467,
1836 0.0114826,
1837 0.0474927,
1838 0.229064,
1839 0.468742,
1840 0.499709,
1841 0,
1842 0,
1843 0,
1844 0,
1845 0,
1846 0,
1847 0 },
1848 { 0, 0, 0, 0, 0, 0.0102398, 0,
1849 0.00963375, 0.0256482, 0.0205175, 0, 0.00772904, 0.0187665, 0.00648284,
1850 0.0302874, 0.0335031, 0, 0.0321651, 0.0210876, 0, 0.0164572,
1851 0.00975459, 0.0144711, 0.0349518, 0.0103055, 0, 0.00774833, 0,
1852 0, 0, 0.00881944, 0.00900243, 0.0166373, 0.0190377, 0.0121791,
1853 0.0290389, 0.0218861, 0.0282682, 0.00957256, 0.0210033, 0.0188598, 0.0160125,
1854 -0.00114259, 0, 0, 0.0232449, 0.0218724, 0.0136043, 0.00759345,
1855 0.0100437, 0.0400415, 0.129409, 0.382832, 0.499742, 0.169623, 0,
1856 0, 0, 0, 0, 0, 0 },
1857 { 0, 0, 0,
1858 0, 0.0217826, 0.00954403,
1859 0.00978454, 0.0314549, 0.00763196,
1860 0.0107647, 0.0477279, 0.049444,
1861 0.00818761, 0.0191171, 0.0187272,
1862 0.00827722, 0.00809931, 0.00807353,
1863 0.00865924, 0.0440351, 0.0188898,
1864 0.0178678, 0.0187399, 0,
1865 0.0353044, 0.0199281, 0,
1866 0.0092765, 0, 0.0341788,
1867 0.0137857, 0, -0.000917367,
1868 0, 0.0222421, 0.0175874,
1869 0.0179275, 0.0111494, 0.0397602,
1870 0.0228989, 0.0195081, 0.0173706,
1871 0, 0.00988815, 0.00821166,
1872 0.0105519, 0, 0.0425583,
1873 0.0290794, 0.0377965, 0.144596,
1874 0.31633, 0.309859, 0.407018,
1875 0, 0, 0,
1876 0, 0, 0,
1877 0, 0 },
1878 { 0, 0, 0, 0, 0.0373188, -0.000800659,
1879 0, 0.0206369, 0.0201636, 0.024405, -0.00164074, 0.0322904,
1880 0.00772539, 0.017094, 0, 0.0102001, 0.0170068, 0.00789728,
1881 0.0285433, 0.017654, 0.0209081, 0.0353514, 0.00991326, 0,
1882 0, 0.0112883, 0.0272436, 0.0221948, 0, 0.0303313,
1883 0, 0.00892256, 0.0112371, 0.00807134, 0.0216746, 0.0174218,
1884 0, 0, 0.00866308, 0.0088484, 0.0820871, 0.0221963,
1885 0.00996698, 0.00907869, 0.00964537, 0.0253436, 0.00854222, 0.0357143,
1886 0.0486163, 0.0954776, 0.259473, 0.149898, 0.052, 0,
1887 0, 0, 0, 0, 0, 0,
1888 0, 0 },
1889 { 0,
1890 0,
1891 0,
1892 0,
1893 0,
1894 0.00925421,
1895 -0.000826828,
1896 0.0114796,
1897 0.0186418,
1898 -0.000892872,
1899 0,
1900 0.0120563,
1901 -0.000841373,
1902 0.0396193,
1903 0.00921667,
1904 0,
1905 -0.000839938,
1906 0,
1907 0.0216088,
1908 0,
1909 0.00909285,
1910 0.00989559,
1911 0.0187188,
1912 0.0143512,
1913 0.0166337,
1914 0,
1915 0,
1916 0.0433543,
1917 0,
1918 0.00886178,
1919 0,
1920 0.0337719,
1921 0.00869666,
1922 0.00954008,
1923 0.0229141,
1924 0.00920752,
1925 0.00859955,
1926 0.016442,
1927 0.0264347,
1928 0.016196,
1929 0.0310447,
1930 0,
1931 0.0232209,
1932 0.0200668,
1933 0.0194759,
1934 0.0252697,
1935 0.00791894,
1936 0.0218616,
1937 0.0886278,
1938 0.107202,
1939 0.120529,
1940 0.017576,
1941 0,
1942 0,
1943 0,
1944 0,
1945 0,
1946 0,
1947 0,
1948 0,
1949 0,
1950 0 },
1951 { 0, 0, 0, 0, 0, 0.0176408, 0.0402815,
1952 0.00900193, 0.00910931, 0, 0.00995597, 0, 0.0193276, 0.02625,
1953 0.0102957, 0.0175661, 0.00946315, 0.0235473, 0.0224084, 0.0174483, 0.0167777,
1954 0.0246835, 0.0131795, 0.0464314, 0.018863, 0.00888001, 0.019547, 0.00803371,
1955 0.0103648, 0.00857505, 0.0340203, 0, 0.0257845, 0, 0.020093,
1956 0.0185917, 0.009653, 0.0361675, 0.0340909, 0.0165112, 0.0213852, -0.000725483,
1957 0.00793697, 0.00797994, 0.0197249, 0.0117026, 0.0102868, 0.0365747, 0.124825,
1958 0.0937213, 0.0546665, 0, 0, 0, 0, 0,
1959 0, 0, 0, 0, 0, 0 },
1960 { 0, 0, 0, 0, 0, 0, 0.0158892,
1961 0.00842231, 0.00791289, 0.0107855, 0.0479872, 0, 0, 0.0337553,
1962 0.0166636, 0.0180018, 0.0818542, 0.00921793, 0, 0, 0.0161146,
1963 0.0104829, 0.0112961, 0.00880284, 0, 0.00894671, 0.0321728, 0.00823045,
1964 0.00836536, 0.0116709, 0.0169351, 0, 0.0504936, 0.0272736, 0.00943627,
1965 0.0163528, 0.00930591, 0.0102789, 0, 0.0395859, 0.0206268, 0.00947595,
1966 0.0454126, 0.0282035, 0.0215823, 0.0425532, 0.0490323, 0.0663051, 0.093707,
1967 0.0430499, 0.204082, 0, 0, 0, 0, 0,
1968 0, 0, 0, 0, 0, 0 },
1969 { 0, 0, 0, 0, 0,
1970 0.0474398, 0.010514, 0.0171634, 0.0472313, 0.0225918,
1971 0.00918337, 0, 0.021168, 0.0214649, 0.0295599,
1972 0, 0.0215697, 0.00883459, 0.0191107, 0.0101801,
1973 0.00683388, 0.0288497, 0.0176917, 0.0115658, 0,
1974 0.0205727, 0.0160846, 0.00916362, 0.00781841, 0.00923558,
1975 0.00764855, 0.0115033, 0.0476958, 0.00914876, 0.0340413,
1976 0.0260466, 0, 0.0203144, 0.0224558, 0,
1977 0.00836209, 0.032166, 0, -0.00201207, -0.000999997,
1978 0.0364455, 0.0197025, 0.0580674, 0.025409, 0.20625,
1979 0, 0, 0, 0, 0,
1980 0, 0, 0, 0, 0,
1981 0, 0 },
1982 { 0, 0, 0, 0, 0, 0, 0.0204216,
1983 0.0178749, 0.0339385, 0.0476663, 0.0214505, 0, 0.010011, 0.00825038,
1984 0, 0.0106965, 0.0298976, 0.0240339, 0.00813984, 0.0264443, 0,
1985 0.0377711, 0.0081673, 0.0264415, 0, -0.00081143, 0.00870287, 0.0430039,
1986 0.00920434, 0.00927802, 0.0102824, 0.00919913, 0.0321815, 0.00880137, 0.0167547,
1987 0.0141489, 0.0282248, 0.0100758, 0.0192092, 0.00911131, 0.0191987, 0.00988181,
1988 0.0094089, 0.037713, 0.0211831, 0, 0.0403058, 0.0343354, 0.0401235,
1989 0, 0, 0, 0, 0, 0, 0,
1990 0, 0, 0, 0, 0, 0 },
1991 { 0, 0, 0, 0, 0, 0, 0.0644633,
1992 0.0346, 0.0371551, 0, 0.0102819, 0.00970625, 0.0103222, 0.010594,
1993 0.00810723, 0.047961, 0.102065, 0.0250843, 0.0208742, 0, 0.0254557,
1994 0, 0.0363382, 0.0186905, 0, 0.008927, 0.0300059, 0.0378217,
1995 0.0508027, 0, 0.0479418, 0.026139, 0.0105877, 0.0322802, 0.0196786,
1996 0.0184422, 0.0167361, 0.0191199, 0.0250955, 0, 0.0225517, 0.0328446,
1997 0.0514563, 0.0103672, 0.0577811, 0.0380105, 0.0258166, 0.0629371, 0,
1998 0, 0, 0, 0, 0, 0, 0,
1999 0, 0, 0, 0, 0, 0 },
2000 { 0, 0, 0, 0, 0, 0, 0.0231183,
2001 0.0173595, 0.0199675, 0, 0.0368378, 0.0089956, 0.00849399, 0.0356221,
2002 0.00735485, 0.00863388, 0.0196979, 0.050897, 0.0251263, 0.0181164, 0.00793501,
2003 0.0370859, 0.00934374, 0.0167843, 0.017893, 0.00876617, 0, 0,
2004 0.0275779, 0.0243458, 0.026663, 0.0332204, 0.0393579, 0.0128048, 0.0185263,
2005 0.00752273, 0.0539919, 0.0286538, 0.017602, 0.0211793, 0.0304102, 0.0377677,
2006 0.0875269, 0.0560436, 0.0186506, 0.0278015, 0.0493056, 0, 0,
2007 0, 0, 0, 0, 0, 0, 0,
2008 0, 0, 0, 0, 0, 0 },
2009 { 0, 0, 0, 0, 0, 0, 0,
2010 0.0160049, 0, 0.0104456, 0.0449646, 0, 0.0203922, 0.0206122,
2011 0.054633, 0.0527813, 0.00926439, 0.0356753, 0.0354575, 0.0316874, 0.00742951,
2012 0.0183808, 0.0190379, 0.0430335, 0.00998737, 0.0183085, 0.0152626, 0.0340017,
2013 0.0507289, 0, 0.0104491, 0.0116689, 0.00989029, 0.05604, 0,
2014 0.010513, 0, 0.0129241, 0.0253465, 0.00932121, 0.062275, 0.0260865,
2015 0.101612, 0.0374271, 0.0134468, 0.0345664, 0, 0, 0,
2016 0, 0, 0, 0, 0, 0, 0,
2017 0, 0, 0, 0, 0, 0 },
2018 { 0,
2019 0,
2020 0,
2021 0,
2022 0,
2023 0,
2024 0,
2025 0.0401984,
2026 0.0239462,
2027 0.0323053,
2028 0.0299613,
2029 0.00808859,
2030 -0.000959072,
2031 0.028936,
2032 0,
2033 0.0360868,
2034 0.0583658,
2035 0.0340565,
2036 0,
2037 0.0185563,
2038 0.0109007,
2039 0,
2040 0.0106144,
2041 0,
2042 0.0538789,
2043 0.0100871,
2044 0.0264126,
2045 0.0181142,
2046 0.00924878,
2047 0.0230294,
2048 0,
2049 0.0237235,
2050 0,
2051 0.0100862,
2052 0.0336292,
2053 0.0208131,
2054 0.042786,
2055 0.0274457,
2056 0.0377862,
2057 0.0349889,
2058 0.072125,
2059 0.146772,
2060 0.0632134,
2061 0.0768693,
2062 0.118474,
2063 0,
2064 0,
2065 0,
2066 0,
2067 0,
2068 0,
2069 0,
2070 0,
2071 0,
2072 0,
2073 0,
2074 0,
2075 0,
2076 0,
2077 0,
2078 0,
2079 0 },
2080 { 0, 0, 0, 0, 0, 0, 0,
2081 0.0465492, 0.0388832, 0.0211017, 0.0219103, 0.0224301, 0.0082513, 0.0485614,
2082 0, 0.0513242, 0, 0.0205069, 0.0318284, 0.0180626, 0.0181801,
2083 0.0312115, 0.00763562, 0.021414, 0.00981205, 0.0211101, 0.0296362, 0,
2084 0.0104593, 0.0100904, 0.0180371, 0.0174462, 0.0241922, 0.0450072, 0.00850443,
2085 0, 0.0186514, 0.0374268, 0.0605321, 0.2021, 0.209671, 0.0577898,
2086 0, 0.0483516, 0, 0, 0, 0, 0,
2087 0, 0, 0, 0, 0, 0, 0,
2088 0, 0, 0, 0, 0, 0 },
2089 { 0, 0, 0, 0, 0, 0, 0, 0,
2090 0.0447205, 0.0353172, 0.0401819, 0.048029, 0.0330707, 0.0328415, 0.0186309, 0.0436886,
2091 0.0277081, 0.0419765, 0.0218716, 0.0285834, 0.0168439, 0.0103291, 0.0385388, 0.021661,
2092 0, 0.0257292, 0.0564093, 0, 0.0214601, 0.0109343, 0, 0.0549645,
2093 0.0105792, 0.0403816, 0.0173734, 0.0182437, 0.0456348, 0.0563464, 0.172808, 0.23266,
2094 0.0694025, 0.0296717, 0.0758808, 0, 0, 0, 0, 0,
2095 0, 0, 0, 0, 0, 0, 0, 0,
2096 0, 0, 0, 0, 0, 0 },
2097 { 0, 0, 0, 0, 0, 0, 0,
2098 0, 0.0234194, 0.00993789, 0.0487651, 0.0231491, 0.0203606, 0.0232817,
2099 0.0218335, 0.0238186, 0.0277103, 0.01239, 0.0181719, 0.0307454, 0.0087968,
2100 0.0112293, 0.02037, 0.00871786, 0.021266, 0.0102832, 0.0190427, 0,
2101 0.0172336, 0.0683699, 0.0352564, 0.00866234, 0.0291314, 0.0111943, 0.0182796,
2102 0.0334819, 0.0346112, 0.313737, 0.389412, 0.0889875, 0.0436609, 0,
2103 0, 0, 0, 0, 0, 0, 0,
2104 0, 0, 0, 0, 0, 0, 0,
2105 0, 0, 0, 0, 0, 0 },
2106 { 0, 0, 0, 0, 0, 0, 0,
2107 0, 0.0985222, 0.0230706, 0.0257939, 0.00966886, 0.0216262, 0.0181988,
2108 0.0174742, 0.0201692, 0.0268839, 0.0126801, 0, 0.0184841, 0,
2109 0.0277268, 0.00982999, 0.0135447, 0.0431554, 0.0553891, 0.0198105, 0.00958762,
2110 0.0308985, 0.0303113, 0.0342678, 0.00921841, 0.0196959, 0.0658576, 0.0100004,
2111 0.193647, 0.28473, 0.542411, 0.204171, 0, 0, 0,
2112 0, 0, 0, 0, 0, 0, 0,
2113 0, 0, 0, 0, 0, 0, 0,
2114 0, 0, 0, 0, 0, 0 },
2115 { 0, 0, 0, 0, 0, 0, 0,
2116 0, 0, 0.03174, 0.0105849, 0.0771186, 0.0279586, 0.0266612,
2117 0.0318761, 0.0559981, 0.0390625, 0.0121277, 0.0320627, 0.0104402, 0.0432972,
2118 -0.00110143, 0.00831363, 0.00758654, 0.027497, 0.0120548, 0.0234965, 0.00981269,
2119 0.0075152, 0.0344042, 0.0575213, 0.0467659, 0.0353293, 0.0774836, 0.121804,
2120 0.280416, 0.482075, 0.287566, 0.0886228, 0, 0, 0,
2121 0, 0, 0, 0, 0, 0, 0,
2122 0, 0, 0, 0, 0, 0, 0,
2123 0, 0, 0, 0, 0, 0 },
2124 { 0, 0, 0, 0, 0, 0, 0,
2125 0, 0, 0.0947666, 0.0470336, 0.0443913, 0.0497682, 0.0266242,
2126 0.0106494, 0, 0.00941109, 0.0125906, 0, 0.0104469, 0.0107968,
2127 0.039149, 0.0413446, 0.0306129, 0.00857951, 0.00911935, 0.0358557, 0,
2128 0.0229856, 0.0215492, 0.0385325, 0.0731537, 0.0260051, 0.15089, 0.420476,
2129 0.753433, 0.444165, 0.286792, 0, 0, 0, 0,
2130 0, 0, 0, 0, 0, 0, 0,
2131 0, 0, 0, 0, 0, 0, 0,
2132 0, 0, 0, 0, 0, 0 },
2133 { 0, 0, 0, 0, 0, 0, 0,
2134 0, 0, 0, 0.12383, 0.0510436, 0.0110302, 0.00967458,
2135 0.0349852, 0.0344665, 0.0105125, 0.0104745, 0.0249594, 0.0187075, 0.00857108,
2136 0.0176199, 0.0384546, 0.0209939, 0.0202804, 0.0102498, 0.018094, 0.017984,
2137 0.0481076, 0.0963358, 0.0563204, 0.0829453, 0.0781776, 0.337149, 0.878505,
2138 0.719033, 0.222222, 0, 0, 0, 0, 0,
2139 0, 0, 0, 0, 0, 0, 0,
2140 0, 0, 0, 0, 0, 0, 0,
2141 0, 0, 0, 0, 0, 0 },
2142 { 0, 0, 0, 0, 0, 0, 0, 0,
2143 0, 0, 0, 0.0588235, 0.0230188, 0.0364726, 0, 0.0119376,
2144 0.0163343, 0.0487118, 0, 0.0367248, 0, 0.0277698, 0.0270057, 0.0212871,
2145 0.0320826, 0, 0.0110496, 0.0271488, 0, 0.0169567, 0.0667761, 0.205392,
2146 0.51, 1.19001, 1.39583, 0, 0, 0, 0, 0,
2147 0, 0, 0, 0, 0, 0, 0, 0,
2148 0, 0, 0, 0, 0, 0, 0, 0,
2149 0, 0, 0, 0, 0, 0 },
2150 { 0, 0, 0, 0, 0, 0, 0, 0,
2151 0, 0, 0, 0.166667, 0.0190713, 0.0228799, 0.0135813, 0,
2152 0.0411091, 0.0132791, 0, 0.0332337, 0.0365666, 0.0218031, 0.0308543, 0.0353304,
2153 0.0472783, 0.0345769, 0.0405673, 0.0216447, 0.0711784, 0.0323605, 0.130971, 0.199812,
2154 0.798479, 2, 0, 0, 0, 0, 0, 0,
2155 0, 0, 0, 0, 0, 0, 0, 0,
2156 0, 0, 0, 0, 0, 0, 0, 0,
2157 0, 0, 0, 0, 0, 0 },
2158 { 0, 0, 0, 0, 0, 0, 0, 0,
2159 0, 0, 0, 0, 0, 0.0334781, 0.0112575, 0.0227968,
2160 0.0277248, 0.0143044, 0.0257028, 0.0131532, 0, 0.013922, 0.0491049, 0.0103723,
2161 0.028983, 0.0598813, 0.0451858, 0.0420784, 0, 0.128964, 0.152036, 0,
2162 0, 0, 0, 0, 0, 0, 0, 0,
2163 0, 0, 0, 0, 0, 0, 0, 0,
2164 0, 0, 0, 0, 0, 0, 0, 0,
2165 0, 0, 0, 0, 0, 0 },
2166 { 0, 0, 0, 0, 0, 0, 0, 0,
2167 0, 0, 0, 0, 0, 0, 0, 0.0162025,
2168 0, 0, 0.0130208, 0, 0.0292517, 0.0137187, 0.0675021, 0.0272235,
2169 0.0277778, 0.0194071, 0.0628173, 0.0737442, 0.0212331, 0.113636, 0, 0,
2170 0, 0, 0, 0, 0, 0, 0, 0,
2171 0, 0, 0, 0, 0, 0, 0, 0,
2172 0, 0, 0, 0, 0, 0, 0, 0,
2173 0, 0, 0, 0, 0, 0 },
2174 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2175 0, 0, 0, 0, 0.226268, 0.051353, 0.0232558, 0, 0.0439695, 0, 0.0276163, 0, 0, 0, 0, 0,
2176 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2177 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
2178 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2179 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2180 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
2181 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2182 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2183 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 } };
2184
2185 for ( int i = 0; i < 62; i++ )
2186 {
2187 for ( int j = 0; j < 62; j++ ) { HisPDF[i][j] = HisPDFtmp[i][j]; }
2188 }
2189}
2190
2191void EvtDtoKKpiPlot::initProbMax() { noProbMax(); }
2192
2193void EvtDtoKKpiPlot::decay( EvtParticle* p ) {
2194
2195loop:
2196 p->initializePhaseSpace( getNDaug(), getDaugs() );
2197
2198 EvtParticle *id1, *id2, *id3;
2199 EvtVector4R pd1, pd2, pd3;
2200 double xmass13, xmass12, xmass23;
2201
2202 id1 = p->getDaug( 0 );
2203 id2 = p->getDaug( 1 );
2204 id3 = p->getDaug( 2 );
2205
2206 pd1 = id1->getP4Lab();
2207 pd2 = id2->getP4Lab();
2208 pd3 = id3->getP4Lab();
2209
2210 // xmass12=(pd1+pd2).mass()*(pd1+pd2).mass(); // M_ksopi
2211 xmass13 = ( pd1 + pd3 ).mass() * ( pd1 + pd3 ).mass(); // M_ksow
2212 xmass23 = ( pd2 + pd3 ).mass() * ( pd2 + pd3 ).mass(); // M_piw
2213
2214 int xbin = FindXBin( xmass23 );
2215 int ybin = FindYBin( xmass13 );
2216
2217 double xratio12 = HisPDF[xbin][ybin] / avm1;
2218 if ( xratio12 <= 0 ) goto loop;
2219
2220 double rd12 = EvtRandom::Flat( 0.0, 1.0 );
2221
2222 if ( rd12 > xratio12 ) goto loop;
2223
2224 return;
2225}
2226
2227int EvtDtoKKpiPlot::FindXBin( double mass2 ) {
2228 if ( mass2 < Xmin ) { return 0; }
2229 else if ( mass2 >= Xmax ) { return 61; }
2230 else { return int( ( mass2 - Xmin ) / Xwid ) + 1; }
2231}
2232
2233int EvtDtoKKpiPlot::FindYBin( double mass2 ) {
2234 if ( mass2 < Ymin ) { return 0; }
2235 else if ( mass2 >= Ymax ) { return 61; }
2236 else { return int( ( mass2 - Ymin ) / Ywid ) + 1; }
2237}
mass
double mass
Definition CosmicGenerator.cxx:128
EvtDecayBase::getParentId
EvtId getParentId()
Definition EvtDecayBase.hh:63
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449
EvtDtoKKpiPlot::FindXBin
int FindXBin(double mass2)
Definition EvtDtoKKpiPlot.cc:2227
EvtDtoKKpiPlot::decay
void decay(EvtParticle *p)
Definition EvtDtoKKpiPlot.cc:2193
EvtDtoKKpiPlot::FindYBin
int FindYBin(double mass2)
Definition EvtDtoKKpiPlot.cc:2233
EvtDtoKKpiPlot::~EvtDtoKKpiPlot
virtual ~EvtDtoKKpiPlot()
Definition EvtDtoKKpiPlot.cc:36
EvtDtoKKpiPlot::clone
EvtDecayBase * clone()
Definition EvtDtoKKpiPlot.cc:40
EvtDtoKKpiPlot::getName
void getName(std::string &name)
Definition EvtDtoKKpiPlot.cc:38
EvtPDL::getSpinType
static EvtSpinType::spintype getSpinType(EvtId i)
Definition EvtPDL.hh:66
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:272
EvtParticle::getP4Lab
EvtVector4R getP4Lab()
Definition EvtParticle.cc:694
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtRandom::Flat
static double Flat()
Definition EvtRandom.cc:69