C++InterfacetoTauola
CD.f
1 
2 c --------------------- subprocesses initiated by CD: 13.08.2015
3 
4  SUBROUTINE cd(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 C
12  ans=0.0d0
13 
14  IF(key.EQ.1) THEN
15  IF(i3.EQ.2 .AND. i4.EQ.1) CALL cd_ud_h(p,h1,h2,ans)
16  IF(i3.EQ.2 .AND. i4.EQ.3) CALL cd_us_h(p,h1,h2,ans)
17  IF(i3.EQ.4 .AND. i4.EQ.1) CAll cd_cd_h(p,h1,h2,ans)
18  IF(i3.EQ.4 .AND. i4.EQ.3) CALL cd_cs_h(p,h1,h2,ans)
19  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
20  CALL cd_ud_h(p,h1,h2,ans1)
21  CALL cd_us_h(p,h1,h2,ans2)
22  CALL cd_cd_h(p,h1,h2,ans3)
23  CALL cd_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 cd_ud_noh(p,h1,h2,ans)
28  IF(i3.EQ.2 .AND. i4.EQ.3) CALL cd_us_noh(p,h1,h2,ans)
29  IF(i3.EQ.4 .AND. i4.EQ.1) CAll cd_cd_noh(p,h1,h2,ans)
30  IF(i3.EQ.4 .AND. i4.EQ.3) CALL cd_cs_noh(p,h1,h2,ans)
31  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
32  CALL cd_ud_noh(p,h1,h2,ans1)
33  CALL cd_us_noh(p,h1,h2,ans2)
34  CALL cd_cd_noh(p,h1,h2,ans3)
35  CALL cd_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 CD
43 
44 
45 
46 C ----- begin subprocesses CD->jjH with Higgs->tautau
47 
48  SUBROUTINE cd_cd_h(P,H1,H2,ANS)
49 C
50 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
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: c d > 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_cd_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_cd_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_cd_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: c d > 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_cd_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 
230 C ----------
231 C BEGIN CODE
232 C ----------
233  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
234  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
235  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
236  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
237  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
238  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
239  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
240  CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,6))
241  CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
242 C Amplitude(s) for diagram number 1
243  CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
244  CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,5))
245  CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,4))
246 C Amplitude(s) for diagram number 2
247  CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_72,amp(2))
248  jamp(1)=-amp(1)
249  jamp(2)=+amp(2)
250 
251  matrix = 0.d0
252  DO i = 1, ncolor
253  ztemp = (0.d0,0.d0)
254  DO j = 1, ncolor
255  ztemp = ztemp + cf(j,i)*jamp(j)
256  ENDDO
257  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
258  ENDDO
259  matrix_cd_cd_h=matrix
260  ENDIF ! CLOSES H1 IF
261  ENDIF ! CLOSES H2 IF
262  END
263 
264 c==========================================
265 
266  SUBROUTINE cd_cs_h(P,H1,H2,ANS)
267 C
268 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
269 C By the MadGraph Development Team
270 C Please visit us at https://launchpad.net/madgraph5
271 C
272 C MadGraph StandAlone Version
273 C
274 C Returns amplitude squared summed/avg over colors
275 C and helicities
276 C for the point in phase space P(0:3,NEXTERNAL)
277 C
278 C Process: c d > c s h WEIGHTED=6
279 C * Decay: h > ta+ ta- WEIGHTED=2
280 C
281  IMPLICIT NONE
282 C
283 C CONSTANTS
284 C
285  INTEGER nexternal
286  parameter(nexternal=6)
287  INTEGER ncomb
288  parameter( ncomb=64)
289 C
290 C ARGUMENTS
291 C
292  REAL*8 p(0:3,nexternal),ans
293  INTEGER h1,h2
294 C
295 C LOCAL VARIABLES
296 C
297  INTEGER nhel(nexternal,ncomb),ntry
298  REAL*8 t
299  REAL*8 matrix_cd_cs_h,matrix
300  INTEGER ihel,iden, i
301  INTEGER jc(nexternal)
302  LOGICAL goodhel(ncomb)
303  DATA ntry/0/
304  DATA goodhel/ncomb*.false./
305  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
306  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
307  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
308  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
309  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
310  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
311  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
312  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
313  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
314  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
315  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
316  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
317  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
318  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
319  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
320  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
321  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
322  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
323  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
324  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
325  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
326  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
327  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
328  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
329  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
330  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
331  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
332  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
333  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
334  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
335  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
336  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
337  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
338  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
339  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
340  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
341  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
342  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
343  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
344  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
345  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
346  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
347  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
348  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
349  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
350  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
351  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
352  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
353  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
354  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
355  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
356  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
357  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
358  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
359  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
360  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
361  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
362  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
363  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
364  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
365  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
366  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
367  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
368  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
369  DATA iden/36/
370 C ----------
371 C BEGIN CODE
372 C ----------
373  DO ihel=1,nexternal
374  jc(ihel) = +1
375  ENDDO
376  ans = 0d0
377  DO ihel=1,ncomb
378  t=matrix_cd_cs_h(p,h1,h2 ,nhel(1,ihel),jc(1))
379  ans=ans+t
380  ENDDO
381  ans=ans/dble(iden)
382  END
383 
384 
385  REAL*8 FUNCTION matrix_cd_cs_h(P,H1,H2,NHEL,IC)
386 C
387 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
388 C By the MadGraph Development Team
389 C Please visit us at https://launchpad.net/madgraph5
390 C
391 C Returns amplitude squared summed/avg over colors
392 C for the point with external lines W(0:6,NEXTERNAL)
393 C
394 C Process: c d > c s h WEIGHTED=6
395 C * Decay: h > ta+ ta- WEIGHTED=2
396 C
397  IMPLICIT NONE
398 C
399 C CONSTANTS
400 C
401  INTEGER ngraphs
402  parameter(ngraphs=1)
403  INTEGER nexternal
404  parameter(nexternal=6)
405  INTEGER nwavefuncs, ncolor
406  parameter(nwavefuncs=7, ncolor=1)
407  REAL*8 zero
408  parameter(zero=0d0)
409  COMPLEX*16 imag1
410  parameter(imag1=(0d0,1d0))
411 C
412 C ARGUMENTS
413 C
414  REAL*8 p(0:3,nexternal)
415  INTEGER nhel(nexternal), ic(nexternal)
416 C
417 C LOCAL VARIABLES
418 C
419  INTEGER i,j
420  COMPLEX*16 ztemp
421  REAL*8 denom(ncolor), cf(ncolor,ncolor)
422  COMPLEX*16 amp(ngraphs), jamp(ncolor)
423  COMPLEX*16 w(18,nwavefuncs)
424  COMPLEX*16 dum0,dum1
425  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
426 C
427 C GLOBAL VARIABLES
428 C
429  include 'coupl.inc'
430 C
431 C COLOR DATA
432 C
433  DATA denom(1)/1/
434  DATA (cf(i, 1),i= 1, 1) / 9/
435 C 1 T(3,2) T(4,1)
436 
437  INTEGER h1,h2
438  REAL*8 matrix
439  matrix_cd_cs_h=0.d0
440  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
441  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
442 
443 C ----------
444 C BEGIN CODE
445 C ----------
446  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
447  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
448  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
449  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
450  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
451  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
452  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
453  CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,6))
454  CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,4))
455 C Amplitude(s) for diagram number 1
456  CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
457  jamp(1)=+amp(1)
458 
459  matrix = 0.d0
460  DO i = 1, ncolor
461  ztemp = (0.d0,0.d0)
462  DO j = 1, ncolor
463  ztemp = ztemp + cf(j,i)*jamp(j)
464  ENDDO
465  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
466  ENDDO
467  matrix_cd_cs_h=matrix
468  ENDIF
469  ENDIF
470  END
471 
472 
473 
474 C===============================================
475 
476  SUBROUTINE cd_ud_h(P,H1,H2,ANS)
477 C
478 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
479 C By the MadGraph Development Team
480 C Please visit us at https://launchpad.net/madgraph5
481 C
482 C MadGraph StandAlone Version
483 C
484 C Returns amplitude squared summed/avg over colors
485 C and helicities
486 C for the point in phase space P(0:3,NEXTERNAL)
487 C
488 C Process: c d > u d h WEIGHTED=6
489 C * Decay: h > ta+ ta- WEIGHTED=2
490 C
491  IMPLICIT NONE
492 C
493 C CONSTANTS
494 C
495  INTEGER nexternal
496  parameter(nexternal=6)
497  INTEGER ncomb
498  parameter( ncomb=64)
499 C
500 C ARGUMENTS
501 C
502  REAL*8 p(0:3,nexternal),ans
503  INTEGER h1,h2
504 C
505 C LOCAL VARIABLES
506 C
507  INTEGER nhel(nexternal,ncomb),ntry
508  REAL*8 t
509  REAL*8 matrix_cd_ud_h
510  INTEGER ihel,iden, i
511  INTEGER jc(nexternal)
512  LOGICAL goodhel(ncomb)
513  DATA ntry/0/
514  DATA goodhel/ncomb*.false./
515  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
516  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
517  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
518  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
519  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
520  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
521  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
522  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
523  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
524  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
525  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
526  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
527  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
528  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
529  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
530  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
531  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
532  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
533  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
534  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
535  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
536  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
537  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
538  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
539  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
540  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
541  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
542  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
543  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
544  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
545  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
546  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
547  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
548  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
549  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
550  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
551  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
552  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
553  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
554  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
555  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
556  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
557  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
558  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
559  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
560  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
561  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
562  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
563  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
564  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
565  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
566  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
567  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
568  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
569  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
570  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
571  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
572  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
573  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
574  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
575  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
576  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
577  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
578  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
579  DATA iden/36/
580 C ----------
581 C BEGIN CODE
582 C ----------
583  DO ihel=1,nexternal
584  jc(ihel) = +1
585  ENDDO
586  ans = 0d0
587  DO ihel=1,ncomb
588  t=matrix_cd_ud_h(p,h1,h2 ,nhel(1,ihel),jc(1))
589  ans=ans+t
590  ENDDO
591  ans=ans/dble(iden)
592  END
593 
594 
595  REAL*8 FUNCTION matrix_cd_ud_h(P,H1,H2,NHEL,IC)
596 C
597 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
598 C By the MadGraph Development Team
599 C Please visit us at https://launchpad.net/madgraph5
600 C
601 C Returns amplitude squared summed/avg over colors
602 C for the point with external lines W(0:6,NEXTERNAL)
603 C
604 C Process: c d > u d h WEIGHTED=6
605 C * Decay: h > ta+ ta- WEIGHTED=2
606 C
607  IMPLICIT NONE
608 C
609 C CONSTANTS
610 C
611  INTEGER ngraphs
612  parameter(ngraphs=1)
613  INTEGER nexternal
614  parameter(nexternal=6)
615  INTEGER nwavefuncs, ncolor
616  parameter(nwavefuncs=7, ncolor=1)
617  REAL*8 zero
618  parameter(zero=0d0)
619  COMPLEX*16 imag1
620  parameter(imag1=(0d0,1d0))
621 C
622 C ARGUMENTS
623 C
624  REAL*8 p(0:3,nexternal)
625  INTEGER nhel(nexternal), ic(nexternal)
626 C
627 C LOCAL VARIABLES
628 C
629  INTEGER i,j
630  COMPLEX*16 ztemp
631  REAL*8 denom(ncolor), cf(ncolor,ncolor)
632  COMPLEX*16 amp(ngraphs), jamp(ncolor)
633  COMPLEX*16 w(18,nwavefuncs)
634  COMPLEX*16 dum0,dum1
635  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
636 C
637 C GLOBAL VARIABLES
638 C
639  include 'coupl.inc'
640 C
641 C COLOR DATA
642 C
643  DATA denom(1)/1/
644  DATA (cf(i, 1),i= 1, 1) / 9/
645 C 1 T(3,2) T(4,1)
646 
647  INTEGER h1,h2
648  REAL*8 matrix
649  matrix_cd_ud_h=0d0
650  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
651  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
652 
653 C ----------
654 C BEGIN CODE
655 C ----------
656  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
657  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
658  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
659  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
660  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
661  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
662  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
663  CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,6))
664  CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,4))
665 C Amplitude(s) for diagram number 1
666  CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
667  jamp(1)=+amp(1)
668 
669  matrix = 0.d0
670  DO i = 1, ncolor
671  ztemp = (0.d0,0.d0)
672  DO j = 1, ncolor
673  ztemp = ztemp + cf(j,i)*jamp(j)
674  ENDDO
675  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
676  ENDDO
677  matrix_cd_ud_h=matrix
678  ENDIF
679  ENDIF
680  END
681 
682 
683 
684 C========================================
685 
686  SUBROUTINE cd_us_h(P,H1,H2,ANS)
687 C
688 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
689 C By the MadGraph Development Team
690 C Please visit us at https://launchpad.net/madgraph5
691 C
692 C MadGraph StandAlone Version
693 C
694 C Returns amplitude squared summed/avg over colors
695 C and helicities
696 C for the point in phase space P(0:3,NEXTERNAL)
697 C
698 C Process: c d > u s h WEIGHTED=6
699 C * Decay: h > ta+ ta- WEIGHTED=2
700 C
701  IMPLICIT NONE
702 C
703 C CONSTANTS
704 C
705  INTEGER nexternal
706  parameter(nexternal=6)
707  INTEGER ncomb
708  parameter( ncomb=64)
709 C
710 C ARGUMENTS
711 C
712  REAL*8 p(0:3,nexternal),ans
713  INTEGER h1,h2
714 C
715 C LOCAL VARIABLES
716 C
717  INTEGER nhel(nexternal,ncomb),ntry
718  REAL*8 t
719  REAL*8 matrix_cd_us_h
720  INTEGER ihel,iden, i
721  INTEGER jc(nexternal)
722  LOGICAL goodhel(ncomb)
723  DATA ntry/0/
724  DATA goodhel/ncomb*.false./
725  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
726  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
727  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
728  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
729  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
730  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
731  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
732  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
733  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
734  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
735  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
736  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
737  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
738  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
739  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
740  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
741  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
742  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
743  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
744  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
745  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
746  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
747  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
748  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
749  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
750  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
751  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
752  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
753  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
754  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
755  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
756  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
757  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
758  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
759  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
760  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
761  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
762  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
763  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
764  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
765  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
766  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
767  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
768  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
769  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
770  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
771  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
772  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
773  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
774  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
775  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
776  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
777  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
778  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
779  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
780  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
781  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
782  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
783  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
784  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
785  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
786  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
787  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
788  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
789  DATA iden/36/
790 C ----------
791 C BEGIN CODE
792 C ----------
793  DO ihel=1,nexternal
794  jc(ihel) = +1
795  ENDDO
796  ans = 0d0
797  DO ihel=1,ncomb
798  t=matrix_cd_us_h(p,h1,h2 ,nhel(1,ihel),jc(1))
799  ans=ans+t
800  ENDDO
801  ans=ans/dble(iden)
802  END
803 
804 
805  REAL*8 FUNCTION matrix_cd_us_h(P,H1,H2,NHEL,IC)
806 C
807 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
808 C By the MadGraph Development Team
809 C Please visit us at https://launchpad.net/madgraph5
810 C
811 C Returns amplitude squared summed/avg over colors
812 C for the point with external lines W(0:6,NEXTERNAL)
813 C
814 C Process: c d > u s h WEIGHTED=6
815 C * Decay: h > ta+ ta- WEIGHTED=2
816 C
817  IMPLICIT NONE
818 C
819 C CONSTANTS
820 C
821  INTEGER ngraphs
822  parameter(ngraphs=1)
823  INTEGER nexternal
824  parameter(nexternal=6)
825  INTEGER nwavefuncs, ncolor
826  parameter(nwavefuncs=7, ncolor=1)
827  REAL*8 zero
828  parameter(zero=0d0)
829  COMPLEX*16 imag1
830  parameter(imag1=(0d0,1d0))
831 C
832 C ARGUMENTS
833 C
834  REAL*8 p(0:3,nexternal)
835  INTEGER nhel(nexternal), ic(nexternal)
836 C
837 C LOCAL VARIABLES
838 C
839  INTEGER i,j
840  COMPLEX*16 ztemp
841  REAL*8 denom(ncolor), cf(ncolor,ncolor)
842  COMPLEX*16 amp(ngraphs), jamp(ncolor)
843  COMPLEX*16 w(18,nwavefuncs)
844  COMPLEX*16 dum0,dum1
845  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
846 C
847 C GLOBAL VARIABLES
848 C
849  include 'coupl.inc'
850 C
851 C COLOR DATA
852 C
853  DATA denom(1)/1/
854  DATA (cf(i, 1),i= 1, 1) / 9/
855 C 1 T(3,2) T(4,1)
856 
857  INTEGER h1,h2
858  REAL*8 matrix
859  matrix_cd_us_h=0d0
860  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
861  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
862 
863 C ----------
864 C BEGIN CODE
865 C ----------
866  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
867  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
868  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
869  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
870  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
871  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
872  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
873  CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,6))
874  CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,4))
875 C Amplitude(s) for diagram number 1
876  CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
877  jamp(1)=+amp(1)
878 
879  matrix = 0.d0
880  DO i = 1, ncolor
881  ztemp = (0.d0,0.d0)
882  DO j = 1, ncolor
883  ztemp = ztemp + cf(j,i)*jamp(j)
884  ENDDO
885  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
886  ENDDO
887  matrix_cd_us_h=matrix
888  ENDIF
889  ENDIF
890  END
891 
892 
893 
894 
895 C --- end of subprocess CD->jjH, H->tautau
896 
897 C ----- begin subprocesses CD->jjtautau, no Higgs
898 
899  SUBROUTINE cd_cd_noh(P,H1,H2,ANS)
900 C
901 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
902 C By the MadGraph Development Team
903 C Please visit us at https://launchpad.net/madgraph5
904 C
905 C MadGraph StandAlone Version
906 C
907 C Returns amplitude squared summed/avg over colors
908 C and helicities
909 C for the point in phase space P(0:3,NEXTERNAL)
910 C
911 C Process: c d > c d ta+ ta- / h QED=4
912 C
913  IMPLICIT NONE
914 C
915 C CONSTANTS
916 C
917  INTEGER nexternal
918  parameter(nexternal=6)
919  INTEGER ncomb
920  parameter( ncomb=64)
921 C
922 C ARGUMENTS
923 C
924  REAL*8 p(0:3,nexternal),ans
925  INTEGER h1,h2
926 C
927 C LOCAL VARIABLES
928 C
929  INTEGER nhel(nexternal,ncomb),ntry
930  REAL*8 t
931  REAL*8 matrix_cd_cd_noh
932  INTEGER ihel,iden, i
933  INTEGER jc(nexternal)
934  LOGICAL goodhel(ncomb)
935  DATA ntry/0/
936  DATA goodhel/ncomb*.false./
937  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
938  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
939  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
940  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
941  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
942  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
943  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
944  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
945  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
946  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
947  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
948  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
949  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
950  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
951  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
952  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
953  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
954  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
955  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
956  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
957  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
958  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
959  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
960  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
961  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
962  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
963  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
964  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
965  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
966  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
967  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
968  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
969  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
970  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
971  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
972  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
973  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
974  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
975  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
976  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
977  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
978  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
979  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
980  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
981  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
982  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
983  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
984  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
985  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
986  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
987  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
988  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
989  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
990  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
991  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
992  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
993  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
994  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
995  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
996  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
997  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
998  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
999  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1000  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1001  DATA iden/36/
1002 C ----------
1003 C BEGIN CODE
1004 C ----------
1005  DO ihel=1,nexternal
1006  jc(ihel) = +1
1007  ENDDO
1008  ans = 0d0
1009  DO ihel=1,ncomb
1010  t=matrix_cd_cd_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1011  ans=ans+t
1012  ENDDO
1013  ans=ans/dble(iden)
1014  END
1015 
1016 
1017  REAL*8 FUNCTION matrix_cd_cd_noh(P,H1,H2,NHEL,IC)
1018 C
1019 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1020 C By the MadGraph Development Team
1021 C Please visit us at https://launchpad.net/madgraph5
1022 C
1023 C Returns amplitude squared summed/avg over colors
1024 C for the point with external lines W(0:6,NEXTERNAL)
1025 C
1026 C Process: c d > c d ta+ ta- / h QED=4
1027 C
1028  IMPLICIT NONE
1029 C
1030 C CONSTANTS
1031 C
1032  INTEGER ngraphs
1033  parameter(ngraphs=43)
1034  INTEGER nexternal
1035  parameter(nexternal=6)
1036  INTEGER nwavefuncs, ncolor
1037  parameter(nwavefuncs=14, ncolor=2)
1038  REAL*8 zero
1039  parameter(zero=0d0)
1040  COMPLEX*16 imag1
1041  parameter(imag1=(0d0,1d0))
1042 C
1043 C ARGUMENTS
1044 C
1045  REAL*8 p(0:3,nexternal)
1046  INTEGER nhel(nexternal), ic(nexternal)
1047 C
1048 C LOCAL VARIABLES
1049 C
1050  INTEGER i,j
1051  COMPLEX*16 ztemp
1052  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1053  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1054  COMPLEX*16 w(18,nwavefuncs)
1055  COMPLEX*16 dum0,dum1
1056  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1057 C
1058 C GLOBAL VARIABLES
1059 C
1060  include 'coupl.inc'
1061 C
1062 C COLOR DATA
1063 C
1064  DATA denom(1)/1/
1065  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1066 C 1 T(3,1) T(4,2)
1067  DATA denom(2)/1/
1068  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1069 C 1 T(3,2) T(4,1)
1070 
1071  INTEGER h1,h2
1072  REAL*8 matrix
1073  matrix_cd_cd_noh=0d0
1074  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1075  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1076 
1077 
1078 C ----------
1079 C BEGIN CODE
1080 C ----------
1081  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1082  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1083  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1084  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1085  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1086  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1087  CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,7))
1088  CALL ffv1p0_3(w(1,2),w(1,4),gc_1,zero,zero,w(1,8))
1089  CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1090 C Amplitude(s) for diagram number 1
1091  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1092  CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1093 C Amplitude(s) for diagram number 2
1094  CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1095  CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,11))
1096 C Amplitude(s) for diagram number 3
1097  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1098 C Amplitude(s) for diagram number 4
1099  CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1100  CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,10))
1101  CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1102 C Amplitude(s) for diagram number 5
1103  CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1104  CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1105 C Amplitude(s) for diagram number 6
1106  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1107 C Amplitude(s) for diagram number 7
1108  CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1109 C Amplitude(s) for diagram number 8
1110  CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1111  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1112  CALL ffv1_2(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
1113 C Amplitude(s) for diagram number 9
1114  CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(9))
1115  CALL ffv1_1(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1116 C Amplitude(s) for diagram number 10
1117  CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_1,amp(10))
1118  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1119 C Amplitude(s) for diagram number 11
1120  CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(11))
1121 C Amplitude(s) for diagram number 12
1122  CALL ffv2_3_0(w(1,2),w(1,13),w(1,7),gc_50,gc_58,amp(12))
1123  CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
1124  CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
1125 C Amplitude(s) for diagram number 13
1126  CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(13))
1127  CALL ffv1_1(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1128 C Amplitude(s) for diagram number 14
1129  CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_1,amp(14))
1130 C Amplitude(s) for diagram number 15
1131  CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(15))
1132 C Amplitude(s) for diagram number 16
1133  CALL ffv2_3_0(w(1,2),w(1,14),w(1,7),gc_50,gc_58,amp(16))
1134  CALL ffv2_3_2(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
1135 C Amplitude(s) for diagram number 17
1136  CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_1,amp(17))
1137  CALL ffv2_3_1(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1138 C Amplitude(s) for diagram number 18
1139  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
1140 C Amplitude(s) for diagram number 19
1141  CALL ffv2_3_0(w(1,14),w(1,4),w(1,7),gc_50,gc_58,amp(19))
1142 C Amplitude(s) for diagram number 20
1143  CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(20))
1144  CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,9))
1145  CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,14))
1146  CALL ffv2_2(w(1,5),w(1,9),gc_108,zero,zero,w(1,10))
1147 C Amplitude(s) for diagram number 21
1148  CALL ffv2_0(w(1,10),w(1,6),w(1,14),gc_108,amp(21))
1149 C Amplitude(s) for diagram number 22
1150  CALL vvv1_0(w(1,12),w(1,9),w(1,14),gc_4,amp(22))
1151 C Amplitude(s) for diagram number 23
1152  CALL vvv1_0(w(1,9),w(1,14),w(1,7),gc_53,amp(23))
1153  CALL ffv2_2(w(1,2),w(1,9),gc_44,zero,zero,w(1,10))
1154 C Amplitude(s) for diagram number 24
1155  CALL ffv1_0(w(1,10),w(1,3),w(1,12),gc_2,amp(24))
1156  CALL ffv2_1(w(1,3),w(1,9),gc_44,zero,zero,w(1,6))
1157 C Amplitude(s) for diagram number 25
1158  CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(25))
1159 C Amplitude(s) for diagram number 26
1160  CALL ffv2_5_0(w(1,10),w(1,3),w(1,7),gc_51,gc_58,amp(26))
1161 C Amplitude(s) for diagram number 27
1162  CALL ffv2_3_0(w(1,2),w(1,6),w(1,7),gc_50,gc_58,amp(27))
1163  CALL ffv2_2(w(1,1),w(1,14),gc_44,zero,zero,w(1,6))
1164 C Amplitude(s) for diagram number 28
1165  CALL ffv1_0(w(1,6),w(1,4),w(1,12),gc_1,amp(28))
1166  CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,10))
1167 C Amplitude(s) for diagram number 29
1168  CALL ffv2_0(w(1,10),w(1,4),w(1,14),gc_44,amp(29))
1169 C Amplitude(s) for diagram number 30
1170  CALL ffv2_3_0(w(1,6),w(1,4),w(1,7),gc_50,gc_58,amp(30))
1171  CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,6))
1172 C Amplitude(s) for diagram number 31
1173  CALL ffv2_0(w(1,6),w(1,4),w(1,14),gc_44,amp(31))
1174  CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,14))
1175 C Amplitude(s) for diagram number 32
1176  CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(32))
1177 C Amplitude(s) for diagram number 33
1178  CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_2,amp(33))
1179 C Amplitude(s) for diagram number 34
1180  CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(34))
1181 C Amplitude(s) for diagram number 35
1182  CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_2,amp(35))
1183  CALL ffv1p0_3(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
1184  CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1185 C Amplitude(s) for diagram number 36
1186  CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(36))
1187 C Amplitude(s) for diagram number 37
1188  CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_11,amp(37))
1189 C Amplitude(s) for diagram number 38
1190  CALL ffv2_5_0(w(1,4),w(1,3),w(1,7),gc_51,gc_58,amp(38))
1191 C Amplitude(s) for diagram number 39
1192  CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(39))
1193  CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
1194 C Amplitude(s) for diagram number 40
1195  CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(40))
1196 C Amplitude(s) for diagram number 41
1197  CALL ffv2_5_0(w(1,10),w(1,3),w(1,11),gc_51,gc_58,amp(41))
1198 C Amplitude(s) for diagram number 42
1199  CALL ffv2_5_0(w(1,8),w(1,3),w(1,7),gc_51,gc_58,amp(42))
1200 C Amplitude(s) for diagram number 43
1201  CALL ffv2_5_0(w(1,6),w(1,3),w(1,11),gc_51,gc_58,amp(43))
1202  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1203  $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1204  $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1205  $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
1206  $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
1207  jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1208  $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1209  $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
1210  $ -1./2.*amp(38)-1./2.*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_cd_cd_noh=matrix
1221  ENDIF ! CLOSES H1 IF
1222  ENDIF ! CLOSES H2 IF
1223  END
1224 
1225 
1226  SUBROUTINE cd_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: c d > 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_cd_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_cd_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_cd_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: c d > 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,2) T(4,1)
1394 
1395  INTEGER h1,h2
1396  REAL*8 matrix
1397  matrix_cd_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 
1402 C ----------
1403 C BEGIN CODE
1404 C ----------
1405  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1406  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1407  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1408  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1409  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1410  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1411  CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,7))
1412  CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,8))
1413  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1414 C Amplitude(s) for diagram number 1
1415  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1416  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1417 C Amplitude(s) for diagram number 2
1418  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1419  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1420 C Amplitude(s) for diagram number 3
1421  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1422  CALL ffv2_2(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
1423 C Amplitude(s) for diagram number 4
1424  CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1425  CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,5))
1426 C Amplitude(s) for diagram number 5
1427  CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1428 C Amplitude(s) for diagram number 6
1429  CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1430 C Amplitude(s) for diagram number 7
1431  CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1432  CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1433 C Amplitude(s) for diagram number 8
1434  CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1435  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1436 C Amplitude(s) for diagram number 9
1437  CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_100,amp(9))
1438 C Amplitude(s) for diagram number 10
1439  CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1440  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1441 C Amplitude(s) for diagram number 11
1442  CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_100,amp(11))
1443  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1444  $ +amp(9)+amp(10)+amp(11)
1445 
1446  matrix = 0.d0
1447  DO i = 1, ncolor
1448  ztemp = (0.d0,0.d0)
1449  DO j = 1, ncolor
1450  ztemp = ztemp + cf(j,i)*jamp(j)
1451  ENDDO
1452  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1453  ENDDO
1454  matrix_cd_cs_noh=matrix
1455  ENDIF ! CLOSES H1 IF
1456  ENDIF ! CLOSES H2 IF
1457  END
1458 
1459  SUBROUTINE cd_ud_noh(P,H1,H2,ANS)
1460 C
1461 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1462 C By the MadGraph Development Team
1463 C Please visit us at https://launchpad.net/madgraph5
1464 C
1465 C MadGraph StandAlone Version
1466 C
1467 C Returns amplitude squared summed/avg over colors
1468 C and helicities
1469 C for the point in phase space P(0:3,NEXTERNAL)
1470 C
1471 C Process: c d > u d ta+ ta- / h QED=4
1472 C
1473  IMPLICIT NONE
1474 C
1475 C CONSTANTS
1476 C
1477  INTEGER nexternal
1478  parameter(nexternal=6)
1479  INTEGER ncomb
1480  parameter( ncomb=64)
1481 C
1482 C ARGUMENTS
1483 C
1484  REAL*8 p(0:3,nexternal),ans
1485  INTEGER h1,h2
1486 C
1487 C LOCAL VARIABLES
1488 C
1489  INTEGER nhel(nexternal,ncomb),ntry
1490  REAL*8 t
1491  REAL*8 matrix_cd_ud_noh
1492  INTEGER ihel,iden, i
1493  INTEGER jc(nexternal)
1494  LOGICAL goodhel(ncomb)
1495  DATA ntry/0/
1496  DATA goodhel/ncomb*.false./
1497  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1498  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1499  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1500  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1501  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1502  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1503  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1504  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1505  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1506  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1507  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1508  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1509  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1510  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1511  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1512  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1513  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1514  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1515  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1516  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1517  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1518  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1519  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1520  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1521  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1522  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1523  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1524  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1525  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1526  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1527  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1528  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1529  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1530  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1531  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1532  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1533  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1534  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1535  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1536  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1537  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1538  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1539  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1540  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1541  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1542  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1543  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1544  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1545  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1546  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1547  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1548  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1549  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1550  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1551  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1552  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1553  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1554  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1555  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1556  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1557  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1558  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1559  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1560  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1561  DATA iden/36/
1562 C ----------
1563 C BEGIN CODE
1564 C ----------
1565  DO ihel=1,nexternal
1566  jc(ihel) = +1
1567  ENDDO
1568  ans = 0d0
1569  DO ihel=1,ncomb
1570  t=matrix_cd_ud_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1571  ans=ans+t
1572  ENDDO
1573  ans=ans/dble(iden)
1574  END
1575 
1576 
1577  REAL*8 FUNCTION matrix_cd_ud_noh(P,H1,H2,NHEL,IC)
1578 C
1579 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1580 C By the MadGraph Development Team
1581 C Please visit us at https://launchpad.net/madgraph5
1582 C
1583 C Returns amplitude squared summed/avg over colors
1584 C for the point with external lines W(0:6,NEXTERNAL)
1585 C
1586 C Process: c d > u d ta+ ta- / h QED=4
1587 C
1588  IMPLICIT NONE
1589 C
1590 C CONSTANTS
1591 C
1592  INTEGER ngraphs
1593  parameter(ngraphs=11)
1594  INTEGER nexternal
1595  parameter(nexternal=6)
1596  INTEGER nwavefuncs, ncolor
1597  parameter(nwavefuncs=10, ncolor=1)
1598  REAL*8 zero
1599  parameter(zero=0d0)
1600  COMPLEX*16 imag1
1601  parameter(imag1=(0d0,1d0))
1602 C
1603 C ARGUMENTS
1604 C
1605  REAL*8 p(0:3,nexternal)
1606  INTEGER nhel(nexternal), ic(nexternal)
1607 C
1608 C LOCAL VARIABLES
1609 C
1610  INTEGER i,j
1611  COMPLEX*16 ztemp
1612  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1613  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1614  COMPLEX*16 w(18,nwavefuncs)
1615  COMPLEX*16 dum0,dum1
1616  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1617 C
1618 C GLOBAL VARIABLES
1619 C
1620  include 'coupl.inc'
1621 C
1622 C COLOR DATA
1623 C
1624  DATA denom(1)/1/
1625  DATA (cf(i, 1),i= 1, 1) / 9/
1626 C 1 T(3,2) T(4,1)
1627 
1628  INTEGER h1,h2
1629  REAL*8 matrix
1630  matrix_cd_ud_noh=0d0
1631  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1632  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
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,4),gc_44,mw,ww,w(1,7))
1644  CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,8))
1645  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1646 C Amplitude(s) for diagram number 1
1647  CALL ffv2_0(w(1,9),w(1,6),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,7),w(1,8),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,7),w(1,8),w(1,10),gc_53,amp(3))
1654  CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1655 C Amplitude(s) for diagram number 4
1656  CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1657  CALL ffv2_1(w(1,3),w(1,7),gc_100,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_1,amp(5))
1660 C Amplitude(s) for diagram number 6
1661  CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1662 C Amplitude(s) for diagram number 7
1663  CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1664  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1665 C Amplitude(s) for diagram number 8
1666  CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1667  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1668 C Amplitude(s) for diagram number 9
1669  CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_44,amp(9))
1670 C Amplitude(s) for diagram number 10
1671  CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1672  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1673 C Amplitude(s) for diagram number 11
1674  CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_44,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_cd_ud_noh=matrix
1687  ENDIF ! CLOSES H1 IF
1688  ENDIF ! CLOSES H2 IF
1689  END
1690 
1691 
1692  SUBROUTINE cd_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: c d > 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_cd_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_cd_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_cd_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: c d > 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,2) T(4,1)
1860 
1861  INTEGER h1,h2
1862  REAL*8 matrix
1863  matrix_cd_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 
1868 C ----------
1869 C BEGIN CODE
1870 C ----------
1871  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1872  CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1873  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1874  CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1875  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1876  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1877  CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,7))
1878  CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,8))
1879  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1880 C Amplitude(s) for diagram number 1
1881  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1882  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1883 C Amplitude(s) for diagram number 2
1884  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1885  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1886 C Amplitude(s) for diagram number 3
1887  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1888  CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1889 C Amplitude(s) for diagram number 4
1890  CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1891  CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,5))
1892 C Amplitude(s) for diagram number 5
1893  CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1894 C Amplitude(s) for diagram number 6
1895  CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1896 C Amplitude(s) for diagram number 7
1897  CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1898  CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1899 C Amplitude(s) for diagram number 8
1900  CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1901  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1902 C Amplitude(s) for diagram number 9
1903  CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_100,amp(9))
1904 C Amplitude(s) for diagram number 10
1905  CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1906  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1907 C Amplitude(s) for diagram number 11
1908  CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_100,amp(11))
1909  jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1910  $ +amp(9)+amp(10)+amp(11)
1911 
1912  matrix = 0.d0
1913  DO i = 1, ncolor
1914  ztemp = (0.d0,0.d0)
1915  DO j = 1, ncolor
1916  ztemp = ztemp + cf(j,i)*jamp(j)
1917  ENDDO
1918  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1919  ENDDO
1920  matrix_cd_us_noh=matrix
1921  ENDIF ! CLOSES H1 IF
1922  ENDIF ! CLOSES H2 IF
1923  END
1924 
1925 
1926 C ------- end subprocesses CD->jjtautau, no Higgs
1927 
1928 
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7