2 COMPLEX FUNCTION f3pi(IFORM,QQ,SA,SB)
27 #if defined(CLEO_TYPECHEK) 43 parameter( crname =
'F3PI' )
46 REAL MRO,GRO,MRP,GRP,MF2,GF2,MF0,GF0,MSG,GSG
47 REAL M1,M2,M3,M1SQ,M2SQ,M3SQ,MPIZ,MPIC
49 REAL F134,F150,F15A,F15B,F167
50 REAL F34A,F34B,F35,F35A,F35B,F36A,F36B
51 COMPLEX BT1,BT2,BT3,BT4,BT5,BT6,BT7
52 COMPLEX FRO1,FRO2,FRP1,FRP2
53 COMPLEX FF21,FF22,FF23,FSG1,FSG2,FSG3,FF01,FF02,FF03
83 bt2 = cmplx(0.12,0.)*cexp(cmplx(0., 0.99*pi))
84 bt3 = cmplx(0.37,0.)*cexp(cmplx(0.,-0.15*pi))
85 bt4 = cmplx(0.87,0.)*cexp(cmplx(0., 0.53*pi))
86 bt5 = cmplx(0.71,0.)*cexp(cmplx(0., 0.56*pi))
87 bt6 = cmplx(2.10,0.)*cexp(cmplx(0., 0.23*pi))
88 bt7 = cmplx(0.77,0.)*cexp(cmplx(0.,-0.54*pi))
90 print *,
' In F3pi: add (rho-pi S-wave) + (rhop-pi S-wave) +' 91 print *,
' (rho-pi D-wave) + (rhop-pi D-wave) +' 92 print *,
' (f2 pi D-wave) + (sigmapi S-wave) + (f0pi S-wave)' 131 IF (iform.EQ.1.OR.iform.EQ.2)
THEN 134 s3 = qq-sa-sb+m1sq+m2sq+m3sq
135 IF (s3.LE.0..OR.s2.LE.0.)
RETURN 140 f134 = -(1./3.)*((s3-m3sq)-(s1-m1sq))
141 f150 = (1./18.)*(qq-m3sq+s3)*(2.*m1sq+2.*m2sq-s3)/s3
145 fro1 = bwigml(s1,mro,gro,m2,m3,1)
146 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
147 fro2 = bwigml(s2,mro,gro,m3,m1,1)
148 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
149 ff23 = bwigml(s3,mf2,gf2,m1,m2,2)
150 fsg3 = bwigml(s3,msg,gsg,m1,m2,0)
151 ff03 = bwigml(s3,mf0,gf0,m1,m2,0)
153 f3pi = bt1*fro1+bt2*frp1+
154 1 bt3*cmplx(f134,0.)*fro2+bt4*cmplx(f134,0.)*frp2+
155 1 bt5*cmplx(f150,0.)*ff23+
156 1 bt6*cmplx(f167,0.)*fsg3+bt7*cmplx(f167,0.)*ff03
159 ELSEIF (idk.EQ.2)
THEN 162 f134 = -(1./3.)*((s3-m3sq)-(s1-m1sq))
163 f15a = -(1./2.)*((s2-m2sq)-(s3-m3sq))
164 f15b = -(1./18.)*(qq-m2sq+s2)*(2.*m1sq+2.*m3sq-s2)/s2
168 fro1 = bwigml(s1,mro,gro,m2,m3,1)
169 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
170 fro2 = bwigml(s2,mro,gro,m3,m1,1)
171 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
172 ff21 = bwigml(s1,mf2,gf2,m2,m3,2)
173 ff22 = bwigml(s2,mf2,gf2,m3,m1,2)
174 fsg2 = bwigml(s2,msg,gsg,m3,m1,0)
175 ff02 = bwigml(s2,mf0,gf0,m3,m1,0)
177 f3pi = bt1*fro1+bt2*frp1+
178 1 bt3*cmplx(f134,0.)*fro2+bt4*cmplx(f134,0.)*frp2
179 1 -bt5*cmplx(f15a,0.)*ff21-bt5*cmplx(f15b,0.)*ff22
180 1 -bt6*cmplx(f167,0.)*fsg2-bt7*cmplx(f167,0.)*ff02
185 ELSE IF (iform.EQ.3)
THEN 188 s2 = qq-sa-sb+m1sq+m2sq+m3sq
189 IF (s1.LE.0..OR.s2.LE.0.)
RETURN 194 f34a = (1./3.)*((s2-m2sq)-(s3-m3sq))
195 f34b = (1./3.)*((s3-m3sq)-(s1-m1sq))
196 f35 =-(1./2.)*((s1-m1sq)-(s2-m2sq))
199 fro1 = bwigml(s1,mro,gro,m2,m3,1)
200 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
201 fro2 = bwigml(s2,mro,gro,m3,m1,1)
202 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
203 ff23 = bwigml(s3,mf2,gf2,m1,m2,2)
206 1 bt3*(cmplx(f34a,0.)*fro1+cmplx(f34b,0.)*fro2)+
207 1 bt4*(cmplx(f34a,0.)*frp1+cmplx(f34b,0.)*frp2)+
208 1 bt5*cmplx(f35,0.)*ff23
211 ELSEIF (idk.EQ.2)
THEN 214 f34a = (1./3.)*((s2-m2sq)-(s3-m3sq))
215 f34b = (1./3.)*((s3-m3sq)-(s1-m1sq))
216 f35a = -(1./18.)*(qq-m1sq+s1)*(2.*m2sq+2.*m3sq-s1)/s1
217 f35b = (1./18.)*(qq-m2sq+s2)*(2.*m3sq+2.*m1sq-s2)/s2
222 fro1 = bwigml(s1,mro,gro,m2,m3,1)
223 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
224 fro2 = bwigml(s2,mro,gro,m3,m1,1)
225 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
226 ff21 = bwigml(s1,mf2,gf2,m2,m3,2)
227 ff22 = bwigml(s2,mf2,gf2,m3,m1,2)
228 fsg1 = bwigml(s1,msg,gsg,m2,m3,0)
229 fsg2 = bwigml(s2,msg,gsg,m3,m1,0)
230 ff01 = bwigml(s1,mf0,gf0,m2,m3,0)
231 ff02 = bwigml(s2,mf0,gf0,m3,m1,0)
234 1 bt3*(cmplx(f34a,0.)*fro1+cmplx(f34b,0.)*fro2)+
235 1 bt4*(cmplx(f34a,0.)*frp1+cmplx(f34b,0.)*frp2)
236 1 -bt5*(cmplx(f35a,0.)*ff21+cmplx(f35b,0.)*ff22)
237 1 -bt6*(cmplx(f36a,0.)*fsg1+cmplx(f36b,0.)*fsg2)
238 1 -bt7*(cmplx(f36a,0.)*ff01+cmplx(f36b,0.)*ff02)
251 COMPLEX FUNCTION bwigml(S,M,G,M1,M2,L)
257 REAL MSQ,W,WGS,MP,MM,QS,QM
264 IF (w.GT.(m1+m2))
THEN 265 qs=sqrt(abs((s -mp)*(s -mm)))/w
266 qm=sqrt(abs((msq -mp)*(msq -mm)))/m
268 wgs=g*(msq/w)*(qs/qm)**ipow
271 bwigml=cmplx(msq,0.)/cmplx(msq-s,-wgs)
276 COMPLEX FUNCTION fa1a1p(XMSQ)
283 REAL XM1,XG1,XM2,XG2,XM1SQ,XM2SQ,GG1,GG2,GF,FG1,FG2
287 IF (ifirst.EQ.0)
THEN 295 bet = cmplx(pkorb(3,17),0.)
297 gg1 = xm1*xg1/(1.3281*0.806)
298 gg2 = xm2*xg2/(1.3281*0.806)
307 f1 = cmplx(-xm1sq,0.0)/cmplx(xmsq-xm1sq,fg1)
308 f2 = cmplx(-xm2sq,0.0)/cmplx(xmsq-xm2sq,fg2)
322 DOUBLE PRECISION MKST,MK,MK1SQ,MK2SQ,C3PI,CKST
323 DOUBLE PRECISION S,WGA1C,WGA1N,WG3PIC,WG3PIN,GKST
327 IF (ifirst.NE.987)
THEN 337 ckst = 4.7621d0**2*c3pi
349 IF (s.GT.mk1sq) gkst = sqrt((s-mk1sq)*(s-mk2sq))/(2.*s)
351 wga1 = sngl(c3pi*(wg3pic+wg3pin)+ckst*gkst)
356 DOUBLE PRECISION FUNCTION wga1c(S)
360 DOUBLE PRECISION S,STH,Q0,Q1,Q2,P0,P1,P2,P3,P4,G1_IM
362 parameter(q0 = 5.80900d0,q1 = -3.00980d0,q2 = 4.57920d0,
363 1 p0 = -13.91400d0,p1 = 27.67900d0,p2 = -13.39300d0,
364 2 p3 = 3.19240d0,p4 = -0.10487d0)
366 parameter(sth = 0.1753d0)
371 ELSEIF((s.GT.sth).AND.(s.LT.0.823d0))
THEN 372 g1_im = q0*(s-sth)**3*(1. + q1*(s-sth) + q2*(s-sth)**2)
374 g1_im = p0 + p1*s + p2*s**2+ p3*s**3 + p4*s**4
381 DOUBLE PRECISION FUNCTION wga1n(S)
385 DOUBLE PRECISION S,STH,Q0,Q1,Q2,P0,P1,P2,P3,P4,G1_IM
387 parameter(q0 = 6.28450d0,q1 = -2.95950d0,q2 = 4.33550d0,
388 1 p0 = -15.41100d0,p1 = 32.08800d0,p2 = -17.66600d0,
389 2 p3 = 4.93550d0,p4 = -0.37498d0)
391 parameter(sth = 0.1676d0)
396 ELSEIF((s.GT.sth).AND.(s.LT.0.823d0))
THEN 397 g1_im = q0*(s-sth)**3*(1. + q1*(s-sth) + q2*(s-sth)**2)
399 g1_im = p0 + p1*s + p2*s**2+ p3*s**3 + p4*s**4