2 implicit real*8 (a-h,o-z)
5 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
6 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
7 common /input/ su,su2,qq2,p1,p2,ngen,iseed,mode,iww,nhit
8 common /param/ pi,alpha,f_max
9 common /cbwgrho/ rhom2,rho1m2,rho2m2,omm2,rhomg,rho1mg,rho2mg
11 common /cbwga1/ a1m2,con
12 common /anom/amrop,gamrop,sig,amrop_2,amropg
13 common /cbwgrho_t/ rho1m2_t,rho1mg_t,beta
15 pi = 3.141592653589793238d0
16 alpha = 1.d0/137.0359895d0
22 coupl = sqrt(6.d0)/fpi**2
42 sgo = 1.55d0/sqrt(2.d0)
55 con = a1g*a1m/gfun8(a1m2)
68 rho1mg_t = rho1m_t*rho1g_t
73 complex*16 function anom_bwg(q1_2,q2_2)
74 implicit real*8 (a-h,o-z)
76 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
77 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
78 common /cbwgrho/ rhom2,rho1m2,rho2m2,omm2,rhomg,rho1mg,rho2mg
80 common /anom/amrop,gamrop,sig,amrop_2,amropg
82 anom_bwg = (dcmplx(1.d0,0.d0)/dcmplx(rhom2-q1_2,-rhomg)
83 1 + dcmplx(sig,0.d0)/dcmplx(amrop_2-q1_2,-amropg) )
84 2 * dcmplx(1.d0,0.d0)/dcmplx(omm2-q2_2,-ommg)
88 complex*16 function bwga1(q1_2)
89 implicit real*8 (a-h,o-z)
91 common /cbwga1/ a1m2,con
94 bwga1 = dcmplx(a1m2,0.d0)/dcmplx(a1m2-q1_2,-ggm)
99 real*8 function gfun8(q1_2)
100 implicit real*8 (a-h,o-z)
102 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
103 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
105 if(q1_2.gt.((rhom+pim)**2))
then 106 gfun8 = q1_2*1.623d0 + 10.38d0 - 9.32d0/q1_2 + 0.65d0/q1_2**2
108 c1 = q1_2 - 9.d0*pim**2
109 gfun8 = 4.1d0 *c1**3 *(1.d0 - 3.3d0*c1 + 5.8d0*c1**2)
115 complex*16 function bwgrho(q1_2)
116 implicit real*8 (a-h,o-z)
118 complex*16 cbw,cbw1,cbw2,cbwo
120 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
121 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
122 common /cbwgrho/ rhom2,rho1m2,rho2m2,omm2,rhomg,rho1mg,rho2mg,ommg
124 c2 = 4.d0*pim**2/q1_2
128 gamrho = rhomg*sqrt(c1*((1.d0-c2)/(c1-c2))**3)
130 gamrho1 = rho1mg*sqrt(c1*((1.d0-c2)/(c1-c2))**3)
132 gamrho2 = rho2mg*sqrt(c1*((1.d0-c2)/(c1-c2))**3)
134 gamom = ommg*sqrt(c1*((1.d0-c2)/(c1-c2))**3)
142 cbw = dcmplx(rhom2,0.d0)/dcmplx(rhom2-q1_2,-gamrho)
143 cbw1 = dcmplx(rho1m2,0.d0)/dcmplx(rho1m2-q1_2,-gamrho1)
144 cbw2 = dcmplx(rho2m2,0.d0)/dcmplx(rho2m2-q1_2,-gamrho2)
145 cbwo = dcmplx(omm2,0.d0)/dcmplx(omm2-q1_2,-gamom)
146 bwgrho = ( cbw *(1.d0+aa*cbwo)/(1.d0+aa)
147 1 + bb1*cbw1+bb2*cbw2)/(1.d0+bb1+bb2)
152 complex*16 function bwgrho_t(q1_2)
153 implicit real*8 (a-h,o-z)
155 complex*16 cbw,cbw1,cbw2,cbwo
157 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
158 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
159 common /cbwgrho/ rhom2,rho1m2,rho2m2,omm2,rhomg,rho1mg,rho2mg
161 common /cbwgrho_t/ rho1m2_t,rho1mg_t,beta
163 c2 = 4.d0*pim**2/q1_2
169 gamrho = rhomg*sqrt(c1*((1.d0-c2)/(c1-c2))**3)
175 gamrho1 = rho1mg_t*sqrt(c1*((1.d0-c2)/(c1-c2))**3)
178 cbw = dcmplx(rhom2,0.d0)/dcmplx(rhom2-q1_2,-gamrho)
179 cbw1 = dcmplx(rho1m2,0.d0)/dcmplx(rho1m2-q1_2,-gamrho1)
181 bwgrho_t = (cbw+beta*cbw1)/(1.d0+beta)
186 complex*16 function bwgf0(q1_2)
187 implicit real*8 (a-h,o-z)
189 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
190 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
194 bwgf0 = dcmplx(f0m2,-f0mg)/dcmplx(f0m2-q1_2,-f0mg)
212 subroutine had1(qq2,q1,q2,q3,q4,hadr)
213 implicit real*8 (a-h,o-z)
215 complex*16 hadr(4),hadr1(4),hadr2(4),hadr3(4),hadr4(4)
216 dimension q1(4),q2(4),q3(4),q4(4)
218 call had2(qq2,q1,q2,q3,q4,hadr1)
219 call had2(qq2,q4,q2,q3,q1,hadr2)
220 call had2(qq2,q1,q3,q2,q4,hadr3)
221 call had2(qq2,q4,q3,q2,q1,hadr4)
224 hadr(i) = hadr1(i)+hadr2(i)+hadr3(i)+hadr4(i)
236 subroutine had3(qq2,q1,q2,q3,q4,hadr)
237 implicit real*8 (a-h,o-z)
239 complex*16 hadr(4),hadr1(4),hadr2(4),hadr3(4)
240 dimension q1(4),q2(4),q3(4),q4(4)
242 call had2(qq2,q1,q2,q3,q4,hadr1)
243 call had2(qq2,q1,q3,q2,q4,hadr2)
244 call had2(qq2,q3,q2,q1,q4,hadr3)
247 hadr(i) = (hadr1(i)+hadr2(i)+hadr3(i))*sqrt(2.d0)
260 subroutine had4(qq2,q1,q2,q3,q4,hadr)
262 implicit real*8 (a-h,o-z)
264 complex*16 hadr(4),hadr1(4),hadr2(4)
265 dimension q1(4),q2(4),q3(4),q4(4)
267 call had2(qq2,q3,q1,q2,q4,hadr1)
268 call had2(qq2,q3,q2,q1,q4,hadr2)
271 hadr(i) = (hadr1(i)+hadr2(i))*sqrt(2.d0)
288 subroutine had2(qq2,q1,q2,q3,q4,hadr)
289 implicit real*8 (a-h,o-z)
291 complex*16 hadr(4),cfac(4),tt(4,4,4),ss(4,4,4,4)
292 complex*16 bwga1,bwgrho,bwgrho_t,bwgf0,c0,c5,c6
293 complex*16 c1_t,c2_t,c3_t,c4_t,anom_bwg
294 dimension q1(4),q2(4),q3(4),q4(4),q2m4(4),q3m1(4),q4m1(4),q3m2(4)
295 dimension q123(4),q124(4),qq(4),q3m4(4),q134(4),q234(4)
296 dimension q2p4(4),q1p3(4),q2p3(4),q1p4(4),q1p2(4),q3p4(4)
298 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
299 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
304 q2m4(i) = q2(i)-q4(i)
305 q3m1(i) = q3(i)-q1(i)
306 q3m4(i) = q3(i)-q4(i)
307 q4m1(i) = q4(i)-q1(i)
308 q3m2(i) = q3(i)-q2(i)
309 q2p4(i) = q2(i)+q4(i)
310 q1p3(i) = q1(i)+q3(i)
311 q1p2(i) = q1(i)+q2(i)
312 q2p3(i) = q2(i)+q3(i)
313 q1p4(i) = q1(i)+q4(i)
314 q3p4(i) = q3(i)+q4(i)
315 q123(i) = q2p3(i)+q1(i)
316 q124(i) = q2p4(i)+q1(i)
317 qq(i) = q123(i) + q4(i)
319 q1_2m4 = q1(1)*q2m4(1)-q1(2)*q2m4(2)-q1(3)*q2m4(3)-q1(4)*q2m4(4)
344 c0 = bwgrho(qq2)*coupl
347 c1_t = bwgrho_t(q2p4_2)
348 c2_t = bwgrho_t(q1p3_2)
349 c3_t = bwgrho_t(q2p3_2)
350 c4_t = bwgrho_t(q1p4_2)
355 tt(1,2,4) = c5*c1_t*gam1
356 tt(2,1,4) = c5*c4_t*gam1
357 tt(2,3,1) = c6*c2_t*gam1
358 tt(1,2,3) = c6*c3_t*gam1
360 ss(3,4,1,2) = bwgrho(q3p4_2)*bwgf0(q1p2_2)*gam2
362 cfac(1) = tt(1,2,3) * (-1.d0 - q1_3m2/qmq4_2 )
363 1 + tt(1,2,4) * ( 1.d0 - q1_2m4/qmq3_2 )
364 2 + tt(2,1,4) * ( 3.d0 + q2_4m1/qmq3_2 )
365 3 + tt(2,3,1) * (-3.d0 - q2_3m1/qmq4_2 )
367 cfac(2) = tt(1,2,3) * (-3.d0 - q1_3m2/qmq4_2 )
368 1 + tt(1,2,4) * ( 3.d0 - q1_2m4/qmq3_2 )
369 2 + tt(2,1,4) * ( 1.d0 + q2_4m1/qmq3_2 )
370 3 + tt(2,3,1) * (-1.d0 - q2_3m1/qmq4_2 )
372 cfac(3) = tt(1,2,3) * ( 1.d0 - q1_3m2/qmq4_2 )
373 1 + tt(1,2,4) * ( 1.d0 + q1_2m4/qmq3_2 )
374 2 + tt(2,1,4) * ( 1.d0 - q2_4m1/qmq3_2 )
375 3 + tt(2,3,1) * ( 1.d0 - q2_3m1/qmq4_2 )
379 1 *(1.d0 -2.d0/qq2*(q_q4*q1_3m2/qmq4_2 +q1p4_3m2) +q1_3m2/qmq4_2 )
381 3 *(-1.d0-2.d0/qq2*(q1_2m4/qmq3_2*q_q3 +q1p3_2m4) +q1_2m4/qmq3_2 )
383 5 *(-1.d0+2.d0/qq2*(q_q3*q2_4m1/qmq3_2 +q2p3_4m1) -q2_4m1/qmq3_2 )
385 7 *(1.d0 -2.d0/qq2*(q2_3m1/qmq4_2*q_q4 +q2p4_3m1) +q2_3m1/qmq4_2 )
386 8 +3.d0*ss(3,4,1,2)/qq2*q1p2_3m4
393 hadr(i) = q1(i) *cfac(1) + q2(i)*cfac(2)
394 1 + q3m4(i)*cfac(3) + qq(i)*cfac(4)
399 fac3 = sgo * 1475.98d0*12.924d0 * 0.266d0 *rhom**2
404 q134(i) = q1p3(i)+q4(i)
405 q234(i) = q2p4(i)+q3(i)
408 q1_134 = q1(1)*q134(1)-q1(2)*q134(2)-q1(3)*q134(3)-q1(4)*q134(4)
423 cfac(1) = anom_bwg(qq2,q134_2) *(q3_134*q24 -q4_134*q23)
424 cfac(2) = anom_bwg(qq2,q234_2) *(q3_234*q14 -q4_234*q13)
425 cfac(3) = anom_bwg(qq2,q134_2) *(q4_134*q12 -q1_134*q24)
426 1 + anom_bwg(qq2,q234_2) *(q4_234*q12 -q2_234*q14)
427 cfac(4) = anom_bwg(qq2,q134_2) *(q1_134*q23 -q3_134*q12)
428 1 + anom_bwg(qq2,q234_2) *(q2_234*q13 -q3_234*q12)
431 hadr(i) = hadr(i) + fac3* (q1(i)*cfac(1) + q2(i)*cfac(2)
432 1 + q3(i)*cfac(3) + q4(i)*cfac(4) )
448 subroutine had2_om(qq2,q1,q2,q3,q4,hadr)
449 implicit real*8 (a-h,o-z)
451 complex*16 hadr(4),cfac(4),tt(4,4,4),ss(4,4,4,4)
452 complex*16 bwga1,bwgrho,bwgrho_t,bwgf0,c0,c5,c6
453 complex*16 c1_t,c2_t,c3_t,c4_t,anom_bwg
454 dimension q1(4),q2(4),q3(4),q4(4),q2m4(4),q3m1(4),q4m1(4),q3m2(4)
455 dimension q123(4),q124(4),qq(4),q3m4(4),q134(4),q234(4)
456 dimension q2p4(4),q1p3(4),q2p3(4),q1p4(4),q1p2(4),q3p4(4)
458 common /had_par/ gam1,gam2,coupl,a1m,a1g,rhom,rhog,rho1m,rho1g
459 1 ,rho2m,rho2g,omm,omg,aa,bb1,bb2,f0m,f0g,pim,sgo
464 q2m4(i) = q2(i)-q4(i)
465 q3m1(i) = q3(i)-q1(i)
466 q3m4(i) = q3(i)-q4(i)
467 q4m1(i) = q4(i)-q1(i)
468 q3m2(i) = q3(i)-q2(i)
469 q2p4(i) = q2(i)+q4(i)
470 q1p3(i) = q1(i)+q3(i)
471 q1p2(i) = q1(i)+q2(i)
472 q2p3(i) = q2(i)+q3(i)
473 q1p4(i) = q1(i)+q4(i)
474 q3p4(i) = q3(i)+q4(i)
475 q123(i) = q2p3(i)+q1(i)
476 q124(i) = q2p4(i)+q1(i)
477 qq(i) = q123(i) + q4(i)
479 q1_2m4 = q1(1)*q2m4(1)-q1(2)*q2m4(2)-q1(3)*q2m4(3)-q1(4)*q2m4(4)
507 fac3 = sgo * 1475.98d0*12.924d0 * 0.266d0 *rhom**2
512 q134(i) = q1p3(i)+q4(i)
513 q234(i) = q2p4(i)+q3(i)
516 q1_134 = q1(1)*q134(1)-q1(2)*q134(2)-q1(3)*q134(3)-q1(4)*q134(4)
531 cfac(1) = anom_bwg(qq2,q134_2) *(q3_134*q24 -q4_134*q23)
532 cfac(2) = anom_bwg(qq2,q234_2) *(q3_234*q14 -q4_234*q13)
533 cfac(3) = anom_bwg(qq2,q134_2) *(q4_134*q12 -q1_134*q24)
534 1 + anom_bwg(qq2,q234_2) *(q4_234*q12 -q2_234*q14)
535 cfac(4) = anom_bwg(qq2,q134_2) *(q1_134*q23 -q3_134*q12)
536 1 + anom_bwg(qq2,q234_2) *(q2_234*q13 -q3_234*q12)
539 hadr(i) = fac3* (q1(i)*cfac(1) + q2(i)*cfac(2)
540 1 + q3(i)*cfac(3) + q4(i)*cfac(4) )