BOSS 7.1.1
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDToKSKmpippip.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: EvtDToKSKmpippip.cc
11// the necessary file: EvtDToKSKmpippip.hh
12//
13// Description: D+ -> KS K- pi+ pi+
14//
15// Modification history:
16//
17// Zhenxuan Li Jun. 15, 2023 Module created
18// Jiahao Wang Nov. 03, 2023 modif
19//------------------------------------------------------------------------
20#include "EvtGenBase/EvtVector4R.hh"
21#include "EvtGenBase/EvtPatches.hh"
22#include "EvtGenBase/EvtParticle.hh"
23#include "EvtGenBase/EvtGenKine.hh"
24#include "EvtGenBase/EvtPDL.hh"
25#include "EvtGenBase/EvtReport.hh"
26#include "EvtGenModels/EvtDToKSKmpippip.hh"
27#include "EvtGenBase/EvtComplex.hh"
28#include "EvtGenBase/EvtDecayTable.hh"
29#include <stdlib.h>
30#include <iostream>
31#include <cmath>
32using namespace std;
33
34EvtDToKSKmpippip::~EvtDToKSKmpippip() {}
35
36void EvtDToKSKmpippip::getName(std::string& model_name){
37 model_name="DToKSKmpippip";
38}
39
43
44void EvtDToKSKmpippip::init(){
45 checkNArg(0);
46 checkNDaug(4);
47 checkSpinParent(EvtSpinType::SCALAR);
48 //cout << "Initializing EvtD0ToKSKppimpi0" << endl;
49mass_Pion = 0.13957;
50mass_Pion_N = 0.134977;
51mass_Eta = 0.547862;
52mD = 1.86966; mPi = 0.13957; mKa = 0.493677;
53math_pi = 3.1415926;
54rD2 = 25.0;
55rRes2 = 9.0;
56GS1 = 0.636619783;
57GS2 = 0.01860182466;
58GS3 = 0.1591549458;
59GS4 = 0.00620060822;
60
61 mK1400 = 1.403; GK1400 = 0.174;
62 m1510 = 1.518; G1510 = 0.073;
63 meta1405 =1.4139; Geta1405 = 0.048;
64 meta1475 = 1.475; Geta1475 = 0.09;
65 mf1285 = 1.2819; Gf1285 = 0.0227;
66 mKstp =0.89176; GKstp = 0.0514;
67 mKst0 =0.89555; GKst0 = 0.0473;
68 ma0_980 = 0.990; Ga0_980 = 0.05;
69
70 rho[0] = 1.0; //2S
71 phi[0] = 0.0;
72
73 phi[1] = 4.1350e-02;
74 phi[2] =-8.3438e-01;
75 phi[3] =-7.9307e-01;
76 phi[4] =-2.2340e+00;
77 phi[5] =-3.5357e-01;
78 phi[6] =-3.5357e-01;
79 phi[7] = 9.1698e+00;
80 phi[8] = 2.1065e+00;
81 rho[1] = 5.5001e-01; //1S
82 rho[2] = 8.9374e-01; //2P
83 rho[3] = 1.8706e+00; //16S
84 rho[4] = 1.2629e+01; //201
85 rho[5] = 3.7866e+00; //201
86 rho[6] = 3.7866e+00; //101
87 rho[7] = 8.6436e-01; //101
88 rho[8] = 2.1213e+00; //101
89
90 modetype[0]=2;
91 modetype[1]=2;
92 modetype[2]=2;
93 modetype[3]=15;
94 modetype[4]=1500;
95 modetype[5]=1200;
96 modetype[6]=1100;
97 modetype[7]=9;
98 modetype[8]=20;
99
100
101 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
102 int EE[4][4][4][4] =
103 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
104 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
105 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
106 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
107 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
108 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
109 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
110 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
111 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
112 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
113 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
114 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
115 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
116 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
117 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
118 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
119 for (int i=0; i<4; i++) {
120 for (int j=0; j<4; j++) {
121 G[i][j] = GG[i][j];
122 for (int k=0; k<4; k++) {
123 for (int l=0; l<4; l++) {
124 E[i][j][k][l] = EE[i][j][k][l];
125 }
126 }
127 }
128 }
129}
130
134
135void EvtDToKSKmpippip::decay( EvtParticle *p){
136 //-----------for max value------------------
137///
138////
139/*
140 double maxprob = 0.0;
141 for(int ir=0;ir<=60000000;ir++){
142 p->initializePhaseSpace(getNDaug(),getDaugs());
143 EvtVector4R ks = p->getDaug(0)->getP4();
144 EvtVector4R km = p->getDaug(1)->getP4();
145 EvtVector4R pip1 = p->getDaug(2)->getP4();
146 EvtVector4R pip2 = p->getDaug(3)->getP4();
147 mother_c=EvtPDL::getStdHep(p->getId());
148 int cc;
149 if(mother_c==411){
150 cc=1;
151 } else if(mother_c==-411){
152 cc=-1;
153 }
154 double value;
155 double KS[4],Km[4],Pip1[4],Pip2[4];
156
157 KS[0] = ks.get(0); Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);
158 KS[1] = cc*ks.get(1); Km[1] = cc*km.get(1); Pip1[1] = cc*pip1.get(1); Pip2[1] = cc*pip2.get(1);
159 KS[2] = cc*ks.get(2); Km[2] = cc*km.get(2); Pip1[2] = cc*pip1.get(2); Pip2[2] = cc*pip2.get(2);
160 KS[3] = cc*ks.get(3); Km[3] = cc*km.get(3); Pip1[3] = cc*pip1.get(3); Pip2[3] = cc*pip2.get(3);
161 double mass1[9] ={ mKstp, mKstp, mKstp, mf1285, m1510, m1510, m1510, meta1405, meta1475};
162 double mass2[9] ={ mKst0, mKst0, mKst0, ma0_980, ma0_980, mKst0, mKstp, ma0_980, mKstp };
163 double width1[9]={ GKstp, GKstp, GKstp, Gf1285, G1510, G1510, G1510, Geta1405, Geta1475};
164 double width2[9]={ GKst0, GKst0, GKst0, Ga0_980, Ga0_980, GKst0, GKstp, Ga0_980, GKstp };
165 int g0[9] ={ 1, 1, 1, 1, 1, 1, 1, 1, 1 };
166 int g1[9] ={ 1, 1, 1, 1, 1, 1, 1, 1, 1 };
167 int g2[9] ={ 0, 1, 2, 0, 1, 0, 0, 0, 0 };
168
169 value=calEva(KS,Km,Pip1,Pip2,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,value);
170 if(value>maxprob) {
171 maxprob=value;
172 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
173 }
174 }
175 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
176 return;
177 //-----------------------------------------------
178*/
179 p->initializePhaseSpace(getNDaug(),getDaugs());
180 EvtVector4R ks = p->getDaug(0)->getP4();
181 EvtVector4R km = p->getDaug(1)->getP4();
182 EvtVector4R pip1= p->getDaug(2)->getP4();
183 EvtVector4R pip2= p->getDaug(3)->getP4();
184
185 mother_c=EvtPDL::getStdHep(p->getId());
186 int cc;
187 if(mother_c==411){
188 cc=1;
189 } else if(mother_c==-411){
190 cc=-1;
191 }
192
193 double KS[4],Km[4],Pip1[4],Pip2[4];
194
195 KS[0] = ks.get(0); Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);
196 KS[1] = cc*ks.get(1); Km[1] = cc*km.get(1); Pip1[1] = cc*pip1.get(1); Pip2[1] = cc*pip2.get(1);
197 KS[2] = cc*ks.get(2); Km[2] = cc*km.get(2); Pip1[2] = cc*pip1.get(2); Pip2[2] = cc*pip2.get(2);
198 KS[3] = cc*ks.get(3); Km[3] = cc*km.get(3); Pip1[3] = cc*pip1.get(3); Pip2[3] = cc*pip2.get(3);
199
200// KS[0] = 0.5844057460 ; Km[0] = 0.6644209662 ; Pip1[0] = 0.3998331690 ; Pip2[0] = 0.2369772689 ;
201// KS[1] =-0.0451579928 ; Km[1] =-0.1621486835 ; Pip1[1] = 0.2182356619 ; Pip2[1] = 0.0516153672 ;
202// KS[2] =-0.2455917598 ; Km[2] = 0.3728754865 ; Pip1[2] =-0.1884855750 ; Pip2[2] =-0.1706218636 ;
203// KS[3] =-0.1776479997 ; Km[3] =-0.1800275423 ; Pip1[3] = 0.2392343060 ; Pip2[3] = 0.0700168961 ;
204
205 double value;
206 double mass1[9] ={ mKstp, mKstp, mKstp, mf1285, m1510, m1510, m1510, meta1405, meta1475};
207 double mass2[9] ={ mKst0, mKst0, mKst0, ma0_980, ma0_980, mKst0, mKstp, ma0_980, mKstp };
208 double width1[9]={ GKstp, GKstp, GKstp, Gf1285, G1510, G1510, G1510, Geta1405, Geta1475};
209 double width2[9]={ GKst0, GKst0, GKst0, Ga0_980, Ga0_980, GKst0, GKstp, Ga0_980, GKstp };
210 int g0[9] ={ 1, 1, 1, 1, 1, 1, 1, 1, 1 };
211 int g1[9] ={ 1, 1, 1, 1, 1, 1, 1, 1, 1 };
212 int g2[9] ={ 0, 1, 2, 0, 1, 0, 0, 0, 0 };
213
214 value=calEva(KS,Km,Pip1,Pip2,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,value);
215 setProb(value);
216
217return;
218
219}
220
221void EvtDToKSKmpippip::Com_Multi(double a1[2], double a2[2], double res[2])
222{
223 res[0] = a1[0]*a2[0]-a1[1]*a2[1];
224 res[1] = a1[1]*a2[0]+a1[0]*a2[1];
225}
226void EvtDToKSKmpippip::Com_Divide(double a1[2], double a2[2], double res[2])
227{
228 double tmp = a2[0]*a2[0]+a2[1]*a2[1];
229 res[0] = (a1[0]*a2[0]+a1[1]*a2[1])/tmp;
230 res[1] = (a1[1]*a2[0]-a1[0]*a2[1])/tmp;
231}
232//------------base---------------------------------
233double EvtDToKSKmpippip::SCADot(double a1[4], double a2[4])
234{
235 double _cal = a1[0]*a2[0]-a1[1]*a2[1]-a1[2]*a2[2]-a1[3]*a2[3];
236 return _cal;
237}
238
239double EvtDToKSKmpippip::Barrier(double mass2, int l, double sa, double sb, double sc, double r2)
240{
241 double F;
242 double tmp = sa+sb-sc;
243 double q = fabs(0.25*tmp*tmp/sa-sb);
244 //if (q < 0) q = 1e-16;
245 double tmp2 = mass2+sb-sc;
246 double q0 = fabs(0.25*tmp2*tmp2/mass2-sb);
247 //if (q0 < 0) q0 = 1e-16;
248 double z = q*r2;
249 double z0 = q0*r2;
250 if (l==1) {
251 F = sqrt((1.0+z0)/(1.0+z));
252 }
253 else if (l==2) {
254 double z2 = z*z; double z02 = z0*z0;
255 F = sqrt((9.0+3.0*z0+z02)/(9.0+3.0*z+z2));
256 } else {
257 F = 1.0;
258 }
259 return F;
260}
261
262double EvtDToKSKmpippip::barrier(int l, double sa, double sb, double sc, double r2)
263{
264 double F;
265 double tmp = sa+sb-sc;
266 double q = 0.25*tmp*tmp/sa-sb;
267 if (q < 0) q = 1e-16;
268 double z = q*r2;
269 if (l == 1) {
270 F = sqrt(2.0*z/(1.0+z));
271 }
272 else if (l == 2) {
273 double z2 = z*z;
274 F = sqrt(13.0*z2/(9.0+3.0*z+z2));
275 } else {
276 F = 1.0;
277 }
278 return F;
279}
280//------------------spin-------------------------------------------
281void EvtDToKSKmpippip::calt1(double daug1[4], double daug2[4], double t1[4])
282{
283 double p, pq, tmp;
284 double pa[4], qa[4];
285 for(int i=0; i<4; i++) {
286 pa[i] = daug1[i] + daug2[i];
287 qa[i] = daug1[i] - daug2[i];
288 }
289 p = SCADot(pa,pa);
290 pq = SCADot(pa,qa);
291 tmp = pq/p;
292 for(int i=0; i<4; i++) {
293 t1[i] = qa[i] - tmp*pa[i];
294 }
295}
296void EvtDToKSKmpippip::calt2(double daug1[4], double daug2[4], double t2[4][4])
297{
298 double p, r;
299 double pa[4], t1[4];
300 calt1(daug1,daug2,t1);
301 r = SCADot(t1,t1)/3.0;
302 for(int i=0; i<4; i++) {
303 pa[i] = daug1[i] + daug2[i];
304 }
305 p = SCADot(pa,pa);
306 for(int i=0; i<4; i++) {
307 for(int j=0; j<4; j++) {
308 t2[i][j] = t1[i]*t1[j] - r*(G[i][j]-pa[i]*pa[j]/p);
309 }
310 }
311}
312//-------------------prop--------------------------------------------
313void EvtDToKSKmpippip::Flatte_rhoab(double sa, double sb, double sc, double rho[2])
314{
315 double q = (sa+sb-sc)*(sa+sb-sc)/(4*sa)-sb;
316 if(q>0) {
317 rho[0]=2* sqrt(q/sa);
318 rho[1]=0;
319 }
320 else if(q<0){
321 rho[0]=0;
322 rho[1]=2*sqrt(-q/sa);
323 }
324}
325void EvtDToKSKmpippip::propagator980(double mass2, double mass, double sx, double *sb, double *sc, double prop[2])
326{
327 double unit[2]={1.0};
328 double ci[2]={0,1};
329 double rho1[2];
330 Flatte_rhoab(sx,sb[0],sc[0],rho1);
331 double rho2[2];
332 Flatte_rhoab(sx,sb[1],sc[1],rho2);
333
334 double gK_f980=0.69465, gPi_f980=0.165;
335 double tmp1[2]={gK_f980,0};
336 double tmp11[2];
337 double tmp2[2]={gPi_f980,0};
338 double tmp22[2];
339 Com_Multi(tmp1,rho1,tmp11);
340 Com_Multi(tmp2,rho2,tmp22);
341 double tmp3[2]={tmp11[0]+tmp22[0],tmp11[1]+tmp22[1]};
342 double tmp31[2];
343 Com_Multi(tmp3, ci,tmp31);
344 double tmp4[2]={mass2-sx-tmp31[0], -1.0*tmp31[1]};
345 Com_Divide( unit,tmp4, prop);
346}
347void EvtDToKSKmpippip::propagatora0980(double mass2, double mass, double sx, double *sb, double *sc, double prop[2])
348{
349 double unit[2]={1.0};
350 double ci[2]={0,1};
351 double rho1[2];
352 Flatte_rhoab(sx,sb[0],sc[0],rho1);
353 double rho2[2];
354 Flatte_rhoab(sx,sb[1],sc[1],rho2);
355
356 double gKK_a980=0.892*0.341, gPiEta_a980=0.341;
357 double tmp1[2]={gKK_a980,0};
358 double tmp11[2];
359 double tmp2[2]={gPiEta_a980,0};
360 double tmp22[2];
361 Com_Multi(tmp1,rho1,tmp11);
362 Com_Multi(tmp2,rho2,tmp22);
363 double tmp3[2]={tmp11[0]+tmp22[0],tmp11[1]+tmp22[1]};
364 double tmp31[2];
365 Com_Multi(tmp3, ci,tmp31);
366 double tmp4[2]={mass2-sx-tmp31[0], -1.0*tmp31[1]};
367 Com_Divide( unit,tmp4, prop);
368}
369//--------------------------------------------------------------------------------------------------------
370void EvtDToKSKmpippip::propagator(double mass2, double mass, double width, double sx, double prop[2])
371{
372 double a[2], b[2];
373 a[0] = 1;
374 a[1] = 0;
375 b[0] = mass2-sx;
376 b[1] = -mass*width;
377 Com_Divide(a,b,prop);
378}
379double EvtDToKSKmpippip::wid(double mass2, double mass, double sa, double sb, double sc, double r2, int l)
380{
381 double widm = 0.;
382 double m = sqrt(sa);
383 double tmp = sb-sc;
384 double tmp1 = sa+tmp;
385 double q = 0.25*tmp1*tmp1/sa-sb;
386 if(q<0) q = 1e-16;
387 double tmp2 = mass2+tmp;
388 double q0 = 0.25*tmp2*tmp2/mass2-sb;
389 if(q0<0) q0 = 1e-16;
390 double z = q*r2;
391 double z0 = q0*r2;
392 double t = q/q0;
393 if(l == 0) {widm = sqrt(t)*mass/m;}
394 else if(l == 1) {widm = t*sqrt(t)*mass/m*(1+z0)/(1+z);}
395 else if(l == 2) {widm = t*t*sqrt(t)*mass/m*(9+3*z0+z0*z0)/(9+3*z+z*z);}
396 return widm;
397}
398double EvtDToKSKmpippip::widl1(double mass2, double mass, double sa, double sb, double sc, double r2)
399{
400 double widm = 0.;
401 double m = sqrt(sa);
402 double tmp = sb-sc;
403 double tmp1 = sa+tmp;
404 double q = 0.25*tmp1*tmp1/sa-sb;
405 if(q<0) q = 1e-16;
406 double tmp2 = mass2+tmp;
407 double q0 = 0.25*tmp2*tmp2/mass2-sb;
408 if(q0<0) q0 = 1e-16;
409 double z = q*r2;
410 double z0 = q0*r2;
411 double F = (1+z0)/(1+z);
412 double t = q/q0;
413 widm = t*sqrt(t)*mass/m*F;
414 return widm;
415}
416void EvtDToKSKmpippip::propagatorRBW(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2])
417{
418 double a[2], b[2];
419 a[0] = 1;
420 a[1] = 0;
421 b[0] = mass2-sa;
422 b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l);
423 Com_Divide(a,b,prop);
424}
425void EvtDToKSKmpippip::propagatorRBWl1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2])
426{
427 double a[2], b[2];
428 a[0] = 1;
429 a[1] = 0;
430 b[0] = mass2-sa;
431 b[1] = -mass*width*widl1(mass2,mass,sa,sb,sc,r2);
432 Com_Divide(a,b,prop);
433}
434//------------GS---used by rho----------------------------
435void EvtDToKSKmpippip::propagatorGS(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2])
436{
437 double a[2], b[2];
438 double tmp = sb-sc;
439 double tmp1 = sa+tmp;
440 double q2 = 0.25*tmp1*tmp1/sa-sb;
441 if(q2<0) q2 = 1e-16;
442
443 double tmp2 = mass2+tmp;
444 double q02 = 0.25*tmp2*tmp2/mass2-sb;
445 if(q02<0) q02 = 1e-16;
446
447 double q = sqrt(q2);
448 double q0 = sqrt(q02);
449 double m = sqrt(sa);
450 double q03 = q0*q02;
451 double tmp3 = log(mass+2*q0)+1.2760418309; // log(mass_2Pion) = 1.2760418309;
452
453 double h = GS1*q/m*(log(m+2*q)+1.2760418309);
454 double h0 = GS1*q0/mass*tmp3;
455 double dh = h0*(0.125/q02-0.5/mass2)+GS3/mass2;
456 double d = GS2/q02*tmp3+GS3*mass/q0-GS4*mass/q03;
457 double f = mass2/q03*(q2*(h-h0)+(mass2-sa)*q02*dh);
458
459 a[0] = 1.0+d*width/mass;
460 a[1] = 0.0;
461 b[0] = mass2-sa+width*f;
462 b[1] = -mass*width*widl1(mass2,mass,sa,sb,sc,r2);
463 Com_Divide(a,b,prop);
464}
465void EvtDToKSKmpippip::rhoab(double sa, double sb, double sc, double res[2]) {
466 double tmp = sa+sb-sc;
467 double q = 0.25*tmp*tmp/sa-sb;
468 if(q>=0) {
469 res[0]=2.0*sqrt(q/sa);
470 res[1]=0.0;
471 } else {
472 res[0]=0.0;
473 res[1]=2.0*sqrt(-q/sa);
474 }
475}
476void EvtDToKSKmpippip::rho4Pi(double sa, double res[2]) {
477 double temp = 1.0-0.3116765584/sa; // 0.3116765584=0.13957*0.13957*16
478 if(temp>=0) {
479 res[0]=sqrt(temp)/(1.0+exp(9.8-3.5*sa));
480 res[1]=0.0;
481 } else {
482 res[0]=0.0;
483 res[1]=sqrt(-temp)/(1.0+exp(9.8-3.5*sa));
484 }
485}
486void EvtDToKSKmpippip::propagatorsigma500(double sa, double sb, double sc, double prop[2]) {//for pipi_s wave
487 double f = 0.5843+1.6663*sa;
488 const double M = 0.9264;
489 const double mass2 = 0.85821696; // M*M
490 const double mpi2d2 = 0.00973989245;
491 double g1 = f*(sa-mpi2d2)/(mass2-mpi2d2)*exp((mass2-sa)/1.082);
492 double rho1s[2], rho1M[2], rho2s[2], rho2M[2], rho1[2], rho2[2];
493 rhoab(sa,sb,sc,rho1s);
494 rhoab(mass2,sb,sc,rho1M);
495 rho4Pi(sa,rho2s);
496 rho4Pi(mass2,rho2M);
497 Com_Divide(rho1s,rho1M,rho1);
498 Com_Divide(rho2s,rho2M,rho2);
499 double a[2], b[2];
500 a[0] = 1.0;
501 a[1] = 0.0;
502 b[0] = mass2-sa+M*(g1*rho1[1]+0.0024*rho2[1]);
503 b[1] = -M*(g1*rho1[0]+0.0024*rho2[0]);
504 Com_Divide(a,b,prop);
505}
506void EvtDToKSKmpippip::KPiSLASS(double sa, double sb, double sc, double prop[2]) {
507 const double m1430 = 1.441;
508 const double sa0 = 2.076481; // m1430*m1430;
509 const double w1430 = 0.193;
510 const double Lass1 = 0.25/sa0;
511 double tmp = sb-sc;
512 double tmp1 = sa0+tmp;
513 double q0 = fabs(Lass1*tmp1*tmp1-sb);
514 //if(q0<0) q0 = 1e-16;
515 double tmp2 = sa+tmp;
516 double qs = fabs(0.25*tmp2*tmp2/sa-sb);
517 double q = sqrt(qs);
518 double width = w1430*q*m1430/sqrt(sa*q0);
519 double temp_R = atan(m1430*width/(sa0-sa));
520 if(temp_R<0) temp_R += math_pi;
521 double deltaR = -109.7*math_pi/180.0 + temp_R;
522 double temp_F = atan(0.226*q/(2.0-3.8194*qs)); // 2.0*0.113 = 0.226; -33.8*0.113 = -3.8194
523 if(temp_F<0) temp_F += math_pi;
524 double deltaF = 0.1*math_pi/180.0 + temp_F;
525 double deltaS = deltaR + 2.0*deltaF;
526 double t1 = 0.96*sin(deltaF);
527 double t2 = sin(deltaR);
528 double CF[2], CS[2];
529 CF[0] = cos(deltaF);
530 CF[1] = sin(deltaF);
531 CS[0] = cos(deltaS);
532 CS[1] = sin(deltaS);
533 prop[0] = t1*CF[0] + t2*CS[0];
534 prop[1] = t1*CF[1] + t2*CS[1];
535
536}
537
538//---------------------------------------------------------------------------------------------------------------------
539double EvtDToKSKmpippip::calEva(
540 double* Km,
541 double* Kp,
542 double* Pip,
543 double* Pi0,
544 double* mass1, double* mass2,
545 double* width1, double* width2,
546 double* amp, double* phase,
547 int* g0, int* g1, int* g2,
548 int* modetype,
549 double Result)
550//r0,r1 just send propagatorGS_Omg(mass1sq,mass1[i],width1[i],srho1,spion1,spi0,rRes2,amp_Omg,phs_Omg,proRho1);
551
552{
553 double cof[2], amp_tmp1[2], amp_tmp2[2], amp_tmp[2], amp_PDF[2], PDF[2];
554 double pKstr0[4],pKstrp[4],pKstrm[4],prhop[4],pphi[4],pK10[4],pK11[4],pK12[4],pK13[4],pK14[4];//pD[4];
555 double pD[4],pKsPip1[4],pKmPip2[4],pKsPip2[4],pKmPip1[4],pKsKmPip1[4],pKsKmPip2[4],pKsKm[4],pKmPipPip[4],pKsPipPip[4];
556
557 for(int i=0;i!=4;i++)
558 {
559 pD[i] = Km[i]+Pip[i]+Kp[i]+Pi0[i];
560 pKsPip1[i] =Km[i]+Pip[i];
561 pKsPip2[i] =Km[i]+Pi0[i];
562 pKmPip1[i] =Kp[i]+Pip[i];
563 pKmPip2[i] =Kp[i]+Pi0[i];
564 pKsKm[i] = Km[i]+Kp[i];
565 pKsKmPip1[i] = pKsPip1[i] + Kp[i];
566 pKsKmPip2[i] = pKsPip2[i] + Kp[i];
567 pKmPipPip[i] = Kp[i]+Pip[i]+Pi0[i];
568 pKsPipPip[i] = Km[i]+Pip[i]+Pi0[i];
569
570 pphi[i] =Km[i]+Kp[i];
571 prhop[i] =Pip[i]+Pi0[i];
572 pKstr0[i] =Km[i]+Pip[i];
573 pKstrp[i] =Kp[i]+Pi0[i];
574 pKstrm[i] =Km[i]+Pi0[i];
575 pK10[i] =pKstrm[i]+Kp[i];
576 pK11[i] =pKstrm[i]+Pip[i];
577 pK12[i] =pKstrp[i]+Km[i];
578 pK13[i] =pKstr0[i]+Pi0[i]; //can delete, no use, equal to pK11
579 pK14[i] =pphi[i]+Pi0[i];
580 }
581 double spi0,spionp,skaon1,skaon2,sphi,srhop,skstr0,skstrp,skstrm,sk10,sk11,sk12,sk13,sk14;//,sD;
582 double sPip2,sPip1,sKs,sKm,sKsPip1,sKmPip2,sD,sKsPip2,sKmPip1,sKsKmPip1,sKsKmPip2,sKsKm,sKmPipPip,sKsPipPip;
583
584 sKs = SCADot(Km,Km);
585 sKm = SCADot(Kp,Kp);
586 sPip1= SCADot(Pip,Pip);
587 sPip2= SCADot(Pi0,Pi0);
588 sKsPip1 = SCADot(pKsPip1,pKsPip1);
589 sKsPip2 = SCADot(pKsPip2,pKsPip2);
590 sKmPip1 = SCADot(pKmPip1,pKmPip1);
591 sKmPip2 = SCADot(pKmPip2,pKmPip2);
592 sKsKmPip1 = SCADot(pKsKmPip1,pKsKmPip1);
593 sKsKmPip2 = SCADot(pKsKmPip2,pKsKmPip2);
594 sKmPipPip = SCADot(pKmPipPip,pKmPipPip);
595 sKsPipPip = SCADot(pKsPipPip,pKsPipPip);
596 sKsKm =SCADot(pKsKm,pKsKm);
597
598
599 skaon1 = SCADot(Km,Km);
600 skaon2 = SCADot(Kp,Kp);
601 spionp = SCADot(Pip,Pip);
602 spi0 = SCADot(Pi0,Pi0);
603
604 sphi = SCADot(pphi,pphi);
605 srhop = SCADot(prhop,prhop);
606 skstr0 = SCADot(pKstr0,pKstr0);
607 skstrp = SCADot(pKstrp,pKstrp);
608 skstrm = SCADot(pKstrm,pKstrm);
609 sk10 = SCADot(pK10,pK10);
610 sk11 = SCADot(pK11,pK11);
611 sk12 = SCADot(pK12,pK12);
612// sk13 = SCADot(pK13,pK13);
613 sk14 = SCADot(pK14,pK14);
614 sD = SCADot(pD,pD);
615//-------------------------------------------------------------------------------------------------
616 double t1phi[4],t1A1[4],t1A2[4],t1rhop[4],t1kstr1[4],t1kstr2[4],t1kstr3[4],t1kstr4[4],t2k11[4][4],t2k12[4][4],t2k13[4][4],t2k14[4][4],t2k21[4][4];
617 double t1KsPip1[4],t1KmPip2[4],t1KsPip2[4],t1KmPip1[4],t2KsKmPip1_kspip[4][4],t2KsKmPip2_kspip[4][4],t2KsKmPip1_kmpip[4][4],t2KsKmPip2_kmpip[4][4],t1KsKmPip1_kspip[4],t1KsKmPip2_kspip[4],t1KsKmPip1_kmpip[4],t1KsKmPip2_kmpip[4],t1KsKmPip1_ksk[4],t1KsKmPip2_ksk[4],t2KmPipPip_kmpip[4][4],t2KmPipPip_kmpip2[4][4],t2KsPipPip_kspip[4][4],t2KsPipPip_kspip2[4][4],t1KmPipPip_kmpip[4],t1KmPipPip_kmpip2[4],t1KsPipPip_kspip[4],t1KsPipPip_kspip2[4];
618 calt1(Km,Pip, t1KsPip1);
619 calt1(Km,Pi0, t1KsPip2);
620 calt1(Kp,Pip, t1KmPip1);
621 calt1(Kp,Pi0, t1KmPip2);
622 calt1(pKsPip1,Kp,t1KsKmPip1_kspip);
623 calt1(pKsPip2,Kp,t1KsKmPip2_kspip);
624 calt1(pKmPip1,Km,t1KsKmPip1_kmpip);
625 calt1(pKmPip2,Km,t1KsKmPip2_kmpip);
626
627 calt1(pKmPip1,Pi0,t1KmPipPip_kmpip);
628 calt1(pKmPip2,Pip,t1KmPipPip_kmpip2);
629 calt1(pKsPip1,Pi0,t1KsPipPip_kspip);
630 calt1(pKsPip2,Pip,t1KsPipPip_kspip2);
631 calt1(pKsKm,Pip,t1KsKmPip1_ksk);
632 calt1(pKsKm,Pi0,t1KsKmPip2_ksk);
633
634 calt2(pKsPip1,Kp,t2KsKmPip1_kspip);
635 calt2(pKsPip2,Kp,t2KsKmPip2_kspip);
636 calt2(pKmPip1,Km,t2KsKmPip1_kmpip);
637 calt2(pKmPip2,Km,t2KsKmPip2_kmpip);
638 calt2(pKmPip2,Pip,t2KmPipPip_kmpip);
639 calt2(pKmPip1,Pi0,t2KmPipPip_kmpip2);
640 calt2(pKsPip2,Pip,t2KsPipPip_kspip);
641 calt2(pKsPip1,Pi0,t2KsPipPip_kspip2);
642
643/* calt1(Km,Kp,t1phi);
644 calt1(Pip,Pi0,t1rhop);
645 calt1(Km,Pip,t1kstr1);
646 calt1(Kp,Pi0,t1kstr2);
647 calt1(Km,Pi0,t1kstr3);
648 calt1(Kp,Pi0,t1kstr4);
649 calt1(pKstr0,Pi0,t1A1);
650 calt1(pKstrm,Pip,t1A2);
651
652 calt2(pKstr0,Pi0,t2k11);
653 calt2(pKstrm,Pip,t2k12);
654 calt2(pKstrm,Kp,t2k13);
655 calt2(pKstrp,Km,t2k14);
656 calt2(prhop,Km,t2k21);
657*/
658 double B_phi =-1.0, B_rho =-1.0, B_rho_2 =-1.0, B_phirho1=-1.0, B_phirho2=-1.0, B_phi_2 =-1.0;
659 double B_Kstp =-1.0, B_Kst0 =-1.0, B_KsKs1 =-1.0, B_KsKs2 =-1.0, B_Kstp_2 =-1.0, B_Kst0_2 =-1.0, B_DA2 =-1.0;
660 double B_Kstm =-1.0, B_DA1 =-1.0, B_A1 =-1.0, B_A2 =-1.0, B_A1_1 =-1.0, B_K13_1 =-1.0, B_A2_1 =-1.0;
661 double B_D_Arho =-1.0, B_Arho =-1.0, B_K11 =-1.0, B_D11 =-1.0, B_Arho1 =-1.0, B_K14_phi=-1.0, B_A10 =-1.0;
662 double B_K11rho =-1.0, B_D_K11rho=-1.0, B_K14_1 =-1.0, B_K14_2 =-1.0, B_D_K14 =-1.0, B_K14_km =-1.0, B_K14_kp =-1.0;
663 double B_K_Kst0pi0=-1.0, B_D_KK =-1.0, B_K_Kstmpip=-1.0, B_D_K11_2=-1.0, B_K11_21 =-1.0, B_K11_22 =-1.0, B_A12 =-1.0;
664
665 PDF[0] = 0;
666 PDF[1] = 0;
667 double mass1sq,mass2sq,tt1,tt2,tmp2;
668 double temp_PDF, tmp1, temp[2],tmp3,temp_PDF1;
669 double pro[2], pro0[2], pro1[2],pro2[2],pro3[2],pro4[2];
670 double pro1b[2],pro2b[2],pro_ksk[2];
671 double pro_kspi1[2],pro_kspi2[2],pro_kpi1[2],pro_kpi2[2],pro_kpi2sw[2],pro_kpi1sw[2],pro_kspi2sw[2],pro_kspi1sw[2];
672 double pro_kspi1_0[2],pro_kspi2_0[2],pro_kpi1_0[2],pro_kpi2_0[2];
673 double ispro_kspi1=0,ispro_kspi2=0,ispro_kpi1=0,ispro_kpi2=0,ispro_kpi2sw=0,ispro_kpi1sw=0,ispro_kspi2sw=0,ispro_kspi1sw=0,ispro_ksk=0;
674 double barr_kskpi1_ksk=-1.0,barr_kskpi2_ksk=-1.0,barr_kskpi1_kspi=-1.0,barr_kskpi2_kspi=-1.0,barr_kskpi1_kpi=-1.0,barr_kskpi2_kpi=-1.0,barr_kspi1=-1.0,barr_kspi2=-1.0,barr_kpi1=-1.0,barr_kpi2=-1.0,barr_kspi1_kpi2=-1.0,barr_kspi2_kpi1=-1.0,barr_ds_kskpi1=-1.0,barr_ds_kskpi2=-1.0;
675 double t1A[4],t1V[4],t1D[4],t2D[4][4],B[3];
676 double sKs2[2]={sKs,mass_Pion*mass_Pion};
677 double sKm2[2]={sKm,mass_Eta*mass_Eta};
678 double pS[4], pV[4], pA[4];
679 double flag[3];
680
681 double sa[3], sb[3], sc[3];
682 int isKstp=0.0, isKst0=0.0, isKstm=0.0, isRho=0.0, isf0=0.0, isKmPip_S=0.0, isKpPi0_S=0.0, isKmPi0_S=0.0, isPiPi_S=0.0, isA0980=0.0;
683 double proRho[2], proKstp[2], proKstm[2], proKst0[2], proPiPi_S[2], proKmPip_S[2], proKpPi0_S[2], proKmPi0_S[2], proA0980[2], prof0[2];
684 for(int i=0; i<9; i++){
685 flag[0] = g0[i]; flag[1] = g1[i];flag[2] = g2[i];
686
687 temp_PDF = 0;
688 temp_PDF1 = 0;
689 cof[0] = amp[i]*cos(phase[i]);
690 cof[1] = amp[i]*sin(phase[i]);
691 mass1sq = mass1[i]*mass1[i];//avoid parameters input reversed
692 mass2sq = mass2[i]*mass2[i];
693
694 if(modetype[i]==2)
695 {
696 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
697 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
698 Com_Multi(pro_kspi1,pro_kpi2,pro1b);
699
700 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
701 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
702 if(barr_kspi1_kpi2<0.0) barr_kspi1_kpi2 = barrier(1,sD,sKsPip1,sKmPip2,rD2);
703 if(g2[i]==0)
704 {
705 for(int a=0; a<4; a++)
706 {
707 temp_PDF += G[a][a]*t1KsPip1[a]*t1KmPip2[a];
708 }
709 tmp1 = barr_kspi1*barr_kpi2*temp_PDF;
710 }
711 if(g2[i]==1)
712 {
713 calt1(pKsPip1,pKmPip2,t1D);
714 for(int a=0; a<4; a++)
715 {
716 for(int j=0; j<4; j++)
717 {
718 tt1=G[a][a]*G[j][j]*pD[a]*t1D[j];
719 for(int k=0; k<4; k++)
720 {
721 tt2=t1KsPip1[k]*G[k][k];
722 for(int l=0; l<4; l++)
723 {
724 temp_PDF += E[a][j][k][l]*t1KmPip2[l]*G[l][l]*tt1*tt2;
725 }
726 }
727 }
728 }
729 tmp1 = barr_kspi1*barr_kpi2*barr_kspi1_kpi2*temp_PDF;
730 }
731 if(g2[i]==2)
732 {
733 calt2(pKsPip1,pKmPip2,t2D);
734 for(int a=0; a<4; a++)
735 {
736 for(int j=0; j<4; j++)
737 {
738 temp_PDF += t2D[a][j]*t1KsPip1[a]*t1KmPip2[j]*G[a][a]*G[j][j];
739 }
740 }
741 B[2] = barrier(2,sD,sKsPip1,sKmPip2,rD2);
742 tmp1 = barr_kspi1*barr_kpi2*B[2]*temp_PDF;
743 }
744 temp_PDF=0;
745 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
746 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
747 Com_Multi(pro_kspi2,pro_kpi1,pro2b);
748
749 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
750 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
751 if(barr_kspi2_kpi1<0.0)barr_kspi2_kpi1 = barrier(1,sD,sKsPip2,sKmPip1,rD2);
752 if(g2[i]==0)
753 {
754 for(int a=0; a<4; a++)
755 {
756 temp_PDF += G[a][a]*t1KsPip2[a]*t1KmPip1[a];
757 }
758 tmp2 =barr_kspi2*barr_kpi1*temp_PDF;
759 }
760 if(g2[i]==1)
761 {
762 calt1(pKsPip2,pKmPip1,t1D);
763 for(int a=0; a<4; a++)
764 {
765 for(int j=0; j<4; j++)
766 {
767 tt1 = G[a][a]*G[j][j]*pD[a]*t1D[j];
768 for(int k=0; k<4; k++)
769 {
770 tt2 = t1KsPip2[k]*G[k][k];
771 for(int l=0; l<4; l++)
772 {
773 temp_PDF += E[a][j][k][l]*t1KmPip1[l]*G[l][l]*tt1*tt2;
774 }
775 }
776 }
777 }
778 tmp2 = barr_kspi2*barr_kpi1*barr_kspi2_kpi1*temp_PDF;
779 }
780 if(g2[i]==2)
781 {
782 calt2(pKsPip2,pKmPip1,t2D);
783 for(int a=0; a<4; a++)
784 {
785 for(int j=0; j<4; j++)
786 {
787 temp_PDF += t2D[a][j]*t1KsPip2[a]*t1KmPip1[j]*G[a][a]*G[j][j];
788 }
789 }
790 B[2] = barrier(2,sD,sKsPip2,sKmPip1,rD2);
791 tmp2 = barr_kspi2*barr_kpi1*B[2]*temp_PDF;
792 }
793 }
794
795 //text
796 // K*1410 K*
797 if(modetype[i]==3)
798 {
799 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
800 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
801 Com_Multi(pro_kspi1,pro_kpi2,pro1b);
802
803 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
804 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
805 if(barr_kspi1_kpi2<0.0) barr_kspi1_kpi2 = barrier(1,sD,sKsPip1,sKmPip2,rD2);
806 if(g2[i]==0)
807 {
808 for(int a=0; a<4; a++)
809 {
810 temp_PDF += G[a][a]*t1KsPip1[a]*t1KmPip2[a];
811 }
812 tmp1 = barr_kspi1*barr_kpi2*temp_PDF;
813 }
814 if(g2[i]==1)
815 {
816 calt1(pKsPip1,pKmPip2,t1D);
817 for(int a=0; a<4; a++)
818 {
819 for(int j=0; j<4; j++)
820 {
821 tt1=G[a][a]*G[j][j]*pD[a]*t1D[j];
822 for(int k=0; k<4; k++)
823 {
824 tt2=t1KsPip1[k]*G[k][k];
825 for(int l=0; l<4; l++)
826 {
827 temp_PDF += E[a][j][k][l]*t1KmPip2[l]*G[l][l]*tt1*tt2;
828 }
829 }
830 }
831 }
832 tmp1 = barr_kspi1*barr_kpi2*barr_kspi1_kpi2*temp_PDF;
833 }
834 if(g2[i]==2)
835 {
836 calt2(pKsPip1,pKmPip2,t2D);
837 for(int a=0; a<4; a++)
838 {
839 for(int j=0; j<4; j++)
840 {
841 temp_PDF += t2D[a][j]*t1KsPip1[a]*t1KmPip2[j]*G[a][a]*G[j][j];
842 }
843 }
844 B[2] = barrier(2,sD,sKsPip1,sKmPip2,rD2);
845 tmp1 = barr_kspi1*barr_kpi2*B[2]*temp_PDF;
846 }
847
848 temp_PDF=0;
849 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
850 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
851 Com_Multi(pro_kspi2,pro_kpi1,pro2b);
852
853 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
854 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
855 if(barr_kspi2_kpi1<0.0)barr_kspi2_kpi1 = barrier(1,sD,sKsPip2,sKmPip1,rD2);
856 if(g2[i]==0)
857 {
858 for(int a=0; a<4; a++)
859 {
860 temp_PDF += G[a][a]*t1KsPip2[a]*t1KmPip1[a];
861 }
862 tmp2 =barr_kspi2*barr_kpi1*temp_PDF;
863 }
864 if(g2[i]==1)
865 {
866 calt1(pKsPip2,pKmPip1,t1D);
867 for(int a=0; a<4; a++)
868 {
869 for(int j=0; j<4; j++)
870 {
871 tt1 = G[a][a]*G[j][j]*pD[a]*t1D[j];
872 for(int k=0; k<4; k++)
873 {
874 tt2 = t1KsPip2[k]*G[k][k];
875 for(int l=0; l<4; l++)
876 {
877 temp_PDF += E[a][j][k][l]*t1KmPip1[l]*G[l][l]*tt1*tt2;
878 }
879 }
880 }
881 }
882 tmp2 = barr_kspi2*barr_kpi1*barr_kspi2_kpi1*temp_PDF;
883 }
884 if(g2[i]==2)
885 {
886 calt2(pKsPip2,pKmPip1,t2D);
887 for(int a=0; a<4; a++)
888 {
889 for(int j=0; j<4; j++)
890 {
891 temp_PDF += t2D[a][j]*t1KsPip2[a]*t1KmPip1[j]*G[a][a]*G[j][j];
892 }
893 }
894 B[2] = barrier(2,sD,sKsPip2,sKmPip1,rD2);
895 tmp2 = barr_kspi2*barr_kpi1*B[2]*temp_PDF;
896 }
897 }
898
899 ///Kstp Kmi_S_wave
900 if(modetype[i]==21)
901 {
902 if(ispro_kspi1==0){
903 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
904 ispro_kspi1=1;
905 }
906 if(ispro_kpi2sw==0){
907 KPiSLASS(sKmPip2,sKm,sPip2,pro_kpi2sw);
908 ispro_kpi2sw=1;
909 }
910 Com_Multi(pro_kspi1,pro_kpi2sw,pro1b);
911 calt1(pKsPip1,pKmPip2,t1D);
912
913 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
914 if(barr_kspi1_kpi2<0.0) barr_kspi1_kpi2 = barrier(1,sD,sKsPip1,sKmPip2,rD2);
915 for(int a=0; a<4; a++)
916 {
917 temp_PDF += t1D[a]*t1KsPip1[a]*G[a][a];
918 }
919 tmp1 = barr_kspi1*barr_kspi1_kpi2*temp_PDF;
920
921 temp_PDF=0;
922 if(ispro_kspi2==0) {
923 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
924 ispro_kspi2=1;
925 }
926 if(ispro_kpi1sw==0){
927 KPiSLASS(sKmPip1,sKm,sPip1,pro_kpi1sw);
928 ispro_kpi1sw=1;
929 }
930 Com_Multi(pro_kspi2,pro_kpi1sw,pro2b);
931
932 calt1(pKsPip2,pKmPip1,t1D);
933 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
934 if(barr_kspi2_kpi1<0.0) barr_kspi2_kpi1 = barrier(1,sD,sKsPip2,sKmPip1,rD2);
935 for(int a=0; a<4; a++)
936 {
937 temp_PDF += t1D[a]*t1KsPip2[a]*G[a][a];
938 }
939 tmp2 = barr_kspi2*barr_kspi2_kpi1*temp_PDF;
940 }
941
942 ///Kst0 Kspi_S_wave
943 if(modetype[i]==22)
944 {
945 if(ispro_kpi2==0){
946 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
947 ispro_kpi2=1;
948 }
949 if(ispro_kspi1sw==0){
950 KPiSLASS(sKsPip1,sKs,sPip1,pro_kspi1sw);
951 ispro_kspi1sw=1;
952 }
953 Com_Multi(pro_kpi2,pro_kspi1sw,pro1b);
954
955 calt1(pKsPip1,pKmPip2,t1D);
956 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
957 if(barr_kspi1_kpi2<0.0) barr_kspi1_kpi2 = barrier(1,sD,sKsPip1,sKmPip2,rD2);
958 for(int a=0; a<4; a++)
959 {
960 temp_PDF += t1D[a]*t1KmPip2[a]*G[a][a];
961 }
962 tmp1 = barr_kpi2*barr_kspi1_kpi2*temp_PDF;
963
964 temp_PDF=0;
965 if(ispro_kpi1==0){
966 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
967 ispro_kpi1=1;
968 }
969 if(ispro_kspi2sw==0){
970 KPiSLASS(sKsPip2,sKs,sPip2,pro_kspi2sw);
971 ispro_kspi2sw=1;
972 }
973 Com_Multi(pro_kpi1,pro_kspi2sw,pro2b);
974
975 calt1(pKsPip2,pKmPip1,t1D);
976 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
977 if(barr_kspi2_kpi1<0.0) barr_kspi2_kpi1 = barrier(1,sD,sKsPip2,sKmPip1,rD2);
978 for(int a=0; a<4; a++)
979 {
980 temp_PDF += t1D[a]*t1KmPip1[a]*G[a][a];
981 }
982 tmp2 = barr_kspi2_kpi1*barr_kpi1*temp_PDF;
983 }
984
985
986 //Kspi_swave kpi_swave
987 if(modetype[i]==1)
988 {
989 if(ispro_kspi1sw==0){
990 KPiSLASS(sKsPip1,sKs,sPip1,pro_kspi1sw);
991 ispro_kspi1sw=1;
992 }
993 if(ispro_kpi2sw==0){
994 KPiSLASS(sKmPip2,sKm,sPip2,pro_kpi2sw);
995 ispro_kpi2sw=1;
996 }
997 Com_Multi(pro_kspi1sw,pro_kpi2sw,pro1b);
998 tmp1 = 1;
999
1000 temp_PDF=0;
1001 if(ispro_kspi2sw==0){
1002 KPiSLASS(sKsPip2,sKs,sPip2,pro_kspi2sw);
1003 ispro_kspi2sw=1;
1004 }
1005 if(ispro_kpi1sw==0){
1006 KPiSLASS(sKmPip1,sKm,sPip1,pro_kpi1sw);
1007 ispro_kpi1sw=1;
1008 }
1009 Com_Multi(pro_kspi2sw,pro_kpi1sw,pro2b);
1010 tmp2 = 1;
1011 }
1012
1013
1014
1015 //Ds -> eta(1475)pi+ eta(1475)->K*Ks,K*0>Km pi+
1016 if(modetype[i]==4)
1017 {
1018 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,1,pro1);
1019 if(ispro_kpi1==0){
1020 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
1021 ispro_kpi1=1;
1022 }
1023 Com_Multi(pro1,pro_kpi1, pro1b);
1024
1025 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
1026 barr_kskpi1_kpi = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
1027 for(int a=0; a<4; a++)
1028 {
1029 temp_PDF += Km[a]*t1KmPip1[a]*G[a][a];
1030 }
1031 tmp1 = barr_kpi1*barr_kskpi1_kpi*temp_PDF;
1032
1033 temp_PDF=0;
1034 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro1);
1035 if(ispro_kpi2==0){
1036 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
1037 ispro_kpi2=1;
1038 }
1039 Com_Multi(pro1,pro_kpi2,pro2b);
1040
1041 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
1042 barr_kskpi2_kpi = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
1043 for(int a=0; a<4; a++)
1044 {
1045 temp_PDF += Km[a]*t1KmPip2[a]*G[a][a];
1046 }
1047 tmp2 = barr_kpi2*barr_kskpi2_kpi*temp_PDF;
1048 }
1049
1050 //Ds -> eta(1475)pi+ eta(1475)->K*Km,K*+>Ks pi+
1051 if(modetype[i]==14)
1052 {
1053 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,1,pro1);
1054 if(ispro_kspi1==0){
1055 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
1056 ispro_kspi1=1;
1057 }
1058 Com_Multi(pro1,pro_kspi1,pro1b);
1059
1060 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
1061 if(barr_kskpi1_kspi<0.0)barr_kskpi1_kspi= barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
1062 for(int a=0; a<4; a++)
1063 {
1064 temp_PDF += Kp[a]*t1KsPip1[a]*G[a][a];
1065 }
1066 tmp1 = barr_kspi1*barr_kskpi1_kspi*temp_PDF;
1067
1068 temp_PDF=0;
1069 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,1,pro1);
1070 if(ispro_kspi2==0){
1071 propagatorRBW(mass1sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
1072 ispro_kspi2=1;
1073 }
1074 Com_Multi(pro1,pro_kspi2,pro2b);
1075
1076 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
1077 if(barr_kskpi2_kspi<0.0)barr_kskpi2_kspi= barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
1078 for(int a=0; a<4; a++)
1079 {
1080 temp_PDF += Kp[a]*t1KsPip2[a]*G[a][a];
1081 }
1082 tmp2 = barr_kspi2*barr_kskpi2_kspi*temp_PDF;
1083 }
1084 if(modetype[i]==8)
1085 {
1086 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsKm,sPip1,rRes2,0,pro1);
1087 if(ispro_ksk==0){
1088 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1089 ispro_ksk=1;
1090 }
1091 Com_Multi(pro1,pro_ksk,pro1b);
1092 tmp1=1.0;
1093
1094 temp_PDF=0;
1095 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsKm,sPip2,rRes2,0,pro1);
1096 if(ispro_ksk==0){
1097 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1098 ispro_ksk=1;
1099 }
1100 Com_Multi(pro1,pro_ksk,pro2b);
1101 tmp2=1.0;
1102 }
1103
1104 //f1285->a0 pip
1105 if(modetype[i]==17)
1106 {
1107 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsKm,sPip1,rRes2,1,pro1);
1108 if(ispro_ksk==0){
1109 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1110 ispro_ksk=1;
1111 }
1112 Com_Multi(pro1,pro_ksk,pro1b);
1113 calt1(pKsKmPip1,Pi0,t1D);
1114
1115 if(barr_kskpi1_ksk<0.0) barr_kskpi1_ksk = barrier(1,sKsKmPip1,sKsKm,sPip1,rRes2);
1116 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1117 for(int a=0; a<4; a++)
1118 {
1119 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_ksk[a];
1120 }
1121 tmp1 = barr_kskpi1_ksk*barr_ds_kskpi1*temp_PDF;
1122
1123 temp_PDF=0;
1124 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsKm,sPip2,rRes2,1,pro1);
1125 if(ispro_ksk==0){
1126 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1127 ispro_ksk=1;
1128 }
1129 Com_Multi(pro1,pro_ksk,pro2b);
1130 calt1(pKsKmPip2,Pip,t1D);
1131
1132 if(barr_kskpi2_ksk<0.0)barr_kskpi2_ksk = barrier(1,sKsKmPip2,sKsKm,sPip2,rRes2);
1133 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1134 for(int a=0; a<4; a++)
1135 {
1136 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip2_ksk[a];
1137 }
1138 tmp2 = barr_kskpi2_ksk*barr_ds_kskpi2*temp_PDF;
1139 }
1140
1141 //1+ Ds -> 1550 pi+ , 1550 -> K-K*+ ,K*+ ->Ks pi+
1142 if(modetype[i]==5){
1143 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,g2[i],pro1);
1144 //pro1[0]=1;pro1[1]=0;
1145 if(ispro_kspi1==0){
1146 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
1147 ispro_kspi1=1;
1148 }
1149 Com_Multi(pro1,pro_kspi1,pro1b);
1150 calt1(pKsKmPip1,Pi0,t1D);
1151
1152 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
1153 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1154 barr_kskpi1_kspi= barrier(2,sKsKmPip1,sKsPip1,sKm,rRes2);
1155 if(g2[i]==0)
1156 {
1157 for(int a=0; a<4; a++)
1158 {
1159 for(int j=0; j<4; j++)
1160 {
1161 temp_PDF += t1D[a]*(pKsKmPip1[a]*pKsKmPip1[j]/sKsKmPip1-G[a][j])*t1KsPip1[j]*G[a][a]*G[j][j];
1162 }
1163 }
1164 tmp1 = barr_kspi1*barr_ds_kskpi1*temp_PDF;
1165 }
1166 if(g2[i]==2)
1167 {
1168 for(int a=0; a<4; a++)
1169 {
1170 for(int j=0; j<4; j++)
1171 {
1172 temp_PDF += t1D[a]*t2KsKmPip1_kspip[a][j]*t1KsPip1[j]*G[a][a]*G[j][j];
1173 }
1174 }
1175 tmp1 = barr_kspi1*barr_ds_kskpi1*barr_kskpi1_kspi*temp_PDF;
1176 }
1177
1178 temp_PDF=0;
1179 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,g2[i],pro2);
1180 //pro2[0]=1;pro2[1]=0;
1181 if(ispro_kspi2==0){
1182 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
1183 ispro_kspi2=1;
1184 }
1185 Com_Multi(pro2,pro_kspi2,pro2b);
1186
1187 calt1(pKsKmPip2,Pip,t1D);
1188 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
1189 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1190 barr_kskpi2_kspi = barrier(2,sKsKmPip2,sKsPip2,sKm,rRes2);
1191 if(g2[i]==0)
1192 {
1193 for(int a=0; a<4; a++)
1194 {
1195 for(int j=0; j<4; j++)
1196 {
1197 temp_PDF += t1D[a]*(pKsKmPip2[a]*pKsKmPip2[j]/sKsKmPip2-G[a][j])*t1KsPip2[j]*G[a][a]*G[j][j];
1198 }
1199 }
1200 tmp2 = barr_kspi2*barr_ds_kskpi2*temp_PDF;
1201 }
1202 if(g2[i]==2)
1203 {
1204 for(int a=0; a<4; a++)
1205 {
1206 for(int j=0; j<4; j++)
1207 {
1208 temp_PDF += t1D[a]*t2KsKmPip2_kspip[a][j]*t1KsPip2[j]*G[a][a]*G[j][j];
1209 }
1210 }
1211 tmp2 = barr_kspi2*barr_ds_kskpi2*barr_kskpi2_kspi*temp_PDF;
1212 }
1213 }
1214
1215 //Ds -> 1550 pi+ , 1550 -> KsK*0bar ,K*0bar ->Km pi+
1216 if(modetype[i]==26){
1217 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,g2[i],pro1);
1218 //pro1[0]=1;pro1[1]=0;
1219 if(ispro_kpi1==0){
1220 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
1221 ispro_kpi1=1;
1222 }
1223 Com_Multi(pro1,pro_kpi1, pro1b);
1224 calt1(pKsKmPip1,Pi0,t1D);
1225
1226 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
1227 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1228 barr_kskpi1_kpi = barrier(2,sKsKmPip1,sKmPip1,sKs,rRes2);
1229 if(g2[i]==0)
1230 {
1231 for(int a=0; a<4; a++)
1232 {
1233 for(int j=0; j<4; j++)
1234 {
1235 temp_PDF += t1D[a]*(pKsKmPip1[a]*pKsKmPip1[j]/sKsKmPip1-G[a][j])*t1KmPip1[j]*G[a][a]*G[j][j];
1236 }
1237 }
1238 tmp1 = barr_kpi1*barr_ds_kskpi1*temp_PDF;
1239 }
1240 if(g2[i]==2)
1241 {
1242 for(int a=0; a<4; a++)
1243 {
1244 for(int j=0; j<4; j++)
1245 {
1246 temp_PDF += t1D[a]*t2KsKmPip1_kmpip[a][j]*t1KmPip1[j]*G[a][a]*G[j][j];
1247 }
1248 }
1249 tmp1 = barr_kpi1*barr_ds_kskpi1*barr_kskpi1_kpi*temp_PDF;
1250 }
1251
1252
1253 temp_PDF=0;
1254 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,g2[i],pro2);
1255 //pro2[0]=1;pro2[1]=0;
1256 if(ispro_kpi2==0){
1257 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
1258 ispro_kpi2=1;
1259 }
1260 Com_Multi(pro2,pro_kpi2, pro2b);
1261
1262 calt1(pKsKmPip2,Pip,t1D);
1263 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
1264 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1265 barr_kskpi2_kpi = barrier(2,sKsKmPip2,sKmPip2,sKs,rRes2);
1266 if(g2[i]==0)
1267 {
1268 for(int a=0; a<4; a++)
1269 {
1270 for(int j=0; j<4; j++)
1271 {
1272 temp_PDF += t1D[a]*(pKsKmPip2[a]*pKsKmPip2[j]/sKsKmPip2-G[a][j])*t1KmPip2[j]*G[a][a]*G[j][j];
1273 }
1274
1275 }
1276 tmp2 = barr_kpi2*barr_ds_kskpi2*temp_PDF;
1277 }
1278 if(g2[i]==2)
1279 {
1280 for(int a=0; a<4; a++)
1281 {
1282 for(int j=0; j<4; j++)
1283 {
1284 temp_PDF += t1D[a]*t2KsKmPip2_kmpip[a][j]*t1KmPip2[j]*G[a][a]*G[j][j];
1285 }
1286 }
1287 tmp2 = barr_kpi2*barr_ds_kskpi2*barr_kskpi2_kpi*temp_PDF;
1288 }
1289 }
1290
1291
1292
1293
1294 //text
1295 //Ds -> eta(1405)pi+ eta(1405)->K*Ks,K*0>Km pi+
1296 if(modetype[i]==7)
1297 {
1298 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,1,pro1);
1299 if(ispro_kpi1==0){
1300 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
1301 ispro_kpi1=1;
1302 }
1303 Com_Multi(pro1,pro_kpi1, pro1b);
1304
1305 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
1306 barr_kskpi1_kpi = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
1307 for(int a=0; a<4; a++)
1308 {
1309 temp_PDF += Km[a]*t1KmPip1[a]*G[a][a];
1310 }
1311 tmp1 = barr_kpi1*barr_kskpi1_kpi*temp_PDF;
1312
1313 temp_PDF=0;
1314 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro1);
1315 if(ispro_kpi2==0){
1316 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
1317 ispro_kpi2=1;
1318 }
1319 Com_Multi(pro1,pro_kpi2,pro2b);
1320
1321 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
1322 barr_kskpi2_kpi = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
1323 for(int a=0; a<4; a++)
1324 {
1325 temp_PDF += Km[a]*t1KmPip2[a]*G[a][a];
1326 }
1327 tmp2 = barr_kpi2*barr_kskpi2_kpi*temp_PDF;
1328 }
1329
1330 //text
1331 //Ds -> eta(1405)pi+ eta(1405)->K*Km,K*+>Ks pi+
1332 if(modetype[i]==10)
1333 {
1334 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,1,pro1);
1335 if(ispro_kspi1==0){
1336 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
1337 ispro_kspi1=1;
1338 }
1339 Com_Multi(pro1,pro_kspi1,pro1b);
1340
1341 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
1342 if(barr_kskpi1_kspi<0.0)barr_kskpi1_kspi= barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
1343 for(int a=0; a<4; a++)
1344 {
1345 temp_PDF += Kp[a]*t1KsPip1[a]*G[a][a];
1346 }
1347 tmp1 = barr_kspi1*barr_kskpi1_kspi*temp_PDF;
1348
1349 temp_PDF=0;
1350 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,1,pro1);
1351 if(ispro_kspi2==0){
1352 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
1353 ispro_kspi2=1;
1354 }
1355 Com_Multi(pro1,pro_kspi2,pro2b);
1356
1357 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
1358 if(barr_kskpi2_kspi<0.0)barr_kskpi2_kspi= barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
1359 for(int a=0; a<4; a++)
1360 {
1361 temp_PDF += Kp[a]*t1KsPip2[a]*G[a][a];
1362 }
1363 tmp2 = barr_kspi2*barr_kskpi2_kspi*temp_PDF;
1364 }
1365
1366 // text
1367 // eta1405->a0 pip
1368 if(modetype[i]==9)
1369 {
1370 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsKm,sPip1,rRes2,0,pro1);
1371 if(ispro_ksk==0){
1372 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1373 ispro_ksk=1;
1374 }
1375 Com_Multi(pro1,pro_ksk,pro1b);
1376 tmp1=1.0;
1377
1378 temp_PDF=0;
1379 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsKm,sPip2,rRes2,0,pro1);
1380 if(ispro_ksk==0){
1381 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1382 ispro_ksk=1;
1383 }
1384 Com_Multi(pro1,pro_ksk,pro2b);
1385 tmp2=1.0;
1386 }
1387
1388 //text
1389 //Ds -> eta(1405)pi+ eta(1405)->kpisw-ks
1390 if(modetype[i]==133)
1391 {
1392 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,0,pro1);
1393 if(ispro_kpi1sw==0){
1394 KPiSLASS(sKmPip1,sKm,sPip1,pro_kpi1sw);
1395 ispro_kpi1sw=1;
1396 }
1397 Com_Multi(pro1,pro_kpi1sw, pro1b);
1398 tmp1= 1.0;
1399
1400 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,0,pro1);
1401 if(ispro_kpi2sw==0){
1402 KPiSLASS(sKmPip2,sKm,sPip2,pro_kpi2sw);
1403 ispro_kpi2sw=1;
1404 }
1405 Com_Multi(pro1,pro_kpi2sw, pro2b);
1406 tmp2= 1.0;
1407 }
1408
1409 //Ds -> eta(1405)pi+ eta(1405)->kpisw-ks
1410 if(modetype[i]==13)
1411 {
1412 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,0,pro1);
1413 if(ispro_kpi1sw==0){
1414 KPiSLASS(sKmPip1,sKm,sPip1,pro_kpi1sw);
1415 ispro_kpi1sw=1;
1416 }
1417 Com_Multi(pro1,pro_kpi1sw, pro1b);
1418 tmp1= 1.0;
1419
1420 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,0,pro1);
1421 if(ispro_kpi2sw==0){
1422 KPiSLASS(sKmPip2,sKm,sPip2,pro_kpi2sw);
1423 ispro_kpi2sw=1;
1424 }
1425 Com_Multi(pro1,pro_kpi2sw, pro2b);
1426 tmp2= 1.0;
1427 }
1428
1429 //text
1430 //Ds -> eta(1405)pi+ eta(1405)->kspisw-k
1431 if(modetype[i]==20)
1432 {
1433 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,0,pro1);
1434 if(ispro_kspi1sw==0){
1435 KPiSLASS(sKsPip1,sKs,sPip1,pro_kspi1sw);
1436 ispro_kspi1sw=1;
1437 }
1438 Com_Multi(pro1,pro_kspi1sw,pro1b);
1439 tmp1= 1.0;
1440
1441 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,0,pro1);
1442 if(ispro_kspi2sw==0){
1443 KPiSLASS(sKsPip2,sKs,sPip2,pro_kspi2sw);
1444 ispro_kspi2sw=1;
1445 }
1446 Com_Multi(pro1,pro_kspi2sw,pro2b);
1447 tmp2= 1.0;
1448 }
1449
1450 //text
1451 //Ds -> eta(1475)pi+ eta(1475)->kpisw-kspisw
1452 if(modetype[i]==19)
1453 {
1454 if(ispro_kspi1sw==0){
1455 KPiSLASS(sKsPip1,sKs,sPip1,pro_kspi1sw);
1456 ispro_kspi1sw=1;
1457 }
1458 if(ispro_kpi1sw==0){
1459 KPiSLASS(sKmPip1,sKm,sPip1,pro_kpi1sw);
1460 ispro_kpi1sw=1;
1461 }
1462 Com_Multi(pro_kpi1sw,pro_kspi1sw,pro1b);
1463 tmp1= 1.0;
1464
1465 if(ispro_kspi2sw==0){
1466 KPiSLASS(sKsPip2,sKs,sPip2,pro_kspi2sw);
1467 ispro_kspi2sw=1;
1468 }
1469 if(ispro_kpi2sw==0){
1470 KPiSLASS(sKmPip2,sKm,sPip2,pro_kpi2sw);
1471 ispro_kpi2sw=1;
1472 }
1473 Com_Multi(pro_kpi2sw,pro_kspi2sw,pro2b);
1474 tmp2= 1.0;
1475 }
1476
1477 //text
1478 // Ds -> f1420 pi+ , f1420 -> K-K*+ ,K*+ ->Ks pi+
1479 if(modetype[i]==11){
1480 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,g2[i],pro1);
1481 if(ispro_kspi1==0){
1482 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
1483 ispro_kspi1=1;
1484 }
1485 Com_Multi(pro1,pro_kspi1,pro1b);
1486 calt1(pKsKmPip1,Pi0,t1D);
1487
1488 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
1489 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1490 barr_kskpi1_kspi= barrier(2,sKsKmPip1,sKsPip1,sKm,rRes2);
1491 if(g2[i]==0)
1492 {
1493 for(int a=0; a<4; a++)
1494 {
1495 for(int j=0; j<4; j++)
1496 {
1497 temp_PDF += t1D[a]*(pKsKmPip1[a]*pKsKmPip1[j]/sKsKmPip1-G[a][j])*t1KsPip1[j]*G[a][a]*G[j][j];
1498 }
1499 }
1500 tmp1 = barr_kspi1*barr_ds_kskpi1*temp_PDF;
1501 }
1502 if(g2[i]==2)
1503 {
1504 for(int a=0; a<4; a++)
1505 {
1506 for(int j=0; j<4; j++)
1507 {
1508 temp_PDF += t1D[a]*t2KsKmPip1_kspip[a][j]*t1KsPip1[j]*G[a][a]*G[j][j];
1509 }
1510 }
1511 tmp1 = barr_kspi1*barr_ds_kskpi1*barr_kskpi1_kspi*temp_PDF;
1512 }
1513
1514
1515 temp_PDF=0;
1516 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,g2[i],pro1);
1517 if(ispro_kspi2==0){
1518 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
1519 ispro_kspi2=1;
1520 }
1521 Com_Multi(pro1,pro_kspi2,pro2b);
1522
1523 calt1(pKsKmPip2,Pip,t1D);
1524 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
1525 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1526 barr_kskpi2_kspi = barrier(2,sKsKmPip2,sKsPip2,sKm,rRes2);
1527 if(g2[i]==0)
1528 {
1529 for(int a=0; a<4; a++)
1530 {
1531 for(int j=0; j<4; j++)
1532 {
1533 temp_PDF += t1D[a]*(pKsKmPip2[a]*pKsKmPip2[j]/sKsKmPip2-G[a][j])*t1KsPip2[j]*G[a][a]*G[j][j];
1534 }
1535 }
1536 tmp2 = barr_kspi2*barr_ds_kskpi2*temp_PDF;
1537 }
1538 if(g2[i]==2)
1539 {
1540 for(int a=0; a<4; a++)
1541 {
1542 for(int j=0; j<4; j++)
1543 {
1544 temp_PDF += t1D[a]*t2KsKmPip2_kspip[a][j]*t1KsPip2[j]*G[a][a]*G[j][j];
1545 }
1546 }
1547 tmp2 = barr_kspi2*barr_ds_kskpi2*barr_kskpi2_kspi*temp_PDF;
1548 }
1549 }
1550 if(modetype[i]==1200){
1551 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,g2[i],pro1);
1552
1553 pro1[0]=1;
1554 pro1[1]=0;
1555 if(ispro_kpi1==0){
1556 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
1557 ispro_kpi1=1;
1558 }
1559 Com_Multi(pro1,pro_kpi1, pro1b);
1560 calt1(pKsKmPip1,Pi0,t1D);
1561
1562 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
1563 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1564 barr_kskpi1_kpi = barrier(2,sKsKmPip1,sKmPip1,sKs,rRes2);
1565 if(g2[i]==0)
1566 {
1567 for(int a=0; a<4; a++)
1568 {
1569 for(int j=0; j<4; j++)
1570 {
1571 temp_PDF += t1D[a]*(pKsKmPip1[a]*pKsKmPip1[j]/sKsKmPip1-G[a][j])*t1KmPip1[j]*G[a][a]*G[j][j]; }
1572 }
1573 tmp1 = barr_kpi1*barr_ds_kskpi1*temp_PDF;
1574 }
1575 if(g2[i]==2)
1576 {
1577 for(int a=0; i<4; i++)
1578 {
1579 for(int j=0; j<4; j++)
1580 {
1581 temp_PDF += t1D[a]*t2KsKmPip1_kmpip[a][j]*t1KmPip1[j]*G[a][a]*G[j][j];
1582 }
1583 }
1584 tmp1 = barr_kpi1*barr_ds_kskpi1*barr_kskpi1_kpi*temp_PDF;
1585 }
1586
1587 temp_PDF=0;
1588 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,g2[i],pro2);
1589 pro2[0]=1;
1590 pro2[1]=0;
1591 if(ispro_kpi2==0){
1592 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
1593 ispro_kpi2=1;
1594 }
1595 Com_Multi(pro2,pro_kpi2, pro2b);
1596
1597 calt1(pKsKmPip2,Pip,t1D);
1598 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
1599 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1600 barr_kskpi2_kpi = barrier(2,sKsKmPip2,sKmPip2,sKs,rRes2);
1601 if(g2[i]==0)
1602 {
1603 for(int a=0; a<4; a++)
1604 {
1605 for(int j=0; j<4; j++)
1606 {
1607 temp_PDF += t1D[a]*(pKsKmPip2[a]*pKsKmPip2[j]/sKsKmPip2-G[a][j])*t1KmPip2[j]*G[a][a]*G[j][j];
1608 }
1609 }
1610 tmp2 = barr_kpi2*barr_ds_kskpi2*temp_PDF;
1611 }
1612 if(g2[i]==2)
1613 {
1614 for(int a=0; a<4; a++)
1615 {
1616 for(int j=0; j<4; j++)
1617 {
1618 temp_PDF += t1D[a]*t2KsKmPip2_kmpip[a][j]*t1KmPip2[j]*G[a][a]*G[j][j];
1619 }
1620 }
1621 tmp2 = barr_kpi2*barr_ds_kskpi2*barr_kskpi2_kpi*temp_PDF;
1622 }
1623 }
1624 //text
1625 // Ds -> f1420(nonResonance) pi+ , f1420 -> K-K*+ ,K*+ ->Ks pi+
1626 if(modetype[i]==1100){
1627 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,g2[i],pro1);
1628 pro1[0]=1;
1629 pro1[1]=0;
1630 if(ispro_kspi1==0){
1631 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
1632 ispro_kspi1=1;
1633 }
1634 Com_Multi(pro1,pro_kspi1,pro1b);
1635 calt1(pKsKmPip1,Pi0,t1D);
1636
1637 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
1638 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1639 barr_kskpi1_kspi= barrier(2,sKsKmPip1,sKsPip1,sKm,rRes2);
1640 if(g2[i]==0)
1641 {
1642 for(int a=0; a<4; a++)
1643 {
1644 for(int j=0; j<4; j++)
1645 {
1646 temp_PDF += t1D[a]*(pKsKmPip1[a]*pKsKmPip1[j]/sKsKmPip1-G[a][j])*t1KsPip1[j]*G[a][a]*G[j][j];
1647 }
1648 }
1649 tmp1 = barr_kspi1*barr_ds_kskpi1*temp_PDF;
1650 }
1651 if(g2[i]==2)
1652 {
1653 for(int a=0; a<4; a++)
1654 {
1655 for(int j=0; j<4; j++)
1656 {
1657 temp_PDF += t1D[a]*t2KsKmPip1_kspip[a][j]*t1KsPip1[j]*G[a][a]*G[j][j];
1658 }
1659 }
1660 tmp1 = barr_kspi1*barr_ds_kskpi1*barr_kskpi1_kspi*temp_PDF;
1661 }
1662
1663
1664 temp_PDF=0;
1665 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,g2[i],pro1);
1666 pro1[0]=1;
1667 pro1[1]=0;
1668 if(ispro_kspi2==0){
1669 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
1670 ispro_kspi2=1;
1671 }
1672 Com_Multi(pro1,pro_kspi2,pro2b);
1673
1674 calt1(pKsKmPip2,Pip,t1D);
1675 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
1676 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1677 barr_kskpi2_kspi = barrier(2,sKsKmPip2,sKsPip2,sKm,rRes2);
1678 if(g2[i]==0)
1679 {
1680 for(int a=0; a<4; a++)
1681 {
1682 for(int j=0; j<4; j++)
1683 {
1684 temp_PDF += t1D[a]*(pKsKmPip2[a]*pKsKmPip2[j]/sKsKmPip2-G[a][j])*t1KsPip2[j]*G[a][a]*G[j][j];
1685 }
1686 }
1687 tmp2 = barr_kspi2*barr_ds_kskpi2*temp_PDF;
1688 }
1689 if(g2[i]==2)
1690 {
1691 for(int a=0; a<4; a++)
1692 {
1693 for(int j=0; j<4; j++)
1694 {
1695 temp_PDF += t1D[a]*t2KsKmPip2_kspip[a][j]*t1KsPip2[j]*G[a][a]*G[j][j];
1696 }
1697 }
1698 tmp2 = barr_kspi2*barr_ds_kskpi2*barr_kskpi2_kspi*temp_PDF;
1699 }
1700 }
1701
1702
1703 //text
1704 //Ds -> 1420 pi+ , 1420 -> KsK*0bar ,K*0bar ->Km pi+
1705 if(modetype[i]==12){
1706 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,g2[i],pro1);
1707 if(ispro_kpi1==0){
1708 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
1709 ispro_kpi1=1;
1710 }
1711 Com_Multi(pro1,pro_kpi1, pro1b);
1712 calt1(pKsKmPip1,Pi0,t1D);
1713
1714 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
1715 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1716 barr_kskpi1_kpi = barrier(2,sKsKmPip1,sKmPip1,sKs,rRes2);
1717 if(g2[i]==0)
1718 {
1719 for(int a=0; a<4; a++)
1720 {
1721 for(int j=0; j<4; j++)
1722 {
1723 temp_PDF += t1D[a]*(pKsKmPip1[a]*pKsKmPip1[j]/sKsKmPip1-G[a][j])*t1KmPip1[j]*G[a][a]*G[j][j]; }
1724 }
1725 tmp1 = barr_kpi1*barr_ds_kskpi1*temp_PDF;
1726 }
1727 if(g2[i]==2)
1728 {
1729 for(int a=0; i<4; i++)
1730 {
1731 for(int j=0; j<4; j++)
1732 {
1733 temp_PDF += t1D[a]*t2KsKmPip1_kmpip[a][j]*t1KmPip1[j]*G[a][a]*G[j][j];
1734 }
1735 }
1736 tmp1 = barr_kpi1*barr_ds_kskpi1*barr_kskpi1_kpi*temp_PDF;
1737 }
1738 temp_PDF=0;
1739 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,g2[i],pro2);
1740 if(ispro_kpi2==0){
1741 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
1742 ispro_kpi2=1;
1743 }
1744 Com_Multi(pro2,pro_kpi2, pro2b);
1745
1746 calt1(pKsKmPip2,Pip,t1D);
1747 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
1748 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1749 barr_kskpi2_kpi = barrier(2,sKsKmPip2,sKmPip2,sKs,rRes2);
1750 if(g2[i]==0)
1751 {
1752 for(int a=0; a<4; a++)
1753 {
1754 for(int j=0; j<4; j++)
1755 {
1756 temp_PDF += t1D[a]*(pKsKmPip2[a]*pKsKmPip2[j]/sKsKmPip2-G[a][j])*t1KmPip2[j]*G[a][a]*G[j][j];
1757 }
1758 }
1759 tmp2 = barr_kpi2*barr_ds_kskpi2*temp_PDF;
1760 }
1761 if(g2[i]==2)
1762 {
1763 for(int a=0; a<4; a++)
1764 {
1765 for(int j=0; j<4; j++)
1766 {
1767 temp_PDF += t1D[a]*t2KsKmPip2_kmpip[a][j]*t1KmPip2[j]*G[a][a]*G[j][j];
1768 }
1769 }
1770 tmp2 = barr_kpi2*barr_ds_kskpi2*barr_kskpi2_kpi*temp_PDF;
1771 }
1772 }
1773 if(modetype[i]==1500){
1774 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsKm,sPip1,rRes2,1,pro1);
1775 pro1[0]=1;
1776 pro1[1]=0;
1777 if(ispro_ksk==0){
1778 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1779 ispro_ksk=1;
1780 }
1781 Com_Multi(pro1,pro_ksk,pro1b);
1782 calt1(pKsKmPip1,Pi0,t1D);
1783
1784 if(barr_kskpi1_ksk<0.0) barr_kskpi1_ksk = barrier(1,sKsKmPip1,sKsKm,sPip1,rRes2);
1785 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1786 for(int a=0; a<4; a++)
1787 {
1788 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_ksk[a];
1789 }
1790 tmp1 = barr_kskpi1_ksk*barr_ds_kskpi1*temp_PDF;
1791
1792 temp_PDF=0;
1793 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsKm,sPip2,rRes2,1,pro1);
1794 pro1[0]=1;
1795 pro1[1]=0;
1796 if(ispro_ksk==0){
1797 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1798 ispro_ksk=1;
1799 }
1800 Com_Multi(pro1,pro_ksk,pro2b);
1801 calt1(pKsKmPip2,Pip,t1D);
1802
1803 if(barr_kskpi2_ksk<0.0) barr_kskpi2_ksk = barrier(1,sKsKmPip2,sKsKm,sPip2,rRes2);
1804 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1805 for(int a=0; a<4; a++)
1806 {
1807 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip2_ksk[a];
1808 }
1809 tmp2 = barr_kskpi2_ksk*barr_ds_kskpi2*temp_PDF;
1810 }
1811
1812
1813 //text
1814 /// f1420->a0 pip
1815 if(modetype[i]==15){
1816 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsKm,sPip1,rRes2,1,pro1);
1817 if(ispro_ksk==0){
1818 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1819 ispro_ksk=1;
1820 }
1821 Com_Multi(pro1,pro_ksk,pro1b);
1822 calt1(pKsKmPip1,Pi0,t1D);
1823
1824 if(barr_kskpi1_ksk<0.0) barr_kskpi1_ksk = barrier(1,sKsKmPip1,sKsKm,sPip1,rRes2);
1825 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1826 for(int a=0; a<4; a++)
1827 {
1828 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_ksk[a];
1829 }
1830 tmp1 = barr_kskpi1_ksk*barr_ds_kskpi1*temp_PDF;
1831
1832 temp_PDF=0;
1833 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsKm,sPip2,rRes2,1,pro1);
1834 if(ispro_ksk==0){
1835 propagatora0980(mass2sq,mass2[i],sKsKm,sKs2,sKm2,pro_ksk);
1836 ispro_ksk=1;
1837 }
1838 Com_Multi(pro1,pro_ksk,pro2b);
1839 calt1(pKsKmPip2,Pip,t1D);
1840
1841 if(barr_kskpi2_ksk<0.0) barr_kskpi2_ksk = barrier(1,sKsKmPip2,sKsKm,sPip2,rRes2);
1842 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1843 for(int a=0; a<4; a++)
1844 {
1845 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip2_ksk[a];
1846 }
1847 tmp2 = barr_kskpi2_ksk*barr_ds_kskpi2*temp_PDF;
1848 }
1849
1850 //text
1851 /// f1420->kspi_sw
1852 if(modetype[i]==24){
1853 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,1,pro1);
1854 if(ispro_kspi1sw==0){
1855 KPiSLASS(sKsPip1,sKs,sPip1,pro_kspi1sw);
1856 ispro_kspi1sw=1;
1857 }
1858 Com_Multi(pro1,pro_kspi1sw,pro1b);
1859 calt1(pKsKmPip1,Pi0,t1D);
1860
1861 if(barr_kskpi1_kspi<0.0) barr_kskpi1_kspi = barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
1862 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1863 for(int a=0; a<4; a++)
1864 {
1865 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kspip[a];
1866 }
1867 tmp1 = barr_kskpi1_kspi*barr_ds_kskpi1*temp_PDF;
1868
1869 temp_PDF=0;
1870 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,1,pro1);
1871 if(ispro_kspi2sw==0){
1872 KPiSLASS(sKsPip2,sKs,sPip2,pro_kspi2sw);
1873 ispro_kspi2sw=1;
1874 }
1875 Com_Multi(pro1,pro_kspi2sw,pro2b);
1876 calt1(pKsKmPip2,Pip,t1D);
1877
1878 if(barr_kskpi2_kspi<0.0) barr_kskpi2_kspi = barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
1879 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1880 for(int a=0; a<4; a++)
1881 {
1882 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip2_kspip[a];
1883 }
1884 tmp2 = barr_kskpi2_kspi*barr_ds_kskpi2*temp_PDF;
1885 }
1886
1887
1888 //text
1889 /// f1420->kpi_sw
1890 if(modetype[i]==27){
1891 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,1,pro1);
1892 if(ispro_kspi1sw==0){
1893 KPiSLASS(sKsPip1,sKm,sPip1,pro_kpi1sw);
1894 ispro_kpi1sw=1;
1895 }
1896 Com_Multi(pro1,pro_kpi1sw,pro1b);
1897 calt1(pKsKmPip1,Pi0,t1D);
1898
1899 barr_kskpi1_kpi = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
1900 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1901 for(int a=0; a<4; a++)
1902 {
1903 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kmpip[a];
1904 }
1905 tmp1 = barr_kskpi1_kpi*barr_ds_kskpi1*temp_PDF;
1906
1907 temp_PDF=0;
1908 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro1);
1909 if(ispro_kpi2sw==0){
1910 KPiSLASS(sKsPip2,sKs,sPip2,pro_kpi2sw);
1911 ispro_kpi2sw=1;
1912 }
1913 Com_Multi(pro1,pro_kpi2sw,pro2b);
1914 calt1(pKsKmPip2,Pip,t1D);
1915
1916 barr_kskpi2_kpi = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
1917 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1918 for(int a=0; a<4; a++)
1919 {
1920 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip2_kmpip[a];
1921 }
1922 tmp2 = barr_kskpi2_kpi*barr_ds_kskpi2*temp_PDF;
1923 }
1924
1925 //Ds -> eta(1475)pi+ eta(1475)->kappa Ks ,K*->Km pi+
1926 if(modetype[i]==244)
1927 {
1928 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,0,pro1);
1929 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,0,pro3);
1930 Com_Multi(pro1,pro3,pro1b);
1931 tmp1= 1.0;
1932
1933 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,0,pro1);
1934 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,0,pro3);
1935 Com_Multi(pro1,pro3,pro2b);
1936 tmp2= 1.0;
1937 }
1938
1939 // Ds -> kst(1410)K- , phi(1680)->K-K*+ ,K*+ ->Ks pi+
1940 if(modetype[i]==280)
1941 {
1942 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPipPip,sKsPip1,sPip2,rRes2,1,pro1);
1943 if (ispro_kspi1==0){
1944 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
1945 ispro_kspi1=1;
1946 }
1947 Com_Multi(pro1,pro_kspi1,pro1b);
1948 //printf(" %.15f %.15f %.15f\n",B[0],B[1],B[2]);
1949
1950 B[0] = barrier(1,sKsPip1,sKs,sPip1,rRes2);
1951 B[1] = barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
1952 B[2] = barrier(1,sD,sKsKmPip1,sPip2,rD2);
1953 for(int a=0; a<4; a++)
1954 {
1955 for(int j=0; j<4; j++)
1956 {
1957 for(int k=0; k<4; k++)
1958 {
1959 for(int l=0; l<4; l++)
1960 {
1961 temp_PDF += E[a][j][k][l]*pKsKmPip1[a]*(pKsPip1[j]-Kp[j])*Pi0[k]*(Km[l]-Pip[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
1962 }
1963 }
1964 }
1965 }
1966 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
1967
1968
1969 temp_PDF=0;
1970 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPipPip,sKsPip2,sPip1,rRes2,1,pro2);
1971 if (ispro_kspi2==0){
1972 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
1973 ispro_kspi2=1;
1974 }
1975 Com_Multi(pro2,pro_kspi2,pro2b);
1976
1977
1978 B[0] = barrier(1,sKsPip2,sKs,sPip2,rRes2);
1979 B[1] = barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
1980 B[2] = barrier(1,sD,sKsKmPip2,sPip1,rD2);
1981
1982 for(int a=0; a<4; a++)
1983 {
1984 for(int j=0; j<4; j++)
1985 {
1986 for(int k=0; k<4; k++)
1987 {
1988 for(int l=0; l<4; l++)
1989 {
1990 temp_PDF += E[a][j][k][l]*pKsKmPip2[a]*(pKsPip2[j]-Kp[j])*Pip[k]*(Km[l]-Pi0[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
1991 }
1992 }
1993 }
1994 }
1995 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
1996 }
1997
1998 //text
1999 // Ds -> kst(1410)KS0 , phi(1680)->K-K*+ ,K*+ ->Km pi+
2000 if(modetype[i]==290){
2001 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPipPip,sKmPip1,sPip2,rRes2,1,pro1);
2002 if (ispro_kpi1==0){
2003 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
2004 ispro_kpi1=1;
2005 }
2006 Com_Multi(pro1,pro_kpi1,pro1b);
2007
2008 B[0] = barrier(1,sKmPip1,sKm,sPip1,rRes2);
2009 B[1] = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
2010 B[2] = barrier(1,sD,sKsKmPip1,sPip2,rD2);
2011 for(int a=0; a<4; a++)
2012 {
2013 for(int j=0; j<4; j++)
2014 {
2015 for(int k=0; k<4; k++)
2016 {
2017 for(int l=0; l<4; l++)
2018 {
2019 temp_PDF += E[a][j][k][l]*pKsKmPip1[a]*(pKmPip1[j]-Km[j])*Pi0[k]*(Kp[l]-Pip[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2020 }
2021 }
2022 }
2023 }
2024 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2025
2026
2027 temp_PDF=0;
2028 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPipPip,sKmPip2,sPip1,rRes2,1,pro2);
2029 if (ispro_kpi2==0){
2030 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
2031 ispro_kpi2=1;
2032 }
2033 Com_Multi(pro2,pro_kpi2,pro2b);
2034 B[0] = barrier(1,sKmPip2,sKm,sPip2,rRes2);
2035 B[1] = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
2036 B[2] = barrier(1,sD,sKsKmPip2,sPip1,rD2);
2037
2038 for(int a=0; a<4; a++)
2039 {
2040 for(int j=0; j<4; j++)
2041 {
2042 for(int k=0; k<4; k++)
2043 {
2044 for(int l=0; l<4; l++)
2045 {
2046 temp_PDF += E[a][j][k][l]*pKsKmPip2[a]*(pKmPip2[j]-Km[j])*Pip[k]*(Kp[l]-Pi0[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2047 }
2048 }
2049 }
2050 }
2051
2052 tmp2 = B[0]*B[1]*B[2]*temp_PDF*(-1);
2053 }
2054 // text
2055 // Ds -> phi(1680)pi+ , phi(1680)->K-K*+ ,K*+ ->Ks pi+
2056 if(modetype[i]==28)
2057 {
2058 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,1,pro1);
2059 if (ispro_kspi1==0){
2060 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
2061 ispro_kspi1=1;
2062 }
2063 Com_Multi(pro1,pro_kspi1,pro1b);
2064 //printf(" %.15f %.15f %.15f\n",B[0],B[1],B[2]);
2065
2066 B[0] = barrier(1,sKsPip1,sKs,sPip1,rRes2);
2067 B[1] = barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
2068 B[2] = barrier(1,sD,sKsKmPip1,sPip2,rD2);
2069 for(int a=0; a<4; a++)
2070 {
2071 for(int j=0; j<4; j++)
2072 {
2073 for(int k=0; k<4; k++)
2074 {
2075 for(int l=0; l<4; l++)
2076 {
2077 temp_PDF += E[a][j][k][l]*pKsKmPip1[a]*(pKsPip1[j]-Kp[j])*Pi0[k]*(Km[l]-Pip[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2078 }
2079 }
2080 }
2081 }
2082 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2083 temp_PDF=0;
2084 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,1,pro2);
2085 if (ispro_kspi2==0){
2086 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
2087 ispro_kspi2=1;
2088 }
2089 Com_Multi(pro2,pro_kspi2,pro2b);
2090
2091
2092 B[0] = barrier(1,sKsPip2,sKs,sPip2,rRes2);
2093 B[1] = barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
2094 B[2] = barrier(1,sD,sKsKmPip2,sPip1,rD2);
2095
2096 for(int a=0; a<4; a++)
2097 {
2098 for(int j=0; j<4; j++)
2099 {
2100 for(int k=0; k<4; k++)
2101 {
2102 for(int l=0; l<4; l++)
2103 {
2104 temp_PDF += E[a][j][k][l]*pKsKmPip2[a]*(pKsPip2[j]-Kp[j])*Pip[k]*(Km[l]-Pi0[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2105 }
2106 }
2107 }
2108 }
2109 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
2110 }
2111
2112 //text
2113 // Ds -> phi(1680)pi+ , phi(1680)->K-K*+ ,K*+ ->Km pi+
2114 if(modetype[i]==29){
2115 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKm,rRes2,1,pro1);
2116 if (ispro_kpi1==0){
2117 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
2118 ispro_kpi1=1;
2119 }
2120 Com_Multi(pro1,pro_kpi1,pro1b);
2121
2122 B[0] = barrier(1,sKmPip1,sKm,sPip1,rRes2);
2123 B[1] = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
2124 B[2] = barrier(1,sD,sKsKmPip1,sPip2,rD2);
2125 for(int a=0; a<4; a++)
2126 {
2127 for(int j=0; j<4; j++)
2128 {
2129 for(int k=0; k<4; k++)
2130 {
2131 for(int l=0; l<4; l++)
2132 {
2133 temp_PDF += E[a][j][k][l]*pKsKmPip1[a]*(pKmPip1[j]-Km[j])*Pi0[k]*(Kp[l]-Pip[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2134 }
2135 }
2136 }
2137 }
2138 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2139
2140
2141 temp_PDF=0;
2142 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro2);
2143 if (ispro_kpi2==0){
2144 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
2145 ispro_kpi2=1;
2146 }
2147 Com_Multi(pro2,pro_kpi2,pro2b);
2148
2149
2150 B[0] = barrier(1,sKmPip2,sKm,sPip2,rRes2);
2151 B[1] = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
2152 B[2] = barrier(1,sD,sKsKmPip2,sPip1,rD2);
2153
2154 for(int a=0; a<4; a++)
2155 {
2156 for(int j=0; j<4; j++)
2157 {
2158 for(int k=0; k<4; k++)
2159 {
2160 for(int l=0; l<4; l++)
2161 {
2162 temp_PDF += E[a][j][k][l]*pKsKmPip2[a]*(pKmPip2[j]-Km[j])*Pip[k]*(Kp[l]-Pi0[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2163 }
2164 }
2165 }
2166 }
2167
2168 tmp2 = B[0]*B[1]*B[2]*temp_PDF*(-1);
2169 }
2170
2171 //---------------------- for K1270*0 pi+ ,(K*+ -> KS pi+) ---------------------------
2172 if(modetype[i]==16){
2173 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPipPip,sKsPip1,sPip2,rRes2,g2[i],pro1);
2174 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kpi1);
2175 Com_Multi(pro1,pro_kpi1,pro1b);
2176
2177 B[0]= barrier(1,sKsPip1,sKs,sPip1,rRes2);
2178 B[1]= barrier(1,sD,sKsPipPip,sKm,rD2);
2179 B[2]= barrier(2,sKsPipPip,sKsPip1,sPip2,rRes2);
2180 if(g2[i]==0)
2181 {
2182 for(int a=0; a<4; a++)
2183 {
2184 for(int j=0; j<4; j++)
2185 {
2186 temp_PDF += t1D[a]*(pKsPipPip[a]*pKsPipPip[j]/sKsPipPip-G[a][j])*t1KsPip1[j]*G[a][a]*G[j][j];
2187 }
2188 }
2189 tmp1 = B[0]*B[1]*temp_PDF;
2190 }
2191 if(g2[i]==2)
2192 {
2193 for(int a=0; a<4; a++)
2194 {
2195 for(int j=0; j<4; j++)
2196 {
2197 temp_PDF += t1D[a]*t2KsPipPip_kspip[a][j]*t1KsPip1[j]*G[a][a]*G[j][j];
2198 }
2199 }
2200 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2201 }
2202
2203
2204 temp_PDF=0;
2205 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPipPip,sKsPip2,sPip1,rRes2,g2[i],pro1);
2206 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kpi2);
2207 Com_Multi(pro1,pro_kpi2,pro2b);
2208
2209 B[0]= barrier(1,sKsPip2,sKs,sPip2,rRes2);
2210 B[1]= barrier(1,sD,sKsPipPip,sKm,rD2);
2211 B[2]= barrier(2,sKsPipPip,sKsPip2,sPip1,rRes2);
2212 if(g2[i]==0)
2213 {
2214 for(int a=0; a<4; a++)
2215 {
2216 for(int j=0; j<4; j++)
2217 {
2218 temp_PDF += t1D[a]*(pKsPipPip[a]*pKsPipPip[j]/sKsPipPip-G[a][j])*t1KsPip2[j]*G[a][a]*G[j][j];
2219 }
2220 }
2221 tmp2 = B[0]*B[1]*temp_PDF;
2222 }
2223 if(g2[i]==2)
2224 {
2225 for(int a=0; a<4; a++)
2226 {
2227 for(int j=0; j<4; j++)
2228 {
2229 temp_PDF += t1D[a]*t2KsPipPip_kspip2[a][j]*t1KsPip2[j]*G[a][a]*G[j][j];
2230 }
2231 }
2232 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
2233 }
2234 }
2235//K*1270tokpiswave:36 is kmpiswave
2236 if(modetype[i]==36){
2237 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPipPip,sKmPip1,sPip2,rRes2,1,pro1);
2238 if(ispro_kspi1sw==0){
2239 KPiSLASS(sKmPip1,sKm,sPip1,pro_kspi1sw);
2240 ispro_kspi1sw=1;
2241 }
2242 Com_Multi(pro1,pro_kspi1sw,pro1b);
2243 calt1(pKmPipPip,Km,t1D);
2244
2245 if(barr_kskpi1_kspi<0.0) barr_kskpi1_kspi = barrier(1,sKmPipPip,sKmPip1,sPip2,rRes2);
2246 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKmPipPip,sKs,rD2);
2247 for(int a=0; a<4; a++)
2248 {
2249 temp_PDF += G[a][a]*t1D[a]*t1KmPipPip_kmpip[a];
2250// temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kspip[a];
2251 }
2252 tmp1 = barr_kskpi1_kspi*barr_ds_kskpi1*temp_PDF;
2253
2254 temp_PDF=0;
2255 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPipPip,sKmPip2,sPip1,rRes2,1,pro1);
2256 if(ispro_kspi2sw==0){
2257 KPiSLASS(sKmPip2,sKm,sPip2,pro_kspi2sw);
2258 ispro_kspi2sw=1;
2259 }
2260 Com_Multi(pro1,pro_kspi2sw,pro2b);
2261 calt1(pKmPipPip,Km,t1D);
2262
2263 if(barr_kskpi2_kspi<0.0) barr_kskpi2_kspi = barrier(1,sKmPipPip,sKmPip2,sPip1,rRes2);
2264 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKmPipPip,sKs,rD2);
2265 for(int a=0; a<4; a++)
2266 {
2267 temp_PDF += G[a][a]*t1D[a]*t1KmPipPip_kmpip2[a];
2268// temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kspip[a];
2269 }
2270 tmp2 = barr_kskpi2_kspi*barr_ds_kskpi2*temp_PDF;
2271 }
2272
2273//K*1270tokpiswave:36 is kmpiswave
2274 if(modetype[i]==46){
2275 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPipPip,sKsPip1,sPip2,rRes2,1,pro1);
2276 if(ispro_kspi1sw==0){
2277 KPiSLASS(sKsPip1,sKs,sPip1,pro_kspi1sw);
2278 ispro_kspi1sw=1;
2279 }
2280 Com_Multi(pro1,pro_kspi1sw,pro1b);
2281 calt1(pKsPipPip,Kp,t1D);
2282
2283 if(barr_kskpi1_kspi<0.0) barr_kskpi1_kspi = barrier(1,sKsPipPip,sKsPip1,sPip2,rRes2);
2284 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsPipPip,sKm,rD2);
2285 for(int a=0; a<4; a++)
2286 {
2287 temp_PDF += G[a][a]*t1D[a]*t1KsPipPip_kspip[a];
2288// temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kspip[a];
2289 }
2290 tmp1 = barr_kskpi1_kspi*barr_ds_kskpi1*temp_PDF;
2291 temp_PDF=0;
2292 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPipPip,sKsPip2,sPip1,rRes2,1,pro1);
2293 if(ispro_kspi2sw==0){
2294 KPiSLASS(sKsPip2,sKs,sPip2,pro_kspi2sw);
2295 ispro_kspi2sw=1;
2296 }
2297 Com_Multi(pro1,pro_kspi2sw,pro2b);
2298 calt1(pKsPipPip,Kp,t1D);
2299
2300 if(barr_kskpi2_kspi<0.0) barr_kskpi2_kspi = barrier(1,sKsPipPip,sKsPip2,sPip1,rRes2);
2301 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsPipPip,sKm,rD2);
2302 for(int a=0; a<4; a++)
2303 {
2304 temp_PDF += G[a][a]*t1D[a]*t1KsPipPip_kspip2[a];
2305// temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kspip[a];
2306 }
2307 tmp2 = barr_kskpi2_kspi*barr_ds_kskpi2*temp_PDF;
2308 }
2309
2310
2311 //text
2312 /// f1420->kpi_sw
2313 if(modetype[i]==46){
2314 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,1,pro1);
2315 if(ispro_kspi1sw==0){
2316 KPiSLASS(sKsPip1,sKm,sPip1,pro_kpi1sw);
2317 ispro_kpi1sw=1;
2318 }
2319 Com_Multi(pro1,pro_kpi1sw,pro1b);
2320 calt1(pKsKmPip1,Pi0,t1D);
2321
2322 barr_kskpi1_kpi = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
2323 if(barr_ds_kskpi1<0.0) barr_ds_kskpi1 = barrier(1,sD,sKsKmPip1,sPip2,rD2);
2324 for(int a=0; a<4; a++)
2325 {
2326 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip1_kmpip[a];
2327 }
2328 tmp1 = barr_kskpi1_kpi*barr_ds_kskpi1*temp_PDF;
2329
2330 temp_PDF=0;
2331 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro1);
2332 if(ispro_kpi2sw==0){
2333 KPiSLASS(sKsPip2,sKs,sPip2,pro_kpi2sw);
2334 ispro_kpi2sw=1;
2335 }
2336 Com_Multi(pro1,pro_kpi2sw,pro2b);
2337 calt1(pKsKmPip2,Pip,t1D);
2338
2339 barr_kskpi2_kpi = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
2340 if(barr_ds_kskpi2<0.0) barr_ds_kskpi2 = barrier(1,sD,sKsKmPip2,sPip1,rD2);
2341 for(int a=0; a<4; a++)
2342 {
2343 temp_PDF += G[a][a]*t1D[a]*t1KsKmPip2_kmpip[a];
2344 }
2345 tmp2 = barr_kskpi2_kpi*barr_ds_kskpi2*temp_PDF;
2346 }
2347
2348 //1+ Ds -> K1270 Ks , K1270 -> pi+ K*0 ,K*0 ->K+ pi-
2349 if(modetype[i]==6){
2350 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPipPip,sKmPip1,sPip2,rRes2,g2[i],pro1);
2351 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
2352 Com_Multi(pro1,pro_kpi1,pro1b);
2353
2354 B[0]= barrier(1,sKmPip1,sKm,sPip1,rRes2);
2355 B[1]= barrier(1,sD,sKmPipPip,sKs,rD2);
2356 B[2]= barrier(2,sKmPipPip,sKmPip1,sPip2,rRes2);
2357 if(g2[i]==0)
2358 {
2359 for(int a=0; a<4; a++)
2360 {
2361 for(int j=0; j<4; j++)
2362 {
2363 temp_PDF += t1D[a]*(pKmPipPip[a]*pKmPipPip[j]/sKmPipPip-G[a][j])*t1KmPip1[j]*G[a][a]*G[j][j];
2364 }
2365 }
2366 tmp1 = B[0]*B[1]*temp_PDF;
2367 }
2368 if(g2[i]==2)
2369 {
2370 for(int a=0; a<4; a++)
2371 {
2372 for(int j=0; j<4; j++)
2373 {
2374 temp_PDF += t1D[a]*t2KmPipPip_kmpip[a][j]*t1KmPip1[j]*G[a][a]*G[j][j];
2375 }
2376 }
2377 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2378 }
2379
2380
2381 temp_PDF=0;
2382 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPipPip,sKmPip2,sPip1,rRes2,g2[i],pro1);
2383 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
2384 Com_Multi(pro1,pro_kpi2,pro2b);
2385
2386 B[0]= barrier(1,sKmPip2,sKm,sPip2,rRes2);
2387 B[1]= barrier(1,sD,sKmPipPip,sKs,rD2);
2388 B[2]= barrier(2,sKmPipPip,sKmPip2,sPip1,rRes2);
2389 if(g2[i]==0)
2390 {
2391 for(int a=0; a<4; a++)
2392 {
2393 for(int j=0; j<4; j++)
2394 {
2395 temp_PDF += t1D[a]*(pKmPipPip[a]*pKmPipPip[j]/sKmPipPip-G[a][j])*t1KmPip2[j]*G[a][a]*G[j][j];
2396 }
2397 }
2398 tmp2 = B[0]*B[1]*temp_PDF;
2399 }
2400 if(g2[i]==2)
2401 {
2402 for(int a=0; a<4; a++)
2403 {
2404 for(int j=0; j<4; j++)
2405 {
2406 temp_PDF += t1D[a]*t2KmPipPip_kmpip2[a][j]*t1KmPip2[j]*G[a][a]*G[j][j];
2407 }
2408 }
2409 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
2410 }
2411 }
2412
2413 ///Kstp Kmpi 1430
2414 if(modetype[i]==30)
2415 {
2416 if(ispro_kspi1==0){
2417 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
2418 ispro_kspi1=1;
2419 }
2420 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,0,pro_kpi2_0);
2421
2422 Com_Multi(pro_kspi1,pro_kpi2_0,pro1b);
2423 calt1(pKsPip1,pKmPip2,t1D);
2424
2425 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
2426 if(barr_kspi1_kpi2<0.0) barr_kspi1_kpi2 = barrier(1,sD,sKsPip1,sKmPip2,rD2);
2427 for(int a=0; a<4; a++)
2428 {
2429 temp_PDF += t1D[a]*t1KsPip1[a]*G[a][a];
2430 }
2431 tmp1 = barr_kspi1*barr_kspi1_kpi2*temp_PDF;
2432
2433 temp_PDF=0;
2434 if(ispro_kspi2==0) {
2435 propagatorRBW(mass1sq,mass1[i],width1[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
2436 ispro_kspi2=1;
2437 }
2438
2439 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,0,pro_kpi1_0);
2440 Com_Multi(pro_kspi2,pro_kpi1_0,pro2b);
2441
2442 calt1(pKsPip2,pKmPip1,t1D);
2443 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
2444 if(barr_kspi2_kpi1<0.0) barr_kspi2_kpi1 = barrier(1,sD,sKsPip2,sKmPip1,rD2);
2445 for(int a=0; a<4; a++)
2446 {
2447 temp_PDF += t1D[a]*t1KsPip2[a]*G[a][a];
2448 }
2449 tmp2 = barr_kspi2*barr_kspi2_kpi1*temp_PDF;
2450 }
2451
2452 ///Kst0 Kspi 1430
2453 if(modetype[i]==31)
2454 {
2455 if(ispro_kpi2==0){
2456 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
2457 ispro_kpi2=1;
2458 }
2459 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,0,pro_kspi1_0);
2460 Com_Multi(pro_kpi2,pro_kspi1_0,pro1b);
2461
2462 calt1(pKsPip1,pKmPip2,t1D);
2463 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
2464 if(barr_kspi1_kpi2<0.0) barr_kspi1_kpi2 = barrier(1,sD,sKsPip1,sKmPip2,rD2);
2465 for(int a=0; a<4; a++)
2466 {
2467 temp_PDF += t1D[a]*t1KmPip2[a]*G[a][a];
2468 }
2469 tmp1 = barr_kpi2*barr_kspi1_kpi2*temp_PDF;
2470
2471 temp_PDF=0;
2472 if(ispro_kpi1==0){
2473 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
2474 ispro_kpi1=1;
2475 }
2476 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,0,pro_kspi2_0);
2477 Com_Multi(pro_kpi1,pro_kspi2_0,pro2b);
2478
2479 calt1(pKsPip2,pKmPip1,t1D);
2480 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
2481 if(barr_kspi2_kpi1<0.0) barr_kspi2_kpi1 = barrier(1,sD,sKsPip2,sKmPip1,rD2);
2482 for(int a=0; a<4; a++)
2483 {
2484 temp_PDF += t1D[a]*t1KmPip1[a]*G[a][a];
2485 }
2486 tmp2 = barr_kspi2_kpi1*barr_kpi1*temp_PDF;
2487 }
2488
2489
2490 //Kspi 1430 kpi 1430
2491 if(modetype[i]==32)
2492 {
2493 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPip1,sKm,sPip1,rRes2,0,pro_kpi1_0);
2494 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,0,pro_kspi2_0);
2495 Com_Multi(pro_kpi1_0,pro_kpi2_0,pro1b);
2496 tmp1 = 1;
2497 temp_PDF=0;
2498 propagatorRBW(mass1sq,mass1[i],width1[i],sKmPip2,sKm,sPip2,rRes2,0,pro_kpi2_0);
2499 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,0,pro_kspi1_0);
2500 Com_Multi(pro_kspi2_0,pro_kpi1_0,pro2b);
2501 tmp2 = 1;
2502 }
2503
2504 if(modetype[i]==1000)
2505 {
2506 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,1,pro1);
2507 pro1[0]=1;
2508 pro1[1]=0;
2509 if (ispro_kspi1==0){
2510 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
2511 ispro_kspi1=1;
2512 }
2513 Com_Multi(pro1,pro_kspi1,pro1b);
2514 //printf(" %.15f %.15f %.15f\n",B[0],B[1],B[2]);
2515
2516 B[0] = barrier(1,sKsPip1,sKs,sPip1,rRes2);
2517 B[1] = barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
2518 B[2] = barrier(1,sD,sKsKmPip1,sPip2,rD2);
2519 for(int a=0; a<4; a++)
2520 {
2521 for(int j=0; j<4; j++)
2522 {
2523 for(int k=0; k<4; k++)
2524 {
2525 for(int l=0; l<4; l++)
2526 {
2527 temp_PDF += E[a][j][k][l]*pKsKmPip1[a]*(pKsPip1[j]-Kp[j])*Pi0[k]*(Km[l]-Pip[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2528 }
2529 }
2530 }
2531 }
2532 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2533
2534
2535 temp_PDF=0;
2536 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,1,pro2);
2537
2538 pro2[0]=1;
2539 pro2[1]=0;
2540 if (ispro_kspi2==0){
2541 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
2542 ispro_kspi2=1;
2543 }
2544 Com_Multi(pro2,pro_kspi2,pro2b);
2545
2546
2547 B[0] = barrier(1,sKsPip2,sKs,sPip2,rRes2);
2548 B[1] = barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
2549 B[2] = barrier(1,sD,sKsKmPip2,sPip1,rD2);
2550
2551 for(int a=0; a<4; a++)
2552 {
2553 for(int j=0; j<4; j++)
2554 {
2555 for(int k=0; k<4; k++)
2556 {
2557 for(int l=0; l<4; l++)
2558 {
2559 temp_PDF += E[a][j][k][l]*pKsKmPip2[a]*(pKsPip2[j]-Kp[j])*Pip[k]*(Km[l]-Pi0[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2560 }
2561 }
2562 }
2563 }
2564 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
2565 }
2566
2567 //text
2568 // Ds -> phi(1680)pi+ , phi(1680)->K-K*+ ,K*+ ->Km pi+
2569 if(modetype[i]==2000){
2570 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKm,rRes2,1,pro1);
2571 if (ispro_kpi1==0){
2572
2573 pro1[0]=1;
2574 pro1[1]=0;
2575 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
2576 ispro_kpi1=1;
2577 }
2578 Com_Multi(pro1,pro_kpi1,pro1b);
2579
2580 B[0] = barrier(1,sKmPip1,sKm,sPip1,rRes2);
2581 B[1] = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
2582 B[2] = barrier(1,sD,sKsKmPip1,sPip2,rD2);
2583 for(int a=0; a<4; a++)
2584 {
2585 for(int j=0; j<4; j++)
2586 {
2587 for(int k=0; k<4; k++)
2588 {
2589 for(int l=0; l<4; l++)
2590 {
2591 temp_PDF += E[a][j][k][l]*pKsKmPip1[a]*(pKmPip1[j]-Km[j])*Pi0[k]*(Kp[l]-Pip[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2592 }
2593 }
2594 }
2595 }
2596 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
2597
2598
2599 temp_PDF=0;
2600 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro2);
2601 pro2[0]=1;
2602 pro2[1]=0;
2603 if (ispro_kpi2==0){
2604 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
2605 ispro_kpi2=1;
2606 }
2607 Com_Multi(pro2,pro_kpi2,pro2b);
2608
2609
2610 B[0] = barrier(1,sKmPip2,sKm,sPip2,rRes2);
2611 B[1] = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
2612 B[2] = barrier(1,sD,sKsKmPip2,sPip1,rD2);
2613
2614 for(int a=0; a<4; a++)
2615 {
2616 for(int j=0; j<4; j++)
2617 {
2618 for(int k=0; k<4; k++)
2619 {
2620 for(int l=0; l<4; l++)
2621 {
2622 temp_PDF += E[a][j][k][l]*pKsKmPip2[a]*(pKmPip2[j]-Km[j])*Pip[k]*(Kp[l]-Pi0[l])*G[a][a]*G[j][j]*G[k][k]*G[l][l];
2623 }
2624 }
2625 }
2626 }
2627
2628 tmp2 = B[0]*B[1]*B[2]*temp_PDF*(-1);
2629 }
2630
2631 if(modetype[i]==1001)
2632 {
2633 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKmPip1,sKs,rRes2,1,pro1);
2634 pro1[0]=1;
2635 pro1[1]=0;
2636 if(ispro_kpi1==0){
2637 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip1,sKm,sPip1,rRes2,1,pro_kpi1);
2638 ispro_kpi1=1;
2639 }
2640 Com_Multi(pro1,pro_kpi1, pro1b);
2641
2642 if(barr_kpi1<0.0) barr_kpi1 = barrier(1,sKmPip1,sKm,sPip1,rRes2);
2643 barr_kskpi1_kpi = barrier(1,sKsKmPip1,sKmPip1,sKs,rRes2);
2644 for(int a=0; a<4; a++)
2645 {
2646 temp_PDF += Km[a]*t1KmPip1[a]*G[a][a];
2647 }
2648 tmp1 = barr_kpi1*barr_kskpi1_kpi*temp_PDF;
2649
2650 temp_PDF=0;
2651 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKmPip2,sKs,rRes2,1,pro1);
2652 pro1[0]=1;
2653 pro1[1]=0;
2654 if(ispro_kpi2==0){
2655 propagatorRBW(mass2sq,mass2[i],width2[i],sKmPip2,sKm,sPip2,rRes2,1,pro_kpi2);
2656 ispro_kpi2=1;
2657 }
2658 Com_Multi(pro1,pro_kpi2,pro2b);
2659
2660 if(barr_kpi2<0.0) barr_kpi2 = barrier(1,sKmPip2,sKm,sPip2,rRes2);
2661 barr_kskpi2_kpi = barrier(1,sKsKmPip2,sKmPip2,sKs,rRes2);
2662 for(int a=0; a<4; a++)
2663 {
2664 temp_PDF += Km[a]*t1KmPip2[a]*G[a][a];
2665 }
2666 tmp2 = barr_kpi2*barr_kskpi2_kpi*temp_PDF;
2667 }
2668
2669 //Ds -> eta(1475)pi+ eta(1475)->K*Km,K*+>Ks pi+
2670 if(modetype[i]==2001)
2671 {
2672 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip1,sKsPip1,sKm,rRes2,1,pro1);
2673 pro1[0]=1;
2674 pro1[1]=0;
2675
2676 if(ispro_kspi1==0){
2677 propagatorRBW(mass2sq,mass2[i],width2[i],sKsPip1,sKs,sPip1,rRes2,1,pro_kspi1);
2678 ispro_kspi1=1;
2679 }
2680 Com_Multi(pro1,pro_kspi1,pro1b);
2681
2682 if(barr_kspi1<0.0) barr_kspi1 = barrier(1,sKsPip1,sKs,sPip1,rRes2);
2683 if(barr_kskpi1_kspi<0.0)barr_kskpi1_kspi= barrier(1,sKsKmPip1,sKsPip1,sKm,rRes2);
2684 for(int a=0; a<4; a++)
2685 {
2686 temp_PDF += Kp[a]*t1KsPip1[a]*G[a][a];
2687 }
2688 tmp1 = barr_kspi1*barr_kskpi1_kspi*temp_PDF;
2689
2690 temp_PDF=0;
2691 propagatorRBW(mass1sq,mass1[i],width1[i],sKsKmPip2,sKsPip2,sKm,rRes2,1,pro1);
2692 pro1[0]=1;
2693 pro1[1]=0;
2694 if(ispro_kspi2==0){
2695 propagatorRBW(mass1sq,mass2[i],width2[i],sKsPip2,sKs,sPip2,rRes2,1,pro_kspi2);
2696 ispro_kspi2=1;
2697 }
2698 Com_Multi(pro1,pro_kspi2,pro2b);
2699
2700 if(barr_kspi2<0.0) barr_kspi2 = barrier(1,sKsPip2,sKs,sPip2,rRes2);
2701 if(barr_kskpi2_kspi<0.0)barr_kskpi2_kspi= barrier(1,sKsKmPip2,sKsPip2,sKm,rRes2);
2702 for(int a=0; a<4; a++)
2703 {
2704 temp_PDF += Kp[a]*t1KsPip2[a]*G[a][a];
2705 }
2706 tmp2 = barr_kspi2*barr_kskpi2_kspi*temp_PDF;
2707 }
2708 if(modetype[i]==1100||modetype[i]==1200||modetype[i]==1500||modetype[i]==1001||modetype[i]==2001||modetype[i]==1000||modetype[i]==2000||modetype[i]==290||modetype[i]==280||modetype[i]==244||modetype[i]==888||modetype[i]==133||modetype[i]==200||modetype[i]==99||modetype[i]==36||modetype[i]==46||modetype[i]==16||modetype[i]==2||modetype[i]==21||modetype[i]==22||modetype[i]==8||modetype[i]==4||modetype[i]==14||modetype[i]==5||modetype[i]==17||modetype[i]==3||modetype[i]==7||modetype[i]==10||modetype[i]==9||modetype[i]==13||modetype[i]==20||modetype[i]==11||modetype[i]==12||modetype[i]==15||modetype[i]==19||modetype[i]==26||modetype[i]==24||modetype[i]==28||modetype[i]==29||modetype[i]==27||modetype[i]==6||modetype[i]==1||modetype[i]==30||modetype[i]==31||modetype[i]==32){
2709 amp_tmp[0] = tmp1*pro1b[0] + tmp2*pro2b[0];
2710 amp_tmp[1] = tmp1*pro1b[1] + tmp2*pro2b[1];
2711 }
2712 if(modetype[i]==1111){
2713 amp_tmp[0] = 1;
2714 amp_tmp[1] = 0;
2715 }
2716
2717
2718
2719 Com_Multi(amp_tmp,cof,amp_PDF);
2720 PDF[0] += amp_PDF[0];
2721 PDF[1] += amp_PDF[1];
2722 }
2723 double value = PDF[0]*PDF[0] + PDF[1]*PDF[1];
2724 if(value <=0) value = 0.00000000001;
2725 Result = value;
2726return Result;
2727}
2728
2729
2730
2731
2732
sin
double sin(const BesAngle a)
Definition BesAngle.h:210
cos
double cos(const BesAngle a)
Definition BesAngle.h:213
mass
double mass
Definition CosmicGenerator.cxx:138
f
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
g1
TF1 * g1
Definition DataBase/tau_mode.c:66
exp
EvtComplex exp(const EvtComplex &c)
Definition EvtComplex.hh:252
unit
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ !Last cell in buffer $ !No of sampling when dividing cell $ !No of function total $ !Flag for random ceel for $ !Flag for type of for WtMax $ !Flag which decides whether vertices are included in the sampling $ entire domain is hyp !Maximum effective eevents per saves r n generator level $ !Flag for chat level in !Latex Output unit
Definition FoamA.h:90
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
CS
TCrossPart * CS
Definition Mcgpj.cxx:51
t
TTree * t
Definition binning.cxx:23
EvtDToKSKmpippip::decay
void decay(EvtParticle *p)
Definition EvtDToKSKmpippip.cc:135
EvtDToKSKmpippip::clone
EvtDecayBase * clone()
Definition EvtDToKSKmpippip.cc:40
EvtDToKSKmpippip::~EvtDToKSKmpippip
virtual ~EvtDToKSKmpippip()
Definition EvtDToKSKmpippip.cc:34
EvtDToKSKmpippip::getName
void getName(std::string &name)
Definition EvtDToKSKmpippip.cc:36
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:520
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:297
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:504
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:65
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:482
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:34
EvtPDL::getStdHep
static int getStdHep(EvtId id)
Definition EvtPDL.hh:56
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:112
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:120
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:84
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1070
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:179
b
const double b
Definition slope.cxx:9