3 SUBROUTINE ssx(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, ans5,ans6,ans7
15 IF(i3.EQ.1 .AND. i4.EQ.-1)
CALL ssx_ddx_h(p,h1,h2,ans)
16 IF(i3.EQ.2 .AND. i4.EQ.-2)
CALL ssx_uux_h(p,h1,h2,ans)
17 IF(i3.EQ.3 .AND. i4.EQ.-3)
CAll ssx_ssx_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.-4)
CALL ssx_ccx_h(p,h1,h2,ans)
19 IF(i3.EQ.2 .AND. i4.EQ.-4)
CALL ssx_ucx_h(p,h1,h2,ans)
20 IF(i3.EQ.4 .AND. i4.EQ.-2)
CALL ssx_cux_h(p,h1,h2,ans)
21 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
22 CALL ssx_ddx_h(p,h1,h2,ans1)
23 CALL ssx_uux_h(p,h1,h2,ans2)
24 CALL ssx_ssx_h(p,h1,h2,ans3)
25 CALL ssx_ccx_h(p,h1,h2,ans4)
26 CALL ssx_ucx_h(p,h1,h2,ans5)
27 CALL ssx_cux_h(p,h1,h2,ans6)
28 ans=ans1+ans2+ans3+ans4+ans5+ans6
30 ELSE IF(key.EQ.0)
THEN
31 IF(i3.EQ.1 .AND. i4.EQ.-1)
CALL ssx_ddx_noh(p,h1,h2,ans)
32 IF(i3.EQ.2 .AND. i4.EQ.-2)
CALL ssx_uux_noh(p,h1,h2,ans)
33 IF(i3.EQ.3 .AND. i4.EQ.-3)
CAll ssx_ssx_noh(p,h1,h2,ans)
34 IF(i3.EQ.4 .AND. i4.EQ.-4)
CALL ssx_ccx_noh(p,h1,h2,ans)
35 IF(i3.EQ.2 .AND. i4.EQ.-4)
CALL ssx_ucx_noh(p,h1,h2,ans)
36 IF(i3.EQ.4 .AND. i4.EQ.-2)
CALL ssx_cux_noh(p,h1,h2,ans)
37 IF(i3.EQ.21.AND.i4.EQ.21)
CALL ssx_gg_noh(p,h1,h2,ans)
38 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
39 CALL ssx_ddx_noh(p,h1,h2,ans1)
40 CALL ssx_uux_noh(p,h1,h2,ans2)
41 CALL ssx_ssx_noh(p,h1,h2,ans3)
42 CALL ssx_ccx_noh(p,h1,h2,ans4)
43 CALL ssx_ucx_noh(p,h1,h2,ans5)
44 CALL ssx_cux_noh(p,h1,h2,ans6)
45 CALL ssx_gg_noh(p,h1,h2,ans7)
46 ans=ans1+ans2+ans3+ans4+ans5+ans6+ans7
50 WRITE(*,*)
'NOT FINISHED'
59 SUBROUTINE ssx_ddx_h(P,H1,H2,ANS)
79 parameter(nexternal=6)
85 real*8 p(0:3,nexternal),ans
90 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
92 real*8 matrix_ssx_ddx_h
95 LOGICAL GOODHEL(NCOMB)
97 DATA goodhel/ncomb*.false./
98 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
99 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
100 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
101 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
102 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
103 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
104 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
105 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
106 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
107 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
108 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
109 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
110 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
111 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
112 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
113 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
114 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
115 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
116 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
117 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
118 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
119 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
120 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
121 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
122 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
123 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
124 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
125 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
126 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
127 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
128 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
129 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
130 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
131 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
132 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
133 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
134 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
135 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
136 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
137 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
138 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
139 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
140 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
141 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
142 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
143 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
144 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
145 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
146 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
147 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
148 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
149 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
150 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
151 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
152 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
153 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
154 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
155 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
156 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
157 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
158 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
159 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
160 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
161 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
171 t=matrix_ssx_ddx_h(p ,h1,h2,nhel(1,ihel),jc(1))
178 REAL*8 FUNCTION matrix_ssx_ddx_h(P,H1,H2,NHEL,IC)
197 parameter(nexternal=6)
198 INTEGER NWAVEFUNCS, NCOLOR
199 parameter(nwavefuncs=7, ncolor=1)
203 parameter(imag1=(0d0,1d0))
207 real*8 p(0:3,nexternal)
208 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
214 real*8 denom(ncolor), cf(ncolor,ncolor)
215 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
216 COMPLEX*16 W(18,NWAVEFUNCS)
218 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
227 DATA (cf(i, 1),i= 1, 1) / 9/
232 matrix_ssx_ddx_h=0.d0
233 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
234 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
239 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
240 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
241 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
242 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
243 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
244 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
245 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
246 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
247 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,2))
249 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_81,amp(1))
256 ztemp = ztemp + cf(j,i)*jamp(j)
258 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
260 matrix_ssx_ddx_h=matrix
267 SUBROUTINE ssx_uux_h(P,H1,H2,ANS)
287 parameter(nexternal=6)
293 real*8 p(0:3,nexternal),ans
298 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
300 real*8 matrix_ssx_uux_h
302 INTEGER JC(NEXTERNAL)
303 LOGICAL GOODHEL(NCOMB)
305 DATA goodhel/ncomb*.false./
306 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
307 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
308 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
309 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
310 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
311 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
312 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
313 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
314 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
315 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
316 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
317 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
318 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
319 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
320 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
321 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
322 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
323 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
324 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
325 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
326 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
327 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
328 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
329 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
330 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
331 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
332 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
333 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
334 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
335 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
336 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
337 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
338 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
339 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
340 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
341 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
342 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
343 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
344 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
345 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
346 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
347 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
348 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
349 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
350 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
351 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
352 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
353 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
354 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
355 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
356 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
357 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
358 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
359 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
360 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
361 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
362 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
363 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
364 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
365 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
366 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
367 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
368 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
369 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
379 t=matrix_ssx_uux_h(p ,h1,h2,nhel(1,ihel),jc(1))
386 REAL*8 FUNCTION matrix_ssx_uux_h(P,H1,H2,NHEL,IC)
405 parameter(nexternal=6)
406 INTEGER NWAVEFUNCS, NCOLOR
407 parameter(nwavefuncs=7, ncolor=2)
411 parameter(imag1=(0d0,1d0))
415 real*8 p(0:3,nexternal)
416 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
422 real*8 denom(ncolor), cf(ncolor,ncolor)
423 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
424 COMPLEX*16 W(18,NWAVEFUNCS)
426 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
435 DATA (cf(i, 1),i= 1, 2) / 9, 3/
438 DATA (cf(i, 2),i= 1, 2) / 3, 9/
443 matrix_ssx_uux_h=0.d0
444 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
445 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
450 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
451 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
452 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
453 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
454 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
455 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
456 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
457 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
458 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
460 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
461 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,5))
462 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,3))
464 CALL vvs1_0(w(1,3),w(1,5),w(1,7),gc_72,amp(2))
472 ztemp = ztemp + cf(j,i)*jamp(j)
474 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
476 matrix_ssx_uux_h=matrix
483 SUBROUTINE ssx_ssx_h(P,H1,H2,ANS)
503 parameter(nexternal=6)
509 real*8 p(0:3,nexternal),ans
514 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
516 real*8 matrix_ssx_ssx_h
518 INTEGER JC(NEXTERNAL)
519 LOGICAL GOODHEL(NCOMB)
521 DATA goodhel/ncomb*.false./
522 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
523 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
524 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
525 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
526 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
527 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
528 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
529 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
530 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
531 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
532 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
533 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
534 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
535 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
536 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
537 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
538 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
539 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
540 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
541 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
542 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
543 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
544 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
545 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
546 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
547 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
548 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
549 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
550 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
551 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
552 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
553 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
554 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
555 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
556 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
557 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
558 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
559 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
560 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
561 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
562 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
563 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
564 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
565 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
566 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
567 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
568 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
569 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
570 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
571 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
572 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
573 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
574 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
575 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
576 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
577 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
578 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
579 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
580 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
581 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
582 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
583 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
584 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
585 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
595 t=matrix_ssx_ssx_h(p ,h1,h2,nhel(1,ihel),jc(1))
602 REAL*8 FUNCTION matrix_ssx_ssx_h(P,H1,H2,NHEL,IC)
621 parameter(nexternal=6)
622 INTEGER NWAVEFUNCS, NCOLOR
623 parameter(nwavefuncs=7, ncolor=2)
627 parameter(imag1=(0d0,1d0))
631 real*8 p(0:3,nexternal)
632 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
638 real*8 denom(ncolor), cf(ncolor,ncolor)
639 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
640 COMPLEX*16 W(18,NWAVEFUNCS)
642 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
651 DATA (cf(i, 1),i= 1, 2) / 9, 3/
654 DATA (cf(i, 2),i= 1, 2) / 3, 9/
659 matrix_ssx_ssx_h=0.d0
660 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
661 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
666 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
667 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
668 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
669 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
670 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
671 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
672 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
673 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
674 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
676 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
677 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
678 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
680 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
688 ztemp = ztemp + cf(j,i)*jamp(j)
690 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
692 matrix_ssx_ssx_h=matrix
699 SUBROUTINE ssx_ccx_h(P,H1,H2,ANS)
719 parameter(nexternal=6)
725 real*8 p(0:3,nexternal),ans
730 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
732 real*8 matrix_ssx_ccx_h
734 INTEGER JC(NEXTERNAL)
735 LOGICAL GOODHEL(NCOMB)
737 DATA goodhel/ncomb*.false./
738 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
739 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
740 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
741 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
742 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
743 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
744 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
745 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
746 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
747 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
748 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
749 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
750 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
751 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
752 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
753 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
754 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
755 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
756 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
757 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
758 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
759 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
760 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
761 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
762 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
763 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
764 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
765 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
766 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
767 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
768 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
769 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
770 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
771 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
772 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
773 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
774 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
775 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
776 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
777 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
778 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
779 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
780 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
781 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
782 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
783 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
784 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
785 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
786 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
787 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
788 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
789 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
790 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
791 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
792 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
793 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
794 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
795 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
796 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
797 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
798 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
799 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
800 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
801 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
811 t=matrix_ssx_ccx_h(p ,h1,h2,nhel(1,ihel),jc(1))
818 REAL*8 FUNCTION matrix_ssx_ccx_h(P,H1,H2,NHEL,IC)
837 parameter(nexternal=6)
838 INTEGER NWAVEFUNCS, NCOLOR
839 parameter(nwavefuncs=7, ncolor=2)
843 parameter(imag1=(0d0,1d0))
847 real*8 p(0:3,nexternal)
848 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
854 real*8 denom(ncolor), cf(ncolor,ncolor)
855 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
856 COMPLEX*16 W(18,NWAVEFUNCS)
858 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
867 DATA (cf(i, 1),i= 1, 2) / 9, 3/
870 DATA (cf(i, 2),i= 1, 2) / 3, 9/
875 matrix_ssx_ccx_h=0.d0
876 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
877 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
882 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
883 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
884 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
885 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
886 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
887 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
888 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
889 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
890 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
892 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
893 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,5))
894 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
896 CALL vvs1_0(w(1,3),w(1,5),w(1,7),gc_72,amp(2))
904 ztemp = ztemp + cf(j,i)*jamp(j)
906 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
909 matrix_ssx_ccx_h=matrix
916 SUBROUTINE ssx_ucx_h(P,H1,H2,ANS)
936 parameter(nexternal=6)
942 real*8 p(0:3,nexternal),ans
947 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
949 real*8 matrix_ssx_ucx_h
951 INTEGER JC(NEXTERNAL)
952 LOGICAL GOODHEL(NCOMB)
954 DATA goodhel/ncomb*.false./
955 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
956 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
957 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
958 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
959 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
960 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
961 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
962 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
963 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
964 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
965 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
966 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
967 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
968 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
969 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
970 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
971 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
972 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
973 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
974 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
975 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
976 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
977 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
978 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
979 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
980 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
981 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
982 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
983 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
984 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
985 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
986 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
987 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
988 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
989 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
990 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
991 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
992 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
993 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
994 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
995 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
996 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
997 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
998 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
999 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1000 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1001 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1002 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1003 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1004 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1005 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1006 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1007 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1008 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1009 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1010 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1011 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1012 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1013 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1014 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1015 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1016 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1017 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1018 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1028 t=matrix_ssx_ucx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1035 REAL*8 FUNCTION matrix_ssx_ucx_h(P,H1,H2,NHEL,IC)
1052 parameter(ngraphs=1)
1054 parameter(nexternal=6)
1055 INTEGER NWAVEFUNCS, NCOLOR
1056 parameter(nwavefuncs=7, ncolor=1)
1060 parameter(imag1=(0d0,1d0))
1064 real*8 p(0:3,nexternal)
1065 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1071 real*8 denom(ncolor), cf(ncolor,ncolor)
1072 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1073 COMPLEX*16 W(18,NWAVEFUNCS)
1074 COMPLEX*16 DUM0,DUM1
1075 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1084 DATA (cf(i, 1),i= 1, 1) / 9/
1089 matrix_ssx_ucx_h=0.d0
1090 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1091 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1096 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1097 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1098 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1099 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1100 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1101 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1102 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1103 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,6))
1104 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
1106 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1113 ztemp = ztemp + cf(j,i)*jamp(j)
1115 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1117 matrix_ssx_ucx_h=matrix
1124 SUBROUTINE ssx_cux_h(P,H1,H2,ANS)
1144 parameter(nexternal=6)
1146 parameter( ncomb=64)
1150 real*8 p(0:3,nexternal),ans
1155 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1157 real*8 matrix_ssx_cux_h
1158 INTEGER IHEL,IDEN, I
1159 INTEGER JC(NEXTERNAL)
1160 LOGICAL GOODHEL(NCOMB)
1162 DATA goodhel/ncomb*.false./
1163 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1164 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1165 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1166 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1167 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1168 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1169 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1170 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1171 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1172 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1173 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1174 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1175 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1176 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1177 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1178 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1179 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1180 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1181 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1182 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1183 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1184 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1185 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1186 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1187 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1188 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1189 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1190 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1191 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1192 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1193 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1194 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1195 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1196 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1197 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1198 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1199 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1200 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1201 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1202 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1203 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1204 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1205 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1206 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1207 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1208 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1209 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1210 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1211 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1212 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1213 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1214 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1215 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1216 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1217 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1218 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1219 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1220 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1221 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1222 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1223 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1224 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1225 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1226 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1236 t=matrix_ssx_cux_h(p ,h1,h2,nhel(1,ihel),jc(1))
1243 REAL*8 FUNCTION matrix_ssx_cux_h(P,H1,H2,NHEL,IC)
1260 parameter(ngraphs=1)
1262 parameter(nexternal=6)
1263 INTEGER NWAVEFUNCS, NCOLOR
1264 parameter(nwavefuncs=7, ncolor=1)
1268 parameter(imag1=(0d0,1d0))
1272 real*8 p(0:3,nexternal)
1273 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1279 real*8 denom(ncolor), cf(ncolor,ncolor)
1280 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1281 COMPLEX*16 W(18,NWAVEFUNCS)
1282 COMPLEX*16 DUM0,DUM1
1283 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1292 DATA (cf(i, 1),i= 1, 1) / 9/
1297 matrix_ssx_cux_h=0.d0
1298 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1299 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1304 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1305 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1306 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1307 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1308 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1309 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1310 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1311 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
1312 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,3))
1314 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1321 ztemp = ztemp + cf(j,i)*jamp(j)
1323 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1325 matrix_ssx_cux_h=matrix
1335 SUBROUTINE ssx_ddx_noh(P,H1,H2,ANS)
1354 parameter(nexternal=6)
1356 parameter( ncomb=64)
1360 real*8 p(0:3,nexternal),ans
1365 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1367 real*8 matrix_ssx_ddx_noh
1368 INTEGER IHEL,IDEN, I
1369 INTEGER JC(NEXTERNAL)
1370 LOGICAL GOODHEL(NCOMB)
1372 DATA goodhel/ncomb*.false./
1373 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1374 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1375 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1376 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1377 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1378 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1379 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1380 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1381 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1382 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1383 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1384 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1385 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1386 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1387 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1388 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1389 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1390 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1391 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1392 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1393 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1394 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1395 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1396 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1397 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1398 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1399 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1400 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1401 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1402 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1403 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1404 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1405 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1406 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1407 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1408 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1409 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1410 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1411 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1412 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1413 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1414 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1415 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1416 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1417 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1418 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1419 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1420 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1421 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1422 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1423 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1424 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1425 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1426 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1427 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1428 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1429 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1430 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1431 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1432 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1433 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1434 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1435 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1436 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1446 t=matrix_ssx_ddx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1453 REAL*8 FUNCTION matrix_ssx_ddx_noh(P,H1,H2,NHEL,IC)
1469 parameter(ngraphs=32)
1471 parameter(nexternal=6)
1472 INTEGER NWAVEFUNCS, NCOLOR
1473 parameter(nwavefuncs=13, ncolor=2)
1477 parameter(imag1=(0d0,1d0))
1481 real*8 p(0:3,nexternal)
1482 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1488 real*8 denom(ncolor), cf(ncolor,ncolor)
1489 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1490 COMPLEX*16 W(18,NWAVEFUNCS)
1491 COMPLEX*16 DUM0,DUM1
1492 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1501 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1504 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1509 matrix_ssx_ddx_noh=0.d0
1510 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1511 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1516 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1517 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1518 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1519 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1520 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1521 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1522 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
1523 CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
1524 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1526 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1527 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1529 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1530 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
1532 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1534 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1535 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
1536 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1538 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1539 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1541 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1543 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1545 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1546 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1547 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
1549 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
1550 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1552 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
1553 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1555 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
1557 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
1558 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1559 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1561 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
1562 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
1564 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_1,amp(14))
1566 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
1568 CALL ffv2_3_0(w(1,6),w(1,3),w(1,7),gc_50,gc_58,amp(16))
1569 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
1571 CALL ffv1_0(w(1,4),w(1,6),w(1,12),gc_1,amp(17))
1572 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1574 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
1576 CALL ffv2_3_0(w(1,4),w(1,6),w(1,7),gc_50,gc_58,amp(19))
1578 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
1579 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
1581 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(21))
1582 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
1584 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_1,amp(22))
1586 CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(23))
1587 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
1589 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_1,amp(24))
1590 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
1591 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1593 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(25))
1595 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(26))
1597 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(27))
1599 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(28))
1600 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
1602 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(29))
1604 CALL ffv2_3_0(w(1,6),w(1,2),w(1,11),gc_50,gc_58,amp(30))
1606 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(31))
1608 CALL ffv2_3_0(w(1,9),w(1,2),w(1,11),gc_50,gc_58,amp(32))
1609 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1610 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1611 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1612 $ -amp(21)-amp(22)-amp(23)-amp(24)+1./6.*amp(25)+1./6.*amp(26)
1613 $ +1./6.*amp(27)+1./6.*amp(28)-amp(29)-amp(30)-amp(31)-amp(32)
1614 jamp(2)=+1./2.*(-amp(13)-amp(14)-amp(15)-amp(16)-amp(25)-amp(26)
1621 ztemp = ztemp + cf(j,i)*jamp(j)
1623 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1625 matrix_ssx_ddx_noh=matrix
1632 SUBROUTINE ssx_uux_noh(P,H1,H2,ANS)
1650 parameter(nexternal=6)
1652 parameter( ncomb=64)
1656 real*8 p(0:3,nexternal),ans
1661 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1663 real*8 matrix_ssx_uux_noh
1664 INTEGER IHEL,IDEN, I
1665 INTEGER JC(NEXTERNAL)
1666 LOGICAL GOODHEL(NCOMB)
1668 DATA goodhel/ncomb*.false./
1669 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1670 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1671 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1672 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1673 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1674 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1675 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1676 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1677 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1678 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1679 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1680 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1681 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1682 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1683 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1684 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1685 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1686 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1687 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1688 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1689 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1690 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1691 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1692 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1693 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1694 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1695 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1696 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1697 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1698 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1699 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1700 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1701 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1702 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1703 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1704 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1705 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1706 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1707 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1708 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1709 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1710 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1711 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1712 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1713 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1714 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1715 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1716 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1717 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1718 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1719 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1720 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1721 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1722 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1723 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1724 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1725 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1726 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1727 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1728 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1729 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1730 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1731 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1732 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1742 t=matrix_ssx_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1749 REAL*8 FUNCTION matrix_ssx_uux_noh(P,H1,H2,NHEL,IC)
1765 parameter(ngraphs=43)
1767 parameter(nexternal=6)
1768 INTEGER NWAVEFUNCS, NCOLOR
1769 parameter(nwavefuncs=14, ncolor=2)
1773 parameter(imag1=(0d0,1d0))
1777 real*8 p(0:3,nexternal)
1778 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1784 real*8 denom(ncolor), cf(ncolor,ncolor)
1785 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1786 COMPLEX*16 W(18,NWAVEFUNCS)
1787 COMPLEX*16 DUM0,DUM1
1788 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1797 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1800 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1805 matrix_ssx_uux_noh=0.d0
1806 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1807 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1812 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1813 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1814 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1815 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1816 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1817 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1818 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
1819 CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
1820 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1822 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1823 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1825 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1826 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
1828 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1830 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1831 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
1832 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1834 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1835 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1837 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1839 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1841 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1842 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1843 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
1845 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
1846 CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
1848 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
1849 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1851 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
1853 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
1854 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1855 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1857 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
1858 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1860 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
1862 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
1864 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
1865 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
1867 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
1868 CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
1870 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
1872 CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
1874 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
1875 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,9))
1876 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,14))
1877 CALL ffv2_1(w(1,6),w(1,9),gc_108,zero,zero,w(1,10))
1879 CALL ffv2_0(w(1,5),w(1,10),w(1,14),gc_108,amp(21))
1881 CALL vvv1_0(w(1,12),w(1,14),w(1,9),gc_4,amp(22))
1883 CALL vvv1_0(w(1,14),w(1,9),w(1,7),gc_53,amp(23))
1884 CALL ffv2_1(w(1,2),w(1,9),gc_101,zero,zero,w(1,10))
1886 CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_2,amp(24))
1887 CALL ffv2_2(w(1,4),w(1,9),gc_101,zero,zero,w(1,5))
1889 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(25))
1891 CALL ffv2_5_0(w(1,4),w(1,10),w(1,7),gc_51,gc_58,amp(26))
1893 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(27))
1894 CALL ffv2_2(w(1,1),w(1,14),gc_101,zero,zero,w(1,5))
1896 CALL ffv1_0(w(1,5),w(1,3),w(1,12),gc_2,amp(28))
1897 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,10))
1899 CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_101,amp(29))
1901 CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(30))
1902 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,5))
1904 CALL ffv2_0(w(1,5),w(1,3),w(1,14),gc_101,amp(31))
1905 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,14))
1907 CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_1,amp(32))
1909 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(33))
1911 CALL ffv2_3_0(w(1,14),w(1,2),w(1,7),gc_50,gc_58,amp(34))
1913 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(35))
1914 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
1915 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1917 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(36))
1919 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
1921 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(38))
1923 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(39))
1924 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
1926 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(40))
1928 CALL ffv2_3_0(w(1,10),w(1,2),w(1,11),gc_50,gc_58,amp(41))
1930 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(42))
1932 CALL ffv2_3_0(w(1,5),w(1,2),w(1,11),gc_50,gc_58,amp(43))
1933 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1934 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1935 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1936 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
1937 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
1938 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1939 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1940 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
1941 $ -1./2.*amp(38)-1./2.*amp(39)
1947 ztemp = ztemp + cf(j,i)*jamp(j)
1949 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1951 matrix_ssx_uux_noh=matrix
1958 SUBROUTINE ssx_ssx_noh(P,H1,H2,ANS)
1977 parameter(nexternal=6)
1979 parameter( ncomb=64)
1983 real*8 p(0:3,nexternal),ans
1988 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1990 real*8 matrix_ssx_ssx_noh
1991 INTEGER IHEL,IDEN, I
1992 INTEGER JC(NEXTERNAL)
1993 LOGICAL GOODHEL(NCOMB)
1995 DATA goodhel/ncomb*.false./
1996 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1997 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1998 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1999 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2000 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2001 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2002 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2003 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2004 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2005 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2006 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2007 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2008 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2009 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2010 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2011 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2012 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2013 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2014 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2015 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2016 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2017 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2018 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2019 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2020 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2021 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2022 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2023 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2024 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2025 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2026 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2027 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2028 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2029 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2030 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2031 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2032 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2033 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2034 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2035 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2036 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2037 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2038 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2039 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2040 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2041 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2042 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2043 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2044 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2045 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2046 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2047 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2048 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2049 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2050 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2051 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2052 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2053 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2054 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2055 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2056 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2057 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2058 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2059 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2069 t=matrix_ssx_ssx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2076 REAL*8 FUNCTION matrix_ssx_ssx_noh(P,H1,H2,NHEL,IC)
2092 parameter(ngraphs=64)
2094 parameter(nexternal=6)
2095 INTEGER NWAVEFUNCS, NCOLOR
2096 parameter(nwavefuncs=16, ncolor=2)
2100 parameter(imag1=(0d0,1d0))
2104 real*8 p(0:3,nexternal)
2105 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2111 real*8 denom(ncolor), cf(ncolor,ncolor)
2112 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2113 COMPLEX*16 W(18,NWAVEFUNCS)
2114 COMPLEX*16 DUM0,DUM1
2115 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2124 DATA (cf(i, 1),i= 1, 2) / 9, 3/
2127 DATA (cf(i, 2),i= 1, 2) / 3, 9/
2132 matrix_ssx_ssx_noh=0.d0
2133 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
2134 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
2139 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2140 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2141 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2142 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2143 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2144 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2145 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
2146 CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
2147 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2149 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2150 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2152 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2153 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
2155 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2157 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2158 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
2159 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2161 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2162 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2164 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2166 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2168 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2169 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2170 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
2172 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
2173 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
2175 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
2176 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2178 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
2180 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
2181 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2182 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2184 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
2185 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
2187 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(14))
2189 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
2191 CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(16))
2192 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
2194 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_1,amp(17))
2195 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
2197 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
2199 CALL ffv2_3_0(w(1,4),w(1,14),w(1,7),gc_50,gc_58,amp(19))
2201 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
2202 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,9))
2203 CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,14))
2204 CALL ffv1_2(w(1,5),w(1,9),gc_3,mta,zero,w(1,10))
2206 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(21))
2207 CALL ffv1_1(w(1,6),w(1,9),gc_3,mta,zero,w(1,13))
2209 CALL ffv1_0(w(1,5),w(1,13),w(1,14),gc_3,amp(22))
2210 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,15))
2212 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(23))
2214 CALL ffv2_4_0(w(1,5),w(1,13),w(1,15),gc_50,gc_59,amp(24))
2215 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,13))
2216 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,10))
2218 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(25))
2219 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,16))
2221 CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
2223 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(27))
2225 CALL ffv2_4_0(w(1,5),w(1,16),w(1,15),gc_50,gc_59,amp(28))
2226 CALL ffv1_1(w(1,2),w(1,9),gc_1,zero,zero,w(1,16))
2228 CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_1,amp(29))
2229 CALL ffv1_2(w(1,4),w(1,9),gc_1,zero,zero,w(1,5))
2231 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(30))
2233 CALL ffv2_3_0(w(1,4),w(1,16),w(1,7),gc_50,gc_58,amp(31))
2235 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(32))
2236 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
2237 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,16))
2239 CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_1,amp(33))
2240 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,9))
2242 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(34))
2244 CALL ffv2_3_0(w(1,4),w(1,16),w(1,7),gc_50,gc_58,amp(35))
2246 CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(36))
2247 CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,9))
2249 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(37))
2250 CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,16))
2252 CALL ffv1_0(w(1,16),w(1,2),w(1,12),gc_1,amp(38))
2254 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(39))
2256 CALL ffv2_3_0(w(1,16),w(1,2),w(1,7),gc_50,gc_58,amp(40))
2257 CALL ffv1_2(w(1,1),w(1,14),gc_1,zero,zero,w(1,16))
2259 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_1,amp(41))
2260 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,9))
2262 CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_1,amp(42))
2264 CALL ffv2_3_0(w(1,16),w(1,3),w(1,7),gc_50,gc_58,amp(43))
2265 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,16))
2267 CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_1,amp(44))
2268 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,14))
2269 CALL ffv1_2(w(1,1),w(1,14),gc_11,zero,zero,w(1,13))
2271 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(45))
2273 CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_11,amp(46))
2275 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(47))
2277 CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_11,amp(48))
2278 CALL ffv2_3_2(w(1,1),w(1,15),gc_50,gc_58,zero,zero,w(1,14))
2280 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(49))
2282 CALL ffv2_3_0(w(1,9),w(1,3),w(1,15),gc_50,gc_58,amp(50))
2284 CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(51))
2286 CALL ffv2_3_0(w(1,16),w(1,3),w(1,15),gc_50,gc_58,amp(52))
2287 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,15))
2289 CALL ffv1_0(w(1,15),w(1,2),w(1,12),gc_1,amp(53))
2291 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_1,amp(54))
2293 CALL ffv2_3_0(w(1,15),w(1,2),w(1,7),gc_50,gc_58,amp(55))
2295 CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_1,amp(56))
2296 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2297 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2299 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(57))
2301 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(58))
2303 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(59))
2305 CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_11,amp(60))
2306 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2308 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(61))
2310 CALL ffv2_3_0(w(1,9),w(1,2),w(1,11),gc_50,gc_58,amp(62))
2312 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(63))
2314 CALL ffv2_3_0(w(1,16),w(1,2),w(1,11),gc_50,gc_58,amp(64))
2315 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2316 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2317 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2318 $ +1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)+1./2.*amp(36)
2319 $ +1./2.*amp(45)+1./2.*amp(46)+1./2.*amp(47)+1./2.*amp(48)
2320 $ -amp(53)-amp(54)-amp(55)-amp(56)+1./6.*amp(57)+1./6.*amp(58)
2321 $ +1./6.*amp(59)+1./6.*amp(60)-amp(61)-amp(62)-amp(63)-amp(64)
2322 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2323 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2324 $ +amp(28)+amp(29)+amp(30)+amp(31)+amp(32)-1./6.*amp(33)
2325 $ -1./6.*amp(34)-1./6.*amp(35)-1./6.*amp(36)+amp(37)+amp(38)
2326 $ +amp(39)+amp(40)+amp(41)+amp(42)+amp(43)+amp(44)-1./6.*amp(45)
2327 $ -1./6.*amp(46)-1./6.*amp(47)-1./6.*amp(48)+amp(49)+amp(50)
2328 $ +amp(51)+amp(52)-1./2.*amp(57)-1./2.*amp(58)-1./2.*amp(59)
2335 ztemp = ztemp + cf(j,i)*jamp(j)
2337 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2339 matrix_ssx_ssx_noh=matrix
2346 SUBROUTINE ssx_ccx_noh(P,H1,H2,ANS)
2365 parameter(nexternal=6)
2367 parameter( ncomb=64)
2371 real*8 p(0:3,nexternal),ans
2376 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2378 real*8 matrix_ssx_ccx_noh
2379 INTEGER IHEL,IDEN, I
2380 INTEGER JC(NEXTERNAL)
2381 LOGICAL GOODHEL(NCOMB)
2383 DATA goodhel/ncomb*.false./
2384 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2385 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2386 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2387 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2388 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2389 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2390 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2391 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2392 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2393 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2394 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2395 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2396 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2397 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2398 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2399 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2400 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2401 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2402 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2403 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2404 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2405 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2406 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2407 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2408 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2409 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2410 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2411 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2412 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2413 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2414 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2415 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2416 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2417 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2418 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2419 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2420 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2421 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2422 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2423 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2424 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2425 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2426 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2427 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2428 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2429 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2430 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2431 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2432 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2433 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2434 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2435 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2436 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2437 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2438 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2439 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2440 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2441 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2442 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2443 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2444 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2445 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2446 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2447 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2457 t=matrix_ssx_ccx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2464 REAL*8 FUNCTION matrix_ssx_ccx_noh(P,H1,H2,NHEL,IC)
2480 parameter(ngraphs=43)
2482 parameter(nexternal=6)
2483 INTEGER NWAVEFUNCS, NCOLOR
2484 parameter(nwavefuncs=14, ncolor=2)
2488 parameter(imag1=(0d0,1d0))
2492 real*8 p(0:3,nexternal)
2493 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2499 real*8 denom(ncolor), cf(ncolor,ncolor)
2500 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2501 COMPLEX*16 W(18,NWAVEFUNCS)
2502 COMPLEX*16 DUM0,DUM1
2503 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2512 DATA (cf(i, 1),i= 1, 2) / 9, 3/
2515 DATA (cf(i, 2),i= 1, 2) / 3, 9/
2520 matrix_ssx_ccx_noh=0.d0
2521 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
2522 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
2527 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2528 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2529 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2530 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2531 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2532 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2533 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
2534 CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
2535 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2537 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2538 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2540 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2541 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
2543 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2545 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2546 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
2547 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2549 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2550 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2552 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2554 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2556 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2557 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2558 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
2560 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
2561 CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
2563 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
2564 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2566 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
2568 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
2569 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2570 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2572 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
2573 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
2575 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
2577 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
2579 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
2580 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
2582 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
2583 CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
2585 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
2587 CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
2589 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
2590 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,9))
2591 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,14))
2592 CALL ffv2_1(w(1,6),w(1,9),gc_108,zero,zero,w(1,10))
2594 CALL ffv2_0(w(1,5),w(1,10),w(1,14),gc_108,amp(21))
2596 CALL vvv1_0(w(1,12),w(1,14),w(1,9),gc_4,amp(22))
2598 CALL vvv1_0(w(1,14),w(1,9),w(1,7),gc_53,amp(23))
2599 CALL ffv2_1(w(1,2),w(1,9),gc_100,zero,zero,w(1,10))
2601 CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_2,amp(24))
2602 CALL ffv2_2(w(1,4),w(1,9),gc_100,zero,zero,w(1,5))
2604 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(25))
2606 CALL ffv2_5_0(w(1,4),w(1,10),w(1,7),gc_51,gc_58,amp(26))
2608 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2609 CALL ffv2_2(w(1,1),w(1,14),gc_100,zero,zero,w(1,5))
2611 CALL ffv1_0(w(1,5),w(1,3),w(1,12),gc_2,amp(28))
2612 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,10))
2614 CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_100,amp(29))
2616 CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(30))
2617 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,5))
2619 CALL ffv2_0(w(1,5),w(1,3),w(1,14),gc_100,amp(31))
2620 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,14))
2622 CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_1,amp(32))
2624 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(33))
2626 CALL ffv2_3_0(w(1,14),w(1,2),w(1,7),gc_50,gc_58,amp(34))
2628 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(35))
2629 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2630 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2632 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(36))
2634 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
2636 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(38))
2638 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(39))
2639 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2641 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(40))
2643 CALL ffv2_3_0(w(1,10),w(1,2),w(1,11),gc_50,gc_58,amp(41))
2645 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(42))
2647 CALL ffv2_3_0(w(1,5),w(1,2),w(1,11),gc_50,gc_58,amp(43))
2648 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2649 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2650 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2651 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
2652 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
2653 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2654 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2655 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
2656 $ -1./2.*amp(38)-1./2.*amp(39)
2662 ztemp = ztemp + cf(j,i)*jamp(j)
2664 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2666 matrix_ssx_ccx_noh=matrix
2673 SUBROUTINE ssx_ucx_noh(P,H1,H2,ANS)
2692 parameter(nexternal=6)
2694 parameter( ncomb=64)
2698 real*8 p(0:3,nexternal),ans
2703 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2705 real*8 matrix_ssx_ucx_noh
2706 INTEGER IHEL,IDEN, I
2707 INTEGER JC(NEXTERNAL)
2708 LOGICAL GOODHEL(NCOMB)
2710 DATA goodhel/ncomb*.false./
2711 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2712 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2713 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2714 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2715 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2716 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2717 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2718 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2719 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2720 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2721 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2722 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2723 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2724 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2725 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2726 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2727 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2728 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2729 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2730 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2731 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2732 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2733 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2734 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2735 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2736 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2737 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2738 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2739 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2740 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2741 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2742 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2743 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2744 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2745 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2746 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2747 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2748 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2749 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2750 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2751 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2752 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2753 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2754 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2755 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2756 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2757 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2758 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2759 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2760 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2761 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2762 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2763 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2764 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2765 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2766 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2767 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2768 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2769 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2770 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2771 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2772 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2773 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2774 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2784 t=matrix_ssx_ucx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2791 REAL*8 FUNCTION matrix_ssx_ucx_noh(P,H1,H2,NHEL,IC)
2807 parameter(ngraphs=11)
2809 parameter(nexternal=6)
2810 INTEGER NWAVEFUNCS, NCOLOR
2811 parameter(nwavefuncs=10, ncolor=1)
2815 parameter(imag1=(0d0,1d0))
2819 real*8 p(0:3,nexternal)
2820 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2826 real*8 denom(ncolor), cf(ncolor,ncolor)
2827 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2828 COMPLEX*16 W(18,NWAVEFUNCS)
2829 COMPLEX*16 DUM0,DUM1
2830 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2839 DATA (cf(i, 1),i= 1, 1) / 9/
2844 matrix_ssx_ucx_noh=0.d0
2845 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
2846 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
2851 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2852 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2853 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2854 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2855 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2856 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2857 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,7))
2858 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
2859 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2861 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2862 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2864 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2865 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2867 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2868 CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
2870 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2871 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
2873 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2875 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2877 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2878 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
2880 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2881 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2883 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_101,amp(9))
2885 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2886 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2888 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_101,amp(11))
2889 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2890 $ +amp(9)+amp(10)+amp(11)
2896 ztemp = ztemp + cf(j,i)*jamp(j)
2898 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2900 matrix_ssx_ucx_noh=matrix
2907 SUBROUTINE ssx_cux_noh(P,H1,H2,ANS)
2926 parameter(nexternal=6)
2928 parameter( ncomb=64)
2932 real*8 p(0:3,nexternal),ans
2937 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2939 real*8 matrix_ssx_cux_noh
2940 INTEGER IHEL,IDEN, I
2941 INTEGER JC(NEXTERNAL)
2942 LOGICAL GOODHEL(NCOMB)
2944 DATA goodhel/ncomb*.false./
2945 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2946 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2947 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2948 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2949 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2950 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2951 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2952 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2953 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2954 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2955 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2956 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2957 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2958 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2959 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2960 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2961 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2962 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2963 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2964 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2965 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2966 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2967 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2968 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2969 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2970 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2971 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2972 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2973 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2974 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2975 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2976 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2977 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2978 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2979 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2980 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2981 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2982 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2983 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2984 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2985 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2986 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2987 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2988 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2989 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2990 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2991 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2992 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2993 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2994 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2995 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2996 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2997 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2998 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2999 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3000 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3001 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3002 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3003 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3004 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3005 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3006 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3007 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3008 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3018 t=matrix_ssx_cux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3025 REAL*8 FUNCTION matrix_ssx_cux_noh(P,H1,H2,NHEL,IC)
3041 parameter(ngraphs=11)
3043 parameter(nexternal=6)
3044 INTEGER NWAVEFUNCS, NCOLOR
3045 parameter(nwavefuncs=10, ncolor=1)
3049 parameter(imag1=(0d0,1d0))
3053 real*8 p(0:3,nexternal)
3054 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
3060 real*8 denom(ncolor), cf(ncolor,ncolor)
3061 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
3062 COMPLEX*16 W(18,NWAVEFUNCS)
3063 COMPLEX*16 DUM0,DUM1
3064 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3073 DATA (cf(i, 1),i= 1, 1) / 9/
3078 matrix_ssx_cux_noh=0.d0
3079 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
3080 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
3085 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3086 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3087 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3088 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
3089 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3090 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3091 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
3092 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,8))
3093 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
3095 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
3096 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
3098 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
3099 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
3101 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
3102 CALL ffv2_1(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
3104 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
3105 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
3107 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
3109 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
3111 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
3112 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
3114 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
3115 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
3117 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
3119 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
3120 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
3122 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
3123 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
3124 $ +amp(9)+amp(10)+amp(11)
3130 ztemp = ztemp + cf(j,i)*jamp(j)
3132 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3134 matrix_ssx_cux_noh=matrix
3140 SUBROUTINE ssx_gg_noh(P,H1,H2,ANS)
3159 parameter(nexternal=6)
3161 parameter( ncomb=64)
3165 real*8 p(0:3,nexternal),ans
3170 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
3172 real*8 matrix_ssx_gg_noh
3173 INTEGER IHEL,IDEN, I
3174 INTEGER JC(NEXTERNAL)
3175 LOGICAL GOODHEL(NCOMB)
3177 DATA goodhel/ncomb*.false./
3178 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
3179 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
3180 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
3181 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
3182 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
3183 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
3184 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
3185 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
3186 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
3187 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
3188 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
3189 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
3190 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
3191 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
3192 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3193 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3194 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3195 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3196 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3197 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3198 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3199 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3200 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3201 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3202 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3203 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3204 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3205 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3206 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3207 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3208 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3209 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3210 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3211 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3212 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3213 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3214 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3215 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3216 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3217 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3218 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3219 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3220 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3221 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3222 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3223 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3224 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3225 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3226 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3227 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3228 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3229 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3230 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3231 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3232 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3233 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3234 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3235 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3236 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3237 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3238 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3239 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3240 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3241 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3251 t=matrix_ssx_gg_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3258 REAL*8 FUNCTION matrix_ssx_gg_noh(P,H1,H2,NHEL,IC)
3274 parameter(ngraphs=16)
3276 parameter(nexternal=6)
3277 INTEGER NWAVEFUNCS, NCOLOR
3278 parameter(nwavefuncs=10, ncolor=2)
3282 parameter(imag1=(0d0,1d0))
3286 real*8 p(0:3,nexternal)
3287 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
3293 real*8 denom(ncolor), cf(ncolor,ncolor)
3294 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
3295 COMPLEX*16 W(18,NWAVEFUNCS)
3296 COMPLEX*16 DUM0,DUM1
3297 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3306 DATA (cf(i, 1),i= 1, 2) / 16, -2/
3309 DATA (cf(i, 2),i= 1, 2) / -2, 16/
3314 matrix_ssx_gg_noh=0.d0
3315 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
3316 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
3321 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3322 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3323 CALL vxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3324 CALL vxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
3325 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3326 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3327 CALL ffv1_2(w(1,1),w(1,3),gc_11,zero,zero,w(1,7))
3328 CALL ffv1_1(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
3329 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
3331 CALL ffv1_0(w(1,7),w(1,8),w(1,9),gc_1,amp(1))
3332 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
3334 CALL ffv2_3_0(w(1,7),w(1,8),w(1,10),gc_50,gc_58,amp(2))
3335 CALL ffv1_2(w(1,7),w(1,4),gc_11,zero,zero,w(1,6))
3337 CALL ffv1_0(w(1,6),w(1,2),w(1,9),gc_1,amp(3))
3339 CALL ffv2_3_0(w(1,6),w(1,2),w(1,10),gc_50,gc_58,amp(4))
3340 CALL ffv1_2(w(1,1),w(1,4),gc_11,zero,zero,w(1,6))
3341 CALL ffv1_1(w(1,2),w(1,3),gc_11,zero,zero,w(1,7))
3343 CALL ffv1_0(w(1,6),w(1,7),w(1,9),gc_1,amp(5))
3345 CALL ffv2_3_0(w(1,6),w(1,7),w(1,10),gc_50,gc_58,amp(6))
3346 CALL ffv1_2(w(1,6),w(1,3),gc_11,zero,zero,w(1,5))
3348 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(7))
3350 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(8))
3351 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,5))
3353 CALL ffv1_0(w(1,5),w(1,7),w(1,4),gc_11,amp(9))
3354 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
3356 CALL ffv1_0(w(1,6),w(1,7),w(1,4),gc_11,amp(10))
3358 CALL ffv1_0(w(1,5),w(1,8),w(1,3),gc_11,amp(11))
3360 CALL ffv1_0(w(1,6),w(1,8),w(1,3),gc_11,amp(12))
3361 CALL vvv1p0_1(w(1,3),w(1,4),gc_10,zero,zero,w(1,8))
3362 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
3364 CALL ffv1_0(w(1,4),w(1,2),w(1,9),gc_1,amp(13))
3366 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(14))
3368 CALL ffv2_3_0(w(1,4),w(1,2),w(1,10),gc_50,gc_58,amp(15))
3370 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(16))
3371 jamp(1)=+amp(5)+amp(6)+amp(7)+amp(8)+amp(9)+amp(10)-imag1
3372 $ *amp(13)-imag1*amp(14)-imag1*amp(15)-imag1*amp(16)
3373 jamp(2)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(11)+amp(12)+imag1
3374 $ *amp(13)+imag1*amp(14)+imag1*amp(15)+imag1*amp(16)
3380 ztemp = ztemp + cf(j,i)*jamp(j)
3382 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3385 matrix_ssx_gg_noh=matrix