C++ Interface to Tauola
DCX.f
1c --------------------- 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
9C I3=0,1 OR 3 I4=0,-2 OR -4
10
11 ans=0.d0
12
13C 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
47C ----- begin subprocesses DCX->DUX with Higgs->tautau
48
49
50 SUBROUTINE dcx_dux_h(P,H1,H2,ANS)
51C
52C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
53C By the MadGraph Development Team
54C Please visit us at https://launchpad.net/madgraph5
55C
56C MadGraph StandAlone Version
57C
58C Returns amplitude squared summed/avg over colors
59C and helicities
60C for the point in phase space P(0:3,NEXTERNAL)
61C
62C Process: d c~ > d u~ h WEIGHTED=6
63C * Decay: h > ta+ ta- WEIGHTED=2
64C Process: d c~ > s c~ h WEIGHTED=6
65C * Decay: h > ta+ ta- WEIGHTED=2
66C
67 IMPLICIT NONE
68C
69C CONSTANTS
70C
71 INTEGER NEXTERNAL
72 parameter(nexternal=6)
73 INTEGER NCOMB
74 parameter( ncomb=64)
75C
76C ARGUMENTS
77C
78 real*8 p(0:3,nexternal),ans
79 INTEGER H1,H2
80C
81C LOCAL VARIABLES
82C
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/
156C ----------
157C BEGIN CODE
158C ----------
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)
172C
173C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
174C By the MadGraph Development Team
175C Please visit us at https://launchpad.net/madgraph5
176C
177C Returns amplitude squared summed/avg over colors
178C for the point with external lines W(0:6,NEXTERNAL)
179C
180C Process: d c~ > d u~ h WEIGHTED=6
181C * Decay: h > ta+ ta- WEIGHTED=2
182C Process: d c~ > s c~ h WEIGHTED=6
183C * Decay: h > ta+ ta- WEIGHTED=2
184C
185 IMPLICIT NONE
186C
187C CONSTANTS
188C
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))
199C
200C ARGUMENTS
201C
202 real*8 p(0:3,nexternal)
203 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
204C
205C LOCAL VARIABLES
206C
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)/
214C
215C GLOBAL VARIABLES
216C
217 include 'coupl.inc'
218C
219C COLOR DATA
220C
221 DATA denom(1)/1/
222 DATA (cf(i, 1),i= 1, 1) / 9/
223C 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
231C ----------
232C BEGIN CODE
233C ----------
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))
243C 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
261C ----- begin subprocesses DCX->DCX with Higgs->tautau
262
263 SUBROUTINE dcx_dcx_h(P,H1,H2,ANS)
264C
265C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
266C By the MadGraph Development Team
267C Please visit us at https://launchpad.net/madgraph5
268C
269C MadGraph StandAlone Version
270C
271C Returns amplitude squared summed/avg over colors
272C and helicities
273C for the point in phase space P(0:3,NEXTERNAL)
274C
275C Process: d c~ > d c~ h WEIGHTED=6
276C * Decay: h > ta+ ta- WEIGHTED=2
277C
278 IMPLICIT NONE
279C
280C CONSTANTS
281C
282 INTEGER NEXTERNAL
283 parameter(nexternal=6)
284 INTEGER NCOMB
285 parameter( ncomb=64)
286C
287C ARGUMENTS
288C
289 real*8 p(0:3,nexternal),ans
290 INTEGER H1,H2
291C
292C LOCAL VARIABLES
293C
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/
367C ----------
368C BEGIN CODE
369C ----------
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)
383C
384C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
385C By the MadGraph Development Team
386C Please visit us at https://launchpad.net/madgraph5
387C
388C Returns amplitude squared summed/avg over colors
389C for the point with external lines W(0:6,NEXTERNAL)
390C
391C Process: d c~ > d c~ h WEIGHTED=6
392C * Decay: h > ta+ ta- WEIGHTED=2
393C
394 IMPLICIT NONE
395C
396C CONSTANTS
397C
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))
408C
409C ARGUMENTS
410C
411 real*8 p(0:3,nexternal)
412 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
413C
414C LOCAL VARIABLES
415C
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)/
423C
424C GLOBAL VARIABLES
425C
426 include 'coupl.inc'
427C
428C COLOR DATA
429C
430 DATA denom(1)/1/
431 DATA (cf(i, 1),i= 1, 2) / 9, 3/
432C 1 T(2,1) T(3,4)
433 DATA denom(2)/1/
434 DATA (cf(i, 2),i= 1, 2) / 3, 9/
435C 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
443C ----------
444C BEGIN CODE
445C ----------
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))
455C 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))
459C 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
477C ----- begin subprocesses DCX->SUX with Higgs->tautau
478
479 SUBROUTINE dcx_sux_h(P,H1,H2,ANS)
480C
481C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
482C By the MadGraph Development Team
483C Please visit us at https://launchpad.net/madgraph5
484C
485C MadGraph StandAlone Version
486C
487C Returns amplitude squared summed/avg over colors
488C and helicities
489C for the point in phase space P(0:3,NEXTERNAL)
490C
491C Process: d c~ > s u~ h WEIGHTED=6
492C * Decay: h > ta+ ta- WEIGHTED=2
493C
494 IMPLICIT NONE
495C
496C CONSTANTS
497C
498 INTEGER NEXTERNAL
499 parameter(nexternal=6)
500 INTEGER NCOMB
501 parameter( ncomb=64)
502C
503C ARGUMENTS
504C
505 real*8 p(0:3,nexternal),ans
506 INTEGER H1,H2
507C
508C LOCAL VARIABLES
509C
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/
583C ----------
584C BEGIN CODE
585C ----------
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)
599C
600C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
601C By the MadGraph Development Team
602C Please visit us at https://launchpad.net/madgraph5
603C
604C Returns amplitude squared summed/avg over colors
605C for the point with external lines W(0:6,NEXTERNAL)
606C
607C Process: d c~ > s u~ h WEIGHTED=6
608C * Decay: h > ta+ ta- WEIGHTED=2
609C
610 IMPLICIT NONE
611C
612C CONSTANTS
613C
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))
624C
625C ARGUMENTS
626C
627 real*8 p(0:3,nexternal)
628 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
629C
630C LOCAL VARIABLES
631C
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)/
639C
640C GLOBAL VARIABLES
641C
642 include 'coupl.inc'
643C
644C COLOR DATA
645C
646 DATA denom(1)/1/
647 DATA (cf(i, 1),i= 1, 1) / 9/
648C 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
656C ----------
657C BEGIN CODE
658C ----------
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))
668C 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
685C ----- begin subprocesses DCX->SCX with Higgs->tautau
686C ----- THE SAME AS DCX->DUX WITH Higgs->tautau
687
688C -------------- NOW PROCESSIS WITHOUT HIGGS
689
690C ----- begin subprocesses DCX->DUX NO Higgs
691
692 SUBROUTINE dcx_dux_noh(P,H1,H2,ANS)
693C
694C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
695C By the MadGraph Development Team
696C Please visit us at https://launchpad.net/madgraph5
697C
698C MadGraph StandAlone Version
699C
700C Returns amplitude squared summed/avg over colors
701C and helicities
702C for the point in phase space P(0:3,NEXTERNAL)
703C
704C Process: d c~ > d u~ ta+ ta- / h QED=4
705C Process: d c~ > s c~ ta+ ta- / h QED=4
706C
707 IMPLICIT NONE
708C
709C CONSTANTS
710C
711 INTEGER NEXTERNAL
712 parameter(nexternal=6)
713 INTEGER NCOMB
714 parameter( ncomb=64)
715C
716C ARGUMENTS
717C
718 real*8 p(0:3,nexternal),ans
719 INTEGER H1,H2
720C
721C LOCAL VARIABLES
722C
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/
796C ----------
797C BEGIN CODE
798C ----------
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)
812C
813C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
814C By the MadGraph Development Team
815C Please visit us at https://launchpad.net/madgraph5
816C
817C Returns amplitude squared summed/avg over colors
818C for the point with external lines W(0:6,NEXTERNAL)
819C
820C Process: d c~ > d u~ ta+ ta- / h QED=4
821C Process: d c~ > s c~ ta+ ta- / h QED=4
822C
823 IMPLICIT NONE
824C
825C CONSTANTS
826C
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))
837C
838C ARGUMENTS
839C
840 real*8 p(0:3,nexternal)
841 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
842C
843C LOCAL VARIABLES
844C
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)/
852C
853C GLOBAL VARIABLES
854C
855 include 'coupl.inc'
856C
857C COLOR DATA
858C
859 DATA denom(1)/1/
860 DATA (cf(i, 1),i= 1, 1) / 9/
861C 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
869C ----------
870C BEGIN CODE
871C ----------
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))
881C 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))
884C 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))
887C 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))
890C 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))
893C Amplitude(s) for diagram number 5
894 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
895C 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))
897C 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))
900C 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))
903C Amplitude(s) for diagram number 9
904 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_44,amp(9))
905C 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))
908C 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
927C ----- begin subprocesses DCX->DCX NO Higgs
928
929 SUBROUTINE dcx_dcx_noh(P,H1,H2,ANS)
930C
931C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
932C By the MadGraph Development Team
933C Please visit us at https://launchpad.net/madgraph5
934C
935C MadGraph StandAlone Version
936C
937C Returns amplitude squared summed/avg over colors
938C and helicities
939C for the point in phase space P(0:3,NEXTERNAL)
940C
941C Process: d c~ > d c~ ta+ ta- / h QED=4
942C
943 IMPLICIT NONE
944C
945C CONSTANTS
946C
947 INTEGER NEXTERNAL
948 parameter(nexternal=6)
949 INTEGER NCOMB
950 parameter( ncomb=64)
951C
952C ARGUMENTS
953C
954 real*8 p(0:3,nexternal),ans
955 INTEGER H1,H2
956C
957C LOCAL VARIABLES
958C
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/
1032C ----------
1033C BEGIN CODE
1034C ----------
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)
1048C
1049C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1050C By the MadGraph Development Team
1051C Please visit us at https://launchpad.net/madgraph5
1052C
1053C Returns amplitude squared summed/avg over colors
1054C for the point with external lines W(0:6,NEXTERNAL)
1055C
1056C Process: d c~ > d c~ ta+ ta- / h QED=4
1057C
1058 IMPLICIT NONE
1059C
1060C CONSTANTS
1061C
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))
1072C
1073C ARGUMENTS
1074C
1075 real*8 p(0:3,nexternal)
1076 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1077C
1078C LOCAL VARIABLES
1079C
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)/
1087C
1088C GLOBAL VARIABLES
1089C
1090 include 'coupl.inc'
1091C
1092C COLOR DATA
1093C
1094 DATA denom(1)/1/
1095 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1096C 1 T(2,1) T(3,4)
1097 DATA denom(2)/1/
1098 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1099C 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
1107C ----------
1108C BEGIN CODE
1109C ----------
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))
1119C 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))
1122C 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))
1125C 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))
1128C 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))
1131C Amplitude(s) for diagram number 5
1132 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_1,amp(5))
1133C 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))
1135C 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))
1140C 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))
1143C 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))
1146C 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))
1148C 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))
1152C 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))
1155C Amplitude(s) for diagram number 13
1156 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1157C 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))
1159C 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))
1162C 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))
1165C Amplitude(s) for diagram number 17
1166 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(17))
1167C 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))
1169C 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))
1173C 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))
1176C Amplitude(s) for diagram number 21
1177 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_2,amp(21))
1178C 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))
1180C 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))
1183C 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))
1186C Amplitude(s) for diagram number 25
1187 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_2,amp(25))
1188C 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))
1190C 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))
1193C 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))
1196C Amplitude(s) for diagram number 29
1197 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_1,amp(29))
1198C 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))
1201C 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))
1205C Amplitude(s) for diagram number 32
1206 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_1,amp(32))
1207C Amplitude(s) for diagram number 33
1208 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1209C 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))
1211C 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))
1214C Amplitude(s) for diagram number 36
1215 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_1,amp(36))
1216C 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))
1218C 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))
1220C 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))
1223C Amplitude(s) for diagram number 40
1224 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_2,amp(40))
1225C Amplitude(s) for diagram number 41
1226 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_44,amp(41))
1227C 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))
1229C 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
1255C ----- begin subprocesses DCX->SUX NO Higgs
1256
1257 SUBROUTINE dcx_sux_noh(P,H1,H2,ANS)
1258C
1259C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1260C By the MadGraph Development Team
1261C Please visit us at https://launchpad.net/madgraph5
1262C
1263C MadGraph StandAlone Version
1264C
1265C Returns amplitude squared summed/avg over colors
1266C and helicities
1267C for the point in phase space P(0:3,NEXTERNAL)
1268C
1269C Process: d c~ > s u~ ta+ ta- / h QED=4
1270C
1271 IMPLICIT NONE
1272C
1273C CONSTANTS
1274C
1275 INTEGER NEXTERNAL
1276 parameter(nexternal=6)
1277 INTEGER NCOMB
1278 parameter( ncomb=64)
1279C
1280C ARGUMENTS
1281C
1282 real*8 p(0:3,nexternal),ans
1283 INTEGER H1,H2
1284C
1285C LOCAL VARIABLES
1286C
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/
1360C ----------
1361C BEGIN CODE
1362C ----------
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)
1376C
1377C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1378C By the MadGraph Development Team
1379C Please visit us at https://launchpad.net/madgraph5
1380C
1381C Returns amplitude squared summed/avg over colors
1382C for the point with external lines W(0:6,NEXTERNAL)
1383C
1384C Process: d c~ > s u~ ta+ ta- / h QED=4
1385C
1386 IMPLICIT NONE
1387C
1388C CONSTANTS
1389C
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))
1400C
1401C ARGUMENTS
1402C
1403 real*8 p(0:3,nexternal)
1404 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1405C
1406C LOCAL VARIABLES
1407C
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)/
1415C
1416C GLOBAL VARIABLES
1417C
1418 include 'coupl.inc'
1419C
1420C COLOR DATA
1421C
1422 DATA denom(1)/1/
1423 DATA (cf(i, 1),i= 1, 1) / 9/
1424C 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
1432C ----------
1433C BEGIN CODE
1434C ----------
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))
1444C 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))
1447C 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))
1450C 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))
1453C 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))
1456C Amplitude(s) for diagram number 5
1457 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1458C 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))
1460C 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))
1463C 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))
1466C Amplitude(s) for diagram number 9
1467 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_44,amp(9))
1468C 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))
1471C 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
1489C ----- begin subprocesses DCX->SCX NO Higgs
1490C IS THE SAME AS DCX->DUX NO Higgs
1491
1492
1493
1494C DONE 26.10.2015 STILL S=2 MISSING