4 SUBROUTINE cs(P,I3,I4,H1,H2,KEY,ANS)
7 INTEGER I3,I4, H1,H2, KEY
8 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
15 IF(i3.EQ.2 .AND. i4.EQ.1)
CALL cs_ud_h(p,h1,h2,ans)
16 IF(i3.EQ.2 .AND. i4.EQ.3)
CALL cs_us_h(p,h1,h2,ans)
17 IF(i3.EQ.4 .AND. i4.EQ.1)
CAll cs_cd_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.3)
CALL cs_cs_h(p,h1,h2,ans)
19 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
20 CALL cs_ud_h(p,h1,h2,ans1)
21 CALL cs_us_h(p,h1,h2,ans2)
22 CALL cs_cd_h(p,h1,h2,ans3)
23 CALL cs_cs_h(p,h1,h2,ans4)
24 ans=ans1+ans2+ans3+ans4
26 ELSE IF(key.EQ.0)
THEN
27 IF(i3.EQ.2 .AND. i4.EQ.1)
CALL cs_ud_noh(p,h1,h2,ans)
28 IF(i3.EQ.2 .AND. i4.EQ.3)
CALL cs_us_noh(p,h1,h2,ans)
29 IF(i3.EQ.4 .AND. i4.EQ.1)
CAll cs_cd_noh(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.3)
CALL cs_cs_noh(p,h1,h2,ans)
31 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
32 CALL cs_ud_noh(p,h1,h2,ans1)
33 CALL cs_us_noh(p,h1,h2,ans2)
34 CALL cs_cd_noh(p,h1,h2,ans3)
35 CALL cs_cs_noh(p,h1,h2,ans4)
36 ans=ans1+ans2+ans3+ans4
39 WRITE(*,*)
'NOT FINISHED'
48 SUBROUTINE cs_cd_h(P,H1,H2,ANS)
68 parameter(nexternal=6)
74 real*8 p(0:3,nexternal),ans
79 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
84 LOGICAL GOODHEL(NCOMB)
86 DATA goodhel/ncomb*.false./
87 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
88 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
89 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
90 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
91 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
92 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
93 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
94 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
95 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
96 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
97 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
98 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
99 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
100 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
101 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
102 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
103 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
104 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
105 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
106 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
107 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
108 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
109 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
110 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
111 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
112 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
113 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
114 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
115 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
116 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
117 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
118 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
119 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
120 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
121 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
122 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
123 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
124 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
125 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
126 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
127 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
128 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
129 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
130 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
131 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
132 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
133 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
134 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
135 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
136 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
137 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
138 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
139 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
140 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
141 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
142 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
143 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
144 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
145 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
146 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
147 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
148 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
149 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
150 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
160 t=matrix_cs_cd_h(p,h1,h2 ,nhel(1,ihel),jc(1))
167 REAL*8 FUNCTION matrix_cs_cd_h(P,H1,H2,NHEL,IC)
186 parameter(nexternal=6)
187 INTEGER NWAVEFUNCS, NCOLOR
188 parameter(nwavefuncs=7, ncolor=1)
192 parameter(imag1=(0d0,1d0))
196 real*8 p(0:3,nexternal)
197 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
203 real*8 denom(ncolor), cf(ncolor,ncolor)
204 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
205 COMPLEX*16 W(18,NWAVEFUNCS)
207 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
216 DATA (cf(i, 1),i= 1, 1) / 9/
222 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
223 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
228 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
229 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
230 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
231 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
232 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
233 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
234 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
235 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,6))
236 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,4))
238 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
245 ztemp = ztemp + cf(j,i)*jamp(j)
247 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
249 matrix_cs_cd_h=matrix
256 SUBROUTINE cs_cs_h(P,H1,H2,ANS)
276 parameter(nexternal=6)
282 real*8 p(0:3,nexternal),ans
287 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
289 real*8 matrix_cs_cs_h
291 INTEGER JC(NEXTERNAL)
292 LOGICAL GOODHEL(NCOMB)
294 DATA goodhel/ncomb*.false./
295 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
296 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
297 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
298 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
299 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
300 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
301 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
302 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
303 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
304 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
305 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
306 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
307 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
308 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
309 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
310 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
311 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
312 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
313 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
314 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
315 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
316 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
317 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
318 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
319 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
320 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
321 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
322 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
323 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
324 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
325 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
326 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
327 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
328 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
329 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
330 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
331 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
332 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
333 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
334 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
335 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
336 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
337 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
338 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
339 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
340 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
341 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
342 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
343 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
344 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
345 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
346 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
347 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
348 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
349 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
350 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
351 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
352 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
353 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
354 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
355 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
356 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
357 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
358 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
368 t=matrix_cs_cs_h(p,h1,h2 ,nhel(1,ihel),jc(1))
375 REAL*8 FUNCTION matrix_cs_cs_h(P,H1,H2,NHEL,IC)
394 parameter(nexternal=6)
395 INTEGER NWAVEFUNCS, NCOLOR
396 parameter(nwavefuncs=7, ncolor=2)
400 parameter(imag1=(0d0,1d0))
404 real*8 p(0:3,nexternal)
405 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
411 real*8 denom(ncolor), cf(ncolor,ncolor)
412 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
413 COMPLEX*16 W(18,NWAVEFUNCS)
415 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
424 DATA (cf(i, 1),i= 1, 2) / 9, 3/
427 DATA (cf(i, 2),i= 1, 2) / 3, 9/
433 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
434 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
441 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
442 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
443 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
444 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
445 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
446 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
447 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
448 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,6))
449 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
451 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
452 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,5))
453 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,4))
455 CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_72,amp(2))
463 ztemp = ztemp + cf(j,i)*jamp(j)
465 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
467 matrix_cs_cs_h=matrix
476 SUBROUTINE cs_ud_h(P,H1,H2,ANS)
496 parameter(nexternal=6)
502 real*8 p(0:3,nexternal),ans
507 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
509 real*8 matrix_cs_ud_h
511 INTEGER JC(NEXTERNAL)
512 LOGICAL GOODHEL(NCOMB)
514 DATA goodhel/ncomb*.false./
515 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
516 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
517 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
518 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
519 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
520 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
521 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
522 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
523 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
524 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
525 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
526 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
527 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
528 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
529 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
530 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
531 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
532 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
533 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
534 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
535 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
536 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
537 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
538 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
539 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
540 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
541 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
542 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
543 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
544 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
545 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
546 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
547 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
548 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
549 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
550 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
551 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
552 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
553 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
554 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
555 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
556 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
557 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
558 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
559 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
560 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
561 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
562 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
563 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
564 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
565 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
566 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
567 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
568 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
569 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
570 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
571 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
572 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
573 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
574 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
575 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
576 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
577 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
578 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
588 t=matrix_cs_ud_h(p,h1,h2 ,nhel(1,ihel),jc(1))
595 REAL*8 FUNCTION matrix_cs_ud_h(P,H1,H2,NHEL,IC)
614 parameter(nexternal=6)
615 INTEGER NWAVEFUNCS, NCOLOR
616 parameter(nwavefuncs=7, ncolor=1)
620 parameter(imag1=(0d0,1d0))
624 real*8 p(0:3,nexternal)
625 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
631 real*8 denom(ncolor), cf(ncolor,ncolor)
632 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
633 COMPLEX*16 W(18,NWAVEFUNCS)
635 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
644 DATA (cf(i, 1),i= 1, 1) / 9/
650 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
651 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
657 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
658 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
659 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
660 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
661 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
662 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
663 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
664 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,6))
665 CALL ffv2_3(w(1,2),w(1,3),gc_101,mw,ww,w(1,4))
667 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
674 ztemp = ztemp + cf(j,i)*jamp(j)
676 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
678 matrix_cs_ud_h=matrix
687 SUBROUTINE cs_us_h(P,H1,H2,ANS)
707 parameter(nexternal=6)
713 real*8 p(0:3,nexternal),ans
718 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
720 real*8 matrix_cs_us_h
722 INTEGER JC(NEXTERNAL)
723 LOGICAL GOODHEL(NCOMB)
725 DATA goodhel/ncomb*.false./
726 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
727 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
728 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
729 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
730 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
731 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
732 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
733 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
734 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
735 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
736 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
737 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
738 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
739 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
740 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
741 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
742 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
743 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
744 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
745 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
746 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
747 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
748 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
749 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
750 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
751 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
752 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
753 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
754 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
755 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
756 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
757 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
758 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
759 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
760 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
761 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
762 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
763 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
764 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
765 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
766 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
767 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
768 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
769 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
770 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
771 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
772 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
773 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
774 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
775 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
776 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
777 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
778 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
779 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
780 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
781 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
782 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
783 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
784 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
785 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
786 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
787 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
788 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
789 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
799 t=matrix_cs_us_h(p,h1,h2 ,nhel(1,ihel),jc(1))
806 REAL*8 FUNCTION matrix_cs_us_h(P,H1,H2,NHEL,IC)
825 parameter(nexternal=6)
826 INTEGER NWAVEFUNCS, NCOLOR
827 parameter(nwavefuncs=7, ncolor=1)
831 parameter(imag1=(0d0,1d0))
835 real*8 p(0:3,nexternal)
836 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
842 real*8 denom(ncolor), cf(ncolor,ncolor)
843 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
844 COMPLEX*16 W(18,NWAVEFUNCS)
846 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
855 DATA (cf(i, 1),i= 1, 1) / 9/
861 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
862 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
868 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
869 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
870 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
871 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
872 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
873 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
874 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
875 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,6))
876 CALL ffv2_3(w(1,2),w(1,3),gc_101,mw,ww,w(1,4))
878 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
885 ztemp = ztemp + cf(j,i)*jamp(j)
887 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
889 matrix_cs_us_h=matrix
901 SUBROUTINE cs_cd_noh(P,H1,H2,ANS)
920 parameter(nexternal=6)
926 real*8 p(0:3,nexternal),ans
931 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
933 real*8 matrix_cs_cd_noh
935 INTEGER JC(NEXTERNAL)
936 LOGICAL GOODHEL(NCOMB)
938 DATA goodhel/ncomb*.false./
939 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
940 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
941 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
942 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
943 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
944 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
945 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
946 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
947 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
948 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
949 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
950 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
951 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
952 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
953 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
954 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
955 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
956 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
957 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
958 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
959 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
960 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
961 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
962 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
963 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
964 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
965 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
966 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
967 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
968 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
969 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
970 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
971 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
972 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
973 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
974 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
975 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
976 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
977 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
978 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
979 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
980 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
981 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
982 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
983 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
984 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
985 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
986 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
987 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
988 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
989 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
990 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
991 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
992 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
993 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
994 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
995 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
996 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
997 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
998 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
999 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1000 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1001 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1002 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1012 t=matrix_cs_cd_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1019 REAL*8 FUNCTION matrix_cs_cd_noh(P,H1,H2,NHEL,IC)
1035 parameter(ngraphs=11)
1037 parameter(nexternal=6)
1038 INTEGER NWAVEFUNCS, NCOLOR
1039 parameter(nwavefuncs=10, ncolor=1)
1043 parameter(imag1=(0d0,1d0))
1047 real*8 p(0:3,nexternal)
1048 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1054 real*8 denom(ncolor), cf(ncolor,ncolor)
1055 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1056 COMPLEX*16 W(18,NWAVEFUNCS)
1057 COMPLEX*16 DUM0,DUM1
1058 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1067 DATA (cf(i, 1),i= 1, 1) / 9/
1072 matrix_cs_cd_noh=0.d0
1073 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1074 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1079 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1080 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1081 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1082 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1083 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1084 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1085 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,7))
1086 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,8))
1087 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1089 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1090 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1092 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1093 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1095 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1096 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1098 CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1099 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,5))
1101 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1103 CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1105 CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1106 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1108 CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1109 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1111 CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_44,amp(9))
1113 CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1114 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1116 CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_44,amp(11))
1117 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1118 $ +amp(9)+amp(10)+amp(11)
1124 ztemp = ztemp + cf(j,i)*jamp(j)
1126 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1128 matrix_cs_cd_noh=matrix
1134 SUBROUTINE cs_cs_noh(P,H1,H2,ANS)
1153 parameter(nexternal=6)
1155 parameter( ncomb=64)
1159 real*8 p(0:3,nexternal),ans
1164 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1166 real*8 matrix_cs_cs_noh
1167 INTEGER IHEL,IDEN, I
1168 INTEGER JC(NEXTERNAL)
1169 LOGICAL GOODHEL(NCOMB)
1171 DATA goodhel/ncomb*.false./
1172 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1173 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1174 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1175 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1176 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1177 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1178 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1179 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1180 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1181 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1182 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1183 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1184 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1185 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1186 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1187 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1188 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1189 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1190 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1191 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1192 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1193 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1194 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1195 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1196 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1197 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1198 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1199 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1200 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1201 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1202 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1203 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1204 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1205 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1206 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1207 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1208 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1209 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1210 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1211 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1212 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1213 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1214 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1215 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1216 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1217 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1218 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1219 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1220 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1221 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1222 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1223 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1224 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1225 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1226 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1227 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1228 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1229 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1230 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1231 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1232 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1233 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1234 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1235 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1245 t=matrix_cs_cs_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1252 REAL*8 FUNCTION matrix_cs_cs_noh(P,H1,H2,NHEL,IC)
1268 parameter(ngraphs=43)
1270 parameter(nexternal=6)
1271 INTEGER NWAVEFUNCS, NCOLOR
1272 parameter(nwavefuncs=14, ncolor=2)
1276 parameter(imag1=(0d0,1d0))
1280 real*8 p(0:3,nexternal)
1281 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1287 real*8 denom(ncolor), cf(ncolor,ncolor)
1288 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1289 COMPLEX*16 W(18,NWAVEFUNCS)
1290 COMPLEX*16 DUM0,DUM1
1291 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1300 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1303 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1308 matrix_cs_cs_noh=0.d0
1309 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1310 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1316 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1317 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1318 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1319 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1320 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1321 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1322 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,7))
1323 CALL ffv1p0_3(w(1,2),w(1,4),gc_1,zero,zero,w(1,8))
1324 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1326 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1327 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1329 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1330 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,11))
1332 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1334 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1335 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,10))
1336 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1338 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1339 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1341 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1343 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1345 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1346 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1347 CALL ffv1_2(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
1349 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(9))
1350 CALL ffv1_1(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1352 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_1,amp(10))
1353 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1355 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(11))
1357 CALL ffv2_3_0(w(1,2),w(1,13),w(1,7),gc_50,gc_58,amp(12))
1358 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
1359 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
1361 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(13))
1362 CALL ffv1_1(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1364 CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_1,amp(14))
1366 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(15))
1368 CALL ffv2_3_0(w(1,2),w(1,14),w(1,7),gc_50,gc_58,amp(16))
1369 CALL ffv2_3_2(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
1371 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_1,amp(17))
1372 CALL ffv2_3_1(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1374 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
1376 CALL ffv2_3_0(w(1,14),w(1,4),w(1,7),gc_50,gc_58,amp(19))
1378 CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(20))
1379 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,9))
1380 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,14))
1381 CALL ffv2_2(w(1,5),w(1,9),gc_108,zero,zero,w(1,10))
1383 CALL ffv2_0(w(1,10),w(1,6),w(1,14),gc_108,amp(21))
1385 CALL vvv1_0(w(1,12),w(1,9),w(1,14),gc_4,amp(22))
1387 CALL vvv1_0(w(1,9),w(1,14),w(1,7),gc_53,amp(23))
1388 CALL ffv2_2(w(1,2),w(1,9),gc_100,zero,zero,w(1,10))
1390 CALL ffv1_0(w(1,10),w(1,3),w(1,12),gc_2,amp(24))
1391 CALL ffv2_1(w(1,3),w(1,9),gc_100,zero,zero,w(1,6))
1393 CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(25))
1395 CALL ffv2_5_0(w(1,10),w(1,3),w(1,7),gc_51,gc_58,amp(26))
1397 CALL ffv2_3_0(w(1,2),w(1,6),w(1,7),gc_50,gc_58,amp(27))
1398 CALL ffv2_2(w(1,1),w(1,14),gc_100,zero,zero,w(1,6))
1400 CALL ffv1_0(w(1,6),w(1,4),w(1,12),gc_1,amp(28))
1401 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,10))
1403 CALL ffv2_0(w(1,10),w(1,4),w(1,14),gc_100,amp(29))
1405 CALL ffv2_3_0(w(1,6),w(1,4),w(1,7),gc_50,gc_58,amp(30))
1406 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,6))
1408 CALL ffv2_0(w(1,6),w(1,4),w(1,14),gc_100,amp(31))
1409 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,14))
1411 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(32))
1413 CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_2,amp(33))
1415 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(34))
1417 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_2,amp(35))
1418 CALL ffv1p0_3(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
1419 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1421 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(36))
1423 CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_11,amp(37))
1425 CALL ffv2_5_0(w(1,4),w(1,3),w(1,7),gc_51,gc_58,amp(38))
1427 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(39))
1428 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
1430 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(40))
1432 CALL ffv2_5_0(w(1,10),w(1,3),w(1,11),gc_51,gc_58,amp(41))
1434 CALL ffv2_5_0(w(1,8),w(1,3),w(1,7),gc_51,gc_58,amp(42))
1436 CALL ffv2_5_0(w(1,6),w(1,3),w(1,11),gc_51,gc_58,amp(43))
1437 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1438 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1439 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1440 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
1441 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
1442 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1443 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1444 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
1445 $ -1./2.*amp(38)-1./2.*amp(39)
1451 ztemp = ztemp + cf(j,i)*jamp(j)
1453 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1455 matrix_cs_cs_noh=matrix
1460 SUBROUTINE cs_ud_noh(P,H1,H2,ANS)
1479 parameter(nexternal=6)
1481 parameter( ncomb=64)
1485 real*8 p(0:3,nexternal),ans
1490 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1492 real*8 matrix_cs_ud_noh
1493 INTEGER IHEL,IDEN, I
1494 INTEGER JC(NEXTERNAL)
1495 LOGICAL GOODHEL(NCOMB)
1497 DATA goodhel/ncomb*.false./
1498 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1499 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1500 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1501 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1502 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1503 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1504 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1505 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1506 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1507 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1508 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1509 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1510 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1511 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1512 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1513 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1514 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1515 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1516 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1517 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1518 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1519 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1520 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1521 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1522 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1523 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1524 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1525 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1526 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1527 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1528 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1529 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1530 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1531 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1532 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1533 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1534 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1535 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1536 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1537 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1538 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1539 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1540 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1541 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1542 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1543 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1544 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1545 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1546 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1547 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1548 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1549 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1550 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1551 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1552 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1553 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1554 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1555 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1556 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1557 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1558 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1559 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1560 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1561 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1571 t=matrix_cs_ud_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1578 REAL*8 FUNCTION matrix_cs_ud_noh(P,H1,H2,NHEL,IC)
1594 parameter(ngraphs=11)
1596 parameter(nexternal=6)
1597 INTEGER NWAVEFUNCS, NCOLOR
1598 parameter(nwavefuncs=10, ncolor=1)
1602 parameter(imag1=(0d0,1d0))
1606 real*8 p(0:3,nexternal)
1607 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1613 real*8 denom(ncolor), cf(ncolor,ncolor)
1614 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1615 COMPLEX*16 W(18,NWAVEFUNCS)
1616 COMPLEX*16 DUM0,DUM1
1617 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1626 DATA (cf(i, 1),i= 1, 1) / 9/
1631 matrix_cs_ud_noh=0.d0
1632 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1633 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1639 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1640 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1641 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1642 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1643 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1644 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1645 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,7))
1646 CALL ffv2_3(w(1,2),w(1,3),gc_101,mw,ww,w(1,8))
1647 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1649 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1650 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1652 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1653 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1655 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1656 CALL ffv2_2(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
1658 CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1659 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,5))
1661 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1663 CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1665 CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1666 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1668 CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1669 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1671 CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_44,amp(9))
1673 CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1674 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1676 CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_44,amp(11))
1677 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1678 $ +amp(9)+amp(10)+amp(11)
1684 ztemp = ztemp + cf(j,i)*jamp(j)
1686 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1688 matrix_cs_ud_noh=matrix
1694 SUBROUTINE cs_us_noh(P,H1,H2,ANS)
1713 parameter(nexternal=6)
1715 parameter( ncomb=64)
1719 real*8 p(0:3,nexternal),ans
1724 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1726 real*8 matrix_cs_us_noh
1727 INTEGER IHEL,IDEN, I
1728 INTEGER JC(NEXTERNAL)
1729 LOGICAL GOODHEL(NCOMB)
1731 DATA goodhel/ncomb*.false./
1732 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1733 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1734 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1735 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1736 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1737 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1738 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1739 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1740 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1741 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1742 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1743 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1744 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1745 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1746 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1747 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1748 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1749 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1750 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1751 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1752 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1753 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1754 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1755 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1756 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1757 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1758 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1759 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1760 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1761 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1762 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1763 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1764 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1765 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1766 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1767 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1768 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1769 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1770 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1771 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1772 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1773 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1774 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1775 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1776 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1777 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1778 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1779 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1780 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1781 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1782 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1783 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1784 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1785 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1786 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1787 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1788 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1789 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1790 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1791 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1792 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1793 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1794 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1795 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1805 t=matrix_cs_us_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1812 REAL*8 FUNCTION matrix_cs_us_noh(P,H1,H2,NHEL,IC)
1828 parameter(ngraphs=11)
1830 parameter(nexternal=6)
1831 INTEGER NWAVEFUNCS, NCOLOR
1832 parameter(nwavefuncs=10, ncolor=1)
1836 parameter(imag1=(0d0,1d0))
1840 real*8 p(0:3,nexternal)
1841 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1847 real*8 denom(ncolor), cf(ncolor,ncolor)
1848 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1849 COMPLEX*16 W(18,NWAVEFUNCS)
1850 COMPLEX*16 DUM0,DUM1
1851 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1860 DATA (cf(i, 1),i= 1, 1) / 9/
1865 matrix_cs_us_noh=0.d0
1866 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1867 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1873 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1874 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1875 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1876 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1877 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1878 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1879 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,7))
1880 CALL ffv2_3(w(1,2),w(1,3),gc_101,mw,ww,w(1,8))
1881 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1883 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1884 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1886 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1887 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1889 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1890 CALL ffv2_2(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
1892 CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1893 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,5))
1895 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1897 CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1899 CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1900 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1902 CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1903 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1905 CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_100,amp(9))
1907 CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1908 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1910 CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_100,amp(11))
1911 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1912 $ +amp(9)+amp(10)+amp(11)
1918 ztemp = ztemp + cf(j,i)*jamp(j)
1920 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1922 matrix_cs_us_noh=matrix