BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDtopipipiPlot.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: EvtDtopipipiPlot.cc
12//
13// Modification history:
14//
15// Liaoyuan Dong August, 2022 Module created
16//
17//------------------------------------------------------------------------
18#include "EvtDtopipipiPlot.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
37
38void EvtDtopipipiPlot::getName( std::string& model_name ) { model_name = "DtopipipiPlot"; }
39
41
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( "pi+" ) ) ||
49 getDaugs()[1] == EvtPDL::getId( std::string( "pi-" ) );
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 << "EvtDtopipipiPlot: the daughter sequence should be pi+ pi+ pi-" << std::endl;
55 abort();
56 }
58
59 Xmin = 0.072919;
60 Xmax = 1.74;
61 Xwid = 0.0333416;
62 Ymin = 0.2;
63 Ymax = 3.1229;
64 Ywid = 0.058458;
65 avm1 = 0.2;
66 double HisPDFtmp[52][52] = {
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, 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,
70 0, 0, 0, 0, 0,
71 0, 0, 0, -0.00150066, 0.0155324,
72 0.00978731, 0, 0, 0, -0.000862562,
73 -0.00027163, 0.00905843, 0.00159547, 0.00674777, 0.00686001,
74 0.00893704, 0.0263911, 0.0177556, 0.0262276, 0.0163004,
75 0.00810661, 0.00266422, 0.0165903, 0.0130485, 0.0090461,
76 0.0138679, 0.00678281, 0.0142944, 0.0124154, 0.00501558,
77 0, 0, 0, 0, 0,
78 0, 0, 0, 0, 0,
79 0, 0 },
80 { 0,
81 0,
82 0,
83 0,
84 0,
85 0,
86 0,
87 0,
88 0.025412,
89 0.018101,
90 0.0105117,
91 -0.000252347,
92 0.000857337,
93 0.00371809,
94 0.0022962,
95 0.00552758,
96 -0.000193383,
97 0.0016975,
98 0.0020368,
99 0.00152886,
100 0.00781011,
101 0.00445376,
102 0.00129068,
103 -0.000225775,
104 0.0105552,
105 0.0142455,
106 0.00369242,
107 0.014687,
108 0.00895646,
109 0.010251,
110 0.0174155,
111 0.0114122,
112 0.0157341,
113 0.00530172,
114 0.00958369,
115 0.00871203,
116 0.0110746,
117 0.00846772,
118 0.004005,
119 0.00648727,
120 0.00412489,
121 0.00510047,
122 0.0077632,
123 0.00452142,
124 0.00439873,
125 -0.00254842,
126 0,
127 0,
128 0,
129 0,
130 0,
131 0 },
132 { 0, 0, 0, 0, 0, 0.2, 0.0187074,
133 0.0151469, 0.0150374, 0.0149382, 0.00852575, 0.0117726, 0.00292745, 0.00594418,
134 0.00210557, 0.00437641, 0.0127742, 0, 0.00188013, -0.000209292, -0.000645717,
135 0.00506688, 0.00161251, 0.00951329, 0.0113891, 0.0154694, 0.011311, 0.0102122,
136 -0.0011918, 0.0115697, 0.0033551, 0.00681684, 0.0142289, 0.00486508, 0.0103071,
137 0.00585104, 0.00654329, 0.00773042, 0.0249313, 0.00100089, 0.00687448, 0.00108067,
138 0.00575115, 0.00836541, 0.00610641, 0.00498157, 0.0101283, -0.00449353, 0,
139 0, 0, 0 },
140 { 0, 0, 0, 0, 0.0357639, 0.0163919, 0.0267385,
141 0.0255677, 0.0161029, 0.0172624, 0.00385577, 0.00721942, 0.00383896, -0.000469486,
142 0.00519038, 0.00132785, 0.00298115, 0.00308214, 0.0060887, 0.00291204, 0.00383603,
143 0.0067211, 0.00317584, 0.00565859, 0.00959031, 0.00772591, 0.0165113, 0.00941265,
144 0.0105124, 0.0193347, 0.012998, 0.00561977, 0.00327781, 0.00355879, 0.0113581,
145 0.0114383, 0.0115121, 0.00388694, 0.00876026, 0.0113876, 0.00312005, 0.0143682,
146 0.0062588, 0.000176478, 0.0138092, 0.011377, 0.00383445, 0.0129828, 0.0272171,
147 0, 0, 0 },
148 { 0, 0, 0, 0, 0.0112905, 0.0116124, 0.0213804,
149 0.0328621, 0.0183141, 0.0246236, 0.00394498, 0.00887383, 0.00415104, -0.000519089,
150 0.00633283, 0.00132945, 0.00350115, 0.00362386, 0.00763285, 0.00291064, 0.00451877,
151 0.00697439, 0.00438741, 0.00652817, 0.0103673, 0.00837873, 0.022633, 0.00731777,
152 0.01128, 0.0180281, 0.0147394, 0.00613579, 0.00345585, 0.00418337, 0.0113416,
153 0.0148126, 0.011185, 0.00484846, 0.0102792, 0.0128383, 0.00288329, 0.0143363,
154 0.00623997, 0.000258956, 0.0151972, 0.0107407, 0.00460596, 0.0105749, 0.0123443,
155 0, 0, 0 },
156 { 0, 0, 0, -0.00280361, 0.00599815, 0.0124128, 0.026802,
157 0.0277765, 0.0203783, 0.0152604, 0.00563037, 0.00539587, 0.00742196, 0.00125699,
158 0.00518185, 0.00362133, 0.00502016, 0.00452122, 0.00639375, 0.00156148, 0.00331946,
159 0.0037184, 0.0030193, 0.00504048, 0.00936035, 0.00428231, 0.0157025, 0.00756299,
160 0.00826151, 0.0135952, 0.0104938, 0.00635791, 0.00549246, 0.00312205, 0.00749632,
161 0.0104138, 0.00878021, 0.00360666, 0.00967957, 0.00906651, 0.00563426, 0.0105905,
162 0.00818554, 0.00765972, 0.0114866, 0.00864717, 0.00973636, 0.00834703, 0.0127143,
163 0, 0, 0 },
164 { 0, 0, 0, 0.00690654, 0.0201675, 0.0178712,
165 0.0108822, 0.0112214, 0.0164782, 0.00218226, 0.00930132, 0.0135348,
166 0.00754681, 0.0178071, 0.00760194, 0.0131149, 0.00212195, 0.00508351,
167 0.00224026, 0.00228926, 0.00408006, 0.00238163, 0.00168934, 0.00388352,
168 0.00723827, 0.00615329, 0.0111326, 0.00904487, 0.00498642, -0.000847188,
169 0.00170438, 0.0126457, 0.0010792, 0.00472998, 0.00522703, 0.0115935,
170 0.0157492, 0.00747376, 0.00691549, 0.00766887, 0.00714786, 0.00519651,
171 0.00828873, 0.0199414, 0.00664845, 0.0152771, 0.0136218, 0.0189482,
172 0.0189, 0, 0, 0 },
173 { 0, 0, 0, 0.00402643, 0.0147216, 0.01715,
174 0.0120051, 0.0136828, 0.0133204, 0.00254233, 0.0104003, 0.0121178,
175 0.00545226, 0.0217436, 0.00735621, 0.0117699, 0.00172737, 0.00469086,
176 0.00251054, 0.00212634, 0.00446061, 0.00189206, 0.00206932, 0.00367898,
177 0.00638883, 0.00674313, 0.0114174, 0.0080861, 0.00369373, -0.000753843,
178 0.00182248, 0.0106247, 0.0010497, 0.00362544, 0.00541468, 0.0136716,
179 0.0148882, 0.00695346, 0.0065737, 0.00548947, 0.00620907, 0.00451994,
180 0.00673862, 0.0194927, 0.00707029, 0.0124406, 0.0129697, 0.0165352,
181 0.0247423, 0, 0, 0 },
182 { 0,
183 0,
184 0,
185 0,
186 -0.000274699,
187 0.00285875,
188 0.00658438,
189 0.0104937,
190 0.0193084,
191 0.00435312,
192 0.002284,
193 -0.00094086,
194 0.00681717,
195 0.00351616,
196 0.00649049,
197 0.00615866,
198 -0.000691215,
199 0.00167221,
200 0.00374339,
201 -0.000538108,
202 0.00450478,
203 0.00609494,
204 0.00136808,
205 -0.00170308,
206 0.00505613,
207 0.00796647,
208 -0.000908463,
209 0.00644384,
210 0.0009644,
211 0.00801776,
212 0.0138419,
213 0.0108086,
214 0.00331586,
215 0.0193975,
216 -0.00128345,
217 0.00831479,
218 0.00567967,
219 0.00248256,
220 0.00798811,
221 0.0113088,
222 0.00398216,
223 0.0146145,
224 0.000559689,
225 0.0136526,
226 0.00900583,
227 0.00765891,
228 0.0113426,
229 0,
230 0,
231 0,
232 0,
233 0 },
234 { 0,
235 0,
236 0,
237 0,
238 -0.000837011,
239 0.00783656,
240 0.0103523,
241 0.0127591,
242 0.014678,
243 0.0184914,
244 0.00536743,
245 0.00480503,
246 0.00397989,
247 0.0141369,
248 0.00160304,
249 0.00385771,
250 0.0148946,
251 0.00584356,
252 0.00591229,
253 -0.000800461,
254 0.00142273,
255 0.000380504,
256 0.00107338,
257 0.00532479,
258 0.0055054,
259 0.000956682,
260 -0.00130562,
261 0.0105865,
262 0.0106789,
263 0.00217481,
264 0.00547282,
265 -0.000182769,
266 -0.000466737,
267 0.0040426,
268 0.0077561,
269 0.0121245,
270 0.00192074,
271 0.00713558,
272 0.00700134,
273 0.00474334,
274 0.00742682,
275 0.00491013,
276 0.00400981,
277 0.00431014,
278 0.00592222,
279 0.00413596,
280 0.00991255,
281 -0.000899015,
282 0.0260504,
283 0,
284 0,
285 0 },
286 { 0,
287 0,
288 0,
289 0,
290 0,
291 0.0122997,
292 0.0169225,
293 0.0184353,
294 0.0164037,
295 0.0117505,
296 0.00480769,
297 0.0126396,
298 0.00399975,
299 0.00669493,
300 0.0118789,
301 0.00342785,
302 0.00141258,
303 0.00272821,
304 -0.000557028,
305 0.00156787,
306 0.00605785,
307 0.00222844,
308 0.000804709,
309 0.00897894,
310 0.00733436,
311 0.00605425,
312 -0.00170722,
313 0.00645383,
314 -0.0024872,
315 0.00638874,
316 0.000761802,
317 0.00587418,
318 0.00361459,
319 0.00200659,
320 0.0038677,
321 0.0114339,
322 0.00660603,
323 0.00888449,
324 0.0164075,
325 0.0104761,
326 0.0036128,
327 0.00721415,
328 0.00894716,
329 0.0200528,
330 0.0127048,
331 0.0114264,
332 0.0298287,
333 0.0122927,
334 0,
335 0,
336 0,
337 0 },
338 { 0, 0, 0, 0, 0, 0.00399151,
339 0.00615773, 0.0177878, 0.0189579, 0.010571, 0.0172364, -0.000927745,
340 0.00426211, 0.0109308, 0.00434394, 0.00467184, 0.0102868, 0.00455592,
341 0.00346033, 0.0011815, 0.00623132, 0.00405261, 0.000927003, 0.00308591,
342 0.00160743, 0.000821355, 0.0018395, 0.00390381, 0.00580204, -0.00146727,
343 0.00466621, 0.00100443, 0.00128368, 0.00196323, 0.00562257, 0.00326293,
344 0.00187777, 0.00104736, 0.00881352, 0.0143722, 0.00304789, 0.0209486,
345 0.00508492, 0.00432553, 0.0032542, 0.0253143, 0.0332302, 0.0152654,
346 0, 0, 0, 0 },
347 { 0,
348 0,
349 0,
350 0,
351 0,
352 0,
353 0.0202926,
354 0.0186309,
355 0.0104131,
356 0.00690488,
357 0.00236005,
358 0.00434038,
359 0.0142355,
360 0.0066333,
361 0.0109639,
362 0.0203775,
363 0.00361271,
364 0.0162799,
365 0.00223694,
366 0.00560054,
367 0.00378602,
368 0.00436651,
369 0.00919364,
370 0.00832168,
371 0.00484681,
372 0.00520677,
373 0.0132129,
374 0.00412628,
375 0.00188165,
376 0.00560339,
377 0.00905057,
378 0.00380511,
379 0.0101017,
380 -0.00105115,
381 -0.00106044,
382 0.0177763,
383 -0.000432789,
384 0.0078263,
385 0.00926294,
386 0.00649598,
387 0.0165854,
388 0.013363,
389 0.00388318,
390 0.027023,
391 0.0204139,
392 0.0245485,
393 0.0147441,
394 0.00703364,
395 0,
396 0,
397 0,
398 0 },
399 { 0,
400 0,
401 0,
402 0,
403 0,
404 0,
405 0.0106125,
406 0.0188372,
407 0.0146971,
408 0.00478025,
409 0.00663613,
410 0,
411 0.00601951,
412 0.0113968,
413 0.0259423,
414 0.00674669,
415 0.00506002,
416 0.00193793,
417 0.00393283,
418 0.00247295,
419 0.00721892,
420 0.00038463,
421 0.00737362,
422 0.00384348,
423 0.00494624,
424 0.00129381,
425 0.0121777,
426 0.00153305,
427 0,
428 0.00309061,
429 0.00096083,
430 0.00635241,
431 0.00437105,
432 0.00271235,
433 0.00323285,
434 0.0152463,
435 0.00790541,
436 0.00197237,
437 0.00842918,
438 0.00351487,
439 0.00691746,
440 0.0108789,
441 0.0061616,
442 0.0240904,
443 0.0129392,
444 0.0129275,
445 -0.000830457,
446 0,
447 0,
448 0,
449 0,
450 0 },
451 { 0, 0, 0, 0, 0, 0,
452 0, 0.0247863, 0.0176282, 0.0175289, 0.00457808, 0.00992452,
453 0.0115541, 0.0135212, 0.00994434, 0.00713981, 0.0111169, 0.0103148,
454 0.00416804, 0.00410405, 0.00111858, 0.00873464, 0.00331374, 0.0106796,
455 0.00201416, 0.00122554, 0.00385263, 0.00661668, 0.00463202, -0.000606177,
456 -0.00120285, 0.00754939, 0.00930348, 0.00377447, 0.00465278, 0.00101872,
457 0.00231525, 0.0112804, 0.00641832, 0.00746406, 0.00930169, 0.0118229,
458 0.00592735, 0.0206328, 0.0224974, 0.0115758, 0, 0,
459 0, 0, 0, 0 },
460 { 0,
461 0,
462 0,
463 0,
464 0,
465 0,
466 0,
467 0,
468 0.0136459,
469 0.00712392,
470 0.0138838,
471 0.00159252,
472 -0.000765476,
473 0.010438,
474 0.0126405,
475 0.0152418,
476 0.0034147,
477 0.00790645,
478 0.00381993,
479 0.00428926,
480 0.00622248,
481 0.00742261,
482 0.00232875,
483 0.00613106,
484 0.00168197,
485 0.00352778,
486 0.00934557,
487 0.00150865,
488 0.00118901,
489 0.00179582,
490 0.00168158,
491 0.00214762,
492 0.00615021,
493 0.00110451,
494 -0.000977385,
495 -0.000972892,
496 0.00159201,
497 -0.00237044,
498 0.00777439,
499 0.00259004,
500 0.00859064,
501 0.00659134,
502 0.0016455,
503 0.00686578,
504 0.00407572,
505 0.0201268,
506 0,
507 0,
508 0,
509 0,
510 0,
511 0 },
512 { 0,
513 0,
514 0,
515 0,
516 0,
517 0,
518 0,
519 0,
520 0.0190655,
521 0.0193906,
522 0.0298653,
523 0.0156586,
524 0.00188532,
525 0.0140794,
526 0.00497449,
527 0.00465134,
528 0.00547378,
529 -0.00155504,
530 0.0112145,
531 0.00138882,
532 0.0100356,
533 0.00871306,
534 0.00347013,
535 0.000565267,
536 -0.00168501,
537 0.00220346,
538 0.00541041,
539 -0.000871444,
540 0.00192547,
541 -0.000892462,
542 0.0015037,
543 0.000526234,
544 0.00204403,
545 0.00432735,
546 0.00207964,
547 -0.00125049,
548 0.00368116,
549 0.00582942,
550 0.00563678,
551 -2.20906e-05,
552 -0.000726604,
553 0.000425052,
554 0.00823579,
555 0.00334216,
556 0.0107216,
557 0.0110852,
558 0,
559 0,
560 0,
561 0,
562 0,
563 0 },
564 { 0,
565 0,
566 0,
567 0,
568 0,
569 0,
570 0,
571 0,
572 0.0102941,
573 0.0227606,
574 0.00958853,
575 0.0112508,
576 0.00437405,
577 0.019164,
578 0.0107376,
579 0.0073221,
580 0.0124903,
581 0.0161441,
582 0.000924579,
583 0.00352839,
584 0.0065253,
585 0.000840571,
586 0.00522096,
587 0.00320728,
588 0.00169485,
589 -0.0014829,
590 -0.000385922,
591 0.00380332,
592 0.00130042,
593 0.003991,
594 -0.00152388,
595 -0.000335786,
596 0.00115229,
597 0.00738576,
598 -0.000475654,
599 0.00231245,
600 -0.00128109,
601 0.00129967,
602 0.00403285,
603 -0.00104735,
604 0.00164014,
605 0.0088376,
606 0.00717689,
607 0.000210979,
608 -0.000834027,
609 -0.00523533,
610 0,
611 0,
612 0,
613 0,
614 0,
615 0 },
616 { 0,
617 0,
618 0,
619 0,
620 0,
621 0,
622 0,
623 0,
624 0,
625 0.00906013,
626 0.00850382,
627 0.00337777,
628 0.00474848,
629 0.00501928,
630 0.00262243,
631 0.00361578,
632 0.00412986,
633 0.00828773,
634 0.00317526,
635 0.000418765,
636 0.00571304,
637 0.0039906,
638 0.00398882,
639 0.00315977,
640 0.00154839,
641 0.0053812,
642 -0.00215534,
643 0.00814652,
644 0.00199931,
645 0.00156904,
646 0.00225004,
647 -0.000790613,
648 0.00182555,
649 0.00237361,
650 0,
651 0.00425822,
652 -0.00124752,
653 -0.00125848,
654 -0.000457286,
655 0.00132635,
656 0.00158466,
657 0.00165806,
658 0.00104607,
659 0.00136616,
660 0.00464344,
661 0,
662 0,
663 0,
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.0154465,
678 0.00895089,
679 0.00729578,
680 0.00496513,
681 0.0110113,
682 0.00404945,
683 0.00868326,
684 0.00803355,
685 0.00121479,
686 0.000957591,
687 -0.00140863,
688 0.0111029,
689 0.00542703,
690 0.00315061,
691 0.0142478,
692 0.00984325,
693 0.00326606,
694 0.00141415,
695 -0.000485355,
696 0.00169373,
697 0.00175355,
698 -0.000214901,
699 0.00382858,
700 -0.000839424,
701 -0.00041114,
702 0.00122857,
703 -0.000428686,
704 -0.000231425,
705 -0.000484801,
706 -0.000677427,
707 -0.00107819,
708 0.00398122,
709 0.00389996,
710 0.00160107,
711 -0.000596897,
712 -0.00473512,
713 0,
714 0,
715 0,
716 0,
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.00557644,
731 0.0117601,
732 -2.014e-05,
733 0.00532312,
734 0.00164892,
735 0.00123443,
736 0.000548319,
737 0.00379273,
738 0.000819126,
739 0.00315643,
740 -0.00112786,
741 0.00447013,
742 0.0039469,
743 -0.000698325,
744 0.00637213,
745 -0.000918547,
746 -0.000877958,
747 -0.000192955,
748 0.00386532,
749 0.00108823,
750 0,
751 0.00216934,
752 -0.000612477,
753 -0.000201423,
754 0,
755 0.00195611,
756 0.00260219,
757 -0.000438483,
758 -0.000499695,
759 -0.00107599,
760 -0.00051163,
761 -0.000468805,
762 -0.00136743,
763 -0.000565784,
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.0113123,
784 0.00768323,
785 0.00389946,
786 0.00981506,
787 0.00978708,
788 0.0127364,
789 0.000633142,
790 0.00285515,
791 0.00344821,
792 0.00683397,
793 -0.000793583,
794 0.00484872,
795 0.00140184,
796 0.00287053,
797 0.00122665,
798 0.0019104,
799 0.00146129,
800 0.00395712,
801 0.00180566,
802 -0.00142892,
803 -0.000177317,
804 -0.000627594,
805 0.00208295,
806 0,
807 -0.000770852,
808 -0.00199946,
809 -0.000717908,
810 -0.000867859,
811 -0.00107504,
812 -0.00128431,
813 -0.000493855,
814 -0.000265389,
815 -0.0104253,
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.0101946,
837 0.00354491,
838 0.00156002,
839 0.00514551,
840 0.00414005,
841 0.00259554,
842 0.00681059,
843 0.018554,
844 0.00333984,
845 0.005891,
846 0.00563848,
847 0.00506756,
848 0.00172031,
849 0.0014715,
850 0.0049191,
851 -0.00172828,
852 0.00219327,
853 -0.000623313,
854 0.000879113,
855 0.00240745,
856 -0.00113367,
857 -0.000812239,
858 0.00148996,
859 -0.000525341,
860 0.00365075,
861 -0.00108374,
862 -0.000710412,
863 0.00358373,
864 -0.00126742,
865 -0.0012819,
866 0,
867 0,
868 0,
869 0,
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.00342704,
889 0.00486663,
890 0.00801254,
891 0.00482659,
892 0.0067724,
893 0.000801564,
894 0.00216852,
895 0.00120544,
896 0.00117051,
897 0.00292921,
898 0.00293786,
899 0.00435708,
900 0.00336549,
901 0.00184459,
902 0.00356918,
903 -0.00143222,
904 0.00271472,
905 -0.000847567,
906 -0.00117673,
907 -0.000591092,
908 -0.000920364,
909 -0.00162241,
910 0.00213705,
911 -0.000244555,
912 -0.000452254,
913 0.000957899,
914 -0.000852445,
915 -0.000524743,
916 0.00174564,
917 0.00208469,
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.01942,
941 0.0104793,
942 0.00238674,
943 -0.00028427,
944 0.00323803,
945 0.00428694,
946 0.0180191,
947 0.00205732,
948 0.004827,
949 0.013351,
950 0.0164632,
951 0.00661175,
952 0.00625597,
953 0.002846,
954 0.00201106,
955 0.000985679,
956 -0.000678121,
957 -0.000691974,
958 0.00196908,
959 0.000656031,
960 -0.00179116,
961 0.000615601,
962 0.00184952,
963 -0.000819999,
964 0.00320967,
965 -0.00193909,
966 -0.00105235,
967 -0.00205237,
968 0.00162034,
969 0.00482908,
970 0,
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.0250022,
994 0.00643193,
995 0.0145664,
996 0.0105645,
997 0.005722,
998 0.0114656,
999 0.0105034,
1000 0.00863349,
1001 0.00663601,
1002 0.00534037,
1003 0.0083902,
1004 -2.58751e-05,
1005 0.00179951,
1006 0.00120073,
1007 -0.000662521,
1008 -0.000201754,
1009 -0.000705186,
1010 0.0016609,
1011 -0.000236756,
1012 0.00330577,
1013 0.00504613,
1014 0.000588122,
1015 0.00174472,
1016 -0.000951208,
1017 0.00253816,
1018 0.000456801,
1019 0.0059189,
1020 -0.0011968,
1021 -0.00077095,
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.00243762,
1046 0.0101053,
1047 0.0137352,
1048 0.00989112,
1049 0.00419056,
1050 0.0167929,
1051 -0.00150345,
1052 0.00391754,
1053 0.00962137,
1054 0.0111504,
1055 0.013032,
1056 0.0091798,
1057 0.00632917,
1058 0.00406742,
1059 0.00171305,
1060 0.0030055,
1061 0.00404168,
1062 -0.00197613,
1063 0.00241273,
1064 0.00871242,
1065 0.00473052,
1066 0.00142909,
1067 0.00582945,
1068 0.00698106,
1069 0.00411218,
1070 0.00363546,
1071 0.00151842,
1072 -0.000477552,
1073 0,
1074 0,
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.0137108,
1099 0.0101587,
1100 0.00512209,
1101 0.00581154,
1102 0.0063175,
1103 0.0148431,
1104 0.00784464,
1105 0.0153501,
1106 0.00301844,
1107 0.0187667,
1108 0.007718,
1109 0.00198616,
1110 0.00344126,
1111 0.000856434,
1112 0.006128,
1113 0.00407622,
1114 0.0034297,
1115 0.000401058,
1116 -0.000712372,
1117 -0.000616223,
1118 -0.000855106,
1119 0.00335076,
1120 -0.000459257,
1121 -0.00100754,
1122 -0.000922785,
1123 0.00190541,
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,
1151 0.0115846,
1152 0.00513737,
1153 0.00809827,
1154 0.00997606,
1155 0.00655055,
1156 0.00941147,
1157 0.0038603,
1158 0.00929101,
1159 0.00884371,
1160 0.0032594,
1161 -0.000438409,
1162 -0.000647537,
1163 0.00152629,
1164 -0.00160462,
1165 0.000699107,
1166 0.00165622,
1167 -0.00124744,
1168 -0.000977481,
1169 0.0014039,
1170 0.00174867,
1171 0.00401864,
1172 0.00399156,
1173 0.00201647,
1174 -0.000601122,
1175 -0.00247459,
1176 0,
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,
1203 -0.00105089,
1204 0.0048501,
1205 0.00845776,
1206 0.00199121,
1207 0.00916815,
1208 0.0118089,
1209 0.0101806,
1210 0.0137388,
1211 0.0176601,
1212 0.00328823,
1213 0.00352053,
1214 0.00419714,
1215 0.00546382,
1216 0,
1217 0.00161903,
1218 -0.00119964,
1219 0.00191271,
1220 -0.000415298,
1221 0.00674503,
1222 0.00439691,
1223 -0.000762733,
1224 -0.000813587,
1225 0.00370975,
1226 0.00708438,
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,
1255 0,
1256 0.0145711,
1257 0.0094525,
1258 0.00818695,
1259 -0.000665325,
1260 0.0047514,
1261 0.00368851,
1262 0.00805664,
1263 0.00446138,
1264 0.00914161,
1265 0.00622046,
1266 0.00112126,
1267 0.000882917,
1268 0.0016548,
1269 0.00429119,
1270 -0.00142295,
1271 -2.06638e-05,
1272 -0.00157297,
1273 -0.000195611,
1274 0.00347694,
1275 -0.000951127,
1276 0.00278539,
1277 0.00166737,
1278 -0.000765855,
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,
1306 0,
1307 0,
1308 0,
1309 0.0124235,
1310 0.00144789,
1311 0.00820455,
1312 0.0031008,
1313 0.00870771,
1314 0.00334443,
1315 0.0121713,
1316 0.00201216,
1317 0.00607529,
1318 0.00177875,
1319 0.00140863,
1320 0.00199694,
1321 0.00242735,
1322 0.00139591,
1323 0.00121349,
1324 0,
1325 -0.00124997,
1326 -0.000215363,
1327 0.00361555,
1328 -0.000206134,
1329 -0.000716385,
1330 0,
1331 0,
1332 0,
1333 0,
1334 0,
1335 0,
1336 0,
1337 0,
1338 0,
1339 0,
1340 0,
1341 0,
1342 0,
1343 0 },
1344 { 0, 0, 0, 0, 0, 0,
1345 0, 0, 0, 0, 0, 0,
1346 0, 0, 0, 0, 0, 0.00128829,
1347 0.00342748, 0.00818716, 0.00565558, 0.00871036, 0.0138593, 0.00145614,
1348 0.00132705, 0.0106219, 0.00618985, 0.00187038, -0.00137881, 0.0018989,
1349 0, 0.00185971, 0.00161027, 0, 0, 0.00198384,
1350 0.0021881, -0.00316355, 0, 0, 0, 0,
1351 0, 0, 0, 0, 0, 0,
1352 0, 0, 0, 0 },
1353 { 0,
1354 0,
1355 0,
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.00148772,
1371 0.00977291,
1372 0.00858047,
1373 0.00392799,
1374 0.00286997,
1375 0.00420178,
1376 0.0063399,
1377 0.00186727,
1378 0.00875324,
1379 0.00169687,
1380 0.00182879,
1381 0.00501749,
1382 0,
1383 -0.000424115,
1384 -0.000195288,
1385 -0.00132863,
1386 0.00201768,
1387 -0.000723473,
1388 -0.000586035,
1389 -0.00211294,
1390 0,
1391 0,
1392 0,
1393 0,
1394 0,
1395 0,
1396 0,
1397 0,
1398 0,
1399 0,
1400 0,
1401 0,
1402 0,
1403 0,
1404 0 },
1405 { 0, 0, 0, 0, 0,
1406 0, 0, 0, 0, 0,
1407 0, 0, 0, 0, 0,
1408 0, 0, 0, 0.00795127, 0.000761631,
1409 0.0130024, 0.00888905, 0.00533266, 0.00165759, -0.000681874,
1410 0.00613109, 0.00421136, 0.0026058, -0.00100585, -0.000440221,
1411 -0.00102298, 0.00181682, 0.00229075, 0.00159399, -0.000890541,
1412 0.0015385, 0, 0, 0, 0,
1413 0, 0, 0, 0, 0,
1414 0, 0, 0, 0, 0,
1415 0, 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.00176391,
1436 0.0016436,
1437 0.00778176,
1438 0.0019574,
1439 0.00373541,
1440 0.00121501,
1441 0.00194797,
1442 0.00156388,
1443 0,
1444 0.00205196,
1445 0,
1446 0.00194588,
1447 -0.000447248,
1448 0.00164912,
1449 0.00378861,
1450 -0.000967688,
1451 -0.000535921,
1452 0,
1453 0,
1454 0,
1455 0,
1456 0,
1457 0,
1458 0,
1459 0,
1460 0,
1461 0,
1462 0,
1463 0,
1464 0,
1465 0,
1466 0,
1467 0 },
1468 { 0,
1469 0,
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.000414658,
1488 0.00158378,
1489 0.00792649,
1490 0.00747993,
1491 0.00355962,
1492 -0.000222224,
1493 -0.000883339,
1494 -0.000675862,
1495 0.00445083,
1496 -0.00039295,
1497 -0.000978694,
1498 0.00126864,
1499 0.00130999,
1500 -0.000655967,
1501 0,
1502 -0.000605564,
1503 -0.00786369,
1504 0,
1505 0,
1506 0,
1507 0,
1508 0,
1509 0,
1510 0,
1511 0,
1512 0,
1513 0,
1514 0,
1515 0,
1516 0,
1517 0,
1518 0,
1519 0 },
1520 { 0,
1521 0,
1522 0,
1523 0,
1524 0,
1525 0,
1526 0,
1527 0,
1528 0,
1529 0,
1530 0,
1531 0,
1532 0,
1533 0,
1534 0,
1535 0,
1536 0,
1537 0,
1538 0,
1539 0,
1540 0.00534462,
1541 0.00337856,
1542 0.00447068,
1543 0.00202555,
1544 0.00262315,
1545 -0.000791716,
1546 0.0017379,
1547 -0.000440789,
1548 -0.000377595,
1549 0.00199316,
1550 -0.000205376,
1551 0.00113595,
1552 0.000595314,
1553 0,
1554 -0.00101937,
1555 0,
1556 0,
1557 0,
1558 0,
1559 0,
1560 0,
1561 0,
1562 0,
1563 0,
1564 0,
1565 0,
1566 0,
1567 0,
1568 0,
1569 0,
1570 0,
1571 0 },
1572 { 0,
1573 0,
1574 0,
1575 0,
1576 0,
1577 0,
1578 0,
1579 0,
1580 0,
1581 0,
1582 0,
1583 0,
1584 0,
1585 0,
1586 0,
1587 0,
1588 0,
1589 0,
1590 0,
1591 0,
1592 -0.00194879,
1593 0.0043191,
1594 0.00349993,
1595 0.00656448,
1596 -0.00022874,
1597 0.00490897,
1598 0.00173958,
1599 0.00211379,
1600 -0.000919604,
1601 0.00122627,
1602 0.00192831,
1603 0.000924648,
1604 -0.000533183,
1605 0,
1606 0,
1607 0,
1608 0,
1609 0,
1610 0,
1611 0,
1612 0,
1613 0,
1614 0,
1615 0,
1616 0,
1617 0,
1618 0,
1619 0,
1620 0,
1621 0,
1622 0,
1623 0 },
1624 { 0,
1625 0,
1626 0,
1627 0,
1628 0,
1629 0,
1630 0,
1631 0,
1632 0,
1633 0,
1634 0,
1635 0,
1636 0,
1637 0,
1638 0,
1639 0,
1640 0,
1641 0,
1642 0,
1643 0,
1644 0,
1645 0.00550371,
1646 0.00377472,
1647 0.00180269,
1648 0.0031602,
1649 -0.000686968,
1650 -0.000478745,
1651 -0.000771646,
1652 -0.00125527,
1653 0.00573005,
1654 0.00298008,
1655 0.00219993,
1656 0.00456167,
1657 0.00815102,
1658 0,
1659 0,
1660 0,
1661 0,
1662 0,
1663 0,
1664 0,
1665 0,
1666 0,
1667 0,
1668 0,
1669 0,
1670 0,
1671 0,
1672 0,
1673 0,
1674 0,
1675 0 },
1676 { 0,
1677 0,
1678 0,
1679 0,
1680 0,
1681 0,
1682 0,
1683 0,
1684 0,
1685 0,
1686 0,
1687 0,
1688 0,
1689 0,
1690 0,
1691 0,
1692 0,
1693 0,
1694 0,
1695 0,
1696 0,
1697 0,
1698 -0.00171234,
1699 0.00450387,
1700 -0.000745589,
1701 0.000928437,
1702 0.00198709,
1703 -0.00020945,
1704 -0.000248579,
1705 0.00112635,
1706 0.00215338,
1707 0.0066882,
1708 0.0106831,
1709 0,
1710 0,
1711 0,
1712 0,
1713 0,
1714 0,
1715 0,
1716 0,
1717 0,
1718 0,
1719 0,
1720 0,
1721 0,
1722 0,
1723 0,
1724 0,
1725 0,
1726 0,
1727 0 },
1728 { 0,
1729 0,
1730 0,
1731 0,
1732 0,
1733 0,
1734 0,
1735 0,
1736 0,
1737 0,
1738 0,
1739 0,
1740 0,
1741 0,
1742 0,
1743 0,
1744 0,
1745 0,
1746 0,
1747 0,
1748 0,
1749 0,
1750 -0.00106871,
1751 -0.000222508,
1752 0.00205684,
1753 -0.00038869,
1754 0.00157722,
1755 -0.000187058,
1756 -0.000198706,
1757 0.00876609,
1758 0.00483817,
1759 0.00529853,
1760 0.0143598,
1761 0,
1762 0,
1763 0,
1764 0,
1765 0,
1766 0,
1767 0,
1768 0,
1769 0,
1770 0,
1771 0,
1772 0,
1773 0,
1774 0,
1775 0,
1776 0,
1777 0,
1778 0,
1779 0 },
1780 { 0, 0, 0, 0, 0, 0, 0, 0,
1781 0, 0, 0, 0, 0, 0, 0, 0,
1782 0, 0, 0, 0, 0, 0, 0, 0.00220188,
1783 -0.00065229, 0.00168199, 0.00391551, 0, 0.00777458, 0.00353025, 0, -0.000506828,
1784 0, 0, 0, 0, 0, 0, 0, 0,
1785 0, 0, 0, 0, 0, 0, 0, 0,
1786 0, 0, 0, 0 },
1787 { 0, 0, 0, 0, 0, 0, 0,
1788 0, 0, 0, 0, 0, 0, 0,
1789 0, 0, 0, 0, 0, 0, 0,
1790 0, 0, 0, 0.00189452, 0.00193445, 0.00152572, -0.000211959,
1791 0.00383985, 0.00874946, 0.00457722, 0.0299003, 0, 0, 0,
1792 0, 0, 0, 0, 0, 0, 0,
1793 0, 0, 0, 0, 0, 0, 0,
1794 0, 0, 0 },
1795 { 0, 0, 0, 0, 0,
1796 0, 0, 0, 0, 0,
1797 0, 0, 0, 0, 0,
1798 0, 0, 0, 0, 0,
1799 0, 0, 0, 0, -0.000708505,
1800 -0.0011835, 0.00362797, -0.00112053, 0.0057435, 0.00834591,
1801 0.0112548, 0, 0, 0, 0,
1802 0, 0, 0, 0, 0,
1803 0, 0, 0, 0, 0,
1804 0, 0, 0, 0, 0,
1805 0, 0 },
1806 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1807 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00433094,
1808 -0.00162684, 0.00649451, 0.00462054, 0.0141132, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1809 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1810 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1811 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0139919,
1812 -0.00152095, 0.0163195, 0.010359, 0.0186132, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1813 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1814 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1815 0, 0, 0, 0, 0, 0, 0, 0.0147059, 0.00529507, 0.0173027, 0.00279661, 0, 0, 0, 0, 0, 0, 0,
1816 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1817 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1818 0, 0, 0, 0, 0, 0, 0, 0, 0.00801884, 0.0144967, 0, 0, 0, 0, 0, 0, 0, 0,
1819 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1820 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1821 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00296133, 0, 0, 0, 0, 0, 0, 0, 0,
1822 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
1823 { 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,
1824 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 } };
1825
1826 for ( int i = 0; i < 52; i++ )
1827 {
1828 for ( int j = 0; j < 52; j++ ) { HisPDF[i][j] = HisPDFtmp[i][j]; }
1829 }
1830}
1831
1833
1835
1836loop:
1838
1839 EvtParticle *id1, *id2, *id3;
1840 EvtVector4R pd1, pd2, pd3;
1841 double xmass13, xmass12, xmass23;
1842
1843 id1 = p->getDaug( 0 );
1844 id2 = p->getDaug( 1 );
1845 id3 = p->getDaug( 2 );
1846
1847 pd1 = id1->getP4Lab();
1848 pd2 = id2->getP4Lab();
1849 pd3 = id3->getP4Lab();
1850
1851 // xmass12=(pd1+pd2).mass()*(pd1+pd2).mass(); // M_ksopi
1852 xmass13 = ( pd1 + pd3 ).mass() * ( pd1 + pd3 ).mass(); // M_ksow
1853 xmass23 = ( pd2 + pd3 ).mass() * ( pd2 + pd3 ).mass(); // M_piw
1854
1855 int xbin = FindXBin( xmass23 );
1856 int ybin = FindYBin( xmass13 );
1857 double xratio12 = HisPDF[xbin][ybin] / avm1;
1858
1859 if ( xratio12 <= 0 ) goto loop;
1860
1861 double rd12 = EvtRandom::Flat( 0.0, 1.0 );
1862 if ( rd12 > xratio12 ) goto loop;
1863
1864 return;
1865}
1866
1867int EvtDtopipipiPlot::FindXBin( double mass2 ) {
1868 if ( mass2 < Xmin ) { return 0; }
1869 else if ( mass2 >= Xmax ) { return 51; }
1870 else { return int( ( mass2 - Xmin ) / Xwid ) + 1; }
1871}
1872
1873int EvtDtopipipiPlot::FindYBin( double mass2 ) {
1874 if ( mass2 < Ymin ) { return 0; }
1875 else if ( mass2 >= Ymax ) { return 51; }
1876 else { return int( ( mass2 - Ymin ) / Ywid ) + 1; }
1877}
double mass
EvtId getParentId()
EvtId * getDaugs()
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
virtual ~EvtDtopipipiPlot()
int FindYBin(double mass2)
EvtDecayBase * clone()
void decay(EvtParticle *p)
int FindXBin(double mass2)
void getName(std::string &name)
static EvtSpinType::spintype getSpinType(EvtId i)
Definition EvtPDL.hh:66
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:272
EvtVector4R getP4Lab()
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
static double Flat()
Definition EvtRandom.cc:69