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