EvtDToKmPipPipPi0.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: EvtDToKmPipPipPi0.cc
//         the necessary file: EvtDToKmPipPipPi0.hh
//
// Description: D+ -> K- pi+ pi+ pi0 (BAM-00739)
//
// Modification history:
//
//    Liaoyuan Dong    Sat Dec 31 13:54:49 2022    Module created
//                     Fri Mar  1 23:43:09 2024    Module updated
//------------------------------------------------------------------------
#include "EvtGenBase/EvtPatches.hh"
#include <fstream>
#include <stdlib.h>
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtGenKine.hh"
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtReport.hh"
#include "EvtGenBase/EvtResonance.hh"
#include "EvtGenBase/EvtResonance2.hh"
#include "EvtGenModels/EvtDToKmPipPipPi0.hh"
#include <string>
#include "EvtGenBase/EvtConst.hh"
#include "EvtGenBase/EvtFlatte.hh"
#include "EvtGenBase/EvtDecayTable.hh"
#include "TMatrix.h"
#include "TMatrixD.h"
#include "TComplex.h"
//#include "width.h"
 
using std::endl;
 
EvtDToKmPipPipPi0::~EvtDToKmPipPipPi0() {}
 
void EvtDToKmPipPipPi0::getName(std::string& model_name){
    model_name="DToKmPipPipPi0";
}
 
EvtDecayBase* EvtDToKmPipPipPi0::clone(){
    return new EvtDToKmPipPipPi0;
}
 
double a1_width[]={
    0.0000000172, 0.0000000487, 0.0000001049, 0.0000001927, 0.0000003185, 
    0.0000004888, 0.0000007099, 0.0000009878, 0.0000013285, 0.0000017378, 
    0.0000022212, 0.0000027845, 0.0000034328, 0.0000041717, 0.0000050061, 
    0.0000059414, 0.0000069824, 0.0000081342, 0.0000094014, 0.0000107890, 
    0.0000123016, 0.0000139437, 0.0000157200, 0.0000176350, 0.0000196931, 
    0.0000218986, 0.0000242559, 0.0000267692, 0.0000294428, 0.0000322808, 
    0.0000352875, 0.0000384668, 0.0000418228, 0.0000453596, 0.0000490811, 
    0.0000529913, 0.0000570941, 0.0000613934, 0.0000658931, 0.0000705969, 
    0.0000755087, 0.0000806324, 0.0000859716, 0.0000915300, 0.0000973115, 
    0.0001033197, 0.0001095583, 0.0001160310, 0.0001227415, 0.0001296933, 
    0.0001368901, 0.0001443355, 0.0001520332, 0.0001599867, 0.0001681995, 
    0.0001766753, 0.0001854177, 0.0001944301, 0.0002037162, 0.0002132795, 
    0.0002231235, 0.0002332517, 0.0002436677, 0.0002543750, 0.0002653771, 
    0.0002766775, 0.0002882799, 0.0003001876, 0.0003124041, 0.0003249331, 
    0.0003377780, 0.0003509423, 0.0003644296, 0.0003782433, 0.0003923870, 
    0.0004068643, 0.0004216785, 0.0004368333, 0.0004523323, 0.0004681788, 
    0.0004843765, 0.0005009290, 0.0005178397, 0.0005351122, 0.0005527502, 
    0.0005707571, 0.0005891365, 0.0006078922, 0.0006270275, 0.0006465462, 
    0.0006664518, 0.0006867480, 0.0007074384, 0.0007285267, 0.0007500165, 
    0.0007719114, 0.0007942151, 0.0008169314, 0.0008400639, 0.0008636163, 
    0.0008875923, 0.0009119957, 0.0009368303, 0.0009620997, 0.0009878078, 
    0.0010139580, 0.0010405550, 0.0010676020, 0.0010951030, 0.0011230610, 
    0.0011514820, 0.0011803670, 0.0012097230, 0.0012395510, 0.0012698570, 
    0.0013006440, 0.0013319160, 0.0013636770, 0.0013959320, 0.0014286840, 
    0.0014619370, 0.0014956950, 0.0015299630, 0.0015647450, 0.0016000440, 
    0.0016358650, 0.0016722130, 0.0017090900, 0.0017465020, 0.0017844530, 
    0.0018229470, 0.0018619890, 0.0019015820, 0.0019417310, 0.0019824400, 
    0.0020237150, 0.0020655590, 0.0021079760, 0.0021509730, 0.0021945510, 
    0.0022387180, 0.0022834760, 0.0023288310, 0.0023747880, 0.0024213500, 
    0.0024685230, 0.0025163120, 0.0025647210, 0.0026137560, 0.0026634200, 
    0.0027137200, 0.0027646590, 0.0028162430, 0.0028684780, 0.0029213670, 
    0.0029749170, 0.0030291320, 0.0030840170, 0.0031395780, 0.0031958200, 
    0.0032527490, 0.0033103700, 0.0033686880, 0.0034277080, 0.0034874370, 
    0.0035478800, 0.0036090420, 0.0036709290, 0.0037335480, 0.0037969030, 
    0.0038610000, 0.0039258460, 0.0039914470, 0.0040578080, 0.0041249360, 
    0.0041928370, 0.0042615160, 0.0043309810, 0.0044012380, 0.0044722930, 
    0.0045441520, 0.0046168220, 0.0046903100, 0.0047646230, 0.0048397670, 
    0.0049157480, 0.0049925750, 0.0050702540, 0.0051487920, 0.0052281970, 
    0.0053084750, 0.0053896340, 0.0054716810, 0.0055546240, 0.0056384710, 
    0.0057232280, 0.0058089050, 0.0058955090, 0.0059830470, 0.0060715290, 
    0.0061609620, 0.0062513540, 0.0063427140, 0.0064350500, 0.0065283720, 
    0.0066226870, 0.0067180050, 0.0068143350, 0.0069116850, 0.0070100650, 
    0.0071094850, 0.0072099530, 0.0073114800, 0.0074140740, 0.0075177470, 
    0.0076225070, 0.0077283660, 0.0078353330, 0.0079434180, 0.0080526330, 
    0.0081629880, 0.0082744940, 0.0083871620, 0.0085010030, 0.0086160290, 
    0.0087322510, 0.0088496800, 0.0089683290, 0.0090882100, 0.0092093350, 
    0.0093317160, 0.0094553660, 0.0095802980, 0.0097065250, 0.0098340590, 
    0.0099629140, 0.0100931000, 0.0102246400, 0.0103575500, 0.0104918200, 
    0.0106274900, 0.0107645700, 0.0109030700, 0.0110430000, 0.0111843800, 
    0.0113272300, 0.0114715700, 0.0116174000, 0.0117647500, 0.0119136300, 
    0.0120640600, 0.0122160600, 0.0123696300, 0.0125248100, 0.0126816100, 
    0.0128400500, 0.0130001400, 0.0131619100, 0.0133253800, 0.0134905500, 
    0.0136574600, 0.0138261300, 0.0139965600, 0.0141688000, 0.0143428500, 
    0.0145187400, 0.0146964900, 0.0148761200, 0.0150576500, 0.0152411200, 
    0.0154265300, 0.0156139200, 0.0158033100, 0.0159947200, 0.0161881800, 
    0.0163837200, 0.0165813500, 0.0167811100, 0.0169830300, 0.0171871200, 
    0.0173934200, 0.0176019600, 0.0178127600, 0.0180258500, 0.0182412700, 
    0.0184590400, 0.0186791900, 0.0189017500, 0.0191267600, 0.0193542500, 
    0.0195842500, 0.0198167900, 0.0200519000, 0.0202896300, 0.0205300000, 
    0.0207730500, 0.0210188200, 0.0212673300, 0.0215186400, 0.0217727700, 
    0.0220297700, 0.0222896700, 0.0225525200, 0.0228183400, 0.0230871900, 
    0.0233591000, 0.0236341200, 0.0239122800, 0.0241936400, 0.0244782300, 
    0.0247661000, 0.0250572900, 0.0253518600, 0.0256498400, 0.0259512800, 
    0.0262562300, 0.0265647400, 0.0268768600, 0.0271926300, 0.0275121100, 
    0.0278353500, 0.0281624000, 0.0284933100, 0.0288281300, 0.0291669200, 
    0.0295097300, 0.0298566200, 0.0302076400, 0.0305628500, 0.0309223000, 
    0.0312860500, 0.0316541600, 0.0320266900, 0.0324036900, 0.0327852300, 
    0.0331713600, 0.0335621500, 0.0339576500, 0.0343579300, 0.0347630500, 
    0.0351730600, 0.0355880400, 0.0360080500, 0.0364331400, 0.0368633900, 
    0.0372988500, 0.0377395900, 0.0381856600, 0.0386371500, 0.0390941000, 
    0.0395565900, 0.0400246700, 0.0404984100, 0.0409778800, 0.0414631300, 
    0.0419542300, 0.0424512400, 0.0429542200, 0.0434632400, 0.0439783500, 
    0.0444996200, 0.0450271000, 0.0455608500, 0.0461009400, 0.0466474000, 
    0.0472003200, 0.0477597200, 0.0483256700, 0.0488982300, 0.0494774300, 
    0.0500633300, 0.0506559700, 0.0512554000, 0.0518616600, 0.0524748000, 
    0.0530948400, 0.0537218300, 0.0543558000, 0.0549967800, 0.0556448000, 
    0.0562998900, 0.0569620700, 0.0576313600, 0.0583077800, 0.0589913500, 
    0.0596820600, 0.0603799400, 0.0610849900, 0.0617972100, 0.0625166000, 
    0.0632431500, 0.0639768500, 0.0647176900, 0.0654656600, 0.0662207300, 
    0.0669828800, 0.0677520900, 0.0685283100, 0.0693115200, 0.0701016800, 
    0.0708987400, 0.0717026500, 0.0725133700, 0.0733308400, 0.0741550000, 
    0.0749858000, 0.0758231500, 0.0766670000, 0.0775172700, 0.0783738800, 
    0.0792367600, 0.0801058200, 0.0809809800, 0.0818621400, 0.0827492100, 
    0.0836421000, 0.0845407100, 0.0854449400, 0.0863546800, 0.0872698400, 
    0.0881902900, 0.0891159400, 0.0900466700, 0.0909823700, 0.0919229100, 
    0.0928682000, 0.0938181000, 0.0947725000, 0.0957312800, 0.0966943100, 
    0.0976614900, 0.0986326700, 0.0996077500, 0.1005866000, 0.1015691000, 
    0.1025551000, 0.1035445000, 0.1045372000, 0.1055331000, 0.1065320000, 
    0.1075338000, 0.1085384000, 0.1095456000, 0.1105555000, 0.1115677000, 
    0.1125823000, 0.1135991000, 0.1146180000, 0.1156389000, 0.1166616000, 
    0.1176862000, 0.1187123000, 0.1197401000, 0.1207692000, 0.1217997000, 
    0.1228315000, 0.1238644000, 0.1248983000, 0.1259332000, 0.1269690000, 
    0.1280055000, 0.1290427000, 0.1300805000, 0.1311188000, 0.1321575000, 
    0.1331966000, 0.1342359000, 0.1352754000, 0.1363150000, 0.1373546000, 
    0.1383942000, 0.1394337000, 0.1404730000, 0.1415120000, 0.1425508000, 
    0.1435891000, 0.1446271000, 0.1456645000, 0.1467014000, 0.1477377000, 
    0.1487733000, 0.1498082000, 0.1508423000, 0.1518756000, 0.1529081000, 
    0.1539396000, 0.1549702000, 0.1559998000, 0.1570284000, 0.1580558000, 
    0.1590822000, 0.1601074000, 0.1611314000, 0.1621541000, 0.1631757000, 
    0.1641959000, 0.1652148000, 0.1662323000, 0.1672484000, 0.1682632000, 
    0.1692765000, 0.1702883000, 0.1712986000, 0.1723075000, 0.1733147000, 
    0.1743205000, 0.1753246000, 0.1763271000, 0.1773281000, 0.1783273000, 
    0.1793250000, 0.1803209000, 0.1813152000, 0.1823078000, 0.1832986000, 
    0.1842877000, 0.1852751000, 0.1862607000, 0.1872445000, 0.1882266000, 
    0.1892068000, 0.1901853000, 0.1911619000, 0.1921368000, 0.1931098000, 
    0.1940809000, 0.1950503000, 0.1960177000, 0.1969834000, 0.1979471000, 
    0.1989090000, 0.1998690000, 0.2008271000, 0.2017834000, 0.2027378000, 
    0.2036903000, 0.2046408000, 0.2055895000, 0.2065363000, 0.2074812000, 
    0.2084242000, 0.2093653000, 0.2103045000, 0.2112418000, 0.2121772000, 
    0.2131106000, 0.2140422000, 0.2149719000, 0.2158996000, 0.2168255000, 
    0.2177494000, 0.2186715000, 0.2195916000, 0.2205098000, 0.2214262000, 
    0.2223406000, 0.2232532000, 0.2241638000, 0.2250726000, 0.2259795000, 
    0.2268844000, 0.2277875000, 0.2286888000, 0.2295881000, 0.2304856000, 
    0.2313812000, 0.2322749000, 0.2331668000, 0.2340568000, 0.2349449000, 
    0.2358312000, 0.2367157000, 0.2375983000, 0.2384790000, 0.2393580000, 
    0.2402351000, 0.2411103000, 0.2419837000, 0.2428554000, 0.2437252000, 
    0.2445931000, 0.2454593000, 0.2463237000, 0.2471863000, 0.2480471000, 
    0.2489061000, 0.2497633000, 0.2506187000, 0.2514724000, 0.2523243000, 
    0.2531744000, 0.2540228000, 0.2548694000, 0.2557143000, 0.2565574000, 
    0.2573988000, 0.2582385000, 0.2590764000, 0.2599126000, 0.2607471000, 
    0.2615799000, 0.2624110000, 0.2632404000, 0.2640681000, 0.2648941000, 
    0.2657184000, 0.2665410000, 0.2673620000, 0.2681813000, 0.2689989000, 
    0.2698149000, 0.2706292000, 0.2714419000, 0.2722529000, 0.2730623000, 
    0.2738701000, 0.2746762000, 0.2754807000, 0.2762836000, 0.2770849000, 
    0.2778846000, 0.2786827000, 0.2794792000, 0.2802741000, 0.2810675000, 
    0.2818592000, 0.2826494000, 0.2834380000, 0.2842251000, 0.2850106000, 
    0.2857946000, 0.2865770000, 0.2873578000, 0.2881372000, 0.2889150000, 
    0.2896913000, 0.2904661000, 0.2912393000, 0.2920111000, 0.2927813000, 
    0.2935501000, 0.2943173000, 0.2950831000, 0.2958474000, 0.2966102000, 
    0.2973715000, 0.2981314000, 0.2988898000, 0.2996468000, 0.3004023000, 
    0.3011564000, 0.3019090000, 0.3026602000, 0.3034099000, 0.3041583000, 
    0.3049052000, 0.3056507000, 0.3063948000, 0.3071374000, 0.3078787000, 
    0.3086186000, 0.3093571000, 0.3100942000, 0.3108299000, 0.3115643000, 
    0.3122973000, 0.3130289000, 0.3137591000, 0.3144880000, 0.3152155000, 
    0.3159417000, 0.3166666000, 0.3173901000, 0.3181122000, 0.3188331000, 
    0.3195526000, 0.3202708000, 0.3209877000, 0.3217033000, 0.3224175000, 
    0.3231305000, 0.3238421000, 0.3245525000, 0.3252616000, 0.3259694000, 
    0.3266759000, 0.3273811000, 0.3280851000, 0.3287878000, 0.3294893000, 
    0.3301894000, 0.3308884000, 0.3315861000, 0.3322825000, 0.3329777000, 
    0.3336716000, 0.3343644000, 0.3350559000, 0.3357461000, 0.3364352000, 
    0.3371230000, 0.3378096000, 0.3384951000, 0.3391793000, 0.3398623000, 
    0.3405441000, 0.3412247000, 0.3419042000, 0.3425824000, 0.3432595000, 
    0.3439354000, 0.3446101000, 0.3452837000, 0.3459561000, 0.3466273000, 
    0.3472974000, 0.3479663000, 0.3486341000, 0.3493008000, 0.3499663000, 
    0.3506306000, 0.3512939000, 0.3519560000, 0.3526169000, 0.3532768000, 
    0.3539355000, 0.3545931000, 0.3552497000, 0.3559051000, 0.3565594000, 
    0.3572126000, 0.3578647000, 0.3585157000, 0.3591656000, 0.3598144000, 
    0.3604622000, 0.3611089000, 0.3617545000, 0.3623990000, 0.3630425000, 
    0.3636849000, 0.3643263000, 0.3649665000, 0.3656058000, 0.3662440000, 
    0.3668811000, 0.3675172000, 0.3681523000, 0.3687863000, 0.3694193000, 
    0.3700512000, 0.3706822000, 0.3713121000, 0.3719409000, 0.3725688000, 
    0.3731957000, 0.3738215000, 0.3744464000, 0.3750702000, 0.3756930000, 
    0.3763149000, 0.3769357000, 0.3775556000, 0.3781745000, 0.3787924000, 
    0.3794093000, 0.3800252000, 0.3806401000, 0.3812541000, 0.3818671000, 
    0.3824792000, 0.3830902000, 0.3837004000, 0.3843095000, 0.3849178000, 
    0.3855250000, 0.3861313000, 0.3867367000, 0.3873411000, 0.3879446000, 
    0.3885472000, 0.3891488000, 0.3897495000, 0.3903492000, 0.3909481000, 
    0.3915460000, 0.3921430000, 0.3927391000, 0.3933343000, 0.3939285000, 
    0.3945219000, 0.3951143000, 0.3957059000, 0.3962966000, 0.3968863000, 
    0.3974752000, 0.3980631000, 0.3986502000, 0.3992364000, 0.3998218000, 
    0.4004062000, 0.4009898000, 0.4015725000, 0.4021543000, 0.4027352000, 
    0.4033153000, 0.4038946000, 0.4044729000, 0.4050504000, 0.4056271000, 
    0.4062029000, 0.4067778000, 0.4073519000, 0.4079252000, 0.4084976000, 
    0.4090692000, 0.4096399000, 0.4102098000, 0.4107789000, 0.4113472000, 
    0.4119146000, 0.4124812000, 0.4130470000, 0.4136119000, 0.4141760000, 
    0.4147394000, 0.4153019000, 0.4158636000, 0.4164245000, 0.4169846000, 
    0.4175439000, 0.4181024000, 0.4186601000, 0.4192170000, 0.4197731000, 
    0.4203284000, 0.4208830000, 0.4214367000, 0.4219897000, 0.4225419000, 
    0.4230933000, 0.4236439000, 0.4241938000, 0.4247429000, 0.4252912000, 
    0.4258388000, 0.4263856000, 0.4269317000, 0.4274769000, 0.4280215000, 
    0.4285652000, 0.4291083000, 0.4296505000, 0.4301921000, 0.4307329000, 
    0.4312729000, 0.4318122000, 0.4323508000, 0.4328886000, 0.4334257000, 
    0.4339621000, 0.4344977000, 0.4350326000, 0.4355668000, 0.4361003000, 
    0.4366331000, 0.4371651000, 0.4376964000, 0.4382270000, 0.4387569000, 
    0.4392861000, 0.4398146000, 0.4403424000, 0.4408695000, 0.4413959000, 
    0.4419216000, 0.4424466000, 0.4429709000, 0.4434945000, 0.4440174000, 
    0.4445396000, 0.4450612000, 0.4455821000, 0.4461023000, 0.4466218000, 
    0.4471407000, 0.4476589000, 0.4481764000, 0.4486932000, 0.4492094000, 
    0.4497249000, 0.4502398000, 0.4507540000, 0.4512675000, 0.4517804000, 
    0.4522926000, 0.4528042000, 0.4533152000, 0.4538255000, 0.4543351000, 
    0.4548442000, 0.4553525000, 0.4558603000, 0.4563674000, 0.4568739000, 
    0.4573797000, 0.4578849000, 0.4583895000, 0.4588935000, 0.4593969000, 
    0.4598996000, 0.4604018000, 0.4609033000, 0.4614042000, 0.4619045000, 
    0.4624042000, 0.4629032000, 0.4634017000, 0.4638996000, 0.4643969000, 
    0.4648936000, 0.4653897000, 0.4658852000, 0.4663802000, 0.4668745000, 
    0.4673683000, 0.4678614000, 0.4683541000, 0.4688461000, 0.4693375000, 
    0.4698284000, 0.4703188000, 0.4708085000, 0.4712977000, 0.4717863000, 
    0.4722744000, 0.4727619000, 0.4732489000, 0.4737353000, 0.4742212000, 
    0.4747065000, 0.4751913000, 0.4756756000, 0.4761593000, 0.4766424000, 
    0.4771251000, 0.4776072000, 0.4780888000, 0.4785699000, 0.4790504000, 
    0.4795304000, 0.4800099000, 0.4804889000, 0.4809674000, 0.4814454000, 
    0.4819229000, 0.4823999000, 0.4828764000, 0.4833524000, 0.4838279000, 
    0.4843029000, 0.4847774000, 0.4852515000, 0.4857250000, 0.4861981000, 
    0.4866708000, 0.4871429000, 0.4876146000, 0.4880858000, 0.4885566000, 
    0.4890269000, 0.4894968000, 0.4899662000, 0.4904351000, 0.4909037000, 
    0.4913718000, 0.4918394000, 0.4923067000, 0.4927735000, 0.4932398000, 
    0.4937058000, 0.4941713000, 0.4946365000, 0.4951012000, 0.4955655000, 
    0.4960295000, 0.4964930000, 0.4969562000, 0.4974189000, 0.4978813000, 
    0.4983434000, 0.4988050000, 0.4992663000, 0.4997272000, 0.5001878000, 
    0.5006480000, 0.5011079000, 0.5015674000, 0.5020266000, 0.5024855000, 
    0.5029440000, 0.5034023000, 0.5038602000, 0.5043178000, 0.5047751000, 
    0.5052321000, 0.5056889000, 0.5061453000, 0.5066015000, 0.5070574000, 
    0.5075131000, 0.5079685000, 0.5084237000, 0.5088786000, 0.5093333000, 
    0.5097878000, 0.5102420000, 0.5106961000, 0.5111499000, 0.5116036000, 
    0.5120571000, 0.5125104000, 0.5129636000, 0.5134166000, 0.5138694000, 
    0.5143221000, 0.5147747000, 0.5152272000, 0.5156796000, 0.5161319000, 
    0.5165842000, 0.5170363000, 0.5174884000, 0.5179405000, 0.5183925000, 
    0.5188445000, 0.5192965000, 0.5197485000, 0.5202005000, 0.5206526000, 
    0.5211048000, 0.5215570000, 0.5220092000, 0.5224616000, 0.5229141000, 
    0.5233668000, 0.5238195000, 0.5242725000, 0.5247257000, 0.5251790000, 
    0.5256326000, 0.5260864000, 0.5265405000, 0.5269949000, 0.5274496000, 
    0.5279047000, 0.5283601000, 0.5288159000, 0.5292721000, 0.5297287000, 
    0.5301858000, 0.5306434000, 0.5311015000, 0.5315601000, 0.5320194000, 
    0.5324792000, 0.5329396000, 0.5334008000, 0.5338626000, 0.5343252000, 
    0.5347885000, 0.5352526000, 0.5357176000, 0.5361834000, 0.5366502000, 
    0.5371179000, 0.5375866000, 0.5380563000, 0.5385271000, 0.5389989000, 
    0.5394720000, 0.5399462000, 0.5404217000, 0.5408984000, 0.5413764000, 
    0.5418558000, 0.5423367000, 0.5428189000, 0.5433026000, 0.5437879000, 
    0.5442747000, 0.5447631000, 0.5452532000, 0.5457450000, 0.5462385000, 
    0.5467337000, 0.5472307000, 0.5477296000, 0.5482303000, 0.5487328000, 
    0.5492373000, 0.5497437000, 0.5502520000, 0.5507623000, 0.5512745000, 
    0.5517887000, 0.5523048000, 0.5528230000, 0.5533430000, 0.5538651000, 
    0.5543890000, 0.5549149000, 0.5554426000, 0.5559722000, 0.5565036000, 
    0.5570367000, 0.5575716000, 0.5581082000, 0.5586464000, 0.5591862000, 
    0.5597274000, 0.5602702000, 0.5608143000, 0.5613597000, 0.5619063000, 
    0.5624541000, 0.5630030000, 0.5635529000, 0.5641037000, 0.5646553000, 
    0.5652078000, 0.5657609000, 0.5663146000, 0.5668689000, 0.5674236000, 
    0.5679787000, 0.5685341000, 0.5690898000, 0.5696456000, 0.5702015000, 
    0.5707575000, 0.5713134000, 0.5718692000, 0.5724249000, 0.5729803000, 
    0.5735355000, 0.5740903000, 0.5746448000, 0.5751989000, 0.5757525000, 
    0.5763056000, 0.5768581000, 0.5774101000, 0.5779615000, 0.5785121000, 
    0.5790622000, 0.5796114000, 0.5801600000, 0.5807078000, 0.5812548000, 
    0.5818009000, 0.5823462000, 0.5828907000, 0.5834343000, 0.5839770000, 
    0.5845188000, 0.5850597000, 0.5855996000, 0.5861386000, 0.5866767000, 
    0.5872137000, 0.5877498000, 0.5882850000, 0.5888191000, 0.5893522000, 
    0.5898844000, 0.5904155000, 0.5909457000, 0.5914748000, 0.5920029000, 
    0.5925300000, 0.5930560000, 0.5935811000, 0.5941051000, 0.5946282000, 
    0.5951502000, 0.5956711000, 0.5961911000, 0.5967100000, 0.5972280000, 
    0.5977449000, 0.5982608000, 0.5987757000, 0.5992896000, 0.5998025000, 
    0.6003144000, 0.6008253000, 0.6013352000, 0.6018441000, 0.6023520000, 
    0.6028589000, 0.6033649000, 0.6038699000, 0.6043739000, 0.6048769000, 
    0.6053790000, 0.6058801000, 0.6063803000, 0.6068795000, 0.6073777000, 
    0.6078750000, 0.6083714000, 0.6088668000, 0.6093613000, 0.6098549000, 
    0.6103476000, 0.6108393000, 0.6113301000, 0.6118201000, 0.6123091000, 
    0.6127972000, 0.6132844000, 0.6137707000, 0.6142562000, 0.6147407000, 
    0.6152244000, 0.6157072000, 0.6161892000, 0.6166703000, 0.6171505000, 
    0.6176298000, 0.6181084000, 0.6185860000, 0.6190629000, 0.6195388000, 
    0.6200140000, 0.6204883000, 0.6209618000, 0.6214345000, 0.6219064000, 
    0.6223774000, 0.6228476000, 0.6233171000, 0.6237857000, 0.6242536000, 
    0.6247206000, 0.6251869000, 0.6256523000, 0.6261170000, 0.6265809000, 
    0.6270441000, 0.6275064000, 0.6279681000, 0.6284289000, 0.6288890000, 
    0.6293483000, 0.6298069000, 0.6302647000, 0.6307218000, 0.6311782000, 
    0.6316338000, 0.6320887000, 0.6325428000, 0.6329963000, 0.6334490000, 
    0.6339009000, 0.6343522000, 0.6348028000, 0.6352526000, 0.6357018000, 
    0.6361502000, 0.6365979000, 0.6370450000, 0.6374913000, 0.6379370000, 
    0.6383820000, 0.6388262000, 0.6392699000, 0.6397128000, 0.6401550000, 
    0.6405966000, 0.6410375000, 0.6414778000, 0.6419174000, 0.6423563000, 
    0.6427946000, 0.6432322000, 0.6436692000, 0.6441055000, 0.6445412000, 
    0.6449762000, 0.6454106000, 0.6458443000, 0.6462774000, 0.6467099000, 
    0.6471418000, 0.6475730000, 0.6480036000, 0.6484336000, 0.6488629000, 
    0.6492917000, 0.6497198000, 0.6501473000, 0.6505742000, 0.6510005000, 
    0.6514262000, 0.6518513000, 0.6522758000, 0.6526997000, 0.6531230000, 
    0.6535457000, 0.6539678000, 0.6543893000, 0.6548103000, 0.6552306000, 
    0.6556504000, 0.6560696000, 0.6564882000, 0.6569062000, 0.6573237000, 
    0.6577406000, 0.6581569000, 0.6585727000, 0.6589879000, 0.6594025000, 
    0.6598166000, 0.6602301000, 0.6606431000, 0.6610555000, 0.6614673000, 
    0.6618787000, 0.6622894000, 0.6626996000, 0.6631093000, 0.6635184000, 
    0.6639270000, 0.6643350000, 0.6647425000, 0.6651495000, 0.6655560000, 
    0.6659619000, 0.6663672000, 0.6667721000, 0.6671764000, 0.6675802000, 
    0.6679835000, 0.6683863000, 0.6687885000, 0.6691902000, 0.6695914000, 
    0.6699921000, 0.6703923000, 0.6707920000, 0.6711912000, 0.6715898000, 
    0.6719880000, 0.6723856000, 0.6727828000, 0.6731794000, 0.6735756000, 
    0.6739712000, 0.6743664000, 0.6747610000, 0.6751552000, 0.6755489000, 
    0.6759421000, 0.6763348000, 0.6767270000, 0.6771187000, 0.6775100000, 
    0.6779008000, 0.6782910000, 0.6786809000, 0.6790702000, 0.6794590000, 
    0.6798474000, 0.6802353000, 0.6806228000, 0.6810097000, 0.6813962000, 
    0.6817823000, 0.6821678000, 0.6825530000, 0.6829376000, 0.6833218000, 
    0.6837055000, 0.6840888000, 0.6844716000, 0.6848539000, 0.6852358000, 
    0.6856172000, 0.6859982000, 0.6863788000, 0.6867588000, 0.6871385000, 
    0.6875177000, 0.6878964000, 0.6882747000, 0.6886526000, 0.6890300000, 
    0.6894070000, 0.6897835000, 0.6901596000, 0.6905353000, 0.6909105000, 
    0.6912853000, 0.6916596000, 0.6920336000, 0.6924071000, 0.6927801000, 
    0.6931528000, 0.6935250000, 0.6938967000, 0.6942681000, 0.6946390000, 
    0.6950095000, 0.6953796000, 0.6957493000, 0.6961185000, 0.6964874000, 
    0.6968558000, 0.6972238000, 0.6975914000, 0.6979585000, 0.6983253000, 
    0.6986916000, 0.6990576000, 0.6994231000, 0.6997882000, 0.7001529000, 
    0.7005172000, 0.7008811000, 0.7012446000, 0.7016077000, 0.7019704000, 
    0.7023326000, 0.7026945000, 0.7030560000, 0.7034171000, 0.7037778000, 
    0.7041381000, 0.7044980000, 0.7048575000, 0.7052166000, 0.7055753000, 
    0.7059336000, 0.7062916000, 0.7066491000, 0.7070063000, 0.7073630000, 
    0.7077194000, 0.7080754000, 0.7084310000, 0.7087863000, 0.7091411000, 
    0.7094956000, 0.7098497000, 0.7102034000, 0.7105567000, 0.7109096000, 
    0.7112622000, 0.7116144000, 0.7119662000, 0.7123177000, 0.7126687000, 
    0.7130194000, 0.7133697000, 0.7137197000, 0.7140693000, 0.7144185000, 
    0.7147673000, 0.7151158000, 0.7154639000, 0.7158117000, 0.7161591000, 
    0.7165061000, 0.7168527000, 0.7171990000, 0.7175450000, 0.7178905000, 
    0.7182357000, 0.7185806000, 0.7189251000, 0.7192692000, 0.7196130000, 
    0.7199564000, 0.7202995000, 0.7206422000, 0.7209845000, 0.7213266000, 
    0.7216682000, 0.7220095000, 0.7223505000, 0.7226911000, 0.7230313000, 
    0.7233712000, 0.7237108000, 0.7240500000, 0.7243889000, 0.7247274000, 
    0.7250656000, 0.7254034000, 0.7257409000, 0.7260780000, 0.7264148000, 
    0.7267513000, 0.7270874000, 0.7274232000, 0.7277587000, 0.7280938000, 
    0.7284286000, 0.7287630000, 0.7290971000, 0.7294309000, 0.7297643000, 
    0.7300974000, 0.7304302000, 0.7307626000, 0.7310947000, 0.7314265000, 
    0.7317579000, 0.7320891000, 0.7324198000, 0.7327503000, 0.7330804000, 
    0.7334102000, 0.7337397000, 0.7340689000, 0.7343977000, 0.7347262000, 
    0.7350544000, 0.7353823000, 0.7357098000, 0.7360370000, 0.7363639000, 
    0.7366905000, 0.7370168000, 0.7373427000, 0.7376683000, 0.7379936000, 
    0.7383186000, 0.7386433000, 0.7389677000, 0.7392917000, 0.7396155000, 
    0.7399389000, 0.7402620000, 0.7405848000, 0.7409073000, 0.7412294000, 
    0.7415513000, 0.7418728000, 0.7421941000, 0.7425150000, 0.7428357000, 
    0.7431560000, 0.7434760000, 0.7437957000, 0.7441151000, 0.7444342000, 
    0.7447530000, 0.7450715000, 0.7453897000, 0.7457076000, 0.7460251000, 
    0.7463424000, 0.7466594000, 0.7469761000, 0.7472925000, 0.7476085000, 
    0.7479243000, 0.7482398000, 0.7485550000, 0.7488699000, 0.7491845000, 
    0.7494988000, 0.7498128000, 0.7501265000, 0.7504399000, 0.7507531000, 
    0.7510659000, 0.7513784000, 0.7516907000, 0.7520027000, 0.7523143000, 
    0.7526257000, 0.7529368000, 0.7532476000, 0.7535581000, 0.7538683000, 
    0.7541783000, 0.7544879000, 0.7547973000, 0.7551064000, 0.7554152000, 
    0.7557237000, 0.7560319000, 0.7563398000, 0.7566475000, 0.7569549000, 
    0.7572620000, 0.7575688000, 0.7578753000, 0.7581816000, 0.7584875000, 
    0.7587932000, 0.7590986000, 0.7594038000, 0.7597086000, 0.7600132000, 
    0.7603175000, 0.7606215000, 0.7609253000, 0.7612288000, 0.7615320000, 
    0.7618349000, 0.7621376000, 0.7624399000, 0.7627420000, 0.7630439000, 
    0.7633454000, 0.7636467000, 0.7639477000, 0.7642485000, 0.7645490000, 
    0.7648492000, 0.7651491000, 0.7654488000, 0.7657482000, 0.7660473000, 
    0.7663462000, 0.7666448000, 0.7669431000, 0.7672412000, 0.7675390000, 
    0.7678366000, 0.7681338000, 0.7684309000, 0.7687276000, 0.7690241000, 
    0.7693203000, 0.7696163000, 0.7699120000, 0.7702074000, 0.7705026000, 
    0.7707975000, 0.7710922000, 0.7713866000, 0.7716807000, 0.7719746000, 
    0.7722682000, 0.7725616000, 0.7728547000, 0.7731475000, 0.7734401000, 
    0.7737325000, 0.7740245000, 0.7743164000, 0.7746079000, 0.7748993000, 
    0.7751903000, 0.7754811000, 0.7757717000, 0.7760620000, 0.7763521000, 
    0.7766419000, 0.7769314000, 0.7772207000, 0.7775098000, 0.7777986000, 
    0.7780871000, 0.7783754000, 0.7786635000, 0.7789513000, 0.7792389000, 
    0.7795262000, 0.7798132000, 0.7801000000, 0.7803866000, 0.7806729000, 
    0.7809590000, 0.7812449000, 0.7815304000, 0.7818158000, 0.7821009000, 
    0.7823858000, 0.7826704000, 0.7829547000, 0.7832389000, 0.7835228000, 
    0.7838064000, 0.7840898000, 0.7843730000, 0.7846559000, 0.7849386000, 
    0.7852210000, 0.7855032000, 0.7857852000, 0.7860669000, 0.7863484000, 
    0.7866297000, 0.7869107000, 0.7871915000, 0.7874720000, 0.7877523000, 
    0.7880324000, 0.7883122000, 0.7885919000, 0.7888712000, 0.7891504000, 
    0.7894293000, 0.7897079000, 0.7899864000, 0.7902646000, 0.7905425000, 
    0.7908203000, 0.7910978000, 0.7913750000, 0.7916521000, 0.7919289000, 
    0.7922055000, 0.7924818000, 0.7927579000, 0.7930338000, 0.7933095000, 
    0.7935849000, 0.7938601000, 0.7941351000, 0.7944099000, 0.7946844000, 
    0.7949587000, 0.7952328000, 0.7955066000, 0.7957803000, 0.7960537000, 
    0.7963268000, 0.7965998000, 0.7968725000, 0.7971450000, 0.7974173000, 
    0.7976893000, 0.7979612000, 0.7982328000, 0.7985041000, 0.7987753000, 
    0.7990463000, 0.7993170000, 0.7995875000, 0.7998578000, 0.8001278000, 
    0.8003977000, 0.8006673000, 0.8009367000, 0.8012059000, 0.8014748000, 
    0.8017436000, 0.8020121000, 0.8022804000, 0.8025485000, 0.8028164000, 
    0.8030841000, 0.8033515000, 0.8036187000, 0.8038857000, 0.8041525000, 
    0.8044191000, 0.8046855000, 0.8049517000, 0.8052176000, 0.8054833000, 
    0.8057488000, 0.8060141000, 0.8062792000, 0.8065441000, 0.8068088000, 
    0.8070732000, 0.8073375000, 0.8076015000, 0.8078653000, 0.8081289000, 
    0.8083923000, 0.8086555000, 0.8089185000, 0.8091813000, 0.8094438000, 
    0.8097062000, 0.8099683000, 0.8102303000, 0.8104920000, 0.8107535000, 
    0.8110148000, 0.8112760000, 0.8115369000, 0.8117976000, 0.8120580000, 
    0.8123183000, 0.8125784000, 0.8128383000, 0.8130980000, 0.8133574000, 
    0.8136167000, 0.8138757000, 0.8141346000, 0.8143933000, 0.8146517000, 
    0.8149099000, 0.8151680000, 0.8154258000, 0.8156835000, 0.8159409000, 
    0.8161982000, 0.8164552000, 0.8167120000, 0.8169687000, 0.8172251000, 
    0.8174813000, 0.8177374000, 0.8179932000, 0.8182488000, 0.8185043000, 
    0.8187595000, 0.8190146000, 0.8192694000, 0.8195241000, 0.8197785000, 
    0.8200328000, 0.8202868000, 0.8205407000, 0.8207944000, 0.8210478000, 
    0.8213011000, 0.8215542000, 0.8218071000, 0.8220598000, 0.8223122000, 
    0.8225645000, 0.8228167000, 0.8230686000, 0.8233203000, 0.8235718000, 
    0.8238232000, 0.8240743000, 0.8243252000, 0.8245760000, 0.8248266000, 
    0.8250769000, 0.8253271000, 0.8255771000, 0.8258269000, 0.8260765000, 
    0.8263259000, 0.8265752000, 0.8268242000, 0.8270731000, 0.8273217000, 
    0.8275702000, 0.8278185000, 0.8280666000, 0.8283145000, 0.8285622000, 
    0.8288097000, 0.8290571000, 0.8293043000, 0.8295512000, 0.8297980000, 
    0.8300446000, 0.8302910000, 0.8305373000, 0.8307833000, 0.8310292000, 
    0.8312748000, 0.8315203000, 0.8317656000, 0.8320107000, 0.8322557000, 
    0.8325004000, 0.8327450000, 0.8329894000, 0.8332336000, 0.8334776000, 
    0.8337214000, 0.8339651000, 0.8342085000, 0.8344518000, 0.8346949000 
};
 
double K1270_width[]={
    0.0000000014, 0.0000000052, 0.0000000127, 0.0000000252, 0.0000000442, 
    0.0000000709, 0.0000001067, 0.0000001528, 0.0000002108, 0.0000002819, 
    0.0000003675, 0.0000004689, 0.0000005876, 0.0000007249, 0.0000008822, 
    0.0000010609, 0.0000012624, 0.0000014882, 0.0000017396, 0.0000020181, 
    0.0000023251, 0.0000026620, 0.0000030303, 0.0000034314, 0.0000038668, 
    0.0000043379, 0.0000048463, 0.0000053933, 0.0000059805, 0.0000066094, 
    0.0000072815, 0.0000079983, 0.0000087612, 0.0000095719, 0.0000104318, 
    0.0000113424, 0.0000123054, 0.0000133223, 0.0000143946, 0.0000155239, 
    0.0000167117, 0.0000179597, 0.0000192695, 0.0000206426, 0.0000220807, 
    0.0000235853, 0.0000251581, 0.0000268008, 0.0000285149, 0.0000303022, 
    0.0000321643, 0.0000341028, 0.0000361195, 0.0000382160, 0.0000403941, 
    0.0000426554, 0.0000450017, 0.0000474347, 0.0000499561, 0.0000525678, 
    0.0000552715, 0.0000580689, 0.0000609619, 0.0000639522, 0.0000670417, 
    0.0000702323, 0.0000735256, 0.0000769237, 0.0000804283, 0.0000840414, 
    0.0000877648, 0.0000916004, 0.0000955503, 0.0000996163, 0.0001038003, 
    0.0001081045, 0.0001125306, 0.0001170808, 0.0001217571, 0.0001265615, 
    0.0001314961, 0.0001365630, 0.0001417641, 0.0001471018, 0.0001525781, 
    0.0001581951, 0.0001639551, 0.0001698602, 0.0001759127, 0.0001821148, 
    0.0001884687, 0.0001949769, 0.0002016415, 0.0002084650, 0.0002154497, 
    0.0002225981, 0.0002299124, 0.0002373952, 0.0002450491, 0.0002528764, 
    0.0002608798, 0.0002690618, 0.0002774250, 0.0002859721, 0.0002947058, 
    0.0003036288, 0.0003127438, 0.0003220537, 0.0003315613, 0.0003412695, 
    0.0003511812, 0.0003612993, 0.0003716269, 0.0003821671, 0.0003929230, 
    0.0004038976, 0.0004150944, 0.0004265164, 0.0004381671, 0.0004500499, 
    0.0004621681, 0.0004745254, 0.0004871253, 0.0004999714, 0.0005130674, 
    0.0005264172, 0.0005400246, 0.0005538936, 0.0005680282, 0.0005824325, 
    0.0005971108, 0.0006120672, 0.0006273063, 0.0006428325, 0.0006586504, 
    0.0006747647, 0.0006911803, 0.0007079021, 0.0007249351, 0.0007422845, 
    0.0007599558, 0.0007779542, 0.0007962855, 0.0008149555, 0.0008339699, 
    0.0008533351, 0.0008730571, 0.0008931425, 0.0009135979, 0.0009344302, 
    0.0009556463, 0.0009772537, 0.0009992598, 0.0010216720, 0.0010444990, 
    0.0010677490, 0.0010914300, 0.0011155510, 0.0011401210, 0.0011651510, 
    0.0011906480, 0.0012166250, 0.0012430910, 0.0012700570, 0.0012975350, 
    0.0013255370, 0.0013540740, 0.0013831600, 0.0014128090, 0.0014430330, 
    0.0014738470, 0.0015052670, 0.0015373090, 0.0015699870, 0.0016033210, 
    0.0016373270, 0.0016720240, 0.0017074330, 0.0017435730, 0.0017804660, 
    0.0018181340, 0.0018566010, 0.0018958920, 0.0019360320, 0.0019770490, 
    0.0020189710, 0.0020618280, 0.0021056510, 0.0021504740, 0.0021963300, 
    0.0022432570, 0.0022912920, 0.0023404760, 0.0023908500, 0.0024424590, 
    0.0024953480, 0.0025495680, 0.0026051670, 0.0026622000, 0.0027207220, 
    0.0027807900, 0.0028424670, 0.0029058130, 0.0029708960, 0.0030377810, 
    0.0031065390, 0.0031772420, 0.0032499620, 0.0033247750, 0.0034017560, 
    0.0034809810, 0.0035625250, 0.0036464650, 0.0037328720, 0.0038218190, 
    0.0039133710, 0.0040075910, 0.0041045360, 0.0042042540, 0.0043067850, 
    0.0044121600, 0.0045203990, 0.0046315090, 0.0047454840, 0.0048623040, 
    0.0049819360, 0.0051043310, 0.0052294280, 0.0053571500, 0.0054874080, 
    0.0056201030, 0.0057551230, 0.0058923500, 0.0060316590, 0.0061729200, 
    0.0063159990, 0.0064607640, 0.0066070800, 0.0067548180, 0.0069038480, 
    0.0070540500, 0.0072053060, 0.0073575050, 0.0075105450, 0.0076643270, 
    0.0078187630, 0.0079737720, 0.0081292780, 0.0082852130, 0.0084415170, 
    0.0085981350, 0.0087550200, 0.0089121270, 0.0090694200, 0.0092268660, 
    0.0093844380, 0.0095421120, 0.0096998670, 0.0098576870, 0.0100155600, 
    0.0101734700, 0.0103314200, 0.0104894000, 0.0106474000, 0.0108054200, 
    0.0109634700, 0.0111215500, 0.0112796500, 0.0114377900, 0.0115959700, 
    0.0117541900, 0.0119124600, 0.0120707900, 0.0122292000, 0.0123876800, 
    0.0125462500, 0.0127049200, 0.0128637100, 0.0130226100, 0.0131816400, 
    0.0133408200, 0.0135001600, 0.0136596700, 0.0138193500, 0.0139792400, 
    0.0141393300, 0.0142996400, 0.0144601800, 0.0146209700, 0.0147820200, 
    0.0149433400, 0.0151049500, 0.0152668600, 0.0154290800, 0.0155916300, 
    0.0157545200, 0.0159177700, 0.0160813900, 0.0162453900, 0.0164097800, 
    0.0165745900, 0.0167398200, 0.0169054800, 0.0170716000, 0.0172381900, 
    0.0174052600, 0.0175728200, 0.0177408900, 0.0179094800, 0.0180786100, 
    0.0182483000, 0.0184185500, 0.0185893800, 0.0187608000, 0.0189328400, 
    0.0191054900, 0.0192787900, 0.0194527400, 0.0196273600, 0.0198026700, 
    0.0199786700, 0.0201553900, 0.0203328300, 0.0205110200, 0.0206899700, 
    0.0208696900, 0.0210502100, 0.0212315300, 0.0214136700, 0.0215966500, 
    0.0217804800, 0.0219651900, 0.0221507800, 0.0223372700, 0.0225246900, 
    0.0227130400, 0.0229023400, 0.0230926100, 0.0232838700, 0.0234761300, 
    0.0236694200, 0.0238637400, 0.0240591200, 0.0242555800, 0.0244531300, 
    0.0246517900, 0.0248515800, 0.0250525200, 0.0252546300, 0.0254579200, 
    0.0256624200, 0.0258681400, 0.0260751100, 0.0262833500, 0.0264928700, 
    0.0267036900, 0.0269158400, 0.0271293400, 0.0273442100, 0.0275604600, 
    0.0277781300, 0.0279972300, 0.0282177800, 0.0284398100, 0.0286633400, 
    0.0288884000, 0.0291149900, 0.0293431600, 0.0295729200, 0.0298042900, 
    0.0300373100, 0.0302719900, 0.0305083600, 0.0307464500, 0.0309862700, 
    0.0312278600, 0.0314712500, 0.0317164500, 0.0319635000, 0.0322124200, 
    0.0324632400, 0.0327159900, 0.0329706900, 0.0332273800, 0.0334860800, 
    0.0337468200, 0.0340096300, 0.0342745400, 0.0345415800, 0.0348107800, 
    0.0350821700, 0.0353557800, 0.0356316400, 0.0359097900, 0.0361902600, 
    0.0364730800, 0.0367582700, 0.0370458900, 0.0373359500, 0.0376284900, 
    0.0379235500, 0.0382211500, 0.0385213500, 0.0388241600, 0.0391296300, 
    0.0394377800, 0.0397486700, 0.0400623100, 0.0403787600, 0.0406980400, 
    0.0410202000, 0.0413452600, 0.0416732700, 0.0420042700, 0.0423382900, 
    0.0426753700, 0.0430155500, 0.0433588700, 0.0437053700, 0.0440550800, 
    0.0444080500, 0.0447643100, 0.0451239000, 0.0454868700, 0.0458532400, 
    0.0462230700, 0.0465963900, 0.0469732300, 0.0473536500, 0.0477376700, 
    0.0481253400, 0.0485166900, 0.0489117600, 0.0493106000, 0.0497132400, 
    0.0501197100, 0.0505300600, 0.0509443200, 0.0513625300, 0.0517847200, 
    0.0522109400, 0.0526412100, 0.0530755700, 0.0535140600, 0.0539567000, 
    0.0544035400, 0.0548545900, 0.0553099000, 0.0557694900, 0.0562334000, 
    0.0567016400, 0.0571742500, 0.0576512600, 0.0581326800, 0.0586185500, 
    0.0591088800, 0.0596036900, 0.0601030200, 0.0606068600, 0.0611152500, 
    0.0616281900, 0.0621457100, 0.0626678100, 0.0631945000, 0.0637258000, 
    0.0642617000, 0.0648022200, 0.0653473700, 0.0658971300, 0.0664515200, 
    0.0670105300, 0.0675741500, 0.0681423900, 0.0687152300, 0.0692926600, 
    0.0698746700, 0.0704612600, 0.0710524000, 0.0716480700, 0.0722482600, 
    0.0728529400, 0.0734620900, 0.0740756900, 0.0746937100, 0.0753161100, 
    0.0759428600, 0.0765739400, 0.0772093000, 0.0778489100, 0.0784927300, 
    0.0791407200, 0.0797928300, 0.0804490200, 0.0811092400, 0.0817734600, 
    0.0824416000, 0.0831136300, 0.0837895000, 0.0844691400, 0.0851525100, 
    0.0858395400, 0.0865301900, 0.0872243800, 0.0879220700, 0.0886231900, 
    0.0893276700, 0.0900354700, 0.0907465100, 0.0914607200, 0.0921780600, 
    0.0928984400, 0.0936218200, 0.0943481100, 0.0950772500, 0.0958091900, 
    0.1023867000, 0.1055505000, 0.1081475000, 0.1104523000, 0.1125711000, 
    0.1145581000, 0.1164458000, 0.1182550000, 0.1200004000, 0.1216926000, 
    0.1233395000, 0.1249474000, 0.1265212000, 0.1280649000, 0.1295817000, 
    0.1310743000, 0.1325452000, 0.1339961000, 0.1354288000, 0.1368448000, 
    0.1382454000, 0.1396316000, 0.1410044000, 0.1423647000, 0.1437133000, 
    0.1450510000, 0.1463782000, 0.1476956000, 0.1490038000, 0.1503031000, 
    0.1515941000, 0.1528770000, 0.1541523000, 0.1554203000, 0.1566812000, 
    0.1579354000, 0.1591832000, 0.1604247000, 0.1616602000, 0.1628899000, 
    0.1641140000, 0.1653326000, 0.1665460000, 0.1677543000, 0.1689576000, 
    0.1701561000, 0.1713499000, 0.1725391000, 0.1737239000, 0.1749043000, 
    0.1760805000, 0.1772526000, 0.1784206000, 0.1795847000, 0.1807449000, 
    0.1819014000, 0.1830541000, 0.1842031000, 0.1853487000, 0.1864907000, 
    0.1876292000, 0.1887644000, 0.1898963000, 0.1910250000, 0.1921504000, 
    0.1932727000, 0.1943919000, 0.1955080000, 0.1966211000, 0.1977313000, 
    0.1988386000, 0.1999430000, 0.2010446000, 0.2021434000, 0.2032395000, 
    0.2043329000, 0.2054236000, 0.2065117000, 0.2075972000, 0.2086801000, 
    0.2097605000, 0.2108384000, 0.2119139000, 0.2129870000, 0.2140576000, 
    0.2151259000, 0.2161919000, 0.2172555000, 0.2183169000, 0.2193761000, 
    0.2204330000, 0.2214878000, 0.2225404000, 0.2235909000, 0.2246392000, 
    0.2256855000, 0.2267297000, 0.2277719000, 0.2288121000, 0.2298503000, 
    0.2308866000, 0.2319210000, 0.2329534000, 0.2339840000, 0.2350126000, 
    0.2360395000, 0.2370646000, 0.2380878000, 0.2391093000, 0.2401291000, 
    0.2411471000, 0.2421634000, 0.2431781000, 0.2441911000, 0.2452024000, 
    0.2462122000, 0.2472203000, 0.2482269000, 0.2492319000, 0.2502354000, 
    0.2512374000, 0.2522379000, 0.2532369000, 0.2542345000, 0.2552306000, 
    0.2562253000, 0.2572187000, 0.2582106000, 0.2592012000, 0.2601905000, 
    0.2611784000, 0.2621651000, 0.2631504000, 0.2641345000, 0.2651174000, 
    0.2660991000, 0.2670795000, 0.2680587000, 0.2690368000, 0.2700137000, 
    0.2709895000, 0.2719642000, 0.2729378000, 0.2739103000, 0.2748817000, 
    0.2758521000, 0.2768215000, 0.2777898000, 0.2787572000, 0.2797235000, 
    0.2806890000, 0.2816534000, 0.2826170000, 0.2835796000, 0.2845414000, 
    0.2855022000, 0.2864623000, 0.2874214000, 0.2883798000, 0.2893373000, 
    0.2902940000, 0.2912500000, 0.2922052000, 0.2931597000, 0.2941134000, 
    0.2950664000, 0.2960187000, 0.2969704000, 0.2979213000, 0.2988716000, 
    0.2998213000, 0.3007704000, 0.3017189000, 0.3026667000, 0.3036140000, 
    0.3045608000, 0.3055070000, 0.3064527000, 0.3073978000, 0.3083425000, 
    0.3092867000, 0.3102304000, 0.3111737000, 0.3121166000, 0.3130590000, 
    0.3140010000, 0.3149426000, 0.3158839000, 0.3168248000, 0.3177653000, 
    0.3187056000, 0.3196455000, 0.3205851000, 0.3215244000, 0.3224634000, 
    0.3234022000, 0.3243408000, 0.3252791000, 0.3262172000, 0.3271551000, 
    0.3280928000, 0.3290303000, 0.3299677000, 0.3309050000, 0.3318421000, 
    0.3327791000, 0.3337160000, 0.3346528000, 0.3355896000, 0.3365263000, 
    0.3374629000, 0.3383995000, 0.3393362000, 0.3402728000, 0.3412094000, 
    0.3421460000, 0.3430827000, 0.3440194000, 0.3449562000, 0.3458931000, 
    0.3468301000, 0.3477672000, 0.3487044000, 0.3496418000, 0.3505793000, 
    0.3515170000, 0.3524548000, 0.3533929000, 0.3543311000, 0.3552696000, 
    0.3562083000, 0.3571472000, 0.3580864000, 0.3590259000, 0.3599657000, 
    0.3609057000, 0.3618461000, 0.3627868000, 0.3637279000, 0.3646693000, 
    0.3656111000, 0.3665532000, 0.3674958000, 0.3684387000, 0.3693821000, 
    0.3703259000, 0.3712701000, 0.3722149000, 0.3731600000, 0.3741057000, 
    0.3750519000, 0.3759986000, 0.3769458000, 0.3778936000, 0.3788419000, 
    0.3797908000, 0.3807402000, 0.3816903000, 0.3826409000, 0.3835922000, 
    0.3845441000, 0.3854966000, 0.3864498000, 0.3874037000, 0.3883583000, 
    0.3893135000, 0.3902695000, 0.3912262000, 0.3921836000, 0.3931418000, 
    0.3941007000, 0.3950604000, 0.3960209000, 0.3969822000, 0.3979443000, 
    0.3989072000, 0.3998710000, 0.4008356000, 0.4018011000, 0.4027675000, 
    0.4037347000, 0.4047029000, 0.4056720000, 0.4066420000, 0.4076129000, 
    0.4085848000, 0.4095576000, 0.4105315000, 0.4115063000, 0.4124821000, 
    0.4134590000, 0.4144368000, 0.4154157000, 0.4163957000, 0.4173767000, 
    0.4183588000, 0.4193420000, 0.4203264000, 0.4213118000, 0.4222983000, 
    0.4232860000, 0.4242749000, 0.4252649000, 0.4262561000, 0.4272484000, 
    0.4282420000, 0.4292368000, 0.4302328000, 0.4312301000, 0.4322286000, 
    0.4332284000, 0.4342295000, 0.4352318000, 0.4362354000, 0.4372404000, 
    0.4382467000, 0.4392543000, 0.4402633000, 0.4412736000, 0.4422853000, 
    0.4432984000, 0.4443129000, 0.4453288000, 0.4463461000, 0.4473648000, 
    0.4483850000, 0.4494067000, 0.4504298000, 0.4514544000, 0.4524804000, 
    0.4535080000, 0.4545371000, 0.4555678000, 0.4565999000, 0.4576337000, 
    0.4586689000, 0.4597058000, 0.4607442000, 0.4617843000, 0.4628259000, 
    0.4638692000, 0.4649141000, 0.4659606000, 0.4670088000, 0.4680586000, 
    0.4691101000, 0.4701633000, 0.4712183000, 0.4722749000, 0.4733332000, 
    0.4743933000, 0.4754551000, 0.4765186000, 0.4775840000, 0.4786511000, 
    0.4797199000, 0.4807906000, 0.4818631000, 0.4829374000, 0.4840135000, 
    0.4850915000, 0.4861713000, 0.4872530000, 0.4883366000, 0.4894220000, 
    0.4905093000, 0.4915985000, 0.4926897000, 0.4937827000, 0.4948777000, 
    0.4959747000, 0.4970736000, 0.4981744000, 0.4992772000, 0.5003820000, 
    0.5014888000, 0.5025976000, 0.5037085000, 0.5048213000, 0.5059362000, 
    0.5070531000, 0.5081720000, 0.5092931000, 0.5104162000, 0.5115413000, 
    0.5126686000, 0.5137980000, 0.5149294000, 0.5160630000, 0.5171987000, 
    0.5183366000, 0.5194766000, 0.5206187000, 0.5217630000, 0.5229095000, 
    0.5240581000, 0.5252090000, 0.5263620000, 0.5275173000, 0.5286748000, 
    0.5298344000, 0.5309964000, 0.5321605000, 0.5333269000, 0.5344956000, 
    0.5356665000, 0.5368397000, 0.5380152000, 0.5391930000, 0.5403731000, 
    0.5415554000, 0.5427401000, 0.5439272000, 0.5451165000, 0.5463082000, 
    0.5475022000, 0.5486986000, 0.5498973000, 0.5510984000, 0.5523019000, 
    0.5535077000, 0.5547159000, 0.5559266000, 0.5571396000, 0.5583551000, 
    0.5595729000, 0.5607932000, 0.5620159000, 0.5632411000, 0.5644686000, 
    0.5656987000, 0.5669312000, 0.5681661000, 0.5694036000, 0.5706435000, 
    0.5718859000, 0.5731307000, 0.5743781000, 0.5756280000, 0.5768804000, 
    0.5781352000, 0.5793926000, 0.5806526000, 0.5819150000, 0.5831800000, 
    0.5844476000, 0.5857176000, 0.5869903000, 0.5882655000, 0.5895432000, 
    0.5908235000, 0.5921064000, 0.5933919000, 0.5946799000, 0.5959706000, 
    0.5972638000, 0.5985596000, 0.5998581000, 0.6011591000, 0.6024627000, 
    0.6037690000, 0.6050779000, 0.6063894000, 0.6077036000, 0.6090203000, 
    0.6103397000, 0.6116618000, 0.6129865000, 0.6143139000, 0.6156439000, 
    0.6169765000, 0.6183119000, 0.6196499000, 0.6209905000, 0.6223339000, 
    0.6236799000, 0.6250286000, 0.6263800000, 0.6277340000, 0.6290908000, 
    0.6304503000, 0.6318124000, 0.6331773000, 0.6345448000, 0.6359151000, 
    0.6372881000, 0.6386638000, 0.6400422000, 0.6414233000, 0.6428071000, 
    0.6441937000, 0.6455830000, 0.6469750000, 0.6483698000, 0.6497673000, 
    0.6511675000, 0.6525705000, 0.6539762000, 0.6553846000, 0.6567958000, 
    0.6582098000, 0.6596264000, 0.6610459000, 0.6624681000, 0.6638930000, 
    0.6653208000, 0.6667512000, 0.6681845000, 0.6696205000, 0.6710592000, 
    0.6725007000, 0.6739450000, 0.6753921000, 0.6768419000, 0.6782945000, 
    0.6797499000, 0.6812080000, 0.6826690000, 0.6841327000, 0.6855992000, 
    0.6870684000, 0.6885404000, 0.6900153000, 0.6914929000, 0.6929732000, 
    0.6944564000, 0.6959423000, 0.6974311000, 0.6989226000, 0.7004169000, 
    0.7019140000, 0.7034139000, 0.7049165000, 0.7064220000, 0.7079302000, 
    0.7094412000, 0.7109550000, 0.7124716000, 0.7139910000, 0.7155132000, 
    0.7170382000, 0.7185659000, 0.7200965000, 0.7216298000, 0.7231660000, 
    0.7247049000, 0.7262466000, 0.7277911000, 0.7293383000, 0.7308884000, 
    0.7324413000, 0.7339969000, 0.7355553000, 0.7371166000, 0.7386806000, 
    0.7402474000, 0.7418169000, 0.7433893000, 0.7449644000, 0.7465424000, 
    0.7481231000, 0.7497066000, 0.7512929000, 0.7528819000, 0.7544737000, 
    0.7560684000, 0.7576657000, 0.7592659000, 0.7608689000, 0.7624746000, 
    0.7640831000, 0.7656943000, 0.7673084000, 0.7689252000, 0.7705448000, 
    0.7721671000, 0.7737922000, 0.7754201000, 0.7770507000, 0.7786841000, 
    0.7803203000, 0.7819592000, 0.7836009000, 0.7852454000, 0.7868926000, 
    0.7885425000, 0.7901952000, 0.7918507000, 0.7935089000, 0.7951698000, 
    0.7968335000, 0.7985000000, 0.8001692000, 0.8018411000, 0.8035158000, 
    0.8051932000, 0.8068733000, 0.8085562000, 0.8102418000, 0.8119301000, 
    0.8136212000, 0.8153150000, 0.8170115000, 0.8187108000, 0.8204127000, 
    0.8221174000, 0.8238248000, 0.8255349000, 0.8272477000, 0.8289633000, 
    0.8306815000, 0.8324025000, 0.8341261000, 0.8358525000, 0.8375815000, 
    0.8393133000, 0.8410477000, 0.8427849000, 0.8445247000, 0.8462672000, 
    0.8480124000, 0.8497603000, 0.8515108000, 0.8532641000, 0.8550200000, 
    0.8567786000, 0.8585399000, 0.8603038000, 0.8620704000, 0.8638396000, 
    0.8656116000, 0.8673861000, 0.8691634000, 0.8709433000, 0.8727258000, 
    0.8745110000, 0.8762988000, 0.8780893000, 0.8798824000, 0.8816782000, 
    0.8834766000, 0.8852776000, 0.8870812000, 0.8888875000, 0.8906964000, 
    0.8925079000, 0.8943221000, 0.8961388000, 0.8979582000, 0.8997801000, 
    0.9016047000, 0.9034319000, 0.9052617000, 0.9070941000, 0.9089291000, 
    0.9107666000, 0.9126068000, 0.9144495000, 0.9162949000, 0.9181428000, 
    0.9199933000, 0.9218463000, 0.9237020000, 0.9255602000, 0.9274210000, 
    0.9292843000, 0.9311502000, 0.9330186000, 0.9348896000, 0.9367632000, 
    0.9386393000, 0.9405179000, 0.9423991000, 0.9442828000, 0.9461691000, 
    0.9480579000, 0.9499492000, 0.9518431000, 0.9537394000, 0.9556383000, 
    0.9575397000, 0.9594436000, 0.9613501000, 0.9632590000, 0.9651704000, 
    0.9670844000, 0.9690008000, 0.9709197000, 0.9728411000, 0.9747650000, 
    0.9766914000, 0.9786203000, 0.9805516000, 0.9824854000, 0.9844217000, 
    0.9863605000, 0.9883017000, 0.9902454000, 0.9921915000, 0.9941401000, 
    0.9960911000, 0.9980446000, 1.0000010000, 1.0019590000, 1.0039200000, 
    1.0058830000, 1.0078490000, 1.0098170000, 1.0117870000, 1.0137600000, 
    1.0157360000, 1.0177130000, 1.0196940000, 1.0216760000, 1.0236610000, 
    1.0256490000, 1.0276380000, 1.0296310000, 1.0316250000, 1.0336220000, 
    1.0356210000, 1.0376230000, 1.0396270000, 1.0416340000, 1.0436430000, 
    1.0456540000, 1.0476670000, 1.0496830000, 1.0517020000, 1.0537220000, 
    1.0557450000, 1.0577710000, 1.0597980000, 1.0618280000, 1.0638610000, 
    1.0658950000, 1.0679320000, 1.0699720000, 1.0720130000, 1.0740570000, 
    1.0761040000, 1.0781520000, 1.0802030000, 1.0822560000, 1.0843120000, 
    1.0863690000, 1.0884300000, 1.0904920000, 1.0925570000, 1.0946240000, 
    1.0966930000, 1.0987640000, 1.1008380000, 1.1029140000, 1.1049920000, 
    1.1070730000, 1.1091560000, 1.1112410000, 1.1133280000, 1.1154180000, 
    1.1175090000, 1.1196030000, 1.1217000000, 1.1237980000, 1.1258990000, 
    1.1280020000, 1.1301070000, 1.1322140000, 1.1343240000, 1.1364360000, 
    1.1385500000, 1.1406660000, 1.1427850000, 1.1449050000, 1.1470280000, 
    1.1491530000, 1.1512800000, 1.1534100000, 1.1555410000, 1.1576750000, 
    1.1598110000, 1.1619490000, 1.1640890000, 1.1662310000, 1.1683760000, 
    1.1705230000, 1.1726710000, 1.1748220000, 1.1769750000, 1.1791310000, 
    1.1812880000, 1.1834480000, 1.1856090000, 1.1877730000, 1.1899390000, 
    1.1921070000, 1.1942770000, 1.1964490000, 1.1986240000, 1.2008000000, 
    1.2029780000, 1.2051590000, 1.2073420000, 1.2095270000, 1.2117130000, 
    1.2139020000, 1.2160930000, 1.2182860000, 1.2204820000, 1.2226790000, 
    1.2248780000, 1.2270790000, 1.2292830000, 1.2314880000, 1.2336960000, 
    1.2359050000, 1.2381170000, 1.2403300000, 1.2425460000, 1.2447640000, 
    1.2469830000, 1.2492050000, 1.2514290000, 1.2536540000, 1.2558820000, 
    1.2581120000, 1.2603430000, 1.2625770000, 1.2648130000, 1.2670510000, 
    1.2692900000, 1.2715320000, 1.2737760000, 1.2760210000, 1.2782690000, 
    1.2805180000, 1.2827700000, 1.2850230000, 1.2872790000, 1.2895360000, 
    1.2917950000, 1.2940560000, 1.2963200000, 1.2985850000, 1.3008520000, 
    1.3031210000, 1.3053920000, 1.3076650000, 1.3099390000, 1.3122160000, 
    1.3144940000, 1.3167750000, 1.3190570000, 1.3213420000, 1.3236280000, 
    1.3259160000, 1.3282060000, 1.3304970000, 1.3327910000, 1.3350870000, 
    1.3373840000, 1.3396830000, 1.3419850000, 1.3442880000, 1.3465920000, 
    1.3488990000, 1.3512080000, 1.3535180000, 1.3558300000, 1.3581440000, 
    1.3604600000, 1.3627780000, 1.3650970000, 1.3674190000, 1.3697420000, 
    1.3720670000, 1.3743940000, 1.3767220000, 1.3790530000, 1.3813850000, 
    1.3837190000, 1.3860540000, 1.3883920000, 1.3907310000, 1.3930720000, 
    1.3954150000, 1.3977600000, 1.4001060000, 1.4024540000, 1.4048040000, 
    1.4071560000, 1.4095090000, 1.4118640000, 1.4142210000, 1.4165800000, 
    1.4189400000, 1.4213020000, 1.4236660000, 1.4260310000, 1.4283990000, 
    1.4307670000, 1.4331380000, 1.4355100000, 1.4378840000, 1.4402600000, 
    1.4426370000, 1.4450160000, 1.4473970000, 1.4497800000, 1.4521640000, 
    1.4545490000, 1.4569370000, 1.4593260000, 1.4617170000, 1.4641090000, 
    1.4665030000, 1.4688990000, 1.4712960000, 1.4736950000, 1.4760960000, 
    1.4784980000, 1.4809020000, 1.4833070000, 1.4857140000, 1.4881230000, 
    1.4905330000, 1.4929450000, 1.4953580000, 1.4977730000, 1.5001900000, 
    1.5026080000, 1.5050280000, 1.5074490000, 1.5098720000, 1.5122970000, 
    1.5147230000, 1.5171500000, 1.5195790000, 1.5220100000, 1.5244420000, 
    1.5268760000, 1.5293110000, 1.5317480000, 1.5341870000, 1.5366270000, 
    1.5390680000, 1.5415110000, 1.5439550000, 1.5464010000, 1.5488490000, 
    1.5512980000, 1.5537480000, 1.5562000000, 1.5586530000, 1.5611080000, 
    1.5635650000, 1.5660230000, 1.5684820000, 1.5709430000, 1.5734050000, 
    1.5758680000, 1.5783340000, 1.5808000000, 1.5832680000, 1.5857380000, 
    1.5882080000, 1.5906810000, 1.5931540000, 1.5956290000, 1.5981060000, 
    1.6005840000, 1.6030630000, 1.6055440000, 1.6080260000, 1.6105100000, 
    1.6129950000, 1.6154810000, 1.6179680000, 1.6204570000, 1.6229480000, 
    1.6254400000, 1.6279330000, 1.6304270000, 1.6329230000, 1.6354200000, 
    1.6379190000, 1.6404190000, 1.6429200000, 1.6454220000, 1.6479260000, 
    1.6504310000, 1.6529380000, 1.6554450000, 1.6579540000, 1.6604650000, 
    1.6629760000, 1.6654890000, 1.6680040000, 1.6705190000, 1.6730360000, 
    1.6755540000, 1.6780730000, 1.6805940000, 1.6831160000, 1.6856390000, 
    1.6881630000, 1.6906890000, 1.6932160000, 1.6957440000, 1.6982730000, 
    1.7008030000, 1.7033350000, 1.7058680000, 1.7084020000, 1.7109380000, 
    1.7134740000, 1.7160120000, 1.7185510000, 1.7210910000, 1.7236330000, 
    1.7261750000, 1.7287190000, 1.7312640000, 1.7338100000, 1.7363570000, 
    1.7389050000, 1.7414550000, 1.7440050000, 1.7465570000, 1.7491100000, 
    1.7516640000, 1.7542190000, 1.7567760000, 1.7593330000, 1.7618920000, 
    1.7644520000, 1.7670120000, 1.7695740000, 1.7721370000, 1.7747010000, 
    1.7772660000, 1.7798330000, 1.7824000000, 1.7849680000, 1.7875380000, 
    1.7901080000, 1.7926800000, 1.7952530000, 1.7978260000, 1.8004010000, 
    1.8029770000, 1.8055540000, 1.8081320000, 1.8107100000, 1.8132900000, 
    1.8158710000, 1.8184530000, 1.8210360000, 1.8236200000, 1.8262050000, 
    1.8287910000, 1.8313780000, 1.8339660000, 1.8365550000, 1.8391450000, 
    1.8417360000, 1.8443270000, 1.8469200000, 1.8495140000, 1.8521090000, 
    1.8547040000, 1.8573010000, 1.8598980000, 1.8624970000, 1.8650960000, 
    1.8676970000, 1.8702980000, 1.8729000000, 1.8755030000, 1.8781070000, 
    1.8807120000, 1.8833180000, 1.8859250000, 1.8885320000, 1.8911410000, 
    1.8937500000, 1.8963600000, 1.8989710000, 1.9015830000, 1.9041960000, 
    1.9068100000, 1.9094240000, 1.9120400000, 1.9146560000, 1.9172730000, 
    1.9198910000, 1.9225100000, 1.9251290000, 1.9277500000, 1.9303710000, 
    1.9329930000, 1.9356160000, 1.9382400000, 1.9408640000, 1.9434890000, 
    1.9461150000, 1.9487420000, 1.9513700000, 1.9539980000, 1.9566270000, 
    1.9592570000, 1.9618880000, 1.9645190000, 1.9671520000, 1.9697850000, 
    1.9724180000, 1.9750530000, 1.9776880000, 1.9803240000, 1.9829610000, 
    1.9855980000, 1.9882360000, 1.9908750000, 1.9935140000, 1.9961550000, 
    1.9987960000, 2.0014370000, 2.0040800000, 2.0067230000, 2.0093660000, 
    2.0120110000, 2.0146560000, 2.0173010000, 2.0199480000, 2.0225950000, 
    2.0252430000, 2.0278910000, 2.0305400000, 2.0331900000, 2.0358400000, 
    2.0384910000, 2.0411430000, 2.0437950000, 2.0464480000, 2.0491010000, 
    2.0517550000, 2.0544100000, 2.0570650000, 2.0597210000, 2.0623780000, 
    2.0650350000, 2.0676930000, 2.0703510000, 2.0730100000, 2.0756690000, 
    2.0783290000, 2.0809900000, 2.0836510000, 2.0863130000, 2.0889750000, 
    2.0916380000, 2.0943010000, 2.0969650000, 2.0996290000, 2.1022940000, 
    2.1049600000, 2.1076260000, 2.1102930000, 2.1129600000, 2.1156270000, 
    2.1182950000, 2.1209640000, 2.1236330000, 2.1263030000, 2.1289730000, 
    2.1316440000, 2.1343150000, 2.1369860000, 2.1396580000, 2.1423310000, 
    2.1450040000, 2.1476780000, 2.1503510000, 2.1530260000, 2.1557010000
};
 
double K1400_width[]={
    0.0000000007, 0.0000000025, 0.0000000061, 0.0000000121, 0.0000000213, 
    0.0000000344, 0.0000000519, 0.0000000747, 0.0000001034, 0.0000001389, 
    0.0000001818, 0.0000002329, 0.0000002930, 0.0000003630, 0.0000004435, 
    0.0000005356, 0.0000006400, 0.0000007576, 0.0000008892, 0.0000010359, 
    0.0000011985, 0.0000013780, 0.0000015753, 0.0000017914, 0.0000020274, 
    0.0000022842, 0.0000025629, 0.0000028646, 0.0000031903, 0.0000035413, 
    0.0000039185, 0.0000043232, 0.0000047566, 0.0000052199, 0.0000057143, 
    0.0000062411, 0.0000068015, 0.0000073970, 0.0000080288, 0.0000086983, 
    0.0000094069, 0.0000101561, 0.0000109474, 0.0000117821, 0.0000126619, 
    0.0000135884, 0.0000145631, 0.0000155876, 0.0000166637, 0.0000177930, 
    0.0000189773, 0.0000202184, 0.0000215182, 0.0000228784, 0.0000243010, 
    0.0000257879, 0.0000273412, 0.0000289629, 0.0000306551, 0.0000324199, 
    0.0000342596, 0.0000361763, 0.0000381724, 0.0000402502, 0.0000424121, 
    0.0000446606, 0.0000469982, 0.0000494275, 0.0000519511, 0.0000545718, 
    0.0000572923, 0.0000601155, 0.0000630443, 0.0000660816, 0.0000692306, 
    0.0000724943, 0.0000758761, 0.0000793791, 0.0000830068, 0.0000867626, 
    0.0000906501, 0.0000946730, 0.0000988349, 0.0001031398, 0.0001075916, 
    0.0001121943, 0.0001169520, 0.0001218691, 0.0001269499, 0.0001321989, 
    0.0001376208, 0.0001432202, 0.0001490021, 0.0001549715, 0.0001611335, 
    0.0001674933, 0.0001740566, 0.0001808287, 0.0001878156, 0.0001950231, 
    0.0002024573, 0.0002101244, 0.0002180310, 0.0002261837, 0.0002345892, 
    0.0002432546, 0.0002521872, 0.0002613945, 0.0002708840, 0.0002806638, 
    0.0002907420, 0.0003011271, 0.0003118276, 0.0003228526, 0.0003342113, 
    0.0003459132, 0.0003579682, 0.0003703863, 0.0003831781, 0.0003963544, 
    0.0004099263, 0.0004239055, 0.0004383037, 0.0004531333, 0.0004684072, 
    0.0004841383, 0.0005003405, 0.0005170277, 0.0005342145, 0.0005519161, 
    0.0005701480, 0.0005889265, 0.0006082683, 0.0006281907, 0.0006487119, 
    0.0006698504, 0.0006916256, 0.0007140575, 0.0007371670, 0.0007609758, 
    0.0007855062, 0.0008107815, 0.0008368260, 0.0008636648, 0.0008913240, 
    0.0009198309, 0.0009492137, 0.0009795018, 0.0010107260, 0.0010429180, 
    0.0010761110, 0.0011103390, 0.0011456390, 0.0011820480, 0.0012196060, 
    0.0012583530, 0.0012983320, 0.0013395870, 0.0013821650, 0.0014261150, 
    0.0014714870, 0.0015183340, 0.0015667120, 0.0016166800, 0.0016682970, 
    0.0017216280, 0.0017767390, 0.0018337000, 0.0018925850, 0.0019534700, 
    0.0020164360, 0.0020815670, 0.0021489530, 0.0022186860, 0.0022908640, 
    0.0023655910, 0.0024429750, 0.0025231290, 0.0026061740, 0.0026922350, 
    0.0027814460, 0.0028739440, 0.0029698790, 0.0030694050, 0.0031726840, 
    0.0032798900, 0.0033912030, 0.0035068140, 0.0036269240, 0.0037517470, 
    0.0038815040, 0.0040164320, 0.0041567790, 0.0043028050, 0.0044547850, 
    0.0046130090, 0.0047777790, 0.0049494150, 0.0051282500, 0.0053146360, 
    0.0055089370, 0.0057115350, 0.0059228300, 0.0061432340, 0.0063731760, 
    0.0066130990, 0.0068634590, 0.0071247200, 0.0073973600, 0.0076818580, 
    0.0079786990, 0.0082883670, 0.0086113370, 0.0089480750, 0.0092990290, 
    0.0096646210, 0.0100452400, 0.0104412500, 0.0108529300, 0.0112805500, 
    0.0117242700, 0.0121842100, 0.0126603900, 0.0131527500, 0.0136611300, 
    0.0141852900, 0.0147248800, 0.0152794700, 0.0158485100, 0.0164314100, 
    0.0170274700, 0.0176359300, 0.0182559700, 0.0188867400, 0.0195273500, 
    0.0201768900, 0.0208344400, 0.0214991000, 0.0221699600, 0.0228461600, 
    0.0235268400, 0.0242112100, 0.0248985100, 0.0255880100, 0.0262790500, 
    0.0269710000, 0.0276632900, 0.0283554000, 0.0290468400, 0.0297371800, 
    0.0304260300, 0.0311130400, 0.0317978800, 0.0324802900, 0.0331600000, 
    0.0338368100, 0.0345105300, 0.0351809800, 0.0358480300, 0.0365115500, 
    0.0371714500, 0.0378276300, 0.0384800200, 0.0391285600, 0.0397732200, 
    0.0404139500, 0.0410507400, 0.0416835500, 0.0423124000, 0.0429372700, 
    0.0435581700, 0.0441751100, 0.0447881100, 0.0453971800, 0.0460023400, 
    0.0466036300, 0.0472010600, 0.0477946700, 0.0483845000, 0.0489705600, 
    0.0495529100, 0.0501315800, 0.0507066000, 0.0512780200, 0.0518458700, 
    0.0524101800, 0.0529710200, 0.0535284000, 0.0540823700, 0.0546329800, 
    0.0551802500, 0.0557242400, 0.0562649800, 0.0568025200, 0.0573368800, 
    0.0578681100, 0.0583962500, 0.0589213300, 0.0594434000, 0.0599624900, 
    0.0604786400, 0.0609918800, 0.0615022500, 0.0620097800, 0.0625145200, 
    0.0630164900, 0.0635157300, 0.0640122700, 0.0645061500, 0.0649973900, 
    0.0654860300, 0.0659721100, 0.0664556400, 0.0669366600, 0.0674152100, 
    0.0678913000, 0.0683649800, 0.0688362600, 0.0693051800, 0.0697717600, 
    0.0702360300, 0.0706980200, 0.0711577500, 0.0716152500, 0.0720705400, 
    0.0725236500, 0.0729746100, 0.0734234300, 0.0738701400, 0.0743147600, 
    0.0747573300, 0.0751978500, 0.0756363500, 0.0760728500, 0.0765073800, 
    0.0769399600, 0.0773706000, 0.0777993300, 0.0782261600, 0.0786511200, 
    0.0790742300, 0.0794955000, 0.0799149600, 0.0803326100, 0.0807484900, 
    0.0811626100, 0.0815749800, 0.0819856300, 0.0823945700, 0.0828018200, 
    0.0832073900, 0.0836113000, 0.0840135700, 0.0844142200, 0.0848132500, 
    0.0852106900, 0.0856065500, 0.0860008400, 0.0863935900, 0.0867848000, 
    0.0871744900, 0.0875626700, 0.0879493600, 0.0883345700, 0.0887183200, 
    0.0891006200, 0.0894814800, 0.0898609100, 0.0902389300, 0.0906155600, 
    0.0909908000, 0.0913646700, 0.0917371700, 0.0921083300, 0.0924781500, 
    0.0928466500, 0.0932138400, 0.0935797200, 0.0939443200, 0.0943076400, 
    0.0946696900, 0.0950304800, 0.0953900400, 0.0957483600, 0.0961054500, 
    0.0964613400, 0.0968160300, 0.0971695200, 0.0975218300, 0.0978729800, 
    0.0982229600, 0.0985718000, 0.0989194900, 0.0992660600, 0.0996115100, 
    0.0999558400, 0.1002991000, 0.1006412000, 0.1009823000, 0.1013223000, 
    0.1016612000, 0.1019991000, 0.1023359000, 0.1026717000, 0.1030064000, 
    0.1033402000, 0.1036729000, 0.1040047000, 0.1043354000, 0.1046652000, 
    0.1049939000, 0.1053218000, 0.1056486000, 0.1059745000, 0.1062995000, 
    0.1066235000, 0.1069466000, 0.1072688000, 0.1075901000, 0.1079104000, 
    0.1082299000, 0.1085485000, 0.1088662000, 0.1091831000, 0.1094990000, 
    0.1098142000, 0.1101284000, 0.1104418000, 0.1107544000, 0.1110662000, 
    0.1113771000, 0.1116873000, 0.1119966000, 0.1123051000, 0.1126128000, 
    0.1129198000, 0.1132259000, 0.1135313000, 0.1138360000, 0.1141398000, 
    0.1144429000, 0.1147453000, 0.1150469000, 0.1153478000, 0.1156480000, 
    0.1159474000, 0.1162461000, 0.1165442000, 0.1168415000, 0.1171381000, 
    0.1174340000, 0.1177292000, 0.1180238000, 0.1183177000, 0.1186109000, 
    0.1189034000, 0.1191953000, 0.1194865000, 0.1197771000, 0.1200670000, 
    0.1203563000, 0.1206449000, 0.1209329000, 0.1212203000, 0.1215071000, 
    0.1217933000, 0.1220788000, 0.1223637000, 0.1226480000, 0.1229318000, 
    0.1232149000, 0.1234974000, 0.1237794000, 0.1240607000, 0.1243415000, 
    0.1246217000, 0.1249013000, 0.1251803000, 0.1254588000, 0.1257367000, 
    0.1260140000, 0.1262908000, 0.1265670000, 0.1268427000, 0.1271178000, 
    0.1273923000, 0.1276663000, 0.1279398000, 0.1282127000, 0.1284850000, 
    0.1287568000, 0.1290281000, 0.1292988000, 0.1295690000, 0.1298387000, 
    0.1301078000, 0.1303764000, 0.1306445000, 0.1309120000, 0.1311790000, 
    0.1314455000, 0.1317114000, 0.1319768000, 0.1322417000, 0.1325060000, 
    0.1327699000, 0.1330332000, 0.1332960000, 0.1335582000, 0.1338199000, 
    0.1340812000, 0.1343418000, 0.1346020000, 0.1348617000, 0.1351208000, 
    0.1355387000, 0.1358629000, 0.1361711000, 0.1364708000, 0.1367648000, 
    0.1370546000, 0.1373411000, 0.1376249000, 0.1379063000, 0.1381858000, 
    0.1384634000, 0.1387394000, 0.1390139000, 0.1392870000, 0.1395588000, 
    0.1398294000, 0.1400988000, 0.1403671000, 0.1406344000, 0.1409007000, 
    0.1411659000, 0.1414302000, 0.1416936000, 0.1419562000, 0.1422178000, 
    0.1424786000, 0.1427386000, 0.1429977000, 0.1432561000, 0.1435136000, 
    0.1437704000, 0.1440265000, 0.1442818000, 0.1445364000, 0.1447902000, 
    0.1450434000, 0.1452958000, 0.1455475000, 0.1457986000, 0.1460490000, 
    0.1462986000, 0.1465477000, 0.1467960000, 0.1470437000, 0.1472908000, 
    0.1475372000, 0.1477830000, 0.1480281000, 0.1482726000, 0.1485165000, 
    0.1487597000, 0.1490024000, 0.1492444000, 0.1494858000, 0.1497266000, 
    0.1499668000, 0.1502065000, 0.1504455000, 0.1506839000, 0.1509218000, 
    0.1511590000, 0.1513957000, 0.1516318000, 0.1518674000, 0.1521023000, 
    0.1523367000, 0.1525706000, 0.1528038000, 0.1530366000, 0.1532687000, 
    0.1535003000, 0.1537314000, 0.1539619000, 0.1541918000, 0.1544213000, 
    0.1546501000, 0.1548785000, 0.1551063000, 0.1553336000, 0.1555603000, 
    0.1557866000, 0.1560122000, 0.1562374000, 0.1564621000, 0.1566862000, 
    0.1569098000, 0.1571329000, 0.1573555000, 0.1575776000, 0.1577992000, 
    0.1580203000, 0.1582409000, 0.1584609000, 0.1586805000, 0.1588996000, 
    0.1591182000, 0.1593363000, 0.1595539000, 0.1597710000, 0.1599876000, 
    0.1602038000, 0.1604194000, 0.1606346000, 0.1608493000, 0.1610636000, 
    0.1612773000, 0.1614906000, 0.1617035000, 0.1619158000, 0.1621277000, 
    0.1623391000, 0.1625501000, 0.1627606000, 0.1629706000, 0.1631802000, 
    0.1633894000, 0.1635981000, 0.1638063000, 0.1640141000, 0.1642214000, 
    0.1644283000, 0.1646347000, 0.1648407000, 0.1650463000, 0.1652514000, 
    0.1654561000, 0.1656604000, 0.1658642000, 0.1660676000, 0.1662706000, 
    0.1664731000, 0.1666752000, 0.1668769000, 0.1670782000, 0.1672790000, 
    0.1674794000, 0.1676794000, 0.1678790000, 0.1680782000, 0.1682769000, 
    0.1684753000, 0.1686732000, 0.1688707000, 0.1690679000, 0.1692646000, 
    0.1694609000, 0.1696568000, 0.1698523000, 0.1700474000, 0.1702421000, 
    0.1704364000, 0.1706304000, 0.1708239000, 0.1710170000, 0.1712098000, 
    0.1714022000, 0.1715941000, 0.1717857000, 0.1719769000, 0.1721677000, 
    0.1723582000, 0.1725482000, 0.1727379000, 0.1729272000, 0.1731162000, 
    0.1733047000, 0.1734929000, 0.1736807000, 0.1738682000, 0.1740552000, 
    0.1742420000, 0.1744283000, 0.1746143000, 0.1747999000, 0.1749851000, 
    0.1751700000, 0.1753546000, 0.1755388000, 0.1757226000, 0.1759060000, 
    0.1760892000, 0.1762719000, 0.1764543000, 0.1766364000, 0.1768181000, 
    0.1769995000, 0.1771805000, 0.1773612000, 0.1775415000, 0.1777215000, 
    0.1779011000, 0.1780804000, 0.1782594000, 0.1784380000, 0.1786163000, 
    0.1787943000, 0.1789719000, 0.1791492000, 0.1793262000, 0.1795028000, 
    0.1796791000, 0.1798551000, 0.1800308000, 0.1802061000, 0.1803811000, 
    0.1805558000, 0.1807301000, 0.1809042000, 0.1810779000, 0.1812513000, 
    0.1814244000, 0.1815971000, 0.1817696000, 0.1819417000, 0.1821136000, 
    0.1822851000, 0.1824563000, 0.1826272000, 0.1827978000, 0.1829681000, 
    0.1831380000, 0.1833077000, 0.1834771000, 0.1836461000, 0.1838149000, 
    0.1839834000, 0.1841515000, 0.1843194000, 0.1844869000, 0.1846542000, 
    0.1848212000, 0.1849879000, 0.1851543000, 0.1853203000, 0.1854862000, 
    0.1856517000, 0.1858169000, 0.1859818000, 0.1861465000, 0.1863108000, 
    0.1864749000, 0.1866387000, 0.1868022000, 0.1869654000, 0.1871284000, 
    0.1872910000, 0.1874534000, 0.1876155000, 0.1877774000, 0.1879389000, 
    0.1881002000, 0.1882612000, 0.1884219000, 0.1885824000, 0.1887426000, 
    0.1889025000, 0.1890621000, 0.1892215000, 0.1893806000, 0.1895395000, 
    0.1896980000, 0.1898563000, 0.1900144000, 0.1901722000, 0.1903297000, 
    0.1904869000, 0.1906439000, 0.1908007000, 0.1909571000, 0.1911133000, 
    0.1912693000, 0.1914250000, 0.1915805000, 0.1917357000, 0.1918906000, 
    0.1920453000, 0.1921997000, 0.1923539000, 0.1925078000, 0.1926615000, 
    0.1928149000, 0.1929681000, 0.1931211000, 0.1932738000, 0.1934262000, 
    0.1935784000, 0.1937304000, 0.1938821000, 0.1940336000, 0.1941848000, 
    0.1943358000, 0.1944865000, 0.1946370000, 0.1947873000, 0.1949374000, 
    0.1950872000, 0.1952367000, 0.1953860000, 0.1955351000, 0.1956840000, 
    0.1958326000, 0.1959810000, 0.1961292000, 0.1962771000, 0.1964248000, 
    0.1965723000, 0.1967196000, 0.1968666000, 0.1970134000, 0.1971600000, 
    0.1973063000, 0.1974524000, 0.1975983000, 0.1977440000, 0.1978894000, 
    0.1980347000, 0.1981797000, 0.1983245000, 0.1984690000, 0.1986134000, 
    0.1987575000, 0.1989015000, 0.1990452000, 0.1991886000, 0.1993319000, 
    0.1994750000, 0.1996178000, 0.1997604000, 0.1999029000, 0.2000451000, 
    0.2001871000, 0.2003289000, 0.2004704000, 0.2006118000, 0.2007530000, 
    0.2008939000, 0.2010347000, 0.2011752000, 0.2013156000, 0.2014557000, 
    0.2015956000, 0.2017353000, 0.2018749000, 0.2020142000, 0.2021533000, 
    0.2022922000, 0.2024310000, 0.2025695000, 0.2027078000, 0.2028459000, 
    0.2029839000, 0.2031216000, 0.2032591000, 0.2033965000, 0.2035336000, 
    0.2036706000, 0.2038073000, 0.2039439000, 0.2040803000, 0.2042165000, 
    0.2043524000, 0.2044883000, 0.2046239000, 0.2047593000, 0.2048945000, 
    0.2050296000, 0.2051644000, 0.2052991000, 0.2054336000, 0.2055679000, 
    0.2057020000, 0.2058360000, 0.2059697000, 0.2061033000, 0.2062367000, 
    0.2063699000, 0.2065029000, 0.2066358000, 0.2067684000, 0.2069009000, 
    0.2070332000, 0.2071653000, 0.2072973000, 0.2074291000, 0.2075607000, 
    0.2076921000, 0.2078233000, 0.2079544000, 0.2080853000, 0.2082161000, 
    0.2083466000, 0.2084770000, 0.2086072000, 0.2087372000, 0.2088671000, 
    0.2089968000, 0.2091263000, 0.2092557000, 0.2093849000, 0.2095139000, 
    0.2096428000, 0.2097715000, 0.2099000000, 0.2100284000, 0.2101566000, 
    0.2102846000, 0.2104125000, 0.2105402000, 0.2106677000, 0.2107951000, 
    0.2109223000, 0.2110494000, 0.2111763000, 0.2113030000, 0.2114296000, 
    0.2115560000, 0.2116823000, 0.2118084000, 0.2119343000, 0.2120601000, 
    0.2121858000, 0.2123112000, 0.2124366000, 0.2125617000, 0.2126868000, 
    0.2128116000, 0.2129363000, 0.2130609000, 0.2131853000, 0.2133095000, 
    0.2134337000, 0.2135576000, 0.2136814000, 0.2138051000, 0.2139286000, 
    0.2140519000, 0.2141751000, 0.2142982000, 0.2144211000, 0.2145439000, 
    0.2146665000, 0.2147890000, 0.2149113000, 0.2150335000, 0.2151555000, 
    0.2152774000, 0.2153992000, 0.2155208000, 0.2156423000, 0.2157636000, 
    0.2158848000, 0.2160059000, 0.2161268000, 0.2162475000, 0.2163682000, 
    0.2164887000, 0.2166090000, 0.2167292000, 0.2168493000, 0.2169693000, 
    0.2170891000, 0.2172088000, 0.2173283000, 0.2174477000, 0.2175670000, 
    0.2176861000, 0.2178051000, 0.2179240000, 0.2180427000, 0.2181613000, 
    0.2182798000, 0.2183982000, 0.2185164000, 0.2186345000, 0.2187524000, 
    0.2188703000, 0.2189880000, 0.2191055000, 0.2192230000, 0.2193403000, 
    0.2194575000, 0.2195745000, 0.2196915000, 0.2198083000, 0.2199250000, 
    0.2200416000, 0.2201580000, 0.2202743000, 0.2203905000, 0.2205066000, 
    0.2206225000, 0.2207384000, 0.2208541000, 0.2209697000, 0.2210851000, 
    0.2212005000, 0.2213157000, 0.2214308000, 0.2215458000, 0.2216607000, 
    0.2217754000, 0.2218901000, 0.2220046000, 0.2221190000, 0.2222333000, 
    0.2223475000, 0.2224615000, 0.2225755000, 0.2226893000, 0.2228030000, 
    0.2229166000, 0.2230301000, 0.2231435000, 0.2232568000, 0.2233699000, 
    0.2234830000, 0.2235959000, 0.2237088000, 0.2238215000, 0.2239341000, 
    0.2240466000, 0.2241590000, 0.2242713000, 0.2243834000, 0.2244955000, 
    0.2246075000, 0.2247193000, 0.2248311000, 0.2249427000, 0.2250543000, 
    0.2251657000, 0.2252770000, 0.2253883000, 0.2254994000, 0.2256104000, 
    0.2257213000, 0.2258321000, 0.2259429000, 0.2260535000, 0.2261640000, 
    0.2262744000, 0.2263847000, 0.2264949000, 0.2266050000, 0.2267151000, 
    0.2268250000, 0.2269348000, 0.2270445000, 0.2271542000, 0.2272637000, 
    0.2273731000, 0.2274825000, 0.2275917000, 0.2277009000, 0.2278099000, 
    0.2279189000, 0.2280277000, 0.2281365000, 0.2282452000, 0.2283538000, 
    0.2284623000, 0.2285707000, 0.2286790000, 0.2287872000, 0.2288954000, 
    0.2290034000, 0.2291114000, 0.2292192000, 0.2293270000, 0.2294347000, 
    0.2295423000, 0.2296498000, 0.2297572000, 0.2298646000, 0.2299718000, 
    0.2300790000, 0.2301861000, 0.2302931000, 0.2304000000, 0.2305068000, 
    0.2306136000, 0.2307202000, 0.2308268000, 0.2309333000, 0.2310397000, 
    0.2311460000, 0.2312523000, 0.2313584000, 0.2314645000, 0.2315705000, 
    0.2316765000, 0.2317823000, 0.2318881000, 0.2319938000, 0.2320994000, 
    0.2322049000, 0.2323104000, 0.2324157000, 0.2325210000, 0.2326263000, 
    0.2327314000, 0.2328365000, 0.2329415000, 0.2330464000, 0.2331513000, 
    0.2332560000, 0.2333607000, 0.2334654000, 0.2335699000, 0.2336744000, 
    0.2337788000, 0.2338831000, 0.2339874000, 0.2340916000, 0.2341957000, 
    0.2342998000, 0.2344038000, 0.2345077000, 0.2346115000, 0.2347153000, 
    0.2348190000, 0.2349227000, 0.2350262000, 0.2351298000, 0.2352332000, 
    0.2353366000, 0.2354399000, 0.2355431000, 0.2356463000, 0.2357494000, 
    0.2358525000, 0.2359555000, 0.2360584000, 0.2361612000, 0.2362640000, 
    0.2363668000, 0.2364694000, 0.2365721000, 0.2366746000, 0.2367771000, 
    0.2368795000, 0.2369819000, 0.2370842000, 0.2371864000, 0.2372886000, 
    0.2373908000, 0.2374928000, 0.2375948000, 0.2376968000, 0.2377987000, 
    0.2379005000, 0.2380023000, 0.2381041000, 0.2382057000, 0.2383074000, 
    0.2384089000, 0.2385104000, 0.2386119000, 0.2387133000, 0.2388146000, 
    0.2389159000, 0.2390172000, 0.2391184000, 0.2392195000, 0.2393206000, 
    0.2394216000, 0.2395226000, 0.2396236000, 0.2397245000, 0.2398253000, 
    0.2399261000, 0.2400268000, 0.2401275000, 0.2402281000, 0.2403287000, 
    0.2404293000, 0.2405298000, 0.2406302000, 0.2407306000, 0.2408310000, 
    0.2409313000, 0.2410316000, 0.2411318000, 0.2412320000, 0.2413321000, 
    0.2414322000, 0.2415323000, 0.2416323000, 0.2417323000, 0.2418322000, 
    0.2419321000, 0.2420319000, 0.2421317000, 0.2422315000, 0.2423312000, 
    0.2424309000, 0.2425305000, 0.2426301000, 0.2427297000, 0.2428292000, 
    0.2429287000, 0.2430281000, 0.2431275000, 0.2432269000, 0.2433262000, 
    0.2434255000, 0.2435248000, 0.2436240000, 0.2437232000, 0.2438224000, 
    0.2439215000, 0.2440206000, 0.2441196000, 0.2442187000, 0.2443176000, 
    0.2444166000, 0.2445155000, 0.2446144000, 0.2447133000, 0.2448121000, 
    0.2449109000, 0.2450097000, 0.2451084000, 0.2452071000, 0.2453058000, 
    0.2454044000, 0.2455031000, 0.2456017000, 0.2457002000, 0.2457988000, 
    0.2458973000, 0.2459957000, 0.2460942000, 0.2461926000, 0.2462910000, 
    0.2463894000, 0.2464878000, 0.2465861000, 0.2466844000, 0.2467827000, 
    0.2468809000, 0.2469792000, 0.2470774000, 0.2471756000, 0.2472737000, 
    0.2473719000, 0.2474700000, 0.2475681000, 0.2476662000, 0.2477642000, 
    0.2478623000, 0.2479603000, 0.2480583000, 0.2481563000, 0.2482542000, 
    0.2483522000, 0.2484501000, 0.2485480000, 0.2486459000, 0.2487437000, 
    0.2488416000, 0.2489394000, 0.2490372000, 0.2491351000, 0.2492328000, 
    0.2493306000, 0.2494284000, 0.2495261000, 0.2496239000, 0.2497216000, 
    0.2498193000, 0.2499170000, 0.2500146000, 0.2501123000, 0.2502100000, 
    0.2503076000, 0.2504052000, 0.2505029000, 0.2506005000, 0.2506981000, 
    0.2507956000, 0.2508932000, 0.2509908000, 0.2510883000, 0.2511859000, 
    0.2512834000, 0.2513810000, 0.2514785000, 0.2515760000, 0.2516735000, 
    0.2517710000, 0.2518685000, 0.2519660000, 0.2520635000, 0.2521610000, 
    0.2522584000, 0.2523559000, 0.2524534000, 0.2525508000, 0.2526483000, 
    0.2527457000, 0.2528432000, 0.2529406000, 0.2530381000, 0.2531355000, 
    0.2532330000, 0.2533304000, 0.2534278000, 0.2535252000, 0.2536227000, 
    0.2537201000, 0.2538175000, 0.2539150000, 0.2540124000, 0.2541098000, 
    0.2542073000, 0.2543047000, 0.2544021000, 0.2544996000, 0.2545970000, 
    0.2546944000, 0.2547919000, 0.2548893000, 0.2549868000, 0.2550842000, 
    0.2551817000, 0.2552791000, 0.2553766000, 0.2554740000, 0.2555715000, 
    0.2556690000, 0.2557665000, 0.2558639000, 0.2559614000, 0.2560589000, 
    0.2561564000, 0.2562539000, 0.2563515000, 0.2564490000, 0.2565465000, 
    0.2566441000, 0.2567416000, 0.2568392000, 0.2569367000, 0.2570343000, 
    0.2571319000, 0.2572295000, 0.2573271000, 0.2574247000, 0.2575223000, 
    0.2576199000, 0.2577176000, 0.2578152000, 0.2579129000, 0.2580106000, 
    0.2581082000, 0.2582059000, 0.2583036000, 0.2584014000, 0.2584991000, 
    0.2585968000, 0.2586946000, 0.2587924000, 0.2588902000, 0.2589879000, 
    0.2590858000, 0.2591836000, 0.2592814000, 0.2593793000, 0.2594771000, 
    0.2595750000, 0.2596729000, 0.2597708000, 0.2598688000, 0.2599667000, 
    0.2600647000, 0.2601627000, 0.2602606000, 0.2603587000, 0.2604567000, 
    0.2605547000, 0.2606528000, 0.2607509000, 0.2608490000, 0.2609471000, 
    0.2610452000, 0.2611434000, 0.2612415000, 0.2613397000, 0.2614379000, 
    0.2615361000, 0.2616344000, 0.2617326000, 0.2618309000, 0.2619292000, 
    0.2620275000, 0.2621259000, 0.2622242000, 0.2623226000, 0.2624210000, 
    0.2625194000, 0.2626179000, 0.2627164000, 0.2628148000, 0.2629133000, 
    0.2630119000, 0.2631104000, 0.2632090000, 0.2633076000, 0.2634062000, 
    0.2635048000, 0.2636035000, 0.2637022000, 0.2638009000, 0.2638996000, 
    0.2639984000, 0.2640971000, 0.2641959000, 0.2642947000, 0.2643936000, 
    0.2644925000, 0.2645913000, 0.2646903000, 0.2647892000, 0.2648882000, 
    0.2649871000, 0.2650861000, 0.2651852000, 0.2652842000, 0.2653833000, 
    0.2654824000, 0.2655816000, 0.2656807000, 0.2657799000, 0.2658791000, 
    0.2659784000, 0.2660776000, 0.2661769000, 0.2662762000, 0.2663755000, 
    0.2664749000, 0.2665743000, 0.2666737000, 0.2667731000, 0.2668726000, 
    0.2669721000, 0.2670716000, 0.2671712000, 0.2672707000, 0.2673703000, 
    0.2674700000, 0.2675696000, 0.2676693000, 0.2677690000, 0.2678687000, 
    0.2679685000, 0.2680683000, 0.2681681000, 0.2682679000, 0.2683678000, 
    0.2684677000, 0.2685676000, 0.2686675000, 0.2687675000, 0.2688675000, 
    0.2689676000, 0.2690676000, 0.2691677000, 0.2692678000, 0.2693680000, 
    0.2694681000, 0.2695683000, 0.2696685000, 0.2697688000, 0.2698691000, 
    0.2699694000, 0.2700697000, 0.2701701000, 0.2702705000, 0.2703709000, 
    0.2704713000, 0.2705718000, 0.2706723000, 0.2707729000, 0.2708734000, 
    0.2709740000, 0.2710746000, 0.2711753000, 0.2712759000, 0.2713766000, 
    0.2714774000, 0.2715781000, 0.2716789000, 0.2717797000, 0.2718806000, 
    0.2719814000, 0.2720823000, 0.2721833000, 0.2722842000, 0.2723852000, 
    0.2724862000, 0.2725873000, 0.2726883000, 0.2727894000, 0.2728905000, 
    0.2729917000, 0.2730929000, 0.2731941000, 0.2732953000, 0.2733966000, 
    0.2734979000, 0.2735992000, 0.2737005000, 0.2738019000, 0.2739033000, 
    0.2740048000, 0.2741062000, 0.2742077000, 0.2743092000, 0.2744108000, 
    0.2745123000, 0.2746139000, 0.2747156000, 0.2748172000, 0.2749189000, 
    0.2750206000, 0.2751223000, 0.2752241000, 0.2753259000, 0.2754277000, 
    0.2755295000, 0.2756314000, 0.2757333000, 0.2758352000, 0.2759372000, 
    0.2760392000, 0.2761412000, 0.2762432000, 0.2763453000, 0.2764474000, 
    0.2765495000, 0.2766516000, 0.2767538000, 0.2768560000, 0.2769582000, 
    0.2770604000, 0.2771627000, 0.2772650000, 0.2773673000, 0.2774697000, 
    0.2775720000, 0.2776745000, 0.2777769000, 0.2778793000, 0.2779818000, 
    0.2780843000, 0.2781868000, 0.2782894000, 0.2783920000, 0.2784946000, 
    0.2785972000, 0.2786999000, 0.2788026000, 0.2789053000, 0.2790080000, 
    0.2791107000, 0.2792135000, 0.2793163000, 0.2794192000, 0.2795220000, 
    0.2796249000, 0.2797278000, 0.2798307000, 0.2799337000, 0.2800366000, 
    0.2801396000, 0.2802427000, 0.2803457000, 0.2804488000, 0.2805519000, 
    0.2806550000, 0.2807581000, 0.2808613000, 0.2809645000, 0.2810677000, 
    0.2811709000, 0.2812742000, 0.2813774000, 0.2814807000, 0.2815841000, 
    0.2816874000, 0.2817908000, 0.2818942000, 0.2819976000, 0.2821010000, 
    0.2822044000, 0.2823079000, 0.2824114000, 0.2825149000, 0.2826185000, 
    0.2827220000, 0.2828256000, 0.2829292000, 0.2830328000, 0.2831365000, 
    0.2832401000, 0.2833438000, 0.2834475000, 0.2835512000, 0.2836550000, 
    0.2837588000, 0.2838625000, 0.2839663000, 0.2840702000, 0.2841740000, 
    0.2842779000, 0.2843817000, 0.2844856000, 0.2845896000, 0.2846935000, 
    0.2847975000, 0.2849014000, 0.2850054000, 0.2851094000, 0.2852135000, 
    0.2853175000, 0.2854216000, 0.2855257000, 0.2856298000, 0.2857339000, 
    0.2858380000, 0.2859422000, 0.2860464000, 0.2861505000, 0.2862548000, 
    0.2863590000, 0.2864632000, 0.2865675000, 0.2866717000, 0.2867760000, 
    0.2868803000, 0.2869847000, 0.2870890000, 0.2871933000, 0.2872977000, 
    0.2874021000, 0.2875065000, 0.2876109000, 0.2877153000, 0.2878198000, 
    0.2879242000, 0.2880287000, 0.2881332000, 0.2882377000, 0.2883422000, 
    0.2884468000, 0.2885513000, 0.2886559000, 0.2887604000, 0.2888650000, 
    0.2889696000, 0.2890742000, 0.2891789000, 0.2892835000, 0.2893882000, 
    0.2894928000, 0.2895975000, 0.2897022000, 0.2898069000, 0.2899116000, 
    0.2900163000, 0.2901211000, 0.2902258000, 0.2903306000, 0.2904354000, 
    0.2905402000, 0.2906449000, 0.2907498000, 0.2908546000, 0.2909594000, 
    0.2910642000, 0.2911691000, 0.2912740000, 0.2913788000, 0.2914837000 
};
 
void EvtDToKmPipPipPi0::init(){
    // check that there are 0 arguments
    checkNArg(0);
    checkNDaug(4);
 
    checkSpinParent(EvtSpinType::SCALAR);
    checkSpinDaughter(0,EvtSpinType::SCALAR);
    checkSpinDaughter(1,EvtSpinType::SCALAR);
    checkSpinDaughter(2,EvtSpinType::SCALAR);
    checkSpinDaughter(3,EvtSpinType::SCALAR);
 
//------------------parameters-----------------
    rho[0] =  1.00;                    phi[0] =  0.00;      
    rho[1] =  1.6407e-01;              phi[1] =  1.4629e+00;
    rho[2] =  5.3805e-01;              phi[2] =  5.6802e-02;
    rho[3] =  9.8491e-01;              phi[3] =  4.5297e-01;
        rho[4] = -7.5361e-01;              phi[4] =  3.9923e+00;
    rho[5] =  7.2022e-01;              phi[5] =  3.6276e-01;
    rho[6] =  1.5187e+00;              phi[6] =  1.1700e+00;
    rho[7] =  3.0364e+00;              phi[7] =  2.9164e+00;
    rho[8] =  5.2613e+00;              phi[8] =  5.0168e+00;
    rho[9] =  1.1521e+01;              phi[9] =  3.8225e+00;
    rho[10]= -7.3292e+00;              phi[10]=  4.8851e+00;
    rho[11]= -8.9717e-01;              phi[11]= -2.7156e+00;
    rho[12]=  7.9479e-01;              phi[12]= -1.9858e+00;
    rho[13]=  4.6589e+00;              phi[13]= -8.5107e-01;
    rho[14]=  6.8280e-01;              phi[14]= -2.1226e+00;
    mD   = 1.86486;
//------------------------new----------------------
    mKst0  = 0.89555;
    mrho   = 0.77511;
    mrho1450   = 1.465;
    mK1270 = 1.289;
    mK1400 = 1.403;
    mK1460 = 1.482;
    mK1650 = 1.672;
    mK1680 = 1.718;
 
        GKst0  = 0.0473;
        Grho   = 0.1491;
        Grho1450   = 0.400;
        GK1270 = 0.116;
        GK1400 = 0.174;
        GK1460 = 0.3356;
        GK1650 = 0.158;
        GK1680 = 0.322;
 
    mass_Pion = 0.13957;
    mass_Pion_N = 0.134977;
    mass_Eta = 0.547862;
    mass_Kaon = 0.493677;
    math_pi = 3.1415926;
 
    rD2 = 5.0;    // 5*5
    rRes1 = 3.0;   // 3*3
    rRes2 = 3.0;   // 3*3
 
    GS1 = 0.636619783;
    GS2 = 0.01860182466;
    GS3 = 0.1591549458;   //  1/(2*math_2pi)
    GS4 = 0.00620060822;  //  mass_Pion2/math_pi
 
    int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
    int EE[4][4][4][4] =
    { { {{0,0,0,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},  {0,0,-1,0}},
          {{0,0,0,0},  {0,0,0,-1}, {0,0,0,0},  {0,1,0,0} },
          {{0,0,0,0},  {0,0,1,0},  {0,-1,0,0}, {0,0,0,0} } },
    { {{0,0,0,0},  {0,0,0,0},  {0,0,0,-1}, {0,0,1,0} },
        {{0,0,0,0},  {0,0,0,0},  {0,0,0,0},  {0,0,0,0} },
        {{0,0,0,1},  {0,0,0,0},  {0,0,0,0},  {-1,0,0,0}},
        {{0,0,-1,0}, {0,0,0,0},  {1,0,0,0},  {0,0,0,0} } },
    { {{0,0,0,0},  {0,0,0,1},  {0,0,0,0},  {0,-1,0,0}},
        {{0,0,0,-1}, {0,0,0,0},  {0,0,0,0},  {1,0,0,0} },
        {{0,0,0,0},  {0,0,0,0},  {0,0,0,0},  {0,0,0,0} },
        {{0,1,0,0},  {-1,0,0,0}, {0,0,0,0},  {0,0,0,0} } },
    { {{0,0,0,0},  {0,0,-1,0}, {0,1,0,0},  {0,0,0,0} },
        {{0,0,1,0},  {0,0,0,0},  {-1,0,0,0}, {0,0,0,0} },
        {{0,-1,0,0}, {1,0,0,0},  {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<4; i++) {
        for (int j=0; j<4; j++) {
            G[i][j] = GG[i][j];
            for (int k=0; k<4; k++) {
                for (int l=0; l<4; l++) {
                    E[i][j][k][l] = EE[i][j][k][l];
                }
            }
        }
    }
}
 
 
 
 
void EvtDToKmPipPipPi0::initProbMax() {
    setProbMax(40000);
}
 
void EvtDToKmPipPipPi0::decay( EvtParticle *p){
//--------------max
/*  
       double maxprob=0,maxq1270;
       for(int ir=0;ir<=10000000;ir++){
          p->initializePhaseSpace(getNDaug(),getDaugs());
             EvtVector4R km   = p->getDaug(0)->getP4();
             EvtVector4R pip1 = p->getDaug(1)->getP4();
             EvtVector4R pip2 = p->getDaug(2)->getP4();
             EvtVector4R pi0  = p->getDaug(3)->getP4();
             int num = EvtPDL::getStdHep(p->getId());
             
        double Pip1[4],Pip2[4],Km[4],Pi0[4];
        Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
        Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
        Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
        Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
 
        if(num>0){
        Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
        Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
        Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
        Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
}
        if(num<0){
        Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
        Km[1] = -km.get(1); Pip1[1] = -pip1.get(1); Pip2[1] = -pip2.get(1);Pi0[1] = -pi0.get(1);
        Km[2] = -km.get(2); Pip1[2] = -pip1.get(2); Pip2[2] = -pip2.get(2);Pi0[2] = -pi0.get(2);
        Km[3] = -km.get(3); Pip1[3] = -pip1.get(3); Pip2[3] = -pip2.get(3);Pi0[3] = -pi0.get(3);
}
 
    double _prob;
        double value;
        int nstates=15;
        int modetype[15]= {1,1,8,4,4,12,11,17,13,12,8,4,9,13,1};
        int g0[15] = {1,1,1,1,1,1,1,1,0,0,1,1,1,0,0};
        int g1[15] = {1,1,1,1,1,1,1,1,1,1,0,0,0,0,1};
        int g2[15] = {0,1,0,0,2,1,1,0,1,1,2,0,0,0,0};
        double mass1[15]  = {mKst0,mKst0,mrho,  mKst0, mKst0, mKst0,  mKst0, mKst0, mrho,  mKst0, mrho,  mKst0,  mKst0, mrho,mKst0};
        double mass2[15]  = {mrho, mrho, mK1270,mK1400,mK1400,mK1460,mK1680, mrho,  mK1460,mK1460,mK1270,mK1400,mK1270, mK1460,mrho};
        double width1[15] = {GKst0,GKst0,Grho,  GKst0, GKst0, GKst0,  GKst0, GKst0, Grho,  GKst0, Grho,  GKst0,  GKst0, Grho,GKst0};
        double width2[15] = {Grho, Grho, GK1270,GK1400,GK1400,GK1460,GK1680, Grho,  GK1460,GK1460,GK1270,GK1400,GK1270, GK1460,Grho};
        calEvaMy(Km,Pip1,Pip2,Pi0,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,nstates,_prob);
    
       if(_prob>maxprob) {
       maxprob=_prob;
       }
    value = _prob;
        setProb(value);
       }
       printf("maxprob = %.10f\n", maxprob);
        return;
*/
//--------------------main--------------
 
       p->initializePhaseSpace(getNDaug(),getDaugs());
         EvtVector4R km   = p->getDaug(0)->getP4();
         EvtVector4R pip1 = p->getDaug(1)->getP4();
         EvtVector4R pip2 = p->getDaug(2)->getP4();
         EvtVector4R pi0  = p->getDaug(3)->getP4();
         int num = EvtPDL::getStdHep(p->getId()); 
         
    double Pip1[4],Pip2[4],Km[4],Pi0[4];    
    //----------------------------------------------------------------------------------------
    
    if(num>0){
    Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
    Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
    Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
    Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
}
        if(num<0){
        Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
        Km[1] = -km.get(1); Pip1[1] = -pip1.get(1); Pip2[1] = -pip2.get(1);Pi0[1] = -pi0.get(1);
        Km[2] = -km.get(2); Pip1[2] = -pip1.get(2); Pip2[2] = -pip2.get(2);Pi0[2] = -pi0.get(2);
        Km[3] = -km.get(3); Pip1[3] = -pip1.get(3); Pip2[3] = -pip2.get(3);Pi0[3] = -pi0.get(3);
}
 
 
    double value;
        int nstates=15;
        int modetype[15]= {1,1,8,4,4,12,11,17,13,12,8,4,9,13,1};
        int g0[15] = {1,1,1,1,1,1,1,1,0,0,1,1,1,0,0};
        int g1[15] = {1,1,1,1,1,1,1,1,1,1,0,0,0,0,1};
        int g2[15] = {0,1,0,0,2,1,1,0,1,1,2,0,0,0,0};
        double mass1[15]  = {mKst0,mKst0,mrho,  mKst0, mKst0, mKst0,  mKst0, mKst0, mrho,  mKst0, mrho,  mKst0,  mKst0, mrho,mKst0};
        double mass2[15]  = {mrho, mrho, mK1270,mK1400,mK1400,mK1460,mK1680, mrho,  mK1460,mK1460,mK1270,mK1400,mK1270, mK1460,mrho};
        double width1[15] = {GKst0,GKst0,Grho,  GKst0, GKst0, GKst0,  GKst0, GKst0, Grho,  GKst0, Grho,  GKst0,  GKst0, Grho,GKst0};
        double width2[15] = {Grho, Grho, GK1270,GK1400,GK1400,GK1460,GK1680, Grho,  GK1460,GK1460,GK1270,GK1400,GK1270, GK1460,Grho};
 
    calEvaMy(Km,Pip1,Pip2,Pi0,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,nstates,value);
    setProb(value);
 
    return;
 
}
 
void EvtDToKmPipPipPi0::KPiSLASS(double sa, double sb, double sc, double prop[2]) {
   const double m1430 = 1.441;
   const double sa0   = 2.076481;   // m1430*m1430;
   const double w1430 = 0.193;
   const double Lass1 = 0.25/sa0;
   double tmp = sb-sc;
   double tmp1 = sa0+tmp;
   double q0 = Lass1*tmp1*tmp1-sb;
   if(q0<0) q0 = 1e-16;
   double tmp2 = sa+tmp;
   double qs = 0.25*tmp2*tmp2/sa-sb;
   double q = sqrt(qs);
   double width = w1430*q*m1430/sqrt(sa*q0);
   double temp_R = atan(m1430*width/(sa0-sa));
   if(temp_R<0) temp_R += math_pi;
   double deltaR = -1.915 + temp_R;
   double temp_F = atan(0.226*q/(2.0-3.8194*qs));  // 2.0*1.07 = 2.14; 1.8*1.07 = 1.926
   if(temp_F<0) temp_F += math_pi;
   double deltaF = 0.002 + temp_F;
   double deltaS = deltaR + 2.0*deltaF;
   double t1 = 0.96*sin(deltaF);
   double t2 = sin(deltaR);
   double CF[2], CS[2];
   CF[0] = cos(deltaF);
   CF[1] = sin(deltaF);
   CS[0] = cos(deltaS);
   CS[1] = sin(deltaS);
   prop[0] = t1*CF[0] + t2*CS[0];
   prop[1] = t1*CF[1] + t2*CS[1];
}
 
 
void EvtDToKmPipPipPi0::Com_Multi(double a1[2], double a2[2], double res[2]){
    res[0] = a1[0]*a2[0]-a1[1]*a2[1];
    res[1] = a1[1]*a2[0]+a1[0]*a2[1];
}
void EvtDToKmPipPipPi0::Com_Divide(double a1[2], double a2[2], double res[2]){
    res[0] = (a1[0]*a2[0]+a1[1]*a2[1])/(a2[0]*a2[0]+a2[1]*a2[1]);
    res[1] = (a1[1]*a2[0]-a1[0]*a2[1])/(a2[0]*a2[0]+a2[1]*a2[1]);
}
double EvtDToKmPipPipPi0::SCADot(double a1[4], double a2[4])
{
    double _cal = a1[0]*a2[0]-a1[1]*a2[1]-a1[2]*a2[2]-a1[3]*a2[3];
    return _cal;
}
double EvtDToKmPipPipPi0::Barrier(int l, double sa, double sb, double sc, double r,double mass)
{
    double q = (sa+sb-sc)*(sa+sb-sc)/(4*sa)-sb;
//  if(q < 0) q = 1e-16;
    if(q < 0) q = -q;
 
    double z;
    z = q*r*r;
    double sa0;
    sa0 = mass*mass;
    double q0 = (sa0+sb-sc)*(sa0+sb-sc)/(4*sa0)-sb;
//  if(q0 < 0) q0 = 1e-16;
    if(q0 < 0) q0 = -q0;
    double z0 = q0*r*r;
    double F = 0.0;
    if(l == 0) F = 1;
    if(l == 1) F = sqrt((1+z0)/(1+z));
    if(l == 2) F = sqrt((9+3*z0+z0*z0)/(9+3*z+z*z));
    //   if(l == 1) F = sqrt((2*r)/(1+z));
    //   if(l == 2) F = sqrt((13*r*r)/(9+3*z+z*z));
    return F;
}
void EvtDToKmPipPipPi0::calt1(double daug1[4], double daug0[4], double t1[4]){
    double p, pq;
    double pa[4], qa[4];
    for(int i=0; i<4; i++){
        pa[i] = daug1[i] + daug0[i];
        qa[i] = daug1[i] - daug0[i];
    }
    p = SCADot(pa,pa);
    pq = SCADot(pa,qa);
    for(int i=0; i<4; i++){
        t1[i] = qa[i] - pq/p*pa[i];
    }
}
void EvtDToKmPipPipPi0::calt2(double daug1[4], double daug0[4], double t2[4][4]){
    double p, r;
    double pa[4], t1[4];
    calt1(daug1,daug0,t1);
    r = SCADot(t1,t1)/3.0;
    for(int i=0; i<4; i++) {
        pa[i] = daug1[i] + daug0[i];
    }
    p = SCADot(pa,pa);
    for(int i=0; i<4; i++) {
        for(int j=0; j<4; j++) {
            t2[i][j] = t1[i]*t1[j] - r*(G[i][j]-pa[i]*pa[j]/p);
        }
    }
}
 
double EvtDToKmPipPipPi0::wid(double mass2, double mass, double sa, double sb, double sc, double r, int l){
    double widm = 0.;
    double m = sqrt(sa);
    double tmp  = sb-sc;
    double tmp1 = sa+tmp;
    double q = 0.25*tmp1*tmp1/sa-sb;
//  if(q<0) q = 1e-16;
    if(q<0) q = -q;
    double tmp2 = mass2+tmp;
    double q0 = 0.25*tmp2*tmp2/mass2-sb;
//  if(q0<0) q0 = 1e-16;
    if(q0<0) q0 = -q0;
    double z = q*r*r;
    double z0 = q0*r*r;
    double t = q/q0;
    if(l == 0)      {widm = sqrt(t)*mass/m;}
    else if(l == 1) {widm = t*sqrt(t)*mass/m*(1+z0)/(1+z);}
    else if(l == 2) {widm = t*t*sqrt(t)*mass/m*(9+3*z0+z0*z0)/(9+3*z+z*z);}
    return widm;
 
}
double EvtDToKmPipPipPi0::widl1(double mass2, double mass, double sa, double sb, double sc, double r)
{
    double widm = 0.;
    double m = sqrt(sa);
    double tmp  = sb-sc;
    double tmp1 = sa+tmp;
    double q = 0.25*tmp1*tmp1/sa-sb;
//  if(q<0) q = 1e-16;
    if(q<0) q = -q;
    double tmp2 = mass2+tmp;
    double q0 = 0.25*tmp2*tmp2/mass2-sb;
//  if(q0<0) q0 = 1e-16;
    if(q0<0) q0 = -q0;
    double z = q*r*r;
    double z0 = q0*r*r;
    double F = (1+z0)/(1+z);
    double t = q/q0;
    widm = t*sqrt(t)*mass/m*F;
    return widm;
 
}
void EvtDToKmPipPipPi0::propagatorNBW(double mass2, double mass, double width, double sa, double sb, double sc, double r, int l, double prop[2]){
   double a[2], b[2];
   a[0] = 1;
   a[1] = 0;
   b[0] = mass2-sa;
   b[1] = -mass*width;
   Com_Divide(a,b,prop);
}
 
void EvtDToKmPipPipPi0::propagatorRBW(double mass2, double mass, double width, double sa, double sb, double sc, double r, int l, double prop[2]){
    double a[2], b[2];
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sa;
    b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r,l);
    Com_Divide(a,b,prop);
}
void EvtDToKmPipPipPi0::propagatorRBW_K1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2])
{
    double a[2], b[2];
    int iii=int(sa*1000/2)-300;//(4000-0)/2000=2
    if(iii<0){ iii=0; }
    double width_K1;
    if(mass<1.35){
        width_K1=K1270_width[iii];
    }else{
        width_K1=K1400_width[iii];
    }
    //printf("width_K1: %.8f\n", width_K1);
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sa;
    //b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l);
    b[1] = -mass*width_K1;
    Com_Divide(a,b,prop);
}
void EvtDToKmPipPipPi0::propagatorRBWl1(double mass2, double mass, double width, double sa, double sb, double sc, double r, double prop[2])
    //firstly code for propagator_omg, now could use for propagatorRBW that l=1, and will be more quick.
{
    double a[2], b[2];
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sa;
    b[1] = -mass*width*widl1(mass2,mass,sa,sb,sc,r);
    Com_Divide(a,b,prop);
}
//------------GS---used by rho----------------------------
void EvtDToKmPipPipPi0::propagatorGS(double mass2, double mass, double width, double sa, double sb, double sc, double r, double prop[2])
{
    double a[2], b[2];
    double tmp  = sb-sc;
    double tmp1 = sa+tmp;
    double q2 = 0.25*tmp1*tmp1/sa-sb;
//  if(q2<0) q2 = 1e-16;
    if(q2<0) q2 = -q2;
 
    double tmp2 = mass2+tmp;
    double q02 = 0.25*tmp2*tmp2/mass2-sb;
//  if(q02<0) q02 = 1e-16;
    if(q02<0) q02 = -q02;
 
    double q  = sqrt(q2);
    double q0 = sqrt(q02);
    double m  = sqrt(sa);
    double q03 = q0*q02;
    double tmp3 = log(mass+2*q0)+1.2760418309;  // log(mass_2Pion) = 1.2760418309;
 
    double h  = GS1*q/m*(log(m+2*q)+1.2760418309);
    double h0 = GS1*q0/mass*tmp3;
    double dh = h0*(0.125/q02-0.5/mass2)+GS3/mass2;
    double d  = GS2/q02*tmp3+GS3*mass/q0-GS4*mass/q03;
    double f  = mass2/q03*(q2*(h-h0)+(mass2-sa)*q02*dh);
 
    a[0] = 1.0+d*width/mass;
    a[1] = 0.0;
    b[0] = mass2-sa+width*f;
    b[1] = -mass*width*widl1(mass2,mass,sa,sb,sc,r);
    Com_Divide(a,b,prop);
}                  
void EvtDToKmPipPipPi0::calEvaMy(
        double* Km,        double* Pip1,     double* Pip2,  double* Pi0,
        double* mass1,     double* mass2,
        double* width1,    double* width2,
        double* amp,       double* phase,
        int* g0,           int*    g1,    int*    g2,
        int* modetype,     int nstates,   double & Result)
 
{
 
    double pKstr1[4],pKstr2[4],pKstr3[4],prho1[4],prho2[4],pK11[4],pK12[4],pK13[4],pA1[4],pD[4];
    double Amp_KPiS1[2], Amp_PiPiS1[2], Amp_KPiS2[2], Amp_PiPiS2[2];
 
    for(int i=0;i!=4;i++){
        pKstr1[i]=Km[i]+Pip1[i];
        pKstr2[i]=Km[i]+Pip2[i];
        pKstr3[i]=Km[i]+Pi0[i];
        prho1[i]=Pip1[i]+Pi0[i];
        prho2[i]=Pip2[i]+Pi0[i];
        pK11[i]=Km[i]+prho1[i];
        pK12[i]=Km[i]+prho2[i];
        pK13[i]=Km[i]+Pip1[i]+Pip2[i];
        pD[i]=Km[i]+prho1[i]+Pip2[i];
    }
    double skstr1,skstr2;
    double skstr3,srho1,srho2,sk11,sk12,sk13,sD;
    double skaon,spion1,spion2,spi0;
    skaon  = SCADot(Km,Km);
    spion1 = SCADot(Pip1,Pip1);
    spion2 = SCADot(Pip2,Pip2);
    spi0   = SCADot(Pi0,Pi0);
    skstr1 = SCADot(pKstr1,pKstr1);
    skstr2 = SCADot(pKstr2,pKstr2);
    skstr3 = SCADot(pKstr3,pKstr3);
    srho1  = SCADot(prho1,prho1);
    srho2  = SCADot(prho2,prho2);
    sk11   = SCADot(pK11,pK11);
    sk12   = SCADot(pK12,pK12);
    sk13   = SCADot(pK13,pK13);
    sD     = SCADot(pD,pD);
    double t1rho1[4],t1rho2[4],t1kstr3[4];
//------------------------
    double t1kstr1[4],t1kstr2[4], t2k21[4][4], t2k22[4][4], t2k31[4][4], t2k32[4][4], t2k11[4][4], t2k12[4][4], t2k13[4][4], t2k14[4][4];
    double t2rho1[4][4], t2rho2[4][4], t2kstr3[4][4], t2kstr1[4][4], t2kstr2[4][4];
    double t1K11[4],t1K12[4],t1K13[4],t1K14[4];
    double t2kk1[4][4],t2kk2[4][4];
 
//--------------------------
    calt1(Pip1,Pi0,t1rho1);
    calt1(Pip2,Pi0,t1rho2);
    calt1(Km,Pip1,t1kstr1);
    calt1(Km,Pip2,t1kstr2);
    calt1(Km,Pi0 ,t1kstr3);
    calt1(pKstr1,Pi0 ,t1K11);
    calt1(pKstr2,Pi0 ,t1K12);
    calt1(pKstr3,Pip1 ,t1K13);
    calt1(pKstr3,Pip2 ,t1K14);
    calt2(Pip2,Pi0,t2rho1);
    calt2(Pip2,Pi0,t2rho2);
    calt2(Km,Pi0,t2kstr3);
    calt2(Km,Pip1,t2kstr1);
    calt2(Km,Pip2,t2kstr2);
    calt2(pKstr3,Pip1,t2k11);
    calt2(pKstr3,Pip2,t2k12);
    calt2(prho1,Km,t2k21);
    calt2(prho2,Km,t2k22);
    calt2(pKstr1,Pi0,t2k31);
    calt2(pKstr2,Pi0,t2k32);
    calt2(pKstr1,Pip2,t2kk1);
    calt2(pKstr2,Pip1,t2kk2);
//----------------------------------
        double cof[2],amp_tmp1[2],amp_PDF[2], PDF[2];
        double amp_tmp2[2];
        double amp_tmp[2];
 
    PDF[0] = 0.0;
    PDF[1] = 0.0;
 
    double temp_PDF1, temp_PDF2, tt1, tt2, tmp1, tmp2;
 
    double temp_PDF3, temp_PDF4;
 
    double  pro_1[2], pro_2[2], pro0_1[2], pro1_1[2], pro0_2[2], pro1_2[2], temp_1[2], temp_2[2];
 
    double t1D[4], t1D_1[4], t1D_2[4];
 
    double t1K1_1[4], t1K1_2[4],t1K1_3[4],t1K1_4[4];
 
    double B1_rho1=-1.0, B1_rho2=-1.0;
    double B1_D_K1_Pi2=-1.0, B1_D_K1_Pi1=-1.0, B2_D_K1_Pi2=-1.0, B2_D_K1_Pi1=-1.0;
    double B1_K1_rho1=-1.0, B1_K1_rho2=-1.0, B2_K1_rho1=-1.0, B2_K1_rho2=-1.0;
    double B1_K1_Kstr3_Pi1=-1.0, B1_K1_Kstr3_Pi2=-1.0, B2_K1_Kstr3_Pi1=-1.0, B2_K1_Kstr3_Pi2=-1.0;
    double B1_K1_Kstr1_Pi0=-1.0, B1_K1_Kstr2_Pi0=-1.0;
    double B1_kstr3=-1.0;
    double B2_K1p_kstr1=-1.0,B2_K1p_kstr2=-1.0;
    double B1_K1p_kstr1=-1.0,B1_K1p_kstr2=-1.0;
    double B1_kstr1_Pi1=-1.0, B1_kstr2_Pi2=-1.0, B1_D_Kstr2_rho1=-1.0, B1_D_Kstr1_rho2=-1.0,B1_D_K1p_Pi0=-1.0;
    double B2_K1_Kstr1=-1.0, B2_K1_Kstr2=-1.0, B2_D_Kstr1_rho2=-1.0, B2_D_Kstr2_rho1=-1.0;
    double B2_rho1=-1.0, B2_rho2=-1.0;
    double B2_kstr1_Pi1=-1.0, B2_kstr2_Pi2=-1.0, B2_kstr3_Pi0=-1.0;;
    double mass1sq, mass2sq;
    int isKst=0, isRho=0, isPiPi_S=0, isKPi_S=0;
    double proRho1[2], proRho2[2], proPiPi1_S[2], proPiPi2_S[2], proKPi_S[2];
 
    double tmp3, tmp4;
    double t2D[4][4], t2D_1[4][4], t2D_2[4][4];
    double pro2_1[2], pro2_2[2], pro3_1[2], pro3_2[2], pro4_1[2], pro4_2[2];
    double t2f2_1[4][4],t2f2_2[4][4];
    double t2G2_1[4][4], t2G2_2[4][4], t2D11_1[4][4],t2D11_2[4][4];
 
//-------------------------------------------------------------------------
    //for(int i=0; i<nstates; i++){
    for(int i=0; i<15; i++){
    //for(int i=2; i<15; i++){
        cof[0] = amp[i]*cos(phase[i]);
        cof[1] = amp[i]*sin(phase[i]);
        mass1sq = mass1[i]*mass1[i];
        mass2sq = mass2[i]*mass2[i];
        temp_PDF1 = 0;
        temp_PDF2 = 0;
        temp_PDF3 = 0;
        temp_PDF4 = 0;
 
        if(modetype[i]==1)
        {
                B1_kstr1_Pi1 = Barrier(1,skstr1,skaon,spion1,rRes2,mass1[i]);
                B1_kstr2_Pi2 = Barrier(1,skstr2,skaon,spion2,rRes2,mass1[i]);
                B1_rho1 = Barrier(1,srho1,spion1,spi0,rRes2,mass2[i]);
                B1_rho2 = Barrier(1,srho2,spion2,spi0,rRes2,mass2[i]);
                if (g2[i]==0)
                {
                        for(int w=0; w<4; w++)
                        {
                                temp_PDF1 += G[w][w]*t1kstr1[w]*t1rho2[w];
                                temp_PDF2 += G[w][w]*t1kstr2[w]*t1rho1[w];
                        }
                        tmp1 = B1_kstr1_Pi1*B1_rho2*temp_PDF1;
                        tmp2 = B1_kstr2_Pi2*B1_rho1*temp_PDF2;
                }
                else if (g2[i]==1)
                {
                        calt1(pKstr1,prho2,t1D_1);
                        calt1(pKstr2,prho1,t1D_2);
                        for(int w=0; w<4; w++)
                        {
                                for(int j=0; j<4; j++)
                                {
                                        for(int k=0; k<4; k++)
                                                {
                                                        for(int l=0; l<4; l++)
                                                        {
                                                                temp_PDF1 += E[w][j][k][l]*pD[w]*t1D_1[j]*t1kstr1[k]*t1rho2[l]*G[w][w]*G[j][j]*G[k][k]*G[l][l];
                                                                temp_PDF2 += E[w][j][k][l]*pD[w]*t1D_2[j]*t1kstr2[k]*t1rho1[l]*G[w][w]*G[j][j]*G[k][k]*G[l][l];
                                                        }
                                                }
                                }
                        }
                        B1_D_Kstr1_rho2 = Barrier(1,sD,skstr1,srho2,rD2,mD);
                        B1_D_Kstr2_rho1 = Barrier(1,sD,skstr2,srho1,rD2,mD);
                        tmp1 = B1_kstr1_Pi1*B1_rho2*B1_D_Kstr1_rho2*temp_PDF1;
                        tmp2 = B1_kstr2_Pi2*B1_rho1*B1_D_Kstr2_rho1*temp_PDF2;
                }
                else if (g2[i]==2)
                {
                        calt2(pKstr1,prho2,t2D_1);
                        calt2(pKstr2,prho1,t2D_2);
                        for(int w=0; w<4; w++)
                                {
                                        for(int j=0; j<4; j++)
                                        {
                                                temp_PDF1 += t2D_1[w][j]*t1kstr1[w]*t1rho2[j]*G[w][w]*G[j][j];
                                                temp_PDF2 += t2D_2[w][j]*t1kstr2[w]*t1rho1[j]*G[w][w]*G[j][j];
                                        }
                                }
                        B2_D_Kstr2_rho1 = Barrier(2,sD,skstr1,srho2,rD2,mD);
                        B2_D_Kstr2_rho1 = Barrier(2,sD,skstr2,srho1,rD2,mD);
                        tmp1 = B1_kstr1_Pi1*B1_rho2*B2_D_Kstr1_rho2*temp_PDF1;
                        tmp2 = B1_kstr2_Pi2*B1_rho1*B2_D_Kstr2_rho1*temp_PDF2;
                }
                if(g0[i]==1){
                        propagatorRBW(mass1sq,mass1[i],width1[i],skstr1,skaon,spion1,rRes2,1,pro0_1);
                        propagatorRBW(mass1sq,mass1[i],width1[i],skstr2,skaon,spion2,rRes2,1,pro0_2);
                }
                else if(g0[i]==0){
                        pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0;
                }
                if(g1[i]==1){
                        propagatorGS(mass2sq,mass2[i],width2[i],srho2,spion2,spi0,rRes2,pro1_1);
                        propagatorGS(mass2sq,mass2[i],width2[i],srho1,spion1,spi0,rRes2,pro1_2);
                }
                else if(g1[i]==0){
                        pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0;
                }
                Com_Multi(pro0_1,pro1_1,pro_1);
                Com_Multi(pro0_2,pro1_2,pro_2);
                amp_tmp1[0] = tmp1*pro_1[0];
                amp_tmp1[1] = tmp1*pro_1[1];
                amp_tmp2[0] = tmp2*pro_2[0];
                amp_tmp2[1] = tmp2*pro_2[1];
}
        else if (modetype[i]==4) {
            B1_D_K1_Pi2 = Barrier(1,sD,sk11,spion2,rD2,mass2[i]);
            B1_D_K1_Pi1 = Barrier(1,sD,sk12,spion1,rD2,mass2[i]);
            B1_kstr3    = Barrier(1,skstr3,skaon,spi0,rRes2,mass1[i]);
            B1_kstr1_Pi1= Barrier(1,skstr1,skaon,spion1,rRes2,mass1[i]);
            B1_kstr2_Pi2= Barrier(1,skstr2,skaon,spion2,rRes2,mass1[i]);
                calt1(pK11,Pip2,t1D_1);
                calt1(pK12,Pip1,t1D_2);
            if (g2[i]==0) {
                for(int w=0; w<4; w++) {
                    for(int j=0; j<4; j++) {
                        tt2 = t1kstr3[j]*G[w][w]*G[j][j];
                        temp_PDF1 += t1D_1[w]*(G[w][j]-pK11[w]*pK11[j]/sk11)*t1kstr1[j]*G[w][w]*G[j][j];
                        temp_PDF2 += t1D_2[w]*(G[w][j]-pK12[w]*pK12[j]/sk12)*t1kstr2[j]*G[w][w]*G[j][j];
                        temp_PDF3 += t1D_1[w]*(G[w][j]-pK11[w]*pK11[j]/sk11)*tt2;
                        temp_PDF4 += t1D_2[w]*(G[w][j]-pK12[w]*pK12[j]/sk12)*tt2;
                    }
                }
                tmp1 = B1_kstr1_Pi1*B1_D_K1_Pi2*temp_PDF1;
                tmp2 = B1_kstr2_Pi2*B1_D_K1_Pi1*temp_PDF2;
                tmp3 = B1_kstr3*B1_D_K1_Pi2*temp_PDF3;
                tmp4 = B1_kstr3*B1_D_K1_Pi1*temp_PDF4;
            }
            else if (g2[i]==2) {
                for(int w=0; w<4; w++) {
                    for(int j=0; j<4; j++) {
                        tt2 = t1kstr3[j]*G[w][w]*G[j][j];
                        temp_PDF1 += t1D_1[w]*t2k31[w][j]*t1kstr1[j]*G[w][w]*G[j][j];
                        temp_PDF2 += t1D_2[w]*t2k32[w][j]*t1kstr2[j]*G[w][w]*G[j][j];
                        temp_PDF3 += t1D_1[w]*t2k11[w][j]*tt2;
                        temp_PDF4 += t1D_2[w]*t2k12[w][j]*tt2;
                    }
                }
                B2_K1_Kstr3_Pi1 = Barrier(2,sk11,skstr3,spion1,rRes2,mass2[i]);
                B2_K1_Kstr3_Pi2 = Barrier(2,sk12,skstr3,spion2,rRes2,mass2[i]);
                B2_K1_Kstr1     = Barrier(2,sk11,skstr1,spi0,rRes2,mass2[i]);
                B2_K1_Kstr2     = Barrier(2,sk12,skstr2,spi0,rRes2,mass2[i]);
                tmp1 = B1_kstr1_Pi1*B2_K1_Kstr1*B1_D_K1_Pi2*temp_PDF1;
                tmp2 = B1_kstr2_Pi2*B2_K1_Kstr2*B1_D_K1_Pi1*temp_PDF2;
                tmp3 = B1_kstr3*B2_K1_Kstr3_Pi1*B1_D_K1_Pi2*temp_PDF3;
                tmp4 = B1_kstr3*B2_K1_Kstr3_Pi2*B1_D_K1_Pi1*temp_PDF4;
            }
            propagatorRBW(mass1sq,mass1[i],width1[i],skstr3,skaon,spi0,rRes2,1,pro2_1);
            propagatorRBW(mass1sq,mass1[i],width1[i],skstr3,skaon,spi0,rRes2,1,pro2_2);
            propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk11,skstr3,spion1,rRes2,g2[i],pro3_1);
            propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk12,skstr3,spion2,rRes2,g2[i],pro3_2);
            if (g0[i]==1)  {
                propagatorRBW(mass1sq,mass1[i],width1[i],skstr1,skaon,spion1,rRes2,1,pro0_1);
                propagatorRBW(mass1sq,mass1[i],width1[i],skstr2,skaon,spion2,rRes2,1,pro0_2);
            } else if (g0[i]==0)  {
                pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0; pro2_1[0] = 1; pro2_1[1] = 0; pro2_2[0] = 1; pro2_2[1] = 0;
            }
            if (g1[i]==1)  {
                propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk11,skstr1,spi0,rRes2,g2[i],pro1_1);
                propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk12,skstr2,spi0,rRes2,g2[i],pro1_2);
            } else if (g1[i]==0)  {
                pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0; pro3_1[0] = 1; pro3_1[1] = 0; pro3_2[0] = 1; pro3_2[1] = 0;
 
            }
            Com_Multi(pro0_1,pro1_1,pro_1);
            Com_Multi(pro2_1,pro3_1,pro4_1);
            Com_Multi(pro0_2,pro1_2,pro_2);
            Com_Multi(pro2_2,pro3_2,pro4_2);
            amp_tmp1[0] = tmp1*pro_1[0]-tmp3*pro4_1[0];
            amp_tmp1[1] = tmp1*pro_1[1]-tmp3*pro4_1[1];
            amp_tmp2[0] = tmp2*pro_2[0]-tmp4*pro4_2[0];
            amp_tmp2[1] = tmp2*pro_2[1]-tmp4*pro4_2[1];
            }
            else if (modetype[i]==8) {
                        B1_rho1     = Barrier(1,srho1,spion1,spi0,rRes2,mass1[i]);
                        B1_rho2     = Barrier(1,srho2,spion2,spi0,rRes2,mass1[i]);
                        B1_D_K1_Pi2 = Barrier(1,sD,sk11,spion2,rD2,mD);
                        B1_D_K1_Pi1 = Barrier(1,sD,sk12,spion1,rD2,mD);
                        calt1(pK11,Pip2,t1D_1);
                        calt1(pK12,Pip1,t1D_2);
                        if(g2[i]==0) {
                                for(int w=0; w<4; w++) {
                                        for(int j=0; j<4; j++) {
                                                temp_PDF1 += t1D_1[w]*(G[w][j]-pK11[w]*pK11[j]/sk11)*t1rho1[j]*G[w][w]*G[j][j];
                                                temp_PDF2 += t1D_2[w]*(G[w][j]-pK12[w]*pK12[j]/sk12)*t1rho2[j]*G[w][w]*G[j][j];
                                        }
                                }
                                tmp1 = B1_rho1*B1_D_K1_Pi2*temp_PDF1;
                                tmp2 = B1_rho2*B1_D_K1_Pi1*temp_PDF2;
                        } else if(g2[i]==2) {
                                for(int w=0; w<4; w++) {
                                        for(int j=0; j<4; j++) {
                                                temp_PDF1 += t1D_1[w]*t2k21[w][j]*t1rho1[j]*G[w][w]*G[j][j];
                                                temp_PDF2 += t1D_2[w]*t2k22[w][j]*t1rho2[j]*G[w][w]*G[j][j];
                                        }
                                }
                                if (B2_K1_rho1<0.0) B2_K1_rho1 = Barrier(2,sk11,srho1,skaon,rRes2,mass2[i]);
                                if (B2_K1_rho2<0.0) B2_K1_rho2 = Barrier(2,sk12,srho2,skaon,rRes2,mass2[i]);
                                tmp1 = B1_rho1*B2_K1_rho1*B1_D_K1_Pi2*temp_PDF1;
                                tmp2 = B1_rho2*B2_K1_rho2*B1_D_K1_Pi1*temp_PDF2;
                        }
                        if(g0[i]==1) {
                                propagatorGS(mass1sq,mass1[i],width1[i],srho1,spion1,spi0,rRes2,pro0_1);
                                propagatorGS(mass1sq,mass1[i],width1[i],srho2,spion2,spi0,rRes2,pro0_2);
                        } else if(g0[i]==0) {
                                 pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0;
                        }
                        if(g1[i]==1) {
                                propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk11,srho1,skaon,rRes2,g2[i],pro1_1);
                                propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk12,srho2,skaon,rRes2,g2[i],pro1_2);
                        } else if(g1[i]==0) {
                                pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0;
                        }
                        
                        Com_Multi(pro0_1,pro1_1,pro_1);
                        Com_Multi(pro0_2,pro1_2,pro_2);
                        amp_tmp1[0] = tmp1*pro_1[0];
                        amp_tmp1[1] = tmp1*pro_1[1];
                        amp_tmp2[0] = tmp2*pro_2[0];
                        amp_tmp2[1] = tmp2*pro_2[1];
                }
            else if (modetype[i]==11) {
            B1_kstr3    = Barrier(1,skstr3,skaon,spi0,rRes2,mass1[i]);
            B1_kstr1_Pi1 = Barrier(1,skstr1,skaon,spion1,rRes2,mass1[i]);
            B1_kstr2_Pi2 = Barrier(1,skstr2,skaon,spion2,rRes2,mass1[i]);
            B1_D_K1_Pi2 = Barrier(1,sD,sk11,spion2,rD2,mD);
            B1_D_K1_Pi1 = Barrier(1,sD,sk12,spion1,rD2,mD);
                                for(int w=0; w<4; w++) {
                                        for(int j=0; j<4; j++) {
                                                for(int k=0; k<4; k++){
                                                      for(int l=0; l<4; l++){
                                                temp_PDF1 += E[w][j][k][l]*pK11[w]*G[w][w]*(pKstr1[j]-Pi0[j])*G[j][j]*Pip2[k]*G[k][k]*(Km[l]-Pip1[l])*G[l][l];
                                                temp_PDF2 += E[w][j][k][l]*pK12[w]*G[w][w]*(pKstr2[j]-Pi0[j])*G[j][j]*Pip1[k]*G[k][k]*(Km[l]-Pip2[l])*G[l][l];
                                                temp_PDF3 += E[w][j][k][l]*pK11[w]*G[w][w]*(pKstr3[j]-Pip1[j])*G[j][j]*Pip2[k]*G[k][k]*(Km[l]-Pi0[l])*G[l][l];
                                                temp_PDF4 += E[w][j][k][l]*pK12[w]*G[w][w]*(pKstr3[j]-Pip2[j])*G[j][j]*Pip1[k]*G[k][k]*(Km[l]-Pi0[l])*G[l][l];
                                                }
                                        }
                                }
                        }
                                B1_K1_Kstr1_Pi0 = Barrier(1,sk11,skstr1,spi0,rRes2,mass2[i]);
                                B1_K1_Kstr2_Pi0 = Barrier(1,sk12,skstr2,spi0,rRes2,mass2[i]);
                                B1_K1_Kstr3_Pi1 = Barrier(1,sk11,skstr3,spion1,rRes2,mass2[i]);
                                B1_K1_Kstr3_Pi1 = Barrier(1,sk12,skstr3,spion2,rRes2,mass2[i]);
                                tmp1 = B1_kstr1_Pi1*B1_K1_Kstr1_Pi0*B1_D_K1_Pi2*temp_PDF1;
                                tmp2 = B1_kstr2_Pi2*B1_K1_Kstr2_Pi0*B1_D_K1_Pi1*temp_PDF2;
                                tmp3 = B1_kstr3*B1_K1_Kstr3_Pi1*B1_D_K1_Pi2*temp_PDF3;
                                tmp4 = B1_kstr3*B1_K1_Kstr3_Pi2*B1_D_K1_Pi1*temp_PDF4;
            propagatorRBW(mass1sq,mass1[i],width1[i],skstr3,skaon,spi0,rRes2,1,pro2_1);
            propagatorRBW(mass1sq,mass1[i],width1[i],skstr3,skaon,spi0,rRes2,1,pro2_2);
            propagatorRBW(mass2sq,mass2[i],width2[i],sk11,skstr3,spion1,rRes2,1,pro3_1);
            propagatorRBW(mass2sq,mass2[i],width2[i],sk12,skstr3,spion2,rRes2,1,pro3_2);
                        if(g0[i]==1) {
                                propagatorRBW(mass1sq,mass1[i],width1[i],skstr1,skaon,spion1,rRes2,1,pro0_1);
                                propagatorRBW(mass1sq,mass1[i],width1[i],skstr2,skaon,spion2,rRes2,1,pro0_2);
                        } else if(g0[i]==0) {
                pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0; pro2_1[0] = 1; pro2_1[1] = 0; pro2_2[0] = 1; pro2_2[1] = 0;
                        }
                        if(g1[i]==1) {
                                propagatorRBW(mass2sq,mass2[i],width2[i],sk11,skstr1,spi0,rRes2,1,pro1_1);
                                propagatorRBW(mass2sq,mass2[i],width2[i],sk12,skstr2,spi0,rRes2,1,pro1_2);
                        } else if(g1[i]==0) {
                pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0; pro3_1[0] = 1; pro3_1[1] = 0; pro3_2[0] = 1; pro3_2[1] = 0;
                        }
            Com_Multi(pro0_1,pro1_1,pro_1);
            Com_Multi(pro2_1,pro3_1,pro4_1);
            Com_Multi(pro0_2,pro1_2,pro_2);
            Com_Multi(pro2_2,pro3_2,pro4_2);
            amp_tmp1[0] = tmp1*pro_1[0]-tmp3*pro4_1[0];
            amp_tmp1[1] = tmp1*pro_1[1]-tmp3*pro4_1[1];
            amp_tmp2[0] = tmp2*pro_2[0]-tmp4*pro4_2[0];
            amp_tmp2[1] = tmp2*pro_2[1]-tmp4*pro4_2[1];
                }
 
 
        else if (modetype[i]==12) {
            B1_K1_Kstr3_Pi1 = Barrier(1,sk11,skstr3,spion1,rRes2,mass2[i]);
            B1_K1_Kstr3_Pi2 = Barrier(1,sk12,skstr3,spion2,rRes2,mass2[i]);
            B1_K1_Kstr1_Pi0 = Barrier(1,sk11,skstr1,spi0,rRes2,mass2[i]);
            B1_K1_Kstr2_Pi0 = Barrier(1,sk12,skstr2,spi0,rRes2,mass2[i]);
            B1_kstr3    = Barrier(1,skstr3,skaon,spi0,rRes2,mass1[i]);
            B1_kstr1_Pi1= Barrier(1,skstr1,skaon,spion1,rRes2,mass1[i]);
            B1_kstr2_Pi2= Barrier(1,skstr2,skaon,spion2,rRes2,mass1[i]);
                for(int w=0; w<4; w++) {
                                temp_PDF1 += G[w][w]*Pi0[w]*t1kstr1[w];
                                temp_PDF2 += G[w][w]*Pi0[w]*t1kstr2[w];
                                temp_PDF3 += G[w][w]*Pip1[w]*t1kstr3[w];
                                temp_PDF4 += G[w][w]*Pip2[w]*t1kstr3[w];
                }
                tmp1 = B1_kstr1_Pi1*B1_K1_Kstr1_Pi0*temp_PDF1;
                tmp2 = B1_kstr2_Pi2*B1_K1_Kstr2_Pi0*temp_PDF2;
                tmp3 = B1_kstr3*B1_K1_Kstr3_Pi1*temp_PDF3;
                tmp4 = B1_kstr3*B1_K1_Kstr3_Pi2*temp_PDF4;
            propagatorRBW(mass1sq,mass1[i],width1[i],skstr3,skaon,spi0,rRes2,1,pro2_1);
            propagatorRBW(mass1sq,mass1[i],width1[i],skstr3,skaon,spi0,rRes2,1,pro2_2);
            propagatorRBW(mass2sq,mass2[i],width2[i],sk11,skstr3,spion1,rRes2,g2[i],pro3_1);
            propagatorRBW(mass2sq,mass2[i],width2[i],sk12,skstr3,spion2,rRes2,g2[i],pro3_2);
            if (g0[i]==1)  {
                propagatorRBW(mass1sq,mass1[i],width1[i],skstr1,skaon,spion1,rRes2,1,pro0_1);
                propagatorRBW(mass1sq,mass1[i],width1[i],skstr2,skaon,spion2,rRes2,1,pro0_2);
            } else if (g0[i]==0)  {
                pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0; pro2_1[0] = 1; pro2_1[1] = 0; pro2_2[0] = 1; pro2_2[1] = 0;
            }
            if (g1[i]==1)  {
                propagatorRBW(mass2sq,mass2[i],width2[i],sk11,skstr1,spi0,rRes2,g2[i],pro1_1);
                propagatorRBW(mass2sq,mass2[i],width2[i],sk12,skstr2,spi0,rRes2,g2[i],pro1_2);
            } else if (g1[i]==0)  {
                pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0; pro3_1[0] = 1; pro3_1[1] = 0; pro3_2[0] = 1; pro3_2[1] = 0;
 
            }
            Com_Multi(pro0_1,pro1_1,pro_1);
            Com_Multi(pro2_1,pro3_1,pro4_1);
            Com_Multi(pro0_2,pro1_2,pro_2);
            Com_Multi(pro2_2,pro3_2,pro4_2);
            amp_tmp1[0] = tmp1*pro_1[0]-tmp3*pro4_1[0];
            amp_tmp1[1] = tmp1*pro_1[1]-tmp3*pro4_1[1];
            amp_tmp2[0] = tmp2*pro_2[0]-tmp4*pro4_2[0];
            amp_tmp2[1] = tmp2*pro_2[1]-tmp4*pro4_2[1];
 
}
      if(modetype[i]==13)
        {       
                B1_K1_rho1 = Barrier(1,sk11,srho1,skaon,rRes2,mass2[i]);
                B1_K1_rho2 = Barrier(1,sk12,srho2,skaon,rRes2,mass2[i]);
                B1_rho1 = Barrier(1,srho1,spion1,spi0,rRes2,mass1[i]);
                B1_rho2 = Barrier(1,srho2,spion2,spi0,rRes2,mass1[i]);
                        for(int w=0; w<4; w++)
                        {       
                                temp_PDF1 += G[w][w]*Km[w]*t1rho1[w];
                                temp_PDF2 += G[w][w]*Km[w]*t1rho2[w];
                        }
                        tmp1 = B1_rho1*B1_K1_rho1*temp_PDF1;
                        tmp2 = B1_rho2*B1_K1_rho2*temp_PDF2;
                if(g1[i]==1){
                        propagatorRBW(mass2sq,mass2[i],width2[i],sk11,skaon,srho1,rRes2,1,pro0_1);
                        propagatorRBW(mass2sq,mass2[i],width2[i],sk12,skaon,srho2,rRes2,1,pro0_2);
                }
                else if(g1[i]==0){
                        pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0;
                }
                if(g0[i]==1){
                        propagatorGS(mass1sq,mass1[i],width1[i],srho1,spion1,spi0,rRes2,pro1_1);
                        propagatorGS(mass1sq,mass1[i],width1[i],srho2,spion2,spi0,rRes2,pro1_2);
                }
                else if(g0[i]==0){
                        pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0;
                }
                Com_Multi(pro0_1,pro1_1,pro_1);
                Com_Multi(pro0_2,pro1_2,pro_2);
                amp_tmp1[0] = tmp1*pro_1[0];
                amp_tmp1[1] = tmp1*pro_1[1];
                amp_tmp2[0] = tmp2*pro_2[0];
                amp_tmp2[1] = tmp2*pro_2[1];
        }
                else if (modetype[i]==17)
                {
                        B1_rho1         = Barrier(1,srho1,spion1,spi0,rRes2,mass2[i]);
                        B1_rho2         = Barrier(1,srho2,spion2,spi0,rRes2,mass2[i]);
                        B1_D_Kstr2_rho1 = Barrier(1,sD,skstr2,srho1,rD2,mD);
                        B1_D_Kstr1_rho2 = Barrier(1,sD,skstr1,srho2,rD2,mD);
                        calt1(pKstr2,prho1,t1D_1);
                        calt1(pKstr1,prho2,t1D_2);
                        for(int w=0; w<4; w++){
                                temp_PDF1 += t1D_1[w]*t1rho1[w]*G[w][w];
                                temp_PDF2 += t1D_2[w]*t1rho2[w]*G[w][w];
                        }
                        tmp1 = temp_PDF1*B1_rho1*B1_D_Kstr2_rho1;
                        tmp2 = temp_PDF2*B1_rho2*B1_D_Kstr1_rho2;
                        if(g1[i]==1) {
                                propagatorGS(mass2sq,mass2[i],width2[i],srho1,spion1,spi0,rRes2,pro_1);
                                propagatorGS(mass2sq,mass2[i],width2[i],srho2,spion2,spi0,rRes2,pro_2);
                        } else if(g1[i]==0) {
                                pro_1[0] = 1; pro_1[1] = 0; pro_2[0] = 1; pro_2[1] = 0;
                        }
                        temp_1[0] = tmp1*pro_1[0];
                        temp_1[1] = tmp1*pro_1[1];
                        temp_2[0] = tmp2*pro_2[0];
                        temp_2[1] = tmp2*pro_2[1];
                        //printf("xx: %.8f\n", temp_2[0]);
                        //printf("xxxx: %.8f\n", temp_2[1]);
                        KPiSLASS(skstr1,skaon,spion1,Amp_KPiS1);
                        KPiSLASS(skstr2,skaon,spion2,Amp_KPiS2);
                        Com_Multi(temp_1,Amp_KPiS2,amp_tmp1);
                        Com_Multi(temp_2,Amp_KPiS1,amp_tmp2);
                }
            else if (modetype[i]==9) {
                B1_kstr1_Pi1 = Barrier(1,skstr1,skaon,spion1,rRes2,mass1[i]);
                B1_kstr2_Pi2 = Barrier(1,skstr2,skaon,spion2,rRes2,mass1[i]);
                B1_D_K1p_Pi0 = Barrier(1,sD,sk13,spi0,rD2,mD);
                        calt1(pK13,Pi0,t1D_1);
                        if(g2[i]==0) {
                                for(int w=0; w<4; w++) {
                                        for(int j=0; j<4; j++) {
                                                temp_PDF1 += t1D_1[w]*(G[w][j]-pK13[w]*pK13[j]/sk13)*t1kstr1[j]*G[w][w]*G[j][j];
                                                temp_PDF2 += t1D_1[w]*(G[w][j]-pK13[w]*pK13[j]/sk13)*t1kstr2[j]*G[w][w]*G[j][j];
                                        }       
                                }       
                                tmp1 = B1_kstr1_Pi1*B1_D_K1p_Pi0*temp_PDF1;
                                tmp2 = B1_kstr2_Pi2*B1_D_K1p_Pi0*temp_PDF2;
                        } else if(g2[i]==2) {
                                for(int w=0; w<4; w++) {
                                        for(int j=0; j<4; j++) {
                                                temp_PDF1 += t1D_1[w]*t2kk1[w][j]*t1kstr1[j]*G[w][w]*G[j][j];
                                                temp_PDF2 += t1D_1[w]*t2kk2[w][j]*t1kstr2[j]*G[w][w]*G[j][j];
                                        }       
                                }       
                                if (B2_K1p_kstr1<0.0) B2_K1p_kstr1 = Barrier(2,sk13,skstr1,spion2,rRes2,mass2[i]);
                                if (B2_K1p_kstr2<0.0) B2_K1p_kstr2 = Barrier(2,sk13,skstr2,spion1,rRes2,mass2[i]);
                                tmp1 = B1_kstr1_Pi1*B2_K1p_kstr1*B1_D_K1p_Pi0*temp_PDF1;
                                tmp2 = B1_kstr2_Pi2*B2_K1p_kstr2*B1_D_K1p_Pi0*temp_PDF2;
                        }       
                        if(g0[i]==1) {
                                propagatorRBW(mass1sq,mass1[i],width1[i],skstr1,skaon,spion1,rRes2,1,pro0_1);
                                propagatorRBW(mass1sq,mass1[i],width1[i],skstr2,skaon,spion2,rRes2,1,pro0_2);
                        } else if(g0[i]==0) {
                                 pro0_1[0] = 1; pro0_1[1] = 0; pro0_2[0] = 1; pro0_2[1] = 0;
                        }        
                        if(g1[i]==1) {
                                propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk13,skstr1,spion2,rRes2,g2[i],pro1_1);
                                propagatorRBW_K1(mass2sq,mass2[i],width2[i],sk13,skstr2,spion1,rRes2,g2[i],pro1_2);
                        } else if(g1[i]==0) {
                                pro1_1[0] = 1; pro1_1[1] = 0; pro1_2[0] = 1; pro1_2[1] = 0;
                        }       
                        Com_Multi(pro0_1,pro1_1,pro_1);
                        Com_Multi(pro0_2,pro1_2,pro_2);
                        amp_tmp1[0] = tmp1*pro_1[0];
                        amp_tmp1[1] = tmp1*pro_1[1];
                        amp_tmp2[0] = tmp2*pro_2[0];
                        amp_tmp2[1] = tmp2*pro_2[1];
                } 
 
        amp_tmp[0] = amp_tmp1[0]+amp_tmp2[0];
        amp_tmp[1] = amp_tmp1[1]+amp_tmp2[1];
    
        Com_Multi(amp_tmp,cof,amp_PDF);
        PDF[0] += amp_PDF[0];
        PDF[1] += amp_PDF[1];
        double value = PDF[0]*PDF[0] + PDF[1]*PDF[1];
        if(value <=0) {value = 1e-20;}
        Result = value;
 
}
}