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