BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtD0toKSpipiPlot.cc
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2//
3// Environment:
4// This software is part of models developed at BES collaboration
5// based on the EvtGen framework. If you use all or part
6// of it, please give an appropriate acknowledgement.
7//
8// Copyright Information: See EvtGen/BesCopyright
9// Copyright (A) 2006 Ping Rong-Gang @IHEP
10//
11// Module: EvtD0toKSpipiPlot.cc
12//
13// Modification history:
14//
15// Liaoyuan Dong August, 2022 Module created
16//
17//------------------------------------------------------------------------
18#include "EvtGenBase/EvtPatches.hh"
19#include "EvtGenBase/EvtParticle.hh"
20#include "EvtGenBase/EvtGenKine.hh"
21#include "EvtGenBase/EvtPDL.hh"
22#include "EvtGenBase/EvtReport.hh"
23#include "EvtGenModels/EvtD0toKSpipiPlot.hh"
24#include "EvtGenBase/EvtRandom.hh"
25#include <stdlib.h>
26#include <string>
27
28#include "TFile.h"
29#include "TApplication.h"
30#include "TROOT.h"
31#include "TH1.h"
32#include "TAxis.h"
33#include "TH2.h"
34using std::endl;
35
36EvtD0toKSpipiPlot::~EvtD0toKSpipiPlot() {}
37
38void EvtD0toKSpipiPlot::getName(std::string& model_name){
39 model_name="D0toKSpipiPlot";
40}
41
42EvtDecayBase* EvtD0toKSpipiPlot::clone(){
43 return new EvtD0toKSpipiPlot;
44}
45
46void EvtD0toKSpipiPlot::init(){
47
48 checkNArg(0);
49
50 bool idN = getDaugs()[0]==EvtPDL::getId(std::string("K_S0"))||getDaugs()[0]==EvtPDL::getId(std::string("K_L0"));
51 bool idKs = getDaugs()[1]==EvtPDL::getId(std::string("pi+"))||getDaugs()[1]==EvtPDL::getId(std::string("pi-")) ;
52 bool idPi = getDaugs()[2]==EvtPDL::getId(std::string("pi-"))||getDaugs()[2]==EvtPDL::getId(std::string("pi+")) ;
53 if(!(idN && idKs && idPi ) ){std::cout<<"EvtD0toKSpipiPlot: the daughter sequence should be K_S0/L0 pi+ pi-"<<std::endl;abort();}
54 EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId());
55
56Xmin = 0.072919;
57Xmax = 3.1229;
58Xwid = 0.0435712;
59Ymin = 0.070375;
60Ymax = 3.1068;
61Ywid = 0.0433775;
62avm1 = 1.89744;
63double HisPDFtmp[72][72] =
64{
65{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
66{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
67{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
68{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
69{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
70{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
71{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
72{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
73{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00665709, 0.0140415, 0.00937031, 0, 0.0335591, 0.0268942, 0.0319302, 0.0336182, 0.082483, 0.0870691, 0.112121, 0, 0, 0, 0, 0, 0, 0},
74{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0482456, 0.0841121, 0.0469484, 0.0234102, 0.0166739, 0.033103, 0.00938078, 0.0034632, 0.00821013, -0.000631524, -0.000642641, 0.0026595, 0.0146251, 0.0194118, 0.0184746, 0.0840415, 0.0687152, 0.057089, 0.0789847, 0.0282787, 0.0280054, 0, 0, 0, 0, 0},
75{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0664898, 0.0652448, 0.0567227, 0.0824385, 0.0992671, 0.0610264, 0.0573868, 0.0573106, 0.0220205, 0.0261934, 0.00844282, 0.00577553, 0.0058175, 0.00297489, 0, 0.00647316, 0.00713415, 0.00989317, 0.0256505, 0.0736834, 0.116204, 0.0670578, 0.0686343, 0.0267318, 0.00918274, 0.0160713, 0.04311, 0, 0, 0, 0},
76{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0521657, 0.0640025, 0.0225772, 0.0829845, 0.0640486, 0.0643463, 0.0918833, 0.0728832, 0.0270935, 0.0350269, 0.0215905, 0.0135926, 0.0131966, 0.011553, 0.00228931, -0.00125326, 0.00596057, -0.00124719, 0.0212871, 0.0265947, 0.0268578, 0.104012, 0.0607158, 0.0337294, 0.00544824, 0.0143701, 0.0103259, -0.000541093, 0.0197782, 0, 0, 0, 0},
77{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.48, 0.0174734, 0.0391642, 0.0233992, 0.0343596, 0.0347096, 0.0447667, 0.0431331, 0.0591761, 0.0630347, 0.0477246, 0.0589755, 0.0535137, 0.0269502, 0.0258217, 0.0358101, 0.00280487, 0.00717157, -0.00122758, 0.00303194, -0.000947307, 0, 0.0118396, 0.0156714, 0.0507447, 0.0770331, 0.0474667, 0.038169, 0.0445178, 0.0274048, 0.0118788, 0.0167624, 0.0144626, 0.0259354, 0, 0, 0, 0},
78{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0347594, 0.0436214, 0.051157, 0.0429697, 0.0247733, 0.0260423, 0.0316134, 0.0247604, 0.039553, 0.0190151, 0.0247182, 0.0507373, 0.0433553, 0.0513938, 0.0397958, 0.0234792, 0.0215416, 0.0193312, 0.00232806, 0.00308476, -0.000714832, 0.00750574, 0.00287595, 0.00927556, 0.0106164, 0.0294285, 0.0933387, 0.056418, 0.0604638, 0.0300376, 0.0263117, 0.0185801, 0.0130549, 0.0178711, 0.0214502, 0.0166943, 0, 0, 0, 0},
79{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0280112, 0.0218703, 0.0154006, 0.0365659, 0.014061, 0.0482109, 0.0247056, 0.032143, 0.0114928, 0.0244444, 0.0456769, 0.0419506, 0.0697134, 0.0386185, 0.0444058, 0.0631189, 0.0427899, 0.0311286, 0.0120179, 0.00208352, 0.00339164, 0.00600556, 0.0126691, -0.00111723, 0.00237435, 0.00971825, 0.0198351, 0.0668728, 0.059098, 0.0301087, 0.0212451, 0.0222507, 0.00582636, 0.0193004, 0.0189585, 0.0104787, 0.018143, 0, 0, 0, 0, 0},
80{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0195625, 0.0334728, 0.0101507, 0.0157114, 0.0219979, 0.0351727, 0.00447182, 0.00700638, 0.0375301, 0.0237089, 0.0242707, 0.0237428, 0.0332346, 0.0368631, 0.032507, 0.032666, 0.019498, 0.0187895, 0.0199937, 0.0102989, 0.00542893, 0.00249367, 0.00323879, -0.000304143, 0.00309895, 0, 0.012236, 0.00541206, 0.0491979, 0.0414711, 0.0323574, 0.0214055, 0.0120577, 0.0046972, 0.00621582, 0.0100278, 0.00861363, 0.0124362, 0.0393549, 0, 0, 0, 0, 0},
81{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.0440366, 0.01841, 0.0107117, 0.0228789, 0.02912, 0.0187117, 0.0188283, 0.015881, 0.0192973, 0.0313393, 0.0166246, 0.0136062, 0.00688693, 0.00687653, 0.0162378, 0.0255552, 0.0387605, 0.0294502, 0.0219534, 0.0342041, 0.0277876, 0.0117542, -0.00187756, 0.00522609, 0.00280581, -0.000376457, -0.000968331, 0.00293612, 0.00639277, 0.0168144, 0.0382846, 0.0604078, 0.0310467, 0.0326504, 0.0279678, 0.0151885, 0.0214212, 0.009829, 0.00262037, 0.0212361, 0.00788974, 0, 0, 0, 0, 0, 0},
82{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.0431732, 0.042595, 0.0114624, 0.0101843, 0.0213717, 0.00894055, 0.013536, 0.0252132, 0.0272387, 0.0103668, 0.0262014, 0.00330917, 0.0261236, 0.00619379, 0.0270658, 0.0201552, 0.0126421, 0.0288047, 0.0184506, 0.0282731, 0.0106444, 0.00314039, 0.00589669, 0.00783229, 0.00325287, 0.00250649, -0.000324501, -0.000299896, 0, 0.00898276, 0.0353091, 0.0570606, 0.0302146, 0.00620741, 0.0114111, 0.0105752, 0.00771583, 0.0106497, 0.00542358, 0.0053125, 0.0200179, 0, 0, 0, 0, 0, 0, 0},
83{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.0222134, 0.00451063, 0.00991126, 0.0196285, 0.00623905, 0.0111046, 0.00620571, 0.00480961, 0.0250631, 0.0171972, 0.0153988, 0.0193861, 0.013, 0.0178049, 0.0278429, 0.0204597, 0.0203379, 0.0202649, 0.020751, 0.0298537, 0.0274692, 0.0186618, 0.00272465, 0.00287571, 0.00675187, -0.000346927, 0.00668447, -0.000317979, 0.00690156, 0.0129432, 0.0369478, 0.0251226, 0.0290798, 0.0407402, 0.0260144, -0.000254573, 0.0177707, 0.00561764, 0.0153272, 0.0323023, 0.0345694, 0, 0, 0, 0, 0, 0, 0, 0},
84{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0317712, 0.0310795, 0.020292, 0.0200376, -0.00177238, 0.020682, 0.0132242, 0.0164988, 0, 0.00656001, 0.00985303, 0.00921659, 0.0112816, 0.0167344, 0.0162036, 0.0324328, 0.0364591, 0.00659146, 0.0308419, 0.0167365, 0.0109779, 0.0119099, 0.0168579, 0.0107088, 0.00287193, -0.000323915, 0.00289249, 0.00890632, 0, 0.00661909, 0.00567849, 0.0191331, 0.0284907, 0.0217304, 0.0386303, 0.0256592, 0.0267499, 0.0169885, 0.0186455, 0.0102709, 0.0276731, 0.0248053, 0, 0, 0, 0, 0, 0, 0, 0},
85{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.00465951, 0.0702822, 0.0447753, 0.0352661, 0.0248955, 0.0211488, 0.0231082, 0.0118607, 0.00946683, 0.0169729, 0.0148534, 0.00951346, 0.00809812, 0.0191495, 0.0176832, 0.0100624, 0.0151274, 0.0163589, 0.0131189, 0.0204952, 0.00948417, 0.0147846, -0.000603759, -0.00108483, 0.0060653, -0.000577155, -0.000294198, 0.00354415, -0.000863502, 0.00554958, 0.017512, 0.0247718, 0.0196635, 0.0299396, 0.00792212, 0.0535733, 0.0125685, 0.0143, 0.0293399, 0.0270971, 0.0123162, 0.00957137, 0, 0, 0, 0, 0, 0, 0, 0, 0},
86{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0865971, 0.12919, 0.0312158, 0.0377967, 0.0145567, 0.00619967, 0.0184486, 0.0163661, 0.0196879, 0.00886918, 0.00568387, 0.0128739, 0.00342042, 0.016514, 0.0302179, 0.00725212, 0.0155085, 0.0198493, 0.0184007, 0.0164522, 0.0117381, 0.00927088, 0.00487471, 0.00517525, 0, 0.00500405, 0.0029352, 0, 0, 0.00889407, 0.0222106, 0.028095, 0.0221151, 0.0164953, 0.0156452, 0.00858472, 0.0207883, 0.0247836, 0.0143985, 0.0159322, 0.0305775, 0.0344795, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
87{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0696162, 0.154043, 0.114228, 0.0551163, 0.0309737, 0.0191352, 0.0385291, 0.0111306, 0.00659411, 0.0336275, 0.0115686, 0.0076806, 0.028192, 0.00694772, 0.0231914, 0.0125673, 0.0152482, 0.0107726, 0.013359, 0.0230453, 0.00557517, 0.0157646, -0.00116608, 0.00619098, -0.000267461, 0.00346956, 0.00295168, -0.000291885, 0.00288836, 0.0080369, 0.0125315, 0.0115649, 0.0457929, 0.0130741, 0.00997373, 0.0223435, 0.00499894, 0.0129204, 0.0187247, 0.0154682, 0.0214145, 0.00931383, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
88{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.550459, 0.233594, 0.101547, 0.0562414, 0.0201483, 0.0441865, 0.0198895, 0.0288849, 0.0246383, 0.00524203, 0.0111237, 0.0135135, 0.0100946, 0.0199063, 0.00630214, 0.0241457, 0.0250312, 0.0158502, 0.0166625, 0.0177019, 0.017145, 0.00993309, 0.00422611, -0.000272863, -0.000790177, 0.00521356, 0.00305057, 0.0112688, -0.000526226, 0.00204006, -0.000231335, 0.014503, 0.0218041, 0.0203275, 0.0139573, 0.0144748, 0.0111991, 0.0251907, 0.0389845, 0.0303032, 0.0321897, 0.031799, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
89{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.31804, 0.462231, 0.230675, 0.10903, 0.0760607, 0.0404051, 0.0469405, 0.0320198, 0.00642959, 0.0186602, 0.0069752, 0.0180583, 0.0216703, 0.00593377, 0.00838789, 0.0186795, 0.0195196, 0.0280292, 0.0110047, 0.02749, 0.0115594, 0.0194028, 0.0028292, 0.00615848, 0.00318525, 0.00250919, 0.00627185, 0.00337776, 0.00765032, 0.00137904, 0.00499126, 0.0173985, 0.0122038, 0.0200898, 0.010377, 0.00883389, 0.0141594, 0.0159646, 0.033311, 0.03605, 0.0156394, 0.021787, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
90{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.314286, 0.738615, 0.662827, 0.196162, 0.132289, 0.0522324, 0.0249344, 0.0432043, 0.0241318, 0.00922015, 0.0218616, 0.0470679, 0.0155021, 0.0192204, 0.0328014, 0.0237527, 0.00878217, 0.0294262, 0.00927312, 0.0127178, 0.010643, 0.0126027, -0.000958031, 0.00715873, -0.000525581, 0.00611947, 0.00244028, 0.00557385, 0, -0.000912307, 0.00229492, 0.00962786, 0.0181747, 0.0248611, 0.0220639, 0.00574051, 0.0186322, 0.0228958, 0.0378924, 0.0285511, 0.0253246, 0.0261138, 0.00980392, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
91{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.866924, 0.621371, 0.331145, 0.170723, 0.0885901, 0.054562, 0.0316905, 0.0346162, 0.0279412, 0.00855344, 0.0206852, 0.0335571, 0.021851, 0.0114343, 0.00607073, 0.0340026, 0.0200892, 0.00613316, 0.0195596, 0.0159075, 0.0114371, -0.00175157, 0.00792463, -0.000240245, 0.00602034, 0.00282451, -0.000593159, 0.00878438, 0.00568302, 0.00577865, 0.00250455, 0.00941139, 0.0103844, 0.0056916, 0.00526885, 0.0122546, 0.0346439, 0.00223482, 0.0225098, 0.0422859, 0.0220061, 0.0259502, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
92{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.529412, 0.260627, 0.599698, 0.414071, 0.132819, 0.0704768, 0.0638888, 0.0526834, 0.0382149, 0.0491349, 0.0288163, 0.0170256, 0.0197696, 0.00849553, 0.0151863, 0.00670989, 0.00349418, 0.00877875, 0.0235989, 0.0162238, 0.0116127, -0.000296778, 0.00244255, 0.00322352, 0.00938696, 0.013315, 0.00501653, 0.00206745, 0.00277726, 0.00785791, 0.00772425, -0.000502062, 0.0107861, 0.00246631, 0.00834297, 0.0134382, 0.0177804, 0.00680223, 0.0282039, 0.0205284, 0.0321745, 0.0154897, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
93{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0949668, 0.325462, 0.433882, 0.383331, 0.192033, 0.081956, 0.0278949, 0.0174577, 0.0202854, 0.0112868, 0.0110913, 0.0241937, 0.00639931, 0.0162666, 0.019322, 0.0099197, 0.0104203, 0.0268816, 0.0233483, 0.00888861, -0.000286063, -0.00054935, -0.000276576, 0.00286216, 0, 0.0123713, -0.000646685, 0.00520266, 0.00211541, 0.00237124, 0.00446584, 0.00659812, 0.00915503, 0.0110205, 0.0155078, 0.00416737, 0.00699636, 0.0151217, 0.0218096, 0.0387199, 0.0407843, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
94{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0919615, 0.270385, 0.632716, 0.277026, 0.169116, 0.0846112, 0.0515777, 0.0301623, 0.0352538, 0.019161, 0.0236401, 0.00992687, 0.00506818, 0.00745408, 0.0112244, 0.0104928, 0.017124, 0.0198827, 0.0092335, 0.00243917, 0.00311102, 0.00252178, 0.0161405, 0.00537439, 0.00931165, 0.00643444, 0.00593156, -0.000277926, 0.00828081, 0.00537097, -0.000638871, -0.000558436, -0.000523722, 0.00303963, 0.0139373, 0.0106944, 0.0245119, 0.0257483, 0.0233867, 0.0281246, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
95{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0259181, 0.0647234, 0.167592, 0.504096, 0.270604, 0.0973994, 0.0802578, 0.0357504, 0.0338436, 0.0390779, 0.00730446, 0.0169072, 0.0120201, 0.00696885, 0.0127019, 0.016837, 0.0158419, 0.00808887, 0.00323731, -0.00117744, 0, 0.00249525, 0.00303478, -0.000257923, -0.000974006, 0.00653244, 0, 0.00468244, 0.0029235, -0.000765462, 0, -0.000502669, 0.00255199, 0.0135758, 0.0142914, 0.0154734, 0.0110644, 0.0220854, 0.0260194, 0.0158371, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
96{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.0420469, 0.067839, 0.155801, 0.408299, 0.311472, 0.0952121, 0.0854928, 0.0426954, 0.0472351, 0.0363388, 0.0189654, 0.0137752, 0.0217982, -0.000864619, 0.0144195, 0.0172023, 0.0151226, 0.00949449, -0.000647373, 0, 0.00298522, 0.00575406, 0.00237401, 0.00563414, 0.00334244, 0.00271553, 0, -0.000551495, 0.00267126, 0.00279805, -0.000280827, 0.00268248, 0.00539638, 0.0201354, 0.00715991, 0.0147071, 0.0256984, 0.0241801, 0.0136786, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
97{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0255061, 0.0600916, 0.0632555, 0.146276, 0.351375, 0.280802, 0.150239, 0.0518022, 0.0308214, 0.0264112, 0.0237854, 0.016946, 0.00236536, 0.0124264, 0.0264236, 0.0158422, 0.0122726, 0.00495155, 0.0108897, -0.000293101, 0.0110721, 0, 0.00324573, 0, 0.00332531, 0, 0.00671498, 0.00251135, 0.00713558, 0, 0.00466496, -0.00047406, -0.000241937, -0.0005298, 0.00350627, 0.0162482, 0.0166161, 0.0060246, 0.0292655, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
98{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00643822, 0.0460297, 0.0527178, 0.168979, 0.324375, 0.235858, 0.111768, 0.0733801, 0.0358306, 0.0175744, 0.0250006, 0.0263039, 0.0366774, 0.00591692, 0.00766176, 0.00540716, 0.00643774, 0.0022853, 0.00437708, 0.00735904, 0.00580082, 0.00676819, 0.00281819, 0.00883046, 0.00760768, 0.0102608, -0.000280281, 0.00270775, -0.000656622, 0, 0.00899496, -0.000480025, 0.00222672, 0.00129486, 0.00238721, 0.0163661, 0.0212252, -0.00148296, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
99{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.019047, 0.0229463, 0.0427084, 0.154597, 0.242385, 0.191279, 0.124368, 0.0572722, 0.0366272, 0.0261736, 0.0222449, 0.0133094, 0.0101925, 0.0132689, 0.00895261, 0.0092906, 0.0105494, 0, -0.00102972, -0.000270034, 0.00274342, 0.00830494, 0.00611775, 0.00804174, 0.00575465, 0.00630444, 0.00699402, 0.00338315, 0.00710927, 0.00627783, 0.00231599, 0.00250542, 0.00444965, 0.0093392, 0.0107216, 0.0179824, 0.0421607, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
100{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.01376, 0.0220924, 0.0238798, 0.0713348, 0.102667, 0.25276, 0.146653, 0.0754997, 0.0322964, 0.0121763, 0.0261575, 0.0147211, 0.00259954, 0.0218871, 0.00233333, 0.00614367, 0.00498323, 0.00606678, 0.00539255, 0, 0, 0.0031289, -0.000529586, 0.00924435, 0.00868669, 0.00230543, 0, 0, -0.000268991, 0, -0.000269554, 0.00161804, -0.00189121, 0.00197102, 0.0115544, 0.0155885, 0.0180045, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
101{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0266696, 0.019719, 0.0227988, 0.0943135, 0.204843, 0.110036, 0.0401929, 0.0260484, 0.0296948, 0.0208731, 0.00669238, 0.00968895, 0.00669996, 0.0126931, 0.00205444, 0.00231361, 0.00334669, 0, -0.000292363, 0.00992995, 0.00207278, 0, 0.00525343, 0.00260067, 0.00531483, -0.000685799, 0.00283505, -0.000259943, 0.00536854, 0.000769969, 0.00152618, 0.00285896, 0.00530952, 0, 0.00868056, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
102{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0246648, 0.0238146, 0.0100666, 0.0481742, 0.0599759, 0.14597, 0.100054, 0.041173, 0.0395047, 0.0115925, 0.012527, 0.00579226, 0.00300786, -0.000953928, 0.00268132, 0.00194135, 0.00261804, 0.00576977, 0.00323769, 0.0081176, 0.00861065, 0.00683228, 0.00549571, 0.00240418, 0.00600316, 0.00249105, 0.00262977, 0, 0, 0.0128301, -0.00121571, 0.00290027, 0, 0.00615531, 0.0186773, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
103{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0295455, 0.0172865, 0.0162673, 0.0264695, 0.0273465, 0.0769191, 0.116368, 0.0934449, 0.0297024, 0.0238043, 0.015884, 0.0061, 0.0122739, 0.00321439, 0.00464015, 0.00506923, 0.00255592, 0.00598962, 0.00824517, 0.00821904, 0.0124875, 0.0133367, -0.000267685, -0.000245579, 0.00594068, 0.00254237, 0.00555178, 0.00562552, 0.00558605, 0.0146251, 0.00771492, 0.00589714, 0, -0.000250171, -0.000914967, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
104{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.000914518, 0.0189659, 0.0155641, 0.0219774, 0.0130828, 0.0218267, 0.0992131, 0.0467706, 0.0171446, 0.0124917, 0.00219422, 0, 0.0027203, 0.00634067, 0.00814849, 0.0132899, 0.00568592, 0.00356134, 0.0108418, 0.00297389, 0.00340212, 0.0044349, 0.0148501, 0.00837351, 0.00270315, 0.00824522, -0.000253934, 0.00233866, -0.000824389, 0.0161455, 0.00250434, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
105{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.0356778, 0.00398527, 0.0401558, 0.0153687, 0.015565, 0.0397546, 0.0427667, 0.0784397, 0.0326612, 0.0237238, 0.00624713, 0.0116745, 0.00980429, 0.00580568, 0.00627916, 0.00585718, 0.00311595, 0.0029994, 0, 0.00266769, 0.00366819, 0.00494056, 0.00985915, 0.0125486, 0.00500476, 0.00263505, 0.00259518, 0.00250188, 0.0065415, 0.0166975, 0.00844484, 0.0132855, -0.00030687, 0.00939542, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
106{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.01593, 0.00672758, 0.00930593, 0.0181723, 0.0166542, 0.0446327, 0.0515939, 0.0439146, 0.00692641, 0.0031473, 0.0062545, 0.00836715, 0.00652286, 0.00265167, 0.0091786, 0.00655072, 0.00918497, 0.0108358, 0.00893231, 0.00460159, 0.00276126, 0.0156746, 0.00766189, 0.00198682, 0.00603986, 0.00766474, 0.00150904, 0.00392221, 0.00651387, 0.00561985, 0.00636521, -0.00204981, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
107{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0648727, 0.00644001, 0.00980392, 0.00739817, 0.0219095, 0.0144736, 0.0291128, 0.0567089, 0.0174434, 0.00993839, 0.00291648, 0.00861031, -0.000543764, 0.00266136, 0.0111344, 0.0161496, 0.00851912, 0.0102413, 0.0130172, 0.00646307, 0.0125816, 0.00825952, 0.00248001, 0.00704836, 0.0116331, 0.00686234, 0.00521266, 0.0112, 0.00863957, -0.000244552, 0.00268468, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
108{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0150545, 0.0175218, 0.00277862, 0.0138633, 0.0070406, 0.00711215, 0.0224095, 0.0227712, 0.00672862, 0.0097044, 0.0128987, -0.000535921, 0.00316704, 0.00853934, 0.00567407, 0.0145124, 0.00705086, 0.00526955, 0.013757, 0.00925532, 0.00476223, 0.0141427, 0.0197044, 0.00620172, 0.0222063, 0.00138602, 0.010233, 0.00566005, 0.00223222, 0.00273963, 0.014265, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
109{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.000435152, 0.0142117, 0.00315209, 0.00972522, 0.0159469, 0.0135021, 0.0122938, 0.00416268, 0.00205024, -0.000512538, 0.00735074, -0.000594068, 0.0134428, 0.00690039, -0.000648465, 0.00609985, 0.0113877, 0.0050996, 0.0172142, 0.00328983, 0.0191053, 0.0102635, 0.00913039, 0.0118038, 0.00669382, 0.0179017, 0.012209, 0.00859542, 0.00620819, -0.000851349, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
110{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0170614, 0.0225141, 0.020126, 0.00624899, 0.00962921, 0.00250487, 0.00543442, 0.0029568, 0.00608116, 0.00236459, 0.00394319, 0.00310485, 0.00762142, 0.00467371, 0.00315239, 0.0107377, 0.0114243, 0.012687, 0.00730633, 0.0116869, 0.0105171, 0.00830238, 0.00482673, 0.0230868, 0.0141066, 0.00602666, 0.0120504, -0.00052991, 0.00684932, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
111{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0127806, 0.00514777, 0.00290073, 0.0160417, 0.0117424, 0.00322849, 0.0072234, 0.0311349, 0.0113657, 0.00888548, 0.00232162, 0, 0.0108722, 0.0031223, 0.00718234, 0.0131875, 0.0158046, 0.0120447, 0.0217175, 0.0103309, 0.0138622, 0.00955613, 0.00934435, 0.0117761, 0.00699558, 0.0158905, 0.0031348, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
112{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0181241, 0.0197183, 0.00606099, 0.00699557, 0.00614893, 0.00339383, 0.0054746, 0.0140354, 0.0214791, 0.0147735, 0.00335442, 0.00566671, 0.00564051, 0.00894947, 0.0133605, 0.0110452, 0.0102459, 0.0127268, 0.0165319, 0.0159548, 0.0059625, 0.00533616, 0.0145233, 0.0115926, 0.0128668, 0.0027315, 0.00388236, -0.00168645, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
113{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00295187, 0.00502366, 0.0125663, 0.00751157, 0.00678689, -0.000729155, 0.0398975, 0.0220088, 0.0360997, 0.0047205, 0.00260361, 0.00275378, 0.0125894, 0.00771866, 0.021691, 0.0217528, 0.0163468, 0.0112058, 0.00790909, 0.0146529, 0.00240104, 0.0103368, 0.00418216, 0.0128487, 0, 0.0125642, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
114{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0199918, 0.0154634, 0.0112614, 0.00289308, 0.00227097, 0.0148894, 0.0107923, 0.0619883, 0.0597212, 0.0187601, 0.00166553, 0.0059313, 0.0129129, 0.0343336, 0.0117161, 0.0319808, 0.0171083, 0.0215434, 0.0158661, 0.0074898, 0.00682951, 0.00862445, 0.0088061, 0.0123945, 0.00318738, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
115{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0188647, 0.00990385, 0.0151969, 0.00620634, 0.00309458, 0.0108622, 0.022087, 0.0626878, 0.0621444, 0.0309275, 0.0100102, 0.011632, 0.0104987, 0.027473, 0.0231617, 0.0227292, 0.0276311, 0.0215664, 0.0182076, 0.0101166, 0.0188944, 0.0191069, 0.00272097, 0.00354929, -0.00209698, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
116{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.020877, 0.00281841, 0.010183, 0.0132576, 0.00950378, 0.0105721, 0.0305541, 0.091212, 0.151017, 0.0258786, 0.0123632, 0.00293484, 0.0302581, 0.023202, 0.0305653, 0.0317708, 0.0223734, 0.0284086, 0.00786618, 0.0120922, 0.025366, 0.0150185, 0.00425779, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
117{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0110763, 0.00241807, 0.0133283, 0.00693999, 0.0199279, 0.0415627, 0.0550772, 0.117899, 0.130702, 0.0416414, 0.00755489, 0.0201174, 0.0244379, 0.0150149, 0.0460953, 0.0237191, 0.0230652, 0.0279088, 0.0150962, 0.0196312, 0.0163489, 0.0077193, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
118{0, 0, 0, 0, 0, 0, 0, 0, -0.0229358, 0.0107774, 0.00192411, 0.00503568, 0.0202753, 0.0396458, 0.0419664, 0.0596123, 0.128126, 0.141565, 0.0333455, 0.0352558, 0.0500001, 0.0260696, 0.0265925, 0.0430366, 0.0237963, 0.0271992, 0.0199606, 0.00803824, 0.0214636, 0.0238036, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
119{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00687079, 0.00808095, 0.0183833, 0.0186556, 0.0482115, 0.0551489, 0.101263, 0.128043, 0.125156, 0.0546276, 0.0597452, 0.0397317, 0.0664489, 0.0365752, 0.0333781, 0.0179711, 0.0215499, 0.0281107, 0.015737, 0.00394378, -0.00222857, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
120{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00559818, 0.00353283, 0.0125729, 0.0573947, 0.0357342, 0.086303, 0.111298, 0.154597, 0.0761374, 0.112503, 0.0954405, 0.05911, 0.0811331, 0.0442798, 0.0295957, 0.0564397, 0.0358336, 0.0169744, 0.00895756, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
121{0, 0, 0, 0, 0, 0, 0, 0, 0.0326923, 0.0065086, 0.00936733, 0.0180777, 0.0361944, 0.0903415, 0.151995, 0.128936, 0.121159, 0.11657, 0.158628, 0.129911, 0.0869573, 0.0374896, 0.0216581, 0.0520988, 0.0205289, 0.0234075, 0.013839, 0.0555556, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
122{0, 0, 0, 0, 0, 0, 0, 0, 0.0293347, 0.0139685, 0.0226661, 0.0552069, 0.0424837, 0.131329, 0.109362, 0.135511, 0.0612981, 0.247144, 0.169952, 0.0947945, 0.106144, 0.0524266, 0.0604242, 0.0434254, 0.0286301, 0.0588995, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
123{0, 0, 0, 0, 0, 0, 0, 0, 0.0137009, 0.0202961, 0.0288702, 0.061196, 0.0683367, 0.110711, 0.0701039, 0.056112, 0.066703, 0.185605, 0.241631, 0.13703, 0.128952, 0.0739922, 0.0572765, 0.0394537, 0.0574397, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
124{0, 0, 0, 0, 0, 0, 0, 0, 0.05228, 0.0312189, 0.0774895, 0.0928928, 0.135765, 0.105098, 0.0492366, 0.0210353, 0.0976445, 0.26425, 0.308278, 0.198294, 0.101087, 0.0561454, 0.0810017, 0.047073, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
125{0, 0, 0, 0, 0, 0, 0, 0, 0.0430778, 0.0421065, 0.093356, 0.0978211, 0.0825034, 0.0275584, 0.0275763, 0.0319886, 0.120972, 0.358682, 0.350457, 0.174366, 0.130494, 0.0821672, 0.0690678, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
126{0, 0, 0, 0, 0, 0, 0, 0, 0.0511346, 0.104765, 0.100829, 0.0613283, 0.0682207, 0.0378536, 0.0281379, 0.0337755, 0.159924, 0.515976, 0.474086, 0.232897, 0.175042, 0.0142627, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
127{0, 0, 0, 0, 0, 0, 0, 0, 0.0868946, 0.138176, 0.0721574, 0.0682048, 0.0487367, 0.0130789, 0.00963483, 0.0317815, 0.231897, 0.469518, 0.454897, 0.21412, 0.137255, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
128{0, 0, 0, 0, 0, 0, 0, 0, 0.143943, 0.141874, 0.0501214, 0.0340479, 0.0270368, 0.0189126, 0.0252699, 0.0405961, 0.185569, 0.682238, 0.491283, 0.513507, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
129{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0671399, 0.0558913, 0.026482, 0.00721053, 0.00775069, 0.0107691, 0.0774007, 0.234518, 0.71368, 1.89744, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
130{0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0207169, 0.0349959, 0.0114188, 0, 0.00752156, 0.041949, 0.0376233, 0.146694, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
131{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.00948451, 0.00767729, 0, 0.0089524, 0.0575039, 0.0570246, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
132{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0697674, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
133{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
134{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
135{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
136{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
137};
138
139for(int i=0;i<72;i++){
140 for (int j=0; j<72; j++) {
141 HisPDF[i][j] = HisPDFtmp[i][j];
142 }
143}
144
145}
146
147void EvtD0toKSpipiPlot::initProbMax(){
148 noProbMax();
149}
150
151void EvtD0toKSpipiPlot::decay( EvtParticle *p ){
152
153loop:
154 p->initializePhaseSpace(getNDaug(),getDaugs());
155
156 EvtParticle *id1,*id2,*id3;
157 EvtVector4R pd1,pd2,pd3;
158 double xmass13,xmass12, xmass23;
159
160 id1 =p->getDaug(0);
161 id2 =p->getDaug(1);
162 id3 =p->getDaug(2);
163
164 pd1 =id1->getP4Lab();
165 pd2 =id2->getP4Lab();
166 pd3 =id3->getP4Lab();
167
168 xmass12=(pd1+pd2).mass()*(pd1+pd2).mass(); // M_ksopi
169 xmass13=(pd1+pd3).mass()*(pd1+pd3).mass(); // M_ksopi
170// xmass23=(pd2+pd3).mass()*(pd2+pd3).mass(); // M_pipi
171
172 int xbin = FindXBin(xmass12);
173 int ybin = FindYBin(xmass13);
174 double xratio12=HisPDF[xbin][ybin]/avm1;
175
176 if(xratio12 <=0) goto loop;
177
178 double rd12=EvtRandom::Flat(0.0, 1.0);
179 if(rd12>xratio12) goto loop;
180
181 return ;
182}
183
184int EvtD0toKSpipiPlot::FindXBin(double mass2) {
185 if (mass2 < Xmin) { return 0; }
186 else if (mass2>=Xmax) { return 71; }
187 else { return int((mass2-Xmin)/Xwid)+1; }
188}
189
190int EvtD0toKSpipiPlot::FindYBin(double mass2) {
191 if (mass2 < Ymin) { return 0; }
192 else if (mass2>=Ymax) { return 71; }
193 else { return int((mass2-Ymin)/Ywid)+1; }
194}
mass
double mass
Definition: CosmicGenerator.cxx:138
EvtD0toKSpipiPlot::decay
void decay(EvtParticle *p)
Definition: EvtD0toKSpipiPlot.cc:151
EvtD0toKSpipiPlot::getName
void getName(std::string &name)
Definition: EvtD0toKSpipiPlot.cc:38
EvtD0toKSpipiPlot::FindYBin
int FindYBin(double mass2)
Definition: EvtD0toKSpipiPlot.cc:190
EvtD0toKSpipiPlot::~EvtD0toKSpipiPlot
virtual ~EvtD0toKSpipiPlot()
Definition: EvtD0toKSpipiPlot.cc:36
EvtD0toKSpipiPlot::FindXBin
int FindXBin(double mass2)
Definition: EvtD0toKSpipiPlot.cc:184
EvtD0toKSpipiPlot::clone
EvtDecayBase * clone()
Definition: EvtD0toKSpipiPlot.cc:42
EvtDecayBase::noProbMax
void noProbMax()
Definition: EvtDecayBase.cc:304
EvtDecayBase::getParentId
EvtId getParentId()
Definition: EvtDecayBase.hh:60
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition: EvtDecayBase.cc:482
EvtPDL::getSpinType
static EvtSpinType::spintype getSpinType(EvtId i)
Definition: EvtPDL.hh:61
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
EvtParticle::getP4Lab
EvtVector4R getP4Lab()
Definition: EvtParticle.cc:684
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:84
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1070
EvtRandom::Flat
static double Flat()
Definition: EvtRandom.cc:74