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