BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDtopipi0pi0Plot.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: EvtDtopipi0pi0Plot.cc
12//
13// Modification history:
14//
15// Liaoyuan Dong August, 2022 Module created
16//
17//------------------------------------------------------------------------
18#include "EvtDtopipi0pi0Plot.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
36EvtDtopipi0pi0Plot::~EvtDtopipi0pi0Plot() {}
37
38void EvtDtopipi0pi0Plot::getName( std::string& model_name ) { model_name = "Dtopipi0pi0Plot"; }
39
40EvtDecayBase* EvtDtopipi0pi0Plot::clone() { return new EvtDtopipi0pi0Plot; }
41
42void EvtDtopipi0pi0Plot::init() {
43
44 checkNArg( 0 );
45
46 bool idN = getDaugs()[0] == EvtPDL::getId( std::string( "pi+" ) ) ||
47 getDaugs()[0] == EvtPDL::getId( std::string( "pi-" ) );
48 bool idKs = getDaugs()[1] == EvtPDL::getId( std::string( "pi0" ) );
49 bool idPi = getDaugs()[2] == EvtPDL::getId( std::string( "pi0" ) );
50 if ( !( idN && idKs && idPi ) )
51 {
52 std::cout << "EvtDtopipi0pi0Plot: the daughter sequence should be pi+ pi0 pi0"
53 << std::endl;
54 abort();
55 }
56 EvtSpinType::spintype parenttype = EvtPDL::getSpinType( getParentId() );
57
58 Xmin = 0.06887;
59 Xmax = 3.0518;
60 Xwid = 0.0596586;
61 Ymin = 0.07037;
62 Ymax = 1.75;
63 Ywid = 0.0335926;
64 avm1 = 0.25;
65 double HisPDFtmp[52][52] = {
66 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
67 0, 0, 0, 0, 0, 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,
69 0, 0, 0, 0, 0, 0, 0,
70 0, 0, 0, 0, 0.126132, 0.0220555, 0.0164693,
71 0.0143236, 0.0615541, 0.0238032, 0.0104177, 0.0514766, 0.00561927, -0.00970829,
72 0.015674, 0.033324, -0.0019576, 0.0179617, 0.00779605, 0.0474916, 0.00439088,
73 0.0191456, 0.00837141, 0.00609675, 0.0840422, 0.0270903, 0.0136279, 0.00854272,
74 0.0278523, 0.00389946, 0.00795048, 0.0450092, 0.0176549, 0.0109168, 0.0398448,
75 0.0144404, 0.0657468, 0 },
76 { 0, 0, 0, 0, 0, 0,
77 0, 0, 0, 0, -0.0122324, -0.00324092,
78 0.0196607, 0.00102101, 0.00378296, 0.00634048, 0.0217871, 0.0081836,
79 0.0140993, 0.00592363, 0.0289073, 0.0116091, 0.015295, 0.0232406,
80 0.016491, 0.0115485, 0.0147406, 0.00684889, 0.00932597, 0.0154916,
81 0.0124573, 0.0166611, 0.0128406, 0.0321747, 0.00527513, 0.011663,
82 0.0212437, 0.000725348, 0.030167, 0.0213951, 0.0147829, 0.0065729,
83 0.0174563, 0.0153622, 0.00786071, 0.0380956, 0.0185521, 0.0102917,
84 0.0206631, 0.0569558, 0, 0 },
85 { 0, 0, 0, 0, 0,
86 0, 0, 0.181193, -0.00514181, -0.00212754,
87 -0.00600245, -0.0047888, 0.00832803, 0.0150015, 0.00962087,
88 0.0126662, 0.0285356, 0.0107688, 0.0164774, 0.0125525,
89 0.0181543, 0.0127544, 0.0180501, 0.0221408, 0.0124838,
90 0.0117366, 0.0103621, 0.00569093, 0.0156799, 0.0104127,
91 0.0134806, 0.0109779, 0.00910387, 0.01754, 0.0138031,
92 0.0108085, 0.0151902, 0.00519298, 0.0126952, 0.0127009,
93 0.0148566, 0.00911647, 0.0103495, 0.0094019, 0.00318659,
94 0.026405, 0.0122733, 0.00718738, 0.0457039, 0,
95 0, 0 },
96 { 0, 0, 0, 0,
97 0, 0, 0.00534517, 0.0171958,
98 -0.00418114, -0.00326283, -0.0103725, -0.00985988,
99 0.0125388, 0.0312826, 0.0177929, 0.0294284,
100 0.0335579, 0.0178685, 0.018411, 0.025134,
101 0.0116659, 0.0108364, 0.0207183, 0.0159755,
102 0.0143117, 0.0141057, 0.00626057, 0.00416359,
103 0.0142275, 0.00944204, 0.00960788, 0.0090878,
104 0.00310312, 0.0050442, 0.0177199, 0.0104969,
105 0.00357364, 0.00960106, -0.00060251, 0.00458286,
106 0.0120471, 0.0169905, 0.00384143, 0,
107 -0.00054573, 0.0156392, 0.0182633, 0,
108 0, 0, 0, 0 },
109 { 0,
110 0,
111 0,
112 0,
113 0,
114 -0.0063761,
115 0.00240996,
116 0.016717,
117 -0.00403196,
118 -0.00315398,
119 -0.0079351,
120 -0.00950213,
121 0.0108979,
122 0.0257102,
123 0.0164158,
124 0.0169603,
125 0.039045,
126 0.0156618,
127 0.0151494,
128 0.0241756,
129 0.0121646,
130 0.014526,
131 0.0206264,
132 0.0169852,
133 0.0141894,
134 0.0114929,
135 0.00528244,
136 0.00348741,
137 0.0119654,
138 0.00944351,
139 0.010706,
140 0.010249,
141 0.00324893,
142 0.00523534,
143 0.0156045,
144 0.012237,
145 0.00412841,
146 0.00945598,
147 -0.000586345,
148 0.00497353,
149 0.00891827,
150 0.016972,
151 0.00367533,
152 0,
153 -0.000588499,
154 0.0145899,
155 0.0285661,
156 0,
157 0,
158 0,
159 0,
160 0 },
161 { 0,
162 0,
163 0,
164 0,
165 0.0214843,
166 0.0152376,
167 0.00361878,
168 0.00659047,
169 0.0015044,
170 0.00383999,
171 0.015594,
172 0.0151121,
173 0.0021722,
174 0.00897155,
175 0.00140851,
176 0.0175562,
177 0.00615938,
178 0.0313451,
179 0.00244723,
180 0.0040116,
181 0.0150819,
182 0.00311473,
183 0.00947754,
184 0,
185 0.0102011,
186 -0.00156026,
187 0.0107358,
188 0.00302003,
189 0.0104483,
190 -0.00228917,
191 -0.000471791,
192 0.0089698,
193 -0.000910742,
194 0.0118354,
195 0.0047521,
196 0.00366432,
197 -0.0004382,
198 0.00986424,
199 0.0137893,
200 0.015136,
201 -0.00208091,
202 0.0104049,
203 0.0257901,
204 0,
205 0.013116,
206 -0.00358994,
207 0.0685484,
208 0,
209 0,
210 0,
211 0,
212 0 },
213 { 0,
214 0,
215 0,
216 0,
217 -0.00226632,
218 0.00771791,
219 -0.00427455,
220 0.0117554,
221 -0.00300342,
222 0.00693718,
223 0.0189155,
224 0.00203017,
225 0.00979611,
226 -0.00473465,
227 0.0173956,
228 0.0131809,
229 0.0128582,
230 0.00145764,
231 0.00549767,
232 0.00423976,
233 0.00889498,
234 0.0105157,
235 0.009075,
236 0.00134348,
237 -0.000548315,
238 0.00165993,
239 0,
240 0.0108027,
241 0.0228641,
242 0.00539549,
243 -0.00174531,
244 -0.000652211,
245 0.00788748,
246 0.00229044,
247 0.01005,
248 0.0107275,
249 0.0084291,
250 0.00547439,
251 0.00902564,
252 -0.000925959,
253 0.0177774,
254 -0.00256084,
255 0.0130949,
256 0,
257 -0.00222004,
258 0,
259 0,
260 0,
261 0,
262 0,
263 0,
264 0 },
265 { 0,
266 0,
267 0,
268 0.00550887,
269 0.00193808,
270 0.00608034,
271 0.0181517,
272 -0.00603111,
273 -0.00472705,
274 0.0131492,
275 0.00704957,
276 0.0114897,
277 0.00189524,
278 0.00536718,
279 0.0189741,
280 0.0298482,
281 0.00440453,
282 0.00249411,
283 0.0229042,
284 0.0121403,
285 0.0107645,
286 0.00982848,
287 0.00466503,
288 0.0182887,
289 -0.00123724,
290 -0.000535069,
291 0.0112898,
292 0.0043009,
293 -0.000457878,
294 0.00370862,
295 0.0168827,
296 0.0104902,
297 0.00383087,
298 0.00505172,
299 0.0046832,
300 0.016317,
301 0.00371058,
302 0.00429451,
303 0,
304 0.00449643,
305 0.00984612,
306 0,
307 0,
308 -0.00110221,
309 0,
310 0,
311 0,
312 0,
313 0,
314 0,
315 0,
316 0 },
317 { 0, 0,
318 0, 0.00559425,
319 0.00784236, -0.00515694,
320 -0.00341012, 0.0056297,
321 0.00161726, 0.00817172,
322 0.0135242, 0.00792446,
323 0.0130368, 0.00220852,
324 0.00213705, 0.0213428,
325 0.00845069, -0.00640989,
326 0.00147561, -0.00334377,
327 0.00396805, -0.00163827,
328 0.015181, -0.00160953,
329 0.00974975, -0.000859794,
330 0.0071086, 0.00382806,
331 -0.00327447, -0.00100732,
332 0.0095671, 0.00985691,
333 0.00516776, -0.000591646,
334 0.0052562, -0.0020691,
335 0.0125954, 0.00530462,
336 0.00524183, -0.00200351,
337 0.00427751, 0.00629911,
338 -0.00116534, 0,
339 0, 0,
340 0, 0,
341 0, 0,
342 0, 0 },
343 { 0,
344 0,
345 0.00924787,
346 0.0220888,
347 0.00547255,
348 -0.00053368,
349 -0.00621351,
350 0.00340835,
351 0.000869906,
352 -0.00453506,
353 -0.00230228,
354 0.00977444,
355 0.00150216,
356 0.0126651,
357 0.00664757,
358 0.00845056,
359 0.0145116,
360 0.0219415,
361 0.00305964,
362 0.0159933,
363 0.00209908,
364 -0.00145665,
365 -0.00114919,
366 0.00413543,
367 0.0037937,
368 0.0197595,
369 -0.000571623,
370 0.00344037,
371 0,
372 0.0110179,
373 0.00531177,
374 0.00501057,
375 -0.001035,
376 0.00452493,
377 0.0103877,
378 0.00574533,
379 0,
380 0.00384947,
381 -0.00148758,
382 0.00478215,
383 0,
384 0,
385 0,
386 0,
387 0,
388 0,
389 0,
390 0,
391 0,
392 0,
393 0,
394 0 },
395 { 0, 0, 0.00740427, 0.00720438,
396 0.00554212, 0.00468727, 0.00749359, 0.0020956,
397 0.00530923, -0.00152604, -0.00157443, 0.0066342,
398 0.00491629, 0.0124885, 0.0119529, 0.016968,
399 0.00564861, 0.00142758, 0.0032944, 0.0125377,
400 0.0107067, -0.0025167, 0.00932271, 0.0151918,
401 0, -0.00107743, -0.00102267, 0.0103105,
402 0.00630324, -0.00215053, 0.00358067, 0.00498285,
403 0.0187526, 0, 0, 0.00485819,
404 0, 0.0109325, 0.00526438, 0.0087764,
405 0, 0, 0, 0,
406 0, 0, 0, 0,
407 0, 0, 0, 0 },
408 { 0,
409 0,
410 -0.00188838,
411 -0.00136398,
412 0.00841484,
413 0.00796693,
414 0.023253,
415 -0.00402091,
416 0.012051,
417 0.0202265,
418 0.0140749,
419 0.0122277,
420 0.01082,
421 0.00468877,
422 0.00300136,
423 0.0257563,
424 0.0326677,
425 0.00203623,
426 0.00314237,
427 0.0101705,
428 0.00789608,
429 -0.000500071,
430 0.00266911,
431 0.00330721,
432 0.00524247,
433 0.0163683,
434 0,
435 0.00357269,
436 -0.000856078,
437 0.00566965,
438 0.00332607,
439 0.0106486,
440 0.00473193,
441 0,
442 0.00403549,
443 0,
444 -0.00104048,
445 -0.0012422,
446 0,
447 0,
448 -0.00278746,
449 0,
450 0,
451 0,
452 0,
453 0,
454 0,
455 0,
456 0,
457 0,
458 0,
459 0 },
460 { 0, -0.0327654, 0.00279954, -0.00156276,
461 -0.00392664, -0.00105142, 0.00949674, 0.00992266,
462 0.00251736, 0.0152795, 0.00306799, -0.00238284,
463 0.0067777, -0.00157091, 0.00526856, 0.0264038,
464 0.00147724, 0.00888431, 0.026411, 0.0198684,
465 -0.00143129, 0.00316219, 0.00514419, 0,
466 0.00560927, 0.00420765, -0.00145472, 0,
467 0.00956539, 0.00241543, 0.00470629, 0.0048465,
468 0.0144799, 0, 0.00570209, 0.00418751,
469 -0.00047658, -0.00110488, 0, 0,
470 0, 0, 0, 0,
471 0, 0, 0, 0,
472 0, 0, 0, 0 },
473 { 0,
474 -0.00753227,
475 0.00397765,
476 0.0123975,
477 0.000839958,
478 -0.0042778,
479 -0.00294332,
480 -0.00318388,
481 -0.00243692,
482 0.00213064,
483 -0.000518138,
484 0.00777234,
485 0.00394029,
486 0,
487 0.00305236,
488 0.00821964,
489 -0.000475594,
490 0.0129634,
491 0.00775832,
492 0.0247704,
493 0.00518487,
494 0.00437631,
495 0.00999549,
496 0.0040436,
497 0.00339696,
498 0.00587269,
499 0.0123544,
500 0,
501 0.0113253,
502 0.00504108,
503 0.00467702,
504 0.00592753,
505 0.00563897,
506 0.0149573,
507 0.00464732,
508 0.00463117,
509 0,
510 0.00761649,
511 0.0151458,
512 0,
513 0,
514 0,
515 0,
516 0,
517 0,
518 0,
519 0,
520 0,
521 0,
522 0,
523 0,
524 0 },
525 { 0,
526 -0.00191676,
527 -0.00423107,
528 0.00475804,
529 -0.00160728,
530 0.0088541,
531 0.00857197,
532 -0.000837,
533 -0.00239001,
534 -0.00049013,
535 -0.000551548,
536 -0.000954402,
537 0.00374396,
538 -0.00129406,
539 0.0189372,
540 -0.00141928,
541 -0.000367202,
542 0.00547228,
543 0.0088018,
544 -0.00043201,
545 0,
546 0.0142185,
547 0.0103194,
548 -0.000443245,
549 0.00500681,
550 0.00566802,
551 0.0117745,
552 0.0128396,
553 -0.00146617,
554 0.00538462,
555 0.0127682,
556 0.00472659,
557 0.00756395,
558 0.0138696,
559 -0.000561553,
560 -0.00167291,
561 0,
562 0,
563 0,
564 0,
565 0,
566 0,
567 0,
568 0,
569 0,
570 0,
571 0,
572 0,
573 0,
574 0,
575 0,
576 0 },
577 { 0,
578 -0.00437213,
579 -0.00102453,
580 0.00473396,
581 -0.0024484,
582 -0.00156082,
583 -0.000590814,
584 -0.00124694,
585 -0.00349041,
586 0.00456603,
587 -0.000421245,
588 0,
589 0.00621102,
590 0.00749433,
591 -0.000960832,
592 0.00314721,
593 0.00359019,
594 0.00284196,
595 0.0137168,
596 0.00562268,
597 -0.00160554,
598 0,
599 0,
600 -0.000470884,
601 0.00336691,
602 0.00459624,
603 -0.00111263,
604 0.0050625,
605 0.0039622,
606 0.00446919,
607 0,
608 0,
609 -0.000451104,
610 0.00630724,
611 0.0259699,
612 -0.0016493,
613 -0.00143629,
614 0,
615 0,
616 0,
617 0,
618 0,
619 0,
620 0,
621 0,
622 0,
623 0,
624 0,
625 0,
626 0,
627 0,
628 0 },
629 { 0,
630 0.0116541,
631 -0.000491055,
632 0,
633 0.0108019,
634 0.00891202,
635 0.0205461,
636 -0.000482761,
637 -0.00115522,
638 -0.00172886,
639 0.00237479,
640 -0.000464065,
641 0.000745853,
642 0.0113383,
643 -0.00227308,
644 0,
645 -0.00154217,
646 -0.000569069,
647 0.00912503,
648 -0.000523669,
649 0,
650 0.0054163,
651 -0.00129641,
652 -0.00213726,
653 0.0056876,
654 0.00560903,
655 0,
656 -0.000929603,
657 0,
658 0,
659 0.00435316,
660 0,
661 -0.000487678,
662 0.0115152,
663 0.00722282,
664 0,
665 0,
666 0,
667 0,
668 0,
669 0,
670 0,
671 0,
672 0,
673 0,
674 0,
675 0,
676 0,
677 0,
678 0,
679 0,
680 0 },
681 { 0,
682 0.0101944,
683 -0.000458053,
684 -0.00136645,
685 0.00293929,
686 -0.00104558,
687 -0.00138759,
688 -0.000995768,
689 0.0041222,
690 -0.000817748,
691 0.00406901,
692 0,
693 0.00431831,
694 -0.000359942,
695 -0.000918891,
696 0.0039494,
697 0.00456461,
698 -0.000445733,
699 -0.000497338,
700 0,
701 0.00431183,
702 0.00336932,
703 0.00508663,
704 0,
705 -0.000925412,
706 0.00551509,
707 0.00448889,
708 -0.000465728,
709 0.00337366,
710 0.00468909,
711 0,
712 -0.00137258,
713 0.00603869,
714 0.00806585,
715 0,
716 -0.00580589,
717 0,
718 0,
719 0,
720 0,
721 0,
722 0,
723 0,
724 0,
725 0,
726 0,
727 0,
728 0,
729 0,
730 0,
731 0,
732 0 },
733 { 0,
734 -0.00237226,
735 0.010259,
736 -0.000895674,
737 0.00411476,
738 0.00512026,
739 -0.00123773,
740 0.00821976,
741 0.00269015,
742 0.00838505,
743 -0.000443918,
744 0.00770911,
745 -0.000451037,
746 -0.000919729,
747 0.00376047,
748 0,
749 -0.00128212,
750 0.00577606,
751 0.00402539,
752 -0.00127397,
753 0.00444049,
754 0,
755 -0.000432523,
756 -0.00189025,
757 -0.000976431,
758 0,
759 0.00498299,
760 0,
761 0.00437452,
762 0,
763 0,
764 0,
765 0,
766 0,
767 0,
768 0,
769 0,
770 0,
771 0,
772 0,
773 0,
774 0,
775 0,
776 0,
777 0,
778 0,
779 0,
780 0,
781 0,
782 0,
783 0,
784 0 },
785 { 0,
786 0,
787 0.00537186,
788 0.00395098,
789 0.0117299,
790 0.00799079,
791 -0.00188147,
792 0.00901231,
793 -0.00115061,
794 -0.0022369,
795 0.00940017,
796 -0.000452578,
797 0.00529762,
798 -0.000898868,
799 -0.00113458,
800 0.00950784,
801 0,
802 -0.00189087,
803 -0.000408221,
804 0,
805 -0.00091378,
806 0.00349189,
807 0,
808 0.00454828,
809 0,
810 -0.000567589,
811 0,
812 0,
813 -0.000869253,
814 0,
815 0.00640011,
816 0,
817 0,
818 0,
819 0,
820 0,
821 0,
822 0,
823 0,
824 0,
825 0,
826 0,
827 0,
828 0,
829 0,
830 0,
831 0,
832 0,
833 0,
834 0,
835 0,
836 0 },
837 { 0,
838 0,
839 0.00875028,
840 0.00507708,
841 0,
842 -0.00136453,
843 0.00372118,
844 0.00202545,
845 0,
846 0,
847 -0.00103392,
848 -0.00142976,
849 0,
850 0.00364562,
851 0.0036454,
852 0.00372184,
853 -0.000416832,
854 -0.00122799,
855 -0.00117455,
856 -0.000848612,
857 0,
858 0.00619048,
859 0,
860 0,
861 0,
862 0.00300899,
863 0,
864 0.00439737,
865 0,
866 -0.000452727,
867 0,
868 0,
869 0.0223262,
870 0,
871 0,
872 0,
873 0,
874 0,
875 0,
876 0,
877 0,
878 0,
879 0,
880 0,
881 0,
882 0,
883 0,
884 0,
885 0,
886 0,
887 0,
888 0 },
889 { 0,
890 0.00813083,
891 0.00984675,
892 0.00418989,
893 0.00419656,
894 -0.000805791,
895 -0.00125829,
896 0.00409336,
897 0.00459685,
898 0.00424519,
899 -0.000454501,
900 0.00829167,
901 0,
902 0.00525152,
903 -0.000449865,
904 0.00472471,
905 -0.000496281,
906 -0.00134592,
907 0,
908 0.013176,
909 0.00972065,
910 0,
911 0.00337632,
912 0.00312468,
913 0,
914 0,
915 0.00465113,
916 0.00523123,
917 0.0096822,
918 0,
919 0,
920 0,
921 0,
922 0,
923 0,
924 0,
925 0,
926 0,
927 0,
928 0,
929 0,
930 0,
931 0,
932 0,
933 0,
934 0,
935 0,
936 0,
937 0,
938 0,
939 0,
940 0 },
941 { 0,
942 0.0078588,
943 0.00438392,
944 0.00364588,
945 0.00309135,
946 -0.00206823,
947 -0.000886876,
948 0.00303936,
949 0.00295383,
950 -0.00250436,
951 -0.00079063,
952 0.00529198,
953 0,
954 0,
955 0,
956 -0.00123548,
957 0.00529033,
958 -0.00133633,
959 0.0055804,
960 -0.0005559,
961 0,
962 0,
963 -0.000420998,
964 0.0047264,
965 0.00372627,
966 0.0050624,
967 0,
968 0.00417586,
969 0.00579403,
970 0.00922432,
971 0,
972 0,
973 0,
974 0,
975 0,
976 0,
977 0,
978 0,
979 0,
980 0,
981 0,
982 0,
983 0,
984 0,
985 0,
986 0,
987 0,
988 0,
989 0,
990 0,
991 0,
992 0 },
993 { 0,
994 0.0134527,
995 -0.000487816,
996 -0.000407827,
997 -0.00164327,
998 0.00714982,
999 0.00834671,
1000 0.0043254,
1001 -0.00150873,
1002 0.00909819,
1003 -0.00107058,
1004 -0.000881873,
1005 -0.000906969,
1006 -0.000801329,
1007 -0.00141174,
1008 0.00361808,
1009 0,
1010 0.00884207,
1011 -0.000526459,
1012 0.00422741,
1013 0.00892857,
1014 0.00334632,
1015 -0.00141564,
1016 0.00534427,
1017 -0.0007912,
1018 0.00382144,
1019 0,
1020 0.00361667,
1021 0.00735391,
1022 0,
1023 0,
1024 0,
1025 0,
1026 0,
1027 0,
1028 0,
1029 0,
1030 0,
1031 0,
1032 0,
1033 0,
1034 0,
1035 0,
1036 0,
1037 0,
1038 0,
1039 0,
1040 0,
1041 0,
1042 0,
1043 0,
1044 0 },
1045 { 0,
1046 -0.00112453,
1047 -0.000921719,
1048 0.00348905,
1049 0.00986824,
1050 0.00323772,
1051 0.00835248,
1052 0.00388933,
1053 0.00499419,
1054 -0.00153203,
1055 0.00456795,
1056 0,
1057 0,
1058 0.00465839,
1059 0,
1060 0.00372416,
1061 0.00465732,
1062 0.00646166,
1063 0,
1064 0.00444293,
1065 0.00290251,
1066 0,
1067 0,
1068 0.00934281,
1069 0.0148351,
1070 -0.000809891,
1071 0.0100435,
1072 0,
1073 -0.00301733,
1074 -0.00375312,
1075 0,
1076 0,
1077 0,
1078 0,
1079 0,
1080 0,
1081 0,
1082 0,
1083 0,
1084 0,
1085 0,
1086 0,
1087 0,
1088 0,
1089 0,
1090 0,
1091 0,
1092 0,
1093 0,
1094 0,
1095 0,
1096 0 },
1097 { 0,
1098 0.0106236,
1099 0.00129889,
1100 0.00448691,
1101 0.00445295,
1102 0.00508553,
1103 -0.00138142,
1104 0,
1105 -0.000445329,
1106 0.00403718,
1107 -0.000866474,
1108 0.00353158,
1109 0,
1110 -0.0014014,
1111 0.00440745,
1112 -0.000460064,
1113 0.00324248,
1114 0.0143277,
1115 0,
1116 -0.00102167,
1117 0.00481517,
1118 0.00437973,
1119 0.0138674,
1120 0.0263436,
1121 0.00353529,
1122 0.00439296,
1123 0.008448,
1124 0,
1125 0,
1126 0,
1127 0,
1128 0,
1129 0,
1130 0,
1131 0,
1132 0,
1133 0,
1134 0,
1135 0,
1136 0,
1137 0,
1138 0,
1139 0,
1140 0,
1141 0,
1142 0,
1143 0,
1144 0,
1145 0,
1146 0,
1147 0,
1148 0 },
1149 { 0,
1150 -0.000615042,
1151 0.00934905,
1152 0.00757931,
1153 -0.00193735,
1154 0,
1155 0,
1156 0,
1157 -0.00162904,
1158 -0.00201206,
1159 0.00316532,
1160 -0.000758207,
1161 -0.000460719,
1162 0.00393524,
1163 0.00274551,
1164 -0.00122072,
1165 -0.00328414,
1166 0.00975402,
1167 0.00479703,
1168 0.00520447,
1169 -0.0012812,
1170 0.003169,
1171 0.00331205,
1172 0.0146124,
1173 0.0111438,
1174 0.00451904,
1175 0,
1176 0.075641,
1177 0,
1178 0,
1179 0,
1180 0,
1181 0,
1182 0,
1183 0,
1184 0,
1185 0,
1186 0,
1187 0,
1188 0,
1189 0,
1190 0,
1191 0,
1192 0,
1193 0,
1194 0,
1195 0,
1196 0,
1197 0,
1198 0,
1199 0,
1200 0 },
1201 { 0,
1202 -0.000516182,
1203 -0.000893867,
1204 0.0121081,
1205 -0.0018224,
1206 -0.000488139,
1207 0.00268232,
1208 0,
1209 0.00402894,
1210 -0.000402761,
1211 0,
1212 0,
1213 0,
1214 -0.000405727,
1215 0.00430231,
1216 0.00434484,
1217 0.00410179,
1218 -0.000975305,
1219 0.0160911,
1220 0.00498146,
1221 0.0173023,
1222 0.00533898,
1223 0.0278033,
1224 -0.000951397,
1225 0.0128432,
1226 0,
1227 0,
1228 0,
1229 0,
1230 0,
1231 0,
1232 0,
1233 0,
1234 0,
1235 0,
1236 0,
1237 0,
1238 0,
1239 0,
1240 0,
1241 0,
1242 0,
1243 0,
1244 0,
1245 0,
1246 0,
1247 0,
1248 0,
1249 0,
1250 0,
1251 0,
1252 0 },
1253 { 0,
1254 -0.00140846,
1255 0.00422503,
1256 -0.00163652,
1257 0.00443282,
1258 0,
1259 -0.000351905,
1260 0.00466989,
1261 0.00342072,
1262 0.00303109,
1263 0,
1264 -0.000763762,
1265 0,
1266 -0.000770398,
1267 0.0202359,
1268 0.00460753,
1269 0.014308,
1270 0.00482234,
1271 -0.00179807,
1272 0.0263769,
1273 0.00670195,
1274 0.00955175,
1275 0.0202865,
1276 0.0155921,
1277 0.0189799,
1278 0.019989,
1279 0,
1280 0,
1281 0,
1282 0,
1283 0,
1284 0,
1285 0,
1286 0,
1287 0,
1288 0,
1289 0,
1290 0,
1291 0,
1292 0,
1293 0,
1294 0,
1295 0,
1296 0,
1297 0,
1298 0,
1299 0,
1300 0,
1301 0,
1302 0,
1303 0,
1304 0 },
1305 { 0, 0,
1306 0, -0.00132663,
1307 0, -0.000857411,
1308 0, -0.000393507,
1309 0, -0.000818886,
1310 0.00722913, -0.000422264,
1311 0, -0.00124678,
1312 0.00794475, 0.0121623,
1313 0.00420758, 0.00262983,
1314 0.00718098, -0.000505903,
1315 0.00895671, 0.0246874,
1316 0.0180739, 0.00844876,
1317 0.0230061, 0,
1318 0, 0,
1319 0, 0,
1320 0, 0,
1321 0, 0,
1322 0, 0,
1323 0, 0,
1324 0, 0,
1325 0, 0,
1326 0, 0,
1327 0, 0,
1328 0, 0,
1329 0, 0,
1330 0, 0 },
1331 { 0,
1332 0,
1333 -0.00043891,
1334 0.0057952,
1335 -0.00153878,
1336 -0.00204123,
1337 -0.00274566,
1338 -0.00166583,
1339 0.00470958,
1340 -0.000736231,
1341 -0.00111919,
1342 -0.00081959,
1343 -0.00271149,
1344 0.00407672,
1345 -0.000947025,
1346 0.00436545,
1347 0.0158883,
1348 0.0121416,
1349 0.0202014,
1350 0.0141576,
1351 0.00915484,
1352 0.00371495,
1353 -0.000838394,
1354 0.0185427,
1355 -0.0330275,
1356 0,
1357 0,
1358 0,
1359 0,
1360 0,
1361 0,
1362 0,
1363 0,
1364 0,
1365 0,
1366 0,
1367 0,
1368 0,
1369 0,
1370 0,
1371 0,
1372 0,
1373 0,
1374 0,
1375 0,
1376 0,
1377 0,
1378 0,
1379 0,
1380 0,
1381 0,
1382 0 },
1383 { 0,
1384 -0.000569017,
1385 -0.00187911,
1386 -0.0013152,
1387 -0.00205846,
1388 -0.00174749,
1389 -0.000859905,
1390 -0.00139486,
1391 0.00919869,
1392 -0.00162995,
1393 -0.0016275,
1394 -0.00136203,
1395 0.00943644,
1396 0.00406865,
1397 0.00304173,
1398 0.00802036,
1399 0.0223577,
1400 0.0243796,
1401 0.0195699,
1402 0.00451811,
1403 0.00637247,
1404 0.00861249,
1405 -0.00157236,
1406 0,
1407 0,
1408 0,
1409 0,
1410 0,
1411 0,
1412 0,
1413 0,
1414 0,
1415 0,
1416 0,
1417 0,
1418 0,
1419 0,
1420 0,
1421 0,
1422 0,
1423 0,
1424 0,
1425 0,
1426 0,
1427 0,
1428 0,
1429 0,
1430 0,
1431 0,
1432 0,
1433 0,
1434 0 },
1435 { 0, 0.00676314, 0, 0,
1436 -0.00151701, -0.00167819, -0.00164177, 0,
1437 0, -0.00126364, 0.00209941, 0.00780059,
1438 0.00523671, 0.00816129, 0.00920645, 0.00914593,
1439 0.0480522, 0.0182587, 0.0174034, 0.0129708,
1440 0, 0.0709505, 0, 0,
1441 0, 0, 0, 0,
1442 0, 0, 0, 0,
1443 0, 0, 0, 0,
1444 0, 0, 0, 0,
1445 0, 0, 0, 0,
1446 0, 0, 0, 0,
1447 0, 0, 0, 0 },
1448 { 0,
1449 -0.00132961,
1450 -0.00105272,
1451 -0.000493242,
1452 -0.000998982,
1453 -0.000963302,
1454 -0.000403533,
1455 0.00321011,
1456 0.00259871,
1457 0.00165038,
1458 0.00945111,
1459 -0.000398643,
1460 0.0215238,
1461 0.00245585,
1462 0.0160536,
1463 0.0330436,
1464 0.038898,
1465 0.0155603,
1466 0.0113935,
1467 0.00682208,
1468 0.0266945,
1469 0.0601852,
1470 0,
1471 0,
1472 0,
1473 0,
1474 0,
1475 0,
1476 0,
1477 0,
1478 0,
1479 0,
1480 0,
1481 0,
1482 0,
1483 0,
1484 0,
1485 0,
1486 0,
1487 0,
1488 0,
1489 0,
1490 0,
1491 0,
1492 0,
1493 0,
1494 0,
1495 0,
1496 0,
1497 0,
1498 0,
1499 0 },
1500 { 0, 0, -0.00219801, -0.00100347, -0.00112403,
1501 0.00820631, 0, 0.000776552, 0, -0.00309955,
1502 0.00344142, 0.0165217, 0.00824275, 0.00318922, 0.0167143,
1503 0.0183902, 0.0209505, 0.0148305, 0.026167, 0.0465092,
1504 0.0393338, 0, 0, 0, 0,
1505 0, 0, 0, 0, 0,
1506 0, 0, 0, 0, 0,
1507 0, 0, 0, 0, 0,
1508 0, 0, 0, 0, 0,
1509 0, 0, 0, 0, 0,
1510 0, 0 },
1511 { 0, -0.00192963,
1512 0.00404023, 0.00429177,
1513 0.00361765, -0.00150455,
1514 0.00368992, 0.00358212,
1515 0.00743376, -0.000420885,
1516 0.0140216, 0.00868896,
1517 0.00296287, 0.0174345,
1518 0.0200781, 0.0511327,
1519 0.03534, 0.0113812,
1520 0.0331168, -0.00423133,
1521 0, 0,
1522 0, 0,
1523 0, 0,
1524 0, 0,
1525 0, 0,
1526 0, 0,
1527 0, 0,
1528 0, 0,
1529 0, 0,
1530 0, 0,
1531 0, 0,
1532 0, 0,
1533 0, 0,
1534 0, 0,
1535 0, 0,
1536 0, 0 },
1537 { 0, -0.0017667, 0.00266838, 0.00821613, 0.00252852, -0.000446041,
1538 0.00859937, 0.00201737, 0.00433766, 0.0037743, 0.0244391, 0.0120824,
1539 0.0348515, 0.0460445, 0.0198451, 0.0172826, 0.0730303, 0.125063,
1540 0.0233494, 0, 0, 0, 0, 0,
1541 0, 0, 0, 0, 0, 0,
1542 0, 0, 0, 0, 0, 0,
1543 0, 0, 0, 0, 0, 0,
1544 0, 0, 0, 0, 0, 0,
1545 0, 0, 0, 0 },
1546 { 0, -0.00103876, 0.00270405, 0.00510568, 0.0100982, 0.0168746, 0.0105157,
1547 0.00277159, 0.00682157, 0.0246823, 0.00884256, 0.00534351, 0.0465627, 0.0387692,
1548 0.071591, 0.0284448, 0.0953723, 0.0658472, 0, 0, 0,
1549 0, 0, 0, 0, 0, 0, 0,
1550 0, 0, 0, 0, 0, 0, 0,
1551 0, 0, 0, 0, 0, 0, 0,
1552 0, 0, 0, 0, 0, 0, 0,
1553 0, 0, 0 },
1554 { 0, 0, 0.0112434, -0.000430343,
1555 0.00425268, -0.00161899, 0.0131034, 0.0092693,
1556 0.00747454, 0.00379318, 0.0193525, 0.0450718,
1557 0.0749774, 0.045105, 0.0726199, 0.0665514,
1558 0.133822, 0, 0, 0,
1559 0, 0, 0, 0,
1560 0, 0, 0, 0,
1561 0, 0, 0, 0,
1562 0, 0, 0, 0,
1563 0, 0, 0, 0,
1564 0, 0, 0, 0,
1565 0, 0, 0, 0,
1566 0, 0, 0, 0 },
1567 { 0, -0.0074107, 0.00855827, 0.00373376, 0.00884158, 0.0110644, 0.00795127,
1568 0.0206009, 0.0167017, 0.0175845, 0.0577236, 0.0370533, 0.0527817, 0.0650344,
1569 0.10323, 0.224505, 0.0510321, 0, 0, 0, 0,
1570 0, 0, 0, 0, 0, 0, 0,
1571 0, 0, 0, 0, 0, 0, 0,
1572 0, 0, 0, 0, 0, 0, 0,
1573 0, 0, 0, 0, 0, 0, 0,
1574 0, 0, 0 },
1575 { 0, 0, 0.0179825, 0.0167169, 0.00888093, -0.000476096,
1576 0.0128328, 0.0112438, 0.0405886, 0.0298025, 0.0406211, 0.0588893,
1577 0.0363251, 0.0611834, 0.0628016, 0.192035, 0, 0,
1578 0, 0, 0, 0, 0, 0,
1579 0, 0, 0, 0, 0, 0,
1580 0, 0, 0, 0, 0, 0,
1581 0, 0, 0, 0, 0, 0,
1582 0, 0, 0, 0, 0, 0,
1583 0, 0, 0, 0 },
1584 { 0, 0, 0.0114356, 0.00720186, 0.0170689, 0.017875, 0.0293027, 0.023858,
1585 0.0125951, 0.0270779, 0.0515326, 0.0715181, 0.0579174, 0.10595, 0.0701754, 0,
1586 0, 0, 0, 0, 0, 0, 0, 0,
1587 0, 0, 0, 0, 0, 0, 0, 0,
1588 0, 0, 0, 0, 0, 0, 0, 0,
1589 0, 0, 0, 0, 0, 0, 0, 0,
1590 0, 0, 0, 0 },
1591 { 0, 0, 0.00974123, 0.0251431, 0.0149583, 0.0333583, 0.0190875,
1592 0.0191277, 0.0457674, 0.0370717, 0.0917013, 0.0913667, 0.105546, 0.0307692,
1593 0, 0, 0, 0, 0, 0, 0,
1594 0, 0, 0, 0, 0, 0, 0,
1595 0, 0, 0, 0, 0, 0, 0,
1596 0, 0, 0, 0, 0, 0, 0,
1597 0, 0, 0, 0, 0, 0, 0,
1598 0, 0, 0 },
1599 { 0, 0, -0.00060854, 0.0338357, 0.0416551, 0.0254439,
1600 0.0576751, 0.0666466, 0.0329934, 0.0596387, 0.0719475, 0.0694444,
1601 0.0710059, 0, 0, 0, 0, 0,
1602 0, 0, 0, 0, 0, 0,
1603 0, 0, 0, 0, 0, 0,
1604 0, 0, 0, 0, 0, 0,
1605 0, 0, 0, 0, 0, 0,
1606 0, 0, 0, 0, 0, 0,
1607 0, 0, 0, 0 },
1608 { 0, 0, 0.0141723, 0.0627466, 0.0553281, 0.0745046, 0.0518118, 0.0624949,
1609 0.0354849, 0.0395422, 0.0723982, 0.0326973, 0, 0, 0, 0,
1610 0, 0, 0, 0, 0, 0, 0, 0,
1611 0, 0, 0, 0, 0, 0, 0, 0,
1612 0, 0, 0, 0, 0, 0, 0, 0,
1613 0, 0, 0, 0, 0, 0, 0, 0,
1614 0, 0, 0, 0 },
1615 { 0, 0, 0.0441052, 0.0414434, 0.0610666, 0.0120298, 0.0320438, 0.0192959,
1616 0.0258374, 0.0776923, 0.0257947, 0, 0, 0, 0, 0,
1617 0, 0, 0, 0, 0, 0, 0, 0,
1618 0, 0, 0, 0, 0, 0, 0, 0,
1619 0, 0, 0, 0, 0, 0, 0, 0,
1620 0, 0, 0, 0, 0, 0, 0, 0,
1621 0, 0, 0, 0 },
1622 { 0, 0, 0, 0.0417211, 0.00655413, 0.0419119, 0.0660837, 0.0144257, 0.0221876,
1623 0.0484496, 0, 0, 0, 0, 0, 0, 0, 0,
1624 0, 0, 0, 0, 0, 0, 0, 0, 0,
1625 0, 0, 0, 0, 0, 0, 0, 0, 0,
1626 0, 0, 0, 0, 0, 0, 0, 0, 0,
1627 0, 0, 0, 0, 0, 0, 0 },
1628 { 0, 0, 0, 0.0515969, 0.042548, 0.046319, 0.0160193, 0.0190285, -0.00564746, 0, 0,
1629 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1630 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1631 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1632 0, 0, 0, 0, 0, 0, 0, 0 },
1633 { 0, 0, 0, 0, 0.046252, 0.0763126, 0, 0, 0.25, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1634 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1635 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1636 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1637 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1638 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1639 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 } };
1640
1641 for ( int i = 0; i < 52; i++ )
1642 {
1643 for ( int j = 0; j < 52; j++ ) { HisPDF[i][j] = HisPDFtmp[i][j]; }
1644 }
1645}
1646
1647void EvtDtopipi0pi0Plot::initProbMax() { noProbMax(); }
1648
1649void EvtDtopipi0pi0Plot::decay( EvtParticle* p ) {
1650
1651loop:
1652 p->initializePhaseSpace( getNDaug(), getDaugs() );
1653
1654 EvtParticle *id1, *id2, *id3;
1655 EvtVector4R pd1, pd2, pd3;
1656 double xmass13, xmass12, xmass23;
1657
1658 id1 = p->getDaug( 0 );
1659 id2 = p->getDaug( 1 );
1660 id3 = p->getDaug( 2 );
1661
1662 pd1 = id1->getP4Lab();
1663 pd2 = id2->getP4Lab();
1664 pd3 = id3->getP4Lab();
1665
1666 // xmass12=(pd1+pd2).mass()*(pd1+pd2).mass(); // M_ksopi
1667 xmass13 = ( pd1 + pd3 ).mass() * ( pd1 + pd3 ).mass(); // M_ksow
1668 xmass23 = ( pd2 + pd3 ).mass() * ( pd2 + pd3 ).mass(); // M_piw
1669
1670 int xbin = FindXBin( xmass23 );
1671 int ybin = FindYBin( xmass13 );
1672 double xratio12 = HisPDF[xbin][ybin] / avm1;
1673
1674 if ( xratio12 <= 0 ) goto loop;
1675
1676 double rd12 = EvtRandom::Flat( 0.0, 1.0 );
1677 if ( rd12 > xratio12 ) goto loop;
1678
1679 return;
1680}
1681
1682int EvtDtopipi0pi0Plot::FindXBin( double mass2 ) {
1683 if ( mass2 < Xmin ) { return 0; }
1684 else if ( mass2 >= Xmax ) { return 51; }
1685 else { return int( ( mass2 - Xmin ) / Xwid ) + 1; }
1686}
1687
1688int EvtDtopipi0pi0Plot::FindYBin( double mass2 ) {
1689 if ( mass2 < Ymin ) { return 0; }
1690 else if ( mass2 >= Ymax ) { return 51; }
1691 else { return int( ( mass2 - Ymin ) / Ywid ) + 1; }
1692}
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
EvtDtopipi0pi0Plot::FindXBin
int FindXBin(double mass2)
Definition EvtDtopipi0pi0Plot.cc:1682
EvtDtopipi0pi0Plot::decay
void decay(EvtParticle *p)
Definition EvtDtopipi0pi0Plot.cc:1649
EvtDtopipi0pi0Plot::FindYBin
int FindYBin(double mass2)
Definition EvtDtopipi0pi0Plot.cc:1688
EvtDtopipi0pi0Plot::clone
EvtDecayBase * clone()
Definition EvtDtopipi0pi0Plot.cc:40
EvtDtopipi0pi0Plot::getName
void getName(std::string &name)
Definition EvtDtopipi0pi0Plot.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