C++ Interface to Tauola
USX.f
1 c --------------------- subprocesses initiated by USbar
2 
3  SUBROUTINE usx(P,I3,I4,H1,H2,KEY,ANS)
4  IMPLICIT NONE
5 
6  INTEGER I3,I4, H1,H2, KEY
7  real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
8 
9 C I3=0,2 OR 4 I4=0,-1 OR -3
10 
11  ans=0.d0
12 
13 C FINAL CSX IS EQUAL TO FINAL UDX
14 
15  IF(key.EQ.1) THEN
16  IF(i3.EQ.2 .AND. i4.EQ.-1) CALL usx_udx_h(p,h1,h2,ans)
17  IF(i3.EQ.2 .AND. i4.EQ.-3) CALL usx_usx_h(p,h1,h2,ans)
18  IF(i3.EQ.4 .AND. i4.EQ.-1) CAll usx_cdx_h(p,h1,h2,ans)
19  IF(i3.EQ.4 .AND. i4.EQ.-3) CALL usx_udx_h(p,h1,h2,ans)
20  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
21  CALL usx_udx_h(p,h1,h2,ans1)
22  CALL usx_usx_h(p,h1,h2,ans2)
23  CALL usx_cdx_h(p,h1,h2,ans3)
24  CALL usx_udx_h(p,h1,h2,ans4)
25  ans=ans1+ans2+ans3+ans4
26  ENDIF
27  ELSE IF(key.EQ.0) THEN
28  IF(i3.EQ.2 .AND. i4.EQ.-1) CALL usx_udx_noh(p,h1,h2,ans)
29  IF(i3.EQ.2 .AND. i4.EQ.-3) CALL usx_usx_noh(p,h1,h2,ans)
30  IF(i3.EQ.4 .AND. i4.EQ.-1) CAll usx_cdx_noh(p,h1,h2,ans)
31  IF(i3.EQ.4 .AND. i4.EQ.-3) CALL usx_udx_noh(p,h1,h2,ans)
32  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
33  CALL usx_udx_noh(p,h1,h2,ans1)
34  CALL usx_usx_noh(p,h1,h2,ans2)
35  CALL usx_cdx_noh(p,h1,h2,ans3)
36  CALL usx_udx_noh(p,h1,h2,ans4)
37  ans=ans1+ans2+ans3+ans4
38  ENDIF
39  ELSE
40  WRITE(*,*) 'NOT FINISHED'
41  stop
42  ENDIF
43  END ! SUBROUTINE USX
44 
45 
46 
47 C ----- begin subprocesses USX->UDX with Higgs->tautau
48 
49 
50  SUBROUTINE usx_udx_h(P,H1,H2,ANS)
51 C
52 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
53 C By the MadGraph Development Team
54 C Please visit us at https://launchpad.net/madgraph5
55 C
56 C MadGraph StandAlone Version
57 C
58 C Returns amplitude squared summed/avg over colors
59 C and helicities
60 C for the point in phase space P(0:3,NEXTERNAL)
61 C
62 C Process: u s~ > u d~ h WEIGHTED=6
63 C * Decay: h > ta+ ta- WEIGHTED=2
64 C Process: u s~ > c s~ h WEIGHTED=6
65 C * Decay: h > ta+ ta- WEIGHTED=2
66 C
67  IMPLICIT NONE
68 C
69 C CONSTANTS
70 C
71  INTEGER NEXTERNAL
72  parameter(nexternal=6)
73  INTEGER NCOMB
74  parameter( ncomb=64)
75 C
76 C ARGUMENTS
77 C
78  real*8 p(0:3,nexternal),ans
79  INTEGER H1,H2
80 C
81 C LOCAL VARIABLES
82 C
83  INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
84  real*8 t
85  real*8 matrix_usx_udx_h
86  INTEGER IHEL,IDEN, I
87  INTEGER JC(NEXTERNAL)
88  LOGICAL GOODHEL(NCOMB)
89  DATA ntry/0/
90  DATA goodhel/ncomb*.false./
91  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
92  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
93  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
94  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
95  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
96  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
97  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
98  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
99  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
100  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
101  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
102  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
103  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
104  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
105  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
106  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
107  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
108  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
109  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
110  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
111  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
112  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
113  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
114  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
115  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
116  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
117  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
118  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
119  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
120  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
121  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
122  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
123  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
124  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
125  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
126  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
127  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
128  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
129  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
130  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
131  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
132  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
133  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
134  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
135  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
136  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
137  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
138  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
139  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
140  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
141  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
142  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
143  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
144  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
145  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
146  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
147  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
148  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
149  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
150  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
151  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
152  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
153  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
154  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
155  DATA iden/36/
156 C ----------
157 C BEGIN CODE
158 C ----------
159  DO ihel=1,nexternal
160  jc(ihel) = +1
161  ENDDO
162  ans = 0d0
163  DO ihel=1,ncomb
164  t=matrix_usx_udx_h(p ,h1,h2,nhel(1,ihel),jc(1))
165  ans=ans+t
166  ENDDO
167  ans=ans/dble(iden)
168  END
169 
170 
171  REAL*8 FUNCTION matrix_usx_udx_h(P,H1,H2,NHEL,IC)
172 C
173 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
174 C By the MadGraph Development Team
175 C Please visit us at https://launchpad.net/madgraph5
176 C
177 C Returns amplitude squared summed/avg over colors
178 C for the point with external lines W(0:6,NEXTERNAL)
179 C
180 C Process: u s~ > u d~ h WEIGHTED=6
181 C * Decay: h > ta+ ta- WEIGHTED=2
182 C Process: u s~ > c s~ h WEIGHTED=6
183 C * Decay: h > ta+ ta- WEIGHTED=2
184 C
185  IMPLICIT NONE
186 C
187 C CONSTANTS
188 C
189  INTEGER NGRAPHS
190  parameter(ngraphs=1)
191  INTEGER NEXTERNAL
192  parameter(nexternal=6)
193  INTEGER NWAVEFUNCS, NCOLOR
194  parameter(nwavefuncs=7, ncolor=1)
195  real*8 zero
196  parameter(zero=0d0)
197  COMPLEX*16 IMAG1
198  parameter(imag1=(0d0,1d0))
199 C
200 C ARGUMENTS
201 C
202  real*8 p(0:3,nexternal)
203  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
204 C
205 C LOCAL VARIABLES
206 C
207  INTEGER I,J
208  COMPLEX*16 ZTEMP
209  real*8 denom(ncolor), cf(ncolor,ncolor)
210  COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
211  COMPLEX*16 W(18,NWAVEFUNCS)
212  COMPLEX*16 DUM0,DUM1
213  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
214 C
215 C GLOBAL VARIABLES
216 C
217  include 'coupl.inc'
218 C
219 C COLOR DATA
220 C
221  DATA denom(1)/1/
222  DATA (cf(i, 1),i= 1, 1) / 9/
223 C 1 T(2,1) T(3,4)
224 
225  INTEGER H1,H2
226  real*8 matrix
227  matrix_usx_udx_h=0.d0
228  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
229  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
230 
231 
232 C ----------
233 C BEGIN CODE
234 C ----------
235  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
236  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
237  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
238  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
239  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
240  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
241  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
242  CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
243  CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
244 C Amplitude(s) for diagram number 1
245  CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
246  jamp(1)=-amp(1)
247 
248  matrix = 0.d0
249  DO i = 1, ncolor
250  ztemp = (0.d0,0.d0)
251  DO j = 1, ncolor
252  ztemp = ztemp + cf(j,i)*jamp(j)
253  ENDDO
254  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
255  ENDDO
256 
257  matrix_usx_udx_h=matrix
258  ENDIF
259  ENDIF
260  END
261 
262 C ----- begin subprocesses USX->USX with Higgs->tautau
263 
264  SUBROUTINE usx_usx_h(P,H1,H2,ANS)
265 C
266 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
267 C By the MadGraph Development Team
268 C Please visit us at https://launchpad.net/madgraph5
269 C
270 C MadGraph StandAlone Version
271 C
272 C Returns amplitude squared summed/avg over colors
273 C and helicities
274 C for the point in phase space P(0:3,NEXTERNAL)
275 C
276 C Process: u s~ > u s~ h WEIGHTED=6
277 C * Decay: h > ta+ ta- WEIGHTED=2
278 C
279  IMPLICIT NONE
280 C
281 C CONSTANTS
282 C
283  INTEGER NEXTERNAL
284  parameter(nexternal=6)
285  INTEGER NCOMB
286  parameter( ncomb=64)
287 C
288 C ARGUMENTS
289 C
290  real*8 p(0:3,nexternal),ans
291  INTEGER H1,H2
292 C
293 C LOCAL VARIABLES
294 C
295  INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
296  real*8 t
297  real*8 matrix_usx_usx_h
298  INTEGER IHEL,IDEN, I
299  INTEGER JC(NEXTERNAL)
300  LOGICAL GOODHEL(NCOMB)
301  DATA ntry/0/
302  DATA goodhel/ncomb*.false./
303  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
304  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
305  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
306  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
307  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
308  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
309  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
310  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
311  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
312  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
313  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
314  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
315  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
316  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
317  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
318  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
319  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
320  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
321  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
322  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
323  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
324  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
325  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
326  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
327  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
328  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
329  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
330  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
331  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
332  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
333  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
334  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
335  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
336  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
337  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
338  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
339  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
340  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
341  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
342  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
343  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
344  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
345  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
346  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
347  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
348  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
349  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
350  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
351  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
352  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
353  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
354  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
355  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
356  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
357  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
358  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
359  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
360  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
361  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
362  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
363  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
364  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
365  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
366  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
367  DATA iden/36/
368 C ----------
369 C BEGIN CODE
370 C ----------
371  DO ihel=1,nexternal
372  jc(ihel) = +1
373  ENDDO
374  ans = 0d0
375  DO ihel=1,ncomb
376  t=matrix_usx_usx_h(p ,h1,h2,nhel(1,ihel),jc(1))
377  ans=ans+t
378  ENDDO
379  ans=ans/dble(iden)
380  END
381 
382 
383  REAL*8 FUNCTION matrix_usx_usx_h(P,H1,H2,NHEL,IC)
384 C
385 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
386 C By the MadGraph Development Team
387 C Please visit us at https://launchpad.net/madgraph5
388 C
389 C Returns amplitude squared summed/avg over colors
390 C for the point with external lines W(0:6,NEXTERNAL)
391 C
392 C Process: u s~ > u s~ h WEIGHTED=6
393 C * Decay: h > ta+ ta- WEIGHTED=2
394 C
395  IMPLICIT NONE
396 C
397 C CONSTANTS
398 C
399  INTEGER NGRAPHS
400  parameter(ngraphs=2)
401  INTEGER NEXTERNAL
402  parameter(nexternal=6)
403  INTEGER NWAVEFUNCS, NCOLOR
404  parameter(nwavefuncs=7, ncolor=2)
405  real*8 zero
406  parameter(zero=0d0)
407  COMPLEX*16 IMAG1
408  parameter(imag1=(0d0,1d0))
409 C
410 C ARGUMENTS
411 C
412  real*8 p(0:3,nexternal)
413  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
414 C
415 C LOCAL VARIABLES
416 C
417  INTEGER I,J
418  COMPLEX*16 ZTEMP
419  real*8 denom(ncolor), cf(ncolor,ncolor)
420  COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
421  COMPLEX*16 W(18,NWAVEFUNCS)
422  COMPLEX*16 DUM0,DUM1
423  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
424 C
425 C GLOBAL VARIABLES
426 C
427  include 'coupl.inc'
428 C
429 C COLOR DATA
430 C
431  DATA denom(1)/1/
432  DATA (cf(i, 1),i= 1, 2) / 9, 3/
433 C 1 T(2,1) T(3,4)
434  DATA denom(2)/1/
435  DATA (cf(i, 2),i= 1, 2) / 3, 9/
436 C 1 T(2,4) T(3,1)
437 
438  INTEGER H1,H2
439  real*8 matrix
440  matrix_usx_usx_h=0.d0
441  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
442  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
443 
444 
445 C ----------
446 C BEGIN CODE
447 C ----------
448  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
449  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
450  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
451  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
452  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
453  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
454  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
455  CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
456  CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,5))
457 C Amplitude(s) for diagram number 1
458  CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_72,amp(1))
459  CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
460  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
461 C Amplitude(s) for diagram number 2
462  CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
463  jamp(1)=-amp(1)
464  jamp(2)=+amp(2)
465 
466  matrix = 0.d0
467  DO i = 1, ncolor
468  ztemp = (0.d0,0.d0)
469  DO j = 1, ncolor
470  ztemp = ztemp + cf(j,i)*jamp(j)
471  ENDDO
472  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
473  ENDDO
474  matrix_usx_usx_h=matrix
475  ENDIF
476  ENDIF
477  END
478 
479 C ----- begin subprocesses USX->CDX with Higgs->tautau
480 
481  SUBROUTINE usx_cdx_h(P,H1,H2,ANS)
482 C
483 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
484 C By the MadGraph Development Team
485 C Please visit us at https://launchpad.net/madgraph5
486 C
487 C MadGraph StandAlone Version
488 C
489 C Returns amplitude squared summed/avg over colors
490 C and helicities
491 C for the point in phase space P(0:3,NEXTERNAL)
492 C
493 C Process: u s~ > c d~ h WEIGHTED=6
494 C * Decay: h > ta+ ta- WEIGHTED=2
495 C
496  IMPLICIT NONE
497 C
498 C CONSTANTS
499 C
500  INTEGER NEXTERNAL
501  parameter(nexternal=6)
502  INTEGER NCOMB
503  parameter( ncomb=64)
504 C
505 C ARGUMENTS
506 C
507  real*8 p(0:3,nexternal),ans
508  INTEGER H1,H2
509 C
510 C LOCAL VARIABLES
511 C
512  INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
513  real*8 t
514  real*8 matrix_usx_cdx_h
515  INTEGER IHEL,IDEN, I
516  INTEGER JC(NEXTERNAL)
517  LOGICAL GOODHEL(NCOMB)
518  DATA ntry/0/
519  DATA goodhel/ncomb*.false./
520  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
521  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
522  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
523  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
524  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
525  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
526  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
527  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
528  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
529  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
530  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
531  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
532  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
533  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
534  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
535  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
536  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
537  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
538  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
539  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
540  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
541  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
542  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
543  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
544  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
545  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
546  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
547  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
548  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
549  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
550  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
551  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
552  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
553  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
554  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
555  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
556  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
557  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
558  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
559  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
560  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
561  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
562  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
563  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
564  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
565  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
566  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
567  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
568  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
569  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
570  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
571  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
572  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
573  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
574  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
575  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
576  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
577  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
578  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
579  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
580  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
581  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
582  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
583  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
584  DATA iden/36/
585 C ----------
586 C BEGIN CODE
587 C ----------
588  DO ihel=1,nexternal
589  jc(ihel) = +1
590  ENDDO
591  ans = 0d0
592  DO ihel=1,ncomb
593  t=matrix_usx_cdx_h(p ,h1,h2,nhel(1,ihel),jc(1))
594  ans=ans+t
595  ENDDO
596  ans=ans/dble(iden)
597  END
598 
599 
600  REAL*8 FUNCTION matrix_usx_cdx_h(P,H1,H2,NHEL,IC)
601 C
602 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
603 C By the MadGraph Development Team
604 C Please visit us at https://launchpad.net/madgraph5
605 C
606 C Returns amplitude squared summed/avg over colors
607 C for the point with external lines W(0:6,NEXTERNAL)
608 C
609 C Process: u s~ > c d~ h WEIGHTED=6
610 C * Decay: h > ta+ ta- WEIGHTED=2
611 C
612  IMPLICIT NONE
613 C
614 C CONSTANTS
615 C
616  INTEGER NGRAPHS
617  parameter(ngraphs=1)
618  INTEGER NEXTERNAL
619  parameter(nexternal=6)
620  INTEGER NWAVEFUNCS, NCOLOR
621  parameter(nwavefuncs=7, ncolor=1)
622  real*8 zero
623  parameter(zero=0d0)
624  COMPLEX*16 IMAG1
625  parameter(imag1=(0d0,1d0))
626 C
627 C ARGUMENTS
628 C
629  real*8 p(0:3,nexternal)
630  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
631 C
632 C LOCAL VARIABLES
633 C
634  INTEGER I,J
635  COMPLEX*16 ZTEMP
636  real*8 denom(ncolor), cf(ncolor,ncolor)
637  COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
638  COMPLEX*16 W(18,NWAVEFUNCS)
639  COMPLEX*16 DUM0,DUM1
640  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
641 C
642 C GLOBAL VARIABLES
643 C
644  include 'coupl.inc'
645 C
646 C COLOR DATA
647 C
648  DATA denom(1)/1/
649  DATA (cf(i, 1),i= 1, 1) / 9/
650 C 1 T(2,1) T(3,4)
651 
652  INTEGER H1,H2
653  real*8 matrix
654  matrix_usx_cdx_h=0.d0
655  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
656  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
657 
658 C ----------
659 C BEGIN CODE
660 C ----------
661  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
662  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
663  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
664  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
665  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
666  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
667  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
668  CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
669  CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,2))
670 C Amplitude(s) for diagram number 1
671  CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
672  jamp(1)=-amp(1)
673 
674  matrix = 0.d0
675  DO i = 1, ncolor
676  ztemp = (0.d0,0.d0)
677  DO j = 1, ncolor
678  ztemp = ztemp + cf(j,i)*jamp(j)
679  ENDDO
680  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
681  ENDDO
682  matrix_usx_cdx_h=matrix
683  ENDIF
684  ENDIF
685  END
686 
687 C ----- begin subprocesses USX->CSX with Higgs->tautau
688 C ----- THE SAME AS USX->UDX WITH Higgs->tautau
689 
690 C -------------- NOW PROCESSIS WITHOUT HIGGS
691 
692 C ----- begin subprocesses USX->UDX NO Higgs
693 
694  SUBROUTINE usx_udx_noh(P,H1,H2,ANS)
695 C
696 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
697 C By the MadGraph Development Team
698 C Please visit us at https://launchpad.net/madgraph5
699 C
700 C MadGraph StandAlone Version
701 C
702 C Returns amplitude squared summed/avg over colors
703 C and helicities
704 C for the point in phase space P(0:3,NEXTERNAL)
705 C
706 C Process: u s~ > u d~ ta+ ta- / h QED=4
707 C Process: u s~ > c s~ ta+ ta- / h QED=4
708 C
709  IMPLICIT NONE
710 C
711 C CONSTANTS
712 C
713  INTEGER NEXTERNAL
714  parameter(nexternal=6)
715  INTEGER NCOMB
716  parameter( ncomb=64)
717 C
718 C ARGUMENTS
719 C
720  real*8 p(0:3,nexternal),ans
721  INTEGER H1,H2
722 C
723 C LOCAL VARIABLES
724 C
725  INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
726  real*8 t
727  real*8 matrix_usx_udx_noh
728  INTEGER IHEL,IDEN, I
729  INTEGER JC(NEXTERNAL)
730  LOGICAL GOODHEL(NCOMB)
731  DATA ntry/0/
732  DATA goodhel/ncomb*.false./
733  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
734  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
735  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
736  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
737  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
738  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
739  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
740  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
741  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
742  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
743  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
744  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
745  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
746  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
747  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
748  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
749  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
750  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
751  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
752  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
753  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
754  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
755  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
756  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
757  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
758  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
759  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
760  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
761  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
762  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
763  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
764  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
765  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
766  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
767  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
768  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
769  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
770  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
771  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
772  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
773  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
774  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
775  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
776  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
777  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
778  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
779  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
780  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
781  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
782  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
783  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
784  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
785  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
786  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
787  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
788  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
789  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
790  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
791  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
792  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
793  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
794  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
795  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
796  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
797  DATA iden/36/
798 C ----------
799 C BEGIN CODE
800 C ----------
801  DO ihel=1,nexternal
802  jc(ihel) = +1
803  ENDDO
804  ans = 0d0
805  DO ihel=1,ncomb
806  t=matrix_usx_udx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
807  ans=ans+t
808  ENDDO
809  ans=ans/dble(iden)
810  END
811 
812 
813  REAL*8 FUNCTION matrix_usx_udx_noh(P,H1,H2,NHEL,IC)
814 C
815 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
816 C By the MadGraph Development Team
817 C Please visit us at https://launchpad.net/madgraph5
818 C
819 C Returns amplitude squared summed/avg over colors
820 C for the point with external lines W(0:6,NEXTERNAL)
821 C
822 C Process: u s~ > u d~ ta+ ta- / h QED=4
823 C Process: u s~ > c s~ ta+ ta- / h QED=4
824 C
825  IMPLICIT NONE
826 C
827 C CONSTANTS
828 C
829  INTEGER NGRAPHS
830  parameter(ngraphs=11)
831  INTEGER NEXTERNAL
832  parameter(nexternal=6)
833  INTEGER NWAVEFUNCS, NCOLOR
834  parameter(nwavefuncs=10, ncolor=1)
835  real*8 zero
836  parameter(zero=0d0)
837  COMPLEX*16 IMAG1
838  parameter(imag1=(0d0,1d0))
839 C
840 C ARGUMENTS
841 C
842  real*8 p(0:3,nexternal)
843  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
844 C
845 C LOCAL VARIABLES
846 C
847  INTEGER I,J
848  COMPLEX*16 ZTEMP
849  real*8 denom(ncolor), cf(ncolor,ncolor)
850  COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
851  COMPLEX*16 W(18,NWAVEFUNCS)
852  COMPLEX*16 DUM0,DUM1
853  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
854 C
855 C GLOBAL VARIABLES
856 C
857  include 'coupl.inc'
858 C
859 C COLOR DATA
860 C
861  DATA denom(1)/1/
862  DATA (cf(i, 1),i= 1, 1) / 9/
863 C 1 T(2,1) T(3,4)
864 
865  INTEGER H1,H2
866  real*8 matrix
867  matrix_usx_udx_noh=0.d0
868  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
869  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
870 
871 C ----------
872 C BEGIN CODE
873 C ----------
874  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
875  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
876  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
877  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
878  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
879  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
880  CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
881  CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
882  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
883 C Amplitude(s) for diagram number 1
884  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
885  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
886 C Amplitude(s) for diagram number 2
887  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
888  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
889 C Amplitude(s) for diagram number 3
890  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
891  CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
892 C Amplitude(s) for diagram number 4
893  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
894  CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
895 C Amplitude(s) for diagram number 5
896  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
897 C Amplitude(s) for diagram number 6
898  CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
899 C Amplitude(s) for diagram number 7
900  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
901  CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
902 C Amplitude(s) for diagram number 8
903  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
904  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
905 C Amplitude(s) for diagram number 9
906  CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
907 C Amplitude(s) for diagram number 10
908  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
909  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
910 C Amplitude(s) for diagram number 11
911  CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
912  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
913  $ -amp(9)-amp(10)-amp(11)
914 
915  matrix = 0.d0
916  DO i = 1, ncolor
917  ztemp = (0.d0,0.d0)
918  DO j = 1, ncolor
919  ztemp = ztemp + cf(j,i)*jamp(j)
920  ENDDO
921  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
922  ENDDO
923  matrix_usx_udx_noh=matrix
924  ENDIF
925  ENDIF
926  END
927 
928 C ----- begin subprocesses USX->USX NO Higgs
929 
930  SUBROUTINE usx_usx_noh(P,H1,H2,ANS)
931 C
932 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
933 C By the MadGraph Development Team
934 C Please visit us at https://launchpad.net/madgraph5
935 C
936 C MadGraph StandAlone Version
937 C
938 C Returns amplitude squared summed/avg over colors
939 C and helicities
940 C for the point in phase space P(0:3,NEXTERNAL)
941 C
942 C Process: u s~ > u s~ ta+ ta- / h QED=4
943 C
944  IMPLICIT NONE
945 C
946 C CONSTANTS
947 C
948  INTEGER NEXTERNAL
949  parameter(nexternal=6)
950  INTEGER NCOMB
951  parameter( ncomb=64)
952 C
953 C ARGUMENTS
954 C
955  real*8 p(0:3,nexternal),ans
956  INTEGER H1,H2
957 C
958 C LOCAL VARIABLES
959 C
960  INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
961  real*8 t
962  real*8 matrix_usx_usx_noh
963  INTEGER IHEL,IDEN, I
964  INTEGER JC(NEXTERNAL)
965  LOGICAL GOODHEL(NCOMB)
966  DATA ntry/0/
967  DATA goodhel/ncomb*.false./
968  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
969  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
970  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
971  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
972  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
973  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
974  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
975  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
976  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
977  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
978  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
979  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
980  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
981  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
982  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
983  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
984  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
985  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
986  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
987  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
988  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
989  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
990  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
991  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
992  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
993  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
994  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
995  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
996  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
997  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
998  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
999  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1000  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1001  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1002  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1003  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1004  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1005  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1006  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1007  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1008  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1009  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1010  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1011  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1012  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1013  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1014  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1015  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1016  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1017  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1018  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1019  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1020  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1021  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1022  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1023  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1024  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1025  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1026  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1027  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1028  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1029  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1030  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1031  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1032  DATA iden/36/
1033 C ----------
1034 C BEGIN CODE
1035 C ----------
1036  DO ihel=1,nexternal
1037  jc(ihel) = +1
1038  ENDDO
1039  ans = 0d0
1040  DO ihel=1,ncomb
1041  t=matrix_usx_usx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1042  ans=ans+t
1043  ENDDO
1044  ans=ans/dble(iden)
1045  END
1046 
1047 
1048  REAL*8 FUNCTION matrix_usx_usx_noh(P,H1,H2,NHEL,IC)
1049 C
1050 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1051 C By the MadGraph Development Team
1052 C Please visit us at https://launchpad.net/madgraph5
1053 C
1054 C Returns amplitude squared summed/avg over colors
1055 C for the point with external lines W(0:6,NEXTERNAL)
1056 C
1057 C Process: u s~ > u s~ ta+ ta- / h QED=4
1058 C
1059  IMPLICIT NONE
1060 C
1061 C CONSTANTS
1062 C
1063  INTEGER NGRAPHS
1064  parameter(ngraphs=43)
1065  INTEGER NEXTERNAL
1066  parameter(nexternal=6)
1067  INTEGER NWAVEFUNCS, NCOLOR
1068  parameter(nwavefuncs=15, ncolor=2)
1069  real*8 zero
1070  parameter(zero=0d0)
1071  COMPLEX*16 IMAG1
1072  parameter(imag1=(0d0,1d0))
1073 C
1074 C ARGUMENTS
1075 C
1076  real*8 p(0:3,nexternal)
1077  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1078 C
1079 C LOCAL VARIABLES
1080 C
1081  INTEGER I,J
1082  COMPLEX*16 ZTEMP
1083  real*8 denom(ncolor), cf(ncolor,ncolor)
1084  COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1085  COMPLEX*16 W(18,NWAVEFUNCS)
1086  COMPLEX*16 DUM0,DUM1
1087  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1088 C
1089 C GLOBAL VARIABLES
1090 C
1091  include 'coupl.inc'
1092 C
1093 C COLOR DATA
1094 C
1095  DATA denom(1)/1/
1096  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1097 C 1 T(2,1) T(3,4)
1098  DATA denom(2)/1/
1099  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1100 C 1 T(2,4) T(3,1)
1101 
1102  INTEGER H1,H2
1103  real*8 matrix
1104  matrix_usx_usx_noh=0.d0
1105  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1106  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1107 
1108 C ----------
1109 C BEGIN CODE
1110 C ----------
1111  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1112  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1113  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1114  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1115  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1116  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1117  CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1118  CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1119  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1120 C Amplitude(s) for diagram number 1
1121  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1122  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1123 C Amplitude(s) for diagram number 2
1124  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1125  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1126 C Amplitude(s) for diagram number 3
1127  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1128  CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,11))
1129 C Amplitude(s) for diagram number 4
1130  CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_1,amp(4))
1131  CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,12))
1132 C Amplitude(s) for diagram number 5
1133  CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(5))
1134 C Amplitude(s) for diagram number 6
1135  CALL ffv2_3_0(w(1,4),w(1,11),w(1,10),gc_50,gc_58,amp(6))
1136 C Amplitude(s) for diagram number 7
1137  CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1138  CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,12))
1139  CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,11))
1140  CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1141 C Amplitude(s) for diagram number 8
1142  CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1143  CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1144 C Amplitude(s) for diagram number 9
1145  CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1146  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,14))
1147 C Amplitude(s) for diagram number 10
1148  CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1149 C Amplitude(s) for diagram number 11
1150  CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1151  CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,13))
1152  CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1153 C Amplitude(s) for diagram number 12
1154  CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1155  CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1156 C Amplitude(s) for diagram number 13
1157  CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1158 C Amplitude(s) for diagram number 14
1159  CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1160 C Amplitude(s) for diagram number 15
1161  CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1162  CALL ffv1_1(w(1,2),w(1,12),gc_1,zero,zero,w(1,15))
1163 C Amplitude(s) for diagram number 16
1164  CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(16))
1165  CALL ffv1_2(w(1,4),w(1,12),gc_1,zero,zero,w(1,5))
1166 C Amplitude(s) for diagram number 17
1167  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(17))
1168 C Amplitude(s) for diagram number 18
1169  CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(18))
1170 C Amplitude(s) for diagram number 19
1171  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(19))
1172  CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1173  CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1174 C Amplitude(s) for diagram number 20
1175  CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(20))
1176  CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1177 C Amplitude(s) for diagram number 21
1178  CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_1,amp(21))
1179 C Amplitude(s) for diagram number 22
1180  CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(22))
1181 C Amplitude(s) for diagram number 23
1182  CALL ffv2_3_0(w(1,12),w(1,2),w(1,10),gc_50,gc_58,amp(23))
1183  CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,12))
1184 C Amplitude(s) for diagram number 24
1185  CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_1,amp(24))
1186  CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,15))
1187 C Amplitude(s) for diagram number 25
1188  CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_1,amp(25))
1189 C Amplitude(s) for diagram number 26
1190  CALL ffv2_3_0(w(1,4),w(1,12),w(1,10),gc_50,gc_58,amp(26))
1191 C Amplitude(s) for diagram number 27
1192  CALL ffv2_3_0(w(1,15),w(1,2),w(1,10),gc_50,gc_58,amp(27))
1193  CALL ffv1_2(w(1,1),w(1,11),gc_2,zero,zero,w(1,15))
1194 C Amplitude(s) for diagram number 28
1195  CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(28))
1196  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,12))
1197 C Amplitude(s) for diagram number 29
1198  CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_2,amp(29))
1199 C Amplitude(s) for diagram number 30
1200  CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(30))
1201  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,15))
1202 C Amplitude(s) for diagram number 31
1203  CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_2,amp(31))
1204  CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1205  CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1206 C Amplitude(s) for diagram number 32
1207  CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_2,amp(32))
1208 C Amplitude(s) for diagram number 33
1209  CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1210 C Amplitude(s) for diagram number 34
1211  CALL ffv2_5_0(w(1,4),w(1,3),w(1,10),gc_51,gc_58,amp(34))
1212 C Amplitude(s) for diagram number 35
1213  CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1214  CALL ffv2_5_2(w(1,1),w(1,14),gc_51,gc_58,zero,zero,w(1,11))
1215 C Amplitude(s) for diagram number 36
1216  CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_2,amp(36))
1217 C Amplitude(s) for diagram number 37
1218  CALL ffv2_5_0(w(1,12),w(1,3),w(1,14),gc_51,gc_58,amp(37))
1219 C Amplitude(s) for diagram number 38
1220  CALL ffv2_5_0(w(1,11),w(1,3),w(1,10),gc_51,gc_58,amp(38))
1221 C Amplitude(s) for diagram number 39
1222  CALL ffv2_5_0(w(1,15),w(1,3),w(1,14),gc_51,gc_58,amp(39))
1223  CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,14))
1224 C Amplitude(s) for diagram number 40
1225  CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_1,amp(40))
1226 C Amplitude(s) for diagram number 41
1227  CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_101,amp(41))
1228 C Amplitude(s) for diagram number 42
1229  CALL ffv2_3_0(w(1,14),w(1,2),w(1,10),gc_50,gc_58,amp(42))
1230 C Amplitude(s) for diagram number 43
1231  CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_101,amp(43))
1232  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1233  $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1234  $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1235  $ -amp(40)-amp(41)-amp(42)-amp(43)
1236  jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1237  $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1238  $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1239  $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1240  $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1241  $ +amp(38)+amp(39)
1242 
1243  matrix = 0.d0
1244  DO i = 1, ncolor
1245  ztemp = (0.d0,0.d0)
1246  DO j = 1, ncolor
1247  ztemp = ztemp + cf(j,i)*jamp(j)
1248  ENDDO
1249  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1250  ENDDO
1251  matrix_usx_usx_noh=matrix
1252  ENDIF
1253  ENDIF
1254  END
1255 
1256 C ----- begin subprocesses USX->CDX NO Higgs
1257 
1258  SUBROUTINE usx_cdx_noh(P,H1,H2,ANS)
1259 C
1260 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1261 C By the MadGraph Development Team
1262 C Please visit us at https://launchpad.net/madgraph5
1263 C
1264 C MadGraph StandAlone Version
1265 C
1266 C Returns amplitude squared summed/avg over colors
1267 C and helicities
1268 C for the point in phase space P(0:3,NEXTERNAL)
1269 C
1270 C Process: u s~ > c d~ ta+ ta- / h QED=4
1271 C
1272  IMPLICIT NONE
1273 C
1274 C CONSTANTS
1275 C
1276  INTEGER NEXTERNAL
1277  parameter(nexternal=6)
1278  INTEGER NCOMB
1279  parameter( ncomb=64)
1280 C
1281 C ARGUMENTS
1282 C
1283  real*8 p(0:3,nexternal),ans
1284  INTEGER H1,H2
1285 C
1286 C LOCAL VARIABLES
1287 C
1288  INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1289  real*8 t
1290  real*8 matrix_usx_cdx_noh
1291  INTEGER IHEL,IDEN, I
1292  INTEGER JC(NEXTERNAL)
1293  LOGICAL GOODHEL(NCOMB)
1294  DATA ntry/0/
1295  DATA goodhel/ncomb*.false./
1296  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1297  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1298  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1299  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1300  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1301  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1302  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1303  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1304  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1305  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1306  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1307  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1308  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1309  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1310  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1311  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1312  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1313  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1314  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1315  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1316  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1317  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1318  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1319  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1320  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1321  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1322  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1323  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1324  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1325  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1326  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1327  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1328  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1329  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1330  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1331  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1332  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1333  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1334  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1335  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1336  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1337  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1338  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1339  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1340  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1341  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1342  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1343  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1344  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1345  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1346  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1347  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1348  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1349  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1350  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1351  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1352  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1353  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1354  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1355  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1356  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1357  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1358  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1359  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1360  DATA iden/36/
1361 C ----------
1362 C BEGIN CODE
1363 C ----------
1364  DO ihel=1,nexternal
1365  jc(ihel) = +1
1366  ENDDO
1367  ans = 0d0
1368  DO ihel=1,ncomb
1369  t=matrix_usx_cdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1370  ans=ans+t
1371  ENDDO
1372 
1373  ans=ans/dble(iden)
1374  END
1375 
1376 
1377  REAL*8 FUNCTION matrix_usx_cdx_noh(P,H1,H2,NHEL,IC)
1378 C
1379 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1380 C By the MadGraph Development Team
1381 C Please visit us at https://launchpad.net/madgraph5
1382 C
1383 C Returns amplitude squared summed/avg over colors
1384 C for the point with external lines W(0:6,NEXTERNAL)
1385 C
1386 C Process: u s~ > c d~ ta+ ta- / h QED=4
1387 C
1388  IMPLICIT NONE
1389 C
1390 C CONSTANTS
1391 C
1392  INTEGER NGRAPHS
1393  parameter(ngraphs=11)
1394  INTEGER NEXTERNAL
1395  parameter(nexternal=6)
1396  INTEGER NWAVEFUNCS, NCOLOR
1397  parameter(nwavefuncs=10, ncolor=1)
1398  real*8 zero
1399  parameter(zero=0d0)
1400  COMPLEX*16 IMAG1
1401  parameter(imag1=(0d0,1d0))
1402 C
1403 C ARGUMENTS
1404 C
1405  real*8 p(0:3,nexternal)
1406  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1407 C
1408 C LOCAL VARIABLES
1409 C
1410  INTEGER I,J
1411  COMPLEX*16 ZTEMP
1412  real*8 denom(ncolor), cf(ncolor,ncolor)
1413  COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1414  COMPLEX*16 W(18,NWAVEFUNCS)
1415  COMPLEX*16 DUM0,DUM1
1416  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1417 C
1418 C GLOBAL VARIABLES
1419 C
1420  include 'coupl.inc'
1421 C
1422 C COLOR DATA
1423 C
1424  DATA denom(1)/1/
1425  DATA (cf(i, 1),i= 1, 1) / 9/
1426 C 1 T(2,1) T(3,4)
1427 
1428  INTEGER H1,H2
1429  real*8 matrix
1430  matrix_usx_cdx_noh=0.d0
1431  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1432  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1433 
1434 C ----------
1435 C BEGIN CODE
1436 C ----------
1437  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1438  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1439  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1440  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1441  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1442  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1443  CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1444  CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1445  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1446 C Amplitude(s) for diagram number 1
1447  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1448  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1449 C Amplitude(s) for diagram number 2
1450  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1451  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1452 C Amplitude(s) for diagram number 3
1453  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1454  CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,6))
1455 C Amplitude(s) for diagram number 4
1456  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1457  CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
1458 C Amplitude(s) for diagram number 5
1459  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1460 C Amplitude(s) for diagram number 6
1461  CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1462 C Amplitude(s) for diagram number 7
1463  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1464  CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1465 C Amplitude(s) for diagram number 8
1466  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1467  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1468 C Amplitude(s) for diagram number 9
1469  CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
1470 C Amplitude(s) for diagram number 10
1471  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1472  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1473 C Amplitude(s) for diagram number 11
1474  CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
1475  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1476  $ -amp(9)-amp(10)-amp(11)
1477 
1478  matrix = 0.d0
1479  DO i = 1, ncolor
1480  ztemp = (0.d0,0.d0)
1481  DO j = 1, ncolor
1482  ztemp = ztemp + cf(j,i)*jamp(j)
1483  ENDDO
1484  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1485  ENDDO
1486  matrix_usx_cdx_noh=matrix
1487  ENDIF
1488  ENDIF
1489  END
1490 
1491 C ----- begin subprocesses USX->CSX NO Higgs
1492 C IS THE SAME AS USX->CSX NO Higgs
1493 
1494 
1495 
1496 C DONE 25.10.2015 STILL S=2 MISSING
Generated by doxygen 1.9.1