4 COMPLEX FUNCTION f3pi(IFORM,QQ,SA,SB)
45 parameter( crname =
'F3PI' )
48 REAL MRO,GRO,MRP,GRP,MF2,GF2,MF0,GF0,MSG,GSG
49 REAL M1,M2,M3,M1SQ,M2SQ,M3SQ,MPIZ,MPIC
51 REAL F134,F150,F15A,F15B,F167
52 REAL F34A,F34B,F35,F35A,F35B,F36A,F36B
53 COMPLEX BT1,BT2,BT3,BT4,BT5,BT6,BT7
54 COMPLEX FRO1,FRO2,FRP1,FRP2
55 COMPLEX FF21,FF22,FF23,FSG1,FSG2,FSG3,FF01,FF02,FF03
85 bt2 = cmplx(0.12,0.)*cexp(cmplx(0., 0.99*pi))
86 bt3 = cmplx(0.37,0.)*cexp(cmplx(0.,-0.15*pi))
87 bt4 = cmplx(0.87,0.)*cexp(cmplx(0., 0.53*pi))
88 bt5 = cmplx(0.71,0.)*cexp(cmplx(0., 0.56*pi))
89 bt6 = cmplx(2.10,0.)*cexp(cmplx(0., 0.23*pi))
90 bt7 = cmplx(0.77,0.)*cexp(cmplx(0.,-0.54*pi))
92 print *,
' In F3pi: add (rho-pi S-wave) + (rhop-pi S-wave) +'
93 print *,
' (rho-pi D-wave) + (rhop-pi D-wave) +'
94 print *,
' (f2 pi D-wave) + (sigmapi S-wave) + (f0pi S-wave)'
133 IF (iform.EQ.1.OR.iform.EQ.2)
THEN
136 s3 = qq-sa-sb+m1sq+m2sq+m3sq
137 IF (s3.LE.0..OR.s2.LE.0.)
RETURN
142 f134 = -(1./3.)*((s3-m3sq)-(s1-m1sq))
143 f150 = (1./18.)*(qq-m3sq+s3)*(2.*m1sq+2.*m2sq-s3)/s3
147 fro1 = bwigml(s1,mro,gro,m2,m3,1)
148 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
149 fro2 = bwigml(s2,mro,gro,m3,m1,1)
150 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
151 ff23 = bwigml(s3,mf2,gf2,m1,m2,2)
152 fsg3 = bwigml(s3,msg,gsg,m1,m2,0)
153 ff03 = bwigml(s3,mf0,gf0,m1,m2,0)
155 f3pi = bt1*fro1+bt2*frp1+
156 1 bt3*cmplx(f134,0.)*fro2+bt4*cmplx(f134,0.)*frp2+
157 1 bt5*cmplx(f150,0.)*ff23+
158 1 bt6*cmplx(f167,0.)*fsg3+bt7*cmplx(f167,0.)*ff03
161 ELSEIF (idk.EQ.2)
THEN
164 f134 = -(1./3.)*((s3-m3sq)-(s1-m1sq))
165 f15a = -(1./2.)*((s2-m2sq)-(s3-m3sq))
166 f15b = -(1./18.)*(qq-m2sq+s2)*(2.*m1sq+2.*m3sq-s2)/s2
170 fro1 = bwigml(s1,mro,gro,m2,m3,1)
171 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
172 fro2 = bwigml(s2,mro,gro,m3,m1,1)
173 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
174 ff21 = bwigml(s1,mf2,gf2,m2,m3,2)
175 ff22 = bwigml(s2,mf2,gf2,m3,m1,2)
176 fsg2 = bwigml(s2,msg,gsg,m3,m1,0)
177 ff02 = bwigml(s2,mf0,gf0,m3,m1,0)
179 f3pi = bt1*fro1+bt2*frp1+
180 1 bt3*cmplx(f134,0.)*fro2+bt4*cmplx(f134,0.)*frp2
181 1 -bt5*cmplx(f15a,0.)*ff21-bt5*cmplx(f15b,0.)*ff22
182 1 -bt6*cmplx(f167,0.)*fsg2-bt7*cmplx(f167,0.)*ff02
187 ELSE IF (iform.EQ.3)
THEN
190 s2 = qq-sa-sb+m1sq+m2sq+m3sq
191 IF (s1.LE.0..OR.s2.LE.0.)
RETURN
196 f34a = (1./3.)*((s2-m2sq)-(s3-m3sq))
197 f34b = (1./3.)*((s3-m3sq)-(s1-m1sq))
198 f35 =-(1./2.)*((s1-m1sq)-(s2-m2sq))
201 fro1 = bwigml(s1,mro,gro,m2,m3,1)
202 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
203 fro2 = bwigml(s2,mro,gro,m3,m1,1)
204 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
205 ff23 = bwigml(s3,mf2,gf2,m1,m2,2)
208 1 bt3*(cmplx(f34a,0.)*fro1+cmplx(f34b,0.)*fro2)+
209 1 bt4*(cmplx(f34a,0.)*frp1+cmplx(f34b,0.)*frp2)+
210 1 bt5*cmplx(f35,0.)*ff23
213 ELSEIF (idk.EQ.2)
THEN
216 f34a = (1./3.)*((s2-m2sq)-(s3-m3sq))
217 f34b = (1./3.)*((s3-m3sq)-(s1-m1sq))
218 f35a = -(1./18.)*(qq-m1sq+s1)*(2.*m2sq+2.*m3sq-s1)/s1
219 f35b = (1./18.)*(qq-m2sq+s2)*(2.*m3sq+2.*m1sq-s2)/s2
224 fro1 = bwigml(s1,mro,gro,m2,m3,1)
225 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
226 fro2 = bwigml(s2,mro,gro,m3,m1,1)
227 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
228 ff21 = bwigml(s1,mf2,gf2,m2,m3,2)
229 ff22 = bwigml(s2,mf2,gf2,m3,m1,2)
230 fsg1 = bwigml(s1,msg,gsg,m2,m3,0)
231 fsg2 = bwigml(s2,msg,gsg,m3,m1,0)
232 ff01 = bwigml(s1,mf0,gf0,m2,m3,0)
233 ff02 = bwigml(s2,mf0,gf0,m3,m1,0)
236 1 bt3*(cmplx(f34a,0.)*fro1+cmplx(f34b,0.)*fro2)+
237 1 bt4*(cmplx(f34a,0.)*frp1+cmplx(f34b,0.)*frp2)
238 1 -bt5*(cmplx(f35a,0.)*ff21+cmplx(f35b,0.)*ff22)
239 1 -bt6*(cmplx(f36a,0.)*fsg1+cmplx(f36b,0.)*fsg2)
240 1 -bt7*(cmplx(f36a,0.)*ff01+cmplx(f36b,0.)*ff02)
253 COMPLEX FUNCTION bwigml(S,M,G,M1,M2,L)
259 REAL MSQ,W,WGS,MP,MM,QS,QM
266 IF (w.GT.(m1+m2))
THEN
267 qs=sqrt(abs((s -mp)*(s -mm)))/w
268 qm=sqrt(abs((msq -mp)*(msq -mm)))/m
270 wgs=g*(msq/w)*(qs/qm)**ipow
273 bwigml=cmplx(msq,0.)/cmplx(msq-s,-wgs)
278 COMPLEX FUNCTION fa1a1p(XMSQ)
285 REAL XM1,XG1,XM2,XG2,XM1SQ,XM2SQ,GG1,GG2,GF,FG1,FG2
289 IF (ifirst.EQ.0)
THEN
297 bet = cmplx(pkorb(3,17),0.)
299 gg1 = xm1*xg1/(1.3281*0.806)
300 gg2 = xm2*xg2/(1.3281*0.806)
309 f1 = cmplx(-xm1sq,0.0)/cmplx(xmsq-xm1sq,fg1)
310 f2 = cmplx(-xm2sq,0.0)/cmplx(xmsq-xm2sq,fg2)
324 DOUBLE PRECISION MKST,MK,MK1SQ,MK2SQ,C3PI,CKST
325 DOUBLE PRECISION S,WGA1C,WGA1N,WG3PIC,WG3PIN,GKST
329 IF (ifirst.NE.987)
THEN
339 ckst = 4.7621d0**2*c3pi
351 IF (s.GT.mk1sq) gkst = sqrt((s-mk1sq)*(s-mk2sq))/(2.*s)
353 wga1 = sngl(c3pi*(wg3pic+wg3pin)+ckst*gkst)
358 DOUBLE PRECISION FUNCTION wga1c(S)
362 DOUBLE PRECISION S,STH,Q0,Q1,Q2,P0,P1,P2,P3,P4,G1_IM
364 parameter(q0 = 5.80900d0,q1 = -3.00980d0,q2 = 4.57920d0,
365 1 p0 = -13.91400d0,p1 = 27.67900d0,p2 = -13.39300d0,
366 2 p3 = 3.19240d0,p4 = -0.10487d0)
368 parameter(sth = 0.1753d0)
373 ELSEIF((s.GT.sth).AND.(s.LT.0.823d0))
THEN
374 g1_im = q0*(s-sth)**3*(1. + q1*(s-sth) + q2*(s-sth)**2)
376 g1_im = p0 + p1*s + p2*s**2+ p3*s**3 + p4*s**4
383 DOUBLE PRECISION FUNCTION wga1n(S)
387 DOUBLE PRECISION S,STH,Q0,Q1,Q2,P0,P1,P2,P3,P4,G1_IM
389 parameter(q0 = 6.28450d0,q1 = -2.95950d0,q2 = 4.33550d0,
390 1 p0 = -15.41100d0,p1 = 32.08800d0,p2 = -17.66600d0,
391 2 p3 = 4.93550d0,p4 = -0.37498d0)
393 parameter(sth = 0.1676d0)
398 ELSEIF((s.GT.sth).AND.(s.LT.0.823d0))
THEN
399 g1_im = q0*(s-sth)**3*(1. + q1*(s-sth) + q2*(s-sth)**2)
401 g1_im = p0 + p1*s + p2*s**2+ p3*s**3 + p4*s**4