3 SUBROUTINE usx(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.2 .AND. i4.EQ.-1)
CALL usx_udx_h(p,h1,h2,ans)
17 IF(i3.EQ.2 .AND. i4.EQ.-3)
CALL usx_usx_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.-1)
CAll usx_cdx_h(p,h1,h2,ans)
19 IF(i3.EQ.4 .AND. i4.EQ.-3)
CALL usx_udx_h(p,h1,h2,ans)
20 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
21 CALL usx_udx_h(p,h1,h2,ans1)
22 CALL usx_usx_h(p,h1,h2,ans2)
23 CALL usx_cdx_h(p,h1,h2,ans3)
24 CALL usx_udx_h(p,h1,h2,ans4)
25 ans=ans1+ans2+ans3+ans4
27 ELSE IF(key.EQ.0)
THEN
28 IF(i3.EQ.2 .AND. i4.EQ.-1)
CALL usx_udx_noh(p,h1,h2,ans)
29 IF(i3.EQ.2 .AND. i4.EQ.-3)
CALL usx_usx_noh(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.-1)
CAll usx_cdx_noh(p,h1,h2,ans)
31 IF(i3.EQ.4 .AND. i4.EQ.-3)
CALL usx_udx_noh(p,h1,h2,ans)
32 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
33 CALL usx_udx_noh(p,h1,h2,ans1)
34 CALL usx_usx_noh(p,h1,h2,ans2)
35 CALL usx_cdx_noh(p,h1,h2,ans3)
36 CALL usx_udx_noh(p,h1,h2,ans4)
37 ans=ans1+ans2+ans3+ans4
40 WRITE(*,*)
'NOT FINISHED'
50 SUBROUTINE usx_udx_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_usx_udx_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_usx_udx_h(p ,h1,h2,nhel(1,ihel),jc(1))
171 REAL*8 FUNCTION matrix_usx_udx_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_usx_udx_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
235 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
236 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
237 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
238 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
239 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
240 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
241 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
242 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
243 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
245 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
252 ztemp = ztemp + cf(j,i)*jamp(j)
254 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
257 matrix_usx_udx_h=matrix
264 SUBROUTINE usx_usx_h(P,H1,H2,ANS)
284 parameter(nexternal=6)
290 real*8 p(0:3,nexternal),ans
295 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
297 real*8 matrix_usx_usx_h
299 INTEGER JC(NEXTERNAL)
300 LOGICAL GOODHEL(NCOMB)
302 DATA goodhel/ncomb*.false./
303 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
304 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
305 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
306 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
307 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
308 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
309 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
310 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
311 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
312 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
313 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
314 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
315 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
316 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
317 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
318 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
319 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
320 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
321 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
322 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
323 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
324 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
325 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
326 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
327 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
328 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
329 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
330 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
331 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
332 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
333 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
334 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
335 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
336 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
337 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
338 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
339 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
340 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
341 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
342 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
343 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
344 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
345 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
346 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
347 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
348 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
349 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
350 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
351 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
352 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
353 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
354 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
355 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
356 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
357 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
358 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
359 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
360 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
361 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
362 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
363 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
364 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
365 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
366 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
376 t=matrix_usx_usx_h(p ,h1,h2,nhel(1,ihel),jc(1))
383 REAL*8 FUNCTION matrix_usx_usx_h(P,H1,H2,NHEL,IC)
402 parameter(nexternal=6)
403 INTEGER NWAVEFUNCS, NCOLOR
404 parameter(nwavefuncs=7, ncolor=2)
408 parameter(imag1=(0d0,1d0))
412 real*8 p(0:3,nexternal)
413 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
419 real*8 denom(ncolor), cf(ncolor,ncolor)
420 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
421 COMPLEX*16 W(18,NWAVEFUNCS)
423 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
432 DATA (cf(i, 1),i= 1, 2) / 9, 3/
435 DATA (cf(i, 2),i= 1, 2) / 3, 9/
440 matrix_usx_usx_h=0.d0
441 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
442 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
448 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
449 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
450 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
451 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
452 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
453 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
454 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
455 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
456 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,5))
458 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_72,amp(1))
459 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
460 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
462 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
470 ztemp = ztemp + cf(j,i)*jamp(j)
472 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
474 matrix_usx_usx_h=matrix
481 SUBROUTINE usx_cdx_h(P,H1,H2,ANS)
501 parameter(nexternal=6)
507 real*8 p(0:3,nexternal),ans
512 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
514 real*8 matrix_usx_cdx_h
516 INTEGER JC(NEXTERNAL)
517 LOGICAL GOODHEL(NCOMB)
519 DATA goodhel/ncomb*.false./
520 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
521 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
522 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
523 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
524 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
525 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
526 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
527 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
528 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
529 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
530 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
531 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
532 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
533 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
534 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
535 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
536 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
537 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
538 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
539 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
540 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
541 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
542 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
543 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
544 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
545 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
546 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
547 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
548 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
549 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
550 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
551 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
552 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
553 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
554 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
555 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
556 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
557 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
558 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
559 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
560 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
561 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
562 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
563 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
564 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
565 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
566 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
567 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
568 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
569 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
570 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
571 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
572 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
573 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
574 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
575 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
576 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
577 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
578 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
579 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
580 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
581 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
582 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
583 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
593 t=matrix_usx_cdx_h(p ,h1,h2,nhel(1,ihel),jc(1))
600 REAL*8 FUNCTION matrix_usx_cdx_h(P,H1,H2,NHEL,IC)
619 parameter(nexternal=6)
620 INTEGER NWAVEFUNCS, NCOLOR
621 parameter(nwavefuncs=7, ncolor=1)
625 parameter(imag1=(0d0,1d0))
629 real*8 p(0:3,nexternal)
630 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
636 real*8 denom(ncolor), cf(ncolor,ncolor)
637 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
638 COMPLEX*16 W(18,NWAVEFUNCS)
640 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
649 DATA (cf(i, 1),i= 1, 1) / 9/
654 matrix_usx_cdx_h=0.d0
655 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
656 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
661 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
662 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
663 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
664 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
665 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
666 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
667 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
668 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
669 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,2))
671 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
678 ztemp = ztemp + cf(j,i)*jamp(j)
680 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
682 matrix_usx_cdx_h=matrix
694 SUBROUTINE usx_udx_noh(P,H1,H2,ANS)
714 parameter(nexternal=6)
720 real*8 p(0:3,nexternal),ans
725 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
727 real*8 matrix_usx_udx_noh
729 INTEGER JC(NEXTERNAL)
730 LOGICAL GOODHEL(NCOMB)
732 DATA goodhel/ncomb*.false./
733 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
734 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
735 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
736 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
737 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
738 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
739 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
740 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
741 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
742 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
743 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
744 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
745 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
746 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
747 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
748 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
749 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
750 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
751 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
752 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
753 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
754 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
755 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
756 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
757 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
758 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
759 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
760 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
761 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
762 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
763 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
764 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
765 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
766 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
767 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
768 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
769 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
770 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
771 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
772 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
773 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
774 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
775 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
776 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
777 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
778 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
779 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
780 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
781 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
782 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
783 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
784 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
785 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
786 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
787 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
788 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
789 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
790 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
791 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
792 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
793 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
794 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
795 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
796 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
806 t=matrix_usx_udx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
813 REAL*8 FUNCTION matrix_usx_udx_noh(P,H1,H2,NHEL,IC)
830 parameter(ngraphs=11)
832 parameter(nexternal=6)
833 INTEGER NWAVEFUNCS, NCOLOR
834 parameter(nwavefuncs=10, ncolor=1)
838 parameter(imag1=(0d0,1d0))
842 real*8 p(0:3,nexternal)
843 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
849 real*8 denom(ncolor), cf(ncolor,ncolor)
850 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
851 COMPLEX*16 W(18,NWAVEFUNCS)
853 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
862 DATA (cf(i, 1),i= 1, 1) / 9/
867 matrix_usx_udx_noh=0.d0
868 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
869 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
874 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
875 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
876 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
877 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
878 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
879 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
880 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
881 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
882 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
884 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
885 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
887 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
888 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
890 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
891 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
893 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
894 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
896 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
898 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
900 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
901 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
903 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
904 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
906 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
908 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
909 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
911 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
912 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
913 $ -amp(9)-amp(10)-amp(11)
919 ztemp = ztemp + cf(j,i)*jamp(j)
921 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
923 matrix_usx_udx_noh=matrix
930 SUBROUTINE usx_usx_noh(P,H1,H2,ANS)
949 parameter(nexternal=6)
955 real*8 p(0:3,nexternal),ans
960 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
962 real*8 matrix_usx_usx_noh
964 INTEGER JC(NEXTERNAL)
965 LOGICAL GOODHEL(NCOMB)
967 DATA goodhel/ncomb*.false./
968 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
969 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
970 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
971 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
972 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
973 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
974 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
975 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
976 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
977 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
978 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
979 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
980 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
981 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
982 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
983 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
984 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
985 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
986 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
987 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
988 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
989 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
990 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
991 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
992 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
993 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
994 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
995 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
996 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
997 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
998 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
999 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1000 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1001 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1002 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1003 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1004 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1005 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1006 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1007 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1008 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1009 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1010 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1011 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1012 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1013 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1014 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1015 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1016 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1017 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1018 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1019 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1020 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1021 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1022 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1023 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1024 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1025 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1026 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1027 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1028 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1029 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1030 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1031 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1041 t=matrix_usx_usx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1048 REAL*8 FUNCTION matrix_usx_usx_noh(P,H1,H2,NHEL,IC)
1064 parameter(ngraphs=43)
1066 parameter(nexternal=6)
1067 INTEGER NWAVEFUNCS, NCOLOR
1068 parameter(nwavefuncs=15, ncolor=2)
1072 parameter(imag1=(0d0,1d0))
1076 real*8 p(0:3,nexternal)
1077 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1083 real*8 denom(ncolor), cf(ncolor,ncolor)
1084 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1085 COMPLEX*16 W(18,NWAVEFUNCS)
1086 COMPLEX*16 DUM0,DUM1
1087 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1096 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1099 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1104 matrix_usx_usx_noh=0.d0
1105 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1106 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1111 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1112 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1113 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1114 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1115 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1116 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1117 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1118 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1119 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1121 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1122 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1124 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1125 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1127 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1128 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,11))
1130 CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_1,amp(4))
1131 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,12))
1133 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(5))
1135 CALL ffv2_3_0(w(1,4),w(1,11),w(1,10),gc_50,gc_58,amp(6))
1137 CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1138 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,12))
1139 CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,11))
1140 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1142 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1143 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1145 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1146 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,14))
1148 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1150 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1151 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,13))
1152 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1154 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1155 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1157 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1159 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1161 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1162 CALL ffv1_1(w(1,2),w(1,12),gc_1,zero,zero,w(1,15))
1164 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(16))
1165 CALL ffv1_2(w(1,4),w(1,12),gc_1,zero,zero,w(1,5))
1167 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(17))
1169 CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(18))
1171 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(19))
1172 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1173 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1175 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(20))
1176 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1178 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_1,amp(21))
1180 CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(22))
1182 CALL ffv2_3_0(w(1,12),w(1,2),w(1,10),gc_50,gc_58,amp(23))
1183 CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,12))
1185 CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_1,amp(24))
1186 CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,15))
1188 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_1,amp(25))
1190 CALL ffv2_3_0(w(1,4),w(1,12),w(1,10),gc_50,gc_58,amp(26))
1192 CALL ffv2_3_0(w(1,15),w(1,2),w(1,10),gc_50,gc_58,amp(27))
1193 CALL ffv1_2(w(1,1),w(1,11),gc_2,zero,zero,w(1,15))
1195 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(28))
1196 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,12))
1198 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_2,amp(29))
1200 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(30))
1201 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,15))
1203 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_2,amp(31))
1204 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1205 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1207 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_2,amp(32))
1209 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1211 CALL ffv2_5_0(w(1,4),w(1,3),w(1,10),gc_51,gc_58,amp(34))
1213 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1214 CALL ffv2_5_2(w(1,1),w(1,14),gc_51,gc_58,zero,zero,w(1,11))
1216 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_2,amp(36))
1218 CALL ffv2_5_0(w(1,12),w(1,3),w(1,14),gc_51,gc_58,amp(37))
1220 CALL ffv2_5_0(w(1,11),w(1,3),w(1,10),gc_51,gc_58,amp(38))
1222 CALL ffv2_5_0(w(1,15),w(1,3),w(1,14),gc_51,gc_58,amp(39))
1223 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,14))
1225 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_1,amp(40))
1227 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_101,amp(41))
1229 CALL ffv2_3_0(w(1,14),w(1,2),w(1,10),gc_50,gc_58,amp(42))
1231 CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_101,amp(43))
1232 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1233 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1234 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1235 $ -amp(40)-amp(41)-amp(42)-amp(43)
1236 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1237 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1238 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1239 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1240 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1247 ztemp = ztemp + cf(j,i)*jamp(j)
1249 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1251 matrix_usx_usx_noh=matrix
1258 SUBROUTINE usx_cdx_noh(P,H1,H2,ANS)
1277 parameter(nexternal=6)
1279 parameter( ncomb=64)
1283 real*8 p(0:3,nexternal),ans
1288 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1290 real*8 matrix_usx_cdx_noh
1291 INTEGER IHEL,IDEN, I
1292 INTEGER JC(NEXTERNAL)
1293 LOGICAL GOODHEL(NCOMB)
1295 DATA goodhel/ncomb*.false./
1296 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1297 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1298 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1299 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1300 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1301 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1302 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1303 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1304 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1305 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1306 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1307 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1308 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1309 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1310 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1311 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1312 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1313 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1314 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1315 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1316 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1317 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1318 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1319 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1320 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1321 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1322 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1323 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1324 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1325 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1326 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1327 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1328 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1329 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1330 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1331 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1332 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1333 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1334 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1335 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1336 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1337 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1338 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1339 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1340 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1341 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1342 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1343 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1344 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1345 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1346 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1347 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1348 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1349 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1350 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1351 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1352 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1353 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1354 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1355 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1356 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1357 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1358 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1359 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1369 t=matrix_usx_cdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1377 REAL*8 FUNCTION matrix_usx_cdx_noh(P,H1,H2,NHEL,IC)
1393 parameter(ngraphs=11)
1395 parameter(nexternal=6)
1396 INTEGER NWAVEFUNCS, NCOLOR
1397 parameter(nwavefuncs=10, ncolor=1)
1401 parameter(imag1=(0d0,1d0))
1405 real*8 p(0:3,nexternal)
1406 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1412 real*8 denom(ncolor), cf(ncolor,ncolor)
1413 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1414 COMPLEX*16 W(18,NWAVEFUNCS)
1415 COMPLEX*16 DUM0,DUM1
1416 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1425 DATA (cf(i, 1),i= 1, 1) / 9/
1430 matrix_usx_cdx_noh=0.d0
1431 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1432 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1437 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1438 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1439 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1440 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1441 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1442 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1443 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1444 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1445 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1447 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1448 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1450 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1451 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1453 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1454 CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,6))
1456 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1457 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
1459 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1461 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1463 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1464 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1466 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1467 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1469 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
1471 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1472 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1474 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
1475 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1476 $ -amp(9)-amp(10)-amp(11)
1482 ztemp = ztemp + cf(j,i)*jamp(j)
1484 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1486 matrix_usx_cdx_noh=matrix