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