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