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