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