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