1 DOUBLE PRECISION FUNCTION gfact(QQ)
5 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
6 * ,ampiz,ampi,amro,gamro,ama1,gama1
7 * ,amk,amkz,amkst,gamkst
9 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
10 * ,ampiz,ampi,amro,gamro,ama1,gama1
11 * ,amk,amkz,amkst,gamkst
12 DOUBLE PRECISION dgamqq,qq,thr,thu,th0,lam
14 DOUBLE PRECISION af4,af5,aa1,aa2,aa4
15 DOUBLE PRECISION f1,f2,f3,f4,f5
16 DOUBLE PRECISION aa,bb,cc,dd,a,b,c,d,x
17 SAVE a,b,c,d,af4,af5,aa,bb,cc,dd,thr,thu,th0
19 CALL ifgfact(1,ii,init)
24 CALL getff3piscal(inum)
38 CALL setff3piscal(inum)
47 bb=1./8.*(10.*aa2-aa4-16.*aa1)
48 aa=(8.*aa1-aa2-4.*bb)/6.
54 d=-9*(-4*f3-f1+4*f2+f4)
55 c=3*(f4+f1-2*f3-11./18.*d)
68 gfact=af4+(af5-af4)*(qq-3*thu)*2/thu
69 ELSEIF (qq.GT.thr)
THEN
70 gfact=a*qq+b+c/qq+d/qq**2
71 ELSEIF(qq.LE.th0)
THEN
75 gfact=aa*x+bb*x**2+cc*x**3
81 DOUBLE PRECISION FUNCTION dgamqq(XQQB)
89 common/precint/ epssq,abs1
90 DOUBLE PRECISION epssq,abs1
91 COMMON /
EXTERNAL/ xqqa
93 EXTERNAL gaus,dgamqqs1
94 DOUBLE PRECISION gaus,dgamqqs1
95 DOUBLE PRECISION xqqb,eps,ups1,downs1
96 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
97 * ,ampiz,ampi,amro,gamro,ama1,gama1
98 * ,amk,amkz,amkst,gamkst
100 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
101 * ,ampiz,ampi,amro,gamro,ama1,gama1
102 * ,amk,amkz,amkst,gamkst
106 ups1=(dsqrt(xqqb) - ampi)**2-abs1
107 downs1=4.d0*ampi**2+abs1
109 dgamqq = gaus(dgamqqs1,downs1,ups1,eps)
116 DOUBLE PRECISION FUNCTION dgamqqs1(S1)
127 EXTERNAL gaus2,ffwid3pi,dgamqqs1s3
128 DOUBLE PRECISION gaus2,ffwid3pi,dgamqqs1s3
129 COMMON /internal/ s1a
131 COMMON /
EXTERNAL/ xqqa
132 DOUBLE PRECISION xqqa
133 common/precint/ epssq,abs1
134 DOUBLE PRECISION epssq,abs1
135 DOUBLE PRECISION eps,ups3,downs3
136 DOUBLE PRECISION xlam,x,y,z
137 DOUBLE PRECISION xampi2
138 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
139 * ,ampiz,ampi,amro,gamro,ama1,gama1
140 * ,amk,amkz,amkst,gamkst
142 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
143 * ,ampiz,ampi,amro,gamro,ama1,gama1
144 * ,amk,amkz,amkst,gamkst
146 xlam(x,y,z) = sqrt(abs((x-y-z)**2 - 4.*y*z))
151 ups3 = (xqqa - ampi**2)**2 -
152 & ( xlam(xqqa,s1,xampi2)
153 & - xlam(s1,xampi2,xampi2) )**2
154 downs3 = (xqqa - ampi**2)**2 -
155 & (xlam(xqqa,s1,xampi2)
156 & + xlam(s1,xampi2,xampi2) )**2
159 downs3 = downs3/4./s1
161 dgamqqs1 = gaus2(dgamqqs1s3,downs3,ups3,eps)
167 DOUBLE PRECISION FUNCTION dgamqqs1s3(XS3)
174 DOUBLE PRECISION ffwid3pi,xs3
175 COMMON /internal/ xs1a
176 DOUBLE PRECISION xs1a
177 COMMON /
EXTERNAL/ xqqa
178 DOUBLE PRECISION xqqa
179 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
180 * ,ampiz,ampi,amro,gamro,ama1,gama1
181 * ,amk,amkz,amkst,gamkst
183 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
184 * ,ampiz,ampi,amro,gamro,ama1,gama1
185 * ,amk,amkz,amkst,gamkst
188 dgamqqs1s3 = ffwid3pi(xqqa,xs1a,xs3)