BOSS 7.1.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 "EvtGenBase/EvtPatches.hh"
19#include "EvtGenBase/EvtParticle.hh"
20#include "EvtGenBase/EvtGenKine.hh"
21#include "EvtGenBase/EvtPDL.hh"
22#include "EvtGenBase/EvtReport.hh"
23#include "EvtGenModels/EvtDtopipi0pi0Plot.hh"
24#include "EvtGenBase/EvtRandom.hh"
25#include <stdlib.h>
26#include <string>
27
28#include "TFile.h"
29#include "TApplication.h"
30#include "TROOT.h"
31#include "TH1.h"
32#include "TAxis.h"
33#include "TH2.h"
34using std::endl;
35
36EvtDtopipi0pi0Plot::~EvtDtopipi0pi0Plot() {}
37
38void EvtDtopipi0pi0Plot::getName(std::string& model_name){
39 model_name="Dtopipi0pi0Plot";
40}
41
42EvtDecayBase* EvtDtopipi0pi0Plot::clone(){
43 return new EvtDtopipi0pi0Plot;
44}
45
46void EvtDtopipi0pi0Plot::init(){
47
48 checkNArg(0);
49
50 bool idN = getDaugs()[0]==EvtPDL::getId(std::string("pi+"))||getDaugs()[0]==EvtPDL::getId(std::string("pi-"));
51 bool idKs = getDaugs()[1]==EvtPDL::getId(std::string("pi0"));
52 bool idPi = getDaugs()[2]==EvtPDL::getId(std::string("pi0"));
53 if(!(idN && idKs && idPi ) ){std::cout<<"EvtDtopipi0pi0Plot: the daughter sequence should be pi+ pi0 pi0"<<std::endl;abort();}
54 EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId());
55
56Xmin = 0.06887;
57Xmax = 3.0518;
58Xwid = 0.0596586;
59Ymin = 0.07037;
60Ymax = 1.75;
61Ywid = 0.0335926;
62avm1 = 0.25;
63double HisPDFtmp[52][52] =
64{
65{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, 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},
66{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.126132, 0.0220555, 0.0164693, 0.0143236, 0.0615541, 0.0238032, 0.0104177, 0.0514766, 0.00561927, -0.00970829, 0.015674, 0.033324, -0.0019576, 0.0179617, 0.00779605, 0.0474916, 0.00439088, 0.0191456, 0.00837141, 0.00609675, 0.0840422, 0.0270903, 0.0136279, 0.00854272, 0.0278523, 0.00389946, 0.00795048, 0.0450092, 0.0176549, 0.0109168, 0.0398448, 0.0144404, 0.0657468, 0},
67{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.0122324, -0.00324092, 0.0196607, 0.00102101, 0.00378296, 0.00634048, 0.0217871, 0.0081836, 0.0140993, 0.00592363, 0.0289073, 0.0116091, 0.015295, 0.0232406, 0.016491, 0.0115485, 0.0147406, 0.00684889, 0.00932597, 0.0154916, 0.0124573, 0.0166611, 0.0128406, 0.0321747, 0.00527513, 0.011663, 0.0212437, 0.000725348, 0.030167, 0.0213951, 0.0147829, 0.0065729, 0.0174563, 0.0153622, 0.00786071, 0.0380956, 0.0185521, 0.0102917, 0.0206631, 0.0569558, 0, 0},
68{0, 0, 0, 0, 0, 0, 0, 0.181193, -0.00514181, -0.00212754, -0.00600245, -0.0047888, 0.00832803, 0.0150015, 0.00962087, 0.0126662, 0.0285356, 0.0107688, 0.0164774, 0.0125525, 0.0181543, 0.0127544, 0.0180501, 0.0221408, 0.0124838, 0.0117366, 0.0103621, 0.00569093, 0.0156799, 0.0104127, 0.0134806, 0.0109779, 0.00910387, 0.01754, 0.0138031, 0.0108085, 0.0151902, 0.00519298, 0.0126952, 0.0127009, 0.0148566, 0.00911647, 0.0103495, 0.0094019, 0.00318659, 0.026405, 0.0122733, 0.00718738, 0.0457039, 0, 0, 0},
69{0, 0, 0, 0, 0, 0, 0.00534517, 0.0171958, -0.00418114, -0.00326283, -0.0103725, -0.00985988, 0.0125388, 0.0312826, 0.0177929, 0.0294284, 0.0335579, 0.0178685, 0.018411, 0.025134, 0.0116659, 0.0108364, 0.0207183, 0.0159755, 0.0143117, 0.0141057, 0.00626057, 0.00416359, 0.0142275, 0.00944204, 0.00960788, 0.0090878, 0.00310312, 0.0050442, 0.0177199, 0.0104969, 0.00357364, 0.00960106, -0.00060251, 0.00458286, 0.0120471, 0.0169905, 0.00384143, 0, -0.00054573, 0.0156392, 0.0182633, 0, 0, 0, 0, 0},
70{0, 0, 0, 0, 0, -0.0063761, 0.00240996, 0.016717, -0.00403196, -0.00315398, -0.0079351, -0.00950213, 0.0108979, 0.0257102, 0.0164158, 0.0169603, 0.039045, 0.0156618, 0.0151494, 0.0241756, 0.0121646, 0.014526, 0.0206264, 0.0169852, 0.0141894, 0.0114929, 0.00528244, 0.00348741, 0.0119654, 0.00944351, 0.010706, 0.010249, 0.00324893, 0.00523534, 0.0156045, 0.012237, 0.00412841, 0.00945598, -0.000586345, 0.00497353, 0.00891827, 0.016972, 0.00367533, 0, -0.000588499, 0.0145899, 0.0285661, 0, 0, 0, 0, 0},
71{0, 0, 0, 0, 0.0214843, 0.0152376, 0.00361878, 0.00659047, 0.0015044, 0.00383999, 0.015594, 0.0151121, 0.0021722, 0.00897155, 0.00140851, 0.0175562, 0.00615938, 0.0313451, 0.00244723, 0.0040116, 0.0150819, 0.00311473, 0.00947754, 0, 0.0102011, -0.00156026, 0.0107358, 0.00302003, 0.0104483, -0.00228917, -0.000471791, 0.0089698, -0.000910742, 0.0118354, 0.0047521, 0.00366432, -0.0004382, 0.00986424, 0.0137893, 0.015136, -0.00208091, 0.0104049, 0.0257901, 0, 0.013116, -0.00358994, 0.0685484, 0, 0, 0, 0, 0},
72{0, 0, 0, 0, -0.00226632, 0.00771791, -0.00427455, 0.0117554, -0.00300342, 0.00693718, 0.0189155, 0.00203017, 0.00979611, -0.00473465, 0.0173956, 0.0131809, 0.0128582, 0.00145764, 0.00549767, 0.00423976, 0.00889498, 0.0105157, 0.009075, 0.00134348, -0.000548315, 0.00165993, 0, 0.0108027, 0.0228641, 0.00539549, -0.00174531, -0.000652211, 0.00788748, 0.00229044, 0.01005, 0.0107275, 0.0084291, 0.00547439, 0.00902564, -0.000925959, 0.0177774, -0.00256084, 0.0130949, 0, -0.00222004, 0, 0, 0, 0, 0, 0, 0},
73{0, 0, 0, 0.00550887, 0.00193808, 0.00608034, 0.0181517, -0.00603111, -0.00472705, 0.0131492, 0.00704957, 0.0114897, 0.00189524, 0.00536718, 0.0189741, 0.0298482, 0.00440453, 0.00249411, 0.0229042, 0.0121403, 0.0107645, 0.00982848, 0.00466503, 0.0182887, -0.00123724, -0.000535069, 0.0112898, 0.0043009, -0.000457878, 0.00370862, 0.0168827, 0.0104902, 0.00383087, 0.00505172, 0.0046832, 0.016317, 0.00371058, 0.00429451, 0, 0.00449643, 0.00984612, 0, 0, -0.00110221, 0, 0, 0, 0, 0, 0, 0, 0},
74{0, 0, 0, 0.00559425, 0.00784236, -0.00515694, -0.00341012, 0.0056297, 0.00161726, 0.00817172, 0.0135242, 0.00792446, 0.0130368, 0.00220852, 0.00213705, 0.0213428, 0.00845069, -0.00640989, 0.00147561, -0.00334377, 0.00396805, -0.00163827, 0.015181, -0.00160953, 0.00974975, -0.000859794, 0.0071086, 0.00382806, -0.00327447, -0.00100732, 0.0095671, 0.00985691, 0.00516776, -0.000591646, 0.0052562, -0.0020691, 0.0125954, 0.00530462, 0.00524183, -0.00200351, 0.00427751, 0.00629911, -0.00116534, 0, 0, 0, 0, 0, 0, 0, 0, 0},
75{0, 0, 0.00924787, 0.0220888, 0.00547255, -0.00053368, -0.00621351, 0.00340835, 0.000869906, -0.00453506, -0.00230228, 0.00977444, 0.00150216, 0.0126651, 0.00664757, 0.00845056, 0.0145116, 0.0219415, 0.00305964, 0.0159933, 0.00209908, -0.00145665, -0.00114919, 0.00413543, 0.0037937, 0.0197595, -0.000571623, 0.00344037, 0, 0.0110179, 0.00531177, 0.00501057, -0.001035, 0.00452493, 0.0103877, 0.00574533, 0, 0.00384947, -0.00148758, 0.00478215, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
76{0, 0, 0.00740427, 0.00720438, 0.00554212, 0.00468727, 0.00749359, 0.0020956, 0.00530923, -0.00152604, -0.00157443, 0.0066342, 0.00491629, 0.0124885, 0.0119529, 0.016968, 0.00564861, 0.00142758, 0.0032944, 0.0125377, 0.0107067, -0.0025167, 0.00932271, 0.0151918, 0, -0.00107743, -0.00102267, 0.0103105, 0.00630324, -0.00215053, 0.00358067, 0.00498285, 0.0187526, 0, 0, 0.00485819, 0, 0.0109325, 0.00526438, 0.0087764, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
77{0, 0, -0.00188838, -0.00136398, 0.00841484, 0.00796693, 0.023253, -0.00402091, 0.012051, 0.0202265, 0.0140749, 0.0122277, 0.01082, 0.00468877, 0.00300136, 0.0257563, 0.0326677, 0.00203623, 0.00314237, 0.0101705, 0.00789608, -0.000500071, 0.00266911, 0.00330721, 0.00524247, 0.0163683, 0, 0.00357269, -0.000856078, 0.00566965, 0.00332607, 0.0106486, 0.00473193, 0, 0.00403549, 0, -0.00104048, -0.0012422, 0, 0, -0.00278746, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
78{0, -0.0327654, 0.00279954, -0.00156276, -0.00392664, -0.00105142, 0.00949674, 0.00992266, 0.00251736, 0.0152795, 0.00306799, -0.00238284, 0.0067777, -0.00157091, 0.00526856, 0.0264038, 0.00147724, 0.00888431, 0.026411, 0.0198684, -0.00143129, 0.00316219, 0.00514419, 0, 0.00560927, 0.00420765, -0.00145472, 0, 0.00956539, 0.00241543, 0.00470629, 0.0048465, 0.0144799, 0, 0.00570209, 0.00418751, -0.00047658, -0.00110488, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
79{0, -0.00753227, 0.00397765, 0.0123975, 0.000839958, -0.0042778, -0.00294332, -0.00318388, -0.00243692, 0.00213064, -0.000518138, 0.00777234, 0.00394029, 0, 0.00305236, 0.00821964, -0.000475594, 0.0129634, 0.00775832, 0.0247704, 0.00518487, 0.00437631, 0.00999549, 0.0040436, 0.00339696, 0.00587269, 0.0123544, 0, 0.0113253, 0.00504108, 0.00467702, 0.00592753, 0.00563897, 0.0149573, 0.00464732, 0.00463117, 0, 0.00761649, 0.0151458, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
80{0, -0.00191676, -0.00423107, 0.00475804, -0.00160728, 0.0088541, 0.00857197, -0.000837, -0.00239001, -0.00049013, -0.000551548, -0.000954402, 0.00374396, -0.00129406, 0.0189372, -0.00141928, -0.000367202, 0.00547228, 0.0088018, -0.00043201, 0, 0.0142185, 0.0103194, -0.000443245, 0.00500681, 0.00566802, 0.0117745, 0.0128396, -0.00146617, 0.00538462, 0.0127682, 0.00472659, 0.00756395, 0.0138696, -0.000561553, -0.00167291, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
81{0, -0.00437213, -0.00102453, 0.00473396, -0.0024484, -0.00156082, -0.000590814, -0.00124694, -0.00349041, 0.00456603, -0.000421245, 0, 0.00621102, 0.00749433, -0.000960832, 0.00314721, 0.00359019, 0.00284196, 0.0137168, 0.00562268, -0.00160554, 0, 0, -0.000470884, 0.00336691, 0.00459624, -0.00111263, 0.0050625, 0.0039622, 0.00446919, 0, 0, -0.000451104, 0.00630724, 0.0259699, -0.0016493, -0.00143629, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
82{0, 0.0116541, -0.000491055, 0, 0.0108019, 0.00891202, 0.0205461, -0.000482761, -0.00115522, -0.00172886, 0.00237479, -0.000464065, 0.000745853, 0.0113383, -0.00227308, 0, -0.00154217, -0.000569069, 0.00912503, -0.000523669, 0, 0.0054163, -0.00129641, -0.00213726, 0.0056876, 0.00560903, 0, -0.000929603, 0, 0, 0.00435316, 0, -0.000487678, 0.0115152, 0.00722282, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
83{0, 0.0101944, -0.000458053, -0.00136645, 0.00293929, -0.00104558, -0.00138759, -0.000995768, 0.0041222, -0.000817748, 0.00406901, 0, 0.00431831, -0.000359942, -0.000918891, 0.0039494, 0.00456461, -0.000445733, -0.000497338, 0, 0.00431183, 0.00336932, 0.00508663, 0, -0.000925412, 0.00551509, 0.00448889, -0.000465728, 0.00337366, 0.00468909, 0, -0.00137258, 0.00603869, 0.00806585, 0, -0.00580589, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
84{0, -0.00237226, 0.010259, -0.000895674, 0.00411476, 0.00512026, -0.00123773, 0.00821976, 0.00269015, 0.00838505, -0.000443918, 0.00770911, -0.000451037, -0.000919729, 0.00376047, 0, -0.00128212, 0.00577606, 0.00402539, -0.00127397, 0.00444049, 0, -0.000432523, -0.00189025, -0.000976431, 0, 0.00498299, 0, 0.00437452, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
85{0, 0, 0.00537186, 0.00395098, 0.0117299, 0.00799079, -0.00188147, 0.00901231, -0.00115061, -0.0022369, 0.00940017, -0.000452578, 0.00529762, -0.000898868, -0.00113458, 0.00950784, 0, -0.00189087, -0.000408221, 0, -0.00091378, 0.00349189, 0, 0.00454828, 0, -0.000567589, 0, 0, -0.000869253, 0, 0.00640011, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
86{0, 0, 0.00875028, 0.00507708, 0, -0.00136453, 0.00372118, 0.00202545, 0, 0, -0.00103392, -0.00142976, 0, 0.00364562, 0.0036454, 0.00372184, -0.000416832, -0.00122799, -0.00117455, -0.000848612, 0, 0.00619048, 0, 0, 0, 0.00300899, 0, 0.00439737, 0, -0.000452727, 0, 0, 0.0223262, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
87{0, 0.00813083, 0.00984675, 0.00418989, 0.00419656, -0.000805791, -0.00125829, 0.00409336, 0.00459685, 0.00424519, -0.000454501, 0.00829167, 0, 0.00525152, -0.000449865, 0.00472471, -0.000496281, -0.00134592, 0, 0.013176, 0.00972065, 0, 0.00337632, 0.00312468, 0, 0, 0.00465113, 0.00523123, 0.0096822, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
88{0, 0.0078588, 0.00438392, 0.00364588, 0.00309135, -0.00206823, -0.000886876, 0.00303936, 0.00295383, -0.00250436, -0.00079063, 0.00529198, 0, 0, 0, -0.00123548, 0.00529033, -0.00133633, 0.0055804, -0.0005559, 0, 0, -0.000420998, 0.0047264, 0.00372627, 0.0050624, 0, 0.00417586, 0.00579403, 0.00922432, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
89{0, 0.0134527, -0.000487816, -0.000407827, -0.00164327, 0.00714982, 0.00834671, 0.0043254, -0.00150873, 0.00909819, -0.00107058, -0.000881873, -0.000906969, -0.000801329, -0.00141174, 0.00361808, 0, 0.00884207, -0.000526459, 0.00422741, 0.00892857, 0.00334632, -0.00141564, 0.00534427, -0.0007912, 0.00382144, 0, 0.00361667, 0.00735391, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
90{0, -0.00112453, -0.000921719, 0.00348905, 0.00986824, 0.00323772, 0.00835248, 0.00388933, 0.00499419, -0.00153203, 0.00456795, 0, 0, 0.00465839, 0, 0.00372416, 0.00465732, 0.00646166, 0, 0.00444293, 0.00290251, 0, 0, 0.00934281, 0.0148351, -0.000809891, 0.0100435, 0, -0.00301733, -0.00375312, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
91{0, 0.0106236, 0.00129889, 0.00448691, 0.00445295, 0.00508553, -0.00138142, 0, -0.000445329, 0.00403718, -0.000866474, 0.00353158, 0, -0.0014014, 0.00440745, -0.000460064, 0.00324248, 0.0143277, 0, -0.00102167, 0.00481517, 0.00437973, 0.0138674, 0.0263436, 0.00353529, 0.00439296, 0.008448, 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},
92{0, -0.000615042, 0.00934905, 0.00757931, -0.00193735, 0, 0, 0, -0.00162904, -0.00201206, 0.00316532, -0.000758207, -0.000460719, 0.00393524, 0.00274551, -0.00122072, -0.00328414, 0.00975402, 0.00479703, 0.00520447, -0.0012812, 0.003169, 0.00331205, 0.0146124, 0.0111438, 0.00451904, 0, 0.075641, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
93{0, -0.000516182, -0.000893867, 0.0121081, -0.0018224, -0.000488139, 0.00268232, 0, 0.00402894, -0.000402761, 0, 0, 0, -0.000405727, 0.00430231, 0.00434484, 0.00410179, -0.000975305, 0.0160911, 0.00498146, 0.0173023, 0.00533898, 0.0278033, -0.000951397, 0.0128432, 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, 0},
94{0, -0.00140846, 0.00422503, -0.00163652, 0.00443282, 0, -0.000351905, 0.00466989, 0.00342072, 0.00303109, 0, -0.000763762, 0, -0.000770398, 0.0202359, 0.00460753, 0.014308, 0.00482234, -0.00179807, 0.0263769, 0.00670195, 0.00955175, 0.0202865, 0.0155921, 0.0189799, 0.019989, 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},
95{0, 0, 0, -0.00132663, 0, -0.000857411, 0, -0.000393507, 0, -0.000818886, 0.00722913, -0.000422264, 0, -0.00124678, 0.00794475, 0.0121623, 0.00420758, 0.00262983, 0.00718098, -0.000505903, 0.00895671, 0.0246874, 0.0180739, 0.00844876, 0.0230061, 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, 0},
96{0, 0, -0.00043891, 0.0057952, -0.00153878, -0.00204123, -0.00274566, -0.00166583, 0.00470958, -0.000736231, -0.00111919, -0.00081959, -0.00271149, 0.00407672, -0.000947025, 0.00436545, 0.0158883, 0.0121416, 0.0202014, 0.0141576, 0.00915484, 0.00371495, -0.000838394, 0.0185427, -0.0330275, 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, 0},
97{0, -0.000569017, -0.00187911, -0.0013152, -0.00205846, -0.00174749, -0.000859905, -0.00139486, 0.00919869, -0.00162995, -0.0016275, -0.00136203, 0.00943644, 0.00406865, 0.00304173, 0.00802036, 0.0223577, 0.0243796, 0.0195699, 0.00451811, 0.00637247, 0.00861249, -0.00157236, 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, 0, 0, 0},
98{0, 0.00676314, 0, 0, -0.00151701, -0.00167819, -0.00164177, 0, 0, -0.00126364, 0.00209941, 0.00780059, 0.00523671, 0.00816129, 0.00920645, 0.00914593, 0.0480522, 0.0182587, 0.0174034, 0.0129708, 0, 0.0709505, 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, 0, 0, 0, 0},
99{0, -0.00132961, -0.00105272, -0.000493242, -0.000998982, -0.000963302, -0.000403533, 0.00321011, 0.00259871, 0.00165038, 0.00945111, -0.000398643, 0.0215238, 0.00245585, 0.0160536, 0.0330436, 0.038898, 0.0155603, 0.0113935, 0.00682208, 0.0266945, 0.0601852, 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, 0, 0, 0, 0},
100{0, 0, -0.00219801, -0.00100347, -0.00112403, 0.00820631, 0, 0.000776552, 0, -0.00309955, 0.00344142, 0.0165217, 0.00824275, 0.00318922, 0.0167143, 0.0183902, 0.0209505, 0.0148305, 0.026167, 0.0465092, 0.0393338, 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, 0, 0, 0, 0, 0},
101{0, -0.00192963, 0.00404023, 0.00429177, 0.00361765, -0.00150455, 0.00368992, 0.00358212, 0.00743376, -0.000420885, 0.0140216, 0.00868896, 0.00296287, 0.0174345, 0.0200781, 0.0511327, 0.03534, 0.0113812, 0.0331168, -0.00423133, 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, 0, 0, 0, 0, 0, 0},
102{0, -0.0017667, 0.00266838, 0.00821613, 0.00252852, -0.000446041, 0.00859937, 0.00201737, 0.00433766, 0.0037743, 0.0244391, 0.0120824, 0.0348515, 0.0460445, 0.0198451, 0.0172826, 0.0730303, 0.125063, 0.0233494, 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, 0, 0, 0, 0, 0, 0, 0},
103{0, -0.00103876, 0.00270405, 0.00510568, 0.0100982, 0.0168746, 0.0105157, 0.00277159, 0.00682157, 0.0246823, 0.00884256, 0.00534351, 0.0465627, 0.0387692, 0.071591, 0.0284448, 0.0953723, 0.0658472, 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, 0, 0, 0, 0, 0, 0, 0, 0},
104{0, 0, 0.0112434, -0.000430343, 0.00425268, -0.00161899, 0.0131034, 0.0092693, 0.00747454, 0.00379318, 0.0193525, 0.0450718, 0.0749774, 0.045105, 0.0726199, 0.0665514, 0.133822, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0},
105{0, -0.0074107, 0.00855827, 0.00373376, 0.00884158, 0.0110644, 0.00795127, 0.0206009, 0.0167017, 0.0175845, 0.0577236, 0.0370533, 0.0527817, 0.0650344, 0.10323, 0.224505, 0.0510321, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0},
106{0, 0, 0.0179825, 0.0167169, 0.00888093, -0.000476096, 0.0128328, 0.0112438, 0.0405886, 0.0298025, 0.0406211, 0.0588893, 0.0363251, 0.0611834, 0.0628016, 0.192035, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
107{0, 0, 0.0114356, 0.00720186, 0.0170689, 0.017875, 0.0293027, 0.023858, 0.0125951, 0.0270779, 0.0515326, 0.0715181, 0.0579174, 0.10595, 0.0701754, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
108{0, 0, 0.00974123, 0.0251431, 0.0149583, 0.0333583, 0.0190875, 0.0191277, 0.0457674, 0.0370717, 0.0917013, 0.0913667, 0.105546, 0.0307692, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
109{0, 0, -0.00060854, 0.0338357, 0.0416551, 0.0254439, 0.0576751, 0.0666466, 0.0329934, 0.0596387, 0.0719475, 0.0694444, 0.0710059, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
110{0, 0, 0.0141723, 0.0627466, 0.0553281, 0.0745046, 0.0518118, 0.0624949, 0.0354849, 0.0395422, 0.0723982, 0.0326973, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
111{0, 0, 0.0441052, 0.0414434, 0.0610666, 0.0120298, 0.0320438, 0.0192959, 0.0258374, 0.0776923, 0.0257947, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
112{0, 0, 0, 0.0417211, 0.00655413, 0.0419119, 0.0660837, 0.0144257, 0.0221876, 0.0484496, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
113{0, 0, 0, 0.0515969, 0.042548, 0.046319, 0.0160193, 0.0190285, -0.00564746, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
114{0, 0, 0, 0, 0.046252, 0.0763126, 0, 0, 0.25, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
115{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, 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},
116{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, 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}
117};
118
119for(int i=0;i<52;i++){
120 for (int j=0; j<52; j++) {
121 HisPDF[i][j] = HisPDFtmp[i][j];
122 }
123}
124
125}
126
127void EvtDtopipi0pi0Plot::initProbMax(){
128 noProbMax();
129}
130
131void EvtDtopipi0pi0Plot::decay( EvtParticle *p ){
132
133loop:
134 p->initializePhaseSpace(getNDaug(),getDaugs());
135
136 EvtParticle *id1,*id2,*id3;
137 EvtVector4R pd1,pd2,pd3;
138 double xmass13,xmass12, xmass23;
139
140 id1 =p->getDaug(0);
141 id2 =p->getDaug(1);
142 id3 =p->getDaug(2);
143
144 pd1 =id1->getP4Lab();
145 pd2 =id2->getP4Lab();
146 pd3 =id3->getP4Lab();
147
148// xmass12=(pd1+pd2).mass()*(pd1+pd2).mass(); // M_ksopi
149 xmass13=(pd1+pd3).mass()*(pd1+pd3).mass(); // M_ksow
150 xmass23=(pd2+pd3).mass()*(pd2+pd3).mass(); // M_piw
151
152 int xbin = FindXBin(xmass23);
153 int ybin = FindYBin(xmass13);
154 double xratio12=HisPDF[xbin][ybin]/avm1;
155
156 if(xratio12 <=0) goto loop;
157
158 double rd12=EvtRandom::Flat(0.0, 1.0);
159 if(rd12>xratio12) goto loop;
160
161 return ;
162}
163
164int EvtDtopipi0pi0Plot::FindXBin(double mass2) {
165 if (mass2 < Xmin) { return 0; }
166 else if (mass2>=Xmax) { return 51; }
167 else { return int((mass2-Xmin)/Xwid)+1; }
168}
169
170int EvtDtopipi0pi0Plot::FindYBin(double mass2) {
171 if (mass2 < Ymin) { return 0; }
172 else if (mass2>=Ymax) { return 51; }
173 else { return int((mass2-Ymin)/Ywid)+1; }
174}
mass
double mass
Definition: CosmicGenerator.cxx:138
EvtDecayBase::noProbMax
void noProbMax()
Definition: EvtDecayBase.cc:304
EvtDecayBase::getParentId
EvtId getParentId()
Definition: EvtDecayBase.hh:60
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition: EvtDecayBase.cc:482
EvtDtopipi0pi0Plot::FindXBin
int FindXBin(double mass2)
Definition: EvtDtopipi0pi0Plot.cc:164
EvtDtopipi0pi0Plot::decay
void decay(EvtParticle *p)
Definition: EvtDtopipi0pi0Plot.cc:131
EvtDtopipi0pi0Plot::FindYBin
int FindYBin(double mass2)
Definition: EvtDtopipi0pi0Plot.cc:170
EvtDtopipi0pi0Plot::clone
EvtDecayBase * clone()
Definition: EvtDtopipi0pi0Plot.cc:42
EvtDtopipi0pi0Plot::getName
void getName(std::string &name)
Definition: EvtDtopipi0pi0Plot.cc:38
EvtPDL::getSpinType
static EvtSpinType::spintype getSpinType(EvtId i)
Definition: EvtPDL.hh:61
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
EvtParticle::getP4Lab
EvtVector4R getP4Lab()
Definition: EvtParticle.cc:684
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:84
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1070
EvtRandom::Flat
static double Flat()
Definition: EvtRandom.cc:74