C++ Interface to Tauola
CD.f
1
2c --------------------- subprocesses initiated by CD: 13.08.2015
3
4 SUBROUTINE cd(P,I3,I4,H1,H2,KEY,ANS)
5 IMPLICIT NONE
6
7 INTEGER I3,I4, H1,H2, KEY
8 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
9
10C I3=0,2 OR 4 I4=0,1 OR 3
11C
12 ans=0.0d0
13
14 IF(key.EQ.1) THEN
15 IF(i3.EQ.2 .AND. i4.EQ.1) CALL cd_ud_h(p,h1,h2,ans)
16 IF(i3.EQ.2 .AND. i4.EQ.3) CALL cd_us_h(p,h1,h2,ans)
17 IF(i3.EQ.4 .AND. i4.EQ.1) CAll cd_cd_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.3) CALL cd_cs_h(p,h1,h2,ans)
19 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
20 CALL cd_ud_h(p,h1,h2,ans1)
21 CALL cd_us_h(p,h1,h2,ans2)
22 CALL cd_cd_h(p,h1,h2,ans3)
23 CALL cd_cs_h(p,h1,h2,ans4)
24 ans=ans1+ans2+ans3+ans4
25 ENDIF
26 ELSE IF(key.EQ.0) THEN
27 IF(i3.EQ.2 .AND. i4.EQ.1) CALL cd_ud_noh(p,h1,h2,ans)
28 IF(i3.EQ.2 .AND. i4.EQ.3) CALL cd_us_noh(p,h1,h2,ans)
29 IF(i3.EQ.4 .AND. i4.EQ.1) CAll cd_cd_noh(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.3) CALL cd_cs_noh(p,h1,h2,ans)
31 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
32 CALL cd_ud_noh(p,h1,h2,ans1)
33 CALL cd_us_noh(p,h1,h2,ans2)
34 CALL cd_cd_noh(p,h1,h2,ans3)
35 CALL cd_cs_noh(p,h1,h2,ans4)
36 ans=ans1+ans2+ans3+ans4
37 ENDIF
38 ELSE
39 WRITE(*,*) "NOT FINISHED"
40 stop
41 ENDIF
42 END ! SUBROUTINE CD
43
44
45
46C ----- begin subprocesses CD->jjH with Higgs->tautau
47
48 SUBROUTINE cd_cd_h(P,H1,H2,ANS)
49C
50C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
51C By the MadGraph Development Team
52C Please visit us at https://launchpad.net/madgraph5
53C
54C MadGraph StandAlone Version
55C
56C Returns amplitude squared summed/avg over colors
57C and helicities
58C for the point in phase space P(0:3,NEXTERNAL)
59C
60C Process: c d > c d h WEIGHTED=6
61C * Decay: h > ta+ ta- WEIGHTED=2
62C
63 IMPLICIT NONE
64C
65C CONSTANTS
66C
67 INTEGER NEXTERNAL
68 parameter(nexternal=6)
69 INTEGER NCOMB
70 parameter( ncomb=64)
71C
72C ARGUMENTS
73C
74 real*8 p(0:3,nexternal),ans
75 INTEGER H1,H2
76C
77C LOCAL VARIABLES
78C
79 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
80 real*8 t
81 real*8 matrix_cd_cd_h
82 INTEGER IHEL,IDEN, I
83 INTEGER JC(NEXTERNAL)
84 LOGICAL GOODHEL(NCOMB)
85 DATA ntry/0/
86 DATA goodhel/ncomb*.false./
87 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
88 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
89 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
90 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
91 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
92 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
93 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
94 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
95 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
96 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
97 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
98 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
99 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
100 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
101 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
102 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
103 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
104 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
105 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
106 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
107 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
108 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
109 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
110 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
111 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
112 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
113 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
114 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
115 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
116 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
117 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
118 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
119 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
120 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
121 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
122 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
123 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
124 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
125 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
126 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
127 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
128 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
129 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
130 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
131 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
132 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
133 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
134 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
135 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
136 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
137 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
138 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
139 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
140 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
141 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
142 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
143 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
144 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
145 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
146 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
147 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
148 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
149 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
150 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
151 DATA iden/36/
152C ----------
153C BEGIN CODE
154C ----------
155 DO ihel=1,nexternal
156 jc(ihel) = +1
157 ENDDO
158 ans = 0d0
159 DO ihel=1,ncomb
160 t=matrix_cd_cd_h(p ,h1,h2,nhel(1,ihel),jc(1))
161 ans=ans+t
162 ENDDO
163 ans=ans/dble(iden)
164 END
165
166
167 REAL*8 FUNCTION matrix_cd_cd_h(P,H1,H2,NHEL,IC)
168C
169C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
170C By the MadGraph Development Team
171C Please visit us at https://launchpad.net/madgraph5
172C
173C Returns amplitude squared summed/avg over colors
174C for the point with external lines W(0:6,NEXTERNAL)
175C
176C Process: c d > c d h WEIGHTED=6
177C * Decay: h > ta+ ta- WEIGHTED=2
178C
179 IMPLICIT NONE
180C
181C CONSTANTS
182C
183 INTEGER NGRAPHS
184 parameter(ngraphs=2)
185 INTEGER NEXTERNAL
186 parameter(nexternal=6)
187 INTEGER NWAVEFUNCS, NCOLOR
188 parameter(nwavefuncs=7, ncolor=2)
189 real*8 zero
190 parameter(zero=0d0)
191 COMPLEX*16 IMAG1
192 parameter(imag1=(0d0,1d0))
193C
194C ARGUMENTS
195C
196 real*8 p(0:3,nexternal)
197 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
198C
199C LOCAL VARIABLES
200C
201 INTEGER I,J
202 COMPLEX*16 ZTEMP
203 real*8 denom(ncolor), cf(ncolor,ncolor)
204 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
205 COMPLEX*16 W(18,NWAVEFUNCS)
206 COMPLEX*16 DUM0,DUM1
207 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
208C
209C GLOBAL VARIABLES
210C
211 include 'coupl.inc'
212C
213C COLOR DATA
214C
215 DATA denom(1)/1/
216 DATA (cf(i, 1),i= 1, 2) / 9, 3/
217C 1 T(3,1) T(4,2)
218 DATA denom(2)/1/
219 DATA (cf(i, 2),i= 1, 2) / 3, 9/
220C 1 T(3,2) T(4,1)
221
222 INTEGER H1,H2
223 real*8 matrix
224 matrix_cd_cd_h=0.d0
225 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
226 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
227
228
229
230C ----------
231C BEGIN CODE
232C ----------
233 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
234 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
235 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
236 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
237 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
238 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
239 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
240 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,6))
241 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
242C Amplitude(s) for diagram number 1
243 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
244 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,5))
245 CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,4))
246C Amplitude(s) for diagram number 2
247 CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_72,amp(2))
248 jamp(1)=-amp(1)
249 jamp(2)=+amp(2)
250
251 matrix = 0.d0
252 DO i = 1, ncolor
253 ztemp = (0.d0,0.d0)
254 DO j = 1, ncolor
255 ztemp = ztemp + cf(j,i)*jamp(j)
256 ENDDO
257 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
258 ENDDO
259 matrix_cd_cd_h=matrix
260 ENDIF ! CLOSES H1 IF
261 ENDIF ! CLOSES H2 IF
262 END
263
264c==========================================
265
266 SUBROUTINE cd_cs_h(P,H1,H2,ANS)
267C
268C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
269C By the MadGraph Development Team
270C Please visit us at https://launchpad.net/madgraph5
271C
272C MadGraph StandAlone Version
273C
274C Returns amplitude squared summed/avg over colors
275C and helicities
276C for the point in phase space P(0:3,NEXTERNAL)
277C
278C Process: c d > c s h WEIGHTED=6
279C * Decay: h > ta+ ta- WEIGHTED=2
280C
281 IMPLICIT NONE
282C
283C CONSTANTS
284C
285 INTEGER NEXTERNAL
286 parameter(nexternal=6)
287 INTEGER NCOMB
288 parameter( ncomb=64)
289C
290C ARGUMENTS
291C
292 real*8 p(0:3,nexternal),ans
293 INTEGER H1,H2
294C
295C LOCAL VARIABLES
296C
297 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
298 real*8 t
299 real*8 matrix_cd_cs_h,matrix
300 INTEGER IHEL,IDEN, I
301 INTEGER JC(NEXTERNAL)
302 LOGICAL GOODHEL(NCOMB)
303 DATA ntry/0/
304 DATA goodhel/ncomb*.false./
305 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
306 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
307 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
308 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
309 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
310 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
311 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
312 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
313 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
314 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
315 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
316 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
317 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
318 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
319 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
320 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
321 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
322 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
323 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
324 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
325 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
326 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
327 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
328 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
329 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
330 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
331 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
332 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
333 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
334 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
335 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
336 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
337 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
338 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
339 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
340 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
341 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
342 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
343 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
344 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
345 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
346 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
347 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
348 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
349 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
350 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
351 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
352 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
353 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
354 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
355 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
356 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
357 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
358 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
359 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
360 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
361 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
362 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
363 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
364 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
365 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
366 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
367 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
368 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
369 DATA iden/36/
370C ----------
371C BEGIN CODE
372C ----------
373 DO ihel=1,nexternal
374 jc(ihel) = +1
375 ENDDO
376 ans = 0d0
377 DO ihel=1,ncomb
378 t=matrix_cd_cs_h(p,h1,h2 ,nhel(1,ihel),jc(1))
379 ans=ans+t
380 ENDDO
381 ans=ans/dble(iden)
382 END
383
384
385 REAL*8 FUNCTION matrix_cd_cs_h(P,H1,H2,NHEL,IC)
386C
387C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
388C By the MadGraph Development Team
389C Please visit us at https://launchpad.net/madgraph5
390C
391C Returns amplitude squared summed/avg over colors
392C for the point with external lines W(0:6,NEXTERNAL)
393C
394C Process: c d > c s h WEIGHTED=6
395C * Decay: h > ta+ ta- WEIGHTED=2
396C
397 IMPLICIT NONE
398C
399C CONSTANTS
400C
401 INTEGER NGRAPHS
402 parameter(ngraphs=1)
403 INTEGER NEXTERNAL
404 parameter(nexternal=6)
405 INTEGER NWAVEFUNCS, NCOLOR
406 parameter(nwavefuncs=7, ncolor=1)
407 real*8 zero
408 parameter(zero=0d0)
409 COMPLEX*16 IMAG1
410 parameter(imag1=(0d0,1d0))
411C
412C ARGUMENTS
413C
414 real*8 p(0:3,nexternal)
415 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
416C
417C LOCAL VARIABLES
418C
419 INTEGER I,J
420 COMPLEX*16 ZTEMP
421 real*8 denom(ncolor), cf(ncolor,ncolor)
422 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
423 COMPLEX*16 W(18,NWAVEFUNCS)
424 COMPLEX*16 DUM0,DUM1
425 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
426C
427C GLOBAL VARIABLES
428C
429 include 'coupl.inc'
430C
431C COLOR DATA
432C
433 DATA denom(1)/1/
434 DATA (cf(i, 1),i= 1, 1) / 9/
435C 1 T(3,2) T(4,1)
436
437 INTEGER H1,H2
438 real*8 matrix
439 matrix_cd_cs_h=0.d0
440 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
441 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
442
443C ----------
444C BEGIN CODE
445C ----------
446 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
447 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
448 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
449 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
450 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
451 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
452 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
453 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,6))
454 CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,4))
455C Amplitude(s) for diagram number 1
456 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
457 jamp(1)=+amp(1)
458
459 matrix = 0.d0
460 DO i = 1, ncolor
461 ztemp = (0.d0,0.d0)
462 DO j = 1, ncolor
463 ztemp = ztemp + cf(j,i)*jamp(j)
464 ENDDO
465 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
466 ENDDO
467 matrix_cd_cs_h=matrix
468 ENDIF
469 ENDIF
470 END
471
472
473
474C===============================================
475
476 SUBROUTINE cd_ud_h(P,H1,H2,ANS)
477C
478C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
479C By the MadGraph Development Team
480C Please visit us at https://launchpad.net/madgraph5
481C
482C MadGraph StandAlone Version
483C
484C Returns amplitude squared summed/avg over colors
485C and helicities
486C for the point in phase space P(0:3,NEXTERNAL)
487C
488C Process: c d > u d h WEIGHTED=6
489C * Decay: h > ta+ ta- WEIGHTED=2
490C
491 IMPLICIT NONE
492C
493C CONSTANTS
494C
495 INTEGER NEXTERNAL
496 parameter(nexternal=6)
497 INTEGER NCOMB
498 parameter( ncomb=64)
499C
500C ARGUMENTS
501C
502 real*8 p(0:3,nexternal),ans
503 INTEGER H1,H2
504C
505C LOCAL VARIABLES
506C
507 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
508 real*8 t
509 real*8 matrix_cd_ud_h
510 INTEGER IHEL,IDEN, I
511 INTEGER JC(NEXTERNAL)
512 LOGICAL GOODHEL(NCOMB)
513 DATA ntry/0/
514 DATA goodhel/ncomb*.false./
515 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
516 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
517 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
518 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
519 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
520 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
521 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
522 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
523 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
524 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
525 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
526 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
527 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
528 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
529 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
530 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
531 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
532 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
533 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
534 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
535 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
536 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
537 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
538 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
539 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
540 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
541 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
542 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
543 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
544 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
545 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
546 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
547 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
548 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
549 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
550 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
551 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
552 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
553 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
554 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
555 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
556 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
557 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
558 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
559 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
560 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
561 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
562 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
563 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
564 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
565 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
566 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
567 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
568 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
569 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
570 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
571 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
572 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
573 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
574 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
575 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
576 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
577 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
578 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
579 DATA iden/36/
580C ----------
581C BEGIN CODE
582C ----------
583 DO ihel=1,nexternal
584 jc(ihel) = +1
585 ENDDO
586 ans = 0d0
587 DO ihel=1,ncomb
588 t=matrix_cd_ud_h(p,h1,h2 ,nhel(1,ihel),jc(1))
589 ans=ans+t
590 ENDDO
591 ans=ans/dble(iden)
592 END
593
594
595 REAL*8 FUNCTION matrix_cd_ud_h(P,H1,H2,NHEL,IC)
596C
597C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
598C By the MadGraph Development Team
599C Please visit us at https://launchpad.net/madgraph5
600C
601C Returns amplitude squared summed/avg over colors
602C for the point with external lines W(0:6,NEXTERNAL)
603C
604C Process: c d > u d h WEIGHTED=6
605C * Decay: h > ta+ ta- WEIGHTED=2
606C
607 IMPLICIT NONE
608C
609C CONSTANTS
610C
611 INTEGER NGRAPHS
612 parameter(ngraphs=1)
613 INTEGER NEXTERNAL
614 parameter(nexternal=6)
615 INTEGER NWAVEFUNCS, NCOLOR
616 parameter(nwavefuncs=7, ncolor=1)
617 real*8 zero
618 parameter(zero=0d0)
619 COMPLEX*16 IMAG1
620 parameter(imag1=(0d0,1d0))
621C
622C ARGUMENTS
623C
624 real*8 p(0:3,nexternal)
625 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
626C
627C LOCAL VARIABLES
628C
629 INTEGER I,J
630 COMPLEX*16 ZTEMP
631 real*8 denom(ncolor), cf(ncolor,ncolor)
632 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
633 COMPLEX*16 W(18,NWAVEFUNCS)
634 COMPLEX*16 DUM0,DUM1
635 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
636C
637C GLOBAL VARIABLES
638C
639 include 'coupl.inc'
640C
641C COLOR DATA
642C
643 DATA denom(1)/1/
644 DATA (cf(i, 1),i= 1, 1) / 9/
645C 1 T(3,2) T(4,1)
646
647 INTEGER H1,H2
648 real*8 matrix
649 matrix_cd_ud_h=0d0
650 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
651 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
652
653C ----------
654C BEGIN CODE
655C ----------
656 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
657 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
658 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
659 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
660 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
661 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
662 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
663 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,6))
664 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,4))
665C Amplitude(s) for diagram number 1
666 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
667 jamp(1)=+amp(1)
668
669 matrix = 0.d0
670 DO i = 1, ncolor
671 ztemp = (0.d0,0.d0)
672 DO j = 1, ncolor
673 ztemp = ztemp + cf(j,i)*jamp(j)
674 ENDDO
675 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
676 ENDDO
677 matrix_cd_ud_h=matrix
678 ENDIF
679 ENDIF
680 END
681
682
683
684C========================================
685
686 SUBROUTINE cd_us_h(P,H1,H2,ANS)
687C
688C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
689C By the MadGraph Development Team
690C Please visit us at https://launchpad.net/madgraph5
691C
692C MadGraph StandAlone Version
693C
694C Returns amplitude squared summed/avg over colors
695C and helicities
696C for the point in phase space P(0:3,NEXTERNAL)
697C
698C Process: c d > u s h WEIGHTED=6
699C * Decay: h > ta+ ta- WEIGHTED=2
700C
701 IMPLICIT NONE
702C
703C CONSTANTS
704C
705 INTEGER NEXTERNAL
706 parameter(nexternal=6)
707 INTEGER NCOMB
708 parameter( ncomb=64)
709C
710C ARGUMENTS
711C
712 real*8 p(0:3,nexternal),ans
713 INTEGER H1,H2
714C
715C LOCAL VARIABLES
716C
717 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
718 real*8 t
719 real*8 matrix_cd_us_h
720 INTEGER IHEL,IDEN, I
721 INTEGER JC(NEXTERNAL)
722 LOGICAL GOODHEL(NCOMB)
723 DATA ntry/0/
724 DATA goodhel/ncomb*.false./
725 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
726 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
727 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
728 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
729 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
730 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
731 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
732 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
733 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
734 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
735 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
736 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
737 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
738 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
739 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
740 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
741 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
742 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
743 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
744 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
745 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
746 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
747 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
748 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
749 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
750 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
751 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
752 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
753 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
754 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
755 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
756 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
757 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
758 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
759 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
760 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
761 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
762 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
763 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
764 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
765 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
766 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
767 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
768 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
769 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
770 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
771 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
772 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
773 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
774 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
775 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
776 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
777 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
778 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
779 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
780 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
781 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
782 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
783 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
784 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
785 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
786 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
787 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
788 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
789 DATA iden/36/
790C ----------
791C BEGIN CODE
792C ----------
793 DO ihel=1,nexternal
794 jc(ihel) = +1
795 ENDDO
796 ans = 0d0
797 DO ihel=1,ncomb
798 t=matrix_cd_us_h(p,h1,h2 ,nhel(1,ihel),jc(1))
799 ans=ans+t
800 ENDDO
801 ans=ans/dble(iden)
802 END
803
804
805 REAL*8 FUNCTION matrix_cd_us_h(P,H1,H2,NHEL,IC)
806C
807C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
808C By the MadGraph Development Team
809C Please visit us at https://launchpad.net/madgraph5
810C
811C Returns amplitude squared summed/avg over colors
812C for the point with external lines W(0:6,NEXTERNAL)
813C
814C Process: c d > u s h WEIGHTED=6
815C * Decay: h > ta+ ta- WEIGHTED=2
816C
817 IMPLICIT NONE
818C
819C CONSTANTS
820C
821 INTEGER NGRAPHS
822 parameter(ngraphs=1)
823 INTEGER NEXTERNAL
824 parameter(nexternal=6)
825 INTEGER NWAVEFUNCS, NCOLOR
826 parameter(nwavefuncs=7, ncolor=1)
827 real*8 zero
828 parameter(zero=0d0)
829 COMPLEX*16 IMAG1
830 parameter(imag1=(0d0,1d0))
831C
832C ARGUMENTS
833C
834 real*8 p(0:3,nexternal)
835 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
836C
837C LOCAL VARIABLES
838C
839 INTEGER I,J
840 COMPLEX*16 ZTEMP
841 real*8 denom(ncolor), cf(ncolor,ncolor)
842 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
843 COMPLEX*16 W(18,NWAVEFUNCS)
844 COMPLEX*16 DUM0,DUM1
845 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
846C
847C GLOBAL VARIABLES
848C
849 include 'coupl.inc'
850C
851C COLOR DATA
852C
853 DATA denom(1)/1/
854 DATA (cf(i, 1),i= 1, 1) / 9/
855C 1 T(3,2) T(4,1)
856
857 INTEGER H1,H2
858 real*8 matrix
859 matrix_cd_us_h=0d0
860 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
861 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
862
863C ----------
864C BEGIN CODE
865C ----------
866 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
867 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
868 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
869 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
870 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
871 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
872 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
873 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,6))
874 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,4))
875C Amplitude(s) for diagram number 1
876 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_72,amp(1))
877 jamp(1)=+amp(1)
878
879 matrix = 0.d0
880 DO i = 1, ncolor
881 ztemp = (0.d0,0.d0)
882 DO j = 1, ncolor
883 ztemp = ztemp + cf(j,i)*jamp(j)
884 ENDDO
885 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
886 ENDDO
887 matrix_cd_us_h=matrix
888 ENDIF
889 ENDIF
890 END
891
892
893
894
895C --- end of subprocess CD->jjH, H->tautau
896
897C ----- begin subprocesses CD->jjtautau, no Higgs
898
899 SUBROUTINE cd_cd_noh(P,H1,H2,ANS)
900C
901C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
902C By the MadGraph Development Team
903C Please visit us at https://launchpad.net/madgraph5
904C
905C MadGraph StandAlone Version
906C
907C Returns amplitude squared summed/avg over colors
908C and helicities
909C for the point in phase space P(0:3,NEXTERNAL)
910C
911C Process: c d > c d ta+ ta- / h QED=4
912C
913 IMPLICIT NONE
914C
915C CONSTANTS
916C
917 INTEGER NEXTERNAL
918 parameter(nexternal=6)
919 INTEGER NCOMB
920 parameter( ncomb=64)
921C
922C ARGUMENTS
923C
924 real*8 p(0:3,nexternal),ans
925 INTEGER H1,H2
926C
927C LOCAL VARIABLES
928C
929 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
930 real*8 t
931 real*8 matrix_cd_cd_noh
932 INTEGER IHEL,IDEN, I
933 INTEGER JC(NEXTERNAL)
934 LOGICAL GOODHEL(NCOMB)
935 DATA ntry/0/
936 DATA goodhel/ncomb*.false./
937 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
938 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
939 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
940 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
941 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
942 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
943 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
944 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
945 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
946 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
947 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
948 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
949 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
950 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
951 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
952 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
953 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
954 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
955 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
956 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
957 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
958 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
959 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
960 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
961 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
962 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
963 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
964 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
965 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
966 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
967 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
968 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
969 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
970 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
971 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
972 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
973 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
974 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
975 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
976 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
977 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
978 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
979 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
980 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
981 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
982 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
983 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
984 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
985 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
986 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
987 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
988 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
989 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
990 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
991 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
992 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
993 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
994 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
995 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
996 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
997 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
998 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
999 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1000 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1001 DATA iden/36/
1002C ----------
1003C BEGIN CODE
1004C ----------
1005 DO ihel=1,nexternal
1006 jc(ihel) = +1
1007 ENDDO
1008 ans = 0d0
1009 DO ihel=1,ncomb
1010 t=matrix_cd_cd_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1011 ans=ans+t
1012 ENDDO
1013 ans=ans/dble(iden)
1014 END
1015
1016
1017 REAL*8 FUNCTION matrix_cd_cd_noh(P,H1,H2,NHEL,IC)
1018C
1019C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1020C By the MadGraph Development Team
1021C Please visit us at https://launchpad.net/madgraph5
1022C
1023C Returns amplitude squared summed/avg over colors
1024C for the point with external lines W(0:6,NEXTERNAL)
1025C
1026C Process: c d > c d ta+ ta- / h QED=4
1027C
1028 IMPLICIT NONE
1029C
1030C CONSTANTS
1031C
1032 INTEGER NGRAPHS
1033 parameter(ngraphs=43)
1034 INTEGER NEXTERNAL
1035 parameter(nexternal=6)
1036 INTEGER NWAVEFUNCS, NCOLOR
1037 parameter(nwavefuncs=14, ncolor=2)
1038 real*8 zero
1039 parameter(zero=0d0)
1040 COMPLEX*16 IMAG1
1041 parameter(imag1=(0d0,1d0))
1042C
1043C ARGUMENTS
1044C
1045 real*8 p(0:3,nexternal)
1046 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1047C
1048C LOCAL VARIABLES
1049C
1050 INTEGER I,J
1051 COMPLEX*16 ZTEMP
1052 real*8 denom(ncolor), cf(ncolor,ncolor)
1053 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1054 COMPLEX*16 W(18,NWAVEFUNCS)
1055 COMPLEX*16 DUM0,DUM1
1056 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1057C
1058C GLOBAL VARIABLES
1059C
1060 include 'coupl.inc'
1061C
1062C COLOR DATA
1063C
1064 DATA denom(1)/1/
1065 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1066C 1 T(3,1) T(4,2)
1067 DATA denom(2)/1/
1068 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1069C 1 T(3,2) T(4,1)
1070
1071 INTEGER H1,H2
1072 real*8 matrix
1073 matrix_cd_cd_noh=0d0
1074 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1075 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1076
1077
1078C ----------
1079C BEGIN CODE
1080C ----------
1081 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1082 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1083 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1084 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1085 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1086 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1087 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,7))
1088 CALL ffv1p0_3(w(1,2),w(1,4),gc_1,zero,zero,w(1,8))
1089 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1090C Amplitude(s) for diagram number 1
1091 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1092 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1093C Amplitude(s) for diagram number 2
1094 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1095 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,11))
1096C Amplitude(s) for diagram number 3
1097 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1098C Amplitude(s) for diagram number 4
1099 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1100 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,10))
1101 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1102C Amplitude(s) for diagram number 5
1103 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1104 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1105C Amplitude(s) for diagram number 6
1106 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1107C Amplitude(s) for diagram number 7
1108 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1109C Amplitude(s) for diagram number 8
1110 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1111 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1112 CALL ffv1_2(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
1113C Amplitude(s) for diagram number 9
1114 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(9))
1115 CALL ffv1_1(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1116C Amplitude(s) for diagram number 10
1117 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_1,amp(10))
1118 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1119C Amplitude(s) for diagram number 11
1120 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(11))
1121C Amplitude(s) for diagram number 12
1122 CALL ffv2_3_0(w(1,2),w(1,13),w(1,7),gc_50,gc_58,amp(12))
1123 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
1124 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
1125C Amplitude(s) for diagram number 13
1126 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(13))
1127 CALL ffv1_1(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1128C Amplitude(s) for diagram number 14
1129 CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_1,amp(14))
1130C Amplitude(s) for diagram number 15
1131 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(15))
1132C Amplitude(s) for diagram number 16
1133 CALL ffv2_3_0(w(1,2),w(1,14),w(1,7),gc_50,gc_58,amp(16))
1134 CALL ffv2_3_2(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
1135C Amplitude(s) for diagram number 17
1136 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_1,amp(17))
1137 CALL ffv2_3_1(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1138C Amplitude(s) for diagram number 18
1139 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
1140C Amplitude(s) for diagram number 19
1141 CALL ffv2_3_0(w(1,14),w(1,4),w(1,7),gc_50,gc_58,amp(19))
1142C Amplitude(s) for diagram number 20
1143 CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(20))
1144 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,9))
1145 CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,14))
1146 CALL ffv2_2(w(1,5),w(1,9),gc_108,zero,zero,w(1,10))
1147C Amplitude(s) for diagram number 21
1148 CALL ffv2_0(w(1,10),w(1,6),w(1,14),gc_108,amp(21))
1149C Amplitude(s) for diagram number 22
1150 CALL vvv1_0(w(1,12),w(1,9),w(1,14),gc_4,amp(22))
1151C Amplitude(s) for diagram number 23
1152 CALL vvv1_0(w(1,9),w(1,14),w(1,7),gc_53,amp(23))
1153 CALL ffv2_2(w(1,2),w(1,9),gc_44,zero,zero,w(1,10))
1154C Amplitude(s) for diagram number 24
1155 CALL ffv1_0(w(1,10),w(1,3),w(1,12),gc_2,amp(24))
1156 CALL ffv2_1(w(1,3),w(1,9),gc_44,zero,zero,w(1,6))
1157C Amplitude(s) for diagram number 25
1158 CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(25))
1159C Amplitude(s) for diagram number 26
1160 CALL ffv2_5_0(w(1,10),w(1,3),w(1,7),gc_51,gc_58,amp(26))
1161C Amplitude(s) for diagram number 27
1162 CALL ffv2_3_0(w(1,2),w(1,6),w(1,7),gc_50,gc_58,amp(27))
1163 CALL ffv2_2(w(1,1),w(1,14),gc_44,zero,zero,w(1,6))
1164C Amplitude(s) for diagram number 28
1165 CALL ffv1_0(w(1,6),w(1,4),w(1,12),gc_1,amp(28))
1166 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,10))
1167C Amplitude(s) for diagram number 29
1168 CALL ffv2_0(w(1,10),w(1,4),w(1,14),gc_44,amp(29))
1169C Amplitude(s) for diagram number 30
1170 CALL ffv2_3_0(w(1,6),w(1,4),w(1,7),gc_50,gc_58,amp(30))
1171 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,6))
1172C Amplitude(s) for diagram number 31
1173 CALL ffv2_0(w(1,6),w(1,4),w(1,14),gc_44,amp(31))
1174 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,14))
1175C Amplitude(s) for diagram number 32
1176 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(32))
1177C Amplitude(s) for diagram number 33
1178 CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_2,amp(33))
1179C Amplitude(s) for diagram number 34
1180 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(34))
1181C Amplitude(s) for diagram number 35
1182 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_2,amp(35))
1183 CALL ffv1p0_3(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
1184 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1185C Amplitude(s) for diagram number 36
1186 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(36))
1187C Amplitude(s) for diagram number 37
1188 CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_11,amp(37))
1189C Amplitude(s) for diagram number 38
1190 CALL ffv2_5_0(w(1,4),w(1,3),w(1,7),gc_51,gc_58,amp(38))
1191C Amplitude(s) for diagram number 39
1192 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(39))
1193 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
1194C Amplitude(s) for diagram number 40
1195 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(40))
1196C Amplitude(s) for diagram number 41
1197 CALL ffv2_5_0(w(1,10),w(1,3),w(1,11),gc_51,gc_58,amp(41))
1198C Amplitude(s) for diagram number 42
1199 CALL ffv2_5_0(w(1,8),w(1,3),w(1,7),gc_51,gc_58,amp(42))
1200C Amplitude(s) for diagram number 43
1201 CALL ffv2_5_0(w(1,6),w(1,3),w(1,11),gc_51,gc_58,amp(43))
1202 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1203 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1204 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1205 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
1206 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
1207 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1208 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1209 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
1210 $ -1./2.*amp(38)-1./2.*amp(39)
1211
1212 matrix = 0.d0
1213 DO i = 1, ncolor
1214 ztemp = (0.d0,0.d0)
1215 DO j = 1, ncolor
1216 ztemp = ztemp + cf(j,i)*jamp(j)
1217 ENDDO
1218 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1219 ENDDO
1220 matrix_cd_cd_noh=matrix
1221 ENDIF ! CLOSES H1 IF
1222 ENDIF ! CLOSES H2 IF
1223 END
1224
1225
1226 SUBROUTINE cd_cs_noh(P,H1,H2,ANS)
1227C
1228C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1229C By the MadGraph Development Team
1230C Please visit us at https://launchpad.net/madgraph5
1231C
1232C MadGraph StandAlone Version
1233C
1234C Returns amplitude squared summed/avg over colors
1235C and helicities
1236C for the point in phase space P(0:3,NEXTERNAL)
1237C
1238C Process: c d > c s ta+ ta- / h QED=4
1239C
1240 IMPLICIT NONE
1241C
1242C CONSTANTS
1243C
1244 INTEGER NEXTERNAL
1245 parameter(nexternal=6)
1246 INTEGER NCOMB
1247 parameter( ncomb=64)
1248C
1249C ARGUMENTS
1250C
1251 real*8 p(0:3,nexternal),ans
1252 INTEGER H1,H2
1253C
1254C LOCAL VARIABLES
1255C
1256 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1257 real*8 t
1258 real*8 matrix_cd_cs_noh
1259 INTEGER IHEL,IDEN, I
1260 INTEGER JC(NEXTERNAL)
1261 LOGICAL GOODHEL(NCOMB)
1262 DATA ntry/0/
1263 DATA goodhel/ncomb*.false./
1264 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1265 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1266 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1267 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1268 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1269 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1270 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1271 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1272 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1273 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1274 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1275 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1276 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1277 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1278 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1279 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1280 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1281 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1282 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1283 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1284 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1285 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1286 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1287 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1288 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1289 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1290 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1291 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1292 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1293 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1294 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1295 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1296 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1297 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1298 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1299 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1300 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1301 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1302 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1303 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1304 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1305 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1306 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1307 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1308 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1309 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1310 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1311 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1312 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1313 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1314 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1315 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1316 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1317 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1318 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1319 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1320 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1321 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1322 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1323 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1324 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1325 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1326 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1327 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1328 DATA iden/36/
1329C ----------
1330C BEGIN CODE
1331C ----------
1332 DO ihel=1,nexternal
1333 jc(ihel) = +1
1334 ENDDO
1335 ans = 0d0
1336 DO ihel=1,ncomb
1337 t=matrix_cd_cs_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1338 ans=ans+t
1339 ENDDO
1340 ans=ans/dble(iden)
1341 END
1342
1343
1344 REAL*8 FUNCTION matrix_cd_cs_noh(P,H1,H2,NHEL,IC)
1345C
1346C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1347C By the MadGraph Development Team
1348C Please visit us at https://launchpad.net/madgraph5
1349C
1350C Returns amplitude squared summed/avg over colors
1351C for the point with external lines W(0:6,NEXTERNAL)
1352C
1353C Process: c d > c s ta+ ta- / h QED=4
1354C
1355 IMPLICIT NONE
1356C
1357C CONSTANTS
1358C
1359 INTEGER NGRAPHS
1360 parameter(ngraphs=11)
1361 INTEGER NEXTERNAL
1362 parameter(nexternal=6)
1363 INTEGER NWAVEFUNCS, NCOLOR
1364 parameter(nwavefuncs=10, ncolor=1)
1365 real*8 zero
1366 parameter(zero=0d0)
1367 COMPLEX*16 IMAG1
1368 parameter(imag1=(0d0,1d0))
1369C
1370C ARGUMENTS
1371C
1372 real*8 p(0:3,nexternal)
1373 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1374C
1375C LOCAL VARIABLES
1376C
1377 INTEGER I,J
1378 COMPLEX*16 ZTEMP
1379 real*8 denom(ncolor), cf(ncolor,ncolor)
1380 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1381 COMPLEX*16 W(18,NWAVEFUNCS)
1382 COMPLEX*16 DUM0,DUM1
1383 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1384C
1385C GLOBAL VARIABLES
1386C
1387 include 'coupl.inc'
1388C
1389C COLOR DATA
1390C
1391 DATA denom(1)/1/
1392 DATA (cf(i, 1),i= 1, 1) / 9/
1393C 1 T(3,2) T(4,1)
1394
1395 INTEGER H1,H2
1396 real*8 matrix
1397 matrix_cd_cs_noh=0d0
1398 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1399 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1400
1401
1402C ----------
1403C BEGIN CODE
1404C ----------
1405 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1406 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1407 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1408 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1409 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1410 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1411 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,7))
1412 CALL ffv2_3(w(1,2),w(1,3),gc_44,mw,ww,w(1,8))
1413 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1414C Amplitude(s) for diagram number 1
1415 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1416 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1417C Amplitude(s) for diagram number 2
1418 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1419 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1420C Amplitude(s) for diagram number 3
1421 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1422 CALL ffv2_2(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
1423C Amplitude(s) for diagram number 4
1424 CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1425 CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,5))
1426C Amplitude(s) for diagram number 5
1427 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1428C Amplitude(s) for diagram number 6
1429 CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1430C Amplitude(s) for diagram number 7
1431 CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1432 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1433C Amplitude(s) for diagram number 8
1434 CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1435 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1436C Amplitude(s) for diagram number 9
1437 CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_100,amp(9))
1438C Amplitude(s) for diagram number 10
1439 CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1440 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1441C Amplitude(s) for diagram number 11
1442 CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_100,amp(11))
1443 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1444 $ +amp(9)+amp(10)+amp(11)
1445
1446 matrix = 0.d0
1447 DO i = 1, ncolor
1448 ztemp = (0.d0,0.d0)
1449 DO j = 1, ncolor
1450 ztemp = ztemp + cf(j,i)*jamp(j)
1451 ENDDO
1452 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1453 ENDDO
1454 matrix_cd_cs_noh=matrix
1455 ENDIF ! CLOSES H1 IF
1456 ENDIF ! CLOSES H2 IF
1457 END
1458
1459 SUBROUTINE cd_ud_noh(P,H1,H2,ANS)
1460C
1461C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1462C By the MadGraph Development Team
1463C Please visit us at https://launchpad.net/madgraph5
1464C
1465C MadGraph StandAlone Version
1466C
1467C Returns amplitude squared summed/avg over colors
1468C and helicities
1469C for the point in phase space P(0:3,NEXTERNAL)
1470C
1471C Process: c d > u d ta+ ta- / h QED=4
1472C
1473 IMPLICIT NONE
1474C
1475C CONSTANTS
1476C
1477 INTEGER NEXTERNAL
1478 parameter(nexternal=6)
1479 INTEGER NCOMB
1480 parameter( ncomb=64)
1481C
1482C ARGUMENTS
1483C
1484 real*8 p(0:3,nexternal),ans
1485 INTEGER H1,H2
1486C
1487C LOCAL VARIABLES
1488C
1489 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1490 real*8 t
1491 real*8 matrix_cd_ud_noh
1492 INTEGER IHEL,IDEN, I
1493 INTEGER JC(NEXTERNAL)
1494 LOGICAL GOODHEL(NCOMB)
1495 DATA ntry/0/
1496 DATA goodhel/ncomb*.false./
1497 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1498 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1499 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1500 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1501 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1502 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1503 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1504 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1505 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1506 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1507 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1508 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1509 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1510 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1511 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1512 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1513 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1514 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1515 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1516 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1517 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1518 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1519 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1520 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1521 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1522 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1523 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1524 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1525 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1526 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1527 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1528 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1529 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1530 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1531 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1532 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1533 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1534 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1535 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1536 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1537 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1538 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1539 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1540 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1541 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1542 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1543 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1544 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1545 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1546 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1547 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1548 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1549 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1550 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1551 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1552 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1553 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1554 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1555 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1556 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1557 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1558 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1559 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1560 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1561 DATA iden/36/
1562C ----------
1563C BEGIN CODE
1564C ----------
1565 DO ihel=1,nexternal
1566 jc(ihel) = +1
1567 ENDDO
1568 ans = 0d0
1569 DO ihel=1,ncomb
1570 t=matrix_cd_ud_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1571 ans=ans+t
1572 ENDDO
1573 ans=ans/dble(iden)
1574 END
1575
1576
1577 REAL*8 FUNCTION matrix_cd_ud_noh(P,H1,H2,NHEL,IC)
1578C
1579C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1580C By the MadGraph Development Team
1581C Please visit us at https://launchpad.net/madgraph5
1582C
1583C Returns amplitude squared summed/avg over colors
1584C for the point with external lines W(0:6,NEXTERNAL)
1585C
1586C Process: c d > u d ta+ ta- / h QED=4
1587C
1588 IMPLICIT NONE
1589C
1590C CONSTANTS
1591C
1592 INTEGER NGRAPHS
1593 parameter(ngraphs=11)
1594 INTEGER NEXTERNAL
1595 parameter(nexternal=6)
1596 INTEGER NWAVEFUNCS, NCOLOR
1597 parameter(nwavefuncs=10, ncolor=1)
1598 real*8 zero
1599 parameter(zero=0d0)
1600 COMPLEX*16 IMAG1
1601 parameter(imag1=(0d0,1d0))
1602C
1603C ARGUMENTS
1604C
1605 real*8 p(0:3,nexternal)
1606 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1607C
1608C LOCAL VARIABLES
1609C
1610 INTEGER I,J
1611 COMPLEX*16 ZTEMP
1612 real*8 denom(ncolor), cf(ncolor,ncolor)
1613 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1614 COMPLEX*16 W(18,NWAVEFUNCS)
1615 COMPLEX*16 DUM0,DUM1
1616 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1617C
1618C GLOBAL VARIABLES
1619C
1620 include 'coupl.inc'
1621C
1622C COLOR DATA
1623C
1624 DATA denom(1)/1/
1625 DATA (cf(i, 1),i= 1, 1) / 9/
1626C 1 T(3,2) T(4,1)
1627
1628 INTEGER H1,H2
1629 real*8 matrix
1630 matrix_cd_ud_noh=0d0
1631 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1632 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1633
1634C ----------
1635C BEGIN CODE
1636C ----------
1637 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1638 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1639 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1640 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1641 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1642 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1643 CALL ffv2_3(w(1,1),w(1,4),gc_44,mw,ww,w(1,7))
1644 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,8))
1645 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1646C Amplitude(s) for diagram number 1
1647 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1648 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1649C Amplitude(s) for diagram number 2
1650 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1651 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1652C Amplitude(s) for diagram number 3
1653 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1654 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1655C Amplitude(s) for diagram number 4
1656 CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1657 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,5))
1658C Amplitude(s) for diagram number 5
1659 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1660C Amplitude(s) for diagram number 6
1661 CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1662C Amplitude(s) for diagram number 7
1663 CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1664 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1665C Amplitude(s) for diagram number 8
1666 CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1667 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1668C Amplitude(s) for diagram number 9
1669 CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_44,amp(9))
1670C Amplitude(s) for diagram number 10
1671 CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1672 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1673C Amplitude(s) for diagram number 11
1674 CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_44,amp(11))
1675 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1676 $ +amp(9)+amp(10)+amp(11)
1677
1678 matrix = 0.d0
1679 DO i = 1, ncolor
1680 ztemp = (0.d0,0.d0)
1681 DO j = 1, ncolor
1682 ztemp = ztemp + cf(j,i)*jamp(j)
1683 ENDDO
1684 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1685 ENDDO
1686 matrix_cd_ud_noh=matrix
1687 ENDIF ! CLOSES H1 IF
1688 ENDIF ! CLOSES H2 IF
1689 END
1690
1691
1692 SUBROUTINE cd_us_noh(P,H1,H2,ANS)
1693C
1694C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1695C By the MadGraph Development Team
1696C Please visit us at https://launchpad.net/madgraph5
1697C
1698C MadGraph StandAlone Version
1699C
1700C Returns amplitude squared summed/avg over colors
1701C and helicities
1702C for the point in phase space P(0:3,NEXTERNAL)
1703C
1704C Process: c d > u s ta+ ta- / h QED=4
1705C
1706 IMPLICIT NONE
1707C
1708C CONSTANTS
1709C
1710 INTEGER NEXTERNAL
1711 parameter(nexternal=6)
1712 INTEGER NCOMB
1713 parameter( ncomb=64)
1714C
1715C ARGUMENTS
1716C
1717 real*8 p(0:3,nexternal),ans
1718 INTEGER H1,H2
1719C
1720C LOCAL VARIABLES
1721C
1722 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1723 real*8 t
1724 real*8 matrix_cd_us_noh
1725 INTEGER IHEL,IDEN, I
1726 INTEGER JC(NEXTERNAL)
1727 LOGICAL GOODHEL(NCOMB)
1728 DATA ntry/0/
1729 DATA goodhel/ncomb*.false./
1730 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1731 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1732 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1733 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1734 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1735 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1736 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1737 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1738 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1739 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1740 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1741 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1742 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1743 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1744 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1745 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1746 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1747 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1748 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1749 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1750 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1751 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1752 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1753 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1754 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1755 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1756 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1757 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1758 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1759 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1760 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1761 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1762 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1763 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1764 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1765 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1766 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1767 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1768 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1769 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1770 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1771 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1772 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1773 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1774 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1775 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1776 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1777 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1778 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1779 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1780 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1781 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1782 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1783 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1784 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1785 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1786 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1787 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1788 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1789 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1790 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1791 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1792 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1793 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1794 DATA iden/36/
1795C ----------
1796C BEGIN CODE
1797C ----------
1798 DO ihel=1,nexternal
1799 jc(ihel) = +1
1800 ENDDO
1801 ans = 0d0
1802 DO ihel=1,ncomb
1803 t=matrix_cd_us_noh(p,h1,h2 ,nhel(1,ihel),jc(1))
1804 ans=ans+t
1805 ENDDO
1806 ans=ans/dble(iden)
1807 END
1808
1809
1810 REAL*8 FUNCTION matrix_cd_us_noh(P,H1,H2,NHEL,IC)
1811C
1812C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1813C By the MadGraph Development Team
1814C Please visit us at https://launchpad.net/madgraph5
1815C
1816C Returns amplitude squared summed/avg over colors
1817C for the point with external lines W(0:6,NEXTERNAL)
1818C
1819C Process: c d > u s ta+ ta- / h QED=4
1820C
1821 IMPLICIT NONE
1822C
1823C CONSTANTS
1824C
1825 INTEGER NGRAPHS
1826 parameter(ngraphs=11)
1827 INTEGER NEXTERNAL
1828 parameter(nexternal=6)
1829 INTEGER NWAVEFUNCS, NCOLOR
1830 parameter(nwavefuncs=10, ncolor=1)
1831 real*8 zero
1832 parameter(zero=0d0)
1833 COMPLEX*16 IMAG1
1834 parameter(imag1=(0d0,1d0))
1835C
1836C ARGUMENTS
1837C
1838 real*8 p(0:3,nexternal)
1839 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1840C
1841C LOCAL VARIABLES
1842C
1843 INTEGER I,J
1844 COMPLEX*16 ZTEMP
1845 real*8 denom(ncolor), cf(ncolor,ncolor)
1846 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1847 COMPLEX*16 W(18,NWAVEFUNCS)
1848 COMPLEX*16 DUM0,DUM1
1849 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1850C
1851C GLOBAL VARIABLES
1852C
1853 include 'coupl.inc'
1854C
1855C COLOR DATA
1856C
1857 DATA denom(1)/1/
1858 DATA (cf(i, 1),i= 1, 1) / 9/
1859C 1 T(3,2) T(4,1)
1860
1861 INTEGER H1,H2
1862 real*8 matrix
1863 matrix_cd_us_noh=0d0
1864 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1865 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1866
1867
1868C ----------
1869C BEGIN CODE
1870C ----------
1871 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1872 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1873 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1874 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1875 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1876 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1877 CALL ffv2_3(w(1,1),w(1,4),gc_100,mw,ww,w(1,7))
1878 CALL ffv2_3(w(1,2),w(1,3),gc_100,mw,ww,w(1,8))
1879 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1880C Amplitude(s) for diagram number 1
1881 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1882 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1883C Amplitude(s) for diagram number 2
1884 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1885 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1886C Amplitude(s) for diagram number 3
1887 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1888 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1889C Amplitude(s) for diagram number 4
1890 CALL ffv1_0(w(1,6),w(1,3),w(1,9),gc_2,amp(4))
1891 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,5))
1892C Amplitude(s) for diagram number 5
1893 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1894C Amplitude(s) for diagram number 6
1895 CALL ffv2_5_0(w(1,6),w(1,3),w(1,10),gc_51,gc_58,amp(6))
1896C Amplitude(s) for diagram number 7
1897 CALL ffv2_3_0(w(1,2),w(1,5),w(1,10),gc_50,gc_58,amp(7))
1898 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1899C Amplitude(s) for diagram number 8
1900 CALL ffv1_0(w(1,5),w(1,4),w(1,9),gc_1,amp(8))
1901 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
1902C Amplitude(s) for diagram number 9
1903 CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_100,amp(9))
1904C Amplitude(s) for diagram number 10
1905 CALL ffv2_3_0(w(1,5),w(1,4),w(1,10),gc_50,gc_58,amp(10))
1906 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1907C Amplitude(s) for diagram number 11
1908 CALL ffv2_0(w(1,5),w(1,4),w(1,8),gc_100,amp(11))
1909 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1910 $ +amp(9)+amp(10)+amp(11)
1911
1912 matrix = 0.d0
1913 DO i = 1, ncolor
1914 ztemp = (0.d0,0.d0)
1915 DO j = 1, ncolor
1916 ztemp = ztemp + cf(j,i)*jamp(j)
1917 ENDDO
1918 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1919 ENDDO
1920 matrix_cd_us_noh=matrix
1921 ENDIF ! CLOSES H1 IF
1922 ENDIF ! CLOSES H2 IF
1923 END
1924
1925
1926C ------- end subprocesses CD->jjtautau, no Higgs
1927
1928