EvtDtopipipiPlot.cc

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/BesCopyright
//      Copyright (A) 2006      Ping Rong-Gang @IHEP
//
// Module: EvtDtopipipiPlot.cc
//
// Modification history:
//
//    Liaoyuan Dong    August, 2022       Module created
//
//------------------------------------------------------------------------
#include "EvtGenBase/EvtPatches.hh"
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtGenKine.hh"
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtReport.hh"
#include "EvtGenModels/EvtDtopipipiPlot.hh"
#include "EvtGenBase/EvtRandom.hh"
#include <stdlib.h>
#include <string>
 
#include "TFile.h"  
#include "TApplication.h"
#include "TROOT.h"
#include "TH1.h"
#include "TAxis.h"
#include "TH2.h"
using std::endl;
 
EvtDtopipipiPlot::~EvtDtopipipiPlot() {}
 
void EvtDtopipipiPlot::getName(std::string& model_name){
        model_name="DtopipipiPlot";     
}
 
EvtDecayBase* EvtDtopipipiPlot::clone(){
        return new EvtDtopipipiPlot;
}
 
void EvtDtopipipiPlot::init(){
 
    checkNArg(0);
 
    bool idN  = getDaugs()[0]==EvtPDL::getId(std::string("pi+"))||getDaugs()[0]==EvtPDL::getId(std::string("pi-"));
    bool idKs = getDaugs()[1]==EvtPDL::getId(std::string("pi+"))||getDaugs()[1]==EvtPDL::getId(std::string("pi-"));
    bool idPi = getDaugs()[2]==EvtPDL::getId(std::string("pi-"))||getDaugs()[2]==EvtPDL::getId(std::string("pi+"));
    if(!(idN && idKs && idPi ) ){std::cout<<"EvtDtopipipiPlot: the daughter sequence should be pi+ pi+ pi-"<<std::endl;abort();}
    EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId());
 
Xmin = 0.072919;
Xmax = 1.74;
Xwid = 0.0333416;
Ymin = 0.2;
Ymax = 3.1229;
Ywid = 0.058458;
avm1  = 0.2;
double HisPDFtmp[52][52] =
{
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00150066, 0.0155324, 0.00978731, 0, 0, 0, -0.000862562, -0.00027163, 0.00905843, 0.00159547, 0.00674777, 0.00686001, 0.00893704, 0.0263911, 0.0177556, 0.0262276, 0.0163004, 0.00810661, 0.00266422, 0.0165903, 0.0130485, 0.0090461, 0.0138679, 0.00678281, 0.0142944, 0.0124154, 0.00501558, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.025412, 0.018101, 0.0105117, -0.000252347, 0.000857337, 0.00371809, 0.0022962, 0.00552758, -0.000193383, 0.0016975, 0.0020368, 0.00152886, 0.00781011, 0.00445376, 0.00129068, -0.000225775, 0.0105552, 0.0142455, 0.00369242, 0.014687, 0.00895646, 0.010251, 0.0174155, 0.0114122, 0.0157341, 0.00530172, 0.00958369, 0.00871203, 0.0110746, 0.00846772, 0.004005, 0.00648727, 0.00412489, 0.00510047, 0.0077632, 0.00452142, 0.00439873, -0.00254842, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.2, 0.0187074, 0.0151469, 0.0150374, 0.0149382, 0.00852575, 0.0117726, 0.00292745, 0.00594418, 0.00210557, 0.00437641, 0.0127742, 0, 0.00188013, -0.000209292, -0.000645717, 0.00506688, 0.00161251, 0.00951329, 0.0113891, 0.0154694, 0.011311, 0.0102122, -0.0011918, 0.0115697, 0.0033551, 0.00681684, 0.0142289, 0.00486508, 0.0103071, 0.00585104, 0.00654329, 0.00773042, 0.0249313, 0.00100089, 0.00687448, 0.00108067, 0.00575115, 0.00836541, 0.00610641, 0.00498157, 0.0101283, -0.00449353, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0357639, 0.0163919, 0.0267385, 0.0255677, 0.0161029, 0.0172624, 0.00385577, 0.00721942, 0.00383896, -0.000469486, 0.00519038, 0.00132785, 0.00298115, 0.00308214, 0.0060887, 0.00291204, 0.00383603, 0.0067211, 0.00317584, 0.00565859, 0.00959031, 0.00772591, 0.0165113, 0.00941265, 0.0105124, 0.0193347, 0.012998, 0.00561977, 0.00327781, 0.00355879, 0.0113581, 0.0114383, 0.0115121, 0.00388694, 0.00876026, 0.0113876, 0.00312005, 0.0143682, 0.0062588, 0.000176478, 0.0138092, 0.011377, 0.00383445, 0.0129828, 0.0272171, 0, 0, 0},
{0, 0, 0, 0, 0.0112905, 0.0116124, 0.0213804, 0.0328621, 0.0183141, 0.0246236, 0.00394498, 0.00887383, 0.00415104, -0.000519089, 0.00633283, 0.00132945, 0.00350115, 0.00362386, 0.00763285, 0.00291064, 0.00451877, 0.00697439, 0.00438741, 0.00652817, 0.0103673, 0.00837873, 0.022633, 0.00731777, 0.01128, 0.0180281, 0.0147394, 0.00613579, 0.00345585, 0.00418337, 0.0113416, 0.0148126, 0.011185, 0.00484846, 0.0102792, 0.0128383, 0.00288329, 0.0143363, 0.00623997, 0.000258956, 0.0151972, 0.0107407, 0.00460596, 0.0105749, 0.0123443, 0, 0, 0},
{0, 0, 0, -0.00280361, 0.00599815, 0.0124128, 0.026802, 0.0277765, 0.0203783, 0.0152604, 0.00563037, 0.00539587, 0.00742196, 0.00125699, 0.00518185, 0.00362133, 0.00502016, 0.00452122, 0.00639375, 0.00156148, 0.00331946, 0.0037184, 0.0030193, 0.00504048, 0.00936035, 0.00428231, 0.0157025, 0.00756299, 0.00826151, 0.0135952, 0.0104938, 0.00635791, 0.00549246, 0.00312205, 0.00749632, 0.0104138, 0.00878021, 0.00360666, 0.00967957, 0.00906651, 0.00563426, 0.0105905, 0.00818554, 0.00765972, 0.0114866, 0.00864717, 0.00973636, 0.00834703, 0.0127143, 0, 0, 0},
{0, 0, 0, 0.00690654, 0.0201675, 0.0178712, 0.0108822, 0.0112214, 0.0164782, 0.00218226, 0.00930132, 0.0135348, 0.00754681, 0.0178071, 0.00760194, 0.0131149, 0.00212195, 0.00508351, 0.00224026, 0.00228926, 0.00408006, 0.00238163, 0.00168934, 0.00388352, 0.00723827, 0.00615329, 0.0111326, 0.00904487, 0.00498642, -0.000847188, 0.00170438, 0.0126457, 0.0010792, 0.00472998, 0.00522703, 0.0115935, 0.0157492, 0.00747376, 0.00691549, 0.00766887, 0.00714786, 0.00519651, 0.00828873, 0.0199414, 0.00664845, 0.0152771, 0.0136218, 0.0189482, 0.0189, 0, 0, 0},
{0, 0, 0, 0.00402643, 0.0147216, 0.01715, 0.0120051, 0.0136828, 0.0133204, 0.00254233, 0.0104003, 0.0121178, 0.00545226, 0.0217436, 0.00735621, 0.0117699, 0.00172737, 0.00469086, 0.00251054, 0.00212634, 0.00446061, 0.00189206, 0.00206932, 0.00367898, 0.00638883, 0.00674313, 0.0114174, 0.0080861, 0.00369373, -0.000753843, 0.00182248, 0.0106247, 0.0010497, 0.00362544, 0.00541468, 0.0136716, 0.0148882, 0.00695346, 0.0065737, 0.00548947, 0.00620907, 0.00451994, 0.00673862, 0.0194927, 0.00707029, 0.0124406, 0.0129697, 0.0165352, 0.0247423, 0, 0, 0},
{0, 0, 0, 0, -0.000274699, 0.00285875, 0.00658438, 0.0104937, 0.0193084, 0.00435312, 0.002284, -0.00094086, 0.00681717, 0.00351616, 0.00649049, 0.00615866, -0.000691215, 0.00167221, 0.00374339, -0.000538108, 0.00450478, 0.00609494, 0.00136808, -0.00170308, 0.00505613, 0.00796647, -0.000908463, 0.00644384, 0.0009644, 0.00801776, 0.0138419, 0.0108086, 0.00331586, 0.0193975, -0.00128345, 0.00831479, 0.00567967, 0.00248256, 0.00798811, 0.0113088, 0.00398216, 0.0146145, 0.000559689, 0.0136526, 0.00900583, 0.00765891, 0.0113426, 0, 0, 0, 0, 0},
{0, 0, 0, 0, -0.000837011, 0.00783656, 0.0103523, 0.0127591, 0.014678, 0.0184914, 0.00536743, 0.00480503, 0.00397989, 0.0141369, 0.00160304, 0.00385771, 0.0148946, 0.00584356, 0.00591229, -0.000800461, 0.00142273, 0.000380504, 0.00107338, 0.00532479, 0.0055054, 0.000956682, -0.00130562, 0.0105865, 0.0106789, 0.00217481, 0.00547282, -0.000182769, -0.000466737, 0.0040426, 0.0077561, 0.0121245, 0.00192074, 0.00713558, 0.00700134, 0.00474334, 0.00742682, 0.00491013, 0.00400981, 0.00431014, 0.00592222, 0.00413596, 0.00991255, -0.000899015, 0.0260504, 0, 0, 0},
{0, 0, 0, 0, 0, 0.0122997, 0.0169225, 0.0184353, 0.0164037, 0.0117505, 0.00480769, 0.0126396, 0.00399975, 0.00669493, 0.0118789, 0.00342785, 0.00141258, 0.00272821, -0.000557028, 0.00156787, 0.00605785, 0.00222844, 0.000804709, 0.00897894, 0.00733436, 0.00605425, -0.00170722, 0.00645383, -0.0024872, 0.00638874, 0.000761802, 0.00587418, 0.00361459, 0.00200659, 0.0038677, 0.0114339, 0.00660603, 0.00888449, 0.0164075, 0.0104761, 0.0036128, 0.00721415, 0.00894716, 0.0200528, 0.0127048, 0.0114264, 0.0298287, 0.0122927, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.00399151, 0.00615773, 0.0177878, 0.0189579, 0.010571, 0.0172364, -0.000927745, 0.00426211, 0.0109308, 0.00434394, 0.00467184, 0.0102868, 0.00455592, 0.00346033, 0.0011815, 0.00623132, 0.00405261, 0.000927003, 0.00308591, 0.00160743, 0.000821355, 0.0018395, 0.00390381, 0.00580204, -0.00146727, 0.00466621, 0.00100443, 0.00128368, 0.00196323, 0.00562257, 0.00326293, 0.00187777, 0.00104736, 0.00881352, 0.0143722, 0.00304789, 0.0209486, 0.00508492, 0.00432553, 0.0032542, 0.0253143, 0.0332302, 0.0152654, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.0202926, 0.0186309, 0.0104131, 0.00690488, 0.00236005, 0.00434038, 0.0142355, 0.0066333, 0.0109639, 0.0203775, 0.00361271, 0.0162799, 0.00223694, 0.00560054, 0.00378602, 0.00436651, 0.00919364, 0.00832168, 0.00484681, 0.00520677, 0.0132129, 0.00412628, 0.00188165, 0.00560339, 0.00905057, 0.00380511, 0.0101017, -0.00105115, -0.00106044, 0.0177763, -0.000432789, 0.0078263, 0.00926294, 0.00649598, 0.0165854, 0.013363, 0.00388318, 0.027023, 0.0204139, 0.0245485, 0.0147441, 0.00703364, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.0106125, 0.0188372, 0.0146971, 0.00478025, 0.00663613, 0, 0.00601951, 0.0113968, 0.0259423, 0.00674669, 0.00506002, 0.00193793, 0.00393283, 0.00247295, 0.00721892, 0.00038463, 0.00737362, 0.00384348, 0.00494624, 0.00129381, 0.0121777, 0.00153305, 0, 0.00309061, 0.00096083, 0.00635241, 0.00437105, 0.00271235, 0.00323285, 0.0152463, 0.00790541, 0.00197237, 0.00842918, 0.00351487, 0.00691746, 0.0108789, 0.0061616, 0.0240904, 0.0129392, 0.0129275, -0.000830457, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0.0247863, 0.0176282, 0.0175289, 0.00457808, 0.00992452, 0.0115541, 0.0135212, 0.00994434, 0.00713981, 0.0111169, 0.0103148, 0.00416804, 0.00410405, 0.00111858, 0.00873464, 0.00331374, 0.0106796, 0.00201416, 0.00122554, 0.00385263, 0.00661668, 0.00463202, -0.000606177, -0.00120285, 0.00754939, 0.00930348, 0.00377447, 0.00465278, 0.00101872, 0.00231525, 0.0112804, 0.00641832, 0.00746406, 0.00930169, 0.0118229, 0.00592735, 0.0206328, 0.0224974, 0.0115758, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.0136459, 0.00712392, 0.0138838, 0.00159252, -0.000765476, 0.010438, 0.0126405, 0.0152418, 0.0034147, 0.00790645, 0.00381993, 0.00428926, 0.00622248, 0.00742261, 0.00232875, 0.00613106, 0.00168197, 0.00352778, 0.00934557, 0.00150865, 0.00118901, 0.00179582, 0.00168158, 0.00214762, 0.00615021, 0.00110451, -0.000977385, -0.000972892, 0.00159201, -0.00237044, 0.00777439, 0.00259004, 0.00859064, 0.00659134, 0.0016455, 0.00686578, 0.00407572, 0.0201268, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.0190655, 0.0193906, 0.0298653, 0.0156586, 0.00188532, 0.0140794, 0.00497449, 0.00465134, 0.00547378, -0.00155504, 0.0112145, 0.00138882, 0.0100356, 0.00871306, 0.00347013, 0.000565267, -0.00168501, 0.00220346, 0.00541041, -0.000871444, 0.00192547, -0.000892462, 0.0015037, 0.000526234, 0.00204403, 0.00432735, 0.00207964, -0.00125049, 0.00368116, 0.00582942, 0.00563678, -2.20906e-05, -0.000726604, 0.000425052, 0.00823579, 0.00334216, 0.0107216, 0.0110852, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.0102941, 0.0227606, 0.00958853, 0.0112508, 0.00437405, 0.019164, 0.0107376, 0.0073221, 0.0124903, 0.0161441, 0.000924579, 0.00352839, 0.0065253, 0.000840571, 0.00522096, 0.00320728, 0.00169485, -0.0014829, -0.000385922, 0.00380332, 0.00130042, 0.003991, -0.00152388, -0.000335786, 0.00115229, 0.00738576, -0.000475654, 0.00231245, -0.00128109, 0.00129967, 0.00403285, -0.00104735, 0.00164014, 0.0088376, 0.00717689, 0.000210979, -0.000834027, -0.00523533, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00906013, 0.00850382, 0.00337777, 0.00474848, 0.00501928, 0.00262243, 0.00361578, 0.00412986, 0.00828773, 0.00317526, 0.000418765, 0.00571304, 0.0039906, 0.00398882, 0.00315977, 0.00154839, 0.0053812, -0.00215534, 0.00814652, 0.00199931, 0.00156904, 0.00225004, -0.000790613, 0.00182555, 0.00237361, 0, 0.00425822, -0.00124752, -0.00125848, -0.000457286, 0.00132635, 0.00158466, 0.00165806, 0.00104607, 0.00136616, 0.00464344, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0154465, 0.00895089, 0.00729578, 0.00496513, 0.0110113, 0.00404945, 0.00868326, 0.00803355, 0.00121479, 0.000957591, -0.00140863, 0.0111029, 0.00542703, 0.00315061, 0.0142478, 0.00984325, 0.00326606, 0.00141415, -0.000485355, 0.00169373, 0.00175355, -0.000214901, 0.00382858, -0.000839424, -0.00041114, 0.00122857, -0.000428686, -0.000231425, -0.000484801, -0.000677427, -0.00107819, 0.00398122, 0.00389996, 0.00160107, -0.000596897, -0.00473512, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00557644, 0.0117601, -2.014e-05, 0.00532312, 0.00164892, 0.00123443, 0.000548319, 0.00379273, 0.000819126, 0.00315643, -0.00112786, 0.00447013, 0.0039469, -0.000698325, 0.00637213, -0.000918547, -0.000877958, -0.000192955, 0.00386532, 0.00108823, 0, 0.00216934, -0.000612477, -0.000201423, 0, 0.00195611, 0.00260219, -0.000438483, -0.000499695, -0.00107599, -0.00051163, -0.000468805, -0.00136743, -0.000565784, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0113123, 0.00768323, 0.00389946, 0.00981506, 0.00978708, 0.0127364, 0.000633142, 0.00285515, 0.00344821, 0.00683397, -0.000793583, 0.00484872, 0.00140184, 0.00287053, 0.00122665, 0.0019104, 0.00146129, 0.00395712, 0.00180566, -0.00142892, -0.000177317, -0.000627594, 0.00208295, 0, -0.000770852, -0.00199946, -0.000717908, -0.000867859, -0.00107504, -0.00128431, -0.000493855, -0.000265389, -0.0104253, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0101946, 0.00354491, 0.00156002, 0.00514551, 0.00414005, 0.00259554, 0.00681059, 0.018554, 0.00333984, 0.005891, 0.00563848, 0.00506756, 0.00172031, 0.0014715, 0.0049191, -0.00172828, 0.00219327, -0.000623313, 0.000879113, 0.00240745, -0.00113367, -0.000812239, 0.00148996, -0.000525341, 0.00365075, -0.00108374, -0.000710412, 0.00358373, -0.00126742, -0.0012819, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00342704, 0.00486663, 0.00801254, 0.00482659, 0.0067724, 0.000801564, 0.00216852, 0.00120544, 0.00117051, 0.00292921, 0.00293786, 0.00435708, 0.00336549, 0.00184459, 0.00356918, -0.00143222, 0.00271472, -0.000847567, -0.00117673, -0.000591092, -0.000920364, -0.00162241, 0.00213705, -0.000244555, -0.000452254, 0.000957899, -0.000852445, -0.000524743, 0.00174564, 0.00208469, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.01942, 0.0104793, 0.00238674, -0.00028427, 0.00323803, 0.00428694, 0.0180191, 0.00205732, 0.004827, 0.013351, 0.0164632, 0.00661175, 0.00625597, 0.002846, 0.00201106, 0.000985679, -0.000678121, -0.000691974, 0.00196908, 0.000656031, -0.00179116, 0.000615601, 0.00184952, -0.000819999, 0.00320967, -0.00193909, -0.00105235, -0.00205237, 0.00162034, 0.00482908, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0250022, 0.00643193, 0.0145664, 0.0105645, 0.005722, 0.0114656, 0.0105034, 0.00863349, 0.00663601, 0.00534037, 0.0083902, -2.58751e-05, 0.00179951, 0.00120073, -0.000662521, -0.000201754, -0.000705186, 0.0016609, -0.000236756, 0.00330577, 0.00504613, 0.000588122, 0.00174472, -0.000951208, 0.00253816, 0.000456801, 0.0059189, -0.0011968, -0.00077095, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00243762, 0.0101053, 0.0137352, 0.00989112, 0.00419056, 0.0167929, -0.00150345, 0.00391754, 0.00962137, 0.0111504, 0.013032, 0.0091798, 0.00632917, 0.00406742, 0.00171305, 0.0030055, 0.00404168, -0.00197613, 0.00241273, 0.00871242, 0.00473052, 0.00142909, 0.00582945, 0.00698106, 0.00411218, 0.00363546, 0.00151842, -0.000477552, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0137108, 0.0101587, 0.00512209, 0.00581154, 0.0063175, 0.0148431, 0.00784464, 0.0153501, 0.00301844, 0.0187667, 0.007718, 0.00198616, 0.00344126, 0.000856434, 0.006128, 0.00407622, 0.0034297, 0.000401058, -0.000712372, -0.000616223, -0.000855106, 0.00335076, -0.000459257, -0.00100754, -0.000922785, 0.00190541, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0115846, 0.00513737, 0.00809827, 0.00997606, 0.00655055, 0.00941147, 0.0038603, 0.00929101, 0.00884371, 0.0032594, -0.000438409, -0.000647537, 0.00152629, -0.00160462, 0.000699107, 0.00165622, -0.00124744, -0.000977481, 0.0014039, 0.00174867, 0.00401864, 0.00399156, 0.00201647, -0.000601122, -0.00247459, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00105089, 0.0048501, 0.00845776, 0.00199121, 0.00916815, 0.0118089, 0.0101806, 0.0137388, 0.0176601, 0.00328823, 0.00352053, 0.00419714, 0.00546382, 0, 0.00161903, -0.00119964, 0.00191271, -0.000415298, 0.00674503, 0.00439691, -0.000762733, -0.000813587, 0.00370975, 0.00708438, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0145711, 0.0094525, 0.00818695, -0.000665325, 0.0047514, 0.00368851, 0.00805664, 0.00446138, 0.00914161, 0.00622046, 0.00112126, 0.000882917, 0.0016548, 0.00429119, -0.00142295, -2.06638e-05, -0.00157297, -0.000195611, 0.00347694, -0.000951127, 0.00278539, 0.00166737, -0.000765855, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0124235, 0.00144789, 0.00820455, 0.0031008, 0.00870771, 0.00334443, 0.0121713, 0.00201216, 0.00607529, 0.00177875, 0.00140863, 0.00199694, 0.00242735, 0.00139591, 0.00121349, 0, -0.00124997, -0.000215363, 0.00361555, -0.000206134, -0.000716385, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00128829, 0.00342748, 0.00818716, 0.00565558, 0.00871036, 0.0138593, 0.00145614, 0.00132705, 0.0106219, 0.00618985, 0.00187038, -0.00137881, 0.0018989, 0, 0.00185971, 0.00161027, 0, 0, 0.00198384, 0.0021881, -0.00316355, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00148772, 0.00977291, 0.00858047, 0.00392799, 0.00286997, 0.00420178, 0.0063399, 0.00186727, 0.00875324, 0.00169687, 0.00182879, 0.00501749, 0, -0.000424115, -0.000195288, -0.00132863, 0.00201768, -0.000723473, -0.000586035, -0.00211294, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00795127, 0.000761631, 0.0130024, 0.00888905, 0.00533266, 0.00165759, -0.000681874, 0.00613109, 0.00421136, 0.0026058, -0.00100585, -0.000440221, -0.00102298, 0.00181682, 0.00229075, 0.00159399, -0.000890541, 0.0015385, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00176391, 0.0016436, 0.00778176, 0.0019574, 0.00373541, 0.00121501, 0.00194797, 0.00156388, 0, 0.00205196, 0, 0.00194588, -0.000447248, 0.00164912, 0.00378861, -0.000967688, -0.000535921, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.000414658, 0.00158378, 0.00792649, 0.00747993, 0.00355962, -0.000222224, -0.000883339, -0.000675862, 0.00445083, -0.00039295, -0.000978694, 0.00126864, 0.00130999, -0.000655967, 0, -0.000605564, -0.00786369, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00534462, 0.00337856, 0.00447068, 0.00202555, 0.00262315, -0.000791716, 0.0017379, -0.000440789, -0.000377595, 0.00199316, -0.000205376, 0.00113595, 0.000595314, 0, -0.00101937, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00194879, 0.0043191, 0.00349993, 0.00656448, -0.00022874, 0.00490897, 0.00173958, 0.00211379, -0.000919604, 0.00122627, 0.00192831, 0.000924648, -0.000533183, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00550371, 0.00377472, 0.00180269, 0.0031602, -0.000686968, -0.000478745, -0.000771646, -0.00125527, 0.00573005, 0.00298008, 0.00219993, 0.00456167, 0.00815102, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00171234, 0.00450387, -0.000745589, 0.000928437, 0.00198709, -0.00020945, -0.000248579, 0.00112635, 0.00215338, 0.0066882, 0.0106831, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00106871, -0.000222508, 0.00205684, -0.00038869, 0.00157722, -0.000187058, -0.000198706, 0.00876609, 0.00483817, 0.00529853, 0.0143598, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00220188, -0.00065229, 0.00168199, 0.00391551, 0, 0.00777458, 0.00353025, 0, -0.000506828, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00189452, 0.00193445, 0.00152572, -0.000211959, 0.00383985, 0.00874946, 0.00457722, 0.0299003, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.000708505, -0.0011835, 0.00362797, -0.00112053, 0.0057435, 0.00834591, 0.0112548, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00433094, -0.00162684, 0.00649451, 0.00462054, 0.0141132, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0139919, -0.00152095, 0.0163195, 0.010359, 0.0186132, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0147059, 0.00529507, 0.0173027, 0.00279661, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00801884, 0.0144967, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00296133, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
};
 
for(int i=0;i<52;i++){
  for (int j=0; j<52; j++) {
     HisPDF[i][j] = HisPDFtmp[i][j];
  }
}
 
}
 
void EvtDtopipipiPlot::initProbMax(){
        noProbMax();
}
 
void EvtDtopipipiPlot::decay( EvtParticle *p ){
 
loop:
        p->initializePhaseSpace(getNDaug(),getDaugs());
 
        EvtParticle *id1,*id2,*id3;
        EvtVector4R pd1,pd2,pd3;
        double xmass13,xmass12, xmass23;
 
        id1 =p->getDaug(0); 
        id2 =p->getDaug(1); 
        id3 =p->getDaug(2);
 
        pd1 =id1->getP4Lab();
        pd2 =id2->getP4Lab();
        pd3 =id3->getP4Lab();
 
//      xmass12=(pd1+pd2).mass()*(pd1+pd2).mass();  // M_ksopi
        xmass13=(pd1+pd3).mass()*(pd1+pd3).mass();  // M_ksow
        xmass23=(pd2+pd3).mass()*(pd2+pd3).mass();  // M_piw
 
                int xbin = FindXBin(xmass23);
                int ybin = FindYBin(xmass13);
                double xratio12=HisPDF[xbin][ybin]/avm1;
 
        if(xratio12 <=0) goto loop;
 
        double rd12=EvtRandom::Flat(0.0, 1.0); 
        if(rd12>xratio12) goto loop; 
 
        return ;
}
 
int EvtDtopipipiPlot::FindXBin(double mass2) {
  if (mass2 < Xmin) { return 0; }
  else if (mass2>=Xmax) { return 51; }
  else { return int((mass2-Xmin)/Xwid)+1; }
}
 
int EvtDtopipipiPlot::FindYBin(double mass2) {
  if (mass2 < Ymin) { return 0; }
  else if (mass2>=Ymax) { return 51; }
  else { return int((mass2-Ymin)/Ywid)+1; }
}