C++InterfacetoTauola
CDX.f
1 c --------------------- subprocesses initiated by CDbar
2 
3  SUBROUTINE cdx(P,I3,I4,H1,H2,KEY,ANS)
4  IMPLICIT NONE
5 
6  INTEGER i3,i4, h1,h2, key
7  REAL*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
8 
9 C I3=0,2 OR 4 I4=0,-1 OR -3
10 
11  ans=0.d0
12 
13 C FINAL CSX IS EQUAL TO FINAL UDX
14 
15  IF(key.EQ.1) THEN
16  IF(i3.EQ.2 .AND. i4.EQ.-1) CALL cdx_udx_h(p,h1,h2,ans)
17  IF(i3.EQ.2 .AND. i4.EQ.-3) CALL cdx_usx_h(p,h1,h2,ans)
18  IF(i3.EQ.4 .AND. i4.EQ.-1) CAll cdx_cdx_h(p,h1,h2,ans)
19  IF(i3.EQ.4 .AND. i4.EQ.-3) CALL cdx_udx_h(p,h1,h2,ans)
20  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
21  CALL cdx_udx_h(p,h1,h2,ans1)
22  CALL cdx_usx_h(p,h1,h2,ans2)
23  CALL cdx_cdx_h(p,h1,h2,ans3)
24  CALL cdx_udx_h(p,h1,h2,ans4)
25  ans=ans1+ans2+ans3+ans4
26  ENDIF
27  ELSE IF(key.EQ.0) THEN
28  IF(i3.EQ.2 .AND. i4.EQ.-1) CALL cdx_udx_noh(p,h1,h2,ans)
29  IF(i3.EQ.2 .AND. i4.EQ.-3) CALL cdx_usx_noh(p,h1,h2,ans)
30  IF(i3.EQ.4 .AND. i4.EQ.-1) CAll cdx_cdx_noh(p,h1,h2,ans)
31  IF(i3.EQ.4 .AND. i4.EQ.-3) CALL cdx_udx_noh(p,h1,h2,ans)
32  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
33  CALL cdx_udx_noh(p,h1,h2,ans1)
34  CALL cdx_usx_noh(p,h1,h2,ans2)
35  CALL cdx_cdx_noh(p,h1,h2,ans3)
36  CALL cdx_udx_noh(p,h1,h2,ans4)
37  ans=ans1+ans2+ans3+ans4
38  ENDIF
39  ELSE
40  WRITE(*,*) 'NOT FINISHED'
41  stop
42  ENDIF
43  END ! SUBROUTINE CDX
44 
45 
46 
47 C ----- begin subprocesses CDX->UDX with Higgs->tautau
48 
49 
50  SUBROUTINE cdx_udx_h(P,H1,H2,ANS)
51 C
52 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
53 C By the MadGraph Development Team
54 C Please visit us at https://launchpad.net/madgraph5
55 C
56 C MadGraph StandAlone Version
57 C
58 C Returns amplitude squared summed/avg over colors
59 C and helicities
60 C for the point in phase space P(0:3,NEXTERNAL)
61 C
62 C Process: c d~ > u d~ h WEIGHTED=6
63 C * Decay: h > ta+ ta- WEIGHTED=2
64 C Process: c d~ > c s~ h WEIGHTED=6
65 C * Decay: h > ta+ ta- WEIGHTED=2
66 C
67  IMPLICIT NONE
68 C
69 C CONSTANTS
70 C
71  INTEGER nexternal
72  parameter(nexternal=6)
73  INTEGER ncomb
74  parameter( ncomb=64)
75 C
76 C ARGUMENTS
77 C
78  REAL*8 p(0:3,nexternal),ans
79  INTEGER h1,h2
80 C
81 C LOCAL VARIABLES
82 C
83  INTEGER nhel(nexternal,ncomb),ntry
84  REAL*8 t
85  REAL*8 matrix_cdx_udx_h
86  INTEGER ihel,iden, i
87  INTEGER jc(nexternal)
88  LOGICAL goodhel(ncomb)
89  DATA ntry/0/
90  DATA goodhel/ncomb*.false./
91  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
92  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
93  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
94  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
95  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
96  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
97  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
98  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
99  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
100  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
101  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
102  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
103  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
104  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
105  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
106  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
107  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
108  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
109  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
110  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
111  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
112  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
113  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
114  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
115  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
116  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
117  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
118  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
119  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
120  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
121  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
122  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
123  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
124  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
125  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
126  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
127  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
128  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
129  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
130  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
131  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
132  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
133  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
134  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
135  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
136  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
137  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
138  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
139  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
140  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
141  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
142  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
143  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
144  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
145  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
146  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
147  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
148  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
149  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
150  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
151  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
152  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
153  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
154  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
155  DATA iden/36/
156 C ----------
157 C BEGIN CODE
158 C ----------
159  DO ihel=1,nexternal
160  jc(ihel) = +1
161  ENDDO
162  ans = 0d0
163  DO ihel=1,ncomb
164  t=matrix_cdx_udx_h(p ,h1,h2,nhel(1,ihel),jc(1))
165  ans=ans+t
166  ENDDO
167  ans=ans/dble(iden)
168  END
169 
170 
171  REAL*8 FUNCTION matrix_cdx_udx_h(P,H1,H2,NHEL,IC)
172 C
173 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
174 C By the MadGraph Development Team
175 C Please visit us at https://launchpad.net/madgraph5
176 C
177 C Returns amplitude squared summed/avg over colors
178 C for the point with external lines W(0:6,NEXTERNAL)
179 C
180 C Process: c d~ > u d~ h WEIGHTED=6
181 C * Decay: h > ta+ ta- WEIGHTED=2
182 C Process: c d~ > c s~ h WEIGHTED=6
183 C * Decay: h > ta+ ta- WEIGHTED=2
184 C
185  IMPLICIT NONE
186 C
187 C CONSTANTS
188 C
189  INTEGER ngraphs
190  parameter(ngraphs=1)
191  INTEGER nexternal
192  parameter(nexternal=6)
193  INTEGER nwavefuncs, ncolor
194  parameter(nwavefuncs=7, ncolor=1)
195  REAL*8 zero
196  parameter(zero=0d0)
197  COMPLEX*16 imag1
198  parameter(imag1=(0d0,1d0))
199 C
200 C ARGUMENTS
201 C
202  REAL*8 p(0:3,nexternal)
203  INTEGER nhel(nexternal), ic(nexternal)
204 C
205 C LOCAL VARIABLES
206 C
207  INTEGER i,j
208  COMPLEX*16 ztemp
209  REAL*8 denom(ncolor), cf(ncolor,ncolor)
210  COMPLEX*16 amp(ngraphs), jamp(ncolor)
211  COMPLEX*16 w(18,nwavefuncs)
212  COMPLEX*16 dum0,dum1
213  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
214 C
215 C GLOBAL VARIABLES
216 C
217  include 'coupl.inc'
218 C
219 C COLOR DATA
220 C
221  DATA denom(1)/1/
222  DATA (cf(i, 1),i= 1, 1) / 9/
223 C 1 T(2,1) T(3,4)
224 
225  INTEGER h1,h2
226  REAL*8 matrix
227  matrix_cdx_udx_h=0.d0
228  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
229  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
230 
231 C ----------
232 C BEGIN CODE
233 C ----------
234  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
235  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
236  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
237  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
238  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
239  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
240  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
241  CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,6))
242  CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
243 C Amplitude(s) for diagram number 1
244  CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
245  jamp(1)=-amp(1)
246 
247  matrix = 0.d0
248  DO i = 1, ncolor
249  ztemp = (0.d0,0.d0)
250  DO j = 1, ncolor
251  ztemp = ztemp + cf(j,i)*jamp(j)
252  ENDDO
253  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
254  ENDDO
255  matrix_cdx_udx_h=matrix
256  ENDIF
257  ENDIF
258  END
259 
260 C ----- begin subprocesses CDX->CDX with Higgs->tautau
261 
262  SUBROUTINE cdx_cdx_h(P,H1,H2,ANS)
263 C
264 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
265 C By the MadGraph Development Team
266 C Please visit us at https://launchpad.net/madgraph5
267 C
268 C MadGraph StandAlone Version
269 C
270 C Returns amplitude squared summed/avg over colors
271 C and helicities
272 C for the point in phase space P(0:3,NEXTERNAL)
273 C
274 C Process: c d~ > c d~ h WEIGHTED=6
275 C * Decay: h > ta+ ta- WEIGHTED=2
276 C
277  IMPLICIT NONE
278 C
279 C CONSTANTS
280 C
281  INTEGER nexternal
282  parameter(nexternal=6)
283  INTEGER ncomb
284  parameter( ncomb=64)
285 C
286 C ARGUMENTS
287 C
288  REAL*8 p(0:3,nexternal),ans
289  INTEGER h1,h2
290 C
291 C LOCAL VARIABLES
292 C
293  INTEGER nhel(nexternal,ncomb),ntry
294  REAL*8 t
295  REAL*8 matrix_cdx_cdx_h
296  INTEGER ihel,iden, i
297  INTEGER jc(nexternal)
298  LOGICAL goodhel(ncomb)
299  DATA ntry/0/
300  DATA goodhel/ncomb*.false./
301  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
302  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
303  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
304  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
305  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
306  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
307  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
308  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
309  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
310  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
311  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
312  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
313  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
314  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
315  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
316  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
317  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
318  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
319  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
320  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
321  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
322  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
323  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
324  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
325  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
326  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
327  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
328  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
329  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
330  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
331  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
332  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
333  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
334  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
335  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
336  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
337  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
338  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
339  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
340  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
341  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
342  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
343  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
344  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
345  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
346  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
347  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
348  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
349  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
350  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
351  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
352  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
353  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
354  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
355  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
356  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
357  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
358  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
359  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
360  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
361  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
362  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
363  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
364  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
365  DATA iden/36/
366 C ----------
367 C BEGIN CODE
368 C ----------
369  DO ihel=1,nexternal
370  jc(ihel) = +1
371  ENDDO
372  ans = 0d0
373  DO ihel=1,ncomb
374  t=matrix_cdx_cdx_h(p ,h1,h2,nhel(1,ihel),jc(1))
375  ans=ans+t
376  ENDDO
377  ans=ans/dble(iden)
378  END
379 
380 
381  REAL*8 FUNCTION matrix_cdx_cdx_h(P,H1,H2,NHEL,IC)
382 C
383 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
384 C By the MadGraph Development Team
385 C Please visit us at https://launchpad.net/madgraph5
386 C
387 C Returns amplitude squared summed/avg over colors
388 C for the point with external lines W(0:6,NEXTERNAL)
389 C
390 C Process: c d~ > c d~ h WEIGHTED=6
391 C * Decay: h > ta+ ta- WEIGHTED=2
392 C
393  IMPLICIT NONE
394 C
395 C CONSTANTS
396 C
397  INTEGER ngraphs
398  parameter(ngraphs=2)
399  INTEGER nexternal
400  parameter(nexternal=6)
401  INTEGER nwavefuncs, ncolor
402  parameter(nwavefuncs=7, ncolor=2)
403  REAL*8 zero
404  parameter(zero=0d0)
405  COMPLEX*16 imag1
406  parameter(imag1=(0d0,1d0))
407 C
408 C ARGUMENTS
409 C
410  REAL*8 p(0:3,nexternal)
411  INTEGER nhel(nexternal), ic(nexternal)
412 C
413 C LOCAL VARIABLES
414 C
415  INTEGER i,j
416  COMPLEX*16 ztemp
417  REAL*8 denom(ncolor), cf(ncolor,ncolor)
418  COMPLEX*16 amp(ngraphs), jamp(ncolor)
419  COMPLEX*16 w(18,nwavefuncs)
420  COMPLEX*16 dum0,dum1
421  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
422 C
423 C GLOBAL VARIABLES
424 C
425  include 'coupl.inc'
426 C
427 C COLOR DATA
428 C
429  DATA denom(1)/1/
430  DATA (cf(i, 1),i= 1, 2) / 9, 3/
431 C 1 T(2,1) T(3,4)
432  DATA denom(2)/1/
433  DATA (cf(i, 2),i= 1, 2) / 3, 9/
434 C 1 T(2,4) T(3,1)
435 
436  INTEGER h1,h2
437  REAL*8 matrix
438  matrix_cdx_cdx_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 oxxxxx(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 ixxxxx(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,2),gc_44,mw,ww,w(1,6))
453  CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,5))
454 C Amplitude(s) for diagram number 1
455  CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_72,amp(1))
456  CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
457  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
458 C Amplitude(s) for diagram number 2
459  CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
460  jamp(1)=-amp(1)
461  jamp(2)=+amp(2)
462 
463  matrix = 0.d0
464  DO i = 1, ncolor
465  ztemp = (0.d0,0.d0)
466  DO j = 1, ncolor
467  ztemp = ztemp + cf(j,i)*jamp(j)
468  ENDDO
469  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
470  ENDDO
471  matrix_cdx_cdx_h=matrix
472  ENDIF
473  ENDIF
474  END
475 
476 C ----- begin subprocesses CDX->USX with Higgs->tautau
477 
478  SUBROUTINE cdx_usx_h(P,H1,H2,ANS)
479 C
480 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
481 C By the MadGraph Development Team
482 C Please visit us at https://launchpad.net/madgraph5
483 C
484 C MadGraph StandAlone Version
485 C
486 C Returns amplitude squared summed/avg over colors
487 C and helicities
488 C for the point in phase space P(0:3,NEXTERNAL)
489 C
490 C Process: c d~ > u s~ h WEIGHTED=6
491 C * Decay: h > ta+ ta- WEIGHTED=2
492 C
493  IMPLICIT NONE
494 C
495 C CONSTANTS
496 C
497  INTEGER nexternal
498  parameter(nexternal=6)
499  INTEGER ncomb
500  parameter( ncomb=64)
501 C
502 C ARGUMENTS
503 C
504  REAL*8 p(0:3,nexternal),ans
505  INTEGER h1,h2
506 C
507 C LOCAL VARIABLES
508 C
509  INTEGER nhel(nexternal,ncomb),ntry
510  REAL*8 t
511  REAL*8 matrix_cdx_usx_h
512  INTEGER ihel,iden, i
513  INTEGER jc(nexternal)
514  LOGICAL goodhel(ncomb)
515  DATA ntry/0/
516  DATA goodhel/ncomb*.false./
517  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
518  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
519  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
520  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
521  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
522  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
523  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
524  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
525  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
526  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
527  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
528  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
529  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
530  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
531  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
532  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
533  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
534  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
535  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
536  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
537  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
538  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
539  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
540  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
541  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
542  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
543  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
544  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
545  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
546  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
547  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
548  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
549  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
550  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
551  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
552  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
553  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
554  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
555  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
556  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
557  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
558  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
559  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
560  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
561  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
562  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
563  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
564  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
565  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
566  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
567  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
568  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
569  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
570  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
571  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
572  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
573  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
574  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
575  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
576  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
577  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
578  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
579  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
580  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
581  DATA iden/36/
582 C ----------
583 C BEGIN CODE
584 C ----------
585  DO ihel=1,nexternal
586  jc(ihel) = +1
587  ENDDO
588  ans = 0d0
589  DO ihel=1,ncomb
590  t=matrix_cdx_usx_h(p ,h1,h2,nhel(1,ihel),jc(1))
591  ans=ans+t
592  ENDDO
593  ans=ans/dble(iden)
594  END
595 
596 
597  REAL*8 FUNCTION matrix_cdx_usx_h(P,H1,H2,NHEL,IC)
598 C
599 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
600 C By the MadGraph Development Team
601 C Please visit us at https://launchpad.net/madgraph5
602 C
603 C Returns amplitude squared summed/avg over colors
604 C for the point with external lines W(0:6,NEXTERNAL)
605 C
606 C Process: c d~ > u s~ h WEIGHTED=6
607 C * Decay: h > ta+ ta- WEIGHTED=2
608 C
609  IMPLICIT NONE
610 C
611 C CONSTANTS
612 C
613  INTEGER ngraphs
614  parameter(ngraphs=1)
615  INTEGER nexternal
616  parameter(nexternal=6)
617  INTEGER nwavefuncs, ncolor
618  parameter(nwavefuncs=7, ncolor=1)
619  REAL*8 zero
620  parameter(zero=0d0)
621  COMPLEX*16 imag1
622  parameter(imag1=(0d0,1d0))
623 C
624 C ARGUMENTS
625 C
626  REAL*8 p(0:3,nexternal)
627  INTEGER nhel(nexternal), ic(nexternal)
628 C
629 C LOCAL VARIABLES
630 C
631  INTEGER i,j
632  COMPLEX*16 ztemp
633  REAL*8 denom(ncolor), cf(ncolor,ncolor)
634  COMPLEX*16 amp(ngraphs), jamp(ncolor)
635  COMPLEX*16 w(18,nwavefuncs)
636  COMPLEX*16 dum0,dum1
637  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
638 C
639 C GLOBAL VARIABLES
640 C
641  include 'coupl.inc'
642 C
643 C COLOR DATA
644 C
645  DATA denom(1)/1/
646  DATA (cf(i, 1),i= 1, 1) / 9/
647 C 1 T(2,1) T(3,4)
648 
649  INTEGER h1,h2
650  REAL*8 matrix
651  matrix_cdx_usx_h=0.d0
652  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
653  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
654 
655 C ----------
656 C BEGIN CODE
657 C ----------
658  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
659  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
660  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
661  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
662  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
663  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
664  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
665  CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,6))
666  CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,2))
667 C Amplitude(s) for diagram number 1
668  CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
669  jamp(1)=-amp(1)
670 
671  matrix = 0.d0
672  DO i = 1, ncolor
673  ztemp = (0.d0,0.d0)
674  DO j = 1, ncolor
675  ztemp = ztemp + cf(j,i)*jamp(j)
676  ENDDO
677  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
678  ENDDO
679  matrix_cdx_usx_h=matrix
680  ENDIF
681  ENDIF
682  END
683 
684 C ----- begin subprocesses CDX->CSX with Higgs->tautau
685 C ----- THE SAME AS CDX->UDX WITH Higgs->tautau
686 
687 C -------------- NOW PROCESSIS WITHOUT HIGGS
688 
689 C ----- begin subprocesses CDX->UDX NO Higgs
690 
691  SUBROUTINE cdx_udx_noh(P,H1,H2,ANS)
692 C
693 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
694 C By the MadGraph Development Team
695 C Please visit us at https://launchpad.net/madgraph5
696 C
697 C MadGraph StandAlone Version
698 C
699 C Returns amplitude squared summed/avg over colors
700 C and helicities
701 C for the point in phase space P(0:3,NEXTERNAL)
702 C
703 C Process: c d~ > u d~ ta+ ta- / h QED=4
704 C Process: c d~ > c s~ ta+ ta- / h QED=4
705 C
706  IMPLICIT NONE
707 C
708 C CONSTANTS
709 C
710  INTEGER nexternal
711  parameter(nexternal=6)
712  INTEGER ncomb
713  parameter( ncomb=64)
714 C
715 C ARGUMENTS
716 C
717  REAL*8 p(0:3,nexternal),ans
718  INTEGER h1,h2
719 C
720 C LOCAL VARIABLES
721 C
722  INTEGER nhel(nexternal,ncomb),ntry
723  REAL*8 t
724  REAL*8 matrix_cdx_udx_noh
725  INTEGER ihel,iden, i
726  INTEGER jc(nexternal)
727  LOGICAL goodhel(ncomb)
728  DATA ntry/0/
729  DATA goodhel/ncomb*.false./
730  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
731  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
732  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
733  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
734  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
735  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
736  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
737  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
738  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
739  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
740  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
741  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
742  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
743  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
744  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
745  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
746  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
747  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
748  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
749  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
750  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
751  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
752  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
753  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
754  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
755  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
756  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
757  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
758  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
759  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
760  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
761  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
762  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
763  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
764  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
765  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
766  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
767  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
768  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
769  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
770  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
771  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
772  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
773  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
774  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
775  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
776  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
777  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
778  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
779  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
780  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
781  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
782  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
783  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
784  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
785  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
786  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
787  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
788  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
789  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
790  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
791  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
792  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
793  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
794  DATA iden/36/
795 C ----------
796 C BEGIN CODE
797 C ----------
798  DO ihel=1,nexternal
799  jc(ihel) = +1
800  ENDDO
801  ans = 0d0
802  DO ihel=1,ncomb
803  t=matrix_cdx_udx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
804  ans=ans+t
805  ENDDO
806  ans=ans/dble(iden)
807  END
808 
809 
810  REAL*8 FUNCTION matrix_cdx_udx_noh(P,H1,H2,NHEL,IC)
811 C
812 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
813 C By the MadGraph Development Team
814 C Please visit us at https://launchpad.net/madgraph5
815 C
816 C Returns amplitude squared summed/avg over colors
817 C for the point with external lines W(0:6,NEXTERNAL)
818 C
819 C Process: c d~ > u d~ ta+ ta- / h QED=4
820 C Process: c d~ > c s~ ta+ ta- / h QED=4
821 C
822  IMPLICIT NONE
823 C
824 C CONSTANTS
825 C
826  INTEGER ngraphs
827  parameter(ngraphs=11)
828  INTEGER nexternal
829  parameter(nexternal=6)
830  INTEGER nwavefuncs, ncolor
831  parameter(nwavefuncs=10, ncolor=1)
832  REAL*8 zero
833  parameter(zero=0d0)
834  COMPLEX*16 imag1
835  parameter(imag1=(0d0,1d0))
836 C
837 C ARGUMENTS
838 C
839  REAL*8 p(0:3,nexternal)
840  INTEGER nhel(nexternal), ic(nexternal)
841 C
842 C LOCAL VARIABLES
843 C
844  INTEGER i,j
845  COMPLEX*16 ztemp
846  REAL*8 denom(ncolor), cf(ncolor,ncolor)
847  COMPLEX*16 amp(ngraphs), jamp(ncolor)
848  COMPLEX*16 w(18,nwavefuncs)
849  COMPLEX*16 dum0,dum1
850  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
851 C
852 C GLOBAL VARIABLES
853 C
854  include 'coupl.inc'
855 C
856 C COLOR DATA
857 C
858  DATA denom(1)/1/
859  DATA (cf(i, 1),i= 1, 1) / 9/
860 C 1 T(2,1) T(3,4)
861 
862  INTEGER h1,h2
863  REAL*8 matrix
864  matrix_cdx_udx_noh=0.d0
865  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
866  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
867 
868 C ----------
869 C BEGIN CODE
870 C ----------
871  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
872  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
873  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
874  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
875  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
876  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
877  CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,7))
878  CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
879  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
880 C Amplitude(s) for diagram number 1
881  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
882  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
883 C Amplitude(s) for diagram number 2
884  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
885  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
886 C Amplitude(s) for diagram number 3
887  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
888  CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
889 C Amplitude(s) for diagram number 4
890  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
891  CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
892 C Amplitude(s) for diagram number 5
893  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
894 C Amplitude(s) for diagram number 6
895  CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
896 C Amplitude(s) for diagram number 7
897  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
898  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
899 C Amplitude(s) for diagram number 8
900  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
901  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
902 C Amplitude(s) for diagram number 9
903  CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_44,amp(9))
904 C Amplitude(s) for diagram number 10
905  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
906  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
907 C Amplitude(s) for diagram number 11
908  CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_44,amp(11))
909  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
910  $ -amp(9)-amp(10)-amp(11)
911 
912  matrix = 0.d0
913  DO i = 1, ncolor
914  ztemp = (0.d0,0.d0)
915  DO j = 1, ncolor
916  ztemp = ztemp + cf(j,i)*jamp(j)
917  ENDDO
918  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
919  ENDDO
920  matrix_cdx_udx_noh=matrix
921  ENDIF
922  ENDIF
923  END
924 
925 C ----- begin subprocesses CDX->CDX NO Higgs
926 
927  SUBROUTINE cdx_cdx_noh(P,H1,H2,ANS)
928 C
929 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
930 C By the MadGraph Development Team
931 C Please visit us at https://launchpad.net/madgraph5
932 C
933 C MadGraph StandAlone Version
934 C
935 C Returns amplitude squared summed/avg over colors
936 C and helicities
937 C for the point in phase space P(0:3,NEXTERNAL)
938 C
939 C Process: c d~ > c d~ ta+ ta- / h QED=4
940 C
941  IMPLICIT NONE
942 C
943 C CONSTANTS
944 C
945  INTEGER nexternal
946  parameter(nexternal=6)
947  INTEGER ncomb
948  parameter( ncomb=64)
949 C
950 C ARGUMENTS
951 C
952  REAL*8 p(0:3,nexternal),ans
953  INTEGER h1,h2
954 C
955 C LOCAL VARIABLES
956 C
957  INTEGER nhel(nexternal,ncomb),ntry
958  REAL*8 t
959  REAL*8 matrix_cdx_cdx_noh
960  INTEGER ihel,iden, i
961  INTEGER jc(nexternal)
962  LOGICAL goodhel(ncomb)
963  DATA ntry/0/
964  DATA goodhel/ncomb*.false./
965  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
966  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
967  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
968  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
969  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
970  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
971  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
972  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
973  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
974  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
975  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
976  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
977  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
978  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
979  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
980  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
981  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
982  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
983  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
984  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
985  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
986  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
987  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
988  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
989  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
990  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
991  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
992  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
993  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
994  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
995  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
996  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
997  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
998  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
999  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1000  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1001  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1002  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1003  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1004  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1005  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1006  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1007  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1008  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1009  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1010  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1011  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1012  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1013  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1014  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1015  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1016  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1017  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1018  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1019  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1020  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1021  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1022  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1023  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1024  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1025  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1026  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1027  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1028  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1029  DATA iden/36/
1030 C ----------
1031 C BEGIN CODE
1032 C ----------
1033  DO ihel=1,nexternal
1034  jc(ihel) = +1
1035  ENDDO
1036  ans = 0d0
1037  DO ihel=1,ncomb
1038  t=matrix_cdx_cdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1039  ans=ans+t
1040  ENDDO
1041  ans=ans/dble(iden)
1042  END
1043 
1044 
1045  REAL*8 FUNCTION matrix_cdx_cdx_noh(P,H1,H2,NHEL,IC)
1046 C
1047 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1048 C By the MadGraph Development Team
1049 C Please visit us at https://launchpad.net/madgraph5
1050 C
1051 C Returns amplitude squared summed/avg over colors
1052 C for the point with external lines W(0:6,NEXTERNAL)
1053 C
1054 C Process: c d~ > c d~ ta+ ta- / h QED=4
1055 C
1056  IMPLICIT NONE
1057 C
1058 C CONSTANTS
1059 C
1060  INTEGER ngraphs
1061  parameter(ngraphs=43)
1062  INTEGER nexternal
1063  parameter(nexternal=6)
1064  INTEGER nwavefuncs, ncolor
1065  parameter(nwavefuncs=15, ncolor=2)
1066  REAL*8 zero
1067  parameter(zero=0d0)
1068  COMPLEX*16 imag1
1069  parameter(imag1=(0d0,1d0))
1070 C
1071 C ARGUMENTS
1072 C
1073  REAL*8 p(0:3,nexternal)
1074  INTEGER nhel(nexternal), ic(nexternal)
1075 C
1076 C LOCAL VARIABLES
1077 C
1078  INTEGER i,j
1079  COMPLEX*16 ztemp
1080  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1081  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1082  COMPLEX*16 w(18,nwavefuncs)
1083  COMPLEX*16 dum0,dum1
1084  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1085 C
1086 C GLOBAL VARIABLES
1087 C
1088  include 'coupl.inc'
1089 C
1090 C COLOR DATA
1091 C
1092  DATA denom(1)/1/
1093  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1094 C 1 T(2,1) T(3,4)
1095  DATA denom(2)/1/
1096  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1097 C 1 T(2,4) T(3,1)
1098 
1099  INTEGER h1,h2
1100  REAL*8 matrix
1101  matrix_cdx_cdx_noh=0.d0
1102  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1103  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1104 
1105 C ----------
1106 C BEGIN CODE
1107 C ----------
1108  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1109  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1110  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1111  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1112  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1113  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1114  CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,7))
1115  CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1116  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1117 C Amplitude(s) for diagram number 1
1118  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1119  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1120 C Amplitude(s) for diagram number 2
1121  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1122  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1123 C Amplitude(s) for diagram number 3
1124  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1125  CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,11))
1126 C Amplitude(s) for diagram number 4
1127  CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_1,amp(4))
1128  CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,12))
1129 C Amplitude(s) for diagram number 5
1130  CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(5))
1131 C Amplitude(s) for diagram number 6
1132  CALL ffv2_3_0(w(1,4),w(1,11),w(1,10),gc_50,gc_58,amp(6))
1133 C Amplitude(s) for diagram number 7
1134  CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1135  CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,12))
1136  CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,11))
1137  CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1138 C Amplitude(s) for diagram number 8
1139  CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1140  CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1141 C Amplitude(s) for diagram number 9
1142  CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1143  CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,14))
1144 C Amplitude(s) for diagram number 10
1145  CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1146 C Amplitude(s) for diagram number 11
1147  CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1148  CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,13))
1149  CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1150 C Amplitude(s) for diagram number 12
1151  CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1152  CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1153 C Amplitude(s) for diagram number 13
1154  CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1155 C Amplitude(s) for diagram number 14
1156  CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1157 C Amplitude(s) for diagram number 15
1158  CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1159  CALL ffv1_1(w(1,2),w(1,12),gc_1,zero,zero,w(1,15))
1160 C Amplitude(s) for diagram number 16
1161  CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(16))
1162  CALL ffv1_2(w(1,4),w(1,12),gc_1,zero,zero,w(1,5))
1163 C Amplitude(s) for diagram number 17
1164  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(17))
1165 C Amplitude(s) for diagram number 18
1166  CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(18))
1167 C Amplitude(s) for diagram number 19
1168  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(19))
1169  CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1170  CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1171 C Amplitude(s) for diagram number 20
1172  CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(20))
1173  CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1174 C Amplitude(s) for diagram number 21
1175  CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_1,amp(21))
1176 C Amplitude(s) for diagram number 22
1177  CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(22))
1178 C Amplitude(s) for diagram number 23
1179  CALL ffv2_3_0(w(1,12),w(1,2),w(1,10),gc_50,gc_58,amp(23))
1180  CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,12))
1181 C Amplitude(s) for diagram number 24
1182  CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_1,amp(24))
1183  CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,15))
1184 C Amplitude(s) for diagram number 25
1185  CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_1,amp(25))
1186 C Amplitude(s) for diagram number 26
1187  CALL ffv2_3_0(w(1,4),w(1,12),w(1,10),gc_50,gc_58,amp(26))
1188 C Amplitude(s) for diagram number 27
1189  CALL ffv2_3_0(w(1,15),w(1,2),w(1,10),gc_50,gc_58,amp(27))
1190  CALL ffv1_2(w(1,1),w(1,11),gc_2,zero,zero,w(1,15))
1191 C Amplitude(s) for diagram number 28
1192  CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(28))
1193  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,12))
1194 C Amplitude(s) for diagram number 29
1195  CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_2,amp(29))
1196 C Amplitude(s) for diagram number 30
1197  CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(30))
1198  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,15))
1199 C Amplitude(s) for diagram number 31
1200  CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_2,amp(31))
1201  CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1202  CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1203 C Amplitude(s) for diagram number 32
1204  CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_2,amp(32))
1205 C Amplitude(s) for diagram number 33
1206  CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1207 C Amplitude(s) for diagram number 34
1208  CALL ffv2_5_0(w(1,4),w(1,3),w(1,10),gc_51,gc_58,amp(34))
1209 C Amplitude(s) for diagram number 35
1210  CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1211  CALL ffv2_5_2(w(1,1),w(1,14),gc_51,gc_58,zero,zero,w(1,11))
1212 C Amplitude(s) for diagram number 36
1213  CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_2,amp(36))
1214 C Amplitude(s) for diagram number 37
1215  CALL ffv2_5_0(w(1,12),w(1,3),w(1,14),gc_51,gc_58,amp(37))
1216 C Amplitude(s) for diagram number 38
1217  CALL ffv2_5_0(w(1,11),w(1,3),w(1,10),gc_51,gc_58,amp(38))
1218 C Amplitude(s) for diagram number 39
1219  CALL ffv2_5_0(w(1,15),w(1,3),w(1,14),gc_51,gc_58,amp(39))
1220  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,14))
1221 C Amplitude(s) for diagram number 40
1222  CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_1,amp(40))
1223 C Amplitude(s) for diagram number 41
1224  CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_44,amp(41))
1225 C Amplitude(s) for diagram number 42
1226  CALL ffv2_3_0(w(1,14),w(1,2),w(1,10),gc_50,gc_58,amp(42))
1227 C Amplitude(s) for diagram number 43
1228  CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_44,amp(43))
1229  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1230  $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1231  $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1232  $ -amp(40)-amp(41)-amp(42)-amp(43)
1233  jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1234  $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1235  $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1236  $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1237  $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1238  $ +amp(38)+amp(39)
1239 
1240  matrix = 0.d0
1241  DO i = 1, ncolor
1242  ztemp = (0.d0,0.d0)
1243  DO j = 1, ncolor
1244  ztemp = ztemp + cf(j,i)*jamp(j)
1245  ENDDO
1246  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1247  ENDDO
1248  matrix_cdx_cdx_noh=matrix
1249  ENDIF
1250  ENDIF
1251  END
1252 
1253 C ----- begin subprocesses CDX->USX NO Higgs
1254 
1255  SUBROUTINE cdx_usx_noh(P,H1,H2,ANS)
1256 C
1257 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1258 C By the MadGraph Development Team
1259 C Please visit us at https://launchpad.net/madgraph5
1260 C
1261 C MadGraph StandAlone Version
1262 C
1263 C Returns amplitude squared summed/avg over colors
1264 C and helicities
1265 C for the point in phase space P(0:3,NEXTERNAL)
1266 C
1267 C Process: c d~ > u s~ ta+ ta- / h QED=4
1268 C
1269  IMPLICIT NONE
1270 C
1271 C CONSTANTS
1272 C
1273  INTEGER nexternal
1274  parameter(nexternal=6)
1275  INTEGER ncomb
1276  parameter( ncomb=64)
1277 C
1278 C ARGUMENTS
1279 C
1280  REAL*8 p(0:3,nexternal),ans
1281  INTEGER h1,h2
1282 C
1283 C LOCAL VARIABLES
1284 C
1285  INTEGER nhel(nexternal,ncomb),ntry
1286  REAL*8 t
1287  REAL*8 matrix_cdx_usx_noh
1288  INTEGER ihel,iden, i
1289  INTEGER jc(nexternal)
1290  LOGICAL goodhel(ncomb)
1291  DATA ntry/0/
1292  DATA goodhel/ncomb*.false./
1293  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1294  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1295  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1296  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1297  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1298  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1299  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1300  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1301  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1302  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1303  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1304  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1305  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1306  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1307  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1308  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1309  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1310  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1311  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1312  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1313  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1314  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1315  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1316  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1317  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1318  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1319  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1320  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1321  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1322  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1323  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1324  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1325  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1326  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1327  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1328  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1329  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1330  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1331  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1332  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1333  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1334  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1335  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1336  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1337  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1338  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1339  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1340  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1341  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1342  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1343  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1344  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1345  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1346  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1347  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1348  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1349  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1350  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1351  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1352  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1353  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1354  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1355  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1356  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1357  DATA iden/36/
1358 C ----------
1359 C BEGIN CODE
1360 C ----------
1361  DO ihel=1,nexternal
1362  jc(ihel) = +1
1363  ENDDO
1364  ans = 0d0
1365  DO ihel=1,ncomb
1366  t=matrix_cdx_usx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1367  ans=ans+t
1368  ENDDO
1369  ans=ans/dble(iden)
1370  END
1371 
1372 
1373  REAL*8 FUNCTION matrix_cdx_usx_noh(P,H1,H2,NHEL,IC)
1374 C
1375 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1376 C By the MadGraph Development Team
1377 C Please visit us at https://launchpad.net/madgraph5
1378 C
1379 C Returns amplitude squared summed/avg over colors
1380 C for the point with external lines W(0:6,NEXTERNAL)
1381 C
1382 C Process: c d~ > u s~ ta+ ta- / h QED=4
1383 C
1384  IMPLICIT NONE
1385 C
1386 C CONSTANTS
1387 C
1388  INTEGER ngraphs
1389  parameter(ngraphs=11)
1390  INTEGER nexternal
1391  parameter(nexternal=6)
1392  INTEGER nwavefuncs, ncolor
1393  parameter(nwavefuncs=10, ncolor=1)
1394  REAL*8 zero
1395  parameter(zero=0d0)
1396  COMPLEX*16 imag1
1397  parameter(imag1=(0d0,1d0))
1398 C
1399 C ARGUMENTS
1400 C
1401  REAL*8 p(0:3,nexternal)
1402  INTEGER nhel(nexternal), ic(nexternal)
1403 C
1404 C LOCAL VARIABLES
1405 C
1406  INTEGER i,j
1407  COMPLEX*16 ztemp
1408  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1409  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1410  COMPLEX*16 w(18,nwavefuncs)
1411  COMPLEX*16 dum0,dum1
1412  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1413 C
1414 C GLOBAL VARIABLES
1415 C
1416  include 'coupl.inc'
1417 C
1418 C COLOR DATA
1419 C
1420  DATA denom(1)/1/
1421  DATA (cf(i, 1),i= 1, 1) / 9/
1422 C 1 T(2,1) T(3,4)
1423 
1424  INTEGER h1,h2
1425  REAL*8 matrix
1426  matrix_cdx_usx_noh=0.d0
1427  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1428  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1429 
1430 C ----------
1431 C BEGIN CODE
1432 C ----------
1433  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1434  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1435  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1436  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1437  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1438  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1439  CALL ffv2_3(w(1,1),w(1,2),gc_44,mw,ww,w(1,7))
1440  CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1441  CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1442 C Amplitude(s) for diagram number 1
1443  CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1444  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1445 C Amplitude(s) for diagram number 2
1446  CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1447  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1448 C Amplitude(s) for diagram number 3
1449  CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1450  CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,6))
1451 C Amplitude(s) for diagram number 4
1452  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1453  CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
1454 C Amplitude(s) for diagram number 5
1455  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1456 C Amplitude(s) for diagram number 6
1457  CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1458 C Amplitude(s) for diagram number 7
1459  CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1460  CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1461 C Amplitude(s) for diagram number 8
1462  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1463  CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1464 C Amplitude(s) for diagram number 9
1465  CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_44,amp(9))
1466 C Amplitude(s) for diagram number 10
1467  CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1468  CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1469 C Amplitude(s) for diagram number 11
1470  CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_44,amp(11))
1471  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1472  $ -amp(9)-amp(10)-amp(11)
1473 
1474  matrix = 0.d0
1475  DO i = 1, ncolor
1476  ztemp = (0.d0,0.d0)
1477  DO j = 1, ncolor
1478  ztemp = ztemp + cf(j,i)*jamp(j)
1479  ENDDO
1480  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1481  ENDDO
1482  matrix_cdx_usx_noh=matrix
1483  ENDIF
1484  ENDIF
1485  END
1486 
1487 C ----- begin subprocesses CDX->CSX NO Higgs
1488 C IS THE SAME AS CDX->UDX NO Higgs
1489 
1490 
1491 
1492 C DONE 26.10.2015 STILL S=2 MISSING
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7