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