C++InterfacetoTauola
DDX.f
1 C --------- begin processes initiated by DDbar
2 
3  SUBROUTINE ddx(P,I3,I4,H1,H2,KEY,ANS)
4  IMPLICIT NONE
5 
6  INTEGER i3,i4, h1,h2, key
7  REAL*8 p(0:3,6), ans, ans1,ans2,ans3,ans4, ans5,ans6
8 
9  REAL*8 pp(0:3,6)
10 
11 C I3=0,1,2,3 OR 4 I4=0,-1,-2,-3, OR -4
12 C or I3=21 AND I4=21 -> GLUONS AS FINAL JETS
13 
14  ans=0.d0
15 
16  ! switch 4-vectors for I3, I4 and provide/(add contribution from) new configuration
17  pp(0:3,1) = p(0:3,1)
18  pp(0:3,2) = p(0:3,2)
19  pp(0:3,5) = p(0:3,5)
20  pp(0:3,6) = p(0:3,6)
21 
22  pp(0:3,3) = p(0:3,4)
23  pp(0:3,4) = p(0:3,3)
24 
25  IF(key.EQ.1) THEN
26  IF(i3.EQ.1 .AND. i4.EQ.-1) CALL ddx_ddx_h(p,h1,h2,ans)
27  IF(i3.EQ.2 .AND. i4.EQ.-2) CALL ddx_uux_h(p,h1,h2,ans)
28  IF(i3.EQ.3 .AND. i4.EQ.-3) CAll ddx_ssx_h(p,h1,h2,ans)
29  IF(i3.EQ.4 .AND. i4.EQ.-4) CALL ddx_ccx_h(p,h1,h2,ans)
30  IF(i3.EQ.2 .AND. i4.EQ.-4) CALL ddx_ucx_h(p,h1,h2,ans)
31  IF(i3.EQ.4 .AND. i4.EQ.-2) CALL ddx_cux_h(p,h1,h2,ans)
32  IF(i3.EQ.-1 .AND. i4.EQ.1) CALL ddx_ddx_h(pp,h1,h2,ans)
33  IF(i3.EQ.-2 .AND. i4.EQ.2) CALL ddx_uux_h(pp,h1,h2,ans)
34  IF(i3.EQ.-3 .AND. i4.EQ.3) CAll ddx_ssx_h(pp,h1,h2,ans)
35  IF(i3.EQ.-4 .AND. i4.EQ.4) CALL ddx_ccx_h(pp,h1,h2,ans)
36  IF(i3.EQ.-2 .AND. i4.EQ.4) CALL ddx_cux_h(pp,h1,h2,ans)
37  IF(i3.EQ.-4 .AND. i4.EQ.2) CALL ddx_ucx_h(pp,h1,h2,ans)
38  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
39  CALL ddx_ddx_h(p,h1,h2,ans1)
40  CALL ddx_uux_h(p,h1,h2,ans2)
41  CALL ddx_ssx_h(p,h1,h2,ans3)
42  CALL ddx_ccx_h(p,h1,h2,ans4)
43  CALL ddx_ucx_h(p,h1,h2,ans5)
44  CALL ddx_cux_h(p,h1,h2,ans6)
45  ans=ans1+ans2+ans3+ans4+ans5+ans6
46  CALL ddx_ddx_h(pp,h1,h2,ans1)
47  CALL ddx_uux_h(pp,h1,h2,ans2)
48  CALL ddx_ssx_h(pp,h1,h2,ans3)
49  CALL ddx_ccx_h(pp,h1,h2,ans4)
50  CALL ddx_ucx_h(pp,h1,h2,ans5)
51  CALL ddx_cux_h(pp,h1,h2,ans6)
52  ans=ans1+ans2+ans3+ans4+ans5+ans6+ans
53  ENDIF
54  ELSE IF(key.EQ.0) THEN
55  IF(i3.EQ.1 .AND. i4.EQ.-1) CALL ddx_ddx_noh(p,h1,h2,ans)
56  IF(i3.EQ.2 .AND. i4.EQ.-2) CALL ddx_uux_noh(p,h1,h2,ans)
57  IF(i3.EQ.3 .AND. i4.EQ.-3) CAll ddx_ssx_noh(p,h1,h2,ans)
58  IF(i3.EQ.4 .AND. i4.EQ.-4) CALL ddx_ccx_noh(p,h1,h2,ans)
59  IF(i3.EQ.2 .AND. i4.EQ.-4) CALL ddx_ucx_noh(p,h1,h2,ans)
60  IF(i3.EQ.4 .AND. i4.EQ.-2) CALL ddx_cux_noh(p,h1,h2,ans)
61 
62  IF(i3.EQ.21 .AND. i4.EQ.21) CALL ddx_gg_noh(p,h1,h2,ans)
63  IF(i3.EQ.-1 .AND. i4.EQ.1) CALL ddx_ddx_noh(pp,h1,h2,ans)
64  IF(i3.EQ.-2 .AND. i4.EQ.2) CALL ddx_uux_noh(pp,h1,h2,ans)
65  IF(i3.EQ.-3 .AND. i4.EQ.3) CAll ddx_ssx_noh(pp,h1,h2,ans)
66  IF(i3.EQ.-4 .AND. i4.EQ.4) CALL ddx_ccx_noh(pp,h1,h2,ans)
67  IF(i3.EQ.-2 .AND. i4.EQ.4) CALL ddx_cux_noh(pp,h1,h2,ans)
68  IF(i3.EQ.-4 .AND. i4.EQ.2) CALL ddx_ucx_noh(pp,h1,h2,ans)
69  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
70  CALL ddx_ddx_noh(p,h1,h2,ans1)
71  CALL ddx_uux_noh(p,h1,h2,ans2)
72  CALL ddx_ssx_noh(p,h1,h2,ans3)
73  CALL ddx_ccx_noh(p,h1,h2,ans4)
74  CALL ddx_ucx_noh(p,h1,h2,ans5)
75  CALL ddx_cux_noh(p,h1,h2,ans6)
76  ans=ans1+ans2+ans3+ans4+ans5+ans6
77  CALL ddx_ddx_noh(pp,h1,h2,ans1)
78  CALL ddx_uux_noh(pp,h1,h2,ans2)
79  CALL ddx_ssx_noh(pp,h1,h2,ans3)
80  CALL ddx_ccx_noh(pp,h1,h2,ans4)
81  CALL ddx_ucx_noh(pp,h1,h2,ans5)
82  CALL ddx_cux_noh(pp,h1,h2,ans6)
83  ans=ans1+ans2+ans3+ans4+ans5+ans6+ans
84  CALL ddx_gg_noh(p,h1,h2,ans1)
85  ans=ans1+ans
86  ENDIF
87  ELSE
88  WRITE(*,*) 'NOT FINISHED'
89  stop
90  ENDIF
91  END SUBROUTINE ddx
92 
93 
94 
95 C ----- begin subprocesses DDX->DDX with Higgs->tautau
96 
97  SUBROUTINE ddx_ddx_h(P,H1,H2,ANS)
98 C
99 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
100 C By the MadGraph Development Team
101 C Please visit us at https://launchpad.net/madgraph5
102 C
103 C MadGraph StandAlone Version
104 C
105 C Returns amplitude squared summed/avg over colors
106 C and helicities
107 C for the point in phase space P(0:3,NEXTERNAL)
108 C
109 C Process: d d~ > d d~ h WEIGHTED=6
110 C * Decay: h > ta+ ta- WEIGHTED=2
111 C
112  IMPLICIT NONE
113 C
114 C CONSTANTS
115 C
116  INTEGER nexternal
117  parameter(nexternal=6)
118  INTEGER ncomb
119  parameter( ncomb=64)
120 C
121 C ARGUMENTS
122 C
123  REAL*8 p(0:3,nexternal),ans
124  INTEGER h1,h2
125 C
126 C LOCAL VARIABLES
127 C
128  INTEGER nhel(nexternal,ncomb),ntry
129  REAL*8 t
130  REAL*8 matrix_ddx_ddx_h
131  INTEGER ihel,iden, i
132  INTEGER jc(nexternal)
133  LOGICAL goodhel(ncomb)
134  DATA ntry/0/
135  DATA goodhel/ncomb*.false./
136  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
137  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
138  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
139  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
140  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
141  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
142  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
143  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
144  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
145  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
146  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
147  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
148  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
149  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
150  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
151  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
152  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
153  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
154  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
155  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
156  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
157  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
158  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
159  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
160  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
161  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
162  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
163  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
164  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
165  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
166  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
167  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
168  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
169  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
170  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
171  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
172  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
173  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
174  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
175  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
176  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
177  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
178  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
179  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
180  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
181  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
182  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
183  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
184  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
185  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
186  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
187  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
188  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
189  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
190  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
191  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
192  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
193  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
194  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
195  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
196  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
197  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
198  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
199  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
200  DATA iden/36/
201 C ----------
202 C BEGIN CODE
203 C ----------
204  DO ihel=1,nexternal
205  jc(ihel) = +1
206  ENDDO
207  ans = 0d0
208  DO ihel=1,ncomb
209  t=matrix_ddx_ddx_h(p ,h1,h2,nhel(1,ihel),jc(1))
210  ans=ans+t
211  ENDDO
212  ans=ans/dble(iden)
213  END
214 
215 
216  REAL*8 FUNCTION matrix_ddx_ddx_h(P,H1,H2,NHEL,IC)
217 C
218 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
219 C By the MadGraph Development Team
220 C Please visit us at https://launchpad.net/madgraph5
221 C
222 C Returns amplitude squared summed/avg over colors
223 C for the point with external lines W(0:6,NEXTERNAL)
224 C
225 C Process: d d~ > d d~ h WEIGHTED=6
226 C * Decay: h > ta+ ta- WEIGHTED=2
227 C
228  IMPLICIT NONE
229 C
230 C CONSTANTS
231 C
232  INTEGER ngraphs
233  parameter(ngraphs=2)
234  INTEGER nexternal
235  parameter(nexternal=6)
236  INTEGER nwavefuncs, ncolor
237  parameter(nwavefuncs=7, ncolor=2)
238  REAL*8 zero
239  parameter(zero=0d0)
240  COMPLEX*16 imag1
241  parameter(imag1=(0d0,1d0))
242 C
243 C ARGUMENTS
244 C
245  REAL*8 p(0:3,nexternal)
246  INTEGER nhel(nexternal), ic(nexternal)
247 C
248 C LOCAL VARIABLES
249 C
250  INTEGER i,j
251  COMPLEX*16 ztemp
252  REAL*8 denom(ncolor), cf(ncolor,ncolor)
253  COMPLEX*16 amp(ngraphs), jamp(ncolor)
254  COMPLEX*16 w(18,nwavefuncs)
255  COMPLEX*16 dum0,dum1
256  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
257 C
258 C GLOBAL VARIABLES
259 C
260  include 'coupl.inc'
261 C
262 C COLOR DATA
263 C
264  DATA denom(1)/1/
265  DATA (cf(i, 1),i= 1, 2) / 9, 3/
266 C 1 T(2,1) T(3,4)
267  DATA denom(2)/1/
268  DATA (cf(i, 2),i= 1, 2) / 3, 9/
269 C 1 T(2,4) T(3,1)
270 
271  INTEGER h1,h2
272  REAL*8 matrix
273  matrix_ddx_ddx_h=0.d0
274  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
275  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
276 
277 
278 C ----------
279 C BEGIN CODE
280 C ----------
281  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
282  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
283  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
284  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
285  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
286  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
287  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
288  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
289  CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
290 C Amplitude(s) for diagram number 1
291  CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
292  CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
293  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
294 C Amplitude(s) for diagram number 2
295  CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
296  jamp(1)=-amp(1)
297  jamp(2)=+amp(2)
298 
299  matrix = 0.d0
300  DO i = 1, ncolor
301  ztemp = (0.d0,0.d0)
302  DO j = 1, ncolor
303  ztemp = ztemp + cf(j,i)*jamp(j)
304  ENDDO
305  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
306  ENDDO
307  matrix_ddx_ddx_h=matrix
308  ENDIF
309  ENDIF
310  END
311 
312 C ----- begin subprocesses DDX->UUX with Higgs->tautau
313 
314  SUBROUTINE ddx_uux_h(P,H1,H2,ANS)
315 C
316 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
317 C By the MadGraph Development Team
318 C Please visit us at https://launchpad.net/madgraph5
319 C
320 C MadGraph StandAlone Version
321 C
322 C Returns amplitude squared summed/avg over colors
323 C and helicities
324 C for the point in phase space P(0:3,NEXTERNAL)
325 C
326 C Process: d d~ > u u~ h WEIGHTED=6
327 C * Decay: h > ta+ ta- WEIGHTED=2
328 C
329  IMPLICIT NONE
330 C
331 C CONSTANTS
332 C
333  INTEGER nexternal
334  parameter(nexternal=6)
335  INTEGER ncomb
336  parameter( ncomb=64)
337 C
338 C ARGUMENTS
339 C
340  REAL*8 p(0:3,nexternal),ans
341  INTEGER h1,h2
342 C
343 C LOCAL VARIABLES
344 C
345  INTEGER nhel(nexternal,ncomb),ntry
346  REAL*8 t
347  REAL*8 matrix_ddx_uux_h
348  INTEGER ihel,iden, i
349  INTEGER jc(nexternal)
350  LOGICAL goodhel(ncomb)
351  DATA ntry/0/
352  DATA goodhel/ncomb*.false./
353  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
354  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
355  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
356  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
357  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
358  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
359  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
360  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
361  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
362  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
363  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
364  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
365  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
366  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
367  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
368  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
369  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
370  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
371  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
372  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
373  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
374  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
375  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
376  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
377  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
378  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
379  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
380  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
381  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
382  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
383  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
384  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
385  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
386  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
387  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
388  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
389  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
390  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
391  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
392  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
393  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
394  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
395  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
396  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
397  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
398  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
399  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
400  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
401  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
402  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
403  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
404  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
405  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
406  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
407  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
408  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
409  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
410  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
411  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
412  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
413  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
414  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
415  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
416  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
417  DATA iden/36/
418 C ----------
419 C BEGIN CODE
420 C ----------
421  DO ihel=1,nexternal
422  jc(ihel) = +1
423  ENDDO
424  ans = 0d0
425  DO ihel=1,ncomb
426  t=matrix_ddx_uux_h(p ,h1,h2,nhel(1,ihel),jc(1))
427  ans=ans+t
428  ENDDO
429  ans=ans/dble(iden)
430  END
431 
432 
433  REAL*8 FUNCTION matrix_ddx_uux_h(P,H1,H2,NHEL,IC)
434 C
435 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
436 C By the MadGraph Development Team
437 C Please visit us at https://launchpad.net/madgraph5
438 C
439 C Returns amplitude squared summed/avg over colors
440 C for the point with external lines W(0:6,NEXTERNAL)
441 C
442 C Process: d d~ > u u~ h WEIGHTED=6
443 C * Decay: h > ta+ ta- WEIGHTED=2
444 C
445  IMPLICIT NONE
446 C
447 C CONSTANTS
448 C
449  INTEGER ngraphs
450  parameter(ngraphs=2)
451  INTEGER nexternal
452  parameter(nexternal=6)
453  INTEGER nwavefuncs, ncolor
454  parameter(nwavefuncs=7, ncolor=2)
455  REAL*8 zero
456  parameter(zero=0d0)
457  COMPLEX*16 imag1
458  parameter(imag1=(0d0,1d0))
459 C
460 C ARGUMENTS
461 C
462  REAL*8 p(0:3,nexternal)
463  INTEGER nhel(nexternal), ic(nexternal)
464 C
465 C LOCAL VARIABLES
466 C
467  INTEGER i,j
468  COMPLEX*16 ztemp
469  REAL*8 denom(ncolor), cf(ncolor,ncolor)
470  COMPLEX*16 amp(ngraphs), jamp(ncolor)
471  COMPLEX*16 w(18,nwavefuncs)
472  COMPLEX*16 dum0,dum1
473  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
474 C
475 C GLOBAL VARIABLES
476 C
477  include 'coupl.inc'
478 C
479 C COLOR DATA
480 C
481  DATA denom(1)/1/
482  DATA (cf(i, 1),i= 1, 2) / 9, 3/
483 C 1 T(2,1) T(3,4)
484  DATA denom(2)/1/
485  DATA (cf(i, 2),i= 1, 2) / 3, 9/
486 C 1 T(2,4) T(3,1)
487 
488  INTEGER h1,h2
489  REAL*8 matrix
490  matrix_ddx_uux_h=0.d0
491  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
492  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
493 
494 C ----------
495 C BEGIN CODE
496 C ----------
497  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
498  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
499  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
500  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
501  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
502  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
503  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
504  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
505  CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
506 C Amplitude(s) for diagram number 1
507  CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
508  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,5))
509  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
510 C Amplitude(s) for diagram number 2
511  CALL vvs1_0(w(1,3),w(1,5),w(1,7),gc_72,amp(2))
512  jamp(1)=-amp(1)
513  jamp(2)=+amp(2)
514 
515  matrix = 0.d0
516  DO i = 1, ncolor
517  ztemp = (0.d0,0.d0)
518  DO j = 1, ncolor
519  ztemp = ztemp + cf(j,i)*jamp(j)
520  ENDDO
521  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
522  ENDDO
523  matrix_ddx_uux_h=matrix
524  ENDIF
525  ENDIF
526  END
527 
528 C ----- begin subprocesses DDX->SSX with Higgs->tautau
529 
530  SUBROUTINE ddx_ssx_h(P,H1,H2,ANS)
531 C
532 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
533 C By the MadGraph Development Team
534 C Please visit us at https://launchpad.net/madgraph5
535 C
536 C MadGraph StandAlone Version
537 C
538 C Returns amplitude squared summed/avg over colors
539 C and helicities
540 C for the point in phase space P(0:3,NEXTERNAL)
541 C
542 C Process: d d~ > s s~ h WEIGHTED=6
543 C * Decay: h > ta+ ta- WEIGHTED=2
544 C
545  IMPLICIT NONE
546 C
547 C CONSTANTS
548 C
549  INTEGER nexternal
550  parameter(nexternal=6)
551  INTEGER ncomb
552  parameter( ncomb=64)
553 C
554 C ARGUMENTS
555 C
556  REAL*8 p(0:3,nexternal),ans
557  INTEGER h1,h2
558 C
559 C LOCAL VARIABLES
560 C
561  INTEGER nhel(nexternal,ncomb),ntry
562  REAL*8 t
563  REAL*8 matrix_ddx_ssx_h
564  INTEGER ihel,iden, i
565  INTEGER jc(nexternal)
566  LOGICAL goodhel(ncomb)
567  DATA ntry/0/
568  DATA goodhel/ncomb*.false./
569  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
570  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
571  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
572  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
573  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
574  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
575  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
576  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
577  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
578  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
579  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
580  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
581  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
582  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
583  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
584  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
585  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
586  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
587  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
588  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
589  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
590  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
591  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
592  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
593  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
594  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
595  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
596  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
597  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
598  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
599  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
600  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
601  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
602  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
603  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
604  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
605  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
606  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
607  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
608  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
609  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
610  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
611  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
612  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
613  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
614  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
615  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
616  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
617  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
618  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
619  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
620  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
621  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
622  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
623  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
624  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
625  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
626  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
627  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
628  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
629  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
630  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
631  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
632  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
633  DATA iden/36/
634 C ----------
635 C BEGIN CODE
636 C ----------
637  DO ihel=1,nexternal
638  jc(ihel) = +1
639  ENDDO
640  ans = 0d0
641  DO ihel=1,ncomb
642  t=matrix_ddx_ssx_h(p , h1, h2, nhel(1,ihel),jc(1))
643  ans=ans+t
644  ENDDO
645  ans=ans/dble(iden)
646  END
647 
648 
649  REAL*8 FUNCTION matrix_ddx_ssx_h(P,H1,H2,NHEL,IC)
650 C
651 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
652 C By the MadGraph Development Team
653 C Please visit us at https://launchpad.net/madgraph5
654 C
655 C Returns amplitude squared summed/avg over colors
656 C for the point with external lines W(0:6,NEXTERNAL)
657 C
658 C Process: d d~ > s s~ h WEIGHTED=6
659 C * Decay: h > ta+ ta- WEIGHTED=2
660 C
661  IMPLICIT NONE
662 C
663 C CONSTANTS
664 C
665  INTEGER ngraphs
666  parameter(ngraphs=1)
667  INTEGER nexternal
668  parameter(nexternal=6)
669  INTEGER nwavefuncs, ncolor
670  parameter(nwavefuncs=7, ncolor=1)
671  REAL*8 zero
672  parameter(zero=0d0)
673  COMPLEX*16 imag1
674  parameter(imag1=(0d0,1d0))
675 C
676 C ARGUMENTS
677 C
678  REAL*8 p(0:3,nexternal)
679  INTEGER nhel(nexternal), ic(nexternal)
680 C
681 C LOCAL VARIABLES
682 C
683  INTEGER i,j
684  COMPLEX*16 ztemp
685  REAL*8 denom(ncolor), cf(ncolor,ncolor)
686  COMPLEX*16 amp(ngraphs), jamp(ncolor)
687  COMPLEX*16 w(18,nwavefuncs)
688  COMPLEX*16 dum0,dum1
689  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
690 C
691 C GLOBAL VARIABLES
692 C
693  include 'coupl.inc'
694 C
695 C COLOR DATA
696 C
697  DATA denom(1)/1/
698  DATA (cf(i, 1),i= 1, 1) / 9/
699 C 1 T(2,1) T(3,4)
700 
701  INTEGER h1,h2
702  REAL*8 matrix
703  matrix_ddx_ssx_h=0.d0
704  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
705  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
706 
707 C ----------
708 C BEGIN CODE
709 C ----------
710  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
711  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
712  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
713  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
714  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
715  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
716  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
717  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
718  CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,2))
719 C Amplitude(s) for diagram number 1
720  CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_81,amp(1))
721  jamp(1)=-amp(1)
722 
723  matrix = 0.d0
724  DO i = 1, ncolor
725  ztemp = (0.d0,0.d0)
726  DO j = 1, ncolor
727  ztemp = ztemp + cf(j,i)*jamp(j)
728  ENDDO
729  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
730  ENDDO
731  matrix_ddx_ssx_h=matrix
732  ENDIF
733  ENDIF
734  END
735 
736 C ----- begin subprocesses DDX->CCX with Higgs->tautau
737 
738  SUBROUTINE ddx_ccx_h(P,H1,H2,ANS)
739 C
740 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
741 C By the MadGraph Development Team
742 C Please visit us at https://launchpad.net/madgraph5
743 C
744 C MadGraph StandAlone Version
745 C
746 C Returns amplitude squared summed/avg over colors
747 C and helicities
748 C for the point in phase space P(0:3,NEXTERNAL)
749 C
750 C Process: d d~ > c c~ h WEIGHTED=6
751 C * Decay: h > ta+ ta- WEIGHTED=2
752 C
753  IMPLICIT NONE
754 C
755 C CONSTANTS
756 C
757  INTEGER nexternal
758  parameter(nexternal=6)
759  INTEGER ncomb
760  parameter( ncomb=64)
761 C
762 C ARGUMENTS
763 C
764  REAL*8 p(0:3,nexternal),ans
765  INTEGER h1,h2
766 C
767 C LOCAL VARIABLES
768 C
769  INTEGER nhel(nexternal,ncomb),ntry
770  REAL*8 t
771  REAL*8 matrix_ddx_ccx_h
772  INTEGER ihel,iden, i
773  INTEGER jc(nexternal)
774  LOGICAL goodhel(ncomb)
775  DATA ntry/0/
776  DATA goodhel/ncomb*.false./
777  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
778  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
779  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
780  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
781  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
782  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
783  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
784  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
785  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
786  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
787  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
788  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
789  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
790  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
791  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
792  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
793  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
794  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
795  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
796  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
797  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
798  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
799  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
800  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
801  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
802  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
803  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
804  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
805  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
806  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
807  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
808  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
809  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
810  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
811  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
812  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
813  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
814  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
815  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
816  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
817  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
818  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
819  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
820  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
821  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
822  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
823  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
824  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
825  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
826  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
827  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
828  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
829  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
830  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
831  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
832  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
833  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
834  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
835  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
836  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
837  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
838  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
839  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
840  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
841  DATA iden/36/
842 C ----------
843 C BEGIN CODE
844 C ----------
845  DO ihel=1,nexternal
846  jc(ihel) = +1
847  ENDDO
848  ans = 0d0
849  DO ihel=1,ncomb
850  t=matrix_ddx_ccx_h(p ,h1,h2,nhel(1,ihel),jc(1))
851  ans=ans+t
852  ENDDO
853  ans=ans/dble(iden)
854  END
855 
856 
857  REAL*8 FUNCTION matrix_ddx_ccx_h(P,H1,H2,NHEL,IC)
858 C
859 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
860 C By the MadGraph Development Team
861 C Please visit us at https://launchpad.net/madgraph5
862 C
863 C Returns amplitude squared summed/avg over colors
864 C for the point with external lines W(0:6,NEXTERNAL)
865 C
866 C Process: d d~ > c c~ h WEIGHTED=6
867 C * Decay: h > ta+ ta- WEIGHTED=2
868 C
869  IMPLICIT NONE
870 C
871 C CONSTANTS
872 C
873  INTEGER ngraphs
874  parameter(ngraphs=2)
875  INTEGER nexternal
876  parameter(nexternal=6)
877  INTEGER nwavefuncs, ncolor
878  parameter(nwavefuncs=7, ncolor=2)
879  REAL*8 zero
880  parameter(zero=0d0)
881  COMPLEX*16 imag1
882  parameter(imag1=(0d0,1d0))
883 C
884 C ARGUMENTS
885 C
886  REAL*8 p(0:3,nexternal)
887  INTEGER nhel(nexternal), ic(nexternal)
888 C
889 C LOCAL VARIABLES
890 C
891  INTEGER i,j
892  COMPLEX*16 ztemp
893  REAL*8 denom(ncolor), cf(ncolor,ncolor)
894  COMPLEX*16 amp(ngraphs), jamp(ncolor)
895  COMPLEX*16 w(18,nwavefuncs)
896  COMPLEX*16 dum0,dum1
897  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
898 C
899 C GLOBAL VARIABLES
900 C
901  include 'coupl.inc'
902 C
903 C COLOR DATA
904 C
905  DATA denom(1)/1/
906  DATA (cf(i, 1),i= 1, 2) / 9, 3/
907 C 1 T(2,1) T(3,4)
908  DATA denom(2)/1/
909  DATA (cf(i, 2),i= 1, 2) / 3, 9/
910 C 1 T(2,4) T(3,1)
911 
912  INTEGER h1,h2
913  REAL*8 matrix
914  matrix_ddx_ccx_h=0.d0
915  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
916  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
917 
918 C ----------
919 C BEGIN CODE
920 C ----------
921  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
922  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
923  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
924  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
925  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
926  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
927  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
928  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
929  CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
930 C Amplitude(s) for diagram number 1
931  CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
932  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,5))
933  CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,3))
934 C Amplitude(s) for diagram number 2
935  CALL vvs1_0(w(1,3),w(1,5),w(1,7),gc_72,amp(2))
936  jamp(1)=-amp(1)
937  jamp(2)=+amp(2)
938 
939  matrix = 0.d0
940  DO i = 1, ncolor
941  ztemp = (0.d0,0.d0)
942  DO j = 1, ncolor
943  ztemp = ztemp + cf(j,i)*jamp(j)
944  ENDDO
945  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
946  ENDDO
947  matrix_ddx_ccx_h=matrix
948  ENDIF
949  ENDIF
950  END
951 
952 C ----- begin subprocesses DDX->UCX with Higgs->tautau
953 
954  SUBROUTINE ddx_ucx_h(P,H1,H2,ANS)
955 C
956 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
957 C By the MadGraph Development Team
958 C Please visit us at https://launchpad.net/madgraph5
959 C
960 C MadGraph StandAlone Version
961 C
962 C Returns amplitude squared summed/avg over colors
963 C and helicities
964 C for the point in phase space P(0:3,NEXTERNAL)
965 C
966 C Process: d d~ > u c~ h WEIGHTED=6
967 C * Decay: h > ta+ ta- WEIGHTED=2
968 C
969  IMPLICIT NONE
970 C
971 C CONSTANTS
972 C
973  INTEGER nexternal
974  parameter(nexternal=6)
975  INTEGER ncomb
976  parameter( ncomb=64)
977 C
978 C ARGUMENTS
979 C
980  REAL*8 p(0:3,nexternal),ans
981  INTEGER h1,h2
982 C
983 C LOCAL VARIABLES
984 C
985  INTEGER nhel(nexternal,ncomb),ntry
986  REAL*8 t
987  REAL*8 matrix_ddx_ucx_h
988  INTEGER ihel,iden, i
989  INTEGER jc(nexternal)
990  LOGICAL goodhel(ncomb)
991  DATA ntry/0/
992  DATA goodhel/ncomb*.false./
993  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
994  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
995  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
996  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
997  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
998  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
999  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1000  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1001  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1002  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1003  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1004  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1005  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1006  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1007  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1008  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1009  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1010  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1011  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1012  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1013  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1014  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1015  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1016  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1017  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1018  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1019  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1020  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1021  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1022  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1023  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1024  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1025  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1026  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1027  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1028  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1029  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1030  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1031  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1032  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1033  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1034  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1035  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1036  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1037  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1038  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1039  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1040  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1041  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1042  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1043  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1044  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1045  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1046  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1047  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1048  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1049  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1050  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1051  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1052  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1053  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1054  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1055  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1056  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1057  DATA iden/36/
1058 C ----------
1059 C BEGIN CODE
1060 C ----------
1061  DO ihel=1,nexternal
1062  jc(ihel) = +1
1063  ENDDO
1064  ans = 0d0
1065  DO ihel=1,ncomb
1066  t=matrix_ddx_ucx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1067  ans=ans+t
1068  ENDDO
1069  ans=ans/dble(iden)
1070  END
1071 
1072 
1073  REAL*8 FUNCTION matrix_ddx_ucx_h(P,H1,H2,NHEL,IC)
1074 C
1075 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1076 C By the MadGraph Development Team
1077 C Please visit us at https://launchpad.net/madgraph5
1078 C
1079 C Returns amplitude squared summed/avg over colors
1080 C for the point with external lines W(0:6,NEXTERNAL)
1081 C
1082 C Process: d d~ > u c~ h WEIGHTED=6
1083 C * Decay: h > ta+ ta- WEIGHTED=2
1084 C
1085  IMPLICIT NONE
1086 C
1087 C CONSTANTS
1088 C
1089  INTEGER ngraphs
1090  parameter(ngraphs=1)
1091  INTEGER nexternal
1092  parameter(nexternal=6)
1093  INTEGER nwavefuncs, ncolor
1094  parameter(nwavefuncs=7, ncolor=1)
1095  REAL*8 zero
1096  parameter(zero=0d0)
1097  COMPLEX*16 imag1
1098  parameter(imag1=(0d0,1d0))
1099 C
1100 C ARGUMENTS
1101 C
1102  REAL*8 p(0:3,nexternal)
1103  INTEGER nhel(nexternal), ic(nexternal)
1104 C
1105 C LOCAL VARIABLES
1106 C
1107  INTEGER i,j
1108  COMPLEX*16 ztemp
1109  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1110  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1111  COMPLEX*16 w(18,nwavefuncs)
1112  COMPLEX*16 dum0,dum1
1113  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1114 C
1115 C GLOBAL VARIABLES
1116 C
1117  include 'coupl.inc'
1118 C
1119 C COLOR DATA
1120 C
1121  DATA denom(1)/1/
1122  DATA (cf(i, 1),i= 1, 1) / 9/
1123 C 1 T(2,4) T(3,1)
1124 
1125  INTEGER h1,h2
1126  REAL*8 matrix
1127  matrix_ddx_ucx_h=0.d0
1128  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1129  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1130 
1131 C ----------
1132 C BEGIN CODE
1133 C ----------
1134  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1135  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1136  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1137  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1138  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1139  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1140  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1141  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
1142  CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,3))
1143 C Amplitude(s) for diagram number 1
1144  CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1145  jamp(1)=+amp(1)
1146 
1147  matrix = 0.d0
1148  DO i = 1, ncolor
1149  ztemp = (0.d0,0.d0)
1150  DO j = 1, ncolor
1151  ztemp = ztemp + cf(j,i)*jamp(j)
1152  ENDDO
1153  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1154  ENDDO
1155  matrix_ddx_ucx_h=matrix
1156  ENDIF
1157  ENDIF
1158  END
1159 
1160 C ------- BEFFINING SUBPROCESS DDX->CUX with H->tautau
1161 
1162  SUBROUTINE ddx_cux_h(P,H1,H2,ANS)
1163 C
1164 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1165 C By the MadGraph Development Team
1166 C Please visit us at https://launchpad.net/madgraph5
1167 C
1168 C MadGraph StandAlone Version
1169 C
1170 C Returns amplitude squared summed/avg over colors
1171 C and helicities
1172 C for the point in phase space P(0:3,NEXTERNAL)
1173 C
1174 C Process: d d~ > c u~ h WEIGHTED=6
1175 C * Decay: h > ta+ ta- WEIGHTED=2
1176 C
1177  IMPLICIT NONE
1178 C
1179 C CONSTANTS
1180 C
1181  INTEGER nexternal
1182  parameter(nexternal=6)
1183  INTEGER ncomb
1184  parameter( ncomb=64)
1185 C
1186 C ARGUMENTS
1187 C
1188  REAL*8 p(0:3,nexternal),ans
1189  INTEGER h1,h2
1190 C
1191 C LOCAL VARIABLES
1192 C
1193  INTEGER nhel(nexternal,ncomb),ntry
1194  REAL*8 t
1195  REAL*8 matrix_ddx_cux_h
1196  INTEGER ihel,iden, i
1197  INTEGER jc(nexternal)
1198  LOGICAL goodhel(ncomb)
1199  DATA ntry/0/
1200  DATA goodhel/ncomb*.false./
1201  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1202  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1203  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1204  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1205  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1206  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1207  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1208  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1209  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1210  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1211  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1212  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1213  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1214  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1215  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1216  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1217  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1218  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1219  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1220  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1221  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1222  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1223  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1224  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1225  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1226  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1227  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1228  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1229  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1230  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1231  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1232  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1233  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1234  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1235  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1236  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1237  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1238  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1239  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1240  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1241  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1242  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1243  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1244  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1245  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1246  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1247  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1248  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1249  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1250  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1251  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1252  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1253  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1254  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1255  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1256  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1257  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1258  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1259  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1260  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1261  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1262  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1263  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1264  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1265  DATA iden/36/
1266 C ----------
1267 C BEGIN CODE
1268 C ----------
1269  DO ihel=1,nexternal
1270  jc(ihel) = +1
1271  ENDDO
1272  ans = 0d0
1273  DO ihel=1,ncomb
1274  t=matrix_ddx_cux_h(p ,h1,h2,nhel(1,ihel),jc(1))
1275  ans=ans+t
1276  ENDDO
1277  ans=ans/dble(iden)
1278  END
1279 
1280 
1281  REAL*8 FUNCTION matrix_ddx_cux_h(P,H1,H2,NHEL,IC)
1282 C
1283 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1284 C By the MadGraph Development Team
1285 C Please visit us at https://launchpad.net/madgraph5
1286 C
1287 C Returns amplitude squared summed/avg over colors
1288 C for the point with external lines W(0:6,NEXTERNAL)
1289 C
1290 C Process: d d~ > c u~ h WEIGHTED=6
1291 C * Decay: h > ta+ ta- WEIGHTED=2
1292 C
1293  IMPLICIT NONE
1294 C
1295 C CONSTANTS
1296 C
1297  INTEGER ngraphs
1298  parameter(ngraphs=1)
1299  INTEGER nexternal
1300  parameter(nexternal=6)
1301  INTEGER nwavefuncs, ncolor
1302  parameter(nwavefuncs=7, ncolor=1)
1303  REAL*8 zero
1304  parameter(zero=0d0)
1305  COMPLEX*16 imag1
1306  parameter(imag1=(0d0,1d0))
1307 C
1308 C ARGUMENTS
1309 C
1310  REAL*8 p(0:3,nexternal)
1311  INTEGER nhel(nexternal), ic(nexternal)
1312 C
1313 C LOCAL VARIABLES
1314 C
1315  INTEGER i,j
1316  COMPLEX*16 ztemp
1317  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1318  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1319  COMPLEX*16 w(18,nwavefuncs)
1320  COMPLEX*16 dum0,dum1
1321  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1322 C
1323 C GLOBAL VARIABLES
1324 C
1325  include 'coupl.inc'
1326 C
1327 C COLOR DATA
1328 C
1329  DATA denom(1)/1/
1330  DATA (cf(i, 1),i= 1, 1) / 9/
1331 C 1 T(2,4) T(3,1)
1332 
1333  INTEGER h1,h2
1334  REAL*8 matrix
1335  matrix_ddx_cux_h=0.d0
1336  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1337  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1338 
1339 C ----------
1340 C BEGIN CODE
1341 C ----------
1342  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1343  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1344  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1345  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1346  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1347  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1348  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1349  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,6))
1350  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
1351 C Amplitude(s) for diagram number 1
1352  CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1353  jamp(1)=+amp(1)
1354 
1355  matrix = 0.d0
1356  DO i = 1, ncolor
1357  ztemp = (0.d0,0.d0)
1358  DO j = 1, ncolor
1359  ztemp = ztemp + cf(j,i)*jamp(j)
1360  ENDDO
1361  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1362  ENDDO
1363  matrix_ddx_cux_h=matrix
1364  ENDIF
1365  ENDIF
1366  END
1367 
1368 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
1369 C ----- begin subprocesses ud->jjtautau, no Higgs
1370 
1371 C-----------SUBPROCESS DDX->DDX NO HIGGS
1372 
1373  SUBROUTINE ddx_ddx_noh(P,H1,H2,ANS)
1374 C
1375 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1376 C By the MadGraph Development Team
1377 C Please visit us at https://launchpad.net/madgraph5
1378 C
1379 C MadGraph StandAlone Version
1380 C
1381 C Returns amplitude squared summed/avg over colors
1382 C and helicities
1383 C for the point in phase space P(0:3,NEXTERNAL)
1384 C
1385 C Process: d d~ > d d~ ta+ ta- / h QED=4
1386 C
1387  IMPLICIT NONE
1388 C
1389 C CONSTANTS
1390 C
1391  INTEGER nexternal
1392  parameter(nexternal=6)
1393  INTEGER ncomb
1394  parameter( ncomb=64)
1395 C
1396 C ARGUMENTS
1397 C
1398  REAL*8 p(0:3,nexternal),ans
1399  INTEGER h1,h2
1400 C
1401 C LOCAL VARIABLES
1402 C
1403  INTEGER nhel(nexternal,ncomb),ntry
1404  REAL*8 t
1405  REAL*8 matrix_ddx_ddx_noh
1406  INTEGER ihel,iden, i
1407  INTEGER jc(nexternal)
1408  LOGICAL goodhel(ncomb)
1409  DATA ntry/0/
1410  DATA goodhel/ncomb*.false./
1411  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1412  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1413  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1414  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1415  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1416  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1417  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1418  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1419  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1420  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1421  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1422  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1423  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1424  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1425  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1426  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1427  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1428  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1429  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1430  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1431  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1432  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1433  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1434  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1435  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1436  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1437  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1438  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1439  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1440  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1441  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1442  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1443  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1444  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1445  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1446  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1447  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1448  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1449  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1450  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1451  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1452  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1453  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1454  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1455  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1456  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1457  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1458  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1459  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1460  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1461  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1462  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1463  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1464  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1465  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1466  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1467  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1468  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1469  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1470  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1471  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1472  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1473  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1474  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1475  DATA iden/36/
1476 C ----------
1477 C BEGIN CODE
1478 C ----------
1479  DO ihel=1,nexternal
1480  jc(ihel) = +1
1481  ENDDO
1482  ans = 0d0
1483  DO ihel=1,ncomb
1484  t=matrix_ddx_ddx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1485  ans=ans+t
1486  ENDDO
1487  ans=ans/dble(iden)
1488  END
1489 
1490 
1491  REAL*8 FUNCTION matrix_ddx_ddx_noh(P,H1,H2,NHEL,IC)
1492 C
1493 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1494 C By the MadGraph Development Team
1495 C Please visit us at https://launchpad.net/madgraph5
1496 C
1497 C Returns amplitude squared summed/avg over colors
1498 C for the point with external lines W(0:6,NEXTERNAL)
1499 C
1500 C Process: d d~ > d d~ ta+ ta- / h QED=4
1501 C
1502  IMPLICIT NONE
1503 C
1504 C CONSTANTS
1505 C
1506  INTEGER ngraphs
1507  parameter(ngraphs=64)
1508  INTEGER nexternal
1509  parameter(nexternal=6)
1510  INTEGER nwavefuncs, ncolor
1511  parameter(nwavefuncs=16, ncolor=2)
1512  REAL*8 zero
1513  parameter(zero=0d0)
1514  COMPLEX*16 imag1
1515  parameter(imag1=(0d0,1d0))
1516 C
1517 C ARGUMENTS
1518 C
1519  REAL*8 p(0:3,nexternal)
1520  INTEGER nhel(nexternal), ic(nexternal)
1521 C
1522 C LOCAL VARIABLES
1523 C
1524  INTEGER i,j
1525  COMPLEX*16 ztemp
1526  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1527  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1528  COMPLEX*16 w(18,nwavefuncs)
1529  COMPLEX*16 dum0,dum1
1530  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1531 C
1532 C GLOBAL VARIABLES
1533 C
1534  include 'coupl.inc'
1535 C
1536 C COLOR DATA
1537 C
1538  DATA denom(1)/1/
1539  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1540 C 1 T(2,1) T(3,4)
1541  DATA denom(2)/1/
1542  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1543 C 1 T(2,4) T(3,1)
1544 
1545  INTEGER h1,h2
1546  REAL*8 matrix
1547  matrix_ddx_ddx_noh=0.d0
1548  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1549  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1550 
1551 C ----------
1552 C BEGIN CODE
1553 C ----------
1554  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1555  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1556  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1557  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1558  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1559  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1560  CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
1561  CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
1562  CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1563 C Amplitude(s) for diagram number 1
1564  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1565  CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1566 C Amplitude(s) for diagram number 2
1567  CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1568  CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
1569 C Amplitude(s) for diagram number 3
1570  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1571 C Amplitude(s) for diagram number 4
1572  CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1573  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
1574  CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1575 C Amplitude(s) for diagram number 5
1576  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1577  CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1578 C Amplitude(s) for diagram number 6
1579  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1580 C Amplitude(s) for diagram number 7
1581  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1582 C Amplitude(s) for diagram number 8
1583  CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1584  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1585  CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
1586 C Amplitude(s) for diagram number 9
1587  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
1588  CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1589 C Amplitude(s) for diagram number 10
1590  CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
1591  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1592 C Amplitude(s) for diagram number 11
1593  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
1594 C Amplitude(s) for diagram number 12
1595  CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
1596  CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1597  CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1598 C Amplitude(s) for diagram number 13
1599  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
1600  CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1601 C Amplitude(s) for diagram number 14
1602  CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(14))
1603 C Amplitude(s) for diagram number 15
1604  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
1605 C Amplitude(s) for diagram number 16
1606  CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(16))
1607  CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
1608 C Amplitude(s) for diagram number 17
1609  CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_1,amp(17))
1610  CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1611 C Amplitude(s) for diagram number 18
1612  CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
1613 C Amplitude(s) for diagram number 19
1614  CALL ffv2_3_0(w(1,4),w(1,14),w(1,7),gc_50,gc_58,amp(19))
1615 C Amplitude(s) for diagram number 20
1616  CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
1617  CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,9))
1618  CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,14))
1619  CALL ffv1_2(w(1,5),w(1,9),gc_3,mta,zero,w(1,10))
1620 C Amplitude(s) for diagram number 21
1621  CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(21))
1622  CALL ffv1_1(w(1,6),w(1,9),gc_3,mta,zero,w(1,13))
1623 C Amplitude(s) for diagram number 22
1624  CALL ffv1_0(w(1,5),w(1,13),w(1,14),gc_3,amp(22))
1625  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,15))
1626 C Amplitude(s) for diagram number 23
1627  CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(23))
1628 C Amplitude(s) for diagram number 24
1629  CALL ffv2_4_0(w(1,5),w(1,13),w(1,15),gc_50,gc_59,amp(24))
1630  CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,13))
1631  CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,10))
1632 C Amplitude(s) for diagram number 25
1633  CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(25))
1634  CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,16))
1635 C Amplitude(s) for diagram number 26
1636  CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
1637 C Amplitude(s) for diagram number 27
1638  CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(27))
1639 C Amplitude(s) for diagram number 28
1640  CALL ffv2_4_0(w(1,5),w(1,16),w(1,15),gc_50,gc_59,amp(28))
1641  CALL ffv1_1(w(1,2),w(1,9),gc_1,zero,zero,w(1,16))
1642 C Amplitude(s) for diagram number 29
1643  CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_1,amp(29))
1644  CALL ffv1_2(w(1,4),w(1,9),gc_1,zero,zero,w(1,5))
1645 C Amplitude(s) for diagram number 30
1646  CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(30))
1647 C Amplitude(s) for diagram number 31
1648  CALL ffv2_3_0(w(1,4),w(1,16),w(1,7),gc_50,gc_58,amp(31))
1649 C Amplitude(s) for diagram number 32
1650  CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(32))
1651  CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1652  CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,16))
1653 C Amplitude(s) for diagram number 33
1654  CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_1,amp(33))
1655  CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,9))
1656 C Amplitude(s) for diagram number 34
1657  CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(34))
1658 C Amplitude(s) for diagram number 35
1659  CALL ffv2_3_0(w(1,4),w(1,16),w(1,7),gc_50,gc_58,amp(35))
1660 C Amplitude(s) for diagram number 36
1661  CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(36))
1662  CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,9))
1663 C Amplitude(s) for diagram number 37
1664  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(37))
1665  CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,16))
1666 C Amplitude(s) for diagram number 38
1667  CALL ffv1_0(w(1,16),w(1,2),w(1,12),gc_1,amp(38))
1668 C Amplitude(s) for diagram number 39
1669  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(39))
1670 C Amplitude(s) for diagram number 40
1671  CALL ffv2_3_0(w(1,16),w(1,2),w(1,7),gc_50,gc_58,amp(40))
1672  CALL ffv1_2(w(1,1),w(1,14),gc_1,zero,zero,w(1,16))
1673 C Amplitude(s) for diagram number 41
1674  CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_1,amp(41))
1675  CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,9))
1676 C Amplitude(s) for diagram number 42
1677  CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_1,amp(42))
1678 C Amplitude(s) for diagram number 43
1679  CALL ffv2_3_0(w(1,16),w(1,3),w(1,7),gc_50,gc_58,amp(43))
1680  CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,16))
1681 C Amplitude(s) for diagram number 44
1682  CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_1,amp(44))
1683  CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,14))
1684  CALL ffv1_2(w(1,1),w(1,14),gc_11,zero,zero,w(1,13))
1685 C Amplitude(s) for diagram number 45
1686  CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(45))
1687 C Amplitude(s) for diagram number 46
1688  CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_11,amp(46))
1689 C Amplitude(s) for diagram number 47
1690  CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(47))
1691 C Amplitude(s) for diagram number 48
1692  CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_11,amp(48))
1693  CALL ffv2_3_2(w(1,1),w(1,15),gc_50,gc_58,zero,zero,w(1,14))
1694 C Amplitude(s) for diagram number 49
1695  CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(49))
1696 C Amplitude(s) for diagram number 50
1697  CALL ffv2_3_0(w(1,9),w(1,3),w(1,15),gc_50,gc_58,amp(50))
1698 C Amplitude(s) for diagram number 51
1699  CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(51))
1700 C Amplitude(s) for diagram number 52
1701  CALL ffv2_3_0(w(1,16),w(1,3),w(1,15),gc_50,gc_58,amp(52))
1702  CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,15))
1703 C Amplitude(s) for diagram number 53
1704  CALL ffv1_0(w(1,15),w(1,2),w(1,12),gc_1,amp(53))
1705 C Amplitude(s) for diagram number 54
1706  CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_1,amp(54))
1707 C Amplitude(s) for diagram number 55
1708  CALL ffv2_3_0(w(1,15),w(1,2),w(1,7),gc_50,gc_58,amp(55))
1709 C Amplitude(s) for diagram number 56
1710  CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_1,amp(56))
1711  CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
1712  CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1713 C Amplitude(s) for diagram number 57
1714  CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(57))
1715 C Amplitude(s) for diagram number 58
1716  CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(58))
1717 C Amplitude(s) for diagram number 59
1718  CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(59))
1719 C Amplitude(s) for diagram number 60
1720  CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_11,amp(60))
1721  CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
1722 C Amplitude(s) for diagram number 61
1723  CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(61))
1724 C Amplitude(s) for diagram number 62
1725  CALL ffv2_3_0(w(1,9),w(1,2),w(1,11),gc_50,gc_58,amp(62))
1726 C Amplitude(s) for diagram number 63
1727  CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(63))
1728 C Amplitude(s) for diagram number 64
1729  CALL ffv2_3_0(w(1,16),w(1,2),w(1,11),gc_50,gc_58,amp(64))
1730  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1731  $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1732  $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1733  $ +1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)+1./2.*amp(36)
1734  $ +1./2.*amp(45)+1./2.*amp(46)+1./2.*amp(47)+1./2.*amp(48)
1735  $ -amp(53)-amp(54)-amp(55)-amp(56)+1./6.*amp(57)+1./6.*amp(58)
1736  $ +1./6.*amp(59)+1./6.*amp(60)-amp(61)-amp(62)-amp(63)-amp(64)
1737  jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1738  $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1739  $ +amp(28)+amp(29)+amp(30)+amp(31)+amp(32)-1./6.*amp(33)
1740  $ -1./6.*amp(34)-1./6.*amp(35)-1./6.*amp(36)+amp(37)+amp(38)
1741  $ +amp(39)+amp(40)+amp(41)+amp(42)+amp(43)+amp(44)-1./6.*amp(45)
1742  $ -1./6.*amp(46)-1./6.*amp(47)-1./6.*amp(48)+amp(49)+amp(50)
1743  $ +amp(51)+amp(52)-1./2.*amp(57)-1./2.*amp(58)-1./2.*amp(59)
1744  $ -1./2.*amp(60)
1745 
1746  matrix = 0.d0
1747  DO i = 1, ncolor
1748  ztemp = (0.d0,0.d0)
1749  DO j = 1, ncolor
1750  ztemp = ztemp + cf(j,i)*jamp(j)
1751  ENDDO
1752  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1753  ENDDO
1754  matrix_ddx_ddx_noh=matrix
1755  ENDIF
1756  ENDIF
1757  END
1758 
1759 C-----------SUBPROCESS DDX->UUX NO HIGGS
1760 
1761  SUBROUTINE ddx_uux_noh(P,H1,H2,ANS)
1762 C
1763 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1764 C By the MadGraph Development Team
1765 C Please visit us at https://launchpad.net/madgraph5
1766 C
1767 C MadGraph StandAlone Version
1768 C
1769 C Returns amplitude squared summed/avg over colors
1770 C and helicities
1771 C for the point in phase space P(0:3,NEXTERNAL)
1772 C
1773 C Process: d d~ > u u~ ta+ ta- / h QED=4
1774 C
1775  IMPLICIT NONE
1776 C
1777 C CONSTANTS
1778 C
1779  INTEGER nexternal
1780  parameter(nexternal=6)
1781  INTEGER ncomb
1782  parameter( ncomb=64)
1783 C
1784 C ARGUMENTS
1785 C
1786  REAL*8 p(0:3,nexternal),ans
1787  INTEGER h1,h2
1788 C
1789 C LOCAL VARIABLES
1790 C
1791  INTEGER nhel(nexternal,ncomb),ntry
1792  REAL*8 t
1793  REAL*8 matrix_ddx_uux_noh
1794  INTEGER ihel,iden, i
1795  INTEGER jc(nexternal)
1796  LOGICAL goodhel(ncomb)
1797  DATA ntry/0/
1798  DATA goodhel/ncomb*.false./
1799  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1800  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1801  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1802  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1803  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1804  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1805  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1806  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1807  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1808  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1809  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1810  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1811  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1812  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1813  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1814  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1815  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1816  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1817  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1818  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1819  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1820  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1821  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1822  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1823  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1824  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1825  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1826  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1827  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1828  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1829  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1830  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1831  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1832  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1833  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1834  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1835  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1836  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1837  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1838  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1839  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1840  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1841  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1842  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1843  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1844  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1845  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1846  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1847  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1848  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1849  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1850  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1851  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1852  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1853  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1854  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1855  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1856  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1857  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1858  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1859  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1860  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1861  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1862  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1863  DATA iden/36/
1864 C ----------
1865 C BEGIN CODE
1866 C ----------
1867  DO ihel=1,nexternal
1868  jc(ihel) = +1
1869  ENDDO
1870  ans = 0d0
1871  DO ihel=1,ncomb
1872  t=matrix_ddx_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1873  ans=ans+t
1874  ENDDO
1875  ans=ans/dble(iden)
1876  END
1877 
1878 
1879  REAL*8 FUNCTION matrix_ddx_uux_noh(P,H1,H2,NHEL,IC)
1880 C
1881 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1882 C By the MadGraph Development Team
1883 C Please visit us at https://launchpad.net/madgraph5
1884 C
1885 C Returns amplitude squared summed/avg over colors
1886 C for the point with external lines W(0:6,NEXTERNAL)
1887 C
1888 C Process: d d~ > u u~ ta+ ta- / h QED=4
1889 C
1890  IMPLICIT NONE
1891 C
1892 C CONSTANTS
1893 C
1894  INTEGER ngraphs
1895  parameter(ngraphs=43)
1896  INTEGER nexternal
1897  parameter(nexternal=6)
1898  INTEGER nwavefuncs, ncolor
1899  parameter(nwavefuncs=14, ncolor=2)
1900  REAL*8 zero
1901  parameter(zero=0d0)
1902  COMPLEX*16 imag1
1903  parameter(imag1=(0d0,1d0))
1904 C
1905 C ARGUMENTS
1906 C
1907  REAL*8 p(0:3,nexternal)
1908  INTEGER nhel(nexternal), ic(nexternal)
1909 C
1910 C LOCAL VARIABLES
1911 C
1912  INTEGER i,j
1913  COMPLEX*16 ztemp
1914  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1915  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1916  COMPLEX*16 w(18,nwavefuncs)
1917  COMPLEX*16 dum0,dum1
1918  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1919 C
1920 C GLOBAL VARIABLES
1921 C
1922  include 'coupl.inc'
1923 C
1924 C COLOR DATA
1925 C
1926  DATA denom(1)/1/
1927  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1928 C 1 T(2,1) T(3,4)
1929  DATA denom(2)/1/
1930  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1931 C 1 T(2,4) T(3,1)
1932 
1933  INTEGER h1,h2
1934  REAL*8 matrix
1935  matrix_ddx_uux_noh=0.d0
1936  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1937  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1938 
1939 C ----------
1940 C BEGIN CODE
1941 C ----------
1942  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1943  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1944  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1945  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1946  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1947  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1948  CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
1949  CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
1950  CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1951 C Amplitude(s) for diagram number 1
1952  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1953  CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1954 C Amplitude(s) for diagram number 2
1955  CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1956  CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
1957 C Amplitude(s) for diagram number 3
1958  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1959 C Amplitude(s) for diagram number 4
1960  CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1961  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
1962  CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1963 C Amplitude(s) for diagram number 5
1964  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1965  CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1966 C Amplitude(s) for diagram number 6
1967  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1968 C Amplitude(s) for diagram number 7
1969  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1970 C Amplitude(s) for diagram number 8
1971  CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1972  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1973  CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
1974 C Amplitude(s) for diagram number 9
1975  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
1976  CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
1977 C Amplitude(s) for diagram number 10
1978  CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
1979  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1980 C Amplitude(s) for diagram number 11
1981  CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
1982 C Amplitude(s) for diagram number 12
1983  CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
1984  CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1985  CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1986 C Amplitude(s) for diagram number 13
1987  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
1988  CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1989 C Amplitude(s) for diagram number 14
1990  CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
1991 C Amplitude(s) for diagram number 15
1992  CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
1993 C Amplitude(s) for diagram number 16
1994  CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
1995  CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
1996 C Amplitude(s) for diagram number 17
1997  CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
1998  CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
1999 C Amplitude(s) for diagram number 18
2000  CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
2001 C Amplitude(s) for diagram number 19
2002  CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
2003 C Amplitude(s) for diagram number 20
2004  CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
2005  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,9))
2006  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,14))
2007  CALL ffv2_1(w(1,6),w(1,9),gc_108,zero,zero,w(1,10))
2008 C Amplitude(s) for diagram number 21
2009  CALL ffv2_0(w(1,5),w(1,10),w(1,14),gc_108,amp(21))
2010 C Amplitude(s) for diagram number 22
2011  CALL vvv1_0(w(1,12),w(1,14),w(1,9),gc_4,amp(22))
2012 C Amplitude(s) for diagram number 23
2013  CALL vvv1_0(w(1,14),w(1,9),w(1,7),gc_53,amp(23))
2014  CALL ffv2_1(w(1,2),w(1,9),gc_100,zero,zero,w(1,10))
2015 C Amplitude(s) for diagram number 24
2016  CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_2,amp(24))
2017  CALL ffv2_2(w(1,4),w(1,9),gc_100,zero,zero,w(1,5))
2018 C Amplitude(s) for diagram number 25
2019  CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(25))
2020 C Amplitude(s) for diagram number 26
2021  CALL ffv2_5_0(w(1,4),w(1,10),w(1,7),gc_51,gc_58,amp(26))
2022 C Amplitude(s) for diagram number 27
2023  CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2024  CALL ffv2_2(w(1,1),w(1,14),gc_100,zero,zero,w(1,5))
2025 C Amplitude(s) for diagram number 28
2026  CALL ffv1_0(w(1,5),w(1,3),w(1,12),gc_2,amp(28))
2027  CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,10))
2028 C Amplitude(s) for diagram number 29
2029  CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_100,amp(29))
2030 C Amplitude(s) for diagram number 30
2031  CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(30))
2032  CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,5))
2033 C Amplitude(s) for diagram number 31
2034  CALL ffv2_0(w(1,5),w(1,3),w(1,14),gc_100,amp(31))
2035  CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,14))
2036 C Amplitude(s) for diagram number 32
2037  CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_1,amp(32))
2038 C Amplitude(s) for diagram number 33
2039  CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(33))
2040 C Amplitude(s) for diagram number 34
2041  CALL ffv2_3_0(w(1,14),w(1,2),w(1,7),gc_50,gc_58,amp(34))
2042 C Amplitude(s) for diagram number 35
2043  CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(35))
2044  CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2045  CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2046 C Amplitude(s) for diagram number 36
2047  CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(36))
2048 C Amplitude(s) for diagram number 37
2049  CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
2050 C Amplitude(s) for diagram number 38
2051  CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(38))
2052 C Amplitude(s) for diagram number 39
2053  CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(39))
2054  CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2055 C Amplitude(s) for diagram number 40
2056  CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(40))
2057 C Amplitude(s) for diagram number 41
2058  CALL ffv2_3_0(w(1,10),w(1,2),w(1,11),gc_50,gc_58,amp(41))
2059 C Amplitude(s) for diagram number 42
2060  CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(42))
2061 C Amplitude(s) for diagram number 43
2062  CALL ffv2_3_0(w(1,5),w(1,2),w(1,11),gc_50,gc_58,amp(43))
2063  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2064  $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2065  $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2066  $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
2067  $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
2068  jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2069  $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2070  $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
2071  $ -1./2.*amp(38)-1./2.*amp(39)
2072 
2073  matrix = 0.d0
2074  DO i = 1, ncolor
2075  ztemp = (0.d0,0.d0)
2076  DO j = 1, ncolor
2077  ztemp = ztemp + cf(j,i)*jamp(j)
2078  ENDDO
2079  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2080  ENDDO
2081  matrix_ddx_uux_noh=matrix
2082  ENDIF
2083  ENDIF
2084  END
2085 
2086 C-----------SUBPROCESS DDX->SSX NO HIGGS
2087 
2088  SUBROUTINE ddx_ssx_noh(P,H1,H2,ANS)
2089 C
2090 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2091 C By the MadGraph Development Team
2092 C Please visit us at https://launchpad.net/madgraph5
2093 C
2094 C MadGraph StandAlone Version
2095 C
2096 C Returns amplitude squared summed/avg over colors
2097 C and helicities
2098 C for the point in phase space P(0:3,NEXTERNAL)
2099 C
2100 C Process: d d~ > s s~ ta+ ta- / h QED=4
2101 C
2102  IMPLICIT NONE
2103 C
2104 C CONSTANTS
2105 C
2106  INTEGER nexternal
2107  parameter(nexternal=6)
2108  INTEGER ncomb
2109  parameter( ncomb=64)
2110 C
2111 C ARGUMENTS
2112 C
2113  REAL*8 p(0:3,nexternal),ans
2114  INTEGER h1,h2
2115 C
2116 C LOCAL VARIABLES
2117 C
2118  INTEGER nhel(nexternal,ncomb),ntry
2119  REAL*8 t
2120  REAL*8 matrix_ddx_ssx_noh
2121  INTEGER ihel,iden, i
2122  INTEGER jc(nexternal)
2123  LOGICAL goodhel(ncomb)
2124  DATA ntry/0/
2125  DATA goodhel/ncomb*.false./
2126  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2127  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2128  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2129  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2130  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2131  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2132  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2133  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2134  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2135  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2136  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2137  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2138  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2139  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2140  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2141  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2142  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2143  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2144  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2145  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2146  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2147  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2148  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2149  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2150  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2151  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2152  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2153  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2154  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2155  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2156  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2157  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2158  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2159  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2160  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2161  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2162  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2163  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2164  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2165  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2166  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2167  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2168  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2169  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2170  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2171  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2172  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2173  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2174  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2175  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2176  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2177  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2178  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2179  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2180  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2181  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2182  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2183  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2184  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2185  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2186  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2187  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2188  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2189  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2190  DATA iden/36/
2191 C ----------
2192 C BEGIN CODE
2193 C ----------
2194  DO ihel=1,nexternal
2195  jc(ihel) = +1
2196  ENDDO
2197  ans = 0d0
2198  DO ihel=1,ncomb
2199  t=matrix_ddx_ssx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2200  ans=ans+t
2201  ENDDO
2202  ans=ans/dble(iden)
2203  END
2204 
2205 
2206  REAL*8 FUNCTION matrix_ddx_ssx_noh(P,H1,H2,NHEL,IC)
2207 C
2208 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2209 C By the MadGraph Development Team
2210 C Please visit us at https://launchpad.net/madgraph5
2211 C
2212 C Returns amplitude squared summed/avg over colors
2213 C for the point with external lines W(0:6,NEXTERNAL)
2214 C
2215 C Process: d d~ > s s~ ta+ ta- / h QED=4
2216 C
2217  IMPLICIT NONE
2218 C
2219 C CONSTANTS
2220 C
2221  INTEGER ngraphs
2222  parameter(ngraphs=32)
2223  INTEGER nexternal
2224  parameter(nexternal=6)
2225  INTEGER nwavefuncs, ncolor
2226  parameter(nwavefuncs=13, ncolor=2)
2227  REAL*8 zero
2228  parameter(zero=0d0)
2229  COMPLEX*16 imag1
2230  parameter(imag1=(0d0,1d0))
2231 C
2232 C ARGUMENTS
2233 C
2234  REAL*8 p(0:3,nexternal)
2235  INTEGER nhel(nexternal), ic(nexternal)
2236 C
2237 C LOCAL VARIABLES
2238 C
2239  INTEGER i,j
2240  COMPLEX*16 ztemp
2241  REAL*8 denom(ncolor), cf(ncolor,ncolor)
2242  COMPLEX*16 amp(ngraphs), jamp(ncolor)
2243  COMPLEX*16 w(18,nwavefuncs)
2244  COMPLEX*16 dum0,dum1
2245  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2246 C
2247 C GLOBAL VARIABLES
2248 C
2249  include 'coupl.inc'
2250 C
2251 C COLOR DATA
2252 C
2253  DATA denom(1)/1/
2254  DATA (cf(i, 1),i= 1, 2) / 9, 3/
2255 C 1 T(2,1) T(3,4)
2256  DATA denom(2)/1/
2257  DATA (cf(i, 2),i= 1, 2) / 3, 9/
2258 C 1 T(2,4) T(3,1)
2259 
2260  INTEGER h1,h2
2261  REAL*8 matrix
2262  matrix_ddx_ssx_noh=0.d0
2263  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2264  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2265 
2266 C ----------
2267 C BEGIN CODE
2268 C ----------
2269  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2270  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2271  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2272  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2273  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2274  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2275  CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
2276  CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
2277  CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2278 C Amplitude(s) for diagram number 1
2279  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2280  CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2281 C Amplitude(s) for diagram number 2
2282  CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2283  CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
2284 C Amplitude(s) for diagram number 3
2285  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2286 C Amplitude(s) for diagram number 4
2287  CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2288  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
2289  CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2290 C Amplitude(s) for diagram number 5
2291  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2292  CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2293 C Amplitude(s) for diagram number 6
2294  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2295 C Amplitude(s) for diagram number 7
2296  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2297 C Amplitude(s) for diagram number 8
2298  CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2299  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2300  CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
2301 C Amplitude(s) for diagram number 9
2302  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
2303  CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
2304 C Amplitude(s) for diagram number 10
2305  CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
2306  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2307 C Amplitude(s) for diagram number 11
2308  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
2309 C Amplitude(s) for diagram number 12
2310  CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
2311  CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2312  CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2313 C Amplitude(s) for diagram number 13
2314  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
2315  CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
2316 C Amplitude(s) for diagram number 14
2317  CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_1,amp(14))
2318 C Amplitude(s) for diagram number 15
2319  CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
2320 C Amplitude(s) for diagram number 16
2321  CALL ffv2_3_0(w(1,6),w(1,3),w(1,7),gc_50,gc_58,amp(16))
2322  CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
2323 C Amplitude(s) for diagram number 17
2324  CALL ffv1_0(w(1,4),w(1,6),w(1,12),gc_1,amp(17))
2325  CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
2326 C Amplitude(s) for diagram number 18
2327  CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
2328 C Amplitude(s) for diagram number 19
2329  CALL ffv2_3_0(w(1,4),w(1,6),w(1,7),gc_50,gc_58,amp(19))
2330 C Amplitude(s) for diagram number 20
2331  CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
2332  CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
2333 C Amplitude(s) for diagram number 21
2334  CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(21))
2335  CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
2336 C Amplitude(s) for diagram number 22
2337  CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_1,amp(22))
2338 C Amplitude(s) for diagram number 23
2339  CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(23))
2340  CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
2341 C Amplitude(s) for diagram number 24
2342  CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_1,amp(24))
2343  CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2344  CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2345 C Amplitude(s) for diagram number 25
2346  CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(25))
2347 C Amplitude(s) for diagram number 26
2348  CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(26))
2349 C Amplitude(s) for diagram number 27
2350  CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2351 C Amplitude(s) for diagram number 28
2352  CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(28))
2353  CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2354 C Amplitude(s) for diagram number 29
2355  CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(29))
2356 C Amplitude(s) for diagram number 30
2357  CALL ffv2_3_0(w(1,6),w(1,2),w(1,11),gc_50,gc_58,amp(30))
2358 C Amplitude(s) for diagram number 31
2359  CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(31))
2360 C Amplitude(s) for diagram number 32
2361  CALL ffv2_3_0(w(1,9),w(1,2),w(1,11),gc_50,gc_58,amp(32))
2362  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2363  $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2364  $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2365  $ -amp(21)-amp(22)-amp(23)-amp(24)+1./6.*amp(25)+1./6.*amp(26)
2366  $ +1./6.*amp(27)+1./6.*amp(28)-amp(29)-amp(30)-amp(31)-amp(32)
2367  jamp(2)=+1./2.*(-amp(13)-amp(14)-amp(15)-amp(16)-amp(25)-amp(26)
2368  $ -amp(27)-amp(28))
2369 
2370  matrix = 0.d0
2371  DO i = 1, ncolor
2372  ztemp = (0.d0,0.d0)
2373  DO j = 1, ncolor
2374  ztemp = ztemp + cf(j,i)*jamp(j)
2375  ENDDO
2376  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2377  ENDDO
2378  matrix_ddx_ssx_noh=matrix
2379  ENDIF
2380  ENDIF
2381  END
2382 
2383 C-----------SUBPROCESS DDX->CCX NO HIGGS
2384 
2385  SUBROUTINE ddx_ccx_noh(P,H1,H2,ANS)
2386 C
2387 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2388 C By the MadGraph Development Team
2389 C Please visit us at https://launchpad.net/madgraph5
2390 C
2391 C MadGraph StandAlone Version
2392 C
2393 C Returns amplitude squared summed/avg over colors
2394 C and helicities
2395 C for the point in phase space P(0:3,NEXTERNAL)
2396 C
2397 C Process: d d~ > c c~ ta+ ta- / h QED=4
2398 C
2399  IMPLICIT NONE
2400 C
2401 C CONSTANTS
2402 C
2403  INTEGER nexternal
2404  parameter(nexternal=6)
2405  INTEGER ncomb
2406  parameter( ncomb=64)
2407 C
2408 C ARGUMENTS
2409 C
2410  REAL*8 p(0:3,nexternal),ans
2411  INTEGER h1,h2
2412 C
2413 C LOCAL VARIABLES
2414 C
2415  INTEGER nhel(nexternal,ncomb),ntry
2416  REAL*8 t
2417  REAL*8 matrix_ddx_ccx_noh
2418  INTEGER ihel,iden, i
2419  INTEGER jc(nexternal)
2420  LOGICAL goodhel(ncomb)
2421  DATA ntry/0/
2422  DATA goodhel/ncomb*.false./
2423  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2424  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2425  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2426  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2427  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2428  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2429  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2430  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2431  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2432  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2433  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2434  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2435  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2436  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2437  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2438  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2439  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2440  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2441  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2442  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2443  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2444  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2445  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2446  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2447  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2448  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2449  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2450  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2451  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2452  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2453  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2454  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2455  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2456  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2457  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2458  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2459  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2460  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2461  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2462  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2463  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2464  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2465  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2466  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2467  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2468  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2469  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2470  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2471  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2472  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2473  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2474  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2475  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2476  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2477  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2478  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2479  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2480  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2481  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2482  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2483  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2484  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2485  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2486  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2487  DATA iden/36/
2488 C ----------
2489 C BEGIN CODE
2490 C ----------
2491  DO ihel=1,nexternal
2492  jc(ihel) = +1
2493  ENDDO
2494  ans = 0d0
2495  DO ihel=1,ncomb
2496  t=matrix_ddx_ccx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2497  ans=ans+t
2498  ENDDO
2499  ans=ans/dble(iden)
2500  END
2501 
2502 
2503  REAL*8 FUNCTION matrix_ddx_ccx_noh(P,H1,H2,NHEL,IC)
2504 C
2505 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2506 C By the MadGraph Development Team
2507 C Please visit us at https://launchpad.net/madgraph5
2508 C
2509 C Returns amplitude squared summed/avg over colors
2510 C for the point with external lines W(0:6,NEXTERNAL)
2511 C
2512 C Process: d d~ > c c~ ta+ ta- / h QED=4
2513 C
2514  IMPLICIT NONE
2515 C
2516 C CONSTANTS
2517 C
2518  INTEGER ngraphs
2519  parameter(ngraphs=43)
2520  INTEGER nexternal
2521  parameter(nexternal=6)
2522  INTEGER nwavefuncs, ncolor
2523  parameter(nwavefuncs=14, ncolor=2)
2524  REAL*8 zero
2525  parameter(zero=0d0)
2526  COMPLEX*16 imag1
2527  parameter(imag1=(0d0,1d0))
2528 C
2529 C ARGUMENTS
2530 C
2531  REAL*8 p(0:3,nexternal)
2532  INTEGER nhel(nexternal), ic(nexternal)
2533 C
2534 C LOCAL VARIABLES
2535 C
2536  INTEGER i,j
2537  COMPLEX*16 ztemp
2538  REAL*8 denom(ncolor), cf(ncolor,ncolor)
2539  COMPLEX*16 amp(ngraphs), jamp(ncolor)
2540  COMPLEX*16 w(18,nwavefuncs)
2541  COMPLEX*16 dum0,dum1
2542  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2543 C
2544 C GLOBAL VARIABLES
2545 C
2546  include 'coupl.inc'
2547 C
2548 C COLOR DATA
2549 C
2550  DATA denom(1)/1/
2551  DATA (cf(i, 1),i= 1, 2) / 9, 3/
2552 C 1 T(2,1) T(3,4)
2553  DATA denom(2)/1/
2554  DATA (cf(i, 2),i= 1, 2) / 3, 9/
2555 C 1 T(2,4) T(3,1)
2556 
2557  INTEGER h1,h2
2558  REAL*8 matrix
2559  matrix_ddx_ccx_noh=0.d0
2560  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2561  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2562 
2563 
2564 C ----------
2565 C BEGIN CODE
2566 C ----------
2567  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2568  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2569  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2570  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2571  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2572  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2573  CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
2574  CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
2575  CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2576 C Amplitude(s) for diagram number 1
2577  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2578  CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2579 C Amplitude(s) for diagram number 2
2580  CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2581  CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
2582 C Amplitude(s) for diagram number 3
2583  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2584 C Amplitude(s) for diagram number 4
2585  CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2586  CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
2587  CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2588 C Amplitude(s) for diagram number 5
2589  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2590  CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2591 C Amplitude(s) for diagram number 6
2592  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2593 C Amplitude(s) for diagram number 7
2594  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2595 C Amplitude(s) for diagram number 8
2596  CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2597  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2598  CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
2599 C Amplitude(s) for diagram number 9
2600  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
2601  CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
2602 C Amplitude(s) for diagram number 10
2603  CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
2604  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2605 C Amplitude(s) for diagram number 11
2606  CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
2607 C Amplitude(s) for diagram number 12
2608  CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
2609  CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2610  CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2611 C Amplitude(s) for diagram number 13
2612  CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
2613  CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
2614 C Amplitude(s) for diagram number 14
2615  CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
2616 C Amplitude(s) for diagram number 15
2617  CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
2618 C Amplitude(s) for diagram number 16
2619  CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
2620  CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
2621 C Amplitude(s) for diagram number 17
2622  CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
2623  CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
2624 C Amplitude(s) for diagram number 18
2625  CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
2626 C Amplitude(s) for diagram number 19
2627  CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
2628 C Amplitude(s) for diagram number 20
2629  CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
2630  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,9))
2631  CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,14))
2632  CALL ffv2_1(w(1,6),w(1,9),gc_108,zero,zero,w(1,10))
2633 C Amplitude(s) for diagram number 21
2634  CALL ffv2_0(w(1,5),w(1,10),w(1,14),gc_108,amp(21))
2635 C Amplitude(s) for diagram number 22
2636  CALL vvv1_0(w(1,12),w(1,14),w(1,9),gc_4,amp(22))
2637 C Amplitude(s) for diagram number 23
2638  CALL vvv1_0(w(1,14),w(1,9),w(1,7),gc_53,amp(23))
2639  CALL ffv2_1(w(1,2),w(1,9),gc_44,zero,zero,w(1,10))
2640 C Amplitude(s) for diagram number 24
2641  CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_2,amp(24))
2642  CALL ffv2_2(w(1,4),w(1,9),gc_44,zero,zero,w(1,5))
2643 C Amplitude(s) for diagram number 25
2644  CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(25))
2645 C Amplitude(s) for diagram number 26
2646  CALL ffv2_5_0(w(1,4),w(1,10),w(1,7),gc_51,gc_58,amp(26))
2647 C Amplitude(s) for diagram number 27
2648  CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2649  CALL ffv2_2(w(1,1),w(1,14),gc_44,zero,zero,w(1,5))
2650 C Amplitude(s) for diagram number 28
2651  CALL ffv1_0(w(1,5),w(1,3),w(1,12),gc_2,amp(28))
2652  CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,10))
2653 C Amplitude(s) for diagram number 29
2654  CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_44,amp(29))
2655 C Amplitude(s) for diagram number 30
2656  CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(30))
2657  CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,5))
2658 C Amplitude(s) for diagram number 31
2659  CALL ffv2_0(w(1,5),w(1,3),w(1,14),gc_44,amp(31))
2660  CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,14))
2661 C Amplitude(s) for diagram number 32
2662  CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_1,amp(32))
2663 C Amplitude(s) for diagram number 33
2664  CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(33))
2665 C Amplitude(s) for diagram number 34
2666  CALL ffv2_3_0(w(1,14),w(1,2),w(1,7),gc_50,gc_58,amp(34))
2667 C Amplitude(s) for diagram number 35
2668  CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(35))
2669  CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2670  CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2671 C Amplitude(s) for diagram number 36
2672  CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(36))
2673 C Amplitude(s) for diagram number 37
2674  CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
2675 C Amplitude(s) for diagram number 38
2676  CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(38))
2677 C Amplitude(s) for diagram number 39
2678  CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(39))
2679  CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2680 C Amplitude(s) for diagram number 40
2681  CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(40))
2682 C Amplitude(s) for diagram number 41
2683  CALL ffv2_3_0(w(1,10),w(1,2),w(1,11),gc_50,gc_58,amp(41))
2684 C Amplitude(s) for diagram number 42
2685  CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(42))
2686 C Amplitude(s) for diagram number 43
2687  CALL ffv2_3_0(w(1,5),w(1,2),w(1,11),gc_50,gc_58,amp(43))
2688  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2689  $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2690  $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2691  $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
2692  $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
2693  jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2694  $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2695  $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
2696  $ -1./2.*amp(38)-1./2.*amp(39)
2697 
2698  matrix = 0.d0
2699  DO i = 1, ncolor
2700  ztemp = (0.d0,0.d0)
2701  DO j = 1, ncolor
2702  ztemp = ztemp + cf(j,i)*jamp(j)
2703  ENDDO
2704  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2705  ENDDO
2706  matrix_ddx_ccx_noh=matrix
2707  ENDIF
2708  ENDIF
2709  END
2710 
2711 C -----------------SUBPROCESS DDX->UCX, NO HIGGS
2712 
2713  SUBROUTINE ddx_ucx_noh(P,H1,H2,ANS)
2714 C
2715 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2716 C By the MadGraph Development Team
2717 C Please visit us at https://launchpad.net/madgraph5
2718 C
2719 C MadGraph StandAlone Version
2720 C
2721 C Returns amplitude squared summed/avg over colors
2722 C and helicities
2723 C for the point in phase space P(0:3,NEXTERNAL)
2724 C
2725 C Process: d d~ > u c~ ta+ ta- / h QED=4
2726 C
2727  IMPLICIT NONE
2728 C
2729 C CONSTANTS
2730 C
2731  INTEGER nexternal
2732  parameter(nexternal=6)
2733  INTEGER ncomb
2734  parameter( ncomb=64)
2735 C
2736 C ARGUMENTS
2737 C
2738  REAL*8 p(0:3,nexternal),ans
2739  INTEGER h1,h2
2740 C
2741 C LOCAL VARIABLES
2742 C
2743  INTEGER nhel(nexternal,ncomb),ntry
2744  REAL*8 t
2745  REAL*8 matrix_ddx_ucx_noh
2746  INTEGER ihel,iden, i
2747  INTEGER jc(nexternal)
2748  LOGICAL goodhel(ncomb)
2749  DATA ntry/0/
2750  DATA goodhel/ncomb*.false./
2751  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2752  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2753  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2754  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2755  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2756  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2757  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2758  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2759  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2760  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2761  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2762  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2763  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2764  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2765  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2766  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2767  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2768  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2769  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2770  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2771  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2772  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2773  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2774  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2775  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2776  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2777  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2778  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2779  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2780  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2781  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2782  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2783  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2784  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2785  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2786  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2787  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2788  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2789  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2790  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2791  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2792  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2793  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2794  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2795  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2796  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2797  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2798  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2799  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2800  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2801  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2802  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2803  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2804  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2805  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2806  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2807  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2808  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2809  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2810  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2811  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2812  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2813  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2814  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2815  DATA iden/36/
2816 C ----------
2817 C BEGIN CODE
2818 C ----------
2819  DO ihel=1,nexternal
2820  jc(ihel) = +1
2821  ENDDO
2822  ans = 0d0
2823  DO ihel=1,ncomb
2824  t=matrix_ddx_ucx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2825  ans=ans+t
2826  ENDDO
2827  ans=ans/dble(iden)
2828  END
2829 
2830 
2831  REAL*8 FUNCTION matrix_ddx_ucx_noh(P,H1,H2,NHEL,IC)
2832 C
2833 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2834 C By the MadGraph Development Team
2835 C Please visit us at https://launchpad.net/madgraph5
2836 C
2837 C Returns amplitude squared summed/avg over colors
2838 C for the point with external lines W(0:6,NEXTERNAL)
2839 C
2840 C Process: d d~ > u c~ ta+ ta- / h QED=4
2841 C
2842  IMPLICIT NONE
2843 C
2844 C CONSTANTS
2845 C
2846  INTEGER ngraphs
2847  parameter(ngraphs=11)
2848  INTEGER nexternal
2849  parameter(nexternal=6)
2850  INTEGER nwavefuncs, ncolor
2851  parameter(nwavefuncs=10, ncolor=1)
2852  REAL*8 zero
2853  parameter(zero=0d0)
2854  COMPLEX*16 imag1
2855  parameter(imag1=(0d0,1d0))
2856 C
2857 C ARGUMENTS
2858 C
2859  REAL*8 p(0:3,nexternal)
2860  INTEGER nhel(nexternal), ic(nexternal)
2861 C
2862 C LOCAL VARIABLES
2863 C
2864  INTEGER i,j
2865  COMPLEX*16 ztemp
2866  REAL*8 denom(ncolor), cf(ncolor,ncolor)
2867  COMPLEX*16 amp(ngraphs), jamp(ncolor)
2868  COMPLEX*16 w(18,nwavefuncs)
2869  COMPLEX*16 dum0,dum1
2870  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2871 C
2872 C GLOBAL VARIABLES
2873 C
2874  include 'coupl.inc'
2875 C
2876 C COLOR DATA
2877 C
2878  DATA denom(1)/1/
2879  DATA (cf(i, 1),i= 1, 1) / 9/
2880 C 1 T(2,4) T(3,1)
2881 
2882  INTEGER h1,h2
2883  REAL*8 matrix
2884  matrix_ddx_ucx_noh=0.d0
2885  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2886  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2887 
2888 C ----------
2889 C BEGIN CODE
2890 C ----------
2891  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2892  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2893  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2894  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2895  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2896  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2897  CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
2898  CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,8))
2899  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2900 C Amplitude(s) for diagram number 1
2901  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2902  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2903 C Amplitude(s) for diagram number 2
2904  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2905  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2906 C Amplitude(s) for diagram number 3
2907  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2908  CALL ffv2_1(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
2909 C Amplitude(s) for diagram number 4
2910  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2911  CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
2912 C Amplitude(s) for diagram number 5
2913  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2914 C Amplitude(s) for diagram number 6
2915  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2916 C Amplitude(s) for diagram number 7
2917  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2918  CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
2919 C Amplitude(s) for diagram number 8
2920  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2921  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2922 C Amplitude(s) for diagram number 9
2923  CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
2924 C Amplitude(s) for diagram number 10
2925  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2926  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2927 C Amplitude(s) for diagram number 11
2928  CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
2929  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2930  $ +amp(9)+amp(10)+amp(11)
2931 
2932  matrix = 0.d0
2933  DO i = 1, ncolor
2934  ztemp = (0.d0,0.d0)
2935  DO j = 1, ncolor
2936  ztemp = ztemp + cf(j,i)*jamp(j)
2937  ENDDO
2938  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2939  ENDDO
2940  matrix_ddx_ucx_noh=matrix
2941  ENDIF
2942  ENDIF
2943  END
2944 
2945 C ----------------SUBPROCESS DDX->CUX, NO HIGGS
2946 
2947  SUBROUTINE ddx_cux_noh(P,H1,H2,ANS)
2948 C
2949 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2950 C By the MadGraph Development Team
2951 C Please visit us at https://launchpad.net/madgraph5
2952 C
2953 C MadGraph StandAlone Version
2954 C
2955 C Returns amplitude squared summed/avg over colors
2956 C and helicities
2957 C for the point in phase space P(0:3,NEXTERNAL)
2958 C
2959 C Process: d d~ > c u~ ta+ ta- / h QED=4
2960 C
2961  IMPLICIT NONE
2962 C
2963 C CONSTANTS
2964 C
2965  INTEGER nexternal
2966  parameter(nexternal=6)
2967  INTEGER ncomb
2968  parameter( ncomb=64)
2969 C
2970 C ARGUMENTS
2971 C
2972  REAL*8 p(0:3,nexternal),ans
2973  INTEGER h1,h2
2974 C
2975 C LOCAL VARIABLES
2976 C
2977  INTEGER nhel(nexternal,ncomb),ntry
2978  REAL*8 t
2979  REAL*8 matrix_ddx_cux_noh
2980  INTEGER ihel,iden, i
2981  INTEGER jc(nexternal)
2982  LOGICAL goodhel(ncomb)
2983  DATA ntry/0/
2984  DATA goodhel/ncomb*.false./
2985  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2986  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2987  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2988  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2989  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2990  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2991  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2992  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2993  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2994  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2995  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2996  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2997  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2998  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2999  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3000  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3001  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3002  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3003  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3004  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3005  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3006  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3007  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3008  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3009  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3010  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3011  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3012  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3013  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3014  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3015  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3016  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3017  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3018  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3019  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3020  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3021  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3022  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3023  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3024  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3025  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3026  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3027  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3028  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3029  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3030  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3031  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3032  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3033  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3034  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3035  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3036  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3037  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3038  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3039  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3040  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3041  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3042  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3043  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3044  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3045  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3046  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3047  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3048  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3049  DATA iden/36/
3050 C ----------
3051 C BEGIN CODE
3052 C ----------
3053  DO ihel=1,nexternal
3054  jc(ihel) = +1
3055  ENDDO
3056  ans = 0d0
3057  DO ihel=1,ncomb
3058  t=matrix_ddx_cux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3059  ans=ans+t
3060  ENDDO
3061  ans=ans/dble(iden)
3062  END
3063 
3064 
3065  REAL*8 FUNCTION matrix_ddx_cux_noh(P,H1,H2,NHEL,IC)
3066 C
3067 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
3068 C By the MadGraph Development Team
3069 C Please visit us at https://launchpad.net/madgraph5
3070 C
3071 C Returns amplitude squared summed/avg over colors
3072 C for the point with external lines W(0:6,NEXTERNAL)
3073 C
3074 C Process: d d~ > c u~ ta+ ta- / h QED=4
3075 C
3076  IMPLICIT NONE
3077 C
3078 C CONSTANTS
3079 C
3080  INTEGER ngraphs
3081  parameter(ngraphs=11)
3082  INTEGER nexternal
3083  parameter(nexternal=6)
3084  INTEGER nwavefuncs, ncolor
3085  parameter(nwavefuncs=10, ncolor=1)
3086  REAL*8 zero
3087  parameter(zero=0d0)
3088  COMPLEX*16 imag1
3089  parameter(imag1=(0d0,1d0))
3090 C
3091 C ARGUMENTS
3092 C
3093  REAL*8 p(0:3,nexternal)
3094  INTEGER nhel(nexternal), ic(nexternal)
3095 C
3096 C LOCAL VARIABLES
3097 C
3098  INTEGER i,j
3099  COMPLEX*16 ztemp
3100  REAL*8 denom(ncolor), cf(ncolor,ncolor)
3101  COMPLEX*16 amp(ngraphs), jamp(ncolor)
3102  COMPLEX*16 w(18,nwavefuncs)
3103  COMPLEX*16 dum0,dum1
3104  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3105 C
3106 C GLOBAL VARIABLES
3107 C
3108  include 'coupl.inc'
3109 C
3110 C COLOR DATA
3111 C
3112  DATA denom(1)/1/
3113  DATA (cf(i, 1),i= 1, 1) / 9/
3114 C 1 T(2,4) T(3,1)
3115 
3116  INTEGER h1,h2
3117  REAL*8 matrix
3118  matrix_ddx_cux_noh=0.d0
3119  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
3120  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
3121 
3122 C ----------
3123 C BEGIN CODE
3124 C ----------
3125  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3126  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3127  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3128  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
3129  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3130  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3131  CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,7))
3132  CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
3133  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
3134 C Amplitude(s) for diagram number 1
3135  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
3136  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
3137 C Amplitude(s) for diagram number 2
3138  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
3139  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
3140 C Amplitude(s) for diagram number 3
3141  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
3142  CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
3143 C Amplitude(s) for diagram number 4
3144  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
3145  CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
3146 C Amplitude(s) for diagram number 5
3147  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
3148 C Amplitude(s) for diagram number 6
3149  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
3150 C Amplitude(s) for diagram number 7
3151  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
3152  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
3153 C Amplitude(s) for diagram number 8
3154  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
3155  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
3156 C Amplitude(s) for diagram number 9
3157  CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_44,amp(9))
3158 C Amplitude(s) for diagram number 10
3159  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
3160  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
3161 C Amplitude(s) for diagram number 11
3162  CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_44,amp(11))
3163  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
3164  $ +amp(9)+amp(10)+amp(11)
3165 
3166  matrix = 0.d0
3167  DO i = 1, ncolor
3168  ztemp = (0.d0,0.d0)
3169  DO j = 1, ncolor
3170  ztemp = ztemp + cf(j,i)*jamp(j)
3171  ENDDO
3172  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3173  ENDDO
3174  matrix_ddx_cux_noh=matrix
3175  ENDIF
3176  ENDIF
3177  END
3178 
3179 C DONE 19.11.2014
3180 
3181 C ------------------------------------------------------
3182 C WITH GLUONS AS JETS ADDED 9.07.2015
3183 
3184 C --- BEGIN SUBROUTINE DDX->GG NO HIGGS
3185 
3186 
3187  SUBROUTINE ddx_gg_noh(P,H1,H2,ANS)
3188 C
3189 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
3190 C By the MadGraph Development Team
3191 C Please visit us at https://launchpad.net/madgraph5
3192 C
3193 C MadGraph StandAlone Version
3194 C
3195 C Returns amplitude squared summed/avg over colors
3196 C and helicities
3197 C for the point in phase space P(0:3,NEXTERNAL)
3198 C
3199 C Process: d d~ > g g ta+ ta- / h QED=4
3200 C Process: s s~ > g g ta+ ta- / h QED=4
3201 C
3202  IMPLICIT NONE
3203 C
3204 C CONSTANTS
3205 C
3206  INTEGER nexternal
3207  parameter(nexternal=6)
3208  INTEGER ncomb
3209  parameter( ncomb=64)
3210 C
3211 C ARGUMENTS
3212 C
3213  REAL*8 p(0:3,nexternal),ans
3214  INTEGER h1,h2
3215 C
3216 C LOCAL VARIABLES
3217 C
3218  INTEGER nhel(nexternal,ncomb),ntry
3219  REAL*8 t
3220  REAL*8 matrix_ddx_gg_noh
3221  INTEGER ihel,iden, i
3222  INTEGER jc(nexternal)
3223  LOGICAL goodhel(ncomb)
3224  DATA ntry/0/
3225  DATA goodhel/ncomb*.false./
3226  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
3227  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
3228  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
3229  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
3230  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
3231  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
3232  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
3233  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
3234  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
3235  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
3236  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
3237  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
3238  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
3239  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
3240  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3241  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3242  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3243  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3244  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3245  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3246  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3247  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3248  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3249  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3250  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3251  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3252  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3253  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3254  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3255  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3256  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3257  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3258  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3259  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3260  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3261  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3262  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3263  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3264  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3265  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3266  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3267  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3268  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3269  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3270  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3271  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3272  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3273  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3274  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3275  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3276  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3277  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3278  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3279  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3280  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3281  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3282  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3283  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3284  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3285  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3286  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3287  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3288  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3289  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3290  DATA iden/72/
3291 C ----------
3292 C BEGIN CODE
3293 C ----------
3294  DO ihel=1,nexternal
3295  jc(ihel) = +1
3296  ENDDO
3297  ans = 0d0
3298  DO ihel=1,ncomb
3299  t=matrix_ddx_gg_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3300  ans=ans+t
3301  ENDDO
3302  ans=ans/dble(iden)
3303  END
3304 
3305 
3306  REAL*8 FUNCTION matrix_ddx_gg_noh(P,H1,H2,NHEL,IC)
3307 C
3308 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
3309 C By the MadGraph Development Team
3310 C Please visit us at https://launchpad.net/madgraph5
3311 C
3312 C Returns amplitude squared summed/avg over colors
3313 C for the point with external lines W(0:6,NEXTERNAL)
3314 C
3315 C Process: d d~ > g g ta+ ta- / h QED=4
3316 C Process: s s~ > g g ta+ ta- / h QED=4
3317 C
3318  IMPLICIT NONE
3319 C
3320 C CONSTANTS
3321 C
3322  INTEGER ngraphs
3323  parameter(ngraphs=16)
3324  INTEGER nexternal
3325  parameter(nexternal=6)
3326  INTEGER nwavefuncs, ncolor
3327  parameter(nwavefuncs=10, ncolor=2)
3328  REAL*8 zero
3329  parameter(zero=0d0)
3330  COMPLEX*16 imag1
3331  parameter(imag1=(0d0,1d0))
3332 C
3333 C ARGUMENTS
3334 C
3335  REAL*8 p(0:3,nexternal)
3336  INTEGER nhel(nexternal), ic(nexternal)
3337 C
3338 C LOCAL VARIABLES
3339 C
3340  INTEGER i,j
3341  COMPLEX*16 ztemp
3342  REAL*8 denom(ncolor), cf(ncolor,ncolor)
3343  COMPLEX*16 amp(ngraphs), jamp(ncolor)
3344  COMPLEX*16 w(18,nwavefuncs)
3345  COMPLEX*16 dum0,dum1
3346  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3347 C
3348 C GLOBAL VARIABLES
3349 C
3350  include 'coupl.inc'
3351 C
3352 C COLOR DATA
3353 C
3354  DATA denom(1)/3/
3355  DATA (cf(i, 1),i= 1, 2) / 16, -2/
3356 C 1 T(3,4,2,1)
3357  DATA denom(2)/3/
3358  DATA (cf(i, 2),i= 1, 2) / -2, 16/
3359 C 1 T(4,3,2,1)
3360 
3361  INTEGER h1,h2
3362  REAL*8 matrix
3363  matrix_ddx_gg_noh=0.d0
3364  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
3365  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
3366 
3367 C ----------
3368 C BEGIN CODE
3369 C ----------
3370  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3371  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3372  CALL vxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3373  CALL vxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
3374  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3375  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3376  CALL vvv1p0_1(w(1,3),w(1,4),gc_10,zero,zero,w(1,7))
3377  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,8))
3378  CALL ffv1_1(w(1,2),w(1,7),gc_11,zero,zero,w(1,9))
3379 C Amplitude(s) for diagram number 1
3380  CALL ffv1_0(w(1,1),w(1,9),w(1,8),gc_1,amp(1))
3381  CALL ffv1_2(w(1,1),w(1,7),gc_11,zero,zero,w(1,10))
3382 C Amplitude(s) for diagram number 2
3383  CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(2))
3384  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
3385 C Amplitude(s) for diagram number 3
3386  CALL ffv2_3_0(w(1,1),w(1,9),w(1,7),gc_50,gc_58,amp(3))
3387 C Amplitude(s) for diagram number 4
3388  CALL ffv2_3_0(w(1,10),w(1,2),w(1,7),gc_50,gc_58,amp(4))
3389  CALL ffv1_1(w(1,2),w(1,3),gc_11,zero,zero,w(1,10))
3390  CALL ffv1_2(w(1,1),w(1,4),gc_11,zero,zero,w(1,9))
3391 C Amplitude(s) for diagram number 5
3392  CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_1,amp(5))
3393 C Amplitude(s) for diagram number 6
3394  CALL ffv2_3_0(w(1,9),w(1,10),w(1,7),gc_50,gc_58,amp(6))
3395  CALL ffv1_1(w(1,10),w(1,4),gc_11,zero,zero,w(1,6))
3396 C Amplitude(s) for diagram number 7
3397  CALL ffv1_0(w(1,1),w(1,6),w(1,8),gc_1,amp(7))
3398 C Amplitude(s) for diagram number 8
3399  CALL ffv2_3_0(w(1,1),w(1,6),w(1,7),gc_50,gc_58,amp(8))
3400  CALL ffv1_2(w(1,1),w(1,3),gc_11,zero,zero,w(1,6))
3401  CALL ffv1_1(w(1,2),w(1,4),gc_11,zero,zero,w(1,10))
3402 C Amplitude(s) for diagram number 9
3403  CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_1,amp(9))
3404 C Amplitude(s) for diagram number 10
3405  CALL ffv2_3_0(w(1,6),w(1,10),w(1,7),gc_50,gc_58,amp(10))
3406  CALL ffv1_2(w(1,6),w(1,4),gc_11,zero,zero,w(1,5))
3407 C Amplitude(s) for diagram number 11
3408  CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(11))
3409 C Amplitude(s) for diagram number 12
3410  CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(12))
3411  CALL ffv1_1(w(1,10),w(1,3),gc_11,zero,zero,w(1,5))
3412 C Amplitude(s) for diagram number 13
3413  CALL ffv1_0(w(1,1),w(1,5),w(1,8),gc_1,amp(13))
3414 C Amplitude(s) for diagram number 14
3415  CALL ffv2_3_0(w(1,1),w(1,5),w(1,7),gc_50,gc_58,amp(14))
3416  CALL ffv1_2(w(1,9),w(1,3),gc_11,zero,zero,w(1,5))
3417 C Amplitude(s) for diagram number 15
3418  CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(15))
3419 C Amplitude(s) for diagram number 16
3420  CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(16))
3421  jamp(1)=-imag1*amp(1)-imag1*amp(2)-imag1*amp(3)-imag1*amp(4)
3422  $ +amp(5)+amp(6)+amp(7)+amp(8)+amp(15)+amp(16)
3423  jamp(2)=+imag1*amp(1)+imag1*amp(2)+imag1*amp(3)+imag1*amp(4)
3424  $ +amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
3425 
3426  matrix = 0.d0
3427  DO i = 1, ncolor
3428  ztemp = (0.d0,0.d0)
3429  DO j = 1, ncolor
3430  ztemp = ztemp + cf(j,i)*jamp(j)
3431  ENDDO
3432  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3433  ENDDO
3434  matrix_ddx_gg_noh=matrix
3435  ENDIF
3436  ENDIF
3437  END
3438 
3439 
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7