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