EvtDtoKKpiPlot.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: EvtDtoKKpiPlot.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/EvtDtoKKpiPlot.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;
 
EvtDtoKKpiPlot::~EvtDtoKKpiPlot() {}
 
void EvtDtoKKpiPlot::getName(std::string& model_name){
        model_name="DtoKKpiPlot";     
}
 
EvtDecayBase* EvtDtoKKpiPlot::clone(){
        return new EvtDtoKKpiPlot;
}
 
void EvtDtoKKpiPlot::init(){
 
    checkNArg(0);
 
    bool idN  = getDaugs()[0]==EvtPDL::getId(std::string("K-"))||getDaugs()[0]==EvtPDL::getId(std::string("K+"));
    bool idKs = getDaugs()[1]==EvtPDL::getId(std::string("K+"))||getDaugs()[1]==EvtPDL::getId(std::string("K-")) ;
    bool idPi = getDaugs()[2]==EvtPDL::getId(std::string("pi+"))||getDaugs()[2]==EvtPDL::getId(std::string("pi-")) ;
    if(!(idN && idKs && idPi ) ){std::cout<<"EvtDtoKKpiPlot: the daughter sequence should be K- K+ pi+"<<std::endl;abort();}
    EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId());
 
Xmin = 0.35;
Xmax = 1.91;
Xwid = 0.026;
Ymin = 0.35;
Ymax = 1.91;
Ywid = 0.026;
avm1  = 4.0;
double HisPDFtmp[62][62] =
{
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0139227, 0.0242858, 0.0102631, 0.019465, 0.00821469, 0.00941803, 0.0196361, 0.00814691, 0.0270675, 0.0188643, 0.0186632, 0.0107977, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0140485, -0.00172283, 0, 0, 0, 0, 0, 0.00962505, 0, 0, -0.000873327, 0, 0.0188526, 0.00761928, 0.0154787, 0, 0.0088296, 0.00889335, 0.0150876, 0.0114478, 0.0180075, 0.0107661, 0.00843499, 0.0130738, 0.0455281, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0134615, 0, 0, 0, 0, 0.0272361, 0.00927996, 0.00887302, 0.00990029, 0.0113715, 0.00920733, -0.000819338, 0.0132542, 0.0195743, 0.00879929, 0.0218734, 0.00831937, 0, 0.0182388, -0.000808949, -0.000794811, 0.0166262, 0, 0.0184469, -0.000680394, 0.00879748, 0, -0.00190808, 0.00985711, -0.00151276, 0.0136869, -0.00110282, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, -0.00224346, -0.00158143, 0, 0.0111645, 0.00963723, 0.00837308, 0, 0, -0.00162554, -0.000811716, -0.00103891, 0, 0.00906974, -0.00080413, -0.000816707, 0, 0, 0.00980291, -0.00183136, 0.0101183, -0.000930091, 0, 0.0185092, 0.0157808, -0.000819735, 0.0196731, -0.000909308, 0, 0.0084228, 0, 0, -0.00270583, -0.000871727, 0, -0.000823917, 0, 0, -0.00290892, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.000883683, 0, 0, 0, -0.00099019, 0, 0, -0.000890028, 0.0107049, 0.00887844, 0, 0.00868849, 0.0078208, 0, -0.00102005, 0.0107161, 0, 0.00602897, 0, 0.028913, -0.00101468, 0, 0.00864125, 0, 0, -0.000888573, -0.000867477, -0.000792421, 0.00912084, 0.0171944, -0.000897794, 0, 0.0104999, 0, 0, -0.00449721, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, -0.00106087, 0, 0, 0, 0, 0, 0, -0.00102738, 0, 0.0200809, 0, 0, -0.000936822, -0.00102728, 0, 0, -0.000944456, 0, -0.000737237, 0.00781316, 0, 0, 0, -0.000795664, 0, -0.00140018, 0, -0.000874954, 0.00909987, -0.00215512, -0.000967849, -0.000805499, 0, 0, 0.00965563, -0.0019969, 0, 0.00854241, 0, 0.0407456, 0.00778027, 0.00917508, 0.0078959, 0.0118948, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, -0.00101511, 0, 0, 0, 0.00886395, 0.00957819, 0, 0, 0, 0.0184697, 0, -0.00101187, 0, -0.00172169, 0.00763349, -0.000895323, 0, -0.00252987, 0.019748, 0.00947828, 0, 0.0192401, 0.00849129, -0.000888861, 0.01189, 0.00896836, 0.00997021, 0.00931974, 0, -0.00170736, 0, 0, -0.000825473, 0, 0.00736716, 0.00858116, 0.0414927, -0.000862929, 0.0108653, 0.00711648, 0, 0.0424946, 0.02056, 0.00888727, 0.0183671, 0.0300872, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, -0.00121692, 0, 0, 0, 0, 0.00980403, 0, 0, 0, -0.000906888, 0.00818758, -0.000931843, 0, -0.000784871, 0, -0.00103353, 0, 0, -0.000859876, 0.00731771, -0.00067335, 0, -0.000905807, -0.000790764, 0.00889752, -0.00070466, 0, 0.0190332, 0, 0, 0.00846915, 0.0110697, 0.00897852, 0.00804933, 0.0201357, 0.00851472, 0, 0.0364704, 0.0211352, 0.0120208, 0.0163991, 0, 0.0236853, 0.0129102, 0.00920625, -0.00154954, 0.0123733, 0.0351667, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, -0.000847885, 0, 0, 0, 0, 0, 0, 0, 0.00934519, 0.00877748, -0.00110685, 0.0100934, 0, 0.0087101, 0.0101893, 0.00770844, 0, 0, 0.01073, 0.00907416, -0.00108313, -0.000690736, 0, 0.0198926, 0, -0.000960642, 0, 0.0114602, -0.000867882, 0, -0.0039687, 0.0108437, 0, 0, -0.00214432, 0.0335963, 0.00767771, 0.00919692, 0.0221316, 0.0167495, 0.00842359, 0, -0.000850634, 0.0139215, 0.010582, 0.00922736, 0.00934376, 0.0640697, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, -0.00139959, 0, 0, 0, 0.0312581, 0.0111713, 0.00825327, -0.000705715, 0.00968094, 0.0079305, 0, 0, 0.0118746, 0, 0, 0.00892758, 0, 0.0102401, 0.0193089, -0.000782156, -0.000885337, 0, 0, 0.00795678, 0, -0.000709818, 0, 0.0155405, 0.00799595, 0, -0.00111399, 0.0105511, 0.0191352, 0, 0.00911081, 0.00939029, 0.0280986, 0.0187336, -0.000793975, 0, 0, 0.0082419, 0.0074632, 0.009078, 0.0357759, 0.0102688, 0.0155454, 0.0498312, 0.017454, 0.0439061, 0.0371816, 0.036918, -0.00125509, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0101589, 0, 0.0395604, 0.0114227, -0.00100364, 0.0213129, 0.0228551, 0, 0, 0, 0.0104772, 0, 0, 0, -0.00091241, 0, -0.0026225, 0.0163946, -0.00064703, -0.000948443, -0.000943019, 0, 0, 0.00841675, 0, 0.00975579, 0, 0, 0, -0.00113992, 0, 0, -0.000855679, 0.0107493, 0.0101579, 0.00928391, 0.00901179, 0, 0.0129731, 0.0198509, 0.0237226, 0.0195224, 0.0268942, 0.0243642, 0.0211605, 0.0359483, 0.0801728, 0.0314633, 0.0332897, 0.0300546, 0.0429153, 0.0455037, 0.0694133, 0.0888124, 0, 0, 0, 0},
{0, 0, 0, 0.0489109, 0.00944428, 0, 0.0173018, 0.0162537, -0.00171801, 0.00829744, 0.0160517, 0, 0.0108203, 0.023278, 0.00844664, 0, 0, -0.00106547, 0.00857296, 0, 0, -0.000827661, -0.000848852, 0, -0.000990982, 0.00832698, -0.00100791, -0.00111151, 0, 0, 0.00859719, 0, -0.000750289, 0.0095587, 0, 0.00928949, 0.0111698, 0.00905906, 0.0154512, -0.000941774, 0.015085, 0.00894646, 0.0198978, 0.0191158, 0.0114518, 0.0109502, 0.0208813, 0.0151289, 0.0259387, 0.0446887, 0.056272, 0.0624139, 0.075904, 0.0216623, 0.092942, 0.0574898, 0.0116131, 0.0468348, 0.042885, 0, 0, 0},
{0, 0, 0, 0.0825589, 0.0221408, 0.0392638, 0.0372548, 0, 0.0281633, 0.0218577, 0, 0.0384219, 0.0113054, 0, 0, 0.0101912, 0.00908258, 0, -0.00165883, -0.000676169, 0, -0.000752131, 0, 0, 0, 0, -0.000812433, 0, 0, -0.000921229, 0, -0.00102465, -0.000717399, 0, -0.00123682, 0, 0.00834694, -0.00118759, 0.00792323, 0.0186343, 0, 0.0301626, 0.0097278, 0, 0.018005, 0.0310391, 0.0104501, 0.0496578, 0.0281955, 0.0165578, 0.0455259, 0.035723, 0.0701144, 0.0881308, 0.0787734, 0.144527, 0.0778648, 0.0562271, 0.0186825, 0, 0, 0},
{0, 0, 0, 0.0698628, 0.16651, 0.105634, 0.0477428, 0.06, 0.0159642, 0.048638, 0.0439435, 0.0187432, -0.000797068, -0.000719579, 0.0207966, 0.0185039, -0.00211158, 0.00953271, 0.0212814, 0, 0.0107911, 0, 0.0104335, 0.00800181, 0.00937909, 0, 0, 0, 0, 0, 0, 0.0102564, 0.0161653, -0.000845697, 0, 0.00823757, 0.0289298, 0.0598826, 0.0111307, 0.00817876, 0.0318672, 0.0770479, 0.0260889, 0.0568727, 0.0695652, 0.0331282, 0.0441718, 0.0661653, 0.0929879, 0.0956522, 0.0896045, 0.0614806, 0.146617, 0.0557448, 0.0764612, 0.15177, 0.11039, 0.0760611, 0.131421, 0, 0, 0},
{0, 0, 0, 0.0844833, 0.0848093, 0.0990991, 0.0516922, 0.0821633, 0.0580325, 0.0579193, 0.057767, 0.0835098, 0.0265978, 0.00721743, -0.000881274, 0.0530418, 0.0379593, 0.00763684, 0.010234, -0.000857086, 0.0197255, 0, 0.00623986, 0.0290557, 0, 0.007842, -0.00105588, -0.000780333, 0.00930577, 0, 0, 0.00649399, 0.00729907, 0.00778127, 0.0308451, 0.0164112, 0.0479435, 0.00776708, 0.0518251, 0.0613207, 0.0714823, 0.07, 0.114109, 0.0487135, 0.0616141, 0.102306, 0.112254, 0.148481, 0.103235, 0.0919618, 0.130377, 0.205457, 0.116028, 0.166836, 0.1459, 0.128874, 0.128328, 0.106844, 0.111886, 0.0611888, 0, 0},
{0, 0, 0, 0.0835314, 0.0221192, 0.112568, 0.0762611, 0.058406, 0.0732344, 0.0692292, 0.055456, 0.0457485, 0.0486681, 0.0481529, 0.0340297, 0.0407047, 0.0280471, 0.0178287, 0.0172723, -0.000836188, 0.00740479, 0.00645997, 0.00971044, 0.00563776, 0.00868042, 0, -0.00153558, -0.000762348, -0.00119712, 0.00938967, 0.00671685, -0.000859723, -0.00099193, 0.00906069, 0.0252852, 0.00772808, 0, 0.0179303, 0.0499147, 0.0057039, 0.0629322, 0.028777, 0.0574818, 0.113005, 0.0967862, 0.0470863, 0.127491, 0.0842513, 0.195622, 0.0993212, 0.185804, 0.12411, 0.146365, 0.170778, 0.197943, 0.259737, 0.107905, 0.154536, 0.177368, 0.0925204, 0, 0},
{0, 0, 0, 0.111431, 0.0654541, 0.0807881, 0.0784929, 0.00795811, 0.0194426, 0.0339753, 0.0677406, 0.0450864, 0.0176048, 0.0196949, -0.000871355, 0.0218863, 0.0161006, 0.0533049, 0.0188992, 0.0378347, 0.00799579, 0.00794461, -0.0009632, 0.0086048, 0.0195132, 0.0211498, -0.000866141, 0, -0.00272947, -0.00176606, 0.00767947, -0.000982462, 0.0104655, 0.0193164, -0.000979882, 0.0171381, 0.00888631, 0.00887055, 0.00893744, 0.0207668, 0.0102099, 0.0718604, 0.052395, 0.0580129, 0.10178, 0.061127, 0.0282911, 0.0639713, 0.0528053, 0.142383, 0.107398, 0.0833401, 0.138912, 0.101017, 0.0642299, 0.0739587, 0.0799474, 0.0678708, 0.109511, 0, 0, 0},
{0, 0, 0, 0.0345705, 0.0293617, 0.0399925, 0.0544398, 0.0200117, 0.00733173, 0.0304692, 0.0426771, 0.0171076, 0.0449947, 0.0242359, 0.0218323, 0.00666036, 0.00687741, 0.0301849, 0.00875201, -0.000985995, 0.00922699, -0.000853764, -0.00083159, 0.00878674, 0, 0, -0.00159877, 0, -0.000724172, -0.000823628, 0, 0.0179337, 0.00984429, 0.034952, 0.0105243, -0.000820197, 0.0176194, 0.0336909, 0.0344162, 0.0406053, 0.0239255, 0.0249303, 0.0207042, 0.0205343, 0.0578339, 0.0451814, 0.0176266, 0.0284754, 0.0276654, 0.0825147, 0.0652409, 0.0672455, 0.119474, 0.062069, 0.0100719, 0.0201064, 0.0101344, 0.0135435, 0.0301031, 0.0376412, 0, 0},
{0, 0, 0, 0.0100569, 0.0093066, 0.00717103, 0.0273388, 0.0070312, 0.0170305, 0.01851, 0.0316539, 0.0186476, 0.00801258, 0.0298422, 0.0162819, 0, 0.00823547, 0, 0.0246154, -0.000997358, 0, 0.0247961, -0.000903635, 0, 0, -0.000988488, 0, 0, 0, 0.00784142, 0.00947002, 0, 0, 0, -0.000766148, 0.0201519, 0.0107802, 0.00888957, 0.0252956, 0.0092967, 0, 0.0165119, 0.0311205, 0.00844633, 0, 0.00759511, 0.00904772, 0.0569211, 0.0219707, 0.0201324, 0.0250272, 0.0178588, -0.00105352, 0.0476843, 0.0251409, 0.0207664, 0.0307567, 0.0220836, 0.0205404, 0.0465057, 0, 0},
{0, 0, 0, 0.00929319, 0, 0.00957557, -0.000775727, 0, 0.0103202, 0, 0, 0, 0.00601536, 0.0330814, 0.01049, -0.00177138, 0.0117269, 0, 0.00818093, 0.00919637, 0.00743863, 0.00665334, 0.00949445, 0, 0, 0, 0, -0.000980402, 0.00944467, 0.0103012, 0, 0.010084, 0.0178529, 0, 0.00904177, -0.000788797, 0, 0.011637, 0.0174232, 0.00806195, 0.0261961, 0, 0.0184518, 0.0206251, 0, 0.0115893, 0.0170932, 0.0312629, 0.0340296, 0.0237025, 0, 0.0275816, 0.0356968, 0.0179104, 0.0328156, 0.0123295, 0.0637828, 0, 0, 0, 0, 0},
{0, 0, 0, 0, -0.00101025, -0.000858735, 0.00870872, 0.0202326, 0, 0, 0, 0, -0.000847364, 0, 0.0389124, 0, 0, 0.0069042, 0.0117152, 0, 0, 0.00749609, 0.0156619, 0.0180838, 0.00901454, 0.00925037, 0, -0.000838362, 0, 0.0104666, 0, 0.0089198, -0.00105854, 0, 0, 0.00941842, 0.0175091, 0.0457886, -0.00157401, 0.0180961, 0.0178693, 0, 0.0258238, 0.0236452, 0.0299774, 0.0208024, 0, 0.0278071, 0.00925399, 0.0437028, 0, 0.0163869, 0.0413868, 0.00910784, 0.0192142, 0.0366326, 0.00841407, -0.0012014, 0.0166667, 0.0436019, 0, 0},
{0, 0, 0, 0.0203152, 0.00949449, 0, -0.000807642, 0.0100947, 0, 0.00798749, 0, 0.0104772, -0.000759592, 0.0285344, -0.000809837, 0.00857392, 0.00926016, 0, 0, 0, 0.0085656, 0, -0.00090472, 0.0101704, 0.00903735, -0.00268799, 0, -0.000837338, 0.00714379, 0.0190639, 0.0199808, 0.0147623, 0.0104001, 0.0297673, 0.0588285, 0, 0.0238077, 0.0195378, 0.0230434, 0.0106557, 0.0297837, 0.0146021, 0, 0.00961409, 0.0100678, 0.0105691, 0.0106211, 0.00962249, -0.00119916, 0.0276335, 0.0223966, 0, 0.0306002, 0.0202501, -0.0010004, 0, 0.0134923, 0.0132922, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0164774, 0.00842688, 0.00911396, 0.00949977, -0.000850423, 0, 0, 0, -0.00105277, -0.000862762, 0.0105271, -0.000807488, 0, 0.00810359, 0, 0, 0, 0.0101164, 0, 0, 0, 0, 0.00842445, 0.0114901, 0.0266553, 0, 0, 0.00915325, 0.0164874, 0.00967876, 0, 0, 0.023019, 0.00997442, 0.0181452, 0.0352483, 0.0235438, 0.00947753, 0.0210715, 0.010443, 0.00770524, 0.0121609, 0.0179744, 0.0279225, 0.0380326, 0.0326071, 0.0259083, 0, 0.0476909, 0.0203499, 0, 0.0108538, 0.0438583, 0.0124589, 0.0649872, 0, 0, 0},
{0, 0, 0, 0, 0, -0.000677822, 0, -0.00196125, -0.00170722, 0.00917402, 0.0168193, 0.0160035, 0.0187455, 0, 0, 0.0105069, 0.00653315, 0.00946888, 0.00939433, -0.000951599, 0.00793651, 0, 0, 0.00800202, 0, 0, 0, 0, 0.0240316, 0, -0.000909594, 0, -0.000833545, 0.00925294, 0.0407852, 0.00822491, 0, 0, -0.000796518, 0.0428913, 0, 0, 0.00864616, 0.0081278, 0.0198069, -0.000900054, 0.0129539, 0.0291967, -0.00107986, 0.0204132, 0.0175669, 0.00921452, 0.020746, 0.00976493, 0, 0, -0.00179696, 0, 0, 0.0412541, 0, 0},
{0, 0, 0, 0.0203955, 0.0106986, 0.00881007, 0.0191054, 0, 0, 0.00945626, 0.0104001, 0.0222121, 0.0085498, 0.00947087, 0.00954733, 0, 0.0116447, 0, 0, 0, 0.0112924, 0, 0.00894737, -0.000928966, -0.000824152, 0, 0.0169718, 0.00829953, -0.0008101, 0.0087907, 0.0212713, 0.00962365, 0.00818169, 0.0186238, 0.0443356, 0.0088327, 0.00907982, 0.0092243, -0.00154134, 0.00745146, 0.0181104, 0, 0.0174172, 0, 0.0188855, 0.0111528, 0.0104965, 0.0187198, 0.00907934, 0, 0.0101033, 0.0132822, 0, 0.0255092, 0.0117876, 0.00884271, 0.0388322, 0.0348326, 0.0268515, 0, 0, 0},
{0, 0, 0, 0.0126706, 0.0172168, 0.0209784, 0.0115417, 0.0198827, 0, -0.000613454, 0, 0, -0.0012482, 0.00666148, 0, 0, -0.000914029, -0.000802899, 0.00751363, 0, 0.0210709, 0, 0.00956299, 0.00755021, 0.0108458, 0, 0.00851648, 0.00970513, 0.00984537, 0.0114522, 0.00811426, 0.014856, 0, 0, 0, 0.0356321, 0.00915935, 0.0079448, 0.0112329, 0, 0.0195858, 0.0361888, 0.0197721, 0.00967225, 0.0190026, 0, 0.0116118, 0.0173056, 0.0256392, 0, 0, 0.0101425, 0, 0.00986239, 0.0204131, 0.0313983, 0.0328703, 0.0150389, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.00925089, -0.000860052, 0.0113342, 0.0200476, 0, 0, 0.0109967, 0, 0.0092987, 0.00817958, 0, 0, -0.000883344, 0.0114325, 0.0287869, 0.0182725, 0.00878291, 0.0194701, -0.000979345, 0.0103969, 0, 0.008714, 0, 0.0102515, -0.000783578, 0.00882312, 0.0181932, 0, -0.00082272, 0.0251304, 0.0243868, 0, 0.0176539, 0.0276354, 0.011233, 0.0169369, 0, 0, 0, -0.000914585, 0, 0.0108333, 0.0105352, 0.0189954, 0.026487, 0.0219047, 0.0214605, 0, 0.0117922, 0.0112767, 0.0208396, 0.0206242, 0, 0, 0, 0},
{0, 0, 0, 0.0215772, 0, 0.0469387, -0.00199938, 0, 0.00806495, 0.0113939, 0.00784219, 0.0102418, 0.00848231, 0, 0, -0.00238053, 0, 0, 0.0108451, 0.0224638, 0.00895348, 0.0104752, 0.00931011, 0, 0.00806999, -0.000835198, 0.00828157, 0, 0, 0.010863, 0, 0.00966184, 0, 0.00880196, 0, 0, 0.00967279, 0, -0.00104673, 0, 0.0193632, 0.0294464, 0.0207633, 0.0180889, 0, 0.0196807, 0.0113285, 0.0207271, 0.00870174, 0.0122184, 0.0084899, 0, 0.0191684, 0.0251521, 0.0244355, 0.0117797, 0.0433962, 0.107924, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0242275, 0, 0, 0.0096565, -0.00155124, 0.0208565, 0, 0, 0.0294917, -0.00180637, 0.00820202, -0.000719485, 0.017647, 0.00939541, 0.00995025, 0, 0.031661, 0.0179949, 0.0088753, 0.00956716, 0, 0.00942728, 0.0182113, 0.019223, 0.0100184, 0, 0.00826353, 0.0109375, 0, 0.0095989, 0.0088364, 0, 0.0177186, 0.0176286, 0, 0.00918478, 0.00979515, 0.0718151, 0.00712146, 0.0388281, 0, 0.0169886, 0.027248, -0.000841679, 0, 0.011249, 0.010187, 0, 0, 0, 0.064453, 0.0347804, 0.14096, -0.0104849, 0, 0, 0, 0},
{0, 0, 0, 0.0304846, 0.00895394, 0, 0, 0.0095143, -0.000817347, 0.0100501, 0.00882872, 0, 0.00960198, 0.00871875, 0.024524, 0.0173022, 0, 0.0126656, 0.0101352, 0.0331152, 0.0228368, 0.0179239, 0, 0.00906672, 0.0325413, 0, 0.00905138, 0, 0.00894021, 0.0172658, 0, 0, 0.00862803, 0, 0.0324659, 0, 0.0155981, 0.0206338, 0.0230017, 0.0108395, -0.000768815, 0.0380435, 0.00846202, 0.0431805, 0.0186941, 0.0185413, 0.0203616, 0.0412492, 0.0285326, 0.0174774, 0.0220788, 0.00855435, 0.0185348, 0.0155923, 0.0847581, 0.17733, 0.384751, 0, 0, 0, 0, 0},
{0, 0, 0, 0.0549451, 0, -0.000809818, 0.0341507, 0.00924517, 0, 0, 0, 0.025915, 0, 0, 0, 0.0100309, 0, 0.00773839, 0, -0.000922284, 0.00824142, 0.00998856, 0.0185869, 0, 0, -0.000930184, 0.00774234, 0.0271866, 0.0261441, 0.00859834, 0.0101072, 0.0106873, 0.0511486, 0.00915786, 0.00992035, 0.00784197, 0.0202428, 0.00991919, 0.0294389, 0, 0.00836111, 0.0168758, 0.0119083, 0.0087859, 0, 0.0168192, 0.00894654, 0, 0.0111521, 0, 0, 0.020193, 0.00745592, 0.043302, 0.33012, 0.756053, 4, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0193347, 0, 0, 0.00877543, 0.00926148, -0.000978804, 0, 0, 0.00616838, 0.0480958, -0.00110401, 0, -0.00094806, 0, 0.0266272, 0.00926523, 0, 0.0470984, 0.0316436, 0.0102558, 0.0137878, 0.00956818, 0, 0.0109435, 0.00795692, 0.00996638, 0, 0.0262011, 0, 0, 0, 0.00821991, 0.00822732, 0.0178063, 0, 0.0131448, 0.0187651, 0.0376607, 0, 0.00832224, 0.035451, 0.0111331, 0.00716381, 0, 0.0200805, 0, 0.0247405, 0.0343828, 0.073858, 0.162882, 0.753372, 0.877119, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.00892485, 0.0087992, 0.00866911, 0, 0, 0.0183706, 0.00892193, 0.0295628, 0, 0.029044, 0, 0.0264391, 0.0130717, -0.000957429, -0.00102977, 0, 0.0217497, 0, 0.0119743, 0, 0.00814037, -0.00216151, 0, 0.00901897, 0.0800344, 0.0311763, 0, -0.000999707, 0.0175548, 0.00776979, 0.00997534, 0.0306324, 0.0163388, 0.00966387, 0, 0.0313085, 0.0438893, 0.0189919, 0.0370669, 0.032723, 0.0253304, 0.0216285, -0.000945404, 0.0126512, 0.0105467, 0.0114826, 0.0474927, 0.229064, 0.468742, 0.499709, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.0102398, 0, 0.00963375, 0.0256482, 0.0205175, 0, 0.00772904, 0.0187665, 0.00648284, 0.0302874, 0.0335031, 0, 0.0321651, 0.0210876, 0, 0.0164572, 0.00975459, 0.0144711, 0.0349518, 0.0103055, 0, 0.00774833, 0, 0, 0, 0.00881944, 0.00900243, 0.0166373, 0.0190377, 0.0121791, 0.0290389, 0.0218861, 0.0282682, 0.00957256, 0.0210033, 0.0188598, 0.0160125, -0.00114259, 0, 0, 0.0232449, 0.0218724, 0.0136043, 0.00759345, 0.0100437, 0.0400415, 0.129409, 0.382832, 0.499742, 0.169623, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0217826, 0.00954403, 0.00978454, 0.0314549, 0.00763196, 0.0107647, 0.0477279, 0.049444, 0.00818761, 0.0191171, 0.0187272, 0.00827722, 0.00809931, 0.00807353, 0.00865924, 0.0440351, 0.0188898, 0.0178678, 0.0187399, 0, 0.0353044, 0.0199281, 0, 0.0092765, 0, 0.0341788, 0.0137857, 0, -0.000917367, 0, 0.0222421, 0.0175874, 0.0179275, 0.0111494, 0.0397602, 0.0228989, 0.0195081, 0.0173706, 0, 0.00988815, 0.00821166, 0.0105519, 0, 0.0425583, 0.0290794, 0.0377965, 0.144596, 0.31633, 0.309859, 0.407018, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0.0373188, -0.000800659, 0, 0.0206369, 0.0201636, 0.024405, -0.00164074, 0.0322904, 0.00772539, 0.017094, 0, 0.0102001, 0.0170068, 0.00789728, 0.0285433, 0.017654, 0.0209081, 0.0353514, 0.00991326, 0, 0, 0.0112883, 0.0272436, 0.0221948, 0, 0.0303313, 0, 0.00892256, 0.0112371, 0.00807134, 0.0216746, 0.0174218, 0, 0, 0.00866308, 0.0088484, 0.0820871, 0.0221963, 0.00996698, 0.00907869, 0.00964537, 0.0253436, 0.00854222, 0.0357143, 0.0486163, 0.0954776, 0.259473, 0.149898, 0.052, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.00925421, -0.000826828, 0.0114796, 0.0186418, -0.000892872, 0, 0.0120563, -0.000841373, 0.0396193, 0.00921667, 0, -0.000839938, 0, 0.0216088, 0, 0.00909285, 0.00989559, 0.0187188, 0.0143512, 0.0166337, 0, 0, 0.0433543, 0, 0.00886178, 0, 0.0337719, 0.00869666, 0.00954008, 0.0229141, 0.00920752, 0.00859955, 0.016442, 0.0264347, 0.016196, 0.0310447, 0, 0.0232209, 0.0200668, 0.0194759, 0.0252697, 0.00791894, 0.0218616, 0.0886278, 0.107202, 0.120529, 0.017576, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.0176408, 0.0402815, 0.00900193, 0.00910931, 0, 0.00995597, 0, 0.0193276, 0.02625, 0.0102957, 0.0175661, 0.00946315, 0.0235473, 0.0224084, 0.0174483, 0.0167777, 0.0246835, 0.0131795, 0.0464314, 0.018863, 0.00888001, 0.019547, 0.00803371, 0.0103648, 0.00857505, 0.0340203, 0, 0.0257845, 0, 0.020093, 0.0185917, 0.009653, 0.0361675, 0.0340909, 0.0165112, 0.0213852, -0.000725483, 0.00793697, 0.00797994, 0.0197249, 0.0117026, 0.0102868, 0.0365747, 0.124825, 0.0937213, 0.0546665, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.0158892, 0.00842231, 0.00791289, 0.0107855, 0.0479872, 0, 0, 0.0337553, 0.0166636, 0.0180018, 0.0818542, 0.00921793, 0, 0, 0.0161146, 0.0104829, 0.0112961, 0.00880284, 0, 0.00894671, 0.0321728, 0.00823045, 0.00836536, 0.0116709, 0.0169351, 0, 0.0504936, 0.0272736, 0.00943627, 0.0163528, 0.00930591, 0.0102789, 0, 0.0395859, 0.0206268, 0.00947595, 0.0454126, 0.0282035, 0.0215823, 0.0425532, 0.0490323, 0.0663051, 0.093707, 0.0430499, 0.204082, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.0474398, 0.010514, 0.0171634, 0.0472313, 0.0225918, 0.00918337, 0, 0.021168, 0.0214649, 0.0295599, 0, 0.0215697, 0.00883459, 0.0191107, 0.0101801, 0.00683388, 0.0288497, 0.0176917, 0.0115658, 0, 0.0205727, 0.0160846, 0.00916362, 0.00781841, 0.00923558, 0.00764855, 0.0115033, 0.0476958, 0.00914876, 0.0340413, 0.0260466, 0, 0.0203144, 0.0224558, 0, 0.00836209, 0.032166, 0, -0.00201207, -0.000999997, 0.0364455, 0.0197025, 0.0580674, 0.025409, 0.20625, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.0204216, 0.0178749, 0.0339385, 0.0476663, 0.0214505, 0, 0.010011, 0.00825038, 0, 0.0106965, 0.0298976, 0.0240339, 0.00813984, 0.0264443, 0, 0.0377711, 0.0081673, 0.0264415, 0, -0.00081143, 0.00870287, 0.0430039, 0.00920434, 0.00927802, 0.0102824, 0.00919913, 0.0321815, 0.00880137, 0.0167547, 0.0141489, 0.0282248, 0.0100758, 0.0192092, 0.00911131, 0.0191987, 0.00988181, 0.0094089, 0.037713, 0.0211831, 0, 0.0403058, 0.0343354, 0.0401235, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.0644633, 0.0346, 0.0371551, 0, 0.0102819, 0.00970625, 0.0103222, 0.010594, 0.00810723, 0.047961, 0.102065, 0.0250843, 0.0208742, 0, 0.0254557, 0, 0.0363382, 0.0186905, 0, 0.008927, 0.0300059, 0.0378217, 0.0508027, 0, 0.0479418, 0.026139, 0.0105877, 0.0322802, 0.0196786, 0.0184422, 0.0167361, 0.0191199, 0.0250955, 0, 0.0225517, 0.0328446, 0.0514563, 0.0103672, 0.0577811, 0.0380105, 0.0258166, 0.0629371, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0.0231183, 0.0173595, 0.0199675, 0, 0.0368378, 0.0089956, 0.00849399, 0.0356221, 0.00735485, 0.00863388, 0.0196979, 0.050897, 0.0251263, 0.0181164, 0.00793501, 0.0370859, 0.00934374, 0.0167843, 0.017893, 0.00876617, 0, 0, 0.0275779, 0.0243458, 0.026663, 0.0332204, 0.0393579, 0.0128048, 0.0185263, 0.00752273, 0.0539919, 0.0286538, 0.017602, 0.0211793, 0.0304102, 0.0377677, 0.0875269, 0.0560436, 0.0186506, 0.0278015, 0.0493056, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0.0160049, 0, 0.0104456, 0.0449646, 0, 0.0203922, 0.0206122, 0.054633, 0.0527813, 0.00926439, 0.0356753, 0.0354575, 0.0316874, 0.00742951, 0.0183808, 0.0190379, 0.0430335, 0.00998737, 0.0183085, 0.0152626, 0.0340017, 0.0507289, 0, 0.0104491, 0.0116689, 0.00989029, 0.05604, 0, 0.010513, 0, 0.0129241, 0.0253465, 0.00932121, 0.062275, 0.0260865, 0.101612, 0.0374271, 0.0134468, 0.0345664, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0.0401984, 0.0239462, 0.0323053, 0.0299613, 0.00808859, -0.000959072, 0.028936, 0, 0.0360868, 0.0583658, 0.0340565, 0, 0.0185563, 0.0109007, 0, 0.0106144, 0, 0.0538789, 0.0100871, 0.0264126, 0.0181142, 0.00924878, 0.0230294, 0, 0.0237235, 0, 0.0100862, 0.0336292, 0.0208131, 0.042786, 0.0274457, 0.0377862, 0.0349889, 0.072125, 0.146772, 0.0632134, 0.0768693, 0.118474, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0.0465492, 0.0388832, 0.0211017, 0.0219103, 0.0224301, 0.0082513, 0.0485614, 0, 0.0513242, 0, 0.0205069, 0.0318284, 0.0180626, 0.0181801, 0.0312115, 0.00763562, 0.021414, 0.00981205, 0.0211101, 0.0296362, 0, 0.0104593, 0.0100904, 0.0180371, 0.0174462, 0.0241922, 0.0450072, 0.00850443, 0, 0.0186514, 0.0374268, 0.0605321, 0.2021, 0.209671, 0.0577898, 0, 0.0483516, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.0447205, 0.0353172, 0.0401819, 0.048029, 0.0330707, 0.0328415, 0.0186309, 0.0436886, 0.0277081, 0.0419765, 0.0218716, 0.0285834, 0.0168439, 0.0103291, 0.0385388, 0.021661, 0, 0.0257292, 0.0564093, 0, 0.0214601, 0.0109343, 0, 0.0549645, 0.0105792, 0.0403816, 0.0173734, 0.0182437, 0.0456348, 0.0563464, 0.172808, 0.23266, 0.0694025, 0.0296717, 0.0758808, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.0234194, 0.00993789, 0.0487651, 0.0231491, 0.0203606, 0.0232817, 0.0218335, 0.0238186, 0.0277103, 0.01239, 0.0181719, 0.0307454, 0.0087968, 0.0112293, 0.02037, 0.00871786, 0.021266, 0.0102832, 0.0190427, 0, 0.0172336, 0.0683699, 0.0352564, 0.00866234, 0.0291314, 0.0111943, 0.0182796, 0.0334819, 0.0346112, 0.313737, 0.389412, 0.0889875, 0.0436609, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0.0985222, 0.0230706, 0.0257939, 0.00966886, 0.0216262, 0.0181988, 0.0174742, 0.0201692, 0.0268839, 0.0126801, 0, 0.0184841, 0, 0.0277268, 0.00982999, 0.0135447, 0.0431554, 0.0553891, 0.0198105, 0.00958762, 0.0308985, 0.0303113, 0.0342678, 0.00921841, 0.0196959, 0.0658576, 0.0100004, 0.193647, 0.28473, 0.542411, 0.204171, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.03174, 0.0105849, 0.0771186, 0.0279586, 0.0266612, 0.0318761, 0.0559981, 0.0390625, 0.0121277, 0.0320627, 0.0104402, 0.0432972, -0.00110143, 0.00831363, 0.00758654, 0.027497, 0.0120548, 0.0234965, 0.00981269, 0.0075152, 0.0344042, 0.0575213, 0.0467659, 0.0353293, 0.0774836, 0.121804, 0.280416, 0.482075, 0.287566, 0.0886228, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0947666, 0.0470336, 0.0443913, 0.0497682, 0.0266242, 0.0106494, 0, 0.00941109, 0.0125906, 0, 0.0104469, 0.0107968, 0.039149, 0.0413446, 0.0306129, 0.00857951, 0.00911935, 0.0358557, 0, 0.0229856, 0.0215492, 0.0385325, 0.0731537, 0.0260051, 0.15089, 0.420476, 0.753433, 0.444165, 0.286792, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.12383, 0.0510436, 0.0110302, 0.00967458, 0.0349852, 0.0344665, 0.0105125, 0.0104745, 0.0249594, 0.0187075, 0.00857108, 0.0176199, 0.0384546, 0.0209939, 0.0202804, 0.0102498, 0.018094, 0.017984, 0.0481076, 0.0963358, 0.0563204, 0.0829453, 0.0781776, 0.337149, 0.878505, 0.719033, 0.222222, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0588235, 0.0230188, 0.0364726, 0, 0.0119376, 0.0163343, 0.0487118, 0, 0.0367248, 0, 0.0277698, 0.0270057, 0.0212871, 0.0320826, 0, 0.0110496, 0.0271488, 0, 0.0169567, 0.0667761, 0.205392, 0.51, 1.19001, 1.39583, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.166667, 0.0190713, 0.0228799, 0.0135813, 0, 0.0411091, 0.0132791, 0, 0.0332337, 0.0365666, 0.0218031, 0.0308543, 0.0353304, 0.0472783, 0.0345769, 0.0405673, 0.0216447, 0.0711784, 0.0323605, 0.130971, 0.199812, 0.798479, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0334781, 0.0112575, 0.0227968, 0.0277248, 0.0143044, 0.0257028, 0.0131532, 0, 0.013922, 0.0491049, 0.0103723, 0.028983, 0.0598813, 0.0451858, 0.0420784, 0, 0.128964, 0.152036, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0162025, 0, 0, 0.0130208, 0, 0.0292517, 0.0137187, 0.0675021, 0.0272235, 0.0277778, 0.0194071, 0.0628173, 0.0737442, 0.0212331, 0.113636, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.226268, 0.051353, 0.0232558, 0, 0.0439695, 0, 0.0276163, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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<62;i++){
  for (int j=0; j<62; j++) {
     HisPDF[i][j] = HisPDFtmp[i][j];
  }
}
 
}
 
void EvtDtoKKpiPlot::initProbMax(){
        noProbMax();
}
 
void EvtDtoKKpiPlot::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 EvtDtoKKpiPlot::FindXBin(double mass2) {
  if (mass2 < Xmin) { return 0; }
  else if (mass2>=Xmax) { return 61; }
  else { return int((mass2-Xmin)/Xwid)+1; }
}
 
int EvtDtoKKpiPlot::FindYBin(double mass2) {
  if (mass2 < Ymin) { return 0; }
  else if (mass2>=Ymax) { return 61; }
  else { return int((mass2-Ymin)/Ywid)+1; }
}