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