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