EvtD0topipietaPlot.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: EvtD0topipietaPlot.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/EvtD0topipietaPlot.hh"
#include "EvtGenBase/EvtRandom.hh"
#include <stdlib.h>
#include <string>
 
#include "TFile.h"  
#include "TApplication.h"
#include "TROOT.h"
#include "TH1.h"
#include "TH2.h"
#include "TAxis.h"
using std::endl;
 
EvtD0topipietaPlot::~EvtD0topipietaPlot() {}
 
void EvtD0topipietaPlot::getName(std::string& model_name){
        model_name="D0topipietaPlot";     
}
 
EvtDecayBase* EvtD0topipietaPlot::clone(){
        return new EvtD0topipietaPlot;
}
 
void EvtD0topipietaPlot::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("eta"));
    if(!(idN && idKs && idPi ) ){std::cout<<"EvtD0topipietaPlot: the daughter sequence should be pi+ pi- eta"<<std::endl;abort();}
    EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId());
 
Xmin = 0.075625;
Xmax = 1.7956;
Xwid = 0.0343995;
Ymin = 0.469225;
Ymax = 3.0625;
Ywid = 0.0518655;
avm1  = 64.8599;
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, 10.9175, 8.27593, 2.90859, 7.20965, 3.42651, 9.24549, 2.76522, 1.40342, 3.18176, 3.93846, 3.91516, 20.6344, 3.58422, 4.35362, 5.7482, 1.41577, 1.75627, 2.39247, 3.92711, 0.627232, 2.4743, -0.146978, 0.984221, -1.93913, -0.646377, -0.00270402, -0.577094, 0.707847, -2.08611, 5.12895, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 3.55907, 4.66306, 0.950407, 1.86677, 1.97746, 3.52687, 2.3963, 6.32973, 5.66186, 1.65555, 4.05733, 7.66306, 3.81898, 3.16844, 7.08957, 6.89748, 6.14049, 2.56748, 4.8378, 3.16366, 3.34193, 1.65233, 2.05854, 0.999099, 1.4982, -0.00876633, -0.21596, 0.271792, -0.217222, 4.63075, 0.42579, -0.369247, 0.197116, 0.0860022, -1.8316, -0.308271, -0.0376458, 1.06267, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 13.8247, 1.98493, 6.44637, 2.03897, 0.584747, 2.04085, 4.12907, 6.49531, 1.62324, 2.83629, 2.76695, 4.02672, 2.6792, 3.0508, 10.8597, 3.66585, 21.62, 2.96017, 2.83761, 6.37106, 2.37234, 11.9749, 7.05564, 4.51659, 1.94883, 2.89088, 0.810742, 0.423342, 3.23828, 1.55276, 0, 1.13574, 1.60283, 0.260828, 0.337017, 0.15641, -0.610435, 0.0861614, -2.34386, 0.302828, -2.80513, 0.135101, 2.21786, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 20.737, 2.9774, 9.66956, 3.05846, 0.877121, 3.06127, 6.19361, 9.74297, 2.43486, 4.25443, 4.15042, 6.04008, 4.0188, 4.57619, 16.2896, 5.49877, 32.43, 4.44025, 4.25642, 9.55659, 3.55851, 17.9623, 10.5835, 6.77489, 2.92324, 4.33632, 1.21611, 0.635013, 4.85742, 2.32914, 0, 1.70361, 2.40424, 0.391241, 0.505525, 0.234615, -0.915652, 0.129242, -3.51579, 0.454242, -4.2077, 0.202651, 3.32679, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 41.474, 5.9548, 19.3391, 6.11692, 1.75424, 6.12254, 12.3872, 19.4859, 4.86973, 8.50887, 8.30084, 12.0802, 8.0376, 9.15239, 32.5791, 10.9975, 64.8599, 8.8805, 8.51284, 19.1132, 7.11701, 35.9247, 21.1669, 13.5498, 5.84648, 8.67264, 2.43223, 1.27003, 9.71485, 4.65829, 0, 3.40723, 4.80848, 0.782483, 1.01105, 0.469229, -1.8313, 0.258484, -7.03158, 0.908484, -8.4154, 0.405303, 6.65357, 0, 0, 0},
{0, 0, 0, 2.59031, 0, 2.65679, 10.2401, 8.52322, 10.2129, 15.129, 8.13256, 10.1182, 6.08062, 3.10212, 3.42579, 1.00428, 0.341389, -0.935525, 0.976943, 1.79233, 2.37185, 3.30414, 4.73653, 6.13036, 6.93277, 5.2688, 2.692, 11.2957, 6.961, 7.67561, 1.11344, 0.977715, 1.63171, 3.59676, 9.15318, 1.97157, 3.30907, 8.5026, 1.04297, 0.499947, 1.74369, 5.55907, 11.3386, 4.02254, 11.5053, 1.8794, 0, 5.51071, 0, 0, 0, 0},
{0, 0, 0, 1.72687, 0, 1.77119, 6.82672, 5.68215, 6.80857, 10.086, 5.4217, 6.74549, 4.05374, 2.06808, 2.28386, 0.669523, 0.227593, -0.623683, 0.651295, 1.19488, 1.58123, 2.20276, 3.15769, 4.0869, 4.62185, 3.51253, 1.79467, 7.53044, 4.64067, 5.11707, 0.742291, 0.65181, 1.0878, 2.39784, 6.10212, 1.31438, 2.20604, 5.6684, 0.695313, 0.333298, 1.16246, 3.70604, 7.55907, 2.6817, 7.6702, 1.25293, 0, 3.67381, 0, 0, 0, 0},
{0, 0, 0, 1.72687, 0, 1.77119, 6.82672, 5.68215, 6.80857, 10.086, 5.4217, 6.74549, 4.05374, 2.06808, 2.28386, 0.669523, 0.227593, -0.623683, 0.651295, 1.19488, 1.58123, 2.20276, 3.15769, 4.0869, 4.62185, 3.51253, 1.79467, 7.53044, 4.64067, 5.11707, 0.742291, 0.65181, 1.0878, 2.39784, 6.10212, 1.31438, 2.20604, 5.6684, 0.695313, 0.333298, 1.16246, 3.70604, 7.55907, 2.6817, 7.6702, 1.25293, 0, 3.67381, 0, 0, 0, 0},
{0, 0, 3.28134, 3.70965, 4.56448, 4.80283, 3.15051, 3.05373, 5.2204, 6.29655, 4.93902, 7.67291, 7.9223, 2.49748, -0.22948, 0.11467, 0.815377, 0.685755, 1.04926, 1.0851, 1.2234, 2.55338, 5.75864, 11.5304, 4.85302, 7.74369, 4.46235, 2.47823, 2.31253, 1.19303, 4.10483, 2.68816, 1.96235, 2.28386, 2.05016, 1.18916, 1.58003, 9.40848, 2.18062, 2.85032, 1.86824, 2.79746, 2.10122, 7.56087, 1.99639, 3.16126, 3.71145, 11.6308, 0, 0, 0, 0},
{0, 1.46055, 2.60212, 1.79712, 2.90589, 2.28194, 4.72487, 4.92182, 8.74099, 18.7329, 7.39901, 10.1338, 3.9737, 1.25208, -0.245896, 0.454465, 2.95695, 0.482418, 0.385033, 1.38346, 1.66366, 6.6702, 1.74243, 2.99784, 5.448, 8.4492, 1.84294, 2.08601, 4.04419, 1.82579, 1.2642, 1.90949, 2.31489, 1.89067, 2.31898, 7.70424, 3.17291, 4.39781, 2.18062, 6.4444, 4.45695, 8.25444, 13.6666, 4.00834, 4.12094, 3.56844, 7.06448, 1.06808, 0, 0, 0, 0},
{0, 4.99279, 13.2006, 2.73475, 2.02503, 2.68025, 1.77686, 2.1884, 13.7024, 5.40242, 7.04655, 7.5824, 3.86736, 3.72847, 1.88346, 0.140282, 0.684553, 0.686957, 0.850619, 3.66066, 1.55913, 3.24549, 2.96844, 4.99639, 6.68455, 2.30362, 8.56628, 2.00716, 4.086, 2.08038, 7.77413, 2.4743, 11.3871, 7.31357, 7.49459, 3.67381, 3.38706, 4.22435, 3.96996, 4.70785, 19.5446, 8.68335, 6.15498, 11.9551, 3.10956, 1.7228, 9.62715, 3.41569, 7.27953, 0, 0, 0},
{0, 4.48738, 2.72757, 4.25323, 1.16997, 2.94308, 34.7525, 6.36644, 4.31478, 11.8924, 2.11247, 3.68025, 6.67291, 21.6988, 2.04479, 2.40589, 1.31037, 9.19343, 0.905257, 3.76691, 0.980276, 3.73535, 3.85446, 6.06808, 8.37632, 5.38887, 2.82884, 7.82913, 4.62185, 1.64085, 2.71055, 1.38606, 2.7234, 6.23294, 7.77593, 3.10956, 4.48565, 10.0155, 3.1671, 30.9031, 41.9031, 18.0967, 4.31002, 15.3476, 7.72757, 20.6988, 1.38706, 0, 0, 0, 0, 0},
{0, 1.92471, 9.90668, 4.94083, 12.844, 5.35122, 6.62218, 12.0597, 4.10212, 11.4014, 6.34201, 5.81718, 8.98738, 5.10392, 1.44398, 4.20965, 4.13905, 0.696034, 2.09531, 1.78013, 3.28314, 2.81988, 7.99279, 2.60009, 4.55913, 4.93977, 4.80523, 1.38749, 5.70965, 3.12185, 2.2485, 7.37752, 3.5541, 4.03284, 3.05733, 6.367, 4.59198, 3.15262, 10.6729, 12.1338, 40.9031, 49.4121, 16.1684, 11.8267, 11.3494, 4.72847, 14.4874, 0, 0, 0, 0, 0},
{0, 7.89841, 7.16006, 12.9175, 6.34761, 9.96649, 5.0865, 6.68395, 5.51507, 3.93984, 4.36915, 9.86734, 3.09387, 3.70604, 3.18134, 3.17561, 1.76999, 3.3241, 3.42606, 0, 8.1702, 6.18515, 3.724, 8.09198, 3.3732, 2.94304, 7.9474, 5.32973, 3.90262, 3.71181, 3.5869, 5.78134, 11.2419, 20.4121, 10.3034, 4.73219, 10.9238, 9.25302, 48.6163, 4.46594, 7.22937, 24.8118, 0, 5.55673, 7.02089, 23.0573, 11.491, 0, 0, 0, 0, 0},
{0, 4.65946, 6.51608, 5.11407, 4.80785, 8.49009, 36.7453, 4.79088, 39.5266, 24.335, 16.534, 7.06568, 7.60122, 2.55553, 4.51251, 1.93971, 3.7541, 5.82793, 9.13436, 5.63676, 8.172, 0, 4.18876, 18.138, 7.46145, 17.4587, 3.78391, 4.34348, 2.07884, 5.62094, 5.01522, 13.8853, 16.0573, 4.56943, 15.0896, 7.41073, 8.34107, 6.45797, 13.2795, 10.7939, 63.1934, 48.3404, 6.97746, 14.4229, 17.8853, 8.30463, 0, 0, 0, 0, 0, 0},
{0, 7.04926, 4.97937, 10.4898, 10.5818, 7.5211, 3.83292, 28.0212, 17.0083, 44.3188, 17.6069, 9.77773, 5.83724, 5.70785, 4.61001, 4.06568, 3.41576, 3.47131, 5.87451, 3.94263, 3.67126, 3.44699, 3.71274, 11.0442, 3.74865, 4.69692, 6.46952, 6.19443, 10.6021, 12.7813, 18.2723, 4.41145, 6.9474, 5.43547, 4.99139, 26.0268, 19.4683, 36.7042, 30.5591, 21.0167, 21.62, 16.2455, 16.353, 27.2437, 8.16246, 0, 0, 0, 0, 0, 0, 0},
{0, 5.88887, 8.23054, 4.94467, 3.47531, 13.6069, 14.3679, 7.83436, 6.90217, 12.0627, 20.7705, 3.76163, 25.1181, 2.58134, 6.35002, 7.20424, 3.37004, 3.81718, 2.4982, 8.51969, 29.3476, 12.4253, 7.45519, 12.3763, 11.2115, 7.97937, 12.1159, 7.92651, 8.54358, 4.19784, 15.4587, 5.79784, 6.10392, 17.5197, 14.6607, 16.5173, 12.8243, 23.1601, 15.8154, 18.3709, 26.3046, 39.6362, 13.0502, 38.4874, 14.3889, 0, 0, 0, 0, 0, 0, 0},
{0, 3.29926, 5.65302, 15.7006, 8.37917, 7.65228, 10.5603, 5.46303, 0, 11.6522, 4.07823, 6.7376, 3.59495, 11.4157, 0, 11.7042, 3.41073, 6.88977, 6.34049, 5.60483, 5.29748, 21.0681, 5.2491, 3.7025, 7.05938, 10.1944, 4.56988, 5.81463, 5.19784, 3.7479, 4.40904, 13.4587, 6.2234, 17.6344, 7.52089, 11.1971, 16.7804, 28.8279, 17.9059, 12.2934, 9.29389, 14.3181, 44.7741, 15.9964, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 9.02683, 11.3739, 3.93495, 6.99639, 5.11734, 17.2269, 5.2186, 5.75624, 10.1326, 12.7365, 5.0833, 3.53832, 4.76161, 3.46115, 5.48565, 2.31047, 3.09983, 12.2437, 4.37634, 10.9498, 5.77419, 7.85543, 6.7246, 38.7024, 7.59784, 12.2103, 3.61085, 7.97937, 4.47847, 4.53905, 4.54049, 13.8494, 9.61287, 11.9597, 20.5173, 52.9067, 18.4707, 21.4707, 30.1003, 31.2813, 19.138, 10.9988, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 10.0286, 10.1601, 4.87752, 7.80042, 19.344, 9.52953, 8.63526, 6.33688, 10.1155, 11.8231, 3.37632, 24.4694, 7.24369, 9.24009, 8.42291, 4.70147, 15.5286, 9.71145, 3.1577, 15.1637, 54.8458, 10.8548, 4.52759, 22.4587, 9.97937, 8.61407, 4.68202, 4.16772, 4.53224, 7.85062, 6.78044, 10.1362, 4.58882, 7.51341, 7.9731, 31.6684, 9.69603, 25.7795, 14.043, 14.1637, 15.9964, 34.0645, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 9.672, 3.2817, 29.0501, 9.63556, 20.8404, 16.0214, 10.6899, 7.90138, 4.70204, 16.8799, 15.0196, 3.88826, 3.89603, 10.8853, 6.96055, 5.48242, 4.68099, 5.75268, 15.4719, 8.7025, 9.94049, 11.9794, 9.96145, 8.26432, 14.8136, 9.35002, 41.4085, 3.81778, 15.0967, 12.1666, 6.30223, 5.05887, 5.0817, 9.10212, 10.6021, 4.3723, 6.55338, 9.7464, 13.3154, 6.16988, 10.9247, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 12.5591, 6.79926, 10.4551, 5.22577, 21.8458, 35.6199, 5.79148, 13.3942, 6.92291, 10.8243, 4.9713, 3.29655, 3.72577, 3.93905, 4.79449, 7.78206, 10.0502, 15.2568, 9.37812, 19.7777, 7.82384, 11.0537, 5.96722, 10.9928, 5.95266, 11.4193, 3.90724, 5.13822, 2.50535, 1.28966, 2.78062, 6.62004, 2.92399, 5.35422, 6.75504, 15.7401, 20.206, 7.81718, 9.25804, 4.06628, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 2.03942, 3.18756, 3.52686, 4.21055, 2.67561, 13.5376, 4.71861, 12.258, 10.2051, 3.5161, 5.27473, 11.344, 6.75384, 3.65302, 3.14022, 7.19567, 11.7084, 5.43624, 14.9597, 9.48242, 8.55613, 8.04479, 9.28043, 3.34355, 7.86824, 8.3315, 2.22416, 2.59207, 3.9211, 3.33688, 2.2216, 1.43545, 1.97487, 13.5609, 4.97667, 11.0967, 5.19712, 2.2344, 2.42651, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 7.48738, 2.10114, 14.5591, 5.32913, 16.6797, 8.47881, 6.71861, 5.09623, 13.206, 7.05014, 7.18812, 3.7501, 4.18062, 5.10212, 2.45249, 16.7741, 7.39495, 15.2625, 18.3458, 16.6117, 22.0949, 13.2066, 20.473, 5.76628, 3.88826, 3.93812, 5.00947, 3.88505, 2.32372, 5.34942, 16.7597, 2.27713, 4.5376, 8.05005, 6.04836, 8.99399, 6.28074, 4.14067, 4.92471, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1.82579, 5.54118, 6.04655, 3.01789, 19.0877, 6.45815, 7.90749, 13.525, 9.8912, 7.77503, 4.84852, 1.95697, 6.15525, 7.98738, 2.29328, 8.79745, 29.4013, 8.54655, 8.11914, 6.36557, 4.49639, 4.1989, 2.80523, 8.30103, 2.97488, 3.18515, 2.05847, 1.71938, 3.40318, 3.9438, 8.48738, 6.84581, 2.5376, 2.80336, 6.16006, 3.63796, 2.586, 9.56267, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1.70965, 3.0842, 18.9748, 3.01789, 2.95374, 4.28056, 23.3744, 10.2401, 10.4014, 7.38346, 4.36317, 2.73234, 2.47881, 5.4193, 3.77653, 13.2006, 7.82552, 4.76161, 3.73168, 2.70095, 7.06327, 2.81538, 1.06688, 2.81778, 0.792092, 1.50821, 1.12182, 1.96649, 2.56538, 1.84168, 1.90081, 15.3404, 3.13706, 5.60212, 7.52683, 2.59712, 0.427026, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 5.06688, 1.83941, 3.45785, 1.83555, 5.34047, 2.79601, 2.91575, 6.10677, 7.27473, 3.62542, 11.8959, 3.98918, 3.00857, 3.31538, 6.25444, 1.88616, 5.03796, 7.08717, 14.4766, 2.45995, 1.65814, 3.53314, 2.40498, 3.01548, 3.8333, 1.28025, 1.6523, 2.60122, 0.608787, 1.3369, 3.70604, 10.197, 1.37428, 1.92591, 1.34348, 3.42291, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 1.79569, 2.29388, 11.9874, 3.28025, 3.85242, 7.11824, 3.24549, 1.6541, 4.10986, 7.37985, 2.50774, 1.81309, 7.63556, 3.2523, 5.99459, 3.83409, 4.06868, 3.50834, 2.89603, 5.12004, 3.94173, -0.334083, 3.87451, 3.4175, 2.13796, 2.5833, 3.31628, 4.32471, 1.87004, 5.34942, 4.77773, 0.896755, 0.681503, 2.06808, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 8.20785, 3.27953, 2.5842, 3.05302, 5.68455, 3.0523, 2.37376, 5.4193, 21.0464, 2.80583, 5.53998, 2.71085, 4.04836, 3.7234, 3.997, 6.80403, 5.42579, 9.28074, 4.78044, 4.02815, 2.83409, 4.38346, 2.49549, 2.66186, 1.71037, 2.37872, 5.35122, 3.63676, 0, 4.18278, 1.21197, 3.24189, 1.01015, -0.0385471, 2.99279, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 1.21306, 5.24369, 0.199505, 1.43941, 2.06172, 3.88242, 9.95514, 3.22507, 0, 6.06448, 9.56267, 4.086, 3.91037, 3.65008, 35.5555, 6.88977, 5.82579, 8.95695, 7.63796, 4.74549, 3.36772, 14.3368, 1.72086, 4.13736, 3.88242, 4.55673, 2.133, 8.30343, 1.17621, 4.02029, 0.844842, 3.85302, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0.529435, 1.88977, 0.511092, 2.01073, 1.27119, 1.17501, 4.9713, 2.0389, 4.40769, 8.92291, 8.21025, 3.28365, 6.82579, 4.31181, 4.83093, 4.20339, 6.13887, 14.9928, 5.33033, 25.0537, 3.12842, 12.3494, 2.79639, 5.24549, 8.61407, 7.85302, 19.9892, 5.65232, 2.1541, 2.81808, 1.62651, 3.77773, 2.49459, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0.435923, 4.70244, 1.49579, 1.4986, 1.39963, 2.38977, 1.38171, 9.56448, 10.1159, 7.62185, 5.62902, 9.31898, 24.3494, 9.25447, 16.0203, 9.19263, 5.28314, 6.22847, 6.83941, 3.30744, 2.0421, 2.93457, 7.12004, 17.3136, 6.88242, 6.39712, 4.68242, 2.19252, 3.55483, 1.78374, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 2.78134, 3.06087, 0.207847, 0.328526, -0.128155, 2.64013, 1.71055, 8.70785, 5.31181, 4.7025, 4.192, 8.25519, 13.9235, 23.0304, 10.1828, 11.2688, 8.69676, 11.7464, 4.79415, 0, 12.2401, 7.75338, 11.5627, 12.5878, 5.35338, 5.83941, 2.20339, 4.32973, 1.9743, 5.96055, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, -0.026299, 0.235348, 1.10212, 0.235949, 0.62579, 8.27232, 1.29926, 3.0688, 3.80824, 4.71181, 7.586, 12.1935, 8.46175, 3.6216, 24.6392, 10.3718, 35.6308, 10.6141, 13.1613, 13.491, 4.66051, 7.95266, 6.99567, 7.70965, 7.76251, 2.84746, 5.28224, 2.09011, 5.37872, 5.13977, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 2.9211, -0.0555673, 0.137964, 0.640668, 1.61407, 1.72847, 1.74049, 5.60392, 5.68314, 13.7109, 16.5424, 9.59318, 9.11181, 15.2819, 26.1021, 29.0268, 6.99794, 7.5851, 9.77773, 7.34107, 7.50965, 6.90893, 8.18134, 38.1217, 12.184, 18.1003, 4.46145, 6.21025, 8.84942, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, -0.163096, -0.217762, -0.0615411, 1.06267, 1.99099, 2.73977, 7.01432, 3.92543, 4.48028, 9.28314, 13.6552, 10.2545, 11.3254, 7.50844, 12.6747, 11.1389, 18.3494, 7.45785, 12.8145, 16.5573, 12.2267, 21.5591, 11.5134, 4.77553, 11.7795, 0, 3.11707, 6.85302, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, -0.646076, -0.335737, -0.401484, 0.340404, 0.895531, 16.8422, 8.672, 7.69263, 9.94353, 7.98821, 33.5663, 62.5591, 14.2831, 9.45041, 6.58278, 23.8835, 10.4444, 7.96055, 9.88616, 6.99519, 11.9937, 7.60929, 3.96772, 6.70845, 3.54539, 7.77773, 3.56808, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, -0.764652, 0.02029, -0.299894, 0.753383, 1.21073, 0.701985, 2.91109, 7.7473, 25.172, 8.24321, 8.39067, 7.77553, 22.4719, 10.4831, 16.9462, 5.5581, 10.2607, 8.85158, 10.8721, 17.1577, 7.50893, 7.03869, 7.97757, 10.0304, 3.83033, 4.9211, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, -0.532525, -2.22407, -0.0896287, 0.264316, 6.45514, 10.6344, 7.836, 8.36737, 22.3002, 23.5914, 15.0258, 28.8279, 8.41336, 6.95646, 0, 7.29459, 9.172, 14.4874, 13.7036, 4.62185, 2.85302, 4.06808, 4.38887, 13.5663, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, -0.145175, -0.359632, -0.0215269, 0.780134, 1.86617, 2.95698, 4.65374, 20.9008, 16.5152, 22.9235, 41.706, 12.902, 26.9175, 8.77863, 6.97989, 4.92255, 8.8333, 7.46592, 8.36914, 11.0406, 16.197, 18.2759, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0.281336, -0.00480715, 1.17471, 2.37094, 1.23941, 7.21145, 17.2464, 16.4606, 11.467, 19.4068, 39.2759, 11.7228, 10.7335, 6.85014, 12.203, 5.31129, 10.4432, 1.9976, 3.36772, 5.35212, 12.5591, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, -1.44093, 1.82913, -0.0932131, 3.4468, 6.9223, 6.7834, 12.3541, 14.9259, 21.9534, 20.9955, 7.39694, 12.4713, 15.206, 25.878, 5.62146, 12.5537, 4.48978, 2.90498, 3.4772, 1.85483, 7.27593, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, -0.259512, 1.32853, 2.02089, 6.80283, 5.7834, 0, 21.9238, 26.4946, 19.5672, 10.3743, 32.163, 8.56628, 6.05732, 4.30105, 3.66592, 6.60926, 3.66366, 9.9193, 6.49099, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.307836, 3.87094, 3.25168, 16.6729, 11.7939, 11.3813, 36.3237, 12.3942, 13.0316, 13.1785, 10.2938, 14.2006, 1.80356, 5.26251, 7.69162, 2.66808, 9.20424, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 2.99639, 3.58278, 8.15589, 26.4175, 10.5788, 15.9233, 8.57665, 13.2933, 31.0573, 4.49579, 6.17261, 4.05822, 15.4694, 7.11106, 9.76691, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 9.49099, 2.55673, 5.77367, 7.00715, 11.8172, 14.795, 16.0516, 19.9916, 16.3954, 15.4181, 7.79475, 11.086, 15.1181, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 3.73535, 27.1326, 18.5295, 6.75536, 11.8236, 27.2383, 5.3328, 37.6127, 5.25804, 0, 3.3404, 0, 7.06448, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 12.5286, 11.6368, 3.9737, 12.405, 10.6027, 9.8324, 11.086, 5.96115, 0, 0, 0, 5.56628, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10.7097, 34.7669, 19.8458, 25.6916, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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 EvtD0topipietaPlot::initProbMax(){
        noProbMax();
}
 
void EvtD0topipietaPlot::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(xmass12);
                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 EvtD0topipietaPlot::FindXBin(double mass2) {
  if (mass2 < Xmin) { return 0; }
  else if (mass2>=Xmax) { return 51; }
  else { return int((mass2-Xmin)/Xwid)+1; }
}
 
int EvtD0topipietaPlot::FindYBin(double mass2) {
  if (mass2 < Ymin) { return 0; }
  else if (mass2>=Ymax) { return 51; }
  else { return int((mass2-Ymin)/Ywid)+1; }
}