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