2 COMPLEX FUNCTION f3pi(IFORM,QQ,SA,SB)
41 parameter( crname =
'F3PI' )
44 REAL MRO,GRO,MRP,GRP,MF2,GF2,MF0,GF0,MSG,GSG
45 REAL M1,M2,M3,M1SQ,M2SQ,M3SQ,MPIZ,MPIC
47 REAL F134,F150,F15A,F15B,F167
48 REAL F34A,F34B,F35,F35A,F35B,F36A,F36B
49 COMPLEX BT1,BT2,BT3,BT4,BT5,BT6,BT7
50 COMPLEX FRO1,FRO2,FRP1,FRP2
51 COMPLEX FF21,FF22,FF23,FSG1,FSG2,FSG3,FF01,FF02,FF03
81 bt2 = cmplx(0.12,0.)*cexp(cmplx(0., 0.99*pi))
82 bt3 = cmplx(0.37,0.)*cexp(cmplx(0.,-0.15*pi))
83 bt4 = cmplx(0.87,0.)*cexp(cmplx(0., 0.53*pi))
84 bt5 = cmplx(0.71,0.)*cexp(cmplx(0., 0.56*pi))
85 bt6 = cmplx(2.10,0.)*cexp(cmplx(0., 0.23*pi))
86 bt7 = cmplx(0.77,0.)*cexp(cmplx(0.,-0.54*pi))
88 print *,
' In F3pi: add (rho-pi S-wave) + (rhop-pi S-wave) +'
89 print *,
' (rho-pi D-wave) + (rhop-pi D-wave) +'
90 print *,
' (f2 pi D-wave) + (sigmapi S-wave) + (f0pi S-wave)'
129 IF (iform.EQ.1.OR.iform.EQ.2)
THEN
132 s3 = qq-sa-sb+m1sq+m2sq+m3sq
133 IF (s3.LE.0..OR.s2.LE.0.)
RETURN
138 f134 = -(1./3.)*((s3-m3sq)-(s1-m1sq))
139 f150 = (1./18.)*(qq-m3sq+s3)*(2.*m1sq+2.*m2sq-s3)/s3
143 fro1 = bwigml(s1,mro,gro,m2,m3,1)
144 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
145 fro2 = bwigml(s2,mro,gro,m3,m1,1)
146 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
147 ff23 = bwigml(s3,mf2,gf2,m1,m2,2)
148 fsg3 = bwigml(s3,msg,gsg,m1,m2,0)
149 ff03 = bwigml(s3,mf0,gf0,m1,m2,0)
151 f3pi = bt1*fro1+bt2*frp1+
152 1 bt3*cmplx(f134,0.)*fro2+bt4*cmplx(f134,0.)*frp2+
153 1 bt5*cmplx(f150,0.)*ff23+
154 1 bt6*cmplx(f167,0.)*fsg3+bt7*cmplx(f167,0.)*ff03
157 ELSEIF (idk.EQ.2)
THEN
160 f134 = -(1./3.)*((s3-m3sq)-(s1-m1sq))
161 f15a = -(1./2.)*((s2-m2sq)-(s3-m3sq))
162 f15b = -(1./18.)*(qq-m2sq+s2)*(2.*m1sq+2.*m3sq-s2)/s2
166 fro1 = bwigml(s1,mro,gro,m2,m3,1)
167 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
168 fro2 = bwigml(s2,mro,gro,m3,m1,1)
169 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
170 ff21 = bwigml(s1,mf2,gf2,m2,m3,2)
171 ff22 = bwigml(s2,mf2,gf2,m3,m1,2)
172 fsg2 = bwigml(s2,msg,gsg,m3,m1,0)
173 ff02 = bwigml(s2,mf0,gf0,m3,m1,0)
175 f3pi = bt1*fro1+bt2*frp1+
176 1 bt3*cmplx(f134,0.)*fro2+bt4*cmplx(f134,0.)*frp2
177 1 -bt5*cmplx(f15a,0.)*ff21-bt5*cmplx(f15b,0.)*ff22
178 1 -bt6*cmplx(f167,0.)*fsg2-bt7*cmplx(f167,0.)*ff02
183 ELSE IF (iform.EQ.3)
THEN
186 s2 = qq-sa-sb+m1sq+m2sq+m3sq
187 IF (s1.LE.0..OR.s2.LE.0.)
RETURN
192 f34a = (1./3.)*((s2-m2sq)-(s3-m3sq))
193 f34b = (1./3.)*((s3-m3sq)-(s1-m1sq))
194 f35 =-(1./2.)*((s1-m1sq)-(s2-m2sq))
197 fro1 = bwigml(s1,mro,gro,m2,m3,1)
198 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
199 fro2 = bwigml(s2,mro,gro,m3,m1,1)
200 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
201 ff23 = bwigml(s3,mf2,gf2,m1,m2,2)
204 1 bt3*(cmplx(f34a,0.)*fro1+cmplx(f34b,0.)*fro2)+
205 1 bt4*(cmplx(f34a,0.)*frp1+cmplx(f34b,0.)*frp2)+
206 1 bt5*cmplx(f35,0.)*ff23
209 ELSEIF (idk.EQ.2)
THEN
212 f34a = (1./3.)*((s2-m2sq)-(s3-m3sq))
213 f34b = (1./3.)*((s3-m3sq)-(s1-m1sq))
214 f35a = -(1./18.)*(qq-m1sq+s1)*(2.*m2sq+2.*m3sq-s1)/s1
215 f35b = (1./18.)*(qq-m2sq+s2)*(2.*m3sq+2.*m1sq-s2)/s2
220 fro1 = bwigml(s1,mro,gro,m2,m3,1)
221 frp1 = bwigml(s1,mrp,grp,m2,m3,1)
222 fro2 = bwigml(s2,mro,gro,m3,m1,1)
223 frp2 = bwigml(s2,mrp,grp,m3,m1,1)
224 ff21 = bwigml(s1,mf2,gf2,m2,m3,2)
225 ff22 = bwigml(s2,mf2,gf2,m3,m1,2)
226 fsg1 = bwigml(s1,msg,gsg,m2,m3,0)
227 fsg2 = bwigml(s2,msg,gsg,m3,m1,0)
228 ff01 = bwigml(s1,mf0,gf0,m2,m3,0)
229 ff02 = bwigml(s2,mf0,gf0,m3,m1,0)
232 1 bt3*(cmplx(f34a,0.)*fro1+cmplx(f34b,0.)*fro2)+
233 1 bt4*(cmplx(f34a,0.)*frp1+cmplx(f34b,0.)*frp2)
234 1 -bt5*(cmplx(f35a,0.)*ff21+cmplx(f35b,0.)*ff22)
235 1 -bt6*(cmplx(f36a,0.)*fsg1+cmplx(f36b,0.)*fsg2)
236 1 -bt7*(cmplx(f36a,0.)*ff01+cmplx(f36b,0.)*ff02)
249 COMPLEX FUNCTION bwigml(S,M,G,M1,M2,L)
255 REAL MSQ,W,WGS,MP,MM,QS,QM
262 IF (w.GT.(m1+m2))
THEN
263 qs=sqrt(abs((s -mp)*(s -mm)))/w
264 qm=sqrt(abs((msq -mp)*(msq -mm)))/m
266 wgs=g*(msq/w)*(qs/qm)**ipow
269 bwigml=cmplx(msq,0.)/cmplx(msq-s,-wgs)
274 COMPLEX FUNCTION fa1a1p(XMSQ)
281 REAL XM1,XG1,XM2,XG2,XM1SQ,XM2SQ,GG1,GG2,GF,FG1,FG2
285 IF (ifirst.EQ.0)
THEN
293 bet = cmplx(pkorb(3,17),0.)
295 gg1 = xm1*xg1/(1.3281*0.806)
296 gg2 = xm2*xg2/(1.3281*0.806)
305 f1 = cmplx(-xm1sq,0.0)/cmplx(xmsq-xm1sq,fg1)
306 f2 = cmplx(-xm2sq,0.0)/cmplx(xmsq-xm2sq,fg2)
320 DOUBLE PRECISION MKST,MK,MK1SQ,MK2SQ,C3PI,CKST
321 DOUBLE PRECISION S,WGA1C,WGA1N,WG3PIC,WG3PIN,GKST
325 IF (ifirst.NE.987)
THEN
335 ckst = 4.7621d0**2*c3pi
347 IF (s.GT.mk1sq) gkst = sqrt((s-mk1sq)*(s-mk2sq))/(2.*s)
349 wga1 = sngl(c3pi*(wg3pic+wg3pin)+ckst*gkst)
354 DOUBLE PRECISION FUNCTION wga1c(S)
358 DOUBLE PRECISION S,STH,Q0,Q1,Q2,P0,P1,P2,P3,P4,G1_IM
360 parameter(q0 = 5.80900d0,q1 = -3.00980d0,q2 = 4.57920d0,
361 1 p0 = -13.91400d0,p1 = 27.67900d0,p2 = -13.39300d0,
362 2 p3 = 3.19240d0,p4 = -0.10487d0)
364 parameter(sth = 0.1753d0)
369 ELSEIF((s.GT.sth).AND.(s.LT.0.823d0))
THEN
370 g1_im = q0*(s-sth)**3*(1. + q1*(s-sth) + q2*(s-sth)**2)
372 g1_im = p0 + p1*s + p2*s**2+ p3*s**3 + p4*s**4
379 DOUBLE PRECISION FUNCTION wga1n(S)
383 DOUBLE PRECISION S,STH,Q0,Q1,Q2,P0,P1,P2,P3,P4,G1_IM
385 parameter(q0 = 6.28450d0,q1 = -2.95950d0,q2 = 4.33550d0,
386 1 p0 = -15.41100d0,p1 = 32.08800d0,p2 = -17.66600d0,
387 2 p3 = 4.93550d0,p4 = -0.37498d0)
389 parameter(sth = 0.1676d0)
394 ELSEIF((s.GT.sth).AND.(s.LT.0.823d0))
THEN
395 g1_im = q0*(s-sth)**3*(1. + q1*(s-sth) + q2*(s-sth)**2)
397 g1_im = p0 + p1*s + p2*s**2+ p3*s**3 + p4*s**4