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