C++InterfacetoTauola
DSX.f
1 
2 c --------------------- subprocesses initiated by d sbar: DSX 21.10.2015
3 
4  SUBROUTINE dsx(P,I3,I4,H1,H2,KEY,ANS)
5  IMPLICIT NONE
6 
7  INTEGER i3,i4, h1,h2, key
8  REAL*8 p(0:3,6), ans, ans1,ans2,ans3,ans4,ans5
9 
10 C I3=0,2,3 OR 4 I4=0,-2,-3 OR -4
11 
12  ans=0.d0
13 
14  IF(key.EQ.1) THEN
15  IF(i3.EQ.2 .AND. i4.EQ.-2) CALL dsx_uux_h(p,h1,h2,ans)
16  IF(i3.EQ.2 .AND. i4.EQ.-4) CALL dsx_ucx_h(p,h1,h2,ans)
17  IF(i3.EQ.1 .AND. i4.EQ.-3) CAll dsx_dsx_h(p,h1,h2,ans)
18  IF(i3.EQ.4 .AND. i4.EQ.-2) CALL dsx_cux_h(p,h1,h2,ans)
19  IF(i3.EQ.4 .AND. i4.EQ.-4) CALL dsx_ccx_h(p,h1,h2,ans)
20  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
21  CALL dsx_uux_h(p,h1,h2,ans1)
22  CALL dsx_ucx_h(p,h1,h2,ans2)
23  CALL dsx_dsx_h(p,h1,h2,ans3)
24  CALL dsx_cux_h(p,h1,h2,ans4)
25  CALL dsx_ccx_h(p,h1,h2,ans5)
26  ans=ans1+ans2+ans3+ans4+ans5
27  ENDIF
28  ELSE IF(key.EQ.0) THEN
29  IF(i3.EQ.2 .AND. i4.EQ.-2) CALL dsx_uux_noh(p,h1,h2,ans)
30  IF(i3.EQ.2 .AND. i4.EQ.-4) CALL dsx_ucx_noh(p,h1,h2,ans)
31  IF(i3.EQ.1 .AND. i4.EQ.-3) CAll dsx_dsx_noh(p,h1,h2,ans)
32  IF(i3.EQ.4 .AND. i4.EQ.-2) CALL dsx_cux_noh(p,h1,h2,ans)
33  IF(i3.EQ.4 .AND. i4.EQ.-4) CALL dsx_ccx_noh(p,h1,h2,ans)
34  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
35  CALL dsx_uux_noh(p,h1,h2,ans1)
36  CALL dsx_ucx_noh(p,h1,h2,ans2)
37  CALL dsx_dsx_noh(p,h1,h2,ans3)
38  CALL dsx_cux_noh(p,h1,h2,ans4)
39  CALL dsx_ccx_noh(p,h1,h2,ans5)
40  ans=ans1+ans2+ans3+ans4+ans5
41 
42  ENDIF
43  ELSE
44  WRITE(*,*) 'NOT FINISHED'
45  stop
46  ENDIF
47  END ! SUBROUTINE DSX
48 
49 
50 
51 C ----- begin subprocesses DSX->jjH with Higgs->tautau
52 
53  SUBROUTINE dsx_uux_h(P,H1,H2,ANS)
54 C
55 C
56 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
57 C By the MadGraph Development Team
58 C Please visit us at https://launchpad.net/madgraph5
59 C
60 C MadGraph StandAlone Version
61 C
62 C Returns amplitude squared summed/avg over colors
63 C and helicities
64 C for the point in phase space P(0:3,NEXTERNAL)
65 C
66 C Process: d s~ > u u~ h WEIGHTED=6
67 C * Decay: h > ta+ ta- WEIGHTED=2
68 C
69  IMPLICIT NONE
70 C
71 C CONSTANTS
72 C
73  INTEGER nexternal
74  parameter(nexternal=6)
75  INTEGER ncomb
76  parameter( ncomb=64)
77 C
78 C ARGUMENTS
79 C
80  REAL*8 p(0:3,nexternal),ans
81  INTEGER h1,h2
82 C
83 C LOCAL VARIABLES
84 C
85  INTEGER nhel(nexternal,ncomb),ntry
86  REAL*8 t
87  REAL*8 matrix_dsx_uux_h
88  INTEGER ihel,iden, i
89  INTEGER jc(nexternal)
90  LOGICAL goodhel(ncomb)
91  DATA ntry/0/
92  DATA goodhel/ncomb*.false./
93  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
94  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
95  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
96  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
97  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
98  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
99  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
100  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
101  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
102  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
103  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
104  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
105  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
106  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
107  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
108  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
109  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
110  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
111  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
112  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
113  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
114  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
115  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
116  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
117  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
118  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
119  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
120  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
121  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
122  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
123  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
124  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
125  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
126  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
127  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
128  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
129  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
130  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
131  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
132  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
133  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
134  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
135  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
136  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
137  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
138  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
139  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
140  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
141  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
142  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
143  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
144  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
145  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
146  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
147  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
148  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
149  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
150  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
151  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
152  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
153  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
154  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
155  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
156  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
157  DATA iden/36/
158 C ----------
159 C BEGIN CODE
160 C ----------
161  DO ihel=1,nexternal
162  jc(ihel) = +1
163  ENDDO
164  ans = 0d0
165  DO ihel=1,ncomb
166  t=matrix_dsx_uux_h(p ,h1,h2,nhel(1,ihel),jc(1))
167  ans=ans+t
168  ENDDO
169  ans=ans/dble(iden)
170  END
171 
172 
173  REAL*8 FUNCTION matrix_dsx_uux_h(P,H1,H2,NHEL,IC)
174 C
175 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
176 C By the MadGraph Development Team
177 C Please visit us at https://launchpad.net/madgraph5
178 C
179 C Returns amplitude squared summed/avg over colors
180 C for the point with external lines W(0:6,NEXTERNAL)
181 C
182 C Process: d s~ > u u~ h WEIGHTED=6
183 C * Decay: h > ta+ ta- WEIGHTED=2
184 C
185  IMPLICIT NONE
186 C
187 C CONSTANTS
188 C
189  INTEGER ngraphs
190  parameter(ngraphs=1)
191  INTEGER nexternal
192  parameter(nexternal=6)
193  INTEGER nwavefuncs, ncolor
194  parameter(nwavefuncs=7, ncolor=1)
195  REAL*8 zero
196  parameter(zero=0d0)
197  COMPLEX*16 imag1
198  parameter(imag1=(0d0,1d0))
199 C
200 C ARGUMENTS
201 C
202  REAL*8 p(0:3,nexternal)
203  INTEGER nhel(nexternal), ic(nexternal)
204 C
205 C LOCAL VARIABLES
206 C
207  INTEGER i,j
208  COMPLEX*16 ztemp
209  REAL*8 denom(ncolor), cf(ncolor,ncolor)
210  COMPLEX*16 amp(ngraphs), jamp(ncolor)
211  COMPLEX*16 w(18,nwavefuncs)
212  COMPLEX*16 dum0,dum1
213  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
214 C
215 C GLOBAL VARIABLES
216 C
217  include 'coupl.inc'
218 C
219 C COLOR DATA
220 C
221  DATA denom(1)/1/
222  DATA (cf(i, 1),i= 1, 1) / 9/
223 C 1 T(2,4) T(3,1)
224 
225  INTEGER h1,h2
226  REAL*8 matrix
227  matrix_dsx_uux_h=0.d0
228  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
229  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
230 
231 C ----------
232 C BEGIN CODE
233 C ----------
234  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
235  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
236  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
237  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
238  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
239  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
240  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
241  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
242  CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,3))
243 C Amplitude(s) for diagram number 1
244  CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
245  jamp(1)=+amp(1)
246 
247  matrix = 0.d0
248  DO i = 1, ncolor
249  ztemp = (0.d0,0.d0)
250  DO j = 1, ncolor
251  ztemp = ztemp + cf(j,i)*jamp(j)
252  ENDDO
253  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
254  ENDDO
255  matrix_dsx_uux_h=matrix
256  ENDIF ! CLOSES H1 IF
257  ENDIF ! CLOSES H2 IF
258  END
259 
260 c==========================================
261 
262  SUBROUTINE dsx_ucx_h(P,H1,H2,ANS)
263 C
264 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
265 C By the MadGraph Development Team
266 C Please visit us at https://launchpad.net/madgraph5
267 C
268 C MadGraph StandAlone Version
269 C
270 C Returns amplitude squared summed/avg over colors
271 C and helicities
272 C for the point in phase space P(0:3,NEXTERNAL)
273 C
274 C Process: d s~ > u c~ h WEIGHTED=6
275 C * Decay: h > ta+ ta- WEIGHTED=2
276 C
277  IMPLICIT NONE
278 C
279 C CONSTANTS
280 C
281  INTEGER nexternal
282  parameter(nexternal=6)
283  INTEGER ncomb
284  parameter( ncomb=64)
285 C
286 C ARGUMENTS
287 C
288  REAL*8 p(0:3,nexternal),ans
289  INTEGER h1,h2
290 C
291 C LOCAL VARIABLES
292 C
293  INTEGER nhel(nexternal,ncomb),ntry
294  REAL*8 t
295  REAL*8 matrix_dsx_ucx_h
296  INTEGER ihel,iden, i
297  INTEGER jc(nexternal)
298  LOGICAL goodhel(ncomb)
299  DATA ntry/0/
300  DATA goodhel/ncomb*.false./
301  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
302  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
303  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
304  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
305  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
306  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
307  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
308  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
309  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
310  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
311  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
312  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
313  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
314  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
315  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
316  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
317  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
318  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
319  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
320  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
321  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
322  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
323  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
324  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
325  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
326  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
327  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
328  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
329  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
330  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
331  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
332  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
333  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
334  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
335  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
336  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
337  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
338  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
339  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
340  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
341  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
342  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
343  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
344  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
345  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
346  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
347  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
348  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
349  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
350  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
351  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
352  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
353  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
354  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
355  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
356  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
357  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
358  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
359  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
360  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
361  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
362  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
363  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
364  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
365  DATA iden/36/
366 C ----------
367 C BEGIN CODE
368 C ----------
369  DO ihel=1,nexternal
370  jc(ihel) = +1
371  ENDDO
372  ans = 0d0
373  DO ihel=1,ncomb
374  t=matrix_dsx_ucx_h(p ,h1,h2,nhel(1,ihel),jc(1))
375  ans=ans+t
376  ENDDO
377  ans=ans/dble(iden)
378  END
379 
380 
381  REAL*8 FUNCTION matrix_dsx_ucx_h(P,H1,H2,NHEL,IC)
382 C
383 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
384 C By the MadGraph Development Team
385 C Please visit us at https://launchpad.net/madgraph5
386 C
387 C Returns amplitude squared summed/avg over colors
388 C for the point with external lines W(0:6,NEXTERNAL)
389 C
390 C Process: d s~ > u c~ h WEIGHTED=6
391 C * Decay: h > ta+ ta- WEIGHTED=2
392 C
393  IMPLICIT NONE
394 C
395 C CONSTANTS
396 C
397  INTEGER ngraphs
398  parameter(ngraphs=1)
399  INTEGER nexternal
400  parameter(nexternal=6)
401  INTEGER nwavefuncs, ncolor
402  parameter(nwavefuncs=7, ncolor=1)
403  REAL*8 zero
404  parameter(zero=0d0)
405  COMPLEX*16 imag1
406  parameter(imag1=(0d0,1d0))
407 C
408 C ARGUMENTS
409 C
410  REAL*8 p(0:3,nexternal)
411  INTEGER nhel(nexternal), ic(nexternal)
412 C
413 C LOCAL VARIABLES
414 C
415  INTEGER i,j
416  COMPLEX*16 ztemp
417  REAL*8 denom(ncolor), cf(ncolor,ncolor)
418  COMPLEX*16 amp(ngraphs), jamp(ncolor)
419  COMPLEX*16 w(18,nwavefuncs)
420  COMPLEX*16 dum0,dum1
421  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
422 C
423 C GLOBAL VARIABLES
424 C
425  include 'coupl.inc'
426 C
427 C COLOR DATA
428 C
429  DATA denom(1)/1/
430  DATA (cf(i, 1),i= 1, 1) / 9/
431 C 1 T(2,4) T(3,1)
432 
433  INTEGER h1,h2
434  REAL*8 matrix
435  matrix_dsx_ucx_h=0.d0
436  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
437  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
438 
439 C ----------
440 C BEGIN CODE
441 C ----------
442  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
443  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
444  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
445  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
446  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
447  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
448  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
449  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
450  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
451 C Amplitude(s) for diagram number 1
452  CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
453  jamp(1)=+amp(1)
454 
455  matrix = 0.d0
456  DO i = 1, ncolor
457  ztemp = (0.d0,0.d0)
458  DO j = 1, ncolor
459  ztemp = ztemp + cf(j,i)*jamp(j)
460  ENDDO
461  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
462  ENDDO
463  matrix_dsx_ucx_h=matrix
464  ENDIF
465  ENDIF
466  END
467 
468 
469 
470 C===============================================
471 
472  SUBROUTINE dsx_dsx_h(P,H1,H2,ANS)
473 C
474 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
475 C By the MadGraph Development Team
476 C Please visit us at https://launchpad.net/madgraph5
477 C
478 C MadGraph StandAlone Version
479 C
480 C Returns amplitude squared summed/avg over colors
481 C and helicities
482 C for the point in phase space P(0:3,NEXTERNAL)
483 C
484 C Process: d s~ > d s~ h WEIGHTED=6
485 C * Decay: h > ta+ ta- WEIGHTED=2
486 C
487  IMPLICIT NONE
488 C
489 C CONSTANTS
490 C
491  INTEGER nexternal
492  parameter(nexternal=6)
493  INTEGER ncomb
494  parameter( ncomb=64)
495 C
496 C ARGUMENTS
497 C
498  REAL*8 p(0:3,nexternal),ans
499  INTEGER h1,h2
500 C
501 C LOCAL VARIABLES
502 C
503  INTEGER nhel(nexternal,ncomb),ntry
504  REAL*8 t
505  REAL*8 matrix_dsx_dsx_h
506  INTEGER ihel,iden, i
507  INTEGER jc(nexternal)
508  LOGICAL goodhel(ncomb)
509  DATA ntry/0/
510  DATA goodhel/ncomb*.false./
511  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
512  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
513  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
514  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
515  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
516  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
517  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
518  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
519  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
520  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
521  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
522  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
523  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
524  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
525  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
526  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
527  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
528  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
529  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
530  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
531  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
532  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
533  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
534  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
535  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
536  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
537  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
538  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
539  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
540  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
541  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
542  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
543  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
544  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
545  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
546  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
547  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
548  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
549  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
550  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
551  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
552  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
553  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
554  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
555  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
556  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
557  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
558  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
559  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
560  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
561  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
562  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
563  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
564  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
565  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
566  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
567  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
568  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
569  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
570  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
571  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
572  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
573  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
574  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
575  DATA iden/36/
576 C ----------
577 C BEGIN CODE
578 C ----------
579  DO ihel=1,nexternal
580  jc(ihel) = +1
581  ENDDO
582  ans = 0d0
583  DO ihel=1,ncomb
584  t=matrix_dsx_dsx_h(p ,h1,h2,nhel(1,ihel),jc(1))
585  ans=ans+t
586  ENDDO
587  ans=ans/dble(iden)
588  END
589 
590 
591  REAL*8 FUNCTION matrix_dsx_dsx_h(P,H1,H2,NHEL,IC)
592 C
593 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
594 C By the MadGraph Development Team
595 C Please visit us at https://launchpad.net/madgraph5
596 C
597 C Returns amplitude squared summed/avg over colors
598 C for the point with external lines W(0:6,NEXTERNAL)
599 C
600 C Process: d s~ > d s~ h WEIGHTED=6
601 C * Decay: h > ta+ ta- WEIGHTED=2
602 C
603  IMPLICIT NONE
604 C
605 C CONSTANTS
606 C
607  INTEGER ngraphs
608  parameter(ngraphs=1)
609  INTEGER nexternal
610  parameter(nexternal=6)
611  INTEGER nwavefuncs, ncolor
612  parameter(nwavefuncs=7, ncolor=1)
613  REAL*8 zero
614  parameter(zero=0d0)
615  COMPLEX*16 imag1
616  parameter(imag1=(0d0,1d0))
617 C
618 C ARGUMENTS
619 C
620  REAL*8 p(0:3,nexternal)
621  INTEGER nhel(nexternal), ic(nexternal)
622 C
623 C LOCAL VARIABLES
624 C
625  INTEGER i,j
626  COMPLEX*16 ztemp
627  REAL*8 denom(ncolor), cf(ncolor,ncolor)
628  COMPLEX*16 amp(ngraphs), jamp(ncolor)
629  COMPLEX*16 w(18,nwavefuncs)
630  COMPLEX*16 dum0,dum1
631  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
632 C
633 C GLOBAL VARIABLES
634 C
635  include 'coupl.inc'
636 C
637 C COLOR DATA
638 C
639  DATA denom(1)/1/
640  DATA (cf(i, 1),i= 1, 1) / 9/
641 C 1 T(2,4) T(3,1)
642 
643  INTEGER h1,h2
644  REAL*8 matrix
645  matrix_dsx_dsx_h=0.d0
646  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
647  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
648 
649 C ----------
650 C BEGIN CODE
651 C ----------
652  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
653  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
654  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
655  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
656  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
657  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
658  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
659  CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,6))
660  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
661 C Amplitude(s) for diagram number 1
662  CALL vvs1_0(w(1,6),w(1,3),w(1,7),gc_81,amp(1))
663  jamp(1)=+amp(1)
664 
665  matrix = 0.d0
666  DO i = 1, ncolor
667  ztemp = (0.d0,0.d0)
668  DO j = 1, ncolor
669  ztemp = ztemp + cf(j,i)*jamp(j)
670  ENDDO
671  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
672  ENDDO
673 
674  matrix_dsx_dsx_h=matrix
675  ENDIF
676  ENDIF
677  END
678 
679 
680 
681 C========================================
682 
683  SUBROUTINE dsx_cux_h(P,H1,H2,ANS)
684 C
685 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
686 C By the MadGraph Development Team
687 C Please visit us at https://launchpad.net/madgraph5
688 C
689 C MadGraph StandAlone Version
690 C
691 C Returns amplitude squared summed/avg over colors
692 C and helicities
693 C for the point in phase space P(0:3,NEXTERNAL)
694 C
695 C Process: d s~ > c u~ h WEIGHTED=6
696 C * Decay: h > ta+ ta- WEIGHTED=2
697 C
698  IMPLICIT NONE
699 C
700 C CONSTANTS
701 C
702  INTEGER nexternal
703  parameter(nexternal=6)
704  INTEGER ncomb
705  parameter( ncomb=64)
706 C
707 C ARGUMENTS
708 C
709  REAL*8 p(0:3,nexternal),ans
710  INTEGER h1,h2
711 C
712 C LOCAL VARIABLES
713 C
714  INTEGER nhel(nexternal,ncomb),ntry
715  REAL*8 t
716  REAL*8 matrix_dsx_cux_h
717  INTEGER ihel,iden, i
718  INTEGER jc(nexternal)
719  LOGICAL goodhel(ncomb)
720  DATA ntry/0/
721  DATA goodhel/ncomb*.false./
722  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
723  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
724  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
725  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
726  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
727  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
728  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
729  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
730  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
731  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
732  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
733  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
734  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
735  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
736  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
737  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
738  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
739  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
740  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
741  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
742  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
743  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
744  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
745  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
746  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
747  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
748  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
749  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
750  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
751  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
752  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
753  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
754  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
755  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
756  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
757  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
758  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
759  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
760  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
761  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
762  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
763  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
764  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
765  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
766  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
767  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
768  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
769  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
770  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
771  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
772  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
773  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
774  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
775  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
776  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
777  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
778  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
779  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
780  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
781  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
782  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
783  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
784  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
785  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
786  DATA iden/36/
787 C ----------
788 C BEGIN CODE
789 C ----------
790  DO ihel=1,nexternal
791  jc(ihel) = +1
792  ENDDO
793  ans = 0d0
794  DO ihel=1,ncomb
795  t=matrix_dsx_cux_h(p ,h1,h2,nhel(1,ihel),jc(1))
796  ans=ans+t
797  ENDDO
798  ans=ans/dble(iden)
799  END
800 
801 
802  REAL*8 FUNCTION matrix_dsx_cux_h(P,H1,H2,NHEL,IC)
803 C
804 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
805 C By the MadGraph Development Team
806 C Please visit us at https://launchpad.net/madgraph5
807 C
808 C Returns amplitude squared summed/avg over colors
809 C for the point with external lines W(0:6,NEXTERNAL)
810 C
811 C Process: d s~ > c u~ h WEIGHTED=6
812 C * Decay: h > ta+ ta- WEIGHTED=2
813 C
814  IMPLICIT NONE
815 C
816 C CONSTANTS
817 C
818  INTEGER ngraphs
819  parameter(ngraphs=1)
820  INTEGER nexternal
821  parameter(nexternal=6)
822  INTEGER nwavefuncs, ncolor
823  parameter(nwavefuncs=7, ncolor=1)
824  REAL*8 zero
825  parameter(zero=0d0)
826  COMPLEX*16 imag1
827  parameter(imag1=(0d0,1d0))
828 C
829 C ARGUMENTS
830 C
831  REAL*8 p(0:3,nexternal)
832  INTEGER nhel(nexternal), ic(nexternal)
833 C
834 C LOCAL VARIABLES
835 C
836  INTEGER i,j
837  COMPLEX*16 ztemp
838  REAL*8 denom(ncolor), cf(ncolor,ncolor)
839  COMPLEX*16 amp(ngraphs), jamp(ncolor)
840  COMPLEX*16 w(18,nwavefuncs)
841  COMPLEX*16 dum0,dum1
842  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
843 C
844 C GLOBAL VARIABLES
845 C
846  include 'coupl.inc'
847 C
848 C COLOR DATA
849 C
850  DATA denom(1)/1/
851  DATA (cf(i, 1),i= 1, 1) / 9/
852 C 1 T(2,4) T(3,1)
853 
854  INTEGER h1,h2
855  REAL*8 matrix
856  matrix_dsx_cux_h=0.d0
857  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
858  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
859 
860 C ----------
861 C BEGIN CODE
862 C ----------
863  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
864  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
865  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
866  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
867  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
868  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
869  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
870  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,6))
871  CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,3))
872 C Amplitude(s) for diagram number 1
873  CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
874  jamp(1)=+amp(1)
875 
876  matrix = 0.d0
877  DO i = 1, ncolor
878  ztemp = (0.d0,0.d0)
879  DO j = 1, ncolor
880  ztemp = ztemp + cf(j,i)*jamp(j)
881  ENDDO
882  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
883  ENDDO
884 
885  matrix_dsx_cux_h=matrix
886  ENDIF
887  ENDIF
888  END
889 
890 
891  SUBROUTINE dsx_ccx_h(P,H1,H2,ANS)
892 C
893 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
894 C By the MadGraph Development Team
895 C Please visit us at https://launchpad.net/madgraph5
896 C
897 C MadGraph StandAlone Version
898 C
899 C Returns amplitude squared summed/avg over colors
900 C and helicities
901 C for the point in phase space P(0:3,NEXTERNAL)
902 C
903 C Process: d s~ > c c~ h WEIGHTED=6
904 C * Decay: h > ta+ ta- WEIGHTED=2
905 C
906  IMPLICIT NONE
907 C
908 C CONSTANTS
909 C
910  INTEGER nexternal
911  parameter(nexternal=6)
912  INTEGER ncomb
913  parameter( ncomb=64)
914 C
915 C ARGUMENTS
916 C
917  REAL*8 p(0:3,nexternal),ans
918  INTEGER h1,h2
919 C
920 C LOCAL VARIABLES
921 C
922  INTEGER nhel(nexternal,ncomb),ntry
923  REAL*8 t
924  REAL*8 matrix_dsx_ccx_h
925  INTEGER ihel,iden, i
926  INTEGER jc(nexternal)
927  LOGICAL goodhel(ncomb)
928  DATA ntry/0/
929  DATA goodhel/ncomb*.false./
930  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
931  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
932  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
933  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
934  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
935  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
936  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
937  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
938  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
939  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
940  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
941  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
942  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
943  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
944  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
945  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
946  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
947  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
948  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
949  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
950  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
951  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
952  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
953  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
954  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
955  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
956  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
957  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
958  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
959  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
960  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
961  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
962  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
963  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
964  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
965  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
966  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
967  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
968  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
969  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
970  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
971  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
972  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
973  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
974  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
975  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
976  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
977  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
978  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
979  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
980  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
981  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
982  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
983  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
984  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
985  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
986  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
987  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
988  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
989  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
990  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
991  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
992  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
993  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
994  DATA iden/36/
995 C ----------
996 C BEGIN CODE
997 C ----------
998  DO ihel=1,nexternal
999  jc(ihel) = +1
1000  ENDDO
1001  ans = 0d0
1002  DO ihel=1,ncomb
1003  t=matrix_dsx_ccx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1004  ans=ans+t
1005  ENDDO
1006  ans=ans/dble(iden)
1007  END
1008 
1009 
1010  REAL*8 FUNCTION matrix_dsx_ccx_h(P,H1,H2,NHEL,IC)
1011 C
1012 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1013 C By the MadGraph Development Team
1014 C Please visit us at https://launchpad.net/madgraph5
1015 C
1016 C Returns amplitude squared summed/avg over colors
1017 C for the point with external lines W(0:6,NEXTERNAL)
1018 C
1019 C Process: d s~ > c c~ h WEIGHTED=6
1020 C * Decay: h > ta+ ta- WEIGHTED=2
1021 C
1022  IMPLICIT NONE
1023 C
1024 C CONSTANTS
1025 C
1026  INTEGER ngraphs
1027  parameter(ngraphs=1)
1028  INTEGER nexternal
1029  parameter(nexternal=6)
1030  INTEGER nwavefuncs, ncolor
1031  parameter(nwavefuncs=7, ncolor=1)
1032  REAL*8 zero
1033  parameter(zero=0d0)
1034  COMPLEX*16 imag1
1035  parameter(imag1=(0d0,1d0))
1036 C
1037 C ARGUMENTS
1038 C
1039  REAL*8 p(0:3,nexternal)
1040  INTEGER nhel(nexternal), ic(nexternal)
1041 C
1042 C LOCAL VARIABLES
1043 C
1044  INTEGER i,j
1045  COMPLEX*16 ztemp
1046  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1047  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1048  COMPLEX*16 w(18,nwavefuncs)
1049  COMPLEX*16 dum0,dum1
1050  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1051 C
1052 C GLOBAL VARIABLES
1053 C
1054  include 'coupl.inc'
1055 C
1056 C COLOR DATA
1057 C
1058  DATA denom(1)/1/
1059  DATA (cf(i, 1),i= 1, 1) / 9/
1060 C 1 T(2,4) T(3,1)
1061 
1062  INTEGER h1,h2
1063  REAL*8 matrix
1064  matrix_dsx_ccx_h=0.d0
1065  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1066  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1067 
1068 C ----------
1069 C BEGIN CODE
1070 C ----------
1071  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1072  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1073  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1074  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1075  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1076  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1077  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1078  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,6))
1079  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
1080 C Amplitude(s) for diagram number 1
1081  CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1082  jamp(1)=+amp(1)
1083 
1084  matrix = 0.d0
1085  DO i = 1, ncolor
1086  ztemp = (0.d0,0.d0)
1087  DO j = 1, ncolor
1088  ztemp = ztemp + cf(j,i)*jamp(j)
1089  ENDDO
1090  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1091  ENDDO
1092 
1093  matrix_dsx_ccx_h=matrix
1094  ENDIF
1095  ENDIF
1096  END
1097 
1098 C --- end of subprocess DSX->jjH, H->tautau
1099 
1100 C ----- begin subprocesses DSX->jjtautau, no Higgs
1101 
1102  SUBROUTINE dsx_uux_noh(P,H1,H2,ANS)
1103 C
1104 
1105 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1106 C By the MadGraph Development Team
1107 C Please visit us at https://launchpad.net/madgraph5
1108 C
1109 C MadGraph StandAlone Version
1110 C
1111 C Returns amplitude squared summed/avg over colors
1112 C and helicities
1113 C for the point in phase space P(0:3,NEXTERNAL)
1114 C
1115 C Process: d s~ > u u~ ta+ ta- / h Qed=4
1116 C
1117  IMPLICIT NONE
1118 C
1119 C CONSTANTS
1120 C
1121  INTEGER nexternal
1122  parameter(nexternal=6)
1123  INTEGER ncomb
1124  parameter( ncomb=64)
1125 C
1126 C ARGUMENTS
1127 C
1128  REAL*8 p(0:3,nexternal),ans
1129  INTEGER h1,h2
1130 C
1131 C LOCAL VARIABLES
1132 C
1133  INTEGER nhel(nexternal,ncomb),ntry
1134  REAL*8 t
1135  REAL*8 matrix_dsx_uux_noh
1136  INTEGER ihel,iden, i
1137  INTEGER jc(nexternal)
1138  LOGICAL goodhel(ncomb)
1139  DATA ntry/0/
1140  DATA goodhel/ncomb*.false./
1141  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1142  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1143  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1144  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1145  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1146  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1147  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1148  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1149  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1150  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1151  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1152  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1153  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1154  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1155  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1156  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1157  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1158  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1159  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1160  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1161  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1162  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1163  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1164  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1165  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1166  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1167  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1168  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1169  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1170  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1171  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1172  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1173  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1174  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1175  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1176  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1177  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1178  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1179  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1180  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1181  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1182  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1183  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1184  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1185  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1186  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1187  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1188  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1189  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1190  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1191  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1192  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1193  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1194  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1195  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1196  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1197  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1198  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1199  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1200  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1201  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1202  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1203  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1204  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1205  DATA iden/36/
1206 C ----------
1207 C BEGIN CODE
1208 C ----------
1209  DO ihel=1,nexternal
1210  jc(ihel) = +1
1211  ENDDO
1212  ans = 0d0
1213  DO ihel=1,ncomb
1214  t=matrix_dsx_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1215  ans=ans+t
1216  ENDDO
1217  ans=ans/dble(iden)
1218  END
1219 
1220 
1221  REAL*8 FUNCTION matrix_dsx_uux_noh(P,H1,H2,NHEL,IC)
1222 C
1223 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1224 C By the MadGraph Development Team
1225 C Please visit us at https://launchpad.net/madgraph5
1226 C
1227 C Returns amplitude squared summed/avg over colors
1228 C for the point with external lines W(0:6,NEXTERNAL)
1229 C
1230 C Process: d s~ > u u~ ta+ ta- / h Qed=4
1231 C
1232  IMPLICIT NONE
1233 C
1234 C CONSTANTS
1235 C
1236  INTEGER ngraphs
1237  parameter(ngraphs=11)
1238  INTEGER nexternal
1239  parameter(nexternal=6)
1240  INTEGER nwavefuncs, ncolor
1241  parameter(nwavefuncs=10, ncolor=1)
1242  REAL*8 zero
1243  parameter(zero=0d0)
1244  COMPLEX*16 imag1
1245  parameter(imag1=(0d0,1d0))
1246 C
1247 C ARGUMENTS
1248 C
1249  REAL*8 p(0:3,nexternal)
1250  INTEGER nhel(nexternal), ic(nexternal)
1251 C
1252 C LOCAL VARIABLES
1253 C
1254  INTEGER i,j
1255  COMPLEX*16 ztemp
1256  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1257  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1258  COMPLEX*16 w(18,nwavefuncs)
1259  COMPLEX*16 dum0,dum1
1260  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1261 C
1262 C GLOBAL VARIABLES
1263 C
1264  include 'coupl.inc'
1265 C
1266 C COLOR DATA
1267 C
1268  DATA denom(1)/1/
1269  DATA (cf(i, 1),i= 1, 1) / 9/
1270 C 1 T(2,4) T(3,1)
1271 
1272  INTEGER h1,h2
1273  REAL*8 matrix
1274  matrix_dsx_uux_noh=0.d0
1275  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1276  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1277 
1278 C ----------
1279 C BEGIN CODE
1280 C ----------
1281  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1282  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1283  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1284  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1285  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1286  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1287  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
1288  CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,8))
1289  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1290 C Amplitude(s) for diagram number 1
1291  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1292  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1293 C Amplitude(s) for diagram number 2
1294  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1295  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1296 C Amplitude(s) for diagram number 3
1297  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1298  CALL ffv2_1(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
1299 C Amplitude(s) for diagram number 4
1300  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1301  CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
1302 C Amplitude(s) for diagram number 5
1303  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
1304 C Amplitude(s) for diagram number 6
1305  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1306 C Amplitude(s) for diagram number 7
1307  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
1308  CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1309 C Amplitude(s) for diagram number 8
1310  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1311  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1312 C Amplitude(s) for diagram number 9
1313  CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
1314 C Amplitude(s) for diagram number 10
1315  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1316  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1317 C Amplitude(s) for diagram number 11
1318  CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
1319  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1320  $ +amp(9)+amp(10)+amp(11)
1321 
1322  matrix = 0.d0
1323  DO i = 1, ncolor
1324  ztemp = (0.d0,0.d0)
1325  DO j = 1, ncolor
1326  ztemp = ztemp + cf(j,i)*jamp(j)
1327  ENDDO
1328  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1329  ENDDO
1330 
1331  matrix_dsx_uux_noh=matrix
1332  ENDIF ! CLOSES H1 IF
1333  ENDIF ! CLOSES H2 IF
1334  END
1335 
1336 
1337  SUBROUTINE dsx_ucx_noh(P,H1,H2,ANS)
1338 C
1339 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1340 C By the MadGraph Development Team
1341 C Please visit us at https://launchpad.net/madgraph5
1342 C
1343 C MadGraph StandAlone Version
1344 C
1345 C Returns amplitude squared summed/avg over colors
1346 C and helicities
1347 C for the point in phase space P(0:3,NEXTERNAL)
1348 C
1349 C Process: d s~ > u c~ ta+ ta- / h Qed=4
1350 C
1351  IMPLICIT NONE
1352 C
1353 C CONSTANTS
1354 C
1355  INTEGER nexternal
1356  parameter(nexternal=6)
1357  INTEGER ncomb
1358  parameter( ncomb=64)
1359 C
1360 C ARGUMENTS
1361 C
1362  REAL*8 p(0:3,nexternal),ans
1363  INTEGER h1,h2
1364 C
1365 C LOCAL VARIABLES
1366 C
1367  INTEGER nhel(nexternal,ncomb),ntry
1368  REAL*8 t
1369  REAL*8 matrix_dsx_ucx_noh
1370  INTEGER ihel,iden, i
1371  INTEGER jc(nexternal)
1372  LOGICAL goodhel(ncomb)
1373  DATA ntry/0/
1374  DATA goodhel/ncomb*.false./
1375  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1376  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1377  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1378  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1379  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1380  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1381  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1382  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1383  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1384  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1385  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1386  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1387  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1388  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1389  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1390  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1391  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1392  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1393  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1394  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1395  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1396  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1397  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1398  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1399  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1400  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1401  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1402  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1403  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1404  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1405  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1406  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1407  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1408  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1409  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1410  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1411  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1412  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1413  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1414  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1415  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1416  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1417  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1418  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1419  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1420  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1421  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1422  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1423  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1424  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1425  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1426  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1427  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1428  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1429  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1430  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1431  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1432  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1433  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1434  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1435  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1436  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1437  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1438  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1439  DATA iden/36/
1440 C ----------
1441 C BEGIN CODE
1442 C ----------
1443  DO ihel=1,nexternal
1444  jc(ihel) = +1
1445  ENDDO
1446  ans = 0d0
1447  DO ihel=1,ncomb
1448  t=matrix_dsx_ucx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1449  ans=ans+t
1450  ENDDO
1451  ans=ans/dble(iden)
1452  END
1453 
1454 
1455  REAL*8 FUNCTION matrix_dsx_ucx_noh(P,H1,H2,NHEL,IC)
1456 C
1457 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1458 C By the MadGraph Development Team
1459 C Please visit us at https://launchpad.net/madgraph5
1460 C
1461 C Returns amplitude squared summed/avg over colors
1462 C for the point with external lines W(0:6,NEXTERNAL)
1463 C
1464 C Process: d s~ > u c~ ta+ ta- / h Qed=4
1465 C
1466  IMPLICIT NONE
1467 C
1468 C CONSTANTS
1469 C
1470  INTEGER ngraphs
1471  parameter(ngraphs=11)
1472  INTEGER nexternal
1473  parameter(nexternal=6)
1474  INTEGER nwavefuncs, ncolor
1475  parameter(nwavefuncs=10, ncolor=1)
1476  REAL*8 zero
1477  parameter(zero=0d0)
1478  COMPLEX*16 imag1
1479  parameter(imag1=(0d0,1d0))
1480 C
1481 C ARGUMENTS
1482 C
1483  REAL*8 p(0:3,nexternal)
1484  INTEGER nhel(nexternal), ic(nexternal)
1485 C
1486 C LOCAL VARIABLES
1487 C
1488  INTEGER i,j
1489  COMPLEX*16 ztemp
1490  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1491  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1492  COMPLEX*16 w(18,nwavefuncs)
1493  COMPLEX*16 dum0,dum1
1494  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1495 C
1496 C GLOBAL VARIABLES
1497 C
1498  include 'coupl.inc'
1499 C
1500 C COLOR DATA
1501 C
1502  DATA denom(1)/1/
1503  DATA (cf(i, 1),i= 1, 1) / 9/
1504 C 1 T(2,4) T(3,1)
1505 
1506  INTEGER h1,h2
1507  REAL*8 matrix
1508  matrix_dsx_ucx_noh=0.d0
1509  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1510  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1511 
1512 C ----------
1513 C BEGIN CODE
1514 C ----------
1515  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1516  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1517  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1518  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1519  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1520  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1521  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
1522  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
1523  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1524 C Amplitude(s) for diagram number 1
1525  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1526  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1527 C Amplitude(s) for diagram number 2
1528  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1529  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1530 C Amplitude(s) for diagram number 3
1531  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1532  CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1533 C Amplitude(s) for diagram number 4
1534  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1535  CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1536 C Amplitude(s) for diagram number 5
1537  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
1538 C Amplitude(s) for diagram number 6
1539  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1540 C Amplitude(s) for diagram number 7
1541  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
1542  CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1543 C Amplitude(s) for diagram number 8
1544  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1545  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1546 C Amplitude(s) for diagram number 9
1547  CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
1548 C Amplitude(s) for diagram number 10
1549  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1550  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1551 C Amplitude(s) for diagram number 11
1552  CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
1553  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1554  $ +amp(9)+amp(10)+amp(11)
1555 
1556  matrix = 0.d0
1557  DO i = 1, ncolor
1558  ztemp = (0.d0,0.d0)
1559  DO j = 1, ncolor
1560  ztemp = ztemp + cf(j,i)*jamp(j)
1561  ENDDO
1562  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1563  ENDDO
1564  matrix_dsx_ucx_noh=matrix
1565  ENDIF ! CLOSES H1 IF
1566  ENDIF ! CLOSES H2 IF
1567  END
1568 
1569  SUBROUTINE dsx_dsx_noh(P,H1,H2,ANS)
1570 C
1571 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1572 C By the MadGraph Development Team
1573 C Please visit us at https://launchpad.net/madgraph5
1574 C
1575 C MadGraph StandAlone Version
1576 C
1577 C Returns amplitude squared summed/avg over colors
1578 C and helicities
1579 C for the point in phase space P(0:3,NEXTERNAL)
1580 C
1581 C Process: d s~ > d s~ ta+ ta- / h Qed=4
1582 C
1583  IMPLICIT NONE
1584 C
1585 C CONSTANTS
1586 C
1587  INTEGER nexternal
1588  parameter(nexternal=6)
1589  INTEGER ncomb
1590  parameter( ncomb=64)
1591 C
1592 C ARGUMENTS
1593 C
1594  REAL*8 p(0:3,nexternal),ans
1595  INTEGER h1,h2
1596 C
1597 C LOCAL VARIABLES
1598 C
1599  INTEGER nhel(nexternal,ncomb),ntry
1600  REAL*8 t
1601  REAL*8 matrix_dsx_dsx_noh
1602  INTEGER ihel,iden, i
1603  INTEGER jc(nexternal)
1604  LOGICAL goodhel(ncomb)
1605  DATA ntry/0/
1606  DATA goodhel/ncomb*.false./
1607  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1608  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1609  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1610  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1611  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1612  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1613  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1614  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1615  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1616  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1617  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1618  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1619  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1620  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1621  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1622  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1623  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1624  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1625  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1626  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1627  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1628  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1629  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1630  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1631  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1632  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1633  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1634  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1635  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1636  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1637  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1638  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1639  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1640  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1641  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1642  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1643  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1644  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1645  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1646  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1647  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1648  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1649  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1650  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1651  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1652  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1653  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1654  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1655  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1656  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1657  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1658  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1659  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1660  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1661  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1662  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1663  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1664  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1665  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1666  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1667  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1668  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1669  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1670  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1671  DATA iden/36/
1672 C ----------
1673 C BEGIN CODE
1674 C ----------
1675  DO ihel=1,nexternal
1676  jc(ihel) = +1
1677  ENDDO
1678  ans = 0d0
1679  DO ihel=1,ncomb
1680  t=matrix_dsx_dsx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1681  ans=ans+t
1682  ENDDO
1683  ans=ans/dble(iden)
1684  END
1685 
1686 
1687  REAL*8 FUNCTION matrix_dsx_dsx_noh(P,H1,H2,NHEL,IC)
1688 C
1689 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1690 C By the MadGraph Development Team
1691 C Please visit us at https://launchpad.net/madgraph5
1692 C
1693 C Returns amplitude squared summed/avg over colors
1694 C for the point with external lines W(0:6,NEXTERNAL)
1695 C
1696 C Process: d s~ > d s~ ta+ ta- / h Qed=4
1697 C
1698  IMPLICIT NONE
1699 C
1700 C CONSTANTS
1701 C
1702  INTEGER ngraphs
1703  parameter(ngraphs=32)
1704  INTEGER nexternal
1705  parameter(nexternal=6)
1706  INTEGER nwavefuncs, ncolor
1707  parameter(nwavefuncs=13, ncolor=2)
1708  REAL*8 zero
1709  parameter(zero=0d0)
1710  COMPLEX*16 imag1
1711  parameter(imag1=(0d0,1d0))
1712 C
1713 C ARGUMENTS
1714 C
1715  REAL*8 p(0:3,nexternal)
1716  INTEGER nhel(nexternal), ic(nexternal)
1717 C
1718 C LOCAL VARIABLES
1719 C
1720  INTEGER i,j
1721  COMPLEX*16 ztemp
1722  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1723  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1724  COMPLEX*16 w(18,nwavefuncs)
1725  COMPLEX*16 dum0,dum1
1726  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1727 C
1728 C GLOBAL VARIABLES
1729 C
1730  include 'coupl.inc'
1731 C
1732 C COLOR DATA
1733 C
1734  DATA denom(1)/1/
1735  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1736 C 1 T(2,1) T(3,4)
1737  DATA denom(2)/1/
1738  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1739 C 1 T(2,4) T(3,1)
1740 
1741  INTEGER h1,h2
1742  REAL*8 matrix
1743  matrix_dsx_dsx_noh=0.d0
1744  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1745  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1746 
1747 C ----------
1748 C BEGIN CODE
1749 C ----------
1750  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1751  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1752  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1753  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1754  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1755  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1756  CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,7))
1757  CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,8))
1758  CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1759 C Amplitude(s) for diagram number 1
1760  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1761  CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1762 C Amplitude(s) for diagram number 2
1763  CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1764  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,11))
1765 C Amplitude(s) for diagram number 3
1766  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1767 C Amplitude(s) for diagram number 4
1768  CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1769  CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,10))
1770  CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1771 C Amplitude(s) for diagram number 5
1772  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1773  CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1774 C Amplitude(s) for diagram number 6
1775  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1776 C Amplitude(s) for diagram number 7
1777  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1778 C Amplitude(s) for diagram number 8
1779  CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1780  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1781  CALL ffv1_1(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
1782 C Amplitude(s) for diagram number 9
1783  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
1784  CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1785 C Amplitude(s) for diagram number 10
1786  CALL ffv1_0(w(1,13),w(1,2),w(1,12),gc_1,amp(10))
1787  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1788 C Amplitude(s) for diagram number 11
1789  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
1790 C Amplitude(s) for diagram number 12
1791  CALL ffv2_3_0(w(1,13),w(1,2),w(1,7),gc_50,gc_58,amp(12))
1792  CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
1793  CALL ffv1_1(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
1794 C Amplitude(s) for diagram number 13
1795  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
1796  CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
1797 C Amplitude(s) for diagram number 14
1798  CALL ffv1_0(w(1,6),w(1,2),w(1,12),gc_1,amp(14))
1799 C Amplitude(s) for diagram number 15
1800  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
1801 C Amplitude(s) for diagram number 16
1802  CALL ffv2_3_0(w(1,6),w(1,2),w(1,7),gc_50,gc_58,amp(16))
1803  CALL ffv2_3_1(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
1804 C Amplitude(s) for diagram number 17
1805  CALL ffv1_0(w(1,4),w(1,6),w(1,12),gc_1,amp(17))
1806  CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1807 C Amplitude(s) for diagram number 18
1808  CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(18))
1809 C Amplitude(s) for diagram number 19
1810  CALL ffv2_3_0(w(1,4),w(1,6),w(1,7),gc_50,gc_58,amp(19))
1811 C Amplitude(s) for diagram number 20
1812  CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(20))
1813  CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
1814 C Amplitude(s) for diagram number 21
1815  CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(21))
1816  CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
1817 C Amplitude(s) for diagram number 22
1818  CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_1,amp(22))
1819 C Amplitude(s) for diagram number 23
1820  CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(23))
1821  CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
1822 C Amplitude(s) for diagram number 24
1823  CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_1,amp(24))
1824  CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,8))
1825  CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1826 C Amplitude(s) for diagram number 25
1827  CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_1,amp(25))
1828 C Amplitude(s) for diagram number 26
1829  CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(26))
1830 C Amplitude(s) for diagram number 27
1831  CALL ffv2_3_0(w(1,4),w(1,3),w(1,7),gc_50,gc_58,amp(27))
1832 C Amplitude(s) for diagram number 28
1833  CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(28))
1834  CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
1835 C Amplitude(s) for diagram number 29
1836  CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_1,amp(29))
1837 C Amplitude(s) for diagram number 30
1838  CALL ffv2_3_0(w(1,6),w(1,3),w(1,11),gc_50,gc_58,amp(30))
1839 C Amplitude(s) for diagram number 31
1840  CALL ffv2_3_0(w(1,8),w(1,3),w(1,7),gc_50,gc_58,amp(31))
1841 C Amplitude(s) for diagram number 32
1842  CALL ffv2_3_0(w(1,9),w(1,3),w(1,11),gc_50,gc_58,amp(32))
1843  jamp(1)=+1./2.*(+amp(13)+amp(14)+amp(15)+amp(16)+amp(25)+amp(26)
1844  $ +amp(27)+amp(28))
1845  jamp(2)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1846  $ +amp(9)+amp(10)+amp(11)+amp(12)-1./6.*amp(13)-1./6.*amp(14)
1847  $ -1./6.*amp(15)-1./6.*amp(16)+amp(17)+amp(18)+amp(19)+amp(20)
1848  $ +amp(21)+amp(22)+amp(23)+amp(24)-1./6.*amp(25)-1./6.*amp(26)
1849  $ -1./6.*amp(27)-1./6.*amp(28)+amp(29)+amp(30)+amp(31)+amp(32)
1850 
1851  matrix = 0.d0
1852  DO i = 1, ncolor
1853  ztemp = (0.d0,0.d0)
1854  DO j = 1, ncolor
1855  ztemp = ztemp + cf(j,i)*jamp(j)
1856  ENDDO
1857  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1858  ENDDO
1859  matrix_dsx_dsx_noh=matrix
1860  ENDIF ! CLOSES H1 IF
1861  ENDIF ! CLOSES H2 IF
1862  END
1863 
1864 
1865  SUBROUTINE dsx_cux_noh(P,H1,H2,ANS)
1866 C
1867 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1868 C By the MadGraph Development Team
1869 C Please visit us at https://launchpad.net/madgraph5
1870 C
1871 C MadGraph StandAlone Version
1872 C
1873 C Returns amplitude squared summed/avg over colors
1874 C and helicities
1875 C for the point in phase space P(0:3,NEXTERNAL)
1876 C
1877 C Process: d s~ > c u~ ta+ ta- / h Qed=4
1878 C
1879  IMPLICIT NONE
1880 C
1881 C CONSTANTS
1882 C
1883  INTEGER nexternal
1884  parameter(nexternal=6)
1885  INTEGER ncomb
1886  parameter( ncomb=64)
1887 C
1888 C ARGUMENTS
1889 C
1890  REAL*8 p(0:3,nexternal),ans
1891  INTEGER h1,h2
1892 C
1893 C LOCAL VARIABLES
1894 C
1895  INTEGER nhel(nexternal,ncomb),ntry
1896  REAL*8 t
1897  REAL*8 matrix_dsx_cux_noh
1898  INTEGER ihel,iden, i
1899  INTEGER jc(nexternal)
1900  LOGICAL goodhel(ncomb)
1901  DATA ntry/0/
1902  DATA goodhel/ncomb*.false./
1903  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1904  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1905  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1906  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1907  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1908  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1909  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1910  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1911  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1912  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1913  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1914  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1915  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1916  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1917  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1918  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1919  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1920  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1921  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1922  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1923  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1924  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1925  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1926  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1927  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1928  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1929  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1930  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1931  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1932  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1933  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1934  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1935  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1936  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1937  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1938  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1939  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1940  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1941  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1942  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1943  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1944  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1945  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1946  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1947  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1948  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1949  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1950  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1951  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1952  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1953  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1954  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1955  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1956  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1957  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1958  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1959  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1960  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1961  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1962  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1963  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1964  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1965  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1966  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1967  DATA iden/36/
1968 C ----------
1969 C BEGIN CODE
1970 C ----------
1971  DO ihel=1,nexternal
1972  jc(ihel) = +1
1973  ENDDO
1974  ans = 0d0
1975  DO ihel=1,ncomb
1976  t=matrix_dsx_cux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1977  ans=ans+t
1978  ENDDO
1979  ans=ans/dble(iden)
1980  END
1981 
1982 
1983  REAL*8 FUNCTION matrix_dsx_cux_noh(P,H1,H2,NHEL,IC)
1984 C
1985 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1986 C By the MadGraph Development Team
1987 C Please visit us at https://launchpad.net/madgraph5
1988 C
1989 C Returns amplitude squared summed/avg over colors
1990 C for the point with external lines W(0:6,NEXTERNAL)
1991 C
1992 C Process: d s~ > c u~ ta+ ta- / h Qed=4
1993 C
1994  IMPLICIT NONE
1995 C
1996 C CONSTANTS
1997 C
1998  INTEGER ngraphs
1999  parameter(ngraphs=11)
2000  INTEGER nexternal
2001  parameter(nexternal=6)
2002  INTEGER nwavefuncs, ncolor
2003  parameter(nwavefuncs=10, ncolor=1)
2004  REAL*8 zero
2005  parameter(zero=0d0)
2006  COMPLEX*16 imag1
2007  parameter(imag1=(0d0,1d0))
2008 C
2009 C ARGUMENTS
2010 C
2011  REAL*8 p(0:3,nexternal)
2012  INTEGER nhel(nexternal), ic(nexternal)
2013 C
2014 C LOCAL VARIABLES
2015 C
2016  INTEGER i,j
2017  COMPLEX*16 ztemp
2018  REAL*8 denom(ncolor), cf(ncolor,ncolor)
2019  COMPLEX*16 amp(ngraphs), jamp(ncolor)
2020  COMPLEX*16 w(18,nwavefuncs)
2021  COMPLEX*16 dum0,dum1
2022  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2023 C
2024 C GLOBAL VARIABLES
2025 C
2026  include 'coupl.inc'
2027 C
2028 C COLOR DATA
2029 C
2030  DATA denom(1)/1/
2031  DATA (cf(i, 1),i= 1, 1) / 9/
2032 C 1 T(2,4) T(3,1)
2033 
2034  INTEGER h1,h2
2035  REAL*8 matrix
2036  matrix_dsx_cux_noh=0.d0
2037  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2038  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2039 
2040 C ----------
2041 C BEGIN CODE
2042 C ----------
2043  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2044  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2045  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2046  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2047  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2048  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2049  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,7))
2050  CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,8))
2051  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2052 C Amplitude(s) for diagram number 1
2053  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2054  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2055 C Amplitude(s) for diagram number 2
2056  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2057  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2058 C Amplitude(s) for diagram number 3
2059  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2060  CALL ffv2_1(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
2061 C Amplitude(s) for diagram number 4
2062  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2063  CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
2064 C Amplitude(s) for diagram number 5
2065  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2066 C Amplitude(s) for diagram number 6
2067  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2068 C Amplitude(s) for diagram number 7
2069  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2070  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
2071 C Amplitude(s) for diagram number 8
2072  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2073  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2074 C Amplitude(s) for diagram number 9
2075  CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_44,amp(9))
2076 C Amplitude(s) for diagram number 10
2077  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2078  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2079 C Amplitude(s) for diagram number 11
2080  CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_44,amp(11))
2081  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2082  $ +amp(9)+amp(10)+amp(11)
2083 
2084  matrix = 0.d0
2085  DO i = 1, ncolor
2086  ztemp = (0.d0,0.d0)
2087  DO j = 1, ncolor
2088  ztemp = ztemp + cf(j,i)*jamp(j)
2089  ENDDO
2090  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2091  ENDDO
2092 
2093  matrix_dsx_cux_noh=matrix
2094  ENDIF ! CLOSES H1 IF
2095  ENDIF ! CLOSES H2 IF
2096  END
2097 
2098 
2099 
2100  SUBROUTINE dsx_ccx_noh(P,H1,H2,ANS)
2101 C
2102 C
2103 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
2104 C By the MadGraph Development Team
2105 C Please visit us at https://launchpad.net/madgraph5
2106 C
2107 C MadGraph StandAlone Version
2108 C
2109 C Returns amplitude squared summed/avg over colors
2110 C and helicities
2111 C for the point in phase space P(0:3,NEXTERNAL)
2112 C
2113 C Process: d s~ > c c~ ta+ ta- / h Qed=4
2114 C
2115  IMPLICIT NONE
2116 C
2117 C CONSTANTS
2118 C
2119  INTEGER nexternal
2120  parameter(nexternal=6)
2121  INTEGER ncomb
2122  parameter( ncomb=64)
2123 C
2124 C ARGUMENTS
2125 C
2126  REAL*8 p(0:3,nexternal),ans
2127  INTEGER h1,h2
2128 C
2129 C LOCAL VARIABLES
2130 C
2131  INTEGER nhel(nexternal,ncomb),ntry
2132  REAL*8 t
2133  REAL*8 matrix_dsx_ccx_noh
2134  INTEGER ihel,iden, i
2135  INTEGER jc(nexternal)
2136  LOGICAL goodhel(ncomb)
2137  DATA ntry/0/
2138  DATA goodhel/ncomb*.false./
2139  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2140  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2141  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2142  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2143  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2144  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2145  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2146  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2147  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2148  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2149  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2150  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2151  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2152  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2153  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2154  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2155  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2156  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2157  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2158  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2159  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2160  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2161  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2162  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2163  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2164  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2165  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2166  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2167  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2168  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2169  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2170  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2171  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2172  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2173  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2174  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2175  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2176  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2177  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2178  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2179  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2180  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2181  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2182  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2183  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2184  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2185  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2186  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2187  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2188  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2189  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2190  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2191  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2192  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2193  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2194  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2195  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2196  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2197  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2198  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2199  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2200  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2201  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2202  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2203  DATA iden/36/
2204 C ----------
2205 C BEGIN CODE
2206 C ----------
2207  DO ihel=1,nexternal
2208  jc(ihel) = +1
2209  ENDDO
2210  ans = 0d0
2211  DO ihel=1,ncomb
2212  t=matrix_dsx_ccx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2213  ans=ans+t
2214  ENDDO
2215  ans=ans/dble(iden)
2216  END
2217 
2218 
2219  REAL*8 FUNCTION matrix_dsx_ccx_noh(P,H1,H2,NHEL,IC)
2220 C
2221 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
2222 C By the MadGraph Development Team
2223 C Please visit us at https://launchpad.net/madgraph5
2224 C
2225 C Returns amplitude squared summed/avg over colors
2226 C for the point with external lines W(0:6,NEXTERNAL)
2227 C
2228 C Process: d s~ > c c~ ta+ ta- / h Qed=4
2229 C
2230  IMPLICIT NONE
2231 C
2232 C CONSTANTS
2233 C
2234  INTEGER ngraphs
2235  parameter(ngraphs=11)
2236  INTEGER nexternal
2237  parameter(nexternal=6)
2238  INTEGER nwavefuncs, ncolor
2239  parameter(nwavefuncs=10, ncolor=1)
2240  REAL*8 zero
2241  parameter(zero=0d0)
2242  COMPLEX*16 imag1
2243  parameter(imag1=(0d0,1d0))
2244 C
2245 C ARGUMENTS
2246 C
2247  REAL*8 p(0:3,nexternal)
2248  INTEGER nhel(nexternal), ic(nexternal)
2249 C
2250 C LOCAL VARIABLES
2251 C
2252  INTEGER i,j
2253  COMPLEX*16 ztemp
2254  REAL*8 denom(ncolor), cf(ncolor,ncolor)
2255  COMPLEX*16 amp(ngraphs), jamp(ncolor)
2256  COMPLEX*16 w(18,nwavefuncs)
2257  COMPLEX*16 dum0,dum1
2258  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2259 C
2260 C GLOBAL VARIABLES
2261 C
2262  include 'coupl.inc'
2263 C
2264 C COLOR DATA
2265 C
2266  DATA denom(1)/1/
2267  DATA (cf(i, 1),i= 1, 1) / 9/
2268 C 1 T(2,4) T(3,1)
2269 
2270  INTEGER h1,h2
2271  REAL*8 matrix
2272  matrix_dsx_ccx_noh=0.d0
2273  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2274  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2275 
2276 C ----------
2277 C BEGIN CODE
2278 C ----------
2279  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2280  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2281  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2282  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2283  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2284  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2285  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,7))
2286  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
2287  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2288 C Amplitude(s) for diagram number 1
2289  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2290  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2291 C Amplitude(s) for diagram number 2
2292  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2293  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2294 C Amplitude(s) for diagram number 3
2295  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2296  CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
2297 C Amplitude(s) for diagram number 4
2298  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2299  CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
2300 C Amplitude(s) for diagram number 5
2301  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2302 C Amplitude(s) for diagram number 6
2303  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2304 C Amplitude(s) for diagram number 7
2305  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2306  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
2307 C Amplitude(s) for diagram number 8
2308  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2309  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2310 C Amplitude(s) for diagram number 9
2311  CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_44,amp(9))
2312 C Amplitude(s) for diagram number 10
2313  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2314  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2315 C Amplitude(s) for diagram number 11
2316  CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_44,amp(11))
2317  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2318  $ +amp(9)+amp(10)+amp(11)
2319 
2320  matrix = 0.d0
2321  DO i = 1, ncolor
2322  ztemp = (0.d0,0.d0)
2323  DO j = 1, ncolor
2324  ztemp = ztemp + cf(j,i)*jamp(j)
2325  ENDDO
2326  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2327  ENDDO
2328  matrix_dsx_ccx_noh=matrix
2329  ENDIF ! CLOSES H1 IF
2330  ENDIF ! CLOSES H2 IF
2331  END
2332 
2333 C ------- end subprocesses DSX->jjtautau, no Higgs
2334 
2335 
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7