EvtDtopipi0etaPlot.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: EvtDtopipi0etaPlot.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/EvtDtopipi0etaPlot.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;
 
EvtDtopipi0etaPlot::~EvtDtopipi0etaPlot() {}
 
void EvtDtopipi0etaPlot::getName(std::string& model_name){
        model_name="Dtopipi0etaPlot";     
}
 
EvtDecayBase* EvtDtopipi0etaPlot::clone(){
        return new EvtDtopipi0etaPlot;
}
 
void EvtDtopipi0etaPlot::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("pi0")) ;
    bool idPi = getDaugs()[2]==EvtPDL::getId(std::string("eta"));
    if(!(idN && idKs && idPi ) ){std::cout<<"EvtDtopipi0etaPlot: the daughter sequence should be pi+ pi0 eta"<<std::endl;abort();}
    EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId());
 
Xmin = 0.0751;
Xmax = 1.796;
Xwid = 0.034418;
Ymin = 0.4692;
Ymax = 3.098;
Ywid = 0.052576;
avm1  = 84.9113;
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, 5.55681, 0, 0, 4.46103, 3.29413, 0, 9.58826, 2.7939, 4.71479, 11.4296, 6.95352, 4.30188, 7.47747, 4.91315, 1.37074, 2.73427, 1.82034, 4.27089, 2.25947, 5.27089, 1.24789, 2.09578, 1.7907, 0, 0.389202, 0.0557514, 0.738381, 4.26268, 1.92207, 0, 11.0493, 13.3667, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 10.7484, 0, 5.83725, 10.8372, 7.11291, 7.19155, 14.1362, 1.93052, 4.66714, 20.3981, 3.25822, 2.29789, 1.65886, 1.89549, 9.23944, 1.11972, 8.95352, 1.01643, 0.997496, 2.68592, 2.50297, 5.19155, 1.90516, 1.0892, 0.952676, 1.31136, -0.107511, 1.56682, 1.43531, 2.04859, 1.23122, -0.17406, 1.93208, 0.651644, 5.10673, 2.41041, 3.78141, 3.34953, 5.60399, 19.97, 11.5568, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 4.59648, 3.36582, 4.24845, 4.1661, 6.0493, 3.59728, 6.33975, 7.22848, 4.48462, 3.97265, 3.02019, 7.13271, 3.01615, 1.70955, 1.56111, 9.92207, 8.5007, 4.42101, 2.72853, 5.58005, 1.07954, 4.19609, 6.35775, 2.47005, 5.19906, 1.92101, 2.57731, 1.9385, 1.37758, 0.275901, -0.720891, 0.828035, -0.158215, 4.65153, 2.0769, 6.3308, 9.70501, 3.29765, 22.0096, 10.1227, 24.5746, 19.6376, 6.60399, 0, 0, 0},
{0, 0, 0, 0, 5.08005, 8.40563, 3.1635, 7.17606, 8.05478, 3.38526, 5.54585, 14.9549, 10.3509, 10.9111, 9.97747, 7.50344, 8.31056, 1.37598, 1.70315, 0.70133, 1.66465, 6.08005, 1.53216, 11.446, 2.32723, 5.78341, 3.91565, 4.14939, 0.995775, 3.67962, 1.92089, 2.71465, 1.54038, 1.15963, 5.95352, 0.485342, 1.98087, -0.0850699, 1.45474, 2.56004, 5.95055, 13.9116, 7.57887, 11.4217, 0, 7.19455, 3.73349, 3.13409, 5.86667, 6.11362, 0, 0},
{0, 0, 0, 0, 10.3352, 17.5418, 0.832316, 10.931, 14.8995, 20.8634, 34.5925, 17.5897, 9.01678, 4.34038, 6.46103, 4.90528, 5.28627, 1.98177, 2.00986, 3.47265, 2.11401, 1.56397, 2.6392, 3.09804, 2.73826, 6.07981, 2.11804, 5.49249, 1.98838, 3.36937, 8.58005, 0, 1.47646, 3.43779, 0.0861977, 0.552348, 0.0736156, 0.70317, 3.80031, 11.4036, 6.33952, 36.9263, 19.773, 52.0221, 6.25531, 13.6225, 6.93169, 1.71252, 1.7754, 5.89202, 0, 0},
{0, 0, 0, 10.684, 0, 1.53709, 3.54288, 12.9909, 9.85845, 8.34476, 9.08279, 5.82606, 18.3906, 5.21972, 4.42923, 1.92887, 2.96479, 3.48462, 2.58779, 1.74249, 1.99399, 3.29789, 5.00047, 4.10626, 3.65117, 3.71854, 6.91315, 12.5089, 3.89707, 8.31878, 2.56479, 1.70955, 1.41315, -0.385289, 0.881691, 1.32019, 4.44601, 1.47426, 21.639, 9.56776, 17.0822, 9.08404, 13.5972, 10.8356, 7.17953, 4.85462, 5.35306, 1.48795, 2.82723, 5.44883, 0, 0},
{0, 0, 13.5883, 5.10399, 1.87872, 8.327, 6.44941, 9.44601, 5.37365, 8.97239, 13.382, 30.3695, 14.6526, 7.03287, 3.67033, 6.33451, 2.78141, 0.709312, 3.5349, 4.10399, 0.331573, 1.56635, 0.540376, 1.58441, 3.87872, 7.0561, 6.61972, 2.66714, 2.7784, 4.31808, 6.14867, 4.23122, 1.01813, 0.964789, -0.0566193, 2.61186, 0.420136, 2.05884, 5.81467, 0, 9.33897, 17.1086, 12.421, 11.4946, 25.0892, 15.2901, 8.754, 2.53339, 1.40223, 0.559628, 11.1615, 0},
{0, 0, 4.15939, 2.12997, 0, 0, 0, 8.65188, 7.95579, 46.4174, 24.4707, 6.24147, 5.68648, 7.37383, 9.08075, 2.96901, 6.96995, -0.300938, 5.6047, 1.24398, 2.69871, 0.545384, 2.18156, 4.0374, 1.18013, 0.95806, 5.54038, 2.46056, 2.48052, 1.57183, 1.14108, -0.404225, 0.468897, 7.95493, -0.118544, 2.2862, 5.40634, 3.74147, 4.58572, 7.97268, 11.7001, 18.3202, 48.8362, 13.649, 6.37681, 9.48263, 1.26229, 4.28944, 2.98157, 2.29554, 0, 0},
{0, 18.2394, 2.60986, 0, 5.18607, 4.05884, 8.74695, 40.0042, 4.85176, 11.3106, 8.54695, 12.788, 7.24695, 2.49249, 4.34338, 1.78627, 0.803568, 0.0818916, 1.12094, -1.31596, 0.489202, 0.500235, 0.907177, 1.23873, 0.780673, 1.28216, 2.69495, 2.17488, 1.00047, 2.0416, -0.190141, 0.702066, 0.766433, 0.750705, 2.00141, 0.897418, 4.17958, 3.32254, 9.39953, 41.9113, 13.7698, 43.4695, 4.35278, 5.94239, 9.14479, 12.125, 3.01737, 1.13686, 14.2916, 0, 0, 0},
{0, 3.22254, 9.03287, 4.24445, 2.63886, 5.525, 7.66385, 17.0021, 20.8038, 10.6736, 38.1465, 5.90074, 9.83177, 3.12535, 2.01315, 0.762207, 0.600282, 2.61972, 0.710681, 1.58005, 0.624726, 2.7784, 0.376338, 1.64789, 0.64331, 1.08529, 0.89468, 2.38134, 2.35364, 1.99249, 0.901174, 1.88545, 0.0401413, 2.66761, 1.427, 1.97711, 20.8455, 8.61565, 3.52535, 13.923, 4.58075, 42.2629, 20.957, 4.92871, 8.71655, 8.67496, 7.83557, 8.61612, 12.3859, 0, 0, 0},
{0, 5.54859, 2.79984, 15.0493, 6.56229, 14.669, 9.79127, 9.57664, 4.27635, 9.06467, 17.1451, 6.74695, 6.87919, 7.91385, 7.52465, 3.10364, 1.87962, 1.45556, 1.14306, 0.741002, 0.780673, -0.184898, 0.296871, 0.262207, 0.542645, 1.19108, 0.508357, 3.09577, 1.65117, 0.911362, 0.96925, 1.24945, -0.0566193, -0.649764, 5.29734, 4.49343, 3.44216, 3.11865, 10.6068, 5.5998, 13.981, 10.967, 14.3736, 23.5761, 11.0182, 6.74479, 2.95493, 11.416, 0, 0, 0, 0},
{0, 4.26338, 3.65775, 0, 7.73533, 26.1793, 13.192, 17.9485, 12.382, 10.2948, 9.86369, 6.52451, 0, 1.30879, 3.32594, 0.619437, 1.28592, 1.60028, 1.13134, 1.04038, 2.28091, 1.46103, 0.961034, 0.545834, 0.492489, 3.52465, 1.34943, 1.56135, 0.790704, 1.40117, 0.593271, 1.01714, 0.184695, 0.0921236, 5.05751, 2.72442, 4.03696, 4.04894, 8.15258, 21.3059, 20.1793, 10.9554, 7.50247, 11.2587, 5.9036, 3.22413, 5.47888, 0, 0, 0, 0, 0},
{0, 3.68357, 0, 4.05681, 6.45329, 25.8634, 9.48347, 5.67285, 5.93108, 32.862, 6.61671, 1.67171, 6.8687, 0.773334, 5.90704, 1.98975, 0.516731, 1.81534, 1.0416, 1.54038, 0.612475, 3.02465, 1.1669, 1.09549, 1.02019, 4.84272, 0.762324, 0.742489, 0.901174, 1.25994, 0.893134, 2.24085, 0.667607, 4.96244, 3.4565, 4.84237, 5.12042, 8.25243, 4.74302, 24.8977, 6.16197, 7.13221, 3.26855, 14.8777, 5.32564, 4.26064, 0, 0, 0, 0, 0, 0},
{0, 2.79127, 3.61017, 14.2491, 4.30628, 6.53537, 6.46651, 5.30509, 17.7887, 6.68347, 3.15981, 4.0804, 6.327, 8.47888, 3.21549, 1.75768, 9.09719, 2.36275, 1.32676, 1.00035, 0.691818, 1.70407, 0.590704, 1.19871, 2.30986, 2.78341, 0.730517, 6.89061, -0.248043, 1.13615, 0.000564003, 3.65399, 5.12113, 3.50892, 3.09934, 11.1341, 1.60763, 7.73404, 17.5836, 5.32066, 1.9198, 14.0371, 11.1355, 0, 17.5803, 2.69226, 0, 0, 0, 0, 0, 0},
{0, 7.9946, 4.94398, 5.0616, 7.87629, 3.07778, 20.1061, 10.2559, 12.5468, 2.47283, 8.22441, 6.28185, 11.4324, 0.524649, 2.79824, 1.85751, 1.50446, 1.71455, 3.04859, -0.237676, 2.35739, 0.409038, 0.190704, 0.958967, 1.50548, 2.07183, 0.868232, 1.19399, 2.07481, 1.08791, 1.4608, 2.26338, 1.56843, 0.772303, 5.61256, 8.59272, 17.9399, 4.14122, 13.0876, 5.41096, 4.7169, 6.2169, 8.85509, 31.9455, 14.9242, 8.52676, 11.877, 0, 0, 0, 0, 0},
{0, 1.75857, 0.415966, 2.92851, 4.36667, 5.41071, 5.3547, 7.98944, 6.6054, 7.44053, 16.877, 3.89883, 1.04517, 2.85399, 3.12751, 2.24945, 1.497, 2.78858, 5.98498, 0.835306, 1.58005, 3.35775, -0.202112, -0.12277, 2.91565, 0, 0.799139, 1.48052, 0.492489, -0.0499214, 0.994273, 2.04859, 3.02876, 8.00689, 2.46718, 4.34423, 7.51373, 8.81162, 6.49116, 10.6446, 9.31565, 13.4653, 11.2887, 9.17596, 0, 29.6075, 0, 0, 0, 0, 0, 0},
{0, 4.11816, 5.2878, 2.45225, 6.50963, 4.15352, 5.9696, 10.9549, 3.35324, 4.92136, 2.47747, 3.96582, 4.08416, 3.21014, 6.87418, 3.32423, 5.99319, 2.77812, 2.47265, 1.28142, 1.83202, 3.09578, -0.320375, 0.448096, 0.829695, 0.693584, -0.0524252, 0.272066, 0.242489, 1.87418, 2.06808, 12.5747, 3.39455, 3.81261, 3.53608, 6.97031, 12.8761, 3.98987, 8.18474, 8.29969, 9.65294, 9.0777, 5.4969, 16.4639, 26.862, 0, 0, 0, 0, 0, 0, 0},
{0, 4.2739, 2.86964, 2.25088, 1.92707, 15.392, 6.21371, 8.0804, 6.74366, 30.6075, 7.35571, 4.38693, 4.83416, 2.87991, 7.65188, 4.66704, 0, 1.04488, 3.78889, 2.7784, 7.19906, 1.03202, -0.087089, 0.147552, 2.63615, 0.365258, 1.69906, -0.10669, 1.25994, 1.19495, 1.80803, 8.7716, 3.91706, 30.3709, 14.8291, 13.4629, 5.14479, 37.185, 9.91902, 25.7869, 12.05, 9.36038, 10.9919, 30.0056, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 6.06843, 5.49934, 12.3562, 18.8359, 6.46197, 10.2709, 7.40814, 17.4228, 34.3052, 10.6202, 4.16828, 2.42911, 8.86322, 4.1338, 4.37812, 1.72778, 3.15439, 3.10876, 3.42136, 0.432019, 0.567723, 0.188397, 0.931377, 0.345801, 1.75897, 0.285916, 9.89061, 2.75035, 4.17512, 4.10535, 1.91725, 6.89026, 12.6418, 11.5527, 7.29421, 10.8446, 4.23738, 13.7669, 7.41756, 20.3613, 13.0613, 14.4024, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 13.2237, 10.0288, 4.84726, 26.4009, 11.9296, 6.21596, 4.92451, 9.87966, 9.20317, 9.97449, 6.63169, 7.70454, 5.14867, 3.44246, 4.09167, 4.26376, 8.80235, 4.93161, 1.82183, 0.886463, 2.38169, 0.54939, -0.121283, 1.08404, 0.675118, 1.73599, 0.0353056, 5.63756, 4.07362, 2.70563, 27.1479, 4.26291, 10.3362, 5.37903, 2.30315, 8.38756, 6.43803, 6.01878, 9.53768, 11.1245, 27.5761, 21.8291, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 8.29186, 2.88948, 5.95953, 6.04809, 4.57533, 7.482, 9.43333, 22.9857, 24.5268, 5.52692, 2.59354, 11.2599, 22.5732, 6.16315, 13.9453, 0.740258, 3.2892, 4.49249, 1.39969, 3.28091, 0.566354, 1.82582, 2.61972, 1.1856, 1.30712, 1.69906, 1.11722, 0.797559, 4.72442, 1.35548, 9.65305, 4.01756, 18.2984, 10.2573, 6.15196, 11.5856, 9.18064, 4.09613, 7.67618, 18.8284, 44.1779, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 28.0178, 8.98087, 7.86404, 5.63535, 10.2904, 8.8867, 4.08659, 16.5411, 52.4624, 12.6758, 3.19085, 13.0246, 5.79002, 2.86937, 3.24789, 2.36937, 8.23944, 10.3981, 1.53193, 0.911362, 0.500353, 1.22244, 0.293897, 0.85227, 0.408451, 1.13615, 3.1651, 2.218, 6.97066, 9.63005, 0, 6.33146, 6.69672, 7.00329, 4.31606, 6.11773, 8.98686, 7.68999, 25.6883, 22.1629, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 5.54413, 14.3742, 10.7998, 29.0725, 14.5523, 12.3577, 9.83005, 6.57183, 14.9378, 4.5374, 3.57594, 3.71273, 4.89331, 7.86596, 3.52805, 1.35897, 1.5431, 10.1601, 2.397, 2.65939, 0.664848, 1.97402, 0.136346, 2.28662, 0.486198, 1.63756, 0.0508614, 8.07395, 4.51549, 4.29005, 10.9759, 19.2977, 10.1191, 7.02641, 1.98987, 8.5439, 6.48404, 6.11667, 12.2272, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 9.22801, 14.0725, 8.77512, 11.4491, 10.5688, 17.3343, 5.36097, 13.2394, 13.6362, 11.2552, 2.18625, 3.29687, 1.86133, 1.36002, 1.91899, 2.15904, 3.24898, 1.39926, 1.04515, 2.95806, 4.71479, 1.56135, 2.35268, 2.42911, -0.113802, 0.807356, 1.14795, 10.3457, 11.6383, 40.6897, 3.93094, 7.04382, 3.27109, 10.2887, 8.12899, 3.57735, 6.14507, 5.70548, 21.3216, 13.4474, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 6.51232, 4.8937, 0, 10.6714, 27.8427, 14.0274, 19.2634, 7.16725, 4.02757, 3.90528, 3.7743, 4.73681, 2.58005, 5.81808, 2.7907, 3.97676, 1.37848, 0.619367, 1.65117, 0.606573, 0.743099, 2.04002, 6.47747, 0, 1.81315, 6.63615, 5.59327, 7.9101, 4.77197, 6.93506, 10.2914, 9.13822, 20.02, 7.89883, 12.2194, 4.96639, 17.2601, 20.1944, 11.5582, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 23.4296, 6.07033, 10.4689, 10.3102, 41.7648, 17.2285, 12.0484, 4.21315, 3.51096, 3.09578, 4.13052, 3.46376, 3.90904, 1.92094, 2.4446, 1.95782, 1.83083, 1.86047, 1.19056, 1.04789, 4.01643, 1.59167, 3.86823, 2.25661, 6.02739, 37.0507, 9.19155, 16.1446, 11.472, 16.3412, 9.70751, 17.9878, 10.605, 8.02357, 11.654, 29.4817, 4.26409, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 7.19399, 5.70681, 8.76883, 7.13271, 14.4535, 5.73326, 4.50571, 4.1277, 7.77089, 6.48498, 2.18075, 3.55268, 2.89707, 2.14071, 4.15368, 5.3338, 4.99249, 19.7469, 9.74695, 5.42958, 2.76854, 4.60376, 3.47653, 0, 6.94038, 12.1089, 0, 64.9592, 24.7748, 7.18905, 8.35505, 17.1715, 5.29033, 11.2012, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 6.02144, 21.2709, 15.7998, 11.2268, 21.7962, 6.54009, 0, 7.84726, 4.60986, 4.09878, 9.36596, 1.93219, 1.65418, 2.23831, 4.25418, 10.8263, 1.21437, 15.7469, 1.93669, 5.01643, 9.12864, 19.0178, 6.80686, 9.70798, 10.275, 10.7919, 9.87747, 20.6371, 15.1485, 26.6321, 7.93065, 8.0169, 2.95981, 5.22825, 4.07324, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 16.1833, 4.88236, 0, 9.67809, 4.24648, 7.7877, 36.2723, 1.75172, 2.24898, 6.88134, 5.37484, 8.46604, 11.7155, 5.22801, 2.21854, 1.42709, 1.22244, 5.27543, 7.15939, 4.70955, 5.06397, 4.65939, 7.06432, 11.1865, 22.6218, 30.6697, 15.8745, 7.02191, 12.6208, 27.5313, 32.273, 13.5014, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 4.60603, 3.83631, 8.36561, 5.21205, 38.5897, 9.38169, 2.17512, 3.64617, 5.00446, 1.49051, 3.74249, 1.67488, 7.49249, 3.8892, 2.80516, 0.939034, 3.8297, 2.23873, 1.67962, 1.6416, 6.70454, 3.28948, 8.63204, 10.1757, 10.7442, 24.6002, 23.3941, 21.1406, 31.7176, 84.9113, 24.4196, 47.5446, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 7.69906, 14.9453, 5.16166, 4.23052, 4.27817, 26.3352, 14.7962, 4.00357, 18.3038, 14.97, 3.29413, 8.04859, 2.6302, 1.96068, 0.834675, 6.90563, 1.16948, 4.11174, 2.14967, 2.44817, 4.36254, 3.69906, 7.18779, 7.76512, 11.2241, 18.3831, 22.4224, 34.9002, 60.8948, 23.3826, 0, 41.4639, 0, 16.8277, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 16.9535, 12.3181, 14.2073, 4.26232, 2.69272, 2.33631, 4.47662, 8.38991, 9.08826, 5.8313, 20.8592, 2.05747, 10.0808, 12.2873, 2.12794, 5.4939, 0.0887327, 1.00571, 0.997277, 5.57778, 4.71822, 8.98157, 7.54335, 12.5725, 34.9167, 41.1068, 23.8856, 9.84202, 9.99162, 14.1706, 47.0357, 10.0986, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 6.47653, 3.3297, 2.30712, 9.00892, 4.19085, 1.81174, 4.69495, 7.42958, 5.34578, 9.62794, 2.44734, 2.9266, 0.402426, 8.36667, 1.12347, 1.02465, 1.59327, 3.15188, 3.48247, -0.0250378, 6.50438, 15.0609, 16.0498, 6.97047, 15.449, 17.0362, 75.4653, 17.3862, 19.6519, 43.4639, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 10.3817, 5.82629, 3.65904, 2.17812, 10.6526, 2.24398, 3.36275, 3.92887, 3.92423, 3.10364, 2.44495, 3.00047, 3.37183, 1.87347, -0.727565, 2.09578, 2.31283, 1.54232, 3.83177, 0.622458, 28.7512, 6.08302, 58.4817, 16.589, 19.1481, 14.525, 78.115, 43.9085, 12.6097, 25.1451, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 4.07981, 2.56225, 13.2394, 5.7784, 3.16554, 3.744, 3.3338, 12.0883, 7.47152, 0.377981, 1.50446, 0.0320593, -0.744834, 3.81056, 0.832301, 0.161436, 0.94781, 5.07718, 6.35981, 10.674, 20.1291, 10.452, 21.1242, 5.95282, 12.9854, 14.8423, 12.2709, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 4.79413, 4.89707, 3.72003, 3.87962, 4.54413, 2.51465, 4.95305, 4.32308, 2.78865, 4.09878, 0.324615, 5.50892, 0.796636, 0.186073, -0.110516, 2.1302, -0.0514854, 8.6054, 6.87829, 4.66784, 17.8441, 8.55657, 9.12549, 34.6211, 25.1765, 45.3202, 7.86322, 20.7005, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 14.461, 4.40657, 2.40657, 8.22301, 2.29725, 3.54896, 4.92094, 2.87333, 6.89331, 3.75446, 2.27662, 6.37418, 1.97747, 11.431, 0.732864, 1.26589, 0.886934, 9.75658, 7.52676, 6.99355, 4.84225, 43.7991, 42.6075, 8.23909, 10.3926, 18.1615, 0, 0, 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.92477, 6.15368, 5.90516, 6.3867, 3.7939, 2.28865, 3.40986, 2.84155, 5.64617, 2.62272, 2.36186, 0.53216, 2.96901, 0.366198, -0.030045, 0.830584, 1.86132, 3.32179, 2.88981, 12.0822, 4.16633, 14.7326, 7.33737, 10.6211, 0, 13.3038, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 21.2545, 15.0643, 3.54515, 2.33881, 6.52441, 4.98462, 15.6047, 23.3352, 3.14009, 2.36554, 0.223006, 1.19427, 0.313303, 1.07629, 2.04073, 2.38404, 1.81071, 3.07289, 4.05478, 7.86737, 3.95352, 10.3083, 9.98568, 5.0169, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 5.80986, 4.18674, 2.95258, 1.9421, 2.93052, 4.18709, 5.17512, 2.01894, 1.90535, 5.84225, 0.884195, 3.39218, 0.653289, 1.05305, 14.7362, 3.36256, 7.70118, 9.39859, 8.19703, 4.05646, 2.01972, 0, 0, 5.57324, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 18.8577, 2.36526, 16.5568, 1.73302, 1.82308, 2.55634, 2.4108, 2.39427, 0.774649, 1.68504, 2.29961, 2.6806, 14.6404, 5.63709, 2.93232, 4.16433, 4.66385, 3.4739, 1.94304, 9.83451, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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.80986, 1.78113, 1.00297, 2.81808, 1.96461, 3.42547, 5.29233, 0.694524, 5.2108, 6.70352, 4.21174, 1.53136, 4.55716, 1.88075, 1.91697, 7.6054, 18.4953, 2.20939, 0, 10.1601, 0, 0, 0, 0, 0, 0, 0, 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.461, 0, 1.25446, 3.00595, 1.2951, 3.68948, 5.67856, 2.30759, 4.69566, 6.43369, 3.9185, 4.90599, 3.02028, 5.88704, 6.30681, 1.74481, 4.23909, 3.85446, 3.73193, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 7.65117, 1.77019, 7.81056, 0.807356, 1.0561, 3.10235, 7.23443, 1.94102, 10.1451, 5.89402, 16.2087, 3.89812, 5.025, 14.8831, 9.57324, 6.83451, 1.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, 0, 0, 0, 0, 6.19155, 10.0958, 5.95352, 2.84225, 2.83887, 3.69606, 1.90993, 2.15653, 8.20869, 7.1601, 3.22027, 2.38169, 19.5432, 4.20204, 0, 14.5418, 4.55681, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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.85775, 7.68404, 4.73498, 17.2737, 1.86432, 0.953052, 6.09719, 2.109, 18.2237, 3.35296, 14.1136, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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.37621, 5.54585, 4.13866, 4.44601, 6.31127, 1.56254, 4.50481, 9.39812, 3.66251, 8.58826, 5.46103, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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.41862, 3.32019, 1.74765, 3.12042, 3.85274, 1.02418, 3.84726, 16.97, 14.0329, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 6.27089, 2.46925, 4.04859, 2.95493, 1.06573, 0, 8.52535, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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.28733, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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 EvtDtopipi0etaPlot::initProbMax(){
        noProbMax();
}
 
void EvtDtopipi0etaPlot::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 EvtDtopipi0etaPlot::FindXBin(double mass2) {
  if (mass2 < Xmin) { return 0; }
  else if (mass2>=Xmax) { return 51; }
  else { return int((mass2-Xmin)/Xwid)+1; }
}
 
int EvtDtopipi0etaPlot::FindYBin(double mass2) {
  if (mass2 < Ymin) { return 0; }
  else if (mass2>=Ymax) { return 51; }
  else { return int((mass2-Ymin)/Ywid)+1; }
}