BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDToKmPipPipPi0.cc
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2// Environment:
3// This software is part of models developed at BES collaboration
4// based on the EvtGen framework. If you use all or part
5// of it, please give an appropriate acknowledgement.
6//
7// Copyright Information: See EvtGen/BesCopyright
8// Copyright (A) 2006 Ping Rong-Gang @IHEP
9//
10// Module: EvtDToKmPipPipPi0.cc
11// the necessary file: EvtDToKmPipPipPi0.hh
12//
13// Description: D+ -> K- pi+ pi+ pi0 (BAM-00739)
14//
15// Modification history:
16//
17// Liaoyuan Dong Sat Dec 31 13:54:49 2022 Module created
18// Fri Mar 1 23:43:09 2024 Module updated
19//------------------------------------------------------------------------
20#include "EvtDToKmPipPipPi0.hh"
21#include "../EvtGenBase/EvtConst.hh"
22#include "../EvtGenBase/EvtDecayTable.hh"
23#include "../EvtGenBase/EvtFlatte.hh"
24#include "../EvtGenBase/EvtGenKine.hh"
25#include "../EvtGenBase/EvtPDL.hh"
26#include "../EvtGenBase/EvtParticle.hh"
27#include "../EvtGenBase/EvtPatches.hh"
28#include "../EvtGenBase/EvtReport.hh"
29#include "../EvtGenBase/EvtResonance.hh"
30#include "../EvtGenBase/EvtResonance2.hh"
31#include "TComplex.h"
32#include "TMatrix.h"
33#include "TMatrixD.h"
34#include <fstream>
35#include <stdlib.h>
36#include <string>
37// #include "width.h"
38
39using std::endl;
40
41EvtDToKmPipPipPi0::~EvtDToKmPipPipPi0() {}
42
43void EvtDToKmPipPipPi0::getName( std::string& model_name ) { model_name = "DToKmPipPipPi0"; }
44
45EvtDecayBase* EvtDToKmPipPipPi0::clone() { return new EvtDToKmPipPipPi0; }
46
47double a1_width[] = {
48 0.0000000172, 0.0000000487, 0.0000001049, 0.0000001927, 0.0000003185, 0.0000004888,
49 0.0000007099, 0.0000009878, 0.0000013285, 0.0000017378, 0.0000022212, 0.0000027845,
50 0.0000034328, 0.0000041717, 0.0000050061, 0.0000059414, 0.0000069824, 0.0000081342,
51 0.0000094014, 0.0000107890, 0.0000123016, 0.0000139437, 0.0000157200, 0.0000176350,
52 0.0000196931, 0.0000218986, 0.0000242559, 0.0000267692, 0.0000294428, 0.0000322808,
53 0.0000352875, 0.0000384668, 0.0000418228, 0.0000453596, 0.0000490811, 0.0000529913,
54 0.0000570941, 0.0000613934, 0.0000658931, 0.0000705969, 0.0000755087, 0.0000806324,
55 0.0000859716, 0.0000915300, 0.0000973115, 0.0001033197, 0.0001095583, 0.0001160310,
56 0.0001227415, 0.0001296933, 0.0001368901, 0.0001443355, 0.0001520332, 0.0001599867,
57 0.0001681995, 0.0001766753, 0.0001854177, 0.0001944301, 0.0002037162, 0.0002132795,
58 0.0002231235, 0.0002332517, 0.0002436677, 0.0002543750, 0.0002653771, 0.0002766775,
59 0.0002882799, 0.0003001876, 0.0003124041, 0.0003249331, 0.0003377780, 0.0003509423,
60 0.0003644296, 0.0003782433, 0.0003923870, 0.0004068643, 0.0004216785, 0.0004368333,
61 0.0004523323, 0.0004681788, 0.0004843765, 0.0005009290, 0.0005178397, 0.0005351122,
62 0.0005527502, 0.0005707571, 0.0005891365, 0.0006078922, 0.0006270275, 0.0006465462,
63 0.0006664518, 0.0006867480, 0.0007074384, 0.0007285267, 0.0007500165, 0.0007719114,
64 0.0007942151, 0.0008169314, 0.0008400639, 0.0008636163, 0.0008875923, 0.0009119957,
65 0.0009368303, 0.0009620997, 0.0009878078, 0.0010139580, 0.0010405550, 0.0010676020,
66 0.0010951030, 0.0011230610, 0.0011514820, 0.0011803670, 0.0012097230, 0.0012395510,
67 0.0012698570, 0.0013006440, 0.0013319160, 0.0013636770, 0.0013959320, 0.0014286840,
68 0.0014619370, 0.0014956950, 0.0015299630, 0.0015647450, 0.0016000440, 0.0016358650,
69 0.0016722130, 0.0017090900, 0.0017465020, 0.0017844530, 0.0018229470, 0.0018619890,
70 0.0019015820, 0.0019417310, 0.0019824400, 0.0020237150, 0.0020655590, 0.0021079760,
71 0.0021509730, 0.0021945510, 0.0022387180, 0.0022834760, 0.0023288310, 0.0023747880,
72 0.0024213500, 0.0024685230, 0.0025163120, 0.0025647210, 0.0026137560, 0.0026634200,
73 0.0027137200, 0.0027646590, 0.0028162430, 0.0028684780, 0.0029213670, 0.0029749170,
74 0.0030291320, 0.0030840170, 0.0031395780, 0.0031958200, 0.0032527490, 0.0033103700,
75 0.0033686880, 0.0034277080, 0.0034874370, 0.0035478800, 0.0036090420, 0.0036709290,
76 0.0037335480, 0.0037969030, 0.0038610000, 0.0039258460, 0.0039914470, 0.0040578080,
77 0.0041249360, 0.0041928370, 0.0042615160, 0.0043309810, 0.0044012380, 0.0044722930,
78 0.0045441520, 0.0046168220, 0.0046903100, 0.0047646230, 0.0048397670, 0.0049157480,
79 0.0049925750, 0.0050702540, 0.0051487920, 0.0052281970, 0.0053084750, 0.0053896340,
80 0.0054716810, 0.0055546240, 0.0056384710, 0.0057232280, 0.0058089050, 0.0058955090,
81 0.0059830470, 0.0060715290, 0.0061609620, 0.0062513540, 0.0063427140, 0.0064350500,
82 0.0065283720, 0.0066226870, 0.0067180050, 0.0068143350, 0.0069116850, 0.0070100650,
83 0.0071094850, 0.0072099530, 0.0073114800, 0.0074140740, 0.0075177470, 0.0076225070,
84 0.0077283660, 0.0078353330, 0.0079434180, 0.0080526330, 0.0081629880, 0.0082744940,
85 0.0083871620, 0.0085010030, 0.0086160290, 0.0087322510, 0.0088496800, 0.0089683290,
86 0.0090882100, 0.0092093350, 0.0093317160, 0.0094553660, 0.0095802980, 0.0097065250,
87 0.0098340590, 0.0099629140, 0.0100931000, 0.0102246400, 0.0103575500, 0.0104918200,
88 0.0106274900, 0.0107645700, 0.0109030700, 0.0110430000, 0.0111843800, 0.0113272300,
89 0.0114715700, 0.0116174000, 0.0117647500, 0.0119136300, 0.0120640600, 0.0122160600,
90 0.0123696300, 0.0125248100, 0.0126816100, 0.0128400500, 0.0130001400, 0.0131619100,
91 0.0133253800, 0.0134905500, 0.0136574600, 0.0138261300, 0.0139965600, 0.0141688000,
92 0.0143428500, 0.0145187400, 0.0146964900, 0.0148761200, 0.0150576500, 0.0152411200,
93 0.0154265300, 0.0156139200, 0.0158033100, 0.0159947200, 0.0161881800, 0.0163837200,
94 0.0165813500, 0.0167811100, 0.0169830300, 0.0171871200, 0.0173934200, 0.0176019600,
95 0.0178127600, 0.0180258500, 0.0182412700, 0.0184590400, 0.0186791900, 0.0189017500,
96 0.0191267600, 0.0193542500, 0.0195842500, 0.0198167900, 0.0200519000, 0.0202896300,
97 0.0205300000, 0.0207730500, 0.0210188200, 0.0212673300, 0.0215186400, 0.0217727700,
98 0.0220297700, 0.0222896700, 0.0225525200, 0.0228183400, 0.0230871900, 0.0233591000,
99 0.0236341200, 0.0239122800, 0.0241936400, 0.0244782300, 0.0247661000, 0.0250572900,
100 0.0253518600, 0.0256498400, 0.0259512800, 0.0262562300, 0.0265647400, 0.0268768600,
101 0.0271926300, 0.0275121100, 0.0278353500, 0.0281624000, 0.0284933100, 0.0288281300,
102 0.0291669200, 0.0295097300, 0.0298566200, 0.0302076400, 0.0305628500, 0.0309223000,
103 0.0312860500, 0.0316541600, 0.0320266900, 0.0324036900, 0.0327852300, 0.0331713600,
104 0.0335621500, 0.0339576500, 0.0343579300, 0.0347630500, 0.0351730600, 0.0355880400,
105 0.0360080500, 0.0364331400, 0.0368633900, 0.0372988500, 0.0377395900, 0.0381856600,
106 0.0386371500, 0.0390941000, 0.0395565900, 0.0400246700, 0.0404984100, 0.0409778800,
107 0.0414631300, 0.0419542300, 0.0424512400, 0.0429542200, 0.0434632400, 0.0439783500,
108 0.0444996200, 0.0450271000, 0.0455608500, 0.0461009400, 0.0466474000, 0.0472003200,
109 0.0477597200, 0.0483256700, 0.0488982300, 0.0494774300, 0.0500633300, 0.0506559700,
110 0.0512554000, 0.0518616600, 0.0524748000, 0.0530948400, 0.0537218300, 0.0543558000,
111 0.0549967800, 0.0556448000, 0.0562998900, 0.0569620700, 0.0576313600, 0.0583077800,
112 0.0589913500, 0.0596820600, 0.0603799400, 0.0610849900, 0.0617972100, 0.0625166000,
113 0.0632431500, 0.0639768500, 0.0647176900, 0.0654656600, 0.0662207300, 0.0669828800,
114 0.0677520900, 0.0685283100, 0.0693115200, 0.0701016800, 0.0708987400, 0.0717026500,
115 0.0725133700, 0.0733308400, 0.0741550000, 0.0749858000, 0.0758231500, 0.0766670000,
116 0.0775172700, 0.0783738800, 0.0792367600, 0.0801058200, 0.0809809800, 0.0818621400,
117 0.0827492100, 0.0836421000, 0.0845407100, 0.0854449400, 0.0863546800, 0.0872698400,
118 0.0881902900, 0.0891159400, 0.0900466700, 0.0909823700, 0.0919229100, 0.0928682000,
119 0.0938181000, 0.0947725000, 0.0957312800, 0.0966943100, 0.0976614900, 0.0986326700,
120 0.0996077500, 0.1005866000, 0.1015691000, 0.1025551000, 0.1035445000, 0.1045372000,
121 0.1055331000, 0.1065320000, 0.1075338000, 0.1085384000, 0.1095456000, 0.1105555000,
122 0.1115677000, 0.1125823000, 0.1135991000, 0.1146180000, 0.1156389000, 0.1166616000,
123 0.1176862000, 0.1187123000, 0.1197401000, 0.1207692000, 0.1217997000, 0.1228315000,
124 0.1238644000, 0.1248983000, 0.1259332000, 0.1269690000, 0.1280055000, 0.1290427000,
125 0.1300805000, 0.1311188000, 0.1321575000, 0.1331966000, 0.1342359000, 0.1352754000,
126 0.1363150000, 0.1373546000, 0.1383942000, 0.1394337000, 0.1404730000, 0.1415120000,
127 0.1425508000, 0.1435891000, 0.1446271000, 0.1456645000, 0.1467014000, 0.1477377000,
128 0.1487733000, 0.1498082000, 0.1508423000, 0.1518756000, 0.1529081000, 0.1539396000,
129 0.1549702000, 0.1559998000, 0.1570284000, 0.1580558000, 0.1590822000, 0.1601074000,
130 0.1611314000, 0.1621541000, 0.1631757000, 0.1641959000, 0.1652148000, 0.1662323000,
131 0.1672484000, 0.1682632000, 0.1692765000, 0.1702883000, 0.1712986000, 0.1723075000,
132 0.1733147000, 0.1743205000, 0.1753246000, 0.1763271000, 0.1773281000, 0.1783273000,
133 0.1793250000, 0.1803209000, 0.1813152000, 0.1823078000, 0.1832986000, 0.1842877000,
134 0.1852751000, 0.1862607000, 0.1872445000, 0.1882266000, 0.1892068000, 0.1901853000,
135 0.1911619000, 0.1921368000, 0.1931098000, 0.1940809000, 0.1950503000, 0.1960177000,
136 0.1969834000, 0.1979471000, 0.1989090000, 0.1998690000, 0.2008271000, 0.2017834000,
137 0.2027378000, 0.2036903000, 0.2046408000, 0.2055895000, 0.2065363000, 0.2074812000,
138 0.2084242000, 0.2093653000, 0.2103045000, 0.2112418000, 0.2121772000, 0.2131106000,
139 0.2140422000, 0.2149719000, 0.2158996000, 0.2168255000, 0.2177494000, 0.2186715000,
140 0.2195916000, 0.2205098000, 0.2214262000, 0.2223406000, 0.2232532000, 0.2241638000,
141 0.2250726000, 0.2259795000, 0.2268844000, 0.2277875000, 0.2286888000, 0.2295881000,
142 0.2304856000, 0.2313812000, 0.2322749000, 0.2331668000, 0.2340568000, 0.2349449000,
143 0.2358312000, 0.2367157000, 0.2375983000, 0.2384790000, 0.2393580000, 0.2402351000,
144 0.2411103000, 0.2419837000, 0.2428554000, 0.2437252000, 0.2445931000, 0.2454593000,
145 0.2463237000, 0.2471863000, 0.2480471000, 0.2489061000, 0.2497633000, 0.2506187000,
146 0.2514724000, 0.2523243000, 0.2531744000, 0.2540228000, 0.2548694000, 0.2557143000,
147 0.2565574000, 0.2573988000, 0.2582385000, 0.2590764000, 0.2599126000, 0.2607471000,
148 0.2615799000, 0.2624110000, 0.2632404000, 0.2640681000, 0.2648941000, 0.2657184000,
149 0.2665410000, 0.2673620000, 0.2681813000, 0.2689989000, 0.2698149000, 0.2706292000,
150 0.2714419000, 0.2722529000, 0.2730623000, 0.2738701000, 0.2746762000, 0.2754807000,
151 0.2762836000, 0.2770849000, 0.2778846000, 0.2786827000, 0.2794792000, 0.2802741000,
152 0.2810675000, 0.2818592000, 0.2826494000, 0.2834380000, 0.2842251000, 0.2850106000,
153 0.2857946000, 0.2865770000, 0.2873578000, 0.2881372000, 0.2889150000, 0.2896913000,
154 0.2904661000, 0.2912393000, 0.2920111000, 0.2927813000, 0.2935501000, 0.2943173000,
155 0.2950831000, 0.2958474000, 0.2966102000, 0.2973715000, 0.2981314000, 0.2988898000,
156 0.2996468000, 0.3004023000, 0.3011564000, 0.3019090000, 0.3026602000, 0.3034099000,
157 0.3041583000, 0.3049052000, 0.3056507000, 0.3063948000, 0.3071374000, 0.3078787000,
158 0.3086186000, 0.3093571000, 0.3100942000, 0.3108299000, 0.3115643000, 0.3122973000,
159 0.3130289000, 0.3137591000, 0.3144880000, 0.3152155000, 0.3159417000, 0.3166666000,
160 0.3173901000, 0.3181122000, 0.3188331000, 0.3195526000, 0.3202708000, 0.3209877000,
161 0.3217033000, 0.3224175000, 0.3231305000, 0.3238421000, 0.3245525000, 0.3252616000,
162 0.3259694000, 0.3266759000, 0.3273811000, 0.3280851000, 0.3287878000, 0.3294893000,
163 0.3301894000, 0.3308884000, 0.3315861000, 0.3322825000, 0.3329777000, 0.3336716000,
164 0.3343644000, 0.3350559000, 0.3357461000, 0.3364352000, 0.3371230000, 0.3378096000,
165 0.3384951000, 0.3391793000, 0.3398623000, 0.3405441000, 0.3412247000, 0.3419042000,
166 0.3425824000, 0.3432595000, 0.3439354000, 0.3446101000, 0.3452837000, 0.3459561000,
167 0.3466273000, 0.3472974000, 0.3479663000, 0.3486341000, 0.3493008000, 0.3499663000,
168 0.3506306000, 0.3512939000, 0.3519560000, 0.3526169000, 0.3532768000, 0.3539355000,
169 0.3545931000, 0.3552497000, 0.3559051000, 0.3565594000, 0.3572126000, 0.3578647000,
170 0.3585157000, 0.3591656000, 0.3598144000, 0.3604622000, 0.3611089000, 0.3617545000,
171 0.3623990000, 0.3630425000, 0.3636849000, 0.3643263000, 0.3649665000, 0.3656058000,
172 0.3662440000, 0.3668811000, 0.3675172000, 0.3681523000, 0.3687863000, 0.3694193000,
173 0.3700512000, 0.3706822000, 0.3713121000, 0.3719409000, 0.3725688000, 0.3731957000,
174 0.3738215000, 0.3744464000, 0.3750702000, 0.3756930000, 0.3763149000, 0.3769357000,
175 0.3775556000, 0.3781745000, 0.3787924000, 0.3794093000, 0.3800252000, 0.3806401000,
176 0.3812541000, 0.3818671000, 0.3824792000, 0.3830902000, 0.3837004000, 0.3843095000,
177 0.3849178000, 0.3855250000, 0.3861313000, 0.3867367000, 0.3873411000, 0.3879446000,
178 0.3885472000, 0.3891488000, 0.3897495000, 0.3903492000, 0.3909481000, 0.3915460000,
179 0.3921430000, 0.3927391000, 0.3933343000, 0.3939285000, 0.3945219000, 0.3951143000,
180 0.3957059000, 0.3962966000, 0.3968863000, 0.3974752000, 0.3980631000, 0.3986502000,
181 0.3992364000, 0.3998218000, 0.4004062000, 0.4009898000, 0.4015725000, 0.4021543000,
182 0.4027352000, 0.4033153000, 0.4038946000, 0.4044729000, 0.4050504000, 0.4056271000,
183 0.4062029000, 0.4067778000, 0.4073519000, 0.4079252000, 0.4084976000, 0.4090692000,
184 0.4096399000, 0.4102098000, 0.4107789000, 0.4113472000, 0.4119146000, 0.4124812000,
185 0.4130470000, 0.4136119000, 0.4141760000, 0.4147394000, 0.4153019000, 0.4158636000,
186 0.4164245000, 0.4169846000, 0.4175439000, 0.4181024000, 0.4186601000, 0.4192170000,
187 0.4197731000, 0.4203284000, 0.4208830000, 0.4214367000, 0.4219897000, 0.4225419000,
188 0.4230933000, 0.4236439000, 0.4241938000, 0.4247429000, 0.4252912000, 0.4258388000,
189 0.4263856000, 0.4269317000, 0.4274769000, 0.4280215000, 0.4285652000, 0.4291083000,
190 0.4296505000, 0.4301921000, 0.4307329000, 0.4312729000, 0.4318122000, 0.4323508000,
191 0.4328886000, 0.4334257000, 0.4339621000, 0.4344977000, 0.4350326000, 0.4355668000,
192 0.4361003000, 0.4366331000, 0.4371651000, 0.4376964000, 0.4382270000, 0.4387569000,
193 0.4392861000, 0.4398146000, 0.4403424000, 0.4408695000, 0.4413959000, 0.4419216000,
194 0.4424466000, 0.4429709000, 0.4434945000, 0.4440174000, 0.4445396000, 0.4450612000,
195 0.4455821000, 0.4461023000, 0.4466218000, 0.4471407000, 0.4476589000, 0.4481764000,
196 0.4486932000, 0.4492094000, 0.4497249000, 0.4502398000, 0.4507540000, 0.4512675000,
197 0.4517804000, 0.4522926000, 0.4528042000, 0.4533152000, 0.4538255000, 0.4543351000,
198 0.4548442000, 0.4553525000, 0.4558603000, 0.4563674000, 0.4568739000, 0.4573797000,
199 0.4578849000, 0.4583895000, 0.4588935000, 0.4593969000, 0.4598996000, 0.4604018000,
200 0.4609033000, 0.4614042000, 0.4619045000, 0.4624042000, 0.4629032000, 0.4634017000,
201 0.4638996000, 0.4643969000, 0.4648936000, 0.4653897000, 0.4658852000, 0.4663802000,
202 0.4668745000, 0.4673683000, 0.4678614000, 0.4683541000, 0.4688461000, 0.4693375000,
203 0.4698284000, 0.4703188000, 0.4708085000, 0.4712977000, 0.4717863000, 0.4722744000,
204 0.4727619000, 0.4732489000, 0.4737353000, 0.4742212000, 0.4747065000, 0.4751913000,
205 0.4756756000, 0.4761593000, 0.4766424000, 0.4771251000, 0.4776072000, 0.4780888000,
206 0.4785699000, 0.4790504000, 0.4795304000, 0.4800099000, 0.4804889000, 0.4809674000,
207 0.4814454000, 0.4819229000, 0.4823999000, 0.4828764000, 0.4833524000, 0.4838279000,
208 0.4843029000, 0.4847774000, 0.4852515000, 0.4857250000, 0.4861981000, 0.4866708000,
209 0.4871429000, 0.4876146000, 0.4880858000, 0.4885566000, 0.4890269000, 0.4894968000,
210 0.4899662000, 0.4904351000, 0.4909037000, 0.4913718000, 0.4918394000, 0.4923067000,
211 0.4927735000, 0.4932398000, 0.4937058000, 0.4941713000, 0.4946365000, 0.4951012000,
212 0.4955655000, 0.4960295000, 0.4964930000, 0.4969562000, 0.4974189000, 0.4978813000,
213 0.4983434000, 0.4988050000, 0.4992663000, 0.4997272000, 0.5001878000, 0.5006480000,
214 0.5011079000, 0.5015674000, 0.5020266000, 0.5024855000, 0.5029440000, 0.5034023000,
215 0.5038602000, 0.5043178000, 0.5047751000, 0.5052321000, 0.5056889000, 0.5061453000,
216 0.5066015000, 0.5070574000, 0.5075131000, 0.5079685000, 0.5084237000, 0.5088786000,
217 0.5093333000, 0.5097878000, 0.5102420000, 0.5106961000, 0.5111499000, 0.5116036000,
218 0.5120571000, 0.5125104000, 0.5129636000, 0.5134166000, 0.5138694000, 0.5143221000,
219 0.5147747000, 0.5152272000, 0.5156796000, 0.5161319000, 0.5165842000, 0.5170363000,
220 0.5174884000, 0.5179405000, 0.5183925000, 0.5188445000, 0.5192965000, 0.5197485000,
221 0.5202005000, 0.5206526000, 0.5211048000, 0.5215570000, 0.5220092000, 0.5224616000,
222 0.5229141000, 0.5233668000, 0.5238195000, 0.5242725000, 0.5247257000, 0.5251790000,
223 0.5256326000, 0.5260864000, 0.5265405000, 0.5269949000, 0.5274496000, 0.5279047000,
224 0.5283601000, 0.5288159000, 0.5292721000, 0.5297287000, 0.5301858000, 0.5306434000,
225 0.5311015000, 0.5315601000, 0.5320194000, 0.5324792000, 0.5329396000, 0.5334008000,
226 0.5338626000, 0.5343252000, 0.5347885000, 0.5352526000, 0.5357176000, 0.5361834000,
227 0.5366502000, 0.5371179000, 0.5375866000, 0.5380563000, 0.5385271000, 0.5389989000,
228 0.5394720000, 0.5399462000, 0.5404217000, 0.5408984000, 0.5413764000, 0.5418558000,
229 0.5423367000, 0.5428189000, 0.5433026000, 0.5437879000, 0.5442747000, 0.5447631000,
230 0.5452532000, 0.5457450000, 0.5462385000, 0.5467337000, 0.5472307000, 0.5477296000,
231 0.5482303000, 0.5487328000, 0.5492373000, 0.5497437000, 0.5502520000, 0.5507623000,
232 0.5512745000, 0.5517887000, 0.5523048000, 0.5528230000, 0.5533430000, 0.5538651000,
233 0.5543890000, 0.5549149000, 0.5554426000, 0.5559722000, 0.5565036000, 0.5570367000,
234 0.5575716000, 0.5581082000, 0.5586464000, 0.5591862000, 0.5597274000, 0.5602702000,
235 0.5608143000, 0.5613597000, 0.5619063000, 0.5624541000, 0.5630030000, 0.5635529000,
236 0.5641037000, 0.5646553000, 0.5652078000, 0.5657609000, 0.5663146000, 0.5668689000,
237 0.5674236000, 0.5679787000, 0.5685341000, 0.5690898000, 0.5696456000, 0.5702015000,
238 0.5707575000, 0.5713134000, 0.5718692000, 0.5724249000, 0.5729803000, 0.5735355000,
239 0.5740903000, 0.5746448000, 0.5751989000, 0.5757525000, 0.5763056000, 0.5768581000,
240 0.5774101000, 0.5779615000, 0.5785121000, 0.5790622000, 0.5796114000, 0.5801600000,
241 0.5807078000, 0.5812548000, 0.5818009000, 0.5823462000, 0.5828907000, 0.5834343000,
242 0.5839770000, 0.5845188000, 0.5850597000, 0.5855996000, 0.5861386000, 0.5866767000,
243 0.5872137000, 0.5877498000, 0.5882850000, 0.5888191000, 0.5893522000, 0.5898844000,
244 0.5904155000, 0.5909457000, 0.5914748000, 0.5920029000, 0.5925300000, 0.5930560000,
245 0.5935811000, 0.5941051000, 0.5946282000, 0.5951502000, 0.5956711000, 0.5961911000,
246 0.5967100000, 0.5972280000, 0.5977449000, 0.5982608000, 0.5987757000, 0.5992896000,
247 0.5998025000, 0.6003144000, 0.6008253000, 0.6013352000, 0.6018441000, 0.6023520000,
248 0.6028589000, 0.6033649000, 0.6038699000, 0.6043739000, 0.6048769000, 0.6053790000,
249 0.6058801000, 0.6063803000, 0.6068795000, 0.6073777000, 0.6078750000, 0.6083714000,
250 0.6088668000, 0.6093613000, 0.6098549000, 0.6103476000, 0.6108393000, 0.6113301000,
251 0.6118201000, 0.6123091000, 0.6127972000, 0.6132844000, 0.6137707000, 0.6142562000,
252 0.6147407000, 0.6152244000, 0.6157072000, 0.6161892000, 0.6166703000, 0.6171505000,
253 0.6176298000, 0.6181084000, 0.6185860000, 0.6190629000, 0.6195388000, 0.6200140000,
254 0.6204883000, 0.6209618000, 0.6214345000, 0.6219064000, 0.6223774000, 0.6228476000,
255 0.6233171000, 0.6237857000, 0.6242536000, 0.6247206000, 0.6251869000, 0.6256523000,
256 0.6261170000, 0.6265809000, 0.6270441000, 0.6275064000, 0.6279681000, 0.6284289000,
257 0.6288890000, 0.6293483000, 0.6298069000, 0.6302647000, 0.6307218000, 0.6311782000,
258 0.6316338000, 0.6320887000, 0.6325428000, 0.6329963000, 0.6334490000, 0.6339009000,
259 0.6343522000, 0.6348028000, 0.6352526000, 0.6357018000, 0.6361502000, 0.6365979000,
260 0.6370450000, 0.6374913000, 0.6379370000, 0.6383820000, 0.6388262000, 0.6392699000,
261 0.6397128000, 0.6401550000, 0.6405966000, 0.6410375000, 0.6414778000, 0.6419174000,
262 0.6423563000, 0.6427946000, 0.6432322000, 0.6436692000, 0.6441055000, 0.6445412000,
263 0.6449762000, 0.6454106000, 0.6458443000, 0.6462774000, 0.6467099000, 0.6471418000,
264 0.6475730000, 0.6480036000, 0.6484336000, 0.6488629000, 0.6492917000, 0.6497198000,
265 0.6501473000, 0.6505742000, 0.6510005000, 0.6514262000, 0.6518513000, 0.6522758000,
266 0.6526997000, 0.6531230000, 0.6535457000, 0.6539678000, 0.6543893000, 0.6548103000,
267 0.6552306000, 0.6556504000, 0.6560696000, 0.6564882000, 0.6569062000, 0.6573237000,
268 0.6577406000, 0.6581569000, 0.6585727000, 0.6589879000, 0.6594025000, 0.6598166000,
269 0.6602301000, 0.6606431000, 0.6610555000, 0.6614673000, 0.6618787000, 0.6622894000,
270 0.6626996000, 0.6631093000, 0.6635184000, 0.6639270000, 0.6643350000, 0.6647425000,
271 0.6651495000, 0.6655560000, 0.6659619000, 0.6663672000, 0.6667721000, 0.6671764000,
272 0.6675802000, 0.6679835000, 0.6683863000, 0.6687885000, 0.6691902000, 0.6695914000,
273 0.6699921000, 0.6703923000, 0.6707920000, 0.6711912000, 0.6715898000, 0.6719880000,
274 0.6723856000, 0.6727828000, 0.6731794000, 0.6735756000, 0.6739712000, 0.6743664000,
275 0.6747610000, 0.6751552000, 0.6755489000, 0.6759421000, 0.6763348000, 0.6767270000,
276 0.6771187000, 0.6775100000, 0.6779008000, 0.6782910000, 0.6786809000, 0.6790702000,
277 0.6794590000, 0.6798474000, 0.6802353000, 0.6806228000, 0.6810097000, 0.6813962000,
278 0.6817823000, 0.6821678000, 0.6825530000, 0.6829376000, 0.6833218000, 0.6837055000,
279 0.6840888000, 0.6844716000, 0.6848539000, 0.6852358000, 0.6856172000, 0.6859982000,
280 0.6863788000, 0.6867588000, 0.6871385000, 0.6875177000, 0.6878964000, 0.6882747000,
281 0.6886526000, 0.6890300000, 0.6894070000, 0.6897835000, 0.6901596000, 0.6905353000,
282 0.6909105000, 0.6912853000, 0.6916596000, 0.6920336000, 0.6924071000, 0.6927801000,
283 0.6931528000, 0.6935250000, 0.6938967000, 0.6942681000, 0.6946390000, 0.6950095000,
284 0.6953796000, 0.6957493000, 0.6961185000, 0.6964874000, 0.6968558000, 0.6972238000,
285 0.6975914000, 0.6979585000, 0.6983253000, 0.6986916000, 0.6990576000, 0.6994231000,
286 0.6997882000, 0.7001529000, 0.7005172000, 0.7008811000, 0.7012446000, 0.7016077000,
287 0.7019704000, 0.7023326000, 0.7026945000, 0.7030560000, 0.7034171000, 0.7037778000,
288 0.7041381000, 0.7044980000, 0.7048575000, 0.7052166000, 0.7055753000, 0.7059336000,
289 0.7062916000, 0.7066491000, 0.7070063000, 0.7073630000, 0.7077194000, 0.7080754000,
290 0.7084310000, 0.7087863000, 0.7091411000, 0.7094956000, 0.7098497000, 0.7102034000,
291 0.7105567000, 0.7109096000, 0.7112622000, 0.7116144000, 0.7119662000, 0.7123177000,
292 0.7126687000, 0.7130194000, 0.7133697000, 0.7137197000, 0.7140693000, 0.7144185000,
293 0.7147673000, 0.7151158000, 0.7154639000, 0.7158117000, 0.7161591000, 0.7165061000,
294 0.7168527000, 0.7171990000, 0.7175450000, 0.7178905000, 0.7182357000, 0.7185806000,
295 0.7189251000, 0.7192692000, 0.7196130000, 0.7199564000, 0.7202995000, 0.7206422000,
296 0.7209845000, 0.7213266000, 0.7216682000, 0.7220095000, 0.7223505000, 0.7226911000,
297 0.7230313000, 0.7233712000, 0.7237108000, 0.7240500000, 0.7243889000, 0.7247274000,
298 0.7250656000, 0.7254034000, 0.7257409000, 0.7260780000, 0.7264148000, 0.7267513000,
299 0.7270874000, 0.7274232000, 0.7277587000, 0.7280938000, 0.7284286000, 0.7287630000,
300 0.7290971000, 0.7294309000, 0.7297643000, 0.7300974000, 0.7304302000, 0.7307626000,
301 0.7310947000, 0.7314265000, 0.7317579000, 0.7320891000, 0.7324198000, 0.7327503000,
302 0.7330804000, 0.7334102000, 0.7337397000, 0.7340689000, 0.7343977000, 0.7347262000,
303 0.7350544000, 0.7353823000, 0.7357098000, 0.7360370000, 0.7363639000, 0.7366905000,
304 0.7370168000, 0.7373427000, 0.7376683000, 0.7379936000, 0.7383186000, 0.7386433000,
305 0.7389677000, 0.7392917000, 0.7396155000, 0.7399389000, 0.7402620000, 0.7405848000,
306 0.7409073000, 0.7412294000, 0.7415513000, 0.7418728000, 0.7421941000, 0.7425150000,
307 0.7428357000, 0.7431560000, 0.7434760000, 0.7437957000, 0.7441151000, 0.7444342000,
308 0.7447530000, 0.7450715000, 0.7453897000, 0.7457076000, 0.7460251000, 0.7463424000,
309 0.7466594000, 0.7469761000, 0.7472925000, 0.7476085000, 0.7479243000, 0.7482398000,
310 0.7485550000, 0.7488699000, 0.7491845000, 0.7494988000, 0.7498128000, 0.7501265000,
311 0.7504399000, 0.7507531000, 0.7510659000, 0.7513784000, 0.7516907000, 0.7520027000,
312 0.7523143000, 0.7526257000, 0.7529368000, 0.7532476000, 0.7535581000, 0.7538683000,
313 0.7541783000, 0.7544879000, 0.7547973000, 0.7551064000, 0.7554152000, 0.7557237000,
314 0.7560319000, 0.7563398000, 0.7566475000, 0.7569549000, 0.7572620000, 0.7575688000,
315 0.7578753000, 0.7581816000, 0.7584875000, 0.7587932000, 0.7590986000, 0.7594038000,
316 0.7597086000, 0.7600132000, 0.7603175000, 0.7606215000, 0.7609253000, 0.7612288000,
317 0.7615320000, 0.7618349000, 0.7621376000, 0.7624399000, 0.7627420000, 0.7630439000,
318 0.7633454000, 0.7636467000, 0.7639477000, 0.7642485000, 0.7645490000, 0.7648492000,
319 0.7651491000, 0.7654488000, 0.7657482000, 0.7660473000, 0.7663462000, 0.7666448000,
320 0.7669431000, 0.7672412000, 0.7675390000, 0.7678366000, 0.7681338000, 0.7684309000,
321 0.7687276000, 0.7690241000, 0.7693203000, 0.7696163000, 0.7699120000, 0.7702074000,
322 0.7705026000, 0.7707975000, 0.7710922000, 0.7713866000, 0.7716807000, 0.7719746000,
323 0.7722682000, 0.7725616000, 0.7728547000, 0.7731475000, 0.7734401000, 0.7737325000,
324 0.7740245000, 0.7743164000, 0.7746079000, 0.7748993000, 0.7751903000, 0.7754811000,
325 0.7757717000, 0.7760620000, 0.7763521000, 0.7766419000, 0.7769314000, 0.7772207000,
326 0.7775098000, 0.7777986000, 0.7780871000, 0.7783754000, 0.7786635000, 0.7789513000,
327 0.7792389000, 0.7795262000, 0.7798132000, 0.7801000000, 0.7803866000, 0.7806729000,
328 0.7809590000, 0.7812449000, 0.7815304000, 0.7818158000, 0.7821009000, 0.7823858000,
329 0.7826704000, 0.7829547000, 0.7832389000, 0.7835228000, 0.7838064000, 0.7840898000,
330 0.7843730000, 0.7846559000, 0.7849386000, 0.7852210000, 0.7855032000, 0.7857852000,
331 0.7860669000, 0.7863484000, 0.7866297000, 0.7869107000, 0.7871915000, 0.7874720000,
332 0.7877523000, 0.7880324000, 0.7883122000, 0.7885919000, 0.7888712000, 0.7891504000,
333 0.7894293000, 0.7897079000, 0.7899864000, 0.7902646000, 0.7905425000, 0.7908203000,
334 0.7910978000, 0.7913750000, 0.7916521000, 0.7919289000, 0.7922055000, 0.7924818000,
335 0.7927579000, 0.7930338000, 0.7933095000, 0.7935849000, 0.7938601000, 0.7941351000,
336 0.7944099000, 0.7946844000, 0.7949587000, 0.7952328000, 0.7955066000, 0.7957803000,
337 0.7960537000, 0.7963268000, 0.7965998000, 0.7968725000, 0.7971450000, 0.7974173000,
338 0.7976893000, 0.7979612000, 0.7982328000, 0.7985041000, 0.7987753000, 0.7990463000,
339 0.7993170000, 0.7995875000, 0.7998578000, 0.8001278000, 0.8003977000, 0.8006673000,
340 0.8009367000, 0.8012059000, 0.8014748000, 0.8017436000, 0.8020121000, 0.8022804000,
341 0.8025485000, 0.8028164000, 0.8030841000, 0.8033515000, 0.8036187000, 0.8038857000,
342 0.8041525000, 0.8044191000, 0.8046855000, 0.8049517000, 0.8052176000, 0.8054833000,
343 0.8057488000, 0.8060141000, 0.8062792000, 0.8065441000, 0.8068088000, 0.8070732000,
344 0.8073375000, 0.8076015000, 0.8078653000, 0.8081289000, 0.8083923000, 0.8086555000,
345 0.8089185000, 0.8091813000, 0.8094438000, 0.8097062000, 0.8099683000, 0.8102303000,
346 0.8104920000, 0.8107535000, 0.8110148000, 0.8112760000, 0.8115369000, 0.8117976000,
347 0.8120580000, 0.8123183000, 0.8125784000, 0.8128383000, 0.8130980000, 0.8133574000,
348 0.8136167000, 0.8138757000, 0.8141346000, 0.8143933000, 0.8146517000, 0.8149099000,
349 0.8151680000, 0.8154258000, 0.8156835000, 0.8159409000, 0.8161982000, 0.8164552000,
350 0.8167120000, 0.8169687000, 0.8172251000, 0.8174813000, 0.8177374000, 0.8179932000,
351 0.8182488000, 0.8185043000, 0.8187595000, 0.8190146000, 0.8192694000, 0.8195241000,
352 0.8197785000, 0.8200328000, 0.8202868000, 0.8205407000, 0.8207944000, 0.8210478000,
353 0.8213011000, 0.8215542000, 0.8218071000, 0.8220598000, 0.8223122000, 0.8225645000,
354 0.8228167000, 0.8230686000, 0.8233203000, 0.8235718000, 0.8238232000, 0.8240743000,
355 0.8243252000, 0.8245760000, 0.8248266000, 0.8250769000, 0.8253271000, 0.8255771000,
356 0.8258269000, 0.8260765000, 0.8263259000, 0.8265752000, 0.8268242000, 0.8270731000,
357 0.8273217000, 0.8275702000, 0.8278185000, 0.8280666000, 0.8283145000, 0.8285622000,
358 0.8288097000, 0.8290571000, 0.8293043000, 0.8295512000, 0.8297980000, 0.8300446000,
359 0.8302910000, 0.8305373000, 0.8307833000, 0.8310292000, 0.8312748000, 0.8315203000,
360 0.8317656000, 0.8320107000, 0.8322557000, 0.8325004000, 0.8327450000, 0.8329894000,
361 0.8332336000, 0.8334776000, 0.8337214000, 0.8339651000, 0.8342085000, 0.8344518000,
362 0.8346949000 };
363
364double K1270_width[] = {
365 0.0000000014, 0.0000000052, 0.0000000127, 0.0000000252, 0.0000000442, 0.0000000709,
366 0.0000001067, 0.0000001528, 0.0000002108, 0.0000002819, 0.0000003675, 0.0000004689,
367 0.0000005876, 0.0000007249, 0.0000008822, 0.0000010609, 0.0000012624, 0.0000014882,
368 0.0000017396, 0.0000020181, 0.0000023251, 0.0000026620, 0.0000030303, 0.0000034314,
369 0.0000038668, 0.0000043379, 0.0000048463, 0.0000053933, 0.0000059805, 0.0000066094,
370 0.0000072815, 0.0000079983, 0.0000087612, 0.0000095719, 0.0000104318, 0.0000113424,
371 0.0000123054, 0.0000133223, 0.0000143946, 0.0000155239, 0.0000167117, 0.0000179597,
372 0.0000192695, 0.0000206426, 0.0000220807, 0.0000235853, 0.0000251581, 0.0000268008,
373 0.0000285149, 0.0000303022, 0.0000321643, 0.0000341028, 0.0000361195, 0.0000382160,
374 0.0000403941, 0.0000426554, 0.0000450017, 0.0000474347, 0.0000499561, 0.0000525678,
375 0.0000552715, 0.0000580689, 0.0000609619, 0.0000639522, 0.0000670417, 0.0000702323,
376 0.0000735256, 0.0000769237, 0.0000804283, 0.0000840414, 0.0000877648, 0.0000916004,
377 0.0000955503, 0.0000996163, 0.0001038003, 0.0001081045, 0.0001125306, 0.0001170808,
378 0.0001217571, 0.0001265615, 0.0001314961, 0.0001365630, 0.0001417641, 0.0001471018,
379 0.0001525781, 0.0001581951, 0.0001639551, 0.0001698602, 0.0001759127, 0.0001821148,
380 0.0001884687, 0.0001949769, 0.0002016415, 0.0002084650, 0.0002154497, 0.0002225981,
381 0.0002299124, 0.0002373952, 0.0002450491, 0.0002528764, 0.0002608798, 0.0002690618,
382 0.0002774250, 0.0002859721, 0.0002947058, 0.0003036288, 0.0003127438, 0.0003220537,
383 0.0003315613, 0.0003412695, 0.0003511812, 0.0003612993, 0.0003716269, 0.0003821671,
384 0.0003929230, 0.0004038976, 0.0004150944, 0.0004265164, 0.0004381671, 0.0004500499,
385 0.0004621681, 0.0004745254, 0.0004871253, 0.0004999714, 0.0005130674, 0.0005264172,
386 0.0005400246, 0.0005538936, 0.0005680282, 0.0005824325, 0.0005971108, 0.0006120672,
387 0.0006273063, 0.0006428325, 0.0006586504, 0.0006747647, 0.0006911803, 0.0007079021,
388 0.0007249351, 0.0007422845, 0.0007599558, 0.0007779542, 0.0007962855, 0.0008149555,
389 0.0008339699, 0.0008533351, 0.0008730571, 0.0008931425, 0.0009135979, 0.0009344302,
390 0.0009556463, 0.0009772537, 0.0009992598, 0.0010216720, 0.0010444990, 0.0010677490,
391 0.0010914300, 0.0011155510, 0.0011401210, 0.0011651510, 0.0011906480, 0.0012166250,
392 0.0012430910, 0.0012700570, 0.0012975350, 0.0013255370, 0.0013540740, 0.0013831600,
393 0.0014128090, 0.0014430330, 0.0014738470, 0.0015052670, 0.0015373090, 0.0015699870,
394 0.0016033210, 0.0016373270, 0.0016720240, 0.0017074330, 0.0017435730, 0.0017804660,
395 0.0018181340, 0.0018566010, 0.0018958920, 0.0019360320, 0.0019770490, 0.0020189710,
396 0.0020618280, 0.0021056510, 0.0021504740, 0.0021963300, 0.0022432570, 0.0022912920,
397 0.0023404760, 0.0023908500, 0.0024424590, 0.0024953480, 0.0025495680, 0.0026051670,
398 0.0026622000, 0.0027207220, 0.0027807900, 0.0028424670, 0.0029058130, 0.0029708960,
399 0.0030377810, 0.0031065390, 0.0031772420, 0.0032499620, 0.0033247750, 0.0034017560,
400 0.0034809810, 0.0035625250, 0.0036464650, 0.0037328720, 0.0038218190, 0.0039133710,
401 0.0040075910, 0.0041045360, 0.0042042540, 0.0043067850, 0.0044121600, 0.0045203990,
402 0.0046315090, 0.0047454840, 0.0048623040, 0.0049819360, 0.0051043310, 0.0052294280,
403 0.0053571500, 0.0054874080, 0.0056201030, 0.0057551230, 0.0058923500, 0.0060316590,
404 0.0061729200, 0.0063159990, 0.0064607640, 0.0066070800, 0.0067548180, 0.0069038480,
405 0.0070540500, 0.0072053060, 0.0073575050, 0.0075105450, 0.0076643270, 0.0078187630,
406 0.0079737720, 0.0081292780, 0.0082852130, 0.0084415170, 0.0085981350, 0.0087550200,
407 0.0089121270, 0.0090694200, 0.0092268660, 0.0093844380, 0.0095421120, 0.0096998670,
408 0.0098576870, 0.0100155600, 0.0101734700, 0.0103314200, 0.0104894000, 0.0106474000,
409 0.0108054200, 0.0109634700, 0.0111215500, 0.0112796500, 0.0114377900, 0.0115959700,
410 0.0117541900, 0.0119124600, 0.0120707900, 0.0122292000, 0.0123876800, 0.0125462500,
411 0.0127049200, 0.0128637100, 0.0130226100, 0.0131816400, 0.0133408200, 0.0135001600,
412 0.0136596700, 0.0138193500, 0.0139792400, 0.0141393300, 0.0142996400, 0.0144601800,
413 0.0146209700, 0.0147820200, 0.0149433400, 0.0151049500, 0.0152668600, 0.0154290800,
414 0.0155916300, 0.0157545200, 0.0159177700, 0.0160813900, 0.0162453900, 0.0164097800,
415 0.0165745900, 0.0167398200, 0.0169054800, 0.0170716000, 0.0172381900, 0.0174052600,
416 0.0175728200, 0.0177408900, 0.0179094800, 0.0180786100, 0.0182483000, 0.0184185500,
417 0.0185893800, 0.0187608000, 0.0189328400, 0.0191054900, 0.0192787900, 0.0194527400,
418 0.0196273600, 0.0198026700, 0.0199786700, 0.0201553900, 0.0203328300, 0.0205110200,
419 0.0206899700, 0.0208696900, 0.0210502100, 0.0212315300, 0.0214136700, 0.0215966500,
420 0.0217804800, 0.0219651900, 0.0221507800, 0.0223372700, 0.0225246900, 0.0227130400,
421 0.0229023400, 0.0230926100, 0.0232838700, 0.0234761300, 0.0236694200, 0.0238637400,
422 0.0240591200, 0.0242555800, 0.0244531300, 0.0246517900, 0.0248515800, 0.0250525200,
423 0.0252546300, 0.0254579200, 0.0256624200, 0.0258681400, 0.0260751100, 0.0262833500,
424 0.0264928700, 0.0267036900, 0.0269158400, 0.0271293400, 0.0273442100, 0.0275604600,
425 0.0277781300, 0.0279972300, 0.0282177800, 0.0284398100, 0.0286633400, 0.0288884000,
426 0.0291149900, 0.0293431600, 0.0295729200, 0.0298042900, 0.0300373100, 0.0302719900,
427 0.0305083600, 0.0307464500, 0.0309862700, 0.0312278600, 0.0314712500, 0.0317164500,
428 0.0319635000, 0.0322124200, 0.0324632400, 0.0327159900, 0.0329706900, 0.0332273800,
429 0.0334860800, 0.0337468200, 0.0340096300, 0.0342745400, 0.0345415800, 0.0348107800,
430 0.0350821700, 0.0353557800, 0.0356316400, 0.0359097900, 0.0361902600, 0.0364730800,
431 0.0367582700, 0.0370458900, 0.0373359500, 0.0376284900, 0.0379235500, 0.0382211500,
432 0.0385213500, 0.0388241600, 0.0391296300, 0.0394377800, 0.0397486700, 0.0400623100,
433 0.0403787600, 0.0406980400, 0.0410202000, 0.0413452600, 0.0416732700, 0.0420042700,
434 0.0423382900, 0.0426753700, 0.0430155500, 0.0433588700, 0.0437053700, 0.0440550800,
435 0.0444080500, 0.0447643100, 0.0451239000, 0.0454868700, 0.0458532400, 0.0462230700,
436 0.0465963900, 0.0469732300, 0.0473536500, 0.0477376700, 0.0481253400, 0.0485166900,
437 0.0489117600, 0.0493106000, 0.0497132400, 0.0501197100, 0.0505300600, 0.0509443200,
438 0.0513625300, 0.0517847200, 0.0522109400, 0.0526412100, 0.0530755700, 0.0535140600,
439 0.0539567000, 0.0544035400, 0.0548545900, 0.0553099000, 0.0557694900, 0.0562334000,
440 0.0567016400, 0.0571742500, 0.0576512600, 0.0581326800, 0.0586185500, 0.0591088800,
441 0.0596036900, 0.0601030200, 0.0606068600, 0.0611152500, 0.0616281900, 0.0621457100,
442 0.0626678100, 0.0631945000, 0.0637258000, 0.0642617000, 0.0648022200, 0.0653473700,
443 0.0658971300, 0.0664515200, 0.0670105300, 0.0675741500, 0.0681423900, 0.0687152300,
444 0.0692926600, 0.0698746700, 0.0704612600, 0.0710524000, 0.0716480700, 0.0722482600,
445 0.0728529400, 0.0734620900, 0.0740756900, 0.0746937100, 0.0753161100, 0.0759428600,
446 0.0765739400, 0.0772093000, 0.0778489100, 0.0784927300, 0.0791407200, 0.0797928300,
447 0.0804490200, 0.0811092400, 0.0817734600, 0.0824416000, 0.0831136300, 0.0837895000,
448 0.0844691400, 0.0851525100, 0.0858395400, 0.0865301900, 0.0872243800, 0.0879220700,
449 0.0886231900, 0.0893276700, 0.0900354700, 0.0907465100, 0.0914607200, 0.0921780600,
450 0.0928984400, 0.0936218200, 0.0943481100, 0.0950772500, 0.0958091900, 0.1023867000,
451 0.1055505000, 0.1081475000, 0.1104523000, 0.1125711000, 0.1145581000, 0.1164458000,
452 0.1182550000, 0.1200004000, 0.1216926000, 0.1233395000, 0.1249474000, 0.1265212000,
453 0.1280649000, 0.1295817000, 0.1310743000, 0.1325452000, 0.1339961000, 0.1354288000,
454 0.1368448000, 0.1382454000, 0.1396316000, 0.1410044000, 0.1423647000, 0.1437133000,
455 0.1450510000, 0.1463782000, 0.1476956000, 0.1490038000, 0.1503031000, 0.1515941000,
456 0.1528770000, 0.1541523000, 0.1554203000, 0.1566812000, 0.1579354000, 0.1591832000,
457 0.1604247000, 0.1616602000, 0.1628899000, 0.1641140000, 0.1653326000, 0.1665460000,
458 0.1677543000, 0.1689576000, 0.1701561000, 0.1713499000, 0.1725391000, 0.1737239000,
459 0.1749043000, 0.1760805000, 0.1772526000, 0.1784206000, 0.1795847000, 0.1807449000,
460 0.1819014000, 0.1830541000, 0.1842031000, 0.1853487000, 0.1864907000, 0.1876292000,
461 0.1887644000, 0.1898963000, 0.1910250000, 0.1921504000, 0.1932727000, 0.1943919000,
462 0.1955080000, 0.1966211000, 0.1977313000, 0.1988386000, 0.1999430000, 0.2010446000,
463 0.2021434000, 0.2032395000, 0.2043329000, 0.2054236000, 0.2065117000, 0.2075972000,
464 0.2086801000, 0.2097605000, 0.2108384000, 0.2119139000, 0.2129870000, 0.2140576000,
465 0.2151259000, 0.2161919000, 0.2172555000, 0.2183169000, 0.2193761000, 0.2204330000,
466 0.2214878000, 0.2225404000, 0.2235909000, 0.2246392000, 0.2256855000, 0.2267297000,
467 0.2277719000, 0.2288121000, 0.2298503000, 0.2308866000, 0.2319210000, 0.2329534000,
468 0.2339840000, 0.2350126000, 0.2360395000, 0.2370646000, 0.2380878000, 0.2391093000,
469 0.2401291000, 0.2411471000, 0.2421634000, 0.2431781000, 0.2441911000, 0.2452024000,
470 0.2462122000, 0.2472203000, 0.2482269000, 0.2492319000, 0.2502354000, 0.2512374000,
471 0.2522379000, 0.2532369000, 0.2542345000, 0.2552306000, 0.2562253000, 0.2572187000,
472 0.2582106000, 0.2592012000, 0.2601905000, 0.2611784000, 0.2621651000, 0.2631504000,
473 0.2641345000, 0.2651174000, 0.2660991000, 0.2670795000, 0.2680587000, 0.2690368000,
474 0.2700137000, 0.2709895000, 0.2719642000, 0.2729378000, 0.2739103000, 0.2748817000,
475 0.2758521000, 0.2768215000, 0.2777898000, 0.2787572000, 0.2797235000, 0.2806890000,
476 0.2816534000, 0.2826170000, 0.2835796000, 0.2845414000, 0.2855022000, 0.2864623000,
477 0.2874214000, 0.2883798000, 0.2893373000, 0.2902940000, 0.2912500000, 0.2922052000,
478 0.2931597000, 0.2941134000, 0.2950664000, 0.2960187000, 0.2969704000, 0.2979213000,
479 0.2988716000, 0.2998213000, 0.3007704000, 0.3017189000, 0.3026667000, 0.3036140000,
480 0.3045608000, 0.3055070000, 0.3064527000, 0.3073978000, 0.3083425000, 0.3092867000,
481 0.3102304000, 0.3111737000, 0.3121166000, 0.3130590000, 0.3140010000, 0.3149426000,
482 0.3158839000, 0.3168248000, 0.3177653000, 0.3187056000, 0.3196455000, 0.3205851000,
483 0.3215244000, 0.3224634000, 0.3234022000, 0.3243408000, 0.3252791000, 0.3262172000,
484 0.3271551000, 0.3280928000, 0.3290303000, 0.3299677000, 0.3309050000, 0.3318421000,
485 0.3327791000, 0.3337160000, 0.3346528000, 0.3355896000, 0.3365263000, 0.3374629000,
486 0.3383995000, 0.3393362000, 0.3402728000, 0.3412094000, 0.3421460000, 0.3430827000,
487 0.3440194000, 0.3449562000, 0.3458931000, 0.3468301000, 0.3477672000, 0.3487044000,
488 0.3496418000, 0.3505793000, 0.3515170000, 0.3524548000, 0.3533929000, 0.3543311000,
489 0.3552696000, 0.3562083000, 0.3571472000, 0.3580864000, 0.3590259000, 0.3599657000,
490 0.3609057000, 0.3618461000, 0.3627868000, 0.3637279000, 0.3646693000, 0.3656111000,
491 0.3665532000, 0.3674958000, 0.3684387000, 0.3693821000, 0.3703259000, 0.3712701000,
492 0.3722149000, 0.3731600000, 0.3741057000, 0.3750519000, 0.3759986000, 0.3769458000,
493 0.3778936000, 0.3788419000, 0.3797908000, 0.3807402000, 0.3816903000, 0.3826409000,
494 0.3835922000, 0.3845441000, 0.3854966000, 0.3864498000, 0.3874037000, 0.3883583000,
495 0.3893135000, 0.3902695000, 0.3912262000, 0.3921836000, 0.3931418000, 0.3941007000,
496 0.3950604000, 0.3960209000, 0.3969822000, 0.3979443000, 0.3989072000, 0.3998710000,
497 0.4008356000, 0.4018011000, 0.4027675000, 0.4037347000, 0.4047029000, 0.4056720000,
498 0.4066420000, 0.4076129000, 0.4085848000, 0.4095576000, 0.4105315000, 0.4115063000,
499 0.4124821000, 0.4134590000, 0.4144368000, 0.4154157000, 0.4163957000, 0.4173767000,
500 0.4183588000, 0.4193420000, 0.4203264000, 0.4213118000, 0.4222983000, 0.4232860000,
501 0.4242749000, 0.4252649000, 0.4262561000, 0.4272484000, 0.4282420000, 0.4292368000,
502 0.4302328000, 0.4312301000, 0.4322286000, 0.4332284000, 0.4342295000, 0.4352318000,
503 0.4362354000, 0.4372404000, 0.4382467000, 0.4392543000, 0.4402633000, 0.4412736000,
504 0.4422853000, 0.4432984000, 0.4443129000, 0.4453288000, 0.4463461000, 0.4473648000,
505 0.4483850000, 0.4494067000, 0.4504298000, 0.4514544000, 0.4524804000, 0.4535080000,
506 0.4545371000, 0.4555678000, 0.4565999000, 0.4576337000, 0.4586689000, 0.4597058000,
507 0.4607442000, 0.4617843000, 0.4628259000, 0.4638692000, 0.4649141000, 0.4659606000,
508 0.4670088000, 0.4680586000, 0.4691101000, 0.4701633000, 0.4712183000, 0.4722749000,
509 0.4733332000, 0.4743933000, 0.4754551000, 0.4765186000, 0.4775840000, 0.4786511000,
510 0.4797199000, 0.4807906000, 0.4818631000, 0.4829374000, 0.4840135000, 0.4850915000,
511 0.4861713000, 0.4872530000, 0.4883366000, 0.4894220000, 0.4905093000, 0.4915985000,
512 0.4926897000, 0.4937827000, 0.4948777000, 0.4959747000, 0.4970736000, 0.4981744000,
513 0.4992772000, 0.5003820000, 0.5014888000, 0.5025976000, 0.5037085000, 0.5048213000,
514 0.5059362000, 0.5070531000, 0.5081720000, 0.5092931000, 0.5104162000, 0.5115413000,
515 0.5126686000, 0.5137980000, 0.5149294000, 0.5160630000, 0.5171987000, 0.5183366000,
516 0.5194766000, 0.5206187000, 0.5217630000, 0.5229095000, 0.5240581000, 0.5252090000,
517 0.5263620000, 0.5275173000, 0.5286748000, 0.5298344000, 0.5309964000, 0.5321605000,
518 0.5333269000, 0.5344956000, 0.5356665000, 0.5368397000, 0.5380152000, 0.5391930000,
519 0.5403731000, 0.5415554000, 0.5427401000, 0.5439272000, 0.5451165000, 0.5463082000,
520 0.5475022000, 0.5486986000, 0.5498973000, 0.5510984000, 0.5523019000, 0.5535077000,
521 0.5547159000, 0.5559266000, 0.5571396000, 0.5583551000, 0.5595729000, 0.5607932000,
522 0.5620159000, 0.5632411000, 0.5644686000, 0.5656987000, 0.5669312000, 0.5681661000,
523 0.5694036000, 0.5706435000, 0.5718859000, 0.5731307000, 0.5743781000, 0.5756280000,
524 0.5768804000, 0.5781352000, 0.5793926000, 0.5806526000, 0.5819150000, 0.5831800000,
525 0.5844476000, 0.5857176000, 0.5869903000, 0.5882655000, 0.5895432000, 0.5908235000,
526 0.5921064000, 0.5933919000, 0.5946799000, 0.5959706000, 0.5972638000, 0.5985596000,
527 0.5998581000, 0.6011591000, 0.6024627000, 0.6037690000, 0.6050779000, 0.6063894000,
528 0.6077036000, 0.6090203000, 0.6103397000, 0.6116618000, 0.6129865000, 0.6143139000,
529 0.6156439000, 0.6169765000, 0.6183119000, 0.6196499000, 0.6209905000, 0.6223339000,
530 0.6236799000, 0.6250286000, 0.6263800000, 0.6277340000, 0.6290908000, 0.6304503000,
531 0.6318124000, 0.6331773000, 0.6345448000, 0.6359151000, 0.6372881000, 0.6386638000,
532 0.6400422000, 0.6414233000, 0.6428071000, 0.6441937000, 0.6455830000, 0.6469750000,
533 0.6483698000, 0.6497673000, 0.6511675000, 0.6525705000, 0.6539762000, 0.6553846000,
534 0.6567958000, 0.6582098000, 0.6596264000, 0.6610459000, 0.6624681000, 0.6638930000,
535 0.6653208000, 0.6667512000, 0.6681845000, 0.6696205000, 0.6710592000, 0.6725007000,
536 0.6739450000, 0.6753921000, 0.6768419000, 0.6782945000, 0.6797499000, 0.6812080000,
537 0.6826690000, 0.6841327000, 0.6855992000, 0.6870684000, 0.6885404000, 0.6900153000,
538 0.6914929000, 0.6929732000, 0.6944564000, 0.6959423000, 0.6974311000, 0.6989226000,
539 0.7004169000, 0.7019140000, 0.7034139000, 0.7049165000, 0.7064220000, 0.7079302000,
540 0.7094412000, 0.7109550000, 0.7124716000, 0.7139910000, 0.7155132000, 0.7170382000,
541 0.7185659000, 0.7200965000, 0.7216298000, 0.7231660000, 0.7247049000, 0.7262466000,
542 0.7277911000, 0.7293383000, 0.7308884000, 0.7324413000, 0.7339969000, 0.7355553000,
543 0.7371166000, 0.7386806000, 0.7402474000, 0.7418169000, 0.7433893000, 0.7449644000,
544 0.7465424000, 0.7481231000, 0.7497066000, 0.7512929000, 0.7528819000, 0.7544737000,
545 0.7560684000, 0.7576657000, 0.7592659000, 0.7608689000, 0.7624746000, 0.7640831000,
546 0.7656943000, 0.7673084000, 0.7689252000, 0.7705448000, 0.7721671000, 0.7737922000,
547 0.7754201000, 0.7770507000, 0.7786841000, 0.7803203000, 0.7819592000, 0.7836009000,
548 0.7852454000, 0.7868926000, 0.7885425000, 0.7901952000, 0.7918507000, 0.7935089000,
549 0.7951698000, 0.7968335000, 0.7985000000, 0.8001692000, 0.8018411000, 0.8035158000,
550 0.8051932000, 0.8068733000, 0.8085562000, 0.8102418000, 0.8119301000, 0.8136212000,
551 0.8153150000, 0.8170115000, 0.8187108000, 0.8204127000, 0.8221174000, 0.8238248000,
552 0.8255349000, 0.8272477000, 0.8289633000, 0.8306815000, 0.8324025000, 0.8341261000,
553 0.8358525000, 0.8375815000, 0.8393133000, 0.8410477000, 0.8427849000, 0.8445247000,
554 0.8462672000, 0.8480124000, 0.8497603000, 0.8515108000, 0.8532641000, 0.8550200000,
555 0.8567786000, 0.8585399000, 0.8603038000, 0.8620704000, 0.8638396000, 0.8656116000,
556 0.8673861000, 0.8691634000, 0.8709433000, 0.8727258000, 0.8745110000, 0.8762988000,
557 0.8780893000, 0.8798824000, 0.8816782000, 0.8834766000, 0.8852776000, 0.8870812000,
558 0.8888875000, 0.8906964000, 0.8925079000, 0.8943221000, 0.8961388000, 0.8979582000,
559 0.8997801000, 0.9016047000, 0.9034319000, 0.9052617000, 0.9070941000, 0.9089291000,
560 0.9107666000, 0.9126068000, 0.9144495000, 0.9162949000, 0.9181428000, 0.9199933000,
561 0.9218463000, 0.9237020000, 0.9255602000, 0.9274210000, 0.9292843000, 0.9311502000,
562 0.9330186000, 0.9348896000, 0.9367632000, 0.9386393000, 0.9405179000, 0.9423991000,
563 0.9442828000, 0.9461691000, 0.9480579000, 0.9499492000, 0.9518431000, 0.9537394000,
564 0.9556383000, 0.9575397000, 0.9594436000, 0.9613501000, 0.9632590000, 0.9651704000,
565 0.9670844000, 0.9690008000, 0.9709197000, 0.9728411000, 0.9747650000, 0.9766914000,
566 0.9786203000, 0.9805516000, 0.9824854000, 0.9844217000, 0.9863605000, 0.9883017000,
567 0.9902454000, 0.9921915000, 0.9941401000, 0.9960911000, 0.9980446000, 1.0000010000,
568 1.0019590000, 1.0039200000, 1.0058830000, 1.0078490000, 1.0098170000, 1.0117870000,
569 1.0137600000, 1.0157360000, 1.0177130000, 1.0196940000, 1.0216760000, 1.0236610000,
570 1.0256490000, 1.0276380000, 1.0296310000, 1.0316250000, 1.0336220000, 1.0356210000,
571 1.0376230000, 1.0396270000, 1.0416340000, 1.0436430000, 1.0456540000, 1.0476670000,
572 1.0496830000, 1.0517020000, 1.0537220000, 1.0557450000, 1.0577710000, 1.0597980000,
573 1.0618280000, 1.0638610000, 1.0658950000, 1.0679320000, 1.0699720000, 1.0720130000,
574 1.0740570000, 1.0761040000, 1.0781520000, 1.0802030000, 1.0822560000, 1.0843120000,
575 1.0863690000, 1.0884300000, 1.0904920000, 1.0925570000, 1.0946240000, 1.0966930000,
576 1.0987640000, 1.1008380000, 1.1029140000, 1.1049920000, 1.1070730000, 1.1091560000,
577 1.1112410000, 1.1133280000, 1.1154180000, 1.1175090000, 1.1196030000, 1.1217000000,
578 1.1237980000, 1.1258990000, 1.1280020000, 1.1301070000, 1.1322140000, 1.1343240000,
579 1.1364360000, 1.1385500000, 1.1406660000, 1.1427850000, 1.1449050000, 1.1470280000,
580 1.1491530000, 1.1512800000, 1.1534100000, 1.1555410000, 1.1576750000, 1.1598110000,
581 1.1619490000, 1.1640890000, 1.1662310000, 1.1683760000, 1.1705230000, 1.1726710000,
582 1.1748220000, 1.1769750000, 1.1791310000, 1.1812880000, 1.1834480000, 1.1856090000,
583 1.1877730000, 1.1899390000, 1.1921070000, 1.1942770000, 1.1964490000, 1.1986240000,
584 1.2008000000, 1.2029780000, 1.2051590000, 1.2073420000, 1.2095270000, 1.2117130000,
585 1.2139020000, 1.2160930000, 1.2182860000, 1.2204820000, 1.2226790000, 1.2248780000,
586 1.2270790000, 1.2292830000, 1.2314880000, 1.2336960000, 1.2359050000, 1.2381170000,
587 1.2403300000, 1.2425460000, 1.2447640000, 1.2469830000, 1.2492050000, 1.2514290000,
588 1.2536540000, 1.2558820000, 1.2581120000, 1.2603430000, 1.2625770000, 1.2648130000,
589 1.2670510000, 1.2692900000, 1.2715320000, 1.2737760000, 1.2760210000, 1.2782690000,
590 1.2805180000, 1.2827700000, 1.2850230000, 1.2872790000, 1.2895360000, 1.2917950000,
591 1.2940560000, 1.2963200000, 1.2985850000, 1.3008520000, 1.3031210000, 1.3053920000,
592 1.3076650000, 1.3099390000, 1.3122160000, 1.3144940000, 1.3167750000, 1.3190570000,
593 1.3213420000, 1.3236280000, 1.3259160000, 1.3282060000, 1.3304970000, 1.3327910000,
594 1.3350870000, 1.3373840000, 1.3396830000, 1.3419850000, 1.3442880000, 1.3465920000,
595 1.3488990000, 1.3512080000, 1.3535180000, 1.3558300000, 1.3581440000, 1.3604600000,
596 1.3627780000, 1.3650970000, 1.3674190000, 1.3697420000, 1.3720670000, 1.3743940000,
597 1.3767220000, 1.3790530000, 1.3813850000, 1.3837190000, 1.3860540000, 1.3883920000,
598 1.3907310000, 1.3930720000, 1.3954150000, 1.3977600000, 1.4001060000, 1.4024540000,
599 1.4048040000, 1.4071560000, 1.4095090000, 1.4118640000, 1.4142210000, 1.4165800000,
600 1.4189400000, 1.4213020000, 1.4236660000, 1.4260310000, 1.4283990000, 1.4307670000,
601 1.4331380000, 1.4355100000, 1.4378840000, 1.4402600000, 1.4426370000, 1.4450160000,
602 1.4473970000, 1.4497800000, 1.4521640000, 1.4545490000, 1.4569370000, 1.4593260000,
603 1.4617170000, 1.4641090000, 1.4665030000, 1.4688990000, 1.4712960000, 1.4736950000,
604 1.4760960000, 1.4784980000, 1.4809020000, 1.4833070000, 1.4857140000, 1.4881230000,
605 1.4905330000, 1.4929450000, 1.4953580000, 1.4977730000, 1.5001900000, 1.5026080000,
606 1.5050280000, 1.5074490000, 1.5098720000, 1.5122970000, 1.5147230000, 1.5171500000,
607 1.5195790000, 1.5220100000, 1.5244420000, 1.5268760000, 1.5293110000, 1.5317480000,
608 1.5341870000, 1.5366270000, 1.5390680000, 1.5415110000, 1.5439550000, 1.5464010000,
609 1.5488490000, 1.5512980000, 1.5537480000, 1.5562000000, 1.5586530000, 1.5611080000,
610 1.5635650000, 1.5660230000, 1.5684820000, 1.5709430000, 1.5734050000, 1.5758680000,
611 1.5783340000, 1.5808000000, 1.5832680000, 1.5857380000, 1.5882080000, 1.5906810000,
612 1.5931540000, 1.5956290000, 1.5981060000, 1.6005840000, 1.6030630000, 1.6055440000,
613 1.6080260000, 1.6105100000, 1.6129950000, 1.6154810000, 1.6179680000, 1.6204570000,
614 1.6229480000, 1.6254400000, 1.6279330000, 1.6304270000, 1.6329230000, 1.6354200000,
615 1.6379190000, 1.6404190000, 1.6429200000, 1.6454220000, 1.6479260000, 1.6504310000,
616 1.6529380000, 1.6554450000, 1.6579540000, 1.6604650000, 1.6629760000, 1.6654890000,
617 1.6680040000, 1.6705190000, 1.6730360000, 1.6755540000, 1.6780730000, 1.6805940000,
618 1.6831160000, 1.6856390000, 1.6881630000, 1.6906890000, 1.6932160000, 1.6957440000,
619 1.6982730000, 1.7008030000, 1.7033350000, 1.7058680000, 1.7084020000, 1.7109380000,
620 1.7134740000, 1.7160120000, 1.7185510000, 1.7210910000, 1.7236330000, 1.7261750000,
621 1.7287190000, 1.7312640000, 1.7338100000, 1.7363570000, 1.7389050000, 1.7414550000,
622 1.7440050000, 1.7465570000, 1.7491100000, 1.7516640000, 1.7542190000, 1.7567760000,
623 1.7593330000, 1.7618920000, 1.7644520000, 1.7670120000, 1.7695740000, 1.7721370000,
624 1.7747010000, 1.7772660000, 1.7798330000, 1.7824000000, 1.7849680000, 1.7875380000,
625 1.7901080000, 1.7926800000, 1.7952530000, 1.7978260000, 1.8004010000, 1.8029770000,
626 1.8055540000, 1.8081320000, 1.8107100000, 1.8132900000, 1.8158710000, 1.8184530000,
627 1.8210360000, 1.8236200000, 1.8262050000, 1.8287910000, 1.8313780000, 1.8339660000,
628 1.8365550000, 1.8391450000, 1.8417360000, 1.8443270000, 1.8469200000, 1.8495140000,
629 1.8521090000, 1.8547040000, 1.8573010000, 1.8598980000, 1.8624970000, 1.8650960000,
630 1.8676970000, 1.8702980000, 1.8729000000, 1.8755030000, 1.8781070000, 1.8807120000,
631 1.8833180000, 1.8859250000, 1.8885320000, 1.8911410000, 1.8937500000, 1.8963600000,
632 1.8989710000, 1.9015830000, 1.9041960000, 1.9068100000, 1.9094240000, 1.9120400000,
633 1.9146560000, 1.9172730000, 1.9198910000, 1.9225100000, 1.9251290000, 1.9277500000,
634 1.9303710000, 1.9329930000, 1.9356160000, 1.9382400000, 1.9408640000, 1.9434890000,
635 1.9461150000, 1.9487420000, 1.9513700000, 1.9539980000, 1.9566270000, 1.9592570000,
636 1.9618880000, 1.9645190000, 1.9671520000, 1.9697850000, 1.9724180000, 1.9750530000,
637 1.9776880000, 1.9803240000, 1.9829610000, 1.9855980000, 1.9882360000, 1.9908750000,
638 1.9935140000, 1.9961550000, 1.9987960000, 2.0014370000, 2.0040800000, 2.0067230000,
639 2.0093660000, 2.0120110000, 2.0146560000, 2.0173010000, 2.0199480000, 2.0225950000,
640 2.0252430000, 2.0278910000, 2.0305400000, 2.0331900000, 2.0358400000, 2.0384910000,
641 2.0411430000, 2.0437950000, 2.0464480000, 2.0491010000, 2.0517550000, 2.0544100000,
642 2.0570650000, 2.0597210000, 2.0623780000, 2.0650350000, 2.0676930000, 2.0703510000,
643 2.0730100000, 2.0756690000, 2.0783290000, 2.0809900000, 2.0836510000, 2.0863130000,
644 2.0889750000, 2.0916380000, 2.0943010000, 2.0969650000, 2.0996290000, 2.1022940000,
645 2.1049600000, 2.1076260000, 2.1102930000, 2.1129600000, 2.1156270000, 2.1182950000,
646 2.1209640000, 2.1236330000, 2.1263030000, 2.1289730000, 2.1316440000, 2.1343150000,
647 2.1369860000, 2.1396580000, 2.1423310000, 2.1450040000, 2.1476780000, 2.1503510000,
648 2.1530260000, 2.1557010000 };
649
650double K1400_width[] = {
651 0.0000000007, 0.0000000025, 0.0000000061, 0.0000000121, 0.0000000213, 0.0000000344,
652 0.0000000519, 0.0000000747, 0.0000001034, 0.0000001389, 0.0000001818, 0.0000002329,
653 0.0000002930, 0.0000003630, 0.0000004435, 0.0000005356, 0.0000006400, 0.0000007576,
654 0.0000008892, 0.0000010359, 0.0000011985, 0.0000013780, 0.0000015753, 0.0000017914,
655 0.0000020274, 0.0000022842, 0.0000025629, 0.0000028646, 0.0000031903, 0.0000035413,
656 0.0000039185, 0.0000043232, 0.0000047566, 0.0000052199, 0.0000057143, 0.0000062411,
657 0.0000068015, 0.0000073970, 0.0000080288, 0.0000086983, 0.0000094069, 0.0000101561,
658 0.0000109474, 0.0000117821, 0.0000126619, 0.0000135884, 0.0000145631, 0.0000155876,
659 0.0000166637, 0.0000177930, 0.0000189773, 0.0000202184, 0.0000215182, 0.0000228784,
660 0.0000243010, 0.0000257879, 0.0000273412, 0.0000289629, 0.0000306551, 0.0000324199,
661 0.0000342596, 0.0000361763, 0.0000381724, 0.0000402502, 0.0000424121, 0.0000446606,
662 0.0000469982, 0.0000494275, 0.0000519511, 0.0000545718, 0.0000572923, 0.0000601155,
663 0.0000630443, 0.0000660816, 0.0000692306, 0.0000724943, 0.0000758761, 0.0000793791,
664 0.0000830068, 0.0000867626, 0.0000906501, 0.0000946730, 0.0000988349, 0.0001031398,
665 0.0001075916, 0.0001121943, 0.0001169520, 0.0001218691, 0.0001269499, 0.0001321989,
666 0.0001376208, 0.0001432202, 0.0001490021, 0.0001549715, 0.0001611335, 0.0001674933,
667 0.0001740566, 0.0001808287, 0.0001878156, 0.0001950231, 0.0002024573, 0.0002101244,
668 0.0002180310, 0.0002261837, 0.0002345892, 0.0002432546, 0.0002521872, 0.0002613945,
669 0.0002708840, 0.0002806638, 0.0002907420, 0.0003011271, 0.0003118276, 0.0003228526,
670 0.0003342113, 0.0003459132, 0.0003579682, 0.0003703863, 0.0003831781, 0.0003963544,
671 0.0004099263, 0.0004239055, 0.0004383037, 0.0004531333, 0.0004684072, 0.0004841383,
672 0.0005003405, 0.0005170277, 0.0005342145, 0.0005519161, 0.0005701480, 0.0005889265,
673 0.0006082683, 0.0006281907, 0.0006487119, 0.0006698504, 0.0006916256, 0.0007140575,
674 0.0007371670, 0.0007609758, 0.0007855062, 0.0008107815, 0.0008368260, 0.0008636648,
675 0.0008913240, 0.0009198309, 0.0009492137, 0.0009795018, 0.0010107260, 0.0010429180,
676 0.0010761110, 0.0011103390, 0.0011456390, 0.0011820480, 0.0012196060, 0.0012583530,
677 0.0012983320, 0.0013395870, 0.0013821650, 0.0014261150, 0.0014714870, 0.0015183340,
678 0.0015667120, 0.0016166800, 0.0016682970, 0.0017216280, 0.0017767390, 0.0018337000,
679 0.0018925850, 0.0019534700, 0.0020164360, 0.0020815670, 0.0021489530, 0.0022186860,
680 0.0022908640, 0.0023655910, 0.0024429750, 0.0025231290, 0.0026061740, 0.0026922350,
681 0.0027814460, 0.0028739440, 0.0029698790, 0.0030694050, 0.0031726840, 0.0032798900,
682 0.0033912030, 0.0035068140, 0.0036269240, 0.0037517470, 0.0038815040, 0.0040164320,
683 0.0041567790, 0.0043028050, 0.0044547850, 0.0046130090, 0.0047777790, 0.0049494150,
684 0.0051282500, 0.0053146360, 0.0055089370, 0.0057115350, 0.0059228300, 0.0061432340,
685 0.0063731760, 0.0066130990, 0.0068634590, 0.0071247200, 0.0073973600, 0.0076818580,
686 0.0079786990, 0.0082883670, 0.0086113370, 0.0089480750, 0.0092990290, 0.0096646210,
687 0.0100452400, 0.0104412500, 0.0108529300, 0.0112805500, 0.0117242700, 0.0121842100,
688 0.0126603900, 0.0131527500, 0.0136611300, 0.0141852900, 0.0147248800, 0.0152794700,
689 0.0158485100, 0.0164314100, 0.0170274700, 0.0176359300, 0.0182559700, 0.0188867400,
690 0.0195273500, 0.0201768900, 0.0208344400, 0.0214991000, 0.0221699600, 0.0228461600,
691 0.0235268400, 0.0242112100, 0.0248985100, 0.0255880100, 0.0262790500, 0.0269710000,
692 0.0276632900, 0.0283554000, 0.0290468400, 0.0297371800, 0.0304260300, 0.0311130400,
693 0.0317978800, 0.0324802900, 0.0331600000, 0.0338368100, 0.0345105300, 0.0351809800,
694 0.0358480300, 0.0365115500, 0.0371714500, 0.0378276300, 0.0384800200, 0.0391285600,
695 0.0397732200, 0.0404139500, 0.0410507400, 0.0416835500, 0.0423124000, 0.0429372700,
696 0.0435581700, 0.0441751100, 0.0447881100, 0.0453971800, 0.0460023400, 0.0466036300,
697 0.0472010600, 0.0477946700, 0.0483845000, 0.0489705600, 0.0495529100, 0.0501315800,
698 0.0507066000, 0.0512780200, 0.0518458700, 0.0524101800, 0.0529710200, 0.0535284000,
699 0.0540823700, 0.0546329800, 0.0551802500, 0.0557242400, 0.0562649800, 0.0568025200,
700 0.0573368800, 0.0578681100, 0.0583962500, 0.0589213300, 0.0594434000, 0.0599624900,
701 0.0604786400, 0.0609918800, 0.0615022500, 0.0620097800, 0.0625145200, 0.0630164900,
702 0.0635157300, 0.0640122700, 0.0645061500, 0.0649973900, 0.0654860300, 0.0659721100,
703 0.0664556400, 0.0669366600, 0.0674152100, 0.0678913000, 0.0683649800, 0.0688362600,
704 0.0693051800, 0.0697717600, 0.0702360300, 0.0706980200, 0.0711577500, 0.0716152500,
705 0.0720705400, 0.0725236500, 0.0729746100, 0.0734234300, 0.0738701400, 0.0743147600,
706 0.0747573300, 0.0751978500, 0.0756363500, 0.0760728500, 0.0765073800, 0.0769399600,
707 0.0773706000, 0.0777993300, 0.0782261600, 0.0786511200, 0.0790742300, 0.0794955000,
708 0.0799149600, 0.0803326100, 0.0807484900, 0.0811626100, 0.0815749800, 0.0819856300,
709 0.0823945700, 0.0828018200, 0.0832073900, 0.0836113000, 0.0840135700, 0.0844142200,
710 0.0848132500, 0.0852106900, 0.0856065500, 0.0860008400, 0.0863935900, 0.0867848000,
711 0.0871744900, 0.0875626700, 0.0879493600, 0.0883345700, 0.0887183200, 0.0891006200,
712 0.0894814800, 0.0898609100, 0.0902389300, 0.0906155600, 0.0909908000, 0.0913646700,
713 0.0917371700, 0.0921083300, 0.0924781500, 0.0928466500, 0.0932138400, 0.0935797200,
714 0.0939443200, 0.0943076400, 0.0946696900, 0.0950304800, 0.0953900400, 0.0957483600,
715 0.0961054500, 0.0964613400, 0.0968160300, 0.0971695200, 0.0975218300, 0.0978729800,
716 0.0982229600, 0.0985718000, 0.0989194900, 0.0992660600, 0.0996115100, 0.0999558400,
717 0.1002991000, 0.1006412000, 0.1009823000, 0.1013223000, 0.1016612000, 0.1019991000,
718 0.1023359000, 0.1026717000, 0.1030064000, 0.1033402000, 0.1036729000, 0.1040047000,
719 0.1043354000, 0.1046652000, 0.1049939000, 0.1053218000, 0.1056486000, 0.1059745000,
720 0.1062995000, 0.1066235000, 0.1069466000, 0.1072688000, 0.1075901000, 0.1079104000,
721 0.1082299000, 0.1085485000, 0.1088662000, 0.1091831000, 0.1094990000, 0.1098142000,
722 0.1101284000, 0.1104418000, 0.1107544000, 0.1110662000, 0.1113771000, 0.1116873000,
723 0.1119966000, 0.1123051000, 0.1126128000, 0.1129198000, 0.1132259000, 0.1135313000,
724 0.1138360000, 0.1141398000, 0.1144429000, 0.1147453000, 0.1150469000, 0.1153478000,
725 0.1156480000, 0.1159474000, 0.1162461000, 0.1165442000, 0.1168415000, 0.1171381000,
726 0.1174340000, 0.1177292000, 0.1180238000, 0.1183177000, 0.1186109000, 0.1189034000,
727 0.1191953000, 0.1194865000, 0.1197771000, 0.1200670000, 0.1203563000, 0.1206449000,
728 0.1209329000, 0.1212203000, 0.1215071000, 0.1217933000, 0.1220788000, 0.1223637000,
729 0.1226480000, 0.1229318000, 0.1232149000, 0.1234974000, 0.1237794000, 0.1240607000,
730 0.1243415000, 0.1246217000, 0.1249013000, 0.1251803000, 0.1254588000, 0.1257367000,
731 0.1260140000, 0.1262908000, 0.1265670000, 0.1268427000, 0.1271178000, 0.1273923000,
732 0.1276663000, 0.1279398000, 0.1282127000, 0.1284850000, 0.1287568000, 0.1290281000,
733 0.1292988000, 0.1295690000, 0.1298387000, 0.1301078000, 0.1303764000, 0.1306445000,
734 0.1309120000, 0.1311790000, 0.1314455000, 0.1317114000, 0.1319768000, 0.1322417000,
735 0.1325060000, 0.1327699000, 0.1330332000, 0.1332960000, 0.1335582000, 0.1338199000,
736 0.1340812000, 0.1343418000, 0.1346020000, 0.1348617000, 0.1351208000, 0.1355387000,
737 0.1358629000, 0.1361711000, 0.1364708000, 0.1367648000, 0.1370546000, 0.1373411000,
738 0.1376249000, 0.1379063000, 0.1381858000, 0.1384634000, 0.1387394000, 0.1390139000,
739 0.1392870000, 0.1395588000, 0.1398294000, 0.1400988000, 0.1403671000, 0.1406344000,
740 0.1409007000, 0.1411659000, 0.1414302000, 0.1416936000, 0.1419562000, 0.1422178000,
741 0.1424786000, 0.1427386000, 0.1429977000, 0.1432561000, 0.1435136000, 0.1437704000,
742 0.1440265000, 0.1442818000, 0.1445364000, 0.1447902000, 0.1450434000, 0.1452958000,
743 0.1455475000, 0.1457986000, 0.1460490000, 0.1462986000, 0.1465477000, 0.1467960000,
744 0.1470437000, 0.1472908000, 0.1475372000, 0.1477830000, 0.1480281000, 0.1482726000,
745 0.1485165000, 0.1487597000, 0.1490024000, 0.1492444000, 0.1494858000, 0.1497266000,
746 0.1499668000, 0.1502065000, 0.1504455000, 0.1506839000, 0.1509218000, 0.1511590000,
747 0.1513957000, 0.1516318000, 0.1518674000, 0.1521023000, 0.1523367000, 0.1525706000,
748 0.1528038000, 0.1530366000, 0.1532687000, 0.1535003000, 0.1537314000, 0.1539619000,
749 0.1541918000, 0.1544213000, 0.1546501000, 0.1548785000, 0.1551063000, 0.1553336000,
750 0.1555603000, 0.1557866000, 0.1560122000, 0.1562374000, 0.1564621000, 0.1566862000,
751 0.1569098000, 0.1571329000, 0.1573555000, 0.1575776000, 0.1577992000, 0.1580203000,
752 0.1582409000, 0.1584609000, 0.1586805000, 0.1588996000, 0.1591182000, 0.1593363000,
753 0.1595539000, 0.1597710000, 0.1599876000, 0.1602038000, 0.1604194000, 0.1606346000,
754 0.1608493000, 0.1610636000, 0.1612773000, 0.1614906000, 0.1617035000, 0.1619158000,
755 0.1621277000, 0.1623391000, 0.1625501000, 0.1627606000, 0.1629706000, 0.1631802000,
756 0.1633894000, 0.1635981000, 0.1638063000, 0.1640141000, 0.1642214000, 0.1644283000,
757 0.1646347000, 0.1648407000, 0.1650463000, 0.1652514000, 0.1654561000, 0.1656604000,
758 0.1658642000, 0.1660676000, 0.1662706000, 0.1664731000, 0.1666752000, 0.1668769000,
759 0.1670782000, 0.1672790000, 0.1674794000, 0.1676794000, 0.1678790000, 0.1680782000,
760 0.1682769000, 0.1684753000, 0.1686732000, 0.1688707000, 0.1690679000, 0.1692646000,
761 0.1694609000, 0.1696568000, 0.1698523000, 0.1700474000, 0.1702421000, 0.1704364000,
762 0.1706304000, 0.1708239000, 0.1710170000, 0.1712098000, 0.1714022000, 0.1715941000,
763 0.1717857000, 0.1719769000, 0.1721677000, 0.1723582000, 0.1725482000, 0.1727379000,
764 0.1729272000, 0.1731162000, 0.1733047000, 0.1734929000, 0.1736807000, 0.1738682000,
765 0.1740552000, 0.1742420000, 0.1744283000, 0.1746143000, 0.1747999000, 0.1749851000,
766 0.1751700000, 0.1753546000, 0.1755388000, 0.1757226000, 0.1759060000, 0.1760892000,
767 0.1762719000, 0.1764543000, 0.1766364000, 0.1768181000, 0.1769995000, 0.1771805000,
768 0.1773612000, 0.1775415000, 0.1777215000, 0.1779011000, 0.1780804000, 0.1782594000,
769 0.1784380000, 0.1786163000, 0.1787943000, 0.1789719000, 0.1791492000, 0.1793262000,
770 0.1795028000, 0.1796791000, 0.1798551000, 0.1800308000, 0.1802061000, 0.1803811000,
771 0.1805558000, 0.1807301000, 0.1809042000, 0.1810779000, 0.1812513000, 0.1814244000,
772 0.1815971000, 0.1817696000, 0.1819417000, 0.1821136000, 0.1822851000, 0.1824563000,
773 0.1826272000, 0.1827978000, 0.1829681000, 0.1831380000, 0.1833077000, 0.1834771000,
774 0.1836461000, 0.1838149000, 0.1839834000, 0.1841515000, 0.1843194000, 0.1844869000,
775 0.1846542000, 0.1848212000, 0.1849879000, 0.1851543000, 0.1853203000, 0.1854862000,
776 0.1856517000, 0.1858169000, 0.1859818000, 0.1861465000, 0.1863108000, 0.1864749000,
777 0.1866387000, 0.1868022000, 0.1869654000, 0.1871284000, 0.1872910000, 0.1874534000,
778 0.1876155000, 0.1877774000, 0.1879389000, 0.1881002000, 0.1882612000, 0.1884219000,
779 0.1885824000, 0.1887426000, 0.1889025000, 0.1890621000, 0.1892215000, 0.1893806000,
780 0.1895395000, 0.1896980000, 0.1898563000, 0.1900144000, 0.1901722000, 0.1903297000,
781 0.1904869000, 0.1906439000, 0.1908007000, 0.1909571000, 0.1911133000, 0.1912693000,
782 0.1914250000, 0.1915805000, 0.1917357000, 0.1918906000, 0.1920453000, 0.1921997000,
783 0.1923539000, 0.1925078000, 0.1926615000, 0.1928149000, 0.1929681000, 0.1931211000,
784 0.1932738000, 0.1934262000, 0.1935784000, 0.1937304000, 0.1938821000, 0.1940336000,
785 0.1941848000, 0.1943358000, 0.1944865000, 0.1946370000, 0.1947873000, 0.1949374000,
786 0.1950872000, 0.1952367000, 0.1953860000, 0.1955351000, 0.1956840000, 0.1958326000,
787 0.1959810000, 0.1961292000, 0.1962771000, 0.1964248000, 0.1965723000, 0.1967196000,
788 0.1968666000, 0.1970134000, 0.1971600000, 0.1973063000, 0.1974524000, 0.1975983000,
789 0.1977440000, 0.1978894000, 0.1980347000, 0.1981797000, 0.1983245000, 0.1984690000,
790 0.1986134000, 0.1987575000, 0.1989015000, 0.1990452000, 0.1991886000, 0.1993319000,
791 0.1994750000, 0.1996178000, 0.1997604000, 0.1999029000, 0.2000451000, 0.2001871000,
792 0.2003289000, 0.2004704000, 0.2006118000, 0.2007530000, 0.2008939000, 0.2010347000,
793 0.2011752000, 0.2013156000, 0.2014557000, 0.2015956000, 0.2017353000, 0.2018749000,
794 0.2020142000, 0.2021533000, 0.2022922000, 0.2024310000, 0.2025695000, 0.2027078000,
795 0.2028459000, 0.2029839000, 0.2031216000, 0.2032591000, 0.2033965000, 0.2035336000,
796 0.2036706000, 0.2038073000, 0.2039439000, 0.2040803000, 0.2042165000, 0.2043524000,
797 0.2044883000, 0.2046239000, 0.2047593000, 0.2048945000, 0.2050296000, 0.2051644000,
798 0.2052991000, 0.2054336000, 0.2055679000, 0.2057020000, 0.2058360000, 0.2059697000,
799 0.2061033000, 0.2062367000, 0.2063699000, 0.2065029000, 0.2066358000, 0.2067684000,
800 0.2069009000, 0.2070332000, 0.2071653000, 0.2072973000, 0.2074291000, 0.2075607000,
801 0.2076921000, 0.2078233000, 0.2079544000, 0.2080853000, 0.2082161000, 0.2083466000,
802 0.2084770000, 0.2086072000, 0.2087372000, 0.2088671000, 0.2089968000, 0.2091263000,
803 0.2092557000, 0.2093849000, 0.2095139000, 0.2096428000, 0.2097715000, 0.2099000000,
804 0.2100284000, 0.2101566000, 0.2102846000, 0.2104125000, 0.2105402000, 0.2106677000,
805 0.2107951000, 0.2109223000, 0.2110494000, 0.2111763000, 0.2113030000, 0.2114296000,
806 0.2115560000, 0.2116823000, 0.2118084000, 0.2119343000, 0.2120601000, 0.2121858000,
807 0.2123112000, 0.2124366000, 0.2125617000, 0.2126868000, 0.2128116000, 0.2129363000,
808 0.2130609000, 0.2131853000, 0.2133095000, 0.2134337000, 0.2135576000, 0.2136814000,
809 0.2138051000, 0.2139286000, 0.2140519000, 0.2141751000, 0.2142982000, 0.2144211000,
810 0.2145439000, 0.2146665000, 0.2147890000, 0.2149113000, 0.2150335000, 0.2151555000,
811 0.2152774000, 0.2153992000, 0.2155208000, 0.2156423000, 0.2157636000, 0.2158848000,
812 0.2160059000, 0.2161268000, 0.2162475000, 0.2163682000, 0.2164887000, 0.2166090000,
813 0.2167292000, 0.2168493000, 0.2169693000, 0.2170891000, 0.2172088000, 0.2173283000,
814 0.2174477000, 0.2175670000, 0.2176861000, 0.2178051000, 0.2179240000, 0.2180427000,
815 0.2181613000, 0.2182798000, 0.2183982000, 0.2185164000, 0.2186345000, 0.2187524000,
816 0.2188703000, 0.2189880000, 0.2191055000, 0.2192230000, 0.2193403000, 0.2194575000,
817 0.2195745000, 0.2196915000, 0.2198083000, 0.2199250000, 0.2200416000, 0.2201580000,
818 0.2202743000, 0.2203905000, 0.2205066000, 0.2206225000, 0.2207384000, 0.2208541000,
819 0.2209697000, 0.2210851000, 0.2212005000, 0.2213157000, 0.2214308000, 0.2215458000,
820 0.2216607000, 0.2217754000, 0.2218901000, 0.2220046000, 0.2221190000, 0.2222333000,
821 0.2223475000, 0.2224615000, 0.2225755000, 0.2226893000, 0.2228030000, 0.2229166000,
822 0.2230301000, 0.2231435000, 0.2232568000, 0.2233699000, 0.2234830000, 0.2235959000,
823 0.2237088000, 0.2238215000, 0.2239341000, 0.2240466000, 0.2241590000, 0.2242713000,
824 0.2243834000, 0.2244955000, 0.2246075000, 0.2247193000, 0.2248311000, 0.2249427000,
825 0.2250543000, 0.2251657000, 0.2252770000, 0.2253883000, 0.2254994000, 0.2256104000,
826 0.2257213000, 0.2258321000, 0.2259429000, 0.2260535000, 0.2261640000, 0.2262744000,
827 0.2263847000, 0.2264949000, 0.2266050000, 0.2267151000, 0.2268250000, 0.2269348000,
828 0.2270445000, 0.2271542000, 0.2272637000, 0.2273731000, 0.2274825000, 0.2275917000,
829 0.2277009000, 0.2278099000, 0.2279189000, 0.2280277000, 0.2281365000, 0.2282452000,
830 0.2283538000, 0.2284623000, 0.2285707000, 0.2286790000, 0.2287872000, 0.2288954000,
831 0.2290034000, 0.2291114000, 0.2292192000, 0.2293270000, 0.2294347000, 0.2295423000,
832 0.2296498000, 0.2297572000, 0.2298646000, 0.2299718000, 0.2300790000, 0.2301861000,
833 0.2302931000, 0.2304000000, 0.2305068000, 0.2306136000, 0.2307202000, 0.2308268000,
834 0.2309333000, 0.2310397000, 0.2311460000, 0.2312523000, 0.2313584000, 0.2314645000,
835 0.2315705000, 0.2316765000, 0.2317823000, 0.2318881000, 0.2319938000, 0.2320994000,
836 0.2322049000, 0.2323104000, 0.2324157000, 0.2325210000, 0.2326263000, 0.2327314000,
837 0.2328365000, 0.2329415000, 0.2330464000, 0.2331513000, 0.2332560000, 0.2333607000,
838 0.2334654000, 0.2335699000, 0.2336744000, 0.2337788000, 0.2338831000, 0.2339874000,
839 0.2340916000, 0.2341957000, 0.2342998000, 0.2344038000, 0.2345077000, 0.2346115000,
840 0.2347153000, 0.2348190000, 0.2349227000, 0.2350262000, 0.2351298000, 0.2352332000,
841 0.2353366000, 0.2354399000, 0.2355431000, 0.2356463000, 0.2357494000, 0.2358525000,
842 0.2359555000, 0.2360584000, 0.2361612000, 0.2362640000, 0.2363668000, 0.2364694000,
843 0.2365721000, 0.2366746000, 0.2367771000, 0.2368795000, 0.2369819000, 0.2370842000,
844 0.2371864000, 0.2372886000, 0.2373908000, 0.2374928000, 0.2375948000, 0.2376968000,
845 0.2377987000, 0.2379005000, 0.2380023000, 0.2381041000, 0.2382057000, 0.2383074000,
846 0.2384089000, 0.2385104000, 0.2386119000, 0.2387133000, 0.2388146000, 0.2389159000,
847 0.2390172000, 0.2391184000, 0.2392195000, 0.2393206000, 0.2394216000, 0.2395226000,
848 0.2396236000, 0.2397245000, 0.2398253000, 0.2399261000, 0.2400268000, 0.2401275000,
849 0.2402281000, 0.2403287000, 0.2404293000, 0.2405298000, 0.2406302000, 0.2407306000,
850 0.2408310000, 0.2409313000, 0.2410316000, 0.2411318000, 0.2412320000, 0.2413321000,
851 0.2414322000, 0.2415323000, 0.2416323000, 0.2417323000, 0.2418322000, 0.2419321000,
852 0.2420319000, 0.2421317000, 0.2422315000, 0.2423312000, 0.2424309000, 0.2425305000,
853 0.2426301000, 0.2427297000, 0.2428292000, 0.2429287000, 0.2430281000, 0.2431275000,
854 0.2432269000, 0.2433262000, 0.2434255000, 0.2435248000, 0.2436240000, 0.2437232000,
855 0.2438224000, 0.2439215000, 0.2440206000, 0.2441196000, 0.2442187000, 0.2443176000,
856 0.2444166000, 0.2445155000, 0.2446144000, 0.2447133000, 0.2448121000, 0.2449109000,
857 0.2450097000, 0.2451084000, 0.2452071000, 0.2453058000, 0.2454044000, 0.2455031000,
858 0.2456017000, 0.2457002000, 0.2457988000, 0.2458973000, 0.2459957000, 0.2460942000,
859 0.2461926000, 0.2462910000, 0.2463894000, 0.2464878000, 0.2465861000, 0.2466844000,
860 0.2467827000, 0.2468809000, 0.2469792000, 0.2470774000, 0.2471756000, 0.2472737000,
861 0.2473719000, 0.2474700000, 0.2475681000, 0.2476662000, 0.2477642000, 0.2478623000,
862 0.2479603000, 0.2480583000, 0.2481563000, 0.2482542000, 0.2483522000, 0.2484501000,
863 0.2485480000, 0.2486459000, 0.2487437000, 0.2488416000, 0.2489394000, 0.2490372000,
864 0.2491351000, 0.2492328000, 0.2493306000, 0.2494284000, 0.2495261000, 0.2496239000,
865 0.2497216000, 0.2498193000, 0.2499170000, 0.2500146000, 0.2501123000, 0.2502100000,
866 0.2503076000, 0.2504052000, 0.2505029000, 0.2506005000, 0.2506981000, 0.2507956000,
867 0.2508932000, 0.2509908000, 0.2510883000, 0.2511859000, 0.2512834000, 0.2513810000,
868 0.2514785000, 0.2515760000, 0.2516735000, 0.2517710000, 0.2518685000, 0.2519660000,
869 0.2520635000, 0.2521610000, 0.2522584000, 0.2523559000, 0.2524534000, 0.2525508000,
870 0.2526483000, 0.2527457000, 0.2528432000, 0.2529406000, 0.2530381000, 0.2531355000,
871 0.2532330000, 0.2533304000, 0.2534278000, 0.2535252000, 0.2536227000, 0.2537201000,
872 0.2538175000, 0.2539150000, 0.2540124000, 0.2541098000, 0.2542073000, 0.2543047000,
873 0.2544021000, 0.2544996000, 0.2545970000, 0.2546944000, 0.2547919000, 0.2548893000,
874 0.2549868000, 0.2550842000, 0.2551817000, 0.2552791000, 0.2553766000, 0.2554740000,
875 0.2555715000, 0.2556690000, 0.2557665000, 0.2558639000, 0.2559614000, 0.2560589000,
876 0.2561564000, 0.2562539000, 0.2563515000, 0.2564490000, 0.2565465000, 0.2566441000,
877 0.2567416000, 0.2568392000, 0.2569367000, 0.2570343000, 0.2571319000, 0.2572295000,
878 0.2573271000, 0.2574247000, 0.2575223000, 0.2576199000, 0.2577176000, 0.2578152000,
879 0.2579129000, 0.2580106000, 0.2581082000, 0.2582059000, 0.2583036000, 0.2584014000,
880 0.2584991000, 0.2585968000, 0.2586946000, 0.2587924000, 0.2588902000, 0.2589879000,
881 0.2590858000, 0.2591836000, 0.2592814000, 0.2593793000, 0.2594771000, 0.2595750000,
882 0.2596729000, 0.2597708000, 0.2598688000, 0.2599667000, 0.2600647000, 0.2601627000,
883 0.2602606000, 0.2603587000, 0.2604567000, 0.2605547000, 0.2606528000, 0.2607509000,
884 0.2608490000, 0.2609471000, 0.2610452000, 0.2611434000, 0.2612415000, 0.2613397000,
885 0.2614379000, 0.2615361000, 0.2616344000, 0.2617326000, 0.2618309000, 0.2619292000,
886 0.2620275000, 0.2621259000, 0.2622242000, 0.2623226000, 0.2624210000, 0.2625194000,
887 0.2626179000, 0.2627164000, 0.2628148000, 0.2629133000, 0.2630119000, 0.2631104000,
888 0.2632090000, 0.2633076000, 0.2634062000, 0.2635048000, 0.2636035000, 0.2637022000,
889 0.2638009000, 0.2638996000, 0.2639984000, 0.2640971000, 0.2641959000, 0.2642947000,
890 0.2643936000, 0.2644925000, 0.2645913000, 0.2646903000, 0.2647892000, 0.2648882000,
891 0.2649871000, 0.2650861000, 0.2651852000, 0.2652842000, 0.2653833000, 0.2654824000,
892 0.2655816000, 0.2656807000, 0.2657799000, 0.2658791000, 0.2659784000, 0.2660776000,
893 0.2661769000, 0.2662762000, 0.2663755000, 0.2664749000, 0.2665743000, 0.2666737000,
894 0.2667731000, 0.2668726000, 0.2669721000, 0.2670716000, 0.2671712000, 0.2672707000,
895 0.2673703000, 0.2674700000, 0.2675696000, 0.2676693000, 0.2677690000, 0.2678687000,
896 0.2679685000, 0.2680683000, 0.2681681000, 0.2682679000, 0.2683678000, 0.2684677000,
897 0.2685676000, 0.2686675000, 0.2687675000, 0.2688675000, 0.2689676000, 0.2690676000,
898 0.2691677000, 0.2692678000, 0.2693680000, 0.2694681000, 0.2695683000, 0.2696685000,
899 0.2697688000, 0.2698691000, 0.2699694000, 0.2700697000, 0.2701701000, 0.2702705000,
900 0.2703709000, 0.2704713000, 0.2705718000, 0.2706723000, 0.2707729000, 0.2708734000,
901 0.2709740000, 0.2710746000, 0.2711753000, 0.2712759000, 0.2713766000, 0.2714774000,
902 0.2715781000, 0.2716789000, 0.2717797000, 0.2718806000, 0.2719814000, 0.2720823000,
903 0.2721833000, 0.2722842000, 0.2723852000, 0.2724862000, 0.2725873000, 0.2726883000,
904 0.2727894000, 0.2728905000, 0.2729917000, 0.2730929000, 0.2731941000, 0.2732953000,
905 0.2733966000, 0.2734979000, 0.2735992000, 0.2737005000, 0.2738019000, 0.2739033000,
906 0.2740048000, 0.2741062000, 0.2742077000, 0.2743092000, 0.2744108000, 0.2745123000,
907 0.2746139000, 0.2747156000, 0.2748172000, 0.2749189000, 0.2750206000, 0.2751223000,
908 0.2752241000, 0.2753259000, 0.2754277000, 0.2755295000, 0.2756314000, 0.2757333000,
909 0.2758352000, 0.2759372000, 0.2760392000, 0.2761412000, 0.2762432000, 0.2763453000,
910 0.2764474000, 0.2765495000, 0.2766516000, 0.2767538000, 0.2768560000, 0.2769582000,
911 0.2770604000, 0.2771627000, 0.2772650000, 0.2773673000, 0.2774697000, 0.2775720000,
912 0.2776745000, 0.2777769000, 0.2778793000, 0.2779818000, 0.2780843000, 0.2781868000,
913 0.2782894000, 0.2783920000, 0.2784946000, 0.2785972000, 0.2786999000, 0.2788026000,
914 0.2789053000, 0.2790080000, 0.2791107000, 0.2792135000, 0.2793163000, 0.2794192000,
915 0.2795220000, 0.2796249000, 0.2797278000, 0.2798307000, 0.2799337000, 0.2800366000,
916 0.2801396000, 0.2802427000, 0.2803457000, 0.2804488000, 0.2805519000, 0.2806550000,
917 0.2807581000, 0.2808613000, 0.2809645000, 0.2810677000, 0.2811709000, 0.2812742000,
918 0.2813774000, 0.2814807000, 0.2815841000, 0.2816874000, 0.2817908000, 0.2818942000,
919 0.2819976000, 0.2821010000, 0.2822044000, 0.2823079000, 0.2824114000, 0.2825149000,
920 0.2826185000, 0.2827220000, 0.2828256000, 0.2829292000, 0.2830328000, 0.2831365000,
921 0.2832401000, 0.2833438000, 0.2834475000, 0.2835512000, 0.2836550000, 0.2837588000,
922 0.2838625000, 0.2839663000, 0.2840702000, 0.2841740000, 0.2842779000, 0.2843817000,
923 0.2844856000, 0.2845896000, 0.2846935000, 0.2847975000, 0.2849014000, 0.2850054000,
924 0.2851094000, 0.2852135000, 0.2853175000, 0.2854216000, 0.2855257000, 0.2856298000,
925 0.2857339000, 0.2858380000, 0.2859422000, 0.2860464000, 0.2861505000, 0.2862548000,
926 0.2863590000, 0.2864632000, 0.2865675000, 0.2866717000, 0.2867760000, 0.2868803000,
927 0.2869847000, 0.2870890000, 0.2871933000, 0.2872977000, 0.2874021000, 0.2875065000,
928 0.2876109000, 0.2877153000, 0.2878198000, 0.2879242000, 0.2880287000, 0.2881332000,
929 0.2882377000, 0.2883422000, 0.2884468000, 0.2885513000, 0.2886559000, 0.2887604000,
930 0.2888650000, 0.2889696000, 0.2890742000, 0.2891789000, 0.2892835000, 0.2893882000,
931 0.2894928000, 0.2895975000, 0.2897022000, 0.2898069000, 0.2899116000, 0.2900163000,
932 0.2901211000, 0.2902258000, 0.2903306000, 0.2904354000, 0.2905402000, 0.2906449000,
933 0.2907498000, 0.2908546000, 0.2909594000, 0.2910642000, 0.2911691000, 0.2912740000,
934 0.2913788000, 0.2914837000 };
935
936void EvtDToKmPipPipPi0::init() {
937 // check that there are 0 arguments
938 checkNArg( 0 );
939 checkNDaug( 4 );
940
941 checkSpinParent( EvtSpinType::SCALAR );
942 checkSpinDaughter( 0, EvtSpinType::SCALAR );
943 checkSpinDaughter( 1, EvtSpinType::SCALAR );
944 checkSpinDaughter( 2, EvtSpinType::SCALAR );
945 checkSpinDaughter( 3, EvtSpinType::SCALAR );
946
947 //------------------parameters-----------------
948 rho[0] = 1.00;
949 phi[0] = 0.00;
950 rho[1] = 1.6407e-01;
951 phi[1] = 1.4629e+00;
952 rho[2] = 5.3805e-01;
953 phi[2] = 5.6802e-02;
954 rho[3] = 9.8491e-01;
955 phi[3] = 4.5297e-01;
956 rho[4] = -7.5361e-01;
957 phi[4] = 3.9923e+00;
958 rho[5] = 7.2022e-01;
959 phi[5] = 3.6276e-01;
960 rho[6] = 1.5187e+00;
961 phi[6] = 1.1700e+00;
962 rho[7] = 3.0364e+00;
963 phi[7] = 2.9164e+00;
964 rho[8] = 5.2613e+00;
965 phi[8] = 5.0168e+00;
966 rho[9] = 1.1521e+01;
967 phi[9] = 3.8225e+00;
968 rho[10] = -7.3292e+00;
969 phi[10] = 4.8851e+00;
970 rho[11] = -8.9717e-01;
971 phi[11] = -2.7156e+00;
972 rho[12] = 7.9479e-01;
973 phi[12] = -1.9858e+00;
974 rho[13] = 4.6589e+00;
975 phi[13] = -8.5107e-01;
976 rho[14] = 6.8280e-01;
977 phi[14] = -2.1226e+00;
978 mD = 1.86486;
979 //------------------------new----------------------
980 mKst0 = 0.89555;
981 mrho = 0.77511;
982 mrho1450 = 1.465;
983 mK1270 = 1.289;
984 mK1400 = 1.403;
985 mK1460 = 1.482;
986 mK1650 = 1.672;
987 mK1680 = 1.718;
988
989 GKst0 = 0.0473;
990 Grho = 0.1491;
991 Grho1450 = 0.400;
992 GK1270 = 0.116;
993 GK1400 = 0.174;
994 GK1460 = 0.3356;
995 GK1650 = 0.158;
996 GK1680 = 0.322;
997
998 mass_Pion = 0.13957;
999 mass_Pion_N = 0.134977;
1000 mass_Eta = 0.547862;
1001 mass_Kaon = 0.493677;
1002 math_pi = 3.1415926;
1003
1004 rD2 = 5.0; // 5*5
1005 rRes1 = 3.0; // 3*3
1006 rRes2 = 3.0; // 3*3
1007
1008 GS1 = 0.636619783;
1009 GS2 = 0.01860182466;
1010 GS3 = 0.1591549458; // 1/(2*math_2pi)
1011 GS4 = 0.00620060822; // mass_Pion2/math_pi
1012
1013 int GG[4][4] = { { 1, 0, 0, 0 }, { 0, -1, 0, 0 }, { 0, 0, -1, 0 }, { 0, 0, 0, -1 } };
1014 int EE[4][4][4][4] = {
1015 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
1016 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, -1, 0 } },
1017 { { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 1, 0, 0 } },
1018 { { 0, 0, 0, 0 }, { 0, 0, 1, 0 }, { 0, -1, 0, 0 }, { 0, 0, 0, 0 } } },
1019 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 1, 0 } },
1020 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
1021 { { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 } },
1022 { { 0, 0, -1, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 } } },
1023 { { { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, -1, 0, 0 } },
1024 { { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 } },
1025 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
1026 { { 0, 1, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } },
1027 { { { 0, 0, 0, 0 }, { 0, 0, -1, 0 }, { 0, 1, 0, 0 }, { 0, 0, 0, 0 } },
1028 { { 0, 0, 1, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 } },
1029 { { 0, -1, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
1030 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } } };
1031 for ( int i = 0; i < 4; i++ )
1032 {
1033 for ( int j = 0; j < 4; j++ )
1034 {
1035 G[i][j] = GG[i][j];
1036 for ( int k = 0; k < 4; k++ )
1037 {
1038 for ( int l = 0; l < 4; l++ ) { E[i][j][k][l] = EE[i][j][k][l]; }
1039 }
1040 }
1041 }
1042}
1043
1044void EvtDToKmPipPipPi0::initProbMax() { setProbMax( 40000 ); }
1045
1046void EvtDToKmPipPipPi0::decay( EvtParticle* p ) {
1047 //--------------max
1048 /*
1049 double maxprob=0,maxq1270;
1050 for(int ir=0;ir<=10000000;ir++){
1051 p->initializePhaseSpace(getNDaug(),getDaugs());
1052 EvtVector4R km = p->getDaug(0)->getP4();
1053 EvtVector4R pip1 = p->getDaug(1)->getP4();
1054 EvtVector4R pip2 = p->getDaug(2)->getP4();
1055 EvtVector4R pi0 = p->getDaug(3)->getP4();
1056 int num = EvtPDL::getStdHep(p->getId());
1057
1058 double Pip1[4],Pip2[4],Km[4],Pi0[4];
1059 Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
1060 Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
1061 Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
1062 Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
1063
1064 if(num>0){
1065 Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
1066 Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
1067 Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
1068 Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
1069 }
1070 if(num<0){
1071 Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
1072 Km[1] = -km.get(1); Pip1[1] = -pip1.get(1); Pip2[1] = -pip2.get(1);Pi0[1] =
1073 -pi0.get(1); Km[2] = -km.get(2); Pip1[2] = -pip1.get(2); Pip2[2] = -pip2.get(2);Pi0[2] =
1074 -pi0.get(2); Km[3] = -km.get(3); Pip1[3] = -pip1.get(3); Pip2[3] = -pip2.get(3);Pi0[3] =
1075 -pi0.get(3);
1076 }
1077
1078 double _prob;
1079 double value;
1080 int nstates=15;
1081 int modetype[15]= {1,1,8,4,4,12,11,17,13,12,8,4,9,13,1};
1082 int g0[15] = {1,1,1,1,1,1,1,1,0,0,1,1,1,0,0};
1083 int g1[15] = {1,1,1,1,1,1,1,1,1,1,0,0,0,0,1};
1084 int g2[15] = {0,1,0,0,2,1,1,0,1,1,2,0,0,0,0};
1085 double mass1[15] = {mKst0,mKst0,mrho, mKst0, mKst0, mKst0, mKst0, mKst0, mrho,
1086 mKst0, mrho, mKst0, mKst0, mrho,mKst0}; double mass2[15] = {mrho, mrho,
1087 mK1270,mK1400,mK1400,mK1460,mK1680, mrho, mK1460,mK1460,mK1270,mK1400,mK1270, mK1460,mrho};
1088 double width1[15] = {GKst0,GKst0,Grho, GKst0, GKst0, GKst0, GKst0, GKst0, Grho, GKst0,
1089 Grho, GKst0, GKst0, Grho,GKst0}; double width2[15] = {Grho, Grho,
1090 GK1270,GK1400,GK1400,GK1460,GK1680, Grho, GK1460,GK1460,GK1270,GK1400,GK1270, GK1460,Grho};
1091 calEvaMy(Km,Pip1,Pip2,Pi0,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,nstates,_prob);
1092
1093 if(_prob>maxprob) {
1094 maxprob=_prob;
1095 }
1096 value = _prob;
1097 setProb(value);
1098 }
1099 printf("maxprob = %.10f\n", maxprob);
1100 return;
1101 */
1102 //--------------------main--------------
1103
1104 p->initializePhaseSpace( getNDaug(), getDaugs() );
1105 EvtVector4R km = p->getDaug( 0 )->getP4();
1106 EvtVector4R pip1 = p->getDaug( 1 )->getP4();
1107 EvtVector4R pip2 = p->getDaug( 2 )->getP4();
1108 EvtVector4R pi0 = p->getDaug( 3 )->getP4();
1109 int num = EvtPDL::getStdHep( p->getId() );
1110
1111 double Pip1[4], Pip2[4], Km[4], Pi0[4];
1112 //----------------------------------------------------------------------------------------
1113
1114 if ( num > 0 )
1115 {
1116 Km[0] = km.get( 0 );
1117 Pip1[0] = pip1.get( 0 );
1118 Pip2[0] = pip2.get( 0 );
1119 Pi0[0] = pi0.get( 0 );
1120 Km[1] = km.get( 1 );
1121 Pip1[1] = pip1.get( 1 );
1122 Pip2[1] = pip2.get( 1 );
1123 Pi0[1] = pi0.get( 1 );
1124 Km[2] = km.get( 2 );
1125 Pip1[2] = pip1.get( 2 );
1126 Pip2[2] = pip2.get( 2 );
1127 Pi0[2] = pi0.get( 2 );
1128 Km[3] = km.get( 3 );
1129 Pip1[3] = pip1.get( 3 );
1130 Pip2[3] = pip2.get( 3 );
1131 Pi0[3] = pi0.get( 3 );
1132 }
1133 if ( num < 0 )
1134 {
1135 Km[0] = km.get( 0 );
1136 Pip1[0] = pip1.get( 0 );
1137 Pip2[0] = pip2.get( 0 );
1138 Pi0[0] = pi0.get( 0 );
1139 Km[1] = -km.get( 1 );
1140 Pip1[1] = -pip1.get( 1 );
1141 Pip2[1] = -pip2.get( 1 );
1142 Pi0[1] = -pi0.get( 1 );
1143 Km[2] = -km.get( 2 );
1144 Pip1[2] = -pip1.get( 2 );
1145 Pip2[2] = -pip2.get( 2 );
1146 Pi0[2] = -pi0.get( 2 );
1147 Km[3] = -km.get( 3 );
1148 Pip1[3] = -pip1.get( 3 );
1149 Pip2[3] = -pip2.get( 3 );
1150 Pi0[3] = -pi0.get( 3 );
1151 }
1152
1153 double value;
1154 int nstates = 15;
1155 int modetype[15] = { 1, 1, 8, 4, 4, 12, 11, 17, 13, 12, 8, 4, 9, 13, 1 };
1156 int g0[15] = { 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0 };
1157 int g1[15] = { 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1 };
1158 int g2[15] = { 0, 1, 0, 0, 2, 1, 1, 0, 1, 1, 2, 0, 0, 0, 0 };
1159 double mass1[15] = { mKst0, mKst0, mrho, mKst0, mKst0, mKst0, mKst0, mKst0,
1160 mrho, mKst0, mrho, mKst0, mKst0, mrho, mKst0 };
1161 double mass2[15] = { mrho, mrho, mK1270, mK1400, mK1400, mK1460, mK1680, mrho,
1162 mK1460, mK1460, mK1270, mK1400, mK1270, mK1460, mrho };
1163 double width1[15] = { GKst0, GKst0, Grho, GKst0, GKst0, GKst0, GKst0, GKst0,
1164 Grho, GKst0, Grho, GKst0, GKst0, Grho, GKst0 };
1165 double width2[15] = { Grho, Grho, GK1270, GK1400, GK1400, GK1460, GK1680, Grho,
1166 GK1460, GK1460, GK1270, GK1400, GK1270, GK1460, Grho };
1167
1168 calEvaMy( Km, Pip1, Pip2, Pi0, mass1, mass2, width1, width2, rho, phi, g0, g1, g2, modetype,
1169 nstates, value );
1170 setProb( value );
1171
1172 return;
1173}
1174
1175void EvtDToKmPipPipPi0::KPiSLASS( double sa, double sb, double sc, double prop[2] ) {
1176 const double m1430 = 1.441;
1177 const double sa0 = 2.076481; // m1430*m1430;
1178 const double w1430 = 0.193;
1179 const double Lass1 = 0.25 / sa0;
1180 double tmp = sb - sc;
1181 double tmp1 = sa0 + tmp;
1182 double q0 = Lass1 * tmp1 * tmp1 - sb;
1183 if ( q0 < 0 ) q0 = 1e-16;
1184 double tmp2 = sa + tmp;
1185 double qs = 0.25 * tmp2 * tmp2 / sa - sb;
1186 double q = sqrt( qs );
1187 double width = w1430 * q * m1430 / sqrt( sa * q0 );
1188 double temp_R = atan( m1430 * width / ( sa0 - sa ) );
1189 if ( temp_R < 0 ) temp_R += math_pi;
1190 double deltaR = -1.915 + temp_R;
1191 double temp_F =
1192 atan( 0.226 * q / ( 2.0 - 3.8194 * qs ) ); // 2.0*1.07 = 2.14; 1.8*1.07 = 1.926
1193 if ( temp_F < 0 ) temp_F += math_pi;
1194 double deltaF = 0.002 + temp_F;
1195 double deltaS = deltaR + 2.0 * deltaF;
1196 double t1 = 0.96 * sin( deltaF );
1197 double t2 = sin( deltaR );
1198 double CF[2], CS[2];
1199 CF[0] = cos( deltaF );
1200 CF[1] = sin( deltaF );
1201 CS[0] = cos( deltaS );
1202 CS[1] = sin( deltaS );
1203 prop[0] = t1 * CF[0] + t2 * CS[0];
1204 prop[1] = t1 * CF[1] + t2 * CS[1];
1205}
1206
1207void EvtDToKmPipPipPi0::Com_Multi( double a1[2], double a2[2], double res[2] ) {
1208 res[0] = a1[0] * a2[0] - a1[1] * a2[1];
1209 res[1] = a1[1] * a2[0] + a1[0] * a2[1];
1210}
1211void EvtDToKmPipPipPi0::Com_Divide( double a1[2], double a2[2], double res[2] ) {
1212 res[0] = ( a1[0] * a2[0] + a1[1] * a2[1] ) / ( a2[0] * a2[0] + a2[1] * a2[1] );
1213 res[1] = ( a1[1] * a2[0] - a1[0] * a2[1] ) / ( a2[0] * a2[0] + a2[1] * a2[1] );
1214}
1215double EvtDToKmPipPipPi0::SCADot( double a1[4], double a2[4] ) {
1216 double _cal = a1[0] * a2[0] - a1[1] * a2[1] - a1[2] * a2[2] - a1[3] * a2[3];
1217 return _cal;
1218}
1219double EvtDToKmPipPipPi0::Barrier( int l, double sa, double sb, double sc, double r,
1220 double mass ) {
1221 double q = ( sa + sb - sc ) * ( sa + sb - sc ) / ( 4 * sa ) - sb;
1222 // if(q < 0) q = 1e-16;
1223 if ( q < 0 ) q = -q;
1224
1225 double z;
1226 z = q * r * r;
1227 double sa0;
1228 sa0 = mass * mass;
1229 double q0 = ( sa0 + sb - sc ) * ( sa0 + sb - sc ) / ( 4 * sa0 ) - sb;
1230 // if(q0 < 0) q0 = 1e-16;
1231 if ( q0 < 0 ) q0 = -q0;
1232 double z0 = q0 * r * r;
1233 double F = 0.0;
1234 if ( l == 0 ) F = 1;
1235 if ( l == 1 ) F = sqrt( ( 1 + z0 ) / ( 1 + z ) );
1236 if ( l == 2 ) F = sqrt( ( 9 + 3 * z0 + z0 * z0 ) / ( 9 + 3 * z + z * z ) );
1237 // if(l == 1) F = sqrt((2*r)/(1+z));
1238 // if(l == 2) F = sqrt((13*r*r)/(9+3*z+z*z));
1239 return F;
1240}
1241void EvtDToKmPipPipPi0::calt1( double daug1[4], double daug0[4], double t1[4] ) {
1242 double p, pq;
1243 double pa[4], qa[4];
1244 for ( int i = 0; i < 4; i++ )
1245 {
1246 pa[i] = daug1[i] + daug0[i];
1247 qa[i] = daug1[i] - daug0[i];
1248 }
1249 p = SCADot( pa, pa );
1250 pq = SCADot( pa, qa );
1251 for ( int i = 0; i < 4; i++ ) { t1[i] = qa[i] - pq / p * pa[i]; }
1252}
1253void EvtDToKmPipPipPi0::calt2( double daug1[4], double daug0[4], double t2[4][4] ) {
1254 double p, r;
1255 double pa[4], t1[4];
1256 calt1( daug1, daug0, t1 );
1257 r = SCADot( t1, t1 ) / 3.0;
1258 for ( int i = 0; i < 4; i++ ) { pa[i] = daug1[i] + daug0[i]; }
1259 p = SCADot( pa, pa );
1260 for ( int i = 0; i < 4; i++ )
1261 {
1262 for ( int j = 0; j < 4; j++ )
1263 { t2[i][j] = t1[i] * t1[j] - r * ( G[i][j] - pa[i] * pa[j] / p ); }
1264 }
1265}
1266
1267double EvtDToKmPipPipPi0::wid( double mass2, double mass, double sa, double sb, double sc,
1268 double r, int l ) {
1269 double widm = 0.;
1270 double m = sqrt( sa );
1271 double tmp = sb - sc;
1272 double tmp1 = sa + tmp;
1273 double q = 0.25 * tmp1 * tmp1 / sa - sb;
1274 // if(q<0) q = 1e-16;
1275 if ( q < 0 ) q = -q;
1276 double tmp2 = mass2 + tmp;
1277 double q0 = 0.25 * tmp2 * tmp2 / mass2 - sb;
1278 // if(q0<0) q0 = 1e-16;
1279 if ( q0 < 0 ) q0 = -q0;
1280 double z = q * r * r;
1281 double z0 = q0 * r * r;
1282 double t = q / q0;
1283 if ( l == 0 ) { widm = sqrt( t ) * mass / m; }
1284 else if ( l == 1 ) { widm = t * sqrt( t ) * mass / m * ( 1 + z0 ) / ( 1 + z ); }
1285 else if ( l == 2 )
1286 { widm = t * t * sqrt( t ) * mass / m * ( 9 + 3 * z0 + z0 * z0 ) / ( 9 + 3 * z + z * z ); }
1287 return widm;
1288}
1289double EvtDToKmPipPipPi0::widl1( double mass2, double mass, double sa, double sb, double sc,
1290 double r ) {
1291 double widm = 0.;
1292 double m = sqrt( sa );
1293 double tmp = sb - sc;
1294 double tmp1 = sa + tmp;
1295 double q = 0.25 * tmp1 * tmp1 / sa - sb;
1296 // if(q<0) q = 1e-16;
1297 if ( q < 0 ) q = -q;
1298 double tmp2 = mass2 + tmp;
1299 double q0 = 0.25 * tmp2 * tmp2 / mass2 - sb;
1300 // if(q0<0) q0 = 1e-16;
1301 if ( q0 < 0 ) q0 = -q0;
1302 double z = q * r * r;
1303 double z0 = q0 * r * r;
1304 double F = ( 1 + z0 ) / ( 1 + z );
1305 double t = q / q0;
1306 widm = t * sqrt( t ) * mass / m * F;
1307 return widm;
1308}
1309void EvtDToKmPipPipPi0::propagatorNBW( double mass2, double mass, double width, double sa,
1310 double sb, double sc, double r, int l,
1311 double prop[2] ) {
1312 double a[2], b[2];
1313 a[0] = 1;
1314 a[1] = 0;
1315 b[0] = mass2 - sa;
1316 b[1] = -mass * width;
1317 Com_Divide( a, b, prop );
1318}
1319
1320void EvtDToKmPipPipPi0::propagatorRBW( double mass2, double mass, double width, double sa,
1321 double sb, double sc, double r, int l,
1322 double prop[2] ) {
1323 double a[2], b[2];
1324 a[0] = 1;
1325 a[1] = 0;
1326 b[0] = mass2 - sa;
1327 b[1] = -mass * width * wid( mass2, mass, sa, sb, sc, r, l );
1328 Com_Divide( a, b, prop );
1329}
1330void EvtDToKmPipPipPi0::propagatorRBW_K1( double mass2, double mass, double width, double sa,
1331 double sb, double sc, double r2, int l,
1332 double prop[2] ) {
1333 double a[2], b[2];
1334 int iii = int( sa * 1000 / 2 ) - 300; //(4000-0)/2000=2
1335 if ( iii < 0 ) { iii = 0; }
1336 double width_K1;
1337 if ( mass < 1.35 ) { width_K1 = K1270_width[iii]; }
1338 else { width_K1 = K1400_width[iii]; }
1339 // printf("width_K1: %.8f\n", width_K1);
1340 a[0] = 1;
1341 a[1] = 0;
1342 b[0] = mass2 - sa;
1343 // b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l);
1344 b[1] = -mass * width_K1;
1345 Com_Divide( a, b, prop );
1346}
1347void EvtDToKmPipPipPi0::propagatorRBWl1( double mass2, double mass, double width, double sa,
1348 double sb, double sc, double r, double prop[2] )
1349// firstly code for propagator_omg, now could use for propagatorRBW that l=1, and will be more
1350// quick.
1351{
1352 double a[2], b[2];
1353 a[0] = 1;
1354 a[1] = 0;
1355 b[0] = mass2 - sa;
1356 b[1] = -mass * width * widl1( mass2, mass, sa, sb, sc, r );
1357 Com_Divide( a, b, prop );
1358}
1359//------------GS---used by rho----------------------------
1360void EvtDToKmPipPipPi0::propagatorGS( double mass2, double mass, double width, double sa,
1361 double sb, double sc, double r, double prop[2] ) {
1362 double a[2], b[2];
1363 double tmp = sb - sc;
1364 double tmp1 = sa + tmp;
1365 double q2 = 0.25 * tmp1 * tmp1 / sa - sb;
1366 // if(q2<0) q2 = 1e-16;
1367 if ( q2 < 0 ) q2 = -q2;
1368
1369 double tmp2 = mass2 + tmp;
1370 double q02 = 0.25 * tmp2 * tmp2 / mass2 - sb;
1371 // if(q02<0) q02 = 1e-16;
1372 if ( q02 < 0 ) q02 = -q02;
1373
1374 double q = sqrt( q2 );
1375 double q0 = sqrt( q02 );
1376 double m = sqrt( sa );
1377 double q03 = q0 * q02;
1378 double tmp3 = log( mass + 2 * q0 ) + 1.2760418309; // log(mass_2Pion) = 1.2760418309;
1379
1380 double h = GS1 * q / m * ( log( m + 2 * q ) + 1.2760418309 );
1381 double h0 = GS1 * q0 / mass * tmp3;
1382 double dh = h0 * ( 0.125 / q02 - 0.5 / mass2 ) + GS3 / mass2;
1383 double d = GS2 / q02 * tmp3 + GS3 * mass / q0 - GS4 * mass / q03;
1384 double f = mass2 / q03 * ( q2 * ( h - h0 ) + ( mass2 - sa ) * q02 * dh );
1385
1386 a[0] = 1.0 + d * width / mass;
1387 a[1] = 0.0;
1388 b[0] = mass2 - sa + width * f;
1389 b[1] = -mass * width * widl1( mass2, mass, sa, sb, sc, r );
1390 Com_Divide( a, b, prop );
1391}
1392void EvtDToKmPipPipPi0::calEvaMy( double* Km, double* Pip1, double* Pip2, double* Pi0,
1393 double* mass1, double* mass2, double* width1, double* width2,
1394 double* amp, double* phase, int* g0, int* g1, int* g2,
1395 int* modetype, int nstates, double& Result )
1396
1397{
1398
1399 double pKstr1[4], pKstr2[4], pKstr3[4], prho1[4], prho2[4], pK11[4], pK12[4], pK13[4],
1400 pA1[4], pD[4];
1401 double Amp_KPiS1[2], Amp_PiPiS1[2], Amp_KPiS2[2], Amp_PiPiS2[2];
1402
1403 for ( int i = 0; i != 4; i++ )
1404 {
1405 pKstr1[i] = Km[i] + Pip1[i];
1406 pKstr2[i] = Km[i] + Pip2[i];
1407 pKstr3[i] = Km[i] + Pi0[i];
1408 prho1[i] = Pip1[i] + Pi0[i];
1409 prho2[i] = Pip2[i] + Pi0[i];
1410 pK11[i] = Km[i] + prho1[i];
1411 pK12[i] = Km[i] + prho2[i];
1412 pK13[i] = Km[i] + Pip1[i] + Pip2[i];
1413 pD[i] = Km[i] + prho1[i] + Pip2[i];
1414 }
1415 double skstr1, skstr2;
1416 double skstr3, srho1, srho2, sk11, sk12, sk13, sD;
1417 double skaon, spion1, spion2, spi0;
1418 skaon = SCADot( Km, Km );
1419 spion1 = SCADot( Pip1, Pip1 );
1420 spion2 = SCADot( Pip2, Pip2 );
1421 spi0 = SCADot( Pi0, Pi0 );
1422 skstr1 = SCADot( pKstr1, pKstr1 );
1423 skstr2 = SCADot( pKstr2, pKstr2 );
1424 skstr3 = SCADot( pKstr3, pKstr3 );
1425 srho1 = SCADot( prho1, prho1 );
1426 srho2 = SCADot( prho2, prho2 );
1427 sk11 = SCADot( pK11, pK11 );
1428 sk12 = SCADot( pK12, pK12 );
1429 sk13 = SCADot( pK13, pK13 );
1430 sD = SCADot( pD, pD );
1431 double t1rho1[4], t1rho2[4], t1kstr3[4];
1432 //------------------------
1433 double t1kstr1[4], t1kstr2[4], t2k21[4][4], t2k22[4][4], t2k31[4][4], t2k32[4][4],
1434 t2k11[4][4], t2k12[4][4], t2k13[4][4], t2k14[4][4];
1435 double t2rho1[4][4], t2rho2[4][4], t2kstr3[4][4], t2kstr1[4][4], t2kstr2[4][4];
1436 double t1K11[4], t1K12[4], t1K13[4], t1K14[4];
1437 double t2kk1[4][4], t2kk2[4][4];
1438
1439 //--------------------------
1440 calt1( Pip1, Pi0, t1rho1 );
1441 calt1( Pip2, Pi0, t1rho2 );
1442 calt1( Km, Pip1, t1kstr1 );
1443 calt1( Km, Pip2, t1kstr2 );
1444 calt1( Km, Pi0, t1kstr3 );
1445 calt1( pKstr1, Pi0, t1K11 );
1446 calt1( pKstr2, Pi0, t1K12 );
1447 calt1( pKstr3, Pip1, t1K13 );
1448 calt1( pKstr3, Pip2, t1K14 );
1449 calt2( Pip2, Pi0, t2rho1 );
1450 calt2( Pip2, Pi0, t2rho2 );
1451 calt2( Km, Pi0, t2kstr3 );
1452 calt2( Km, Pip1, t2kstr1 );
1453 calt2( Km, Pip2, t2kstr2 );
1454 calt2( pKstr3, Pip1, t2k11 );
1455 calt2( pKstr3, Pip2, t2k12 );
1456 calt2( prho1, Km, t2k21 );
1457 calt2( prho2, Km, t2k22 );
1458 calt2( pKstr1, Pi0, t2k31 );
1459 calt2( pKstr2, Pi0, t2k32 );
1460 calt2( pKstr1, Pip2, t2kk1 );
1461 calt2( pKstr2, Pip1, t2kk2 );
1462 //----------------------------------
1463 double cof[2], amp_tmp1[2], amp_PDF[2], PDF[2];
1464 double amp_tmp2[2];
1465 double amp_tmp[2];
1466
1467 PDF[0] = 0.0;
1468 PDF[1] = 0.0;
1469
1470 double temp_PDF1, temp_PDF2, tt1, tt2, tmp1, tmp2;
1471
1472 double temp_PDF3, temp_PDF4;
1473
1474 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];
1475
1476 double t1D[4], t1D_1[4], t1D_2[4];
1477
1478 double t1K1_1[4], t1K1_2[4], t1K1_3[4], t1K1_4[4];
1479
1480 double B1_rho1 = -1.0, B1_rho2 = -1.0;
1481 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;
1482 double B1_K1_rho1 = -1.0, B1_K1_rho2 = -1.0, B2_K1_rho1 = -1.0, B2_K1_rho2 = -1.0;
1483 double B1_K1_Kstr3_Pi1 = -1.0, B1_K1_Kstr3_Pi2 = -1.0, B2_K1_Kstr3_Pi1 = -1.0,
1484 B2_K1_Kstr3_Pi2 = -1.0;
1485 double B1_K1_Kstr1_Pi0 = -1.0, B1_K1_Kstr2_Pi0 = -1.0;
1486 double B1_kstr3 = -1.0;
1487 double B2_K1p_kstr1 = -1.0, B2_K1p_kstr2 = -1.0;
1488 double B1_K1p_kstr1 = -1.0, B1_K1p_kstr2 = -1.0;
1489 double B1_kstr1_Pi1 = -1.0, B1_kstr2_Pi2 = -1.0, B1_D_Kstr2_rho1 = -1.0,
1490 B1_D_Kstr1_rho2 = -1.0, B1_D_K1p_Pi0 = -1.0;
1491 double B2_K1_Kstr1 = -1.0, B2_K1_Kstr2 = -1.0, B2_D_Kstr1_rho2 = -1.0,
1492 B2_D_Kstr2_rho1 = -1.0;
1493 double B2_rho1 = -1.0, B2_rho2 = -1.0;
1494 double B2_kstr1_Pi1 = -1.0, B2_kstr2_Pi2 = -1.0, B2_kstr3_Pi0 = -1.0;
1495 ;
1496 double mass1sq, mass2sq;
1497 int isKst = 0, isRho = 0, isPiPi_S = 0, isKPi_S = 0;
1498 double proRho1[2], proRho2[2], proPiPi1_S[2], proPiPi2_S[2], proKPi_S[2];
1499
1500 double tmp3, tmp4;
1501 double t2D[4][4], t2D_1[4][4], t2D_2[4][4];
1502 double pro2_1[2], pro2_2[2], pro3_1[2], pro3_2[2], pro4_1[2], pro4_2[2];
1503 double t2f2_1[4][4], t2f2_2[4][4];
1504 double t2G2_1[4][4], t2G2_2[4][4], t2D11_1[4][4], t2D11_2[4][4];
1505
1506 //-------------------------------------------------------------------------
1507 // for(int i=0; i<nstates; i++){
1508 for ( int i = 0; i < 15; i++ )
1509 {
1510 // for(int i=2; i<15; i++){
1511 cof[0] = amp[i] * cos( phase[i] );
1512 cof[1] = amp[i] * sin( phase[i] );
1513 mass1sq = mass1[i] * mass1[i];
1514 mass2sq = mass2[i] * mass2[i];
1515 temp_PDF1 = 0;
1516 temp_PDF2 = 0;
1517 temp_PDF3 = 0;
1518 temp_PDF4 = 0;
1519
1520 if ( modetype[i] == 1 )
1521 {
1522 B1_kstr1_Pi1 = Barrier( 1, skstr1, skaon, spion1, rRes2, mass1[i] );
1523 B1_kstr2_Pi2 = Barrier( 1, skstr2, skaon, spion2, rRes2, mass1[i] );
1524 B1_rho1 = Barrier( 1, srho1, spion1, spi0, rRes2, mass2[i] );
1525 B1_rho2 = Barrier( 1, srho2, spion2, spi0, rRes2, mass2[i] );
1526 if ( g2[i] == 0 )
1527 {
1528 for ( int w = 0; w < 4; w++ )
1529 {
1530 temp_PDF1 += G[w][w] * t1kstr1[w] * t1rho2[w];
1531 temp_PDF2 += G[w][w] * t1kstr2[w] * t1rho1[w];
1532 }
1533 tmp1 = B1_kstr1_Pi1 * B1_rho2 * temp_PDF1;
1534 tmp2 = B1_kstr2_Pi2 * B1_rho1 * temp_PDF2;
1535 }
1536 else if ( g2[i] == 1 )
1537 {
1538 calt1( pKstr1, prho2, t1D_1 );
1539 calt1( pKstr2, prho1, t1D_2 );
1540 for ( int w = 0; w < 4; w++ )
1541 {
1542 for ( int j = 0; j < 4; j++ )
1543 {
1544 for ( int k = 0; k < 4; k++ )
1545 {
1546 for ( int l = 0; l < 4; l++ )
1547 {
1548 temp_PDF1 += E[w][j][k][l] * pD[w] * t1D_1[j] * t1kstr1[k] * t1rho2[l] *
1549 G[w][w] * G[j][j] * G[k][k] * G[l][l];
1550 temp_PDF2 += E[w][j][k][l] * pD[w] * t1D_2[j] * t1kstr2[k] * t1rho1[l] *
1551 G[w][w] * G[j][j] * G[k][k] * G[l][l];
1552 }
1553 }
1554 }
1555 }
1556 B1_D_Kstr1_rho2 = Barrier( 1, sD, skstr1, srho2, rD2, mD );
1557 B1_D_Kstr2_rho1 = Barrier( 1, sD, skstr2, srho1, rD2, mD );
1558 tmp1 = B1_kstr1_Pi1 * B1_rho2 * B1_D_Kstr1_rho2 * temp_PDF1;
1559 tmp2 = B1_kstr2_Pi2 * B1_rho1 * B1_D_Kstr2_rho1 * temp_PDF2;
1560 }
1561 else if ( g2[i] == 2 )
1562 {
1563 calt2( pKstr1, prho2, t2D_1 );
1564 calt2( pKstr2, prho1, t2D_2 );
1565 for ( int w = 0; w < 4; w++ )
1566 {
1567 for ( int j = 0; j < 4; j++ )
1568 {
1569 temp_PDF1 += t2D_1[w][j] * t1kstr1[w] * t1rho2[j] * G[w][w] * G[j][j];
1570 temp_PDF2 += t2D_2[w][j] * t1kstr2[w] * t1rho1[j] * G[w][w] * G[j][j];
1571 }
1572 }
1573 B2_D_Kstr2_rho1 = Barrier( 2, sD, skstr1, srho2, rD2, mD );
1574 B2_D_Kstr2_rho1 = Barrier( 2, sD, skstr2, srho1, rD2, mD );
1575 tmp1 = B1_kstr1_Pi1 * B1_rho2 * B2_D_Kstr1_rho2 * temp_PDF1;
1576 tmp2 = B1_kstr2_Pi2 * B1_rho1 * B2_D_Kstr2_rho1 * temp_PDF2;
1577 }
1578 if ( g0[i] == 1 )
1579 {
1580 propagatorRBW( mass1sq, mass1[i], width1[i], skstr1, skaon, spion1, rRes2, 1, pro0_1 );
1581 propagatorRBW( mass1sq, mass1[i], width1[i], skstr2, skaon, spion2, rRes2, 1, pro0_2 );
1582 }
1583 else if ( g0[i] == 0 )
1584 {
1585 pro0_1[0] = 1;
1586 pro0_1[1] = 0;
1587 pro0_2[0] = 1;
1588 pro0_2[1] = 0;
1589 }
1590 if ( g1[i] == 1 )
1591 {
1592 propagatorGS( mass2sq, mass2[i], width2[i], srho2, spion2, spi0, rRes2, pro1_1 );
1593 propagatorGS( mass2sq, mass2[i], width2[i], srho1, spion1, spi0, rRes2, pro1_2 );
1594 }
1595 else if ( g1[i] == 0 )
1596 {
1597 pro1_1[0] = 1;
1598 pro1_1[1] = 0;
1599 pro1_2[0] = 1;
1600 pro1_2[1] = 0;
1601 }
1602 Com_Multi( pro0_1, pro1_1, pro_1 );
1603 Com_Multi( pro0_2, pro1_2, pro_2 );
1604 amp_tmp1[0] = tmp1 * pro_1[0];
1605 amp_tmp1[1] = tmp1 * pro_1[1];
1606 amp_tmp2[0] = tmp2 * pro_2[0];
1607 amp_tmp2[1] = tmp2 * pro_2[1];
1608 }
1609 else if ( modetype[i] == 4 )
1610 {
1611 B1_D_K1_Pi2 = Barrier( 1, sD, sk11, spion2, rD2, mass2[i] );
1612 B1_D_K1_Pi1 = Barrier( 1, sD, sk12, spion1, rD2, mass2[i] );
1613 B1_kstr3 = Barrier( 1, skstr3, skaon, spi0, rRes2, mass1[i] );
1614 B1_kstr1_Pi1 = Barrier( 1, skstr1, skaon, spion1, rRes2, mass1[i] );
1615 B1_kstr2_Pi2 = Barrier( 1, skstr2, skaon, spion2, rRes2, mass1[i] );
1616 calt1( pK11, Pip2, t1D_1 );
1617 calt1( pK12, Pip1, t1D_2 );
1618 if ( g2[i] == 0 )
1619 {
1620 for ( int w = 0; w < 4; w++ )
1621 {
1622 for ( int j = 0; j < 4; j++ )
1623 {
1624 tt2 = t1kstr3[j] * G[w][w] * G[j][j];
1625 temp_PDF1 += t1D_1[w] * ( G[w][j] - pK11[w] * pK11[j] / sk11 ) * t1kstr1[j] *
1626 G[w][w] * G[j][j];
1627 temp_PDF2 += t1D_2[w] * ( G[w][j] - pK12[w] * pK12[j] / sk12 ) * t1kstr2[j] *
1628 G[w][w] * G[j][j];
1629 temp_PDF3 += t1D_1[w] * ( G[w][j] - pK11[w] * pK11[j] / sk11 ) * tt2;
1630 temp_PDF4 += t1D_2[w] * ( G[w][j] - pK12[w] * pK12[j] / sk12 ) * tt2;
1631 }
1632 }
1633 tmp1 = B1_kstr1_Pi1 * B1_D_K1_Pi2 * temp_PDF1;
1634 tmp2 = B1_kstr2_Pi2 * B1_D_K1_Pi1 * temp_PDF2;
1635 tmp3 = B1_kstr3 * B1_D_K1_Pi2 * temp_PDF3;
1636 tmp4 = B1_kstr3 * B1_D_K1_Pi1 * temp_PDF4;
1637 }
1638 else if ( g2[i] == 2 )
1639 {
1640 for ( int w = 0; w < 4; w++ )
1641 {
1642 for ( int j = 0; j < 4; j++ )
1643 {
1644 tt2 = t1kstr3[j] * G[w][w] * G[j][j];
1645 temp_PDF1 += t1D_1[w] * t2k31[w][j] * t1kstr1[j] * G[w][w] * G[j][j];
1646 temp_PDF2 += t1D_2[w] * t2k32[w][j] * t1kstr2[j] * G[w][w] * G[j][j];
1647 temp_PDF3 += t1D_1[w] * t2k11[w][j] * tt2;
1648 temp_PDF4 += t1D_2[w] * t2k12[w][j] * tt2;
1649 }
1650 }
1651 B2_K1_Kstr3_Pi1 = Barrier( 2, sk11, skstr3, spion1, rRes2, mass2[i] );
1652 B2_K1_Kstr3_Pi2 = Barrier( 2, sk12, skstr3, spion2, rRes2, mass2[i] );
1653 B2_K1_Kstr1 = Barrier( 2, sk11, skstr1, spi0, rRes2, mass2[i] );
1654 B2_K1_Kstr2 = Barrier( 2, sk12, skstr2, spi0, rRes2, mass2[i] );
1655 tmp1 = B1_kstr1_Pi1 * B2_K1_Kstr1 * B1_D_K1_Pi2 * temp_PDF1;
1656 tmp2 = B1_kstr2_Pi2 * B2_K1_Kstr2 * B1_D_K1_Pi1 * temp_PDF2;
1657 tmp3 = B1_kstr3 * B2_K1_Kstr3_Pi1 * B1_D_K1_Pi2 * temp_PDF3;
1658 tmp4 = B1_kstr3 * B2_K1_Kstr3_Pi2 * B1_D_K1_Pi1 * temp_PDF4;
1659 }
1660 propagatorRBW( mass1sq, mass1[i], width1[i], skstr3, skaon, spi0, rRes2, 1, pro2_1 );
1661 propagatorRBW( mass1sq, mass1[i], width1[i], skstr3, skaon, spi0, rRes2, 1, pro2_2 );
1662 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk11, skstr3, spion1, rRes2, g2[i],
1663 pro3_1 );
1664 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk12, skstr3, spion2, rRes2, g2[i],
1665 pro3_2 );
1666 if ( g0[i] == 1 )
1667 {
1668 propagatorRBW( mass1sq, mass1[i], width1[i], skstr1, skaon, spion1, rRes2, 1, pro0_1 );
1669 propagatorRBW( mass1sq, mass1[i], width1[i], skstr2, skaon, spion2, rRes2, 1, pro0_2 );
1670 }
1671 else if ( g0[i] == 0 )
1672 {
1673 pro0_1[0] = 1;
1674 pro0_1[1] = 0;
1675 pro0_2[0] = 1;
1676 pro0_2[1] = 0;
1677 pro2_1[0] = 1;
1678 pro2_1[1] = 0;
1679 pro2_2[0] = 1;
1680 pro2_2[1] = 0;
1681 }
1682 if ( g1[i] == 1 )
1683 {
1684 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk11, skstr1, spi0, rRes2, g2[i],
1685 pro1_1 );
1686 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk12, skstr2, spi0, rRes2, g2[i],
1687 pro1_2 );
1688 }
1689 else if ( g1[i] == 0 )
1690 {
1691 pro1_1[0] = 1;
1692 pro1_1[1] = 0;
1693 pro1_2[0] = 1;
1694 pro1_2[1] = 0;
1695 pro3_1[0] = 1;
1696 pro3_1[1] = 0;
1697 pro3_2[0] = 1;
1698 pro3_2[1] = 0;
1699 }
1700 Com_Multi( pro0_1, pro1_1, pro_1 );
1701 Com_Multi( pro2_1, pro3_1, pro4_1 );
1702 Com_Multi( pro0_2, pro1_2, pro_2 );
1703 Com_Multi( pro2_2, pro3_2, pro4_2 );
1704 amp_tmp1[0] = tmp1 * pro_1[0] - tmp3 * pro4_1[0];
1705 amp_tmp1[1] = tmp1 * pro_1[1] - tmp3 * pro4_1[1];
1706 amp_tmp2[0] = tmp2 * pro_2[0] - tmp4 * pro4_2[0];
1707 amp_tmp2[1] = tmp2 * pro_2[1] - tmp4 * pro4_2[1];
1708 }
1709 else if ( modetype[i] == 8 )
1710 {
1711 B1_rho1 = Barrier( 1, srho1, spion1, spi0, rRes2, mass1[i] );
1712 B1_rho2 = Barrier( 1, srho2, spion2, spi0, rRes2, mass1[i] );
1713 B1_D_K1_Pi2 = Barrier( 1, sD, sk11, spion2, rD2, mD );
1714 B1_D_K1_Pi1 = Barrier( 1, sD, sk12, spion1, rD2, mD );
1715 calt1( pK11, Pip2, t1D_1 );
1716 calt1( pK12, Pip1, t1D_2 );
1717 if ( g2[i] == 0 )
1718 {
1719 for ( int w = 0; w < 4; w++ )
1720 {
1721 for ( int j = 0; j < 4; j++ )
1722 {
1723 temp_PDF1 += t1D_1[w] * ( G[w][j] - pK11[w] * pK11[j] / sk11 ) * t1rho1[j] *
1724 G[w][w] * G[j][j];
1725 temp_PDF2 += t1D_2[w] * ( G[w][j] - pK12[w] * pK12[j] / sk12 ) * t1rho2[j] *
1726 G[w][w] * G[j][j];
1727 }
1728 }
1729 tmp1 = B1_rho1 * B1_D_K1_Pi2 * temp_PDF1;
1730 tmp2 = B1_rho2 * B1_D_K1_Pi1 * temp_PDF2;
1731 }
1732 else if ( g2[i] == 2 )
1733 {
1734 for ( int w = 0; w < 4; w++ )
1735 {
1736 for ( int j = 0; j < 4; j++ )
1737 {
1738 temp_PDF1 += t1D_1[w] * t2k21[w][j] * t1rho1[j] * G[w][w] * G[j][j];
1739 temp_PDF2 += t1D_2[w] * t2k22[w][j] * t1rho2[j] * G[w][w] * G[j][j];
1740 }
1741 }
1742 if ( B2_K1_rho1 < 0.0 ) B2_K1_rho1 = Barrier( 2, sk11, srho1, skaon, rRes2, mass2[i] );
1743 if ( B2_K1_rho2 < 0.0 ) B2_K1_rho2 = Barrier( 2, sk12, srho2, skaon, rRes2, mass2[i] );
1744 tmp1 = B1_rho1 * B2_K1_rho1 * B1_D_K1_Pi2 * temp_PDF1;
1745 tmp2 = B1_rho2 * B2_K1_rho2 * B1_D_K1_Pi1 * temp_PDF2;
1746 }
1747 if ( g0[i] == 1 )
1748 {
1749 propagatorGS( mass1sq, mass1[i], width1[i], srho1, spion1, spi0, rRes2, pro0_1 );
1750 propagatorGS( mass1sq, mass1[i], width1[i], srho2, spion2, spi0, rRes2, pro0_2 );
1751 }
1752 else if ( g0[i] == 0 )
1753 {
1754 pro0_1[0] = 1;
1755 pro0_1[1] = 0;
1756 pro0_2[0] = 1;
1757 pro0_2[1] = 0;
1758 }
1759 if ( g1[i] == 1 )
1760 {
1761 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk11, srho1, skaon, rRes2, g2[i],
1762 pro1_1 );
1763 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk12, srho2, skaon, rRes2, g2[i],
1764 pro1_2 );
1765 }
1766 else if ( g1[i] == 0 )
1767 {
1768 pro1_1[0] = 1;
1769 pro1_1[1] = 0;
1770 pro1_2[0] = 1;
1771 pro1_2[1] = 0;
1772 }
1773
1774 Com_Multi( pro0_1, pro1_1, pro_1 );
1775 Com_Multi( pro0_2, pro1_2, pro_2 );
1776 amp_tmp1[0] = tmp1 * pro_1[0];
1777 amp_tmp1[1] = tmp1 * pro_1[1];
1778 amp_tmp2[0] = tmp2 * pro_2[0];
1779 amp_tmp2[1] = tmp2 * pro_2[1];
1780 }
1781 else if ( modetype[i] == 11 )
1782 {
1783 B1_kstr3 = Barrier( 1, skstr3, skaon, spi0, rRes2, mass1[i] );
1784 B1_kstr1_Pi1 = Barrier( 1, skstr1, skaon, spion1, rRes2, mass1[i] );
1785 B1_kstr2_Pi2 = Barrier( 1, skstr2, skaon, spion2, rRes2, mass1[i] );
1786 B1_D_K1_Pi2 = Barrier( 1, sD, sk11, spion2, rD2, mD );
1787 B1_D_K1_Pi1 = Barrier( 1, sD, sk12, spion1, rD2, mD );
1788 for ( int w = 0; w < 4; w++ )
1789 {
1790 for ( int j = 0; j < 4; j++ )
1791 {
1792 for ( int k = 0; k < 4; k++ )
1793 {
1794 for ( int l = 0; l < 4; l++ )
1795 {
1796 temp_PDF1 += E[w][j][k][l] * pK11[w] * G[w][w] * ( pKstr1[j] - Pi0[j] ) *
1797 G[j][j] * Pip2[k] * G[k][k] * ( Km[l] - Pip1[l] ) * G[l][l];
1798 temp_PDF2 += E[w][j][k][l] * pK12[w] * G[w][w] * ( pKstr2[j] - Pi0[j] ) *
1799 G[j][j] * Pip1[k] * G[k][k] * ( Km[l] - Pip2[l] ) * G[l][l];
1800 temp_PDF3 += E[w][j][k][l] * pK11[w] * G[w][w] * ( pKstr3[j] - Pip1[j] ) *
1801 G[j][j] * Pip2[k] * G[k][k] * ( Km[l] - Pi0[l] ) * G[l][l];
1802 temp_PDF4 += E[w][j][k][l] * pK12[w] * G[w][w] * ( pKstr3[j] - Pip2[j] ) *
1803 G[j][j] * Pip1[k] * G[k][k] * ( Km[l] - Pi0[l] ) * G[l][l];
1804 }
1805 }
1806 }
1807 }
1808 B1_K1_Kstr1_Pi0 = Barrier( 1, sk11, skstr1, spi0, rRes2, mass2[i] );
1809 B1_K1_Kstr2_Pi0 = Barrier( 1, sk12, skstr2, spi0, rRes2, mass2[i] );
1810 B1_K1_Kstr3_Pi1 = Barrier( 1, sk11, skstr3, spion1, rRes2, mass2[i] );
1811 B1_K1_Kstr3_Pi1 = Barrier( 1, sk12, skstr3, spion2, rRes2, mass2[i] );
1812 tmp1 = B1_kstr1_Pi1 * B1_K1_Kstr1_Pi0 * B1_D_K1_Pi2 * temp_PDF1;
1813 tmp2 = B1_kstr2_Pi2 * B1_K1_Kstr2_Pi0 * B1_D_K1_Pi1 * temp_PDF2;
1814 tmp3 = B1_kstr3 * B1_K1_Kstr3_Pi1 * B1_D_K1_Pi2 * temp_PDF3;
1815 tmp4 = B1_kstr3 * B1_K1_Kstr3_Pi2 * B1_D_K1_Pi1 * temp_PDF4;
1816 propagatorRBW( mass1sq, mass1[i], width1[i], skstr3, skaon, spi0, rRes2, 1, pro2_1 );
1817 propagatorRBW( mass1sq, mass1[i], width1[i], skstr3, skaon, spi0, rRes2, 1, pro2_2 );
1818 propagatorRBW( mass2sq, mass2[i], width2[i], sk11, skstr3, spion1, rRes2, 1, pro3_1 );
1819 propagatorRBW( mass2sq, mass2[i], width2[i], sk12, skstr3, spion2, rRes2, 1, pro3_2 );
1820 if ( g0[i] == 1 )
1821 {
1822 propagatorRBW( mass1sq, mass1[i], width1[i], skstr1, skaon, spion1, rRes2, 1, pro0_1 );
1823 propagatorRBW( mass1sq, mass1[i], width1[i], skstr2, skaon, spion2, rRes2, 1, pro0_2 );
1824 }
1825 else if ( g0[i] == 0 )
1826 {
1827 pro0_1[0] = 1;
1828 pro0_1[1] = 0;
1829 pro0_2[0] = 1;
1830 pro0_2[1] = 0;
1831 pro2_1[0] = 1;
1832 pro2_1[1] = 0;
1833 pro2_2[0] = 1;
1834 pro2_2[1] = 0;
1835 }
1836 if ( g1[i] == 1 )
1837 {
1838 propagatorRBW( mass2sq, mass2[i], width2[i], sk11, skstr1, spi0, rRes2, 1, pro1_1 );
1839 propagatorRBW( mass2sq, mass2[i], width2[i], sk12, skstr2, spi0, rRes2, 1, pro1_2 );
1840 }
1841 else if ( g1[i] == 0 )
1842 {
1843 pro1_1[0] = 1;
1844 pro1_1[1] = 0;
1845 pro1_2[0] = 1;
1846 pro1_2[1] = 0;
1847 pro3_1[0] = 1;
1848 pro3_1[1] = 0;
1849 pro3_2[0] = 1;
1850 pro3_2[1] = 0;
1851 }
1852 Com_Multi( pro0_1, pro1_1, pro_1 );
1853 Com_Multi( pro2_1, pro3_1, pro4_1 );
1854 Com_Multi( pro0_2, pro1_2, pro_2 );
1855 Com_Multi( pro2_2, pro3_2, pro4_2 );
1856 amp_tmp1[0] = tmp1 * pro_1[0] - tmp3 * pro4_1[0];
1857 amp_tmp1[1] = tmp1 * pro_1[1] - tmp3 * pro4_1[1];
1858 amp_tmp2[0] = tmp2 * pro_2[0] - tmp4 * pro4_2[0];
1859 amp_tmp2[1] = tmp2 * pro_2[1] - tmp4 * pro4_2[1];
1860 }
1861
1862 else if ( modetype[i] == 12 )
1863 {
1864 B1_K1_Kstr3_Pi1 = Barrier( 1, sk11, skstr3, spion1, rRes2, mass2[i] );
1865 B1_K1_Kstr3_Pi2 = Barrier( 1, sk12, skstr3, spion2, rRes2, mass2[i] );
1866 B1_K1_Kstr1_Pi0 = Barrier( 1, sk11, skstr1, spi0, rRes2, mass2[i] );
1867 B1_K1_Kstr2_Pi0 = Barrier( 1, sk12, skstr2, spi0, rRes2, mass2[i] );
1868 B1_kstr3 = Barrier( 1, skstr3, skaon, spi0, rRes2, mass1[i] );
1869 B1_kstr1_Pi1 = Barrier( 1, skstr1, skaon, spion1, rRes2, mass1[i] );
1870 B1_kstr2_Pi2 = Barrier( 1, skstr2, skaon, spion2, rRes2, mass1[i] );
1871 for ( int w = 0; w < 4; w++ )
1872 {
1873 temp_PDF1 += G[w][w] * Pi0[w] * t1kstr1[w];
1874 temp_PDF2 += G[w][w] * Pi0[w] * t1kstr2[w];
1875 temp_PDF3 += G[w][w] * Pip1[w] * t1kstr3[w];
1876 temp_PDF4 += G[w][w] * Pip2[w] * t1kstr3[w];
1877 }
1878 tmp1 = B1_kstr1_Pi1 * B1_K1_Kstr1_Pi0 * temp_PDF1;
1879 tmp2 = B1_kstr2_Pi2 * B1_K1_Kstr2_Pi0 * temp_PDF2;
1880 tmp3 = B1_kstr3 * B1_K1_Kstr3_Pi1 * temp_PDF3;
1881 tmp4 = B1_kstr3 * B1_K1_Kstr3_Pi2 * temp_PDF4;
1882 propagatorRBW( mass1sq, mass1[i], width1[i], skstr3, skaon, spi0, rRes2, 1, pro2_1 );
1883 propagatorRBW( mass1sq, mass1[i], width1[i], skstr3, skaon, spi0, rRes2, 1, pro2_2 );
1884 propagatorRBW( mass2sq, mass2[i], width2[i], sk11, skstr3, spion1, rRes2, g2[i],
1885 pro3_1 );
1886 propagatorRBW( mass2sq, mass2[i], width2[i], sk12, skstr3, spion2, rRes2, g2[i],
1887 pro3_2 );
1888 if ( g0[i] == 1 )
1889 {
1890 propagatorRBW( mass1sq, mass1[i], width1[i], skstr1, skaon, spion1, rRes2, 1, pro0_1 );
1891 propagatorRBW( mass1sq, mass1[i], width1[i], skstr2, skaon, spion2, rRes2, 1, pro0_2 );
1892 }
1893 else if ( g0[i] == 0 )
1894 {
1895 pro0_1[0] = 1;
1896 pro0_1[1] = 0;
1897 pro0_2[0] = 1;
1898 pro0_2[1] = 0;
1899 pro2_1[0] = 1;
1900 pro2_1[1] = 0;
1901 pro2_2[0] = 1;
1902 pro2_2[1] = 0;
1903 }
1904 if ( g1[i] == 1 )
1905 {
1906 propagatorRBW( mass2sq, mass2[i], width2[i], sk11, skstr1, spi0, rRes2, g2[i],
1907 pro1_1 );
1908 propagatorRBW( mass2sq, mass2[i], width2[i], sk12, skstr2, spi0, rRes2, g2[i],
1909 pro1_2 );
1910 }
1911 else if ( g1[i] == 0 )
1912 {
1913 pro1_1[0] = 1;
1914 pro1_1[1] = 0;
1915 pro1_2[0] = 1;
1916 pro1_2[1] = 0;
1917 pro3_1[0] = 1;
1918 pro3_1[1] = 0;
1919 pro3_2[0] = 1;
1920 pro3_2[1] = 0;
1921 }
1922 Com_Multi( pro0_1, pro1_1, pro_1 );
1923 Com_Multi( pro2_1, pro3_1, pro4_1 );
1924 Com_Multi( pro0_2, pro1_2, pro_2 );
1925 Com_Multi( pro2_2, pro3_2, pro4_2 );
1926 amp_tmp1[0] = tmp1 * pro_1[0] - tmp3 * pro4_1[0];
1927 amp_tmp1[1] = tmp1 * pro_1[1] - tmp3 * pro4_1[1];
1928 amp_tmp2[0] = tmp2 * pro_2[0] - tmp4 * pro4_2[0];
1929 amp_tmp2[1] = tmp2 * pro_2[1] - tmp4 * pro4_2[1];
1930 }
1931 if ( modetype[i] == 13 )
1932 {
1933 B1_K1_rho1 = Barrier( 1, sk11, srho1, skaon, rRes2, mass2[i] );
1934 B1_K1_rho2 = Barrier( 1, sk12, srho2, skaon, rRes2, mass2[i] );
1935 B1_rho1 = Barrier( 1, srho1, spion1, spi0, rRes2, mass1[i] );
1936 B1_rho2 = Barrier( 1, srho2, spion2, spi0, rRes2, mass1[i] );
1937 for ( int w = 0; w < 4; w++ )
1938 {
1939 temp_PDF1 += G[w][w] * Km[w] * t1rho1[w];
1940 temp_PDF2 += G[w][w] * Km[w] * t1rho2[w];
1941 }
1942 tmp1 = B1_rho1 * B1_K1_rho1 * temp_PDF1;
1943 tmp2 = B1_rho2 * B1_K1_rho2 * temp_PDF2;
1944 if ( g1[i] == 1 )
1945 {
1946 propagatorRBW( mass2sq, mass2[i], width2[i], sk11, skaon, srho1, rRes2, 1, pro0_1 );
1947 propagatorRBW( mass2sq, mass2[i], width2[i], sk12, skaon, srho2, rRes2, 1, pro0_2 );
1948 }
1949 else if ( g1[i] == 0 )
1950 {
1951 pro0_1[0] = 1;
1952 pro0_1[1] = 0;
1953 pro0_2[0] = 1;
1954 pro0_2[1] = 0;
1955 }
1956 if ( g0[i] == 1 )
1957 {
1958 propagatorGS( mass1sq, mass1[i], width1[i], srho1, spion1, spi0, rRes2, pro1_1 );
1959 propagatorGS( mass1sq, mass1[i], width1[i], srho2, spion2, spi0, rRes2, pro1_2 );
1960 }
1961 else if ( g0[i] == 0 )
1962 {
1963 pro1_1[0] = 1;
1964 pro1_1[1] = 0;
1965 pro1_2[0] = 1;
1966 pro1_2[1] = 0;
1967 }
1968 Com_Multi( pro0_1, pro1_1, pro_1 );
1969 Com_Multi( pro0_2, pro1_2, pro_2 );
1970 amp_tmp1[0] = tmp1 * pro_1[0];
1971 amp_tmp1[1] = tmp1 * pro_1[1];
1972 amp_tmp2[0] = tmp2 * pro_2[0];
1973 amp_tmp2[1] = tmp2 * pro_2[1];
1974 }
1975 else if ( modetype[i] == 17 )
1976 {
1977 B1_rho1 = Barrier( 1, srho1, spion1, spi0, rRes2, mass2[i] );
1978 B1_rho2 = Barrier( 1, srho2, spion2, spi0, rRes2, mass2[i] );
1979 B1_D_Kstr2_rho1 = Barrier( 1, sD, skstr2, srho1, rD2, mD );
1980 B1_D_Kstr1_rho2 = Barrier( 1, sD, skstr1, srho2, rD2, mD );
1981 calt1( pKstr2, prho1, t1D_1 );
1982 calt1( pKstr1, prho2, t1D_2 );
1983 for ( int w = 0; w < 4; w++ )
1984 {
1985 temp_PDF1 += t1D_1[w] * t1rho1[w] * G[w][w];
1986 temp_PDF2 += t1D_2[w] * t1rho2[w] * G[w][w];
1987 }
1988 tmp1 = temp_PDF1 * B1_rho1 * B1_D_Kstr2_rho1;
1989 tmp2 = temp_PDF2 * B1_rho2 * B1_D_Kstr1_rho2;
1990 if ( g1[i] == 1 )
1991 {
1992 propagatorGS( mass2sq, mass2[i], width2[i], srho1, spion1, spi0, rRes2, pro_1 );
1993 propagatorGS( mass2sq, mass2[i], width2[i], srho2, spion2, spi0, rRes2, pro_2 );
1994 }
1995 else if ( g1[i] == 0 )
1996 {
1997 pro_1[0] = 1;
1998 pro_1[1] = 0;
1999 pro_2[0] = 1;
2000 pro_2[1] = 0;
2001 }
2002 temp_1[0] = tmp1 * pro_1[0];
2003 temp_1[1] = tmp1 * pro_1[1];
2004 temp_2[0] = tmp2 * pro_2[0];
2005 temp_2[1] = tmp2 * pro_2[1];
2006 // printf("xx: %.8f\n", temp_2[0]);
2007 // printf("xxxx: %.8f\n", temp_2[1]);
2008 KPiSLASS( skstr1, skaon, spion1, Amp_KPiS1 );
2009 KPiSLASS( skstr2, skaon, spion2, Amp_KPiS2 );
2010 Com_Multi( temp_1, Amp_KPiS2, amp_tmp1 );
2011 Com_Multi( temp_2, Amp_KPiS1, amp_tmp2 );
2012 }
2013 else if ( modetype[i] == 9 )
2014 {
2015 B1_kstr1_Pi1 = Barrier( 1, skstr1, skaon, spion1, rRes2, mass1[i] );
2016 B1_kstr2_Pi2 = Barrier( 1, skstr2, skaon, spion2, rRes2, mass1[i] );
2017 B1_D_K1p_Pi0 = Barrier( 1, sD, sk13, spi0, rD2, mD );
2018 calt1( pK13, Pi0, t1D_1 );
2019 if ( g2[i] == 0 )
2020 {
2021 for ( int w = 0; w < 4; w++ )
2022 {
2023 for ( int j = 0; j < 4; j++ )
2024 {
2025 temp_PDF1 += t1D_1[w] * ( G[w][j] - pK13[w] * pK13[j] / sk13 ) * t1kstr1[j] *
2026 G[w][w] * G[j][j];
2027 temp_PDF2 += t1D_1[w] * ( G[w][j] - pK13[w] * pK13[j] / sk13 ) * t1kstr2[j] *
2028 G[w][w] * G[j][j];
2029 }
2030 }
2031 tmp1 = B1_kstr1_Pi1 * B1_D_K1p_Pi0 * temp_PDF1;
2032 tmp2 = B1_kstr2_Pi2 * B1_D_K1p_Pi0 * temp_PDF2;
2033 }
2034 else if ( g2[i] == 2 )
2035 {
2036 for ( int w = 0; w < 4; w++ )
2037 {
2038 for ( int j = 0; j < 4; j++ )
2039 {
2040 temp_PDF1 += t1D_1[w] * t2kk1[w][j] * t1kstr1[j] * G[w][w] * G[j][j];
2041 temp_PDF2 += t1D_1[w] * t2kk2[w][j] * t1kstr2[j] * G[w][w] * G[j][j];
2042 }
2043 }
2044 if ( B2_K1p_kstr1 < 0.0 )
2045 B2_K1p_kstr1 = Barrier( 2, sk13, skstr1, spion2, rRes2, mass2[i] );
2046 if ( B2_K1p_kstr2 < 0.0 )
2047 B2_K1p_kstr2 = Barrier( 2, sk13, skstr2, spion1, rRes2, mass2[i] );
2048 tmp1 = B1_kstr1_Pi1 * B2_K1p_kstr1 * B1_D_K1p_Pi0 * temp_PDF1;
2049 tmp2 = B1_kstr2_Pi2 * B2_K1p_kstr2 * B1_D_K1p_Pi0 * temp_PDF2;
2050 }
2051 if ( g0[i] == 1 )
2052 {
2053 propagatorRBW( mass1sq, mass1[i], width1[i], skstr1, skaon, spion1, rRes2, 1, pro0_1 );
2054 propagatorRBW( mass1sq, mass1[i], width1[i], skstr2, skaon, spion2, rRes2, 1, pro0_2 );
2055 }
2056 else if ( g0[i] == 0 )
2057 {
2058 pro0_1[0] = 1;
2059 pro0_1[1] = 0;
2060 pro0_2[0] = 1;
2061 pro0_2[1] = 0;
2062 }
2063 if ( g1[i] == 1 )
2064 {
2065 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk13, skstr1, spion2, rRes2, g2[i],
2066 pro1_1 );
2067 propagatorRBW_K1( mass2sq, mass2[i], width2[i], sk13, skstr2, spion1, rRes2, g2[i],
2068 pro1_2 );
2069 }
2070 else if ( g1[i] == 0 )
2071 {
2072 pro1_1[0] = 1;
2073 pro1_1[1] = 0;
2074 pro1_2[0] = 1;
2075 pro1_2[1] = 0;
2076 }
2077 Com_Multi( pro0_1, pro1_1, pro_1 );
2078 Com_Multi( pro0_2, pro1_2, pro_2 );
2079 amp_tmp1[0] = tmp1 * pro_1[0];
2080 amp_tmp1[1] = tmp1 * pro_1[1];
2081 amp_tmp2[0] = tmp2 * pro_2[0];
2082 amp_tmp2[1] = tmp2 * pro_2[1];
2083 }
2084
2085 amp_tmp[0] = amp_tmp1[0] + amp_tmp2[0];
2086 amp_tmp[1] = amp_tmp1[1] + amp_tmp2[1];
2087
2088 Com_Multi( amp_tmp, cof, amp_PDF );
2089 PDF[0] += amp_PDF[0];
2090 PDF[1] += amp_PDF[1];
2091 double value = PDF[0] * PDF[0] + PDF[1] * PDF[1];
2092 if ( value <= 0 ) { value = 1e-20; }
2093 Result = value;
2094 }
2095}
mass
double mass
Definition CosmicGenerator.cxx:128
f
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
g1
TF1 * g1
Definition DataBase/tau_mode.c:58
K1270_width
double K1270_width[]
Definition EvtDToKmPipPipPi0.cc:364
K1400_width
double K1400_width[]
Definition EvtDToKmPipPipPi0.cc:650
a1_width
double a1_width[]
Definition EvtDToKmPipPipPi0.cc:47
w
double w
Definition EvtPhokharaDef.hh:28
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
CS
TCrossPart * CS
Definition Mcgpj.cxx:51
EvtDToKmPipPipPi0::decay
void decay(EvtParticle *p)
Definition EvtDToKmPipPipPi0.cc:1046
EvtDToKmPipPipPi0::getName
void getName(std::string &name)
Definition EvtDToKmPipPipPi0.cc:43
EvtDToKmPipPipPi0::clone
EvtDecayBase * clone()
Definition EvtDToKmPipPipPi0.cc:45
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:492
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:261
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:40
EvtPDL::getStdHep
static int getStdHep(EvtId id)
Definition EvtPDL.hh:61
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:100
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:171
t
int t()
Definition t.c:1