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