3 SUBROUTINE dcx(P,I3,I4,H1,H2,KEY,ANS)
6 INTEGER I3,I4, H1,H2, KEY
7 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
16 IF(i3.EQ.1 .AND. i4.EQ.-2)
CALL dcx_dux_h(p,h1,h2,ans)
17 IF(i3.EQ.3 .AND. i4.EQ.-2)
CALL dcx_sux_h(p,h1,h2,ans)
18 IF(i3.EQ.1 .AND. i4.EQ.-4)
CAll dcx_dcx_h(p,h1,h2,ans)
19 IF(i3.EQ.3 .AND. i4.EQ.-4)
CALL dcx_dux_h(p,h1,h2,ans)
20 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
21 CALL dcx_dux_h(p,h1,h2,ans1)
22 CALL dcx_sux_h(p,h1,h2,ans2)
23 CALL dcx_dcx_h(p,h1,h2,ans3)
24 CALL dcx_dux_h(p,h1,h2,ans4)
25 ans=ans1+ans2+ans3+ans4
27 ELSE IF(key.EQ.0)
THEN
28 IF(i3.EQ.1 .AND. i4.EQ.-2)
CALL dcx_dux_noh(p,h1,h2,ans)
29 IF(i3.EQ.3 .AND. i4.EQ.-2)
CALL dcx_sux_noh(p,h1,h2,ans)
30 IF(i3.EQ.1 .AND. i4.EQ.-4)
CAll dcx_dcx_noh(p,h1,h2,ans)
31 IF(i3.EQ.3 .AND. i4.EQ.-4)
CALL dcx_dux_noh(p,h1,h2,ans)
32 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
33 CALL dcx_dux_noh(p,h1,h2,ans1)
34 CALL dcx_sux_noh(p,h1,h2,ans2)
35 CALL dcx_dcx_noh(p,h1,h2,ans3)
36 CALL dcx_dux_noh(p,h1,h2,ans4)
37 ans=ans1+ans2+ans3+ans4
40 WRITE(*,*)
'NOT FINISHED'
50 SUBROUTINE dcx_dux_h(P,H1,H2,ANS)
72 parameter(nexternal=6)
78 real*8 p(0:3,nexternal),ans
83 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
85 real*8 matrix_dcx_dux_h
88 LOGICAL GOODHEL(NCOMB)
90 DATA goodhel/ncomb*.false./
91 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
92 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
93 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
94 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
95 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
96 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
97 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
98 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
99 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
100 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
101 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
102 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
103 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
104 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
105 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
106 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
107 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
108 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
109 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
110 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
111 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
112 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
113 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
114 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
115 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
116 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
117 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
118 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
119 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
120 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
121 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
122 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
123 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
124 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
125 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
126 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
127 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
128 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
129 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
130 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
131 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
132 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
133 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
134 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
135 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
136 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
137 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
138 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
139 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
140 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
141 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
142 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
143 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
144 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
145 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
146 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
147 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
148 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
149 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
150 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
151 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
152 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
153 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
154 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
164 t=matrix_dcx_dux_h(p ,h1,h2,nhel(1,ihel),jc(1))
171 REAL*8 FUNCTION matrix_dcx_dux_h(P,H1,H2,NHEL,IC)
192 parameter(nexternal=6)
193 INTEGER NWAVEFUNCS, NCOLOR
194 parameter(nwavefuncs=7, ncolor=1)
198 parameter(imag1=(0d0,1d0))
202 real*8 p(0:3,nexternal)
203 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
209 real*8 denom(ncolor), cf(ncolor,ncolor)
210 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
211 COMPLEX*16 W(18,NWAVEFUNCS)
213 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
222 DATA (cf(i, 1),i= 1, 1) / 9/
227 matrix_dcx_dux_h=0.d0
228 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
229 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
234 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
235 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
236 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
237 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
238 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
239 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
240 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
241 CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,6))
242 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
244 CALL vvs1_0(w(1,2),w(1,6),w(1,7),gc_72,amp(1))
251 ztemp = ztemp + cf(j,i)*jamp(j)
253 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
256 matrix_dcx_dux_h=matrix
263 SUBROUTINE dcx_dcx_h(P,H1,H2,ANS)
283 parameter(nexternal=6)
289 real*8 p(0:3,nexternal),ans
294 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
296 real*8 matrix_dcx_dcx_h
298 INTEGER JC(NEXTERNAL)
299 LOGICAL GOODHEL(NCOMB)
301 DATA goodhel/ncomb*.false./
302 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
303 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
304 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
305 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
306 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
307 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
308 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
309 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
310 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
311 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
312 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
313 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
314 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
315 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
316 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
317 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
318 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
319 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
320 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
321 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
322 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
323 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
324 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
325 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
326 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
327 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
328 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
329 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
330 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
331 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
332 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
333 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
334 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
335 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
336 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
337 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
338 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
339 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
340 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
341 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
342 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
343 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
344 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
345 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
346 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
347 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
348 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
349 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
350 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
351 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
352 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
353 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
354 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
355 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
356 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
357 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
358 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
359 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
360 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
361 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
362 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
363 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
364 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
365 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
375 t=matrix_dcx_dcx_h(p ,h1,h2,nhel(1,ihel),jc(1))
382 REAL*8 FUNCTION matrix_dcx_dcx_h(P,H1,H2,NHEL,IC)
401 parameter(nexternal=6)
402 INTEGER NWAVEFUNCS, NCOLOR
403 parameter(nwavefuncs=7, ncolor=2)
407 parameter(imag1=(0d0,1d0))
411 real*8 p(0:3,nexternal)
412 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
418 real*8 denom(ncolor), cf(ncolor,ncolor)
419 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
420 COMPLEX*16 W(18,NWAVEFUNCS)
422 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
431 DATA (cf(i, 1),i= 1, 2) / 9, 3/
434 DATA (cf(i, 2),i= 1, 2) / 3, 9/
439 matrix_dcx_dcx_h=0.d0
440 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
441 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
446 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
447 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
448 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
449 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
450 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
451 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
452 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
453 CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,6))
454 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,5))
456 CALL vvs1_0(w(1,5),w(1,6),w(1,7),gc_72,amp(1))
457 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
458 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,3))
460 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
468 ztemp = ztemp + cf(j,i)*jamp(j)
470 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
472 matrix_dcx_dcx_h=matrix
479 SUBROUTINE dcx_sux_h(P,H1,H2,ANS)
499 parameter(nexternal=6)
505 real*8 p(0:3,nexternal),ans
510 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
512 real*8 matrix_dcx_sux_h
514 INTEGER JC(NEXTERNAL)
515 LOGICAL GOODHEL(NCOMB)
517 DATA goodhel/ncomb*.false./
518 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
519 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
520 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
521 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
522 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
523 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
524 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
525 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
526 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
527 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
528 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
529 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
530 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
531 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
532 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
533 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
534 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
535 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
536 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
537 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
538 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
539 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
540 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
541 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
542 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
543 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
544 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
545 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
546 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
547 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
548 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
549 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
550 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
551 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
552 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
553 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
554 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
555 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
556 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
557 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
558 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
559 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
560 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
561 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
562 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
563 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
564 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
565 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
566 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
567 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
568 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
569 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
570 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
571 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
572 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
573 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
574 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
575 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
576 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
577 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
578 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
579 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
580 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
581 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
591 t=matrix_dcx_sux_h(p ,h1,h2,nhel(1,ihel),jc(1))
598 REAL*8 FUNCTION matrix_dcx_sux_h(P,H1,H2,NHEL,IC)
617 parameter(nexternal=6)
618 INTEGER NWAVEFUNCS, NCOLOR
619 parameter(nwavefuncs=7, ncolor=1)
623 parameter(imag1=(0d0,1d0))
627 real*8 p(0:3,nexternal)
628 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
634 real*8 denom(ncolor), cf(ncolor,ncolor)
635 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
636 COMPLEX*16 W(18,NWAVEFUNCS)
638 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
647 DATA (cf(i, 1),i= 1, 1) / 9/
652 matrix_dcx_sux_h=0.d0
653 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
654 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
659 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
660 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
661 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
662 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
663 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
664 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
665 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
666 CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,6))
667 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,2))
669 CALL vvs1_0(w(1,2),w(1,6),w(1,7),gc_72,amp(1))
676 ztemp = ztemp + cf(j,i)*jamp(j)
678 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
680 matrix_dcx_sux_h=matrix
692 SUBROUTINE dcx_dux_noh(P,H1,H2,ANS)
712 parameter(nexternal=6)
718 real*8 p(0:3,nexternal),ans
723 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
725 real*8 matrix_dcx_dux_noh
727 INTEGER JC(NEXTERNAL)
728 LOGICAL GOODHEL(NCOMB)
730 DATA goodhel/ncomb*.false./
731 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
732 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
733 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
734 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
735 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
736 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
737 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
738 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
739 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
740 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
741 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
742 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
743 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
744 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
745 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
746 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
747 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
748 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
749 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
750 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
751 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
752 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
753 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
754 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
755 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
756 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
757 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
758 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
759 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
760 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
761 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
762 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
763 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
764 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
765 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
766 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
767 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
768 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
769 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
770 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
771 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
772 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
773 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
774 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
775 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
776 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
777 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
778 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
779 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
780 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
781 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
782 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
783 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
784 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
785 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
786 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
787 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
788 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
789 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
790 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
791 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
792 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
793 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
794 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
804 t=matrix_dcx_dux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
811 REAL*8 FUNCTION matrix_dcx_dux_noh(P,H1,H2,NHEL,IC)
828 parameter(ngraphs=11)
830 parameter(nexternal=6)
831 INTEGER NWAVEFUNCS, NCOLOR
832 parameter(nwavefuncs=10, ncolor=1)
836 parameter(imag1=(0d0,1d0))
840 real*8 p(0:3,nexternal)
841 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
847 real*8 denom(ncolor), cf(ncolor,ncolor)
848 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
849 COMPLEX*16 W(18,NWAVEFUNCS)
851 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
860 DATA (cf(i, 1),i= 1, 1) / 9/
865 matrix_dcx_dux_noh=0.d0
866 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
867 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
872 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
873 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
874 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
875 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
876 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
877 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
878 CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,7))
879 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
880 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
882 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
883 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
885 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
886 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
888 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
889 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
891 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
892 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
894 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
896 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
898 CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(7))
899 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
901 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(8))
902 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,3))
904 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_44,amp(9))
906 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
907 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
909 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_44,amp(11))
910 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
911 $ -amp(9)-amp(10)-amp(11)
917 ztemp = ztemp + cf(j,i)*jamp(j)
919 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
922 matrix_dcx_dux_noh=matrix
929 SUBROUTINE dcx_dcx_noh(P,H1,H2,ANS)
948 parameter(nexternal=6)
954 real*8 p(0:3,nexternal),ans
959 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
961 real*8 matrix_dcx_dcx_noh
963 INTEGER JC(NEXTERNAL)
964 LOGICAL GOODHEL(NCOMB)
966 DATA goodhel/ncomb*.false./
967 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
968 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
969 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
970 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
971 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
972 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
973 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
974 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
975 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
976 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
977 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
978 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
979 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
980 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
981 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
982 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
983 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
984 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
985 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
986 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
987 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
988 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
989 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
990 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
991 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
992 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
993 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
994 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
995 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
996 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
997 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
998 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
999 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1000 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1001 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1002 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1003 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1004 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1005 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1006 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1007 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1008 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1009 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1010 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1011 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1012 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1013 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1014 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1015 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1016 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1017 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1018 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1019 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1020 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1021 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1022 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1023 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1024 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1025 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1026 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1027 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1028 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1029 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1030 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1040 t=matrix_dcx_dcx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1047 REAL*8 FUNCTION matrix_dcx_dcx_noh(P,H1,H2,NHEL,IC)
1063 parameter(ngraphs=43)
1065 parameter(nexternal=6)
1066 INTEGER NWAVEFUNCS, NCOLOR
1067 parameter(nwavefuncs=15, ncolor=2)
1071 parameter(imag1=(0d0,1d0))
1075 real*8 p(0:3,nexternal)
1076 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1082 real*8 denom(ncolor), cf(ncolor,ncolor)
1083 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1084 COMPLEX*16 W(18,NWAVEFUNCS)
1085 COMPLEX*16 DUM0,DUM1
1086 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1095 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1098 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1103 matrix_dcx_dcx_noh=0.d0
1104 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1105 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1110 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1111 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1112 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1113 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1114 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1115 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1116 CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,7))
1117 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1118 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1120 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1121 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1123 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1124 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1126 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1127 CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,11))
1129 CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_2,amp(4))
1130 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,12))
1132 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_1,amp(5))
1134 CALL ffv2_5_0(w(1,4),w(1,11),w(1,10),gc_51,gc_58,amp(6))
1136 CALL ffv2_3_0(w(1,12),w(1,3),w(1,10),gc_50,gc_58,amp(7))
1137 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,12))
1138 CALL ffv1p0_3(w(1,4),w(1,2),gc_2,zero,zero,w(1,11))
1139 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1141 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1142 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1144 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1145 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,14))
1147 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1149 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1150 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,13))
1151 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1153 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1154 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1156 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1158 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1160 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1161 CALL ffv1_1(w(1,2),w(1,12),gc_2,zero,zero,w(1,15))
1163 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_2,amp(16))
1164 CALL ffv1_2(w(1,4),w(1,12),gc_2,zero,zero,w(1,5))
1166 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(17))
1168 CALL ffv2_5_0(w(1,4),w(1,15),w(1,10),gc_51,gc_58,amp(18))
1170 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(19))
1171 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1172 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1174 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_2,amp(20))
1175 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1177 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_2,amp(21))
1179 CALL ffv2_5_0(w(1,4),w(1,15),w(1,10),gc_51,gc_58,amp(22))
1181 CALL ffv2_5_0(w(1,12),w(1,2),w(1,10),gc_51,gc_58,amp(23))
1182 CALL ffv2_5_1(w(1,2),w(1,13),gc_51,gc_58,zero,zero,w(1,12))
1184 CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_2,amp(24))
1185 CALL ffv2_5_2(w(1,4),w(1,13),gc_51,gc_58,zero,zero,w(1,15))
1187 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_2,amp(25))
1189 CALL ffv2_5_0(w(1,4),w(1,12),w(1,10),gc_51,gc_58,amp(26))
1191 CALL ffv2_5_0(w(1,15),w(1,2),w(1,10),gc_51,gc_58,amp(27))
1192 CALL ffv1_2(w(1,1),w(1,11),gc_1,zero,zero,w(1,15))
1194 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_1,amp(28))
1195 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,12))
1197 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_1,amp(29))
1199 CALL ffv2_3_0(w(1,15),w(1,3),w(1,10),gc_50,gc_58,amp(30))
1200 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,15))
1202 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_1,amp(31))
1203 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1204 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1206 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_1,amp(32))
1208 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1210 CALL ffv2_3_0(w(1,4),w(1,3),w(1,10),gc_50,gc_58,amp(34))
1212 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1213 CALL ffv2_3_2(w(1,1),w(1,14),gc_50,gc_58,zero,zero,w(1,11))
1215 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_1,amp(36))
1217 CALL ffv2_3_0(w(1,12),w(1,3),w(1,14),gc_50,gc_58,amp(37))
1219 CALL ffv2_3_0(w(1,11),w(1,3),w(1,10),gc_50,gc_58,amp(38))
1221 CALL ffv2_3_0(w(1,15),w(1,3),w(1,14),gc_50,gc_58,amp(39))
1222 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,14))
1224 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_2,amp(40))
1226 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_44,amp(41))
1228 CALL ffv2_5_0(w(1,14),w(1,2),w(1,10),gc_51,gc_58,amp(42))
1230 CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_44,amp(43))
1231 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1232 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1233 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1234 $ -amp(40)-amp(41)-amp(42)-amp(43)
1235 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1236 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1237 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1238 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1239 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1246 ztemp = ztemp + cf(j,i)*jamp(j)
1248 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1250 matrix_dcx_dcx_noh=matrix
1257 SUBROUTINE dcx_sux_noh(P,H1,H2,ANS)
1276 parameter(nexternal=6)
1278 parameter( ncomb=64)
1282 real*8 p(0:3,nexternal),ans
1287 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1289 real*8 matrix_dcx_sux_noh
1290 INTEGER IHEL,IDEN, I
1291 INTEGER JC(NEXTERNAL)
1292 LOGICAL GOODHEL(NCOMB)
1294 DATA goodhel/ncomb*.false./
1295 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1296 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1297 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1298 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1299 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1300 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1301 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1302 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1303 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1304 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1305 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1306 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1307 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1308 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1309 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1310 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1311 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1312 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1313 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1314 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1315 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1316 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1317 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1318 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1319 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1320 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1321 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1322 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1323 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1324 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1325 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1326 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1327 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1328 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1329 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1330 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1331 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1332 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1333 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1334 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1335 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1336 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1337 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1338 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1339 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1340 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1341 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1342 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1343 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1344 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1345 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1346 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1347 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1348 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1349 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1350 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1351 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1352 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1353 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1354 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1355 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1356 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1357 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1358 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1368 t=matrix_dcx_sux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1375 REAL*8 FUNCTION matrix_dcx_sux_noh(P,H1,H2,NHEL,IC)
1391 parameter(ngraphs=11)
1393 parameter(nexternal=6)
1394 INTEGER NWAVEFUNCS, NCOLOR
1395 parameter(nwavefuncs=10, ncolor=1)
1399 parameter(imag1=(0d0,1d0))
1403 real*8 p(0:3,nexternal)
1404 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1410 real*8 denom(ncolor), cf(ncolor,ncolor)
1411 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1412 COMPLEX*16 W(18,NWAVEFUNCS)
1413 COMPLEX*16 DUM0,DUM1
1414 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1423 DATA (cf(i, 1),i= 1, 1) / 9/
1428 matrix_dcx_sux_noh=0.d0
1429 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1430 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1435 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1436 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1437 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1438 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1439 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1440 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1441 CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,7))
1442 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1443 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1445 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1446 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1448 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1449 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1451 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1452 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,6))
1454 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1455 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
1457 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1459 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1461 CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(7))
1462 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1464 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(8))
1465 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,3))
1467 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_44,amp(9))
1469 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
1470 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1472 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_44,amp(11))
1473 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1474 $ -amp(9)-amp(10)-amp(11)
1480 ztemp = ztemp + cf(j,i)*jamp(j)
1482 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1484 matrix_dcx_sux_noh=matrix