C++ Interface to Tauola
DC.f
1
2c --------------------- subprocesses initiated by DC: 27.10.2015
3
4 SUBROUTINE dc(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
11
12 ans=0.d0
13
14 IF(key.EQ.1) THEN
15 IF(i3.EQ.2 .AND. i4.EQ.1) CALL dc_ud_h(p,h1,h2,ans)
16 IF(i3.EQ.2 .AND. i4.EQ.3) CALL dc_us_h(p,h1,h2,ans)
17 IF(i3.EQ.4 .AND. i4.EQ.1) CAll dc_cd_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.3) CALL dc_cs_h(p,h1,h2,ans)
19 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
20 CALL dc_ud_h(p,h1,h2,ans1)
21 CALL dc_us_h(p,h1,h2,ans2)
22 CALL dc_cd_h(p,h1,h2,ans3)
23 CALL dc_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 dc_ud_noh(p,h1,h2,ans)
28 IF(i3.EQ.2 .AND. i4.EQ.3) CALL dc_us_noh(p,h1,h2,ans)
29 IF(i3.EQ.4 .AND. i4.EQ.1) CAll dc_cd_noh(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.3) CALL dc_cs_noh(p,h1,h2,ans)
31 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
32 CALL dc_ud_noh(p,h1,h2,ans1)
33 CALL dc_us_noh(p,h1,h2,ans2)
34 CALL dc_cd_noh(p,h1,h2,ans3)
35 CALL dc_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 DC
43
44
45
46C ----- begin subprocesses DC->jjH with Higgs->tautau
47
48 SUBROUTINE dc_cd_h(P,H1,H2,ANS)
49C
50C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
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: d c > 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_dc_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_dc_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_dc_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: d c > 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_dc_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
229C ----------
230C BEGIN CODE
231C ----------
232 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
233 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
234 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
235 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
236 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
237 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
238 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
239 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,6))
240 CALL ffv2_3(w(1,2),w(1,4),gc_44,mw,ww,w(1,5))
241C Amplitude(s) for diagram number 1
242 CALL vvs1_0(w(1,5),w(1,6),w(1,7),gc_72,amp(1))
243 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
244 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,4))
245C Amplitude(s) for diagram number 2
246 CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_81,amp(2))
247 jamp(1)=-amp(1)
248 jamp(2)=+amp(2)
249
250 matrix = 0.d0
251 DO i = 1, ncolor
252 ztemp = (0.d0,0.d0)
253 DO j = 1, ncolor
254 ztemp = ztemp + cf(j,i)*jamp(j)
255 ENDDO
256 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
257 ENDDO
258 matrix_dc_cd_h=matrix
259 ENDIF ! CLOSES H1 IF
260 ENDIF ! CLOSES H2 IF
261 END
262
263c==========================================
264
265 SUBROUTINE dc_cs_h(P,H1,H2,ANS)
266C
267C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
268C By the MadGraph Development Team
269C Please visit us at https://launchpad.net/madgraph5
270C
271C MadGraph StandAlone Version
272C
273C Returns amplitude squared summed/avg over colors
274C and helicities
275C for the point in phase space P(0:3,NEXTERNAL)
276C
277C Process: d c > c s h WEIGHTED=6
278C * Decay: h > ta+ ta- WEIGHTED=2
279C
280 IMPLICIT NONE
281C
282C CONSTANTS
283C
284 INTEGER NEXTERNAL
285 parameter(nexternal=6)
286 INTEGER NCOMB
287 parameter( ncomb=64)
288C
289C ARGUMENTS
290C
291 real*8 p(0:3,nexternal),ans
292 INTEGER H1,H2
293C
294C LOCAL VARIABLES
295C
296 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
297 real*8 t
298 real*8 matrix_dc_cs_h
299 INTEGER IHEL,IDEN, I
300 INTEGER JC(NEXTERNAL)
301 LOGICAL GOODHEL(NCOMB)
302 DATA ntry/0/
303 DATA goodhel/ncomb*.false./
304 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
305 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
306 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
307 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
308 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
309 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
310 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
311 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
312 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
313 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
314 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
315 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
316 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
317 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
318 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
319 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
320 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
321 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
322 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
323 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
324 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
325 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
326 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
327 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
328 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
329 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
330 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
331 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
332 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
333 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
334 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
335 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
336 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
337 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
338 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
339 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
340 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
341 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
342 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
343 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
344 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
345 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
346 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
347 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
348 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
349 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
350 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
351 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
352 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
353 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
354 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
355 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
356 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
357 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
358 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
359 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
360 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
361 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
362 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
363 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
364 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
365 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
366 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
367 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
368 DATA iden/36/
369C ----------
370C BEGIN CODE
371C ----------
372 DO ihel=1,nexternal
373 jc(ihel) = +1
374 ENDDO
375 ans = 0d0
376 DO ihel=1,ncomb
377 t=matrix_dc_cs_h(p ,h1,h2,nhel(1,ihel),jc(1))
378 ans=ans+t
379 ENDDO
380 ans=ans/dble(iden)
381 END
382
383
384 REAL*8 FUNCTION matrix_dc_cs_h(P,H1,H2,NHEL,IC)
385C
386C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
387C By the MadGraph Development Team
388C Please visit us at https://launchpad.net/madgraph5
389C
390C Returns amplitude squared summed/avg over colors
391C for the point with external lines W(0:6,NEXTERNAL)
392C
393C Process: d c > c s h WEIGHTED=6
394C * Decay: h > ta+ ta- WEIGHTED=2
395C
396 IMPLICIT NONE
397C
398C CONSTANTS
399C
400 INTEGER NGRAPHS
401 parameter(ngraphs=1)
402 INTEGER NEXTERNAL
403 parameter(nexternal=6)
404 INTEGER NWAVEFUNCS, NCOLOR
405 parameter(nwavefuncs=7, ncolor=1)
406 real*8 zero
407 parameter(zero=0d0)
408 COMPLEX*16 IMAG1
409 parameter(imag1=(0d0,1d0))
410C
411C ARGUMENTS
412C
413 real*8 p(0:3,nexternal)
414 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
415C
416C LOCAL VARIABLES
417C
418 INTEGER I,J
419 COMPLEX*16 ZTEMP
420 real*8 denom(ncolor), cf(ncolor,ncolor)
421 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
422 COMPLEX*16 W(18,NWAVEFUNCS)
423 COMPLEX*16 DUM0,DUM1
424 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
425C
426C GLOBAL VARIABLES
427C
428 include 'coupl.inc'
429C
430C COLOR DATA
431C
432 DATA denom(1)/1/
433 DATA (cf(i, 1),i= 1, 1) / 9/
434C 1 T(3,1) T(4,2)
435
436 INTEGER H1,H2
437 real*8 matrix
438 matrix_dc_cs_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 ixxxxx(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 oxxxxx(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,3),gc_44,mw,ww,w(1,6))
453 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,3))
454C Amplitude(s) for diagram number 1
455 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
456 jamp(1)=-amp(1)
457
458 matrix = 0.d0
459 DO i = 1, ncolor
460 ztemp = (0.d0,0.d0)
461 DO j = 1, ncolor
462 ztemp = ztemp + cf(j,i)*jamp(j)
463 ENDDO
464 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
465 ENDDO
466 matrix_dc_cs_h=matrix
467 ENDIF
468 ENDIF
469 END
470
471
472
473C===============================================
474
475 SUBROUTINE dc_ud_h(P,H1,H2,ANS)
476C
477C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
478C By the MadGraph Development Team
479C Please visit us at https://launchpad.net/madgraph5
480C
481C MadGraph StandAlone Version
482C
483C Returns amplitude squared summed/avg over colors
484C and helicities
485C for the point in phase space P(0:3,NEXTERNAL)
486C
487C Process: d c > u d h WEIGHTED=6
488C * Decay: h > ta+ ta- WEIGHTED=2
489C
490 IMPLICIT NONE
491C
492C CONSTANTS
493C
494 INTEGER NEXTERNAL
495 parameter(nexternal=6)
496 INTEGER NCOMB
497 parameter( ncomb=64)
498C
499C ARGUMENTS
500C
501 real*8 p(0:3,nexternal),ans
502 INTEGER H1,H2
503C
504C LOCAL VARIABLES
505C
506 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
507 real*8 t
508 real*8 matrix_dc_ud_h
509 INTEGER IHEL,IDEN, I
510 INTEGER JC(NEXTERNAL)
511 LOGICAL GOODHEL(NCOMB)
512 DATA ntry/0/
513 DATA goodhel/ncomb*.false./
514 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
515 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
516 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
517 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
518 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
519 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
520 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
521 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
522 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
523 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
524 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
525 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
526 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
527 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
528 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
529 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
530 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
531 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
532 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
533 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
534 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
535 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
536 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
537 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
538 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
539 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
540 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
541 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
542 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
543 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
544 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
545 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
546 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
547 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
548 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
549 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
550 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
551 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
552 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
553 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
554 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
555 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
556 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
557 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
558 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
559 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
560 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
561 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
562 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
563 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
564 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
565 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
566 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
567 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
568 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
569 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
570 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
571 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
572 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
573 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
574 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
575 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
576 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
577 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
578 DATA iden/36/
579C ----------
580C BEGIN CODE
581C ----------
582 DO ihel=1,nexternal
583 jc(ihel) = +1
584 ENDDO
585 ans = 0d0
586 DO ihel=1,ncomb
587 t=matrix_dc_ud_h(p ,h1,h2,nhel(1,ihel),jc(1))
588 ans=ans+t
589 ENDDO
590 ans=ans/dble(iden)
591 END
592
593
594 REAL*8 FUNCTION matrix_dc_ud_h(P,H1,H2,NHEL,IC)
595C
596C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
597C By the MadGraph Development Team
598C Please visit us at https://launchpad.net/madgraph5
599C
600C Returns amplitude squared summed/avg over colors
601C for the point with external lines W(0:6,NEXTERNAL)
602C
603C Process: d c > u d h WEIGHTED=6
604C * Decay: h > ta+ ta- WEIGHTED=2
605C
606 IMPLICIT NONE
607C
608C CONSTANTS
609C
610 INTEGER NGRAPHS
611 parameter(ngraphs=1)
612 INTEGER NEXTERNAL
613 parameter(nexternal=6)
614 INTEGER NWAVEFUNCS, NCOLOR
615 parameter(nwavefuncs=7, ncolor=1)
616 real*8 zero
617 parameter(zero=0d0)
618 COMPLEX*16 IMAG1
619 parameter(imag1=(0d0,1d0))
620C
621C ARGUMENTS
622C
623 real*8 p(0:3,nexternal)
624 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
625C
626C LOCAL VARIABLES
627C
628 INTEGER I,J
629 COMPLEX*16 ZTEMP
630 real*8 denom(ncolor), cf(ncolor,ncolor)
631 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
632 COMPLEX*16 W(18,NWAVEFUNCS)
633 COMPLEX*16 DUM0,DUM1
634 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
635C
636C GLOBAL VARIABLES
637C
638 include 'coupl.inc'
639C
640C COLOR DATA
641C
642 DATA denom(1)/1/
643 DATA (cf(i, 1),i= 1, 1) / 9/
644C 1 T(3,1) T(4,2)
645
646 INTEGER H1,H2
647 REAL *8 MATRIX
648 matrix_dc_ud_h=0d0
649 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
650 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
651
652C ----------
653C BEGIN CODE
654C ----------
655 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
656 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
657 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
658 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
659 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
660 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
661 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
662 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
663 CALL ffv2_3(w(1,2),w(1,4),gc_44,mw,ww,w(1,3))
664C Amplitude(s) for diagram number 1
665 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
666 jamp(1)=-amp(1)
667
668 matrix = 0.d0
669 DO i = 1, ncolor
670 ztemp = (0.d0,0.d0)
671 DO j = 1, ncolor
672 ztemp = ztemp + cf(j,i)*jamp(j)
673 ENDDO
674 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
675 ENDDO
676 matrix_dc_ud_h=matrix
677 ENDIF
678 ENDIF
679 END
680
681
682
683C========================================
684
685 SUBROUTINE dc_us_h(P,H1,H2,ANS)
686C
687C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
688C By the MadGraph Development Team
689C Please visit us at https://launchpad.net/madgraph5
690C
691C MadGraph StandAlone Version
692C
693C Returns amplitude squared summed/avg over colors
694C and helicities
695C for the point in phase space P(0:3,NEXTERNAL)
696C
697C Process: d c > u s h WEIGHTED=6
698C * Decay: h > ta+ ta- WEIGHTED=2
699C
700 IMPLICIT NONE
701C
702C CONSTANTS
703C
704 INTEGER NEXTERNAL
705 parameter(nexternal=6)
706 INTEGER NCOMB
707 parameter( ncomb=64)
708C
709C ARGUMENTS
710C
711 real*8 p(0:3,nexternal),ans
712 INTEGER H1,H2
713C
714C LOCAL VARIABLES
715C
716 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
717 real*8 t
718 real*8 matrix_dc_us_h
719 INTEGER IHEL,IDEN, I
720 INTEGER JC(NEXTERNAL)
721 LOGICAL GOODHEL(NCOMB)
722 DATA ntry/0/
723 DATA goodhel/ncomb*.false./
724 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
725 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
726 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
727 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
728 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
729 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
730 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
731 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
732 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
733 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
734 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
735 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
736 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
737 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
738 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
739 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
740 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
741 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
742 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
743 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
744 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
745 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
746 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
747 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
748 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
749 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
750 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
751 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
752 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
753 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
754 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
755 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
756 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
757 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
758 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
759 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
760 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
761 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
762 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
763 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
764 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
765 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
766 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
767 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
768 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
769 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
770 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
771 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
772 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
773 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
774 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
775 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
776 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
777 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
778 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
779 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
780 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
781 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
782 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
783 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
784 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
785 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
786 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
787 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
788 DATA iden/36/
789C ----------
790C BEGIN CODE
791C ----------
792 DO ihel=1,nexternal
793 jc(ihel) = +1
794 ENDDO
795 ans = 0d0
796 DO ihel=1,ncomb
797 t=matrix_dc_us_h(p ,h1,h2,nhel(1,ihel),jc(1))
798 ans=ans+t
799 ENDDO
800 ans=ans/dble(iden)
801 END
802
803
804 REAL*8 FUNCTION matrix_dc_us_h(P,H1,H2,NHEL,IC)
805C
806C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
807C By the MadGraph Development Team
808C Please visit us at https://launchpad.net/madgraph5
809C
810C Returns amplitude squared summed/avg over colors
811C for the point with external lines W(0:6,NEXTERNAL)
812C
813C Process: d c > u s h WEIGHTED=6
814C * Decay: h > ta+ ta- WEIGHTED=2
815C
816 IMPLICIT NONE
817C
818C CONSTANTS
819C
820 INTEGER NGRAPHS
821 parameter(ngraphs=1)
822 INTEGER NEXTERNAL
823 parameter(nexternal=6)
824 INTEGER NWAVEFUNCS, NCOLOR
825 parameter(nwavefuncs=7, ncolor=1)
826 real*8 zero
827 parameter(zero=0d0)
828 COMPLEX*16 IMAG1
829 parameter(imag1=(0d0,1d0))
830C
831C ARGUMENTS
832C
833 real*8 p(0:3,nexternal)
834 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
835C
836C LOCAL VARIABLES
837C
838 INTEGER I,J
839 COMPLEX*16 ZTEMP
840 real*8 denom(ncolor), cf(ncolor,ncolor)
841 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
842 COMPLEX*16 W(18,NWAVEFUNCS)
843 COMPLEX*16 DUM0,DUM1
844 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
845C
846C GLOBAL VARIABLES
847C
848 include 'coupl.inc'
849C
850C COLOR DATA
851C
852 DATA denom(1)/1/
853 DATA (cf(i, 1),i= 1, 1) / 9/
854C 1 T(3,1) T(4,2)
855
856 INTEGER H1,H2
857 real*8 matrix
858 matrix_dc_us_h=0d0
859 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
860 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
861
862C ----------
863C BEGIN CODE
864C ----------
865 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
866 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
867 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
868 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
869 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
870 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
871 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
872 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
873 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,3))
874C Amplitude(s) for diagram number 1
875 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
876 jamp(1)=-amp(1)
877
878 matrix = 0.d0
879 DO i = 1, ncolor
880 ztemp = (0.d0,0.d0)
881 DO j = 1, ncolor
882 ztemp = ztemp + cf(j,i)*jamp(j)
883 ENDDO
884 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
885 ENDDO
886 matrix_dc_us_h=matrix
887 ENDIF
888 ENDIF
889 END
890
891
892
893
894C --- end of subprocess DC->jjH, H->tautau
895
896C ----- begin subprocesses DC->jjtautau, no Higgs
897
898 SUBROUTINE dc_cd_noh(P,H1,H2,ANS)
899C
900C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
901C By the MadGraph Development Team
902C Please visit us at https://launchpad.net/madgraph5
903C
904C MadGraph StandAlone Version
905C
906C Returns amplitude squared summed/avg over colors
907C and helicities
908C for the point in phase space P(0:3,NEXTERNAL)
909C
910C Process: d c > c d ta+ ta- / h QED=4
911C
912 IMPLICIT NONE
913C
914C CONSTANTS
915C
916 INTEGER NEXTERNAL
917 parameter(nexternal=6)
918 INTEGER NCOMB
919 parameter( ncomb=64)
920C
921C ARGUMENTS
922C
923 real*8 p(0:3,nexternal),ans
924 INTEGER H1,H2
925C
926C LOCAL VARIABLES
927C
928 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
929 real*8 t
930 real*8 matrix_dc_cd_noh
931 INTEGER IHEL,IDEN, I
932 INTEGER JC(NEXTERNAL)
933 LOGICAL GOODHEL(NCOMB)
934 DATA ntry/0/
935 DATA goodhel/ncomb*.false./
936 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
937 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
938 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
939 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
940 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
941 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
942 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
943 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
944 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
945 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
946 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
947 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
948 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
949 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
950 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
951 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
952 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
953 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
954 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
955 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
956 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
957 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
958 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
959 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
960 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
961 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
962 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
963 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
964 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
965 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
966 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
967 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
968 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
969 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
970 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
971 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
972 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
973 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
974 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
975 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
976 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
977 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
978 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
979 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
980 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
981 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
982 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
983 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
984 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
985 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
986 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
987 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
988 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
989 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
990 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
991 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
992 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
993 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
994 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
995 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
996 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
997 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
998 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
999 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1000 DATA iden/36/
1001C ----------
1002C BEGIN CODE
1003C ----------
1004 DO ihel=1,nexternal
1005 jc(ihel) = +1
1006 ENDDO
1007 ans = 0d0
1008 DO ihel=1,ncomb
1009 t=matrix_dc_cd_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1010 ans=ans+t
1011 ENDDO
1012 ans=ans/dble(iden)
1013 END
1014
1015
1016 REAL*8 FUNCTION matrix_dc_cd_noh(P,H1,H2,NHEL,IC)
1017C
1018C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1019C By the MadGraph Development Team
1020C Please visit us at https://launchpad.net/madgraph5
1021C
1022C Returns amplitude squared summed/avg over colors
1023C for the point with external lines W(0:6,NEXTERNAL)
1024C
1025C Process: d c > c d ta+ ta- / h QED=4
1026C
1027 IMPLICIT NONE
1028C
1029C CONSTANTS
1030C
1031 INTEGER NGRAPHS
1032 parameter(ngraphs=43)
1033 INTEGER NEXTERNAL
1034 parameter(nexternal=6)
1035 INTEGER NWAVEFUNCS, NCOLOR
1036 parameter(nwavefuncs=15, ncolor=2)
1037 real*8 zero
1038 parameter(zero=0d0)
1039 COMPLEX*16 IMAG1
1040 parameter(imag1=(0d0,1d0))
1041C
1042C ARGUMENTS
1043C
1044 real*8 p(0:3,nexternal)
1045 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1046C
1047C LOCAL VARIABLES
1048C
1049 INTEGER I,J
1050 COMPLEX*16 ZTEMP
1051 real*8 denom(ncolor), cf(ncolor,ncolor)
1052 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1053 COMPLEX*16 W(18,NWAVEFUNCS)
1054 COMPLEX*16 DUM0,DUM1
1055 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1056C
1057C GLOBAL VARIABLES
1058C
1059 include 'coupl.inc'
1060C
1061C COLOR DATA
1062C
1063 DATA denom(1)/1/
1064 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1065C 1 T(3,1) T(4,2)
1066 DATA denom(2)/1/
1067 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1068C 1 T(3,2) T(4,1)
1069
1070 INTEGER H1,H2
1071 real*8 matrix
1072 matrix_dc_cd_noh=0d0
1073 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1074 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1075
1076
1077C ----------
1078C BEGIN CODE
1079C ----------
1080 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1081 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1082 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1083 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1084 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1085 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1086 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,7))
1087 CALL ffv2_3(w(1,2),w(1,4),gc_44,mw,ww,w(1,8))
1088 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1089C Amplitude(s) for diagram number 1
1090 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1091 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1092C Amplitude(s) for diagram number 2
1093 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1094 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1095C Amplitude(s) for diagram number 3
1096 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1097 CALL ffv2_2(w(1,2),w(1,7),gc_44,zero,zero,w(1,11))
1098C Amplitude(s) for diagram number 4
1099 CALL ffv1_0(w(1,11),w(1,4),w(1,9),gc_1,amp(4))
1100 CALL ffv2_1(w(1,4),w(1,7),gc_44,zero,zero,w(1,12))
1101C Amplitude(s) for diagram number 5
1102 CALL ffv1_0(w(1,2),w(1,12),w(1,9),gc_2,amp(5))
1103C Amplitude(s) for diagram number 6
1104 CALL ffv2_3_0(w(1,11),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1105C Amplitude(s) for diagram number 7
1106 CALL ffv2_5_0(w(1,2),w(1,12),w(1,10),gc_51,gc_58,amp(7))
1107 CALL ffv1p0_3(w(1,1),w(1,4),gc_1,zero,zero,w(1,12))
1108 CALL ffv1p0_3(w(1,2),w(1,3),gc_2,zero,zero,w(1,11))
1109 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1110C Amplitude(s) for diagram number 8
1111 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1112 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1113C Amplitude(s) for diagram number 9
1114 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1115 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,14))
1116C Amplitude(s) for diagram number 10
1117 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1118C Amplitude(s) for diagram number 11
1119 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1120 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,13))
1121 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1122C Amplitude(s) for diagram number 12
1123 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1124 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1125C Amplitude(s) for diagram number 13
1126 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1127C Amplitude(s) for diagram number 14
1128 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1129C Amplitude(s) for diagram number 15
1130 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1131 CALL ffv1_2(w(1,2),w(1,12),gc_2,zero,zero,w(1,15))
1132C Amplitude(s) for diagram number 16
1133 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(16))
1134 CALL ffv1_1(w(1,3),w(1,12),gc_2,zero,zero,w(1,5))
1135C Amplitude(s) for diagram number 17
1136 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(17))
1137C Amplitude(s) for diagram number 18
1138 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(18))
1139C Amplitude(s) for diagram number 19
1140 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(19))
1141 CALL ffv1p0_3(w(1,1),w(1,4),gc_11,zero,zero,w(1,5))
1142 CALL ffv1_2(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1143C Amplitude(s) for diagram number 20
1144 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(20))
1145 CALL ffv1_1(w(1,3),w(1,5),gc_11,zero,zero,w(1,12))
1146C Amplitude(s) for diagram number 21
1147 CALL ffv1_0(w(1,2),w(1,12),w(1,9),gc_2,amp(21))
1148C Amplitude(s) for diagram number 22
1149 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(22))
1150C Amplitude(s) for diagram number 23
1151 CALL ffv2_5_0(w(1,2),w(1,12),w(1,10),gc_51,gc_58,amp(23))
1152 CALL ffv2_5_2(w(1,2),w(1,13),gc_51,gc_58,zero,zero,w(1,12))
1153C Amplitude(s) for diagram number 24
1154 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(24))
1155 CALL ffv2_5_1(w(1,3),w(1,13),gc_51,gc_58,zero,zero,w(1,15))
1156C Amplitude(s) for diagram number 25
1157 CALL ffv1_0(w(1,2),w(1,15),w(1,9),gc_2,amp(25))
1158C Amplitude(s) for diagram number 26
1159 CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(26))
1160C Amplitude(s) for diagram number 27
1161 CALL ffv2_5_0(w(1,2),w(1,15),w(1,10),gc_51,gc_58,amp(27))
1162 CALL ffv1_2(w(1,1),w(1,11),gc_1,zero,zero,w(1,15))
1163C Amplitude(s) for diagram number 28
1164 CALL ffv1_0(w(1,15),w(1,4),w(1,9),gc_1,amp(28))
1165 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,12))
1166C Amplitude(s) for diagram number 29
1167 CALL ffv1_0(w(1,12),w(1,4),w(1,11),gc_1,amp(29))
1168C Amplitude(s) for diagram number 30
1169 CALL ffv2_3_0(w(1,15),w(1,4),w(1,10),gc_50,gc_58,amp(30))
1170 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,15))
1171C Amplitude(s) for diagram number 31
1172 CALL ffv1_0(w(1,15),w(1,4),w(1,11),gc_1,amp(31))
1173 CALL ffv1p0_3(w(1,2),w(1,3),gc_11,zero,zero,w(1,11))
1174 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,2))
1175C Amplitude(s) for diagram number 32
1176 CALL ffv1_0(w(1,2),w(1,4),w(1,9),gc_1,amp(32))
1177C Amplitude(s) for diagram number 33
1178 CALL ffv1_0(w(1,12),w(1,4),w(1,11),gc_11,amp(33))
1179C Amplitude(s) for diagram number 34
1180 CALL ffv2_3_0(w(1,2),w(1,4),w(1,10),gc_50,gc_58,amp(34))
1181C Amplitude(s) for diagram number 35
1182 CALL ffv1_0(w(1,15),w(1,4),w(1,11),gc_11,amp(35))
1183 CALL ffv2_3_2(w(1,1),w(1,14),gc_50,gc_58,zero,zero,w(1,11))
1184C Amplitude(s) for diagram number 36
1185 CALL ffv1_0(w(1,11),w(1,4),w(1,9),gc_1,amp(36))
1186C Amplitude(s) for diagram number 37
1187 CALL ffv2_3_0(w(1,12),w(1,4),w(1,14),gc_50,gc_58,amp(37))
1188C Amplitude(s) for diagram number 38
1189 CALL ffv2_3_0(w(1,11),w(1,4),w(1,10),gc_50,gc_58,amp(38))
1190C Amplitude(s) for diagram number 39
1191 CALL ffv2_3_0(w(1,15),w(1,4),w(1,14),gc_50,gc_58,amp(39))
1192 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,14))
1193C Amplitude(s) for diagram number 40
1194 CALL ffv1_0(w(1,14),w(1,3),w(1,9),gc_2,amp(40))
1195C Amplitude(s) for diagram number 41
1196 CALL ffv2_0(w(1,12),w(1,3),w(1,8),gc_44,amp(41))
1197C Amplitude(s) for diagram number 42
1198 CALL ffv2_5_0(w(1,14),w(1,3),w(1,10),gc_51,gc_58,amp(42))
1199C Amplitude(s) for diagram number 43
1200 CALL ffv2_0(w(1,15),w(1,3),w(1,8),gc_44,amp(43))
1201 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1202 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1203 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1204 $ -amp(40)-amp(41)-amp(42)-amp(43)
1205 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1206 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1207 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1208 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1209 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1210 $ +amp(38)+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_dc_cd_noh=matrix
1221 ENDIF ! CLOSES H1 IF
1222 ENDIF ! CLOSES H2 IF
1223 END
1224
1225
1226 SUBROUTINE dc_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: d c > 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_dc_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_dc_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_dc_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: d c > 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,1) T(4,2)
1394
1395 INTEGER H1,H2
1396 real*8 matrix
1397 matrix_dc_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
1401C ----------
1402C BEGIN CODE
1403C ----------
1404 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1405 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1406 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1407 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1408 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1409 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1410 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,7))
1411 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,8))
1412 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1413C Amplitude(s) for diagram number 1
1414 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1415 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1416C Amplitude(s) for diagram number 2
1417 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1418 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1419C Amplitude(s) for diagram number 3
1420 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1421 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1422C Amplitude(s) for diagram number 4
1423 CALL ffv1_0(w(1,6),w(1,4),w(1,9),gc_1,amp(4))
1424 CALL ffv2_1(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1425C Amplitude(s) for diagram number 5
1426 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(5))
1427C Amplitude(s) for diagram number 6
1428 CALL ffv2_3_0(w(1,6),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1429C Amplitude(s) for diagram number 7
1430 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(7))
1431 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
1432C Amplitude(s) for diagram number 8
1433 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1434 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1435C Amplitude(s) for diagram number 9
1436 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_44,amp(9))
1437C Amplitude(s) for diagram number 10
1438 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1439 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1440C Amplitude(s) for diagram number 11
1441 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_44,amp(11))
1442 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1443 $ -amp(9)-amp(10)-amp(11)
1444
1445 matrix = 0.d0
1446 DO i = 1, ncolor
1447 ztemp = (0.d0,0.d0)
1448 DO j = 1, ncolor
1449 ztemp = ztemp + cf(j,i)*jamp(j)
1450 ENDDO
1451 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1452 ENDDO
1453 matrix_dc_cs_noh=matrix
1454 ENDIF ! CLOSES H1 IF
1455 ENDIF ! CLOSES H2 IF
1456 END
1457
1458 SUBROUTINE dc_ud_noh(P,H1,H2,ANS)
1459C
1460C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1461C By the MadGraph Development Team
1462C Please visit us at https://launchpad.net/madgraph5
1463C
1464C MadGraph StandAlone Version
1465C
1466C Returns amplitude squared summed/avg over colors
1467C and helicities
1468C for the point in phase space P(0:3,NEXTERNAL)
1469C
1470C Process: d c > u d ta+ ta- / h QED=4
1471C
1472 IMPLICIT NONE
1473C
1474C CONSTANTS
1475C
1476 INTEGER NEXTERNAL
1477 parameter(nexternal=6)
1478 INTEGER NCOMB
1479 parameter( ncomb=64)
1480C
1481C ARGUMENTS
1482C
1483 real*8 p(0:3,nexternal),ans
1484 INTEGER H1,H2
1485C
1486C LOCAL VARIABLES
1487C
1488 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1489 real*8 t
1490 real*8 matrix_dc_ud_noh
1491 INTEGER IHEL,IDEN, I
1492 INTEGER JC(NEXTERNAL)
1493 LOGICAL GOODHEL(NCOMB)
1494 DATA ntry/0/
1495 DATA goodhel/ncomb*.false./
1496 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1497 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1498 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1499 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1500 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1501 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1502 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1503 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1504 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1505 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1506 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1507 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1508 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1509 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1510 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1511 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1512 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1513 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1514 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1515 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1516 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1517 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1518 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1519 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1520 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1521 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1522 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1523 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1524 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1525 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1526 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1527 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1528 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1529 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1530 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1531 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1532 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1533 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1534 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1535 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1536 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1537 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1538 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1539 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1540 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1541 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1542 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1543 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1544 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1545 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1546 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1547 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1548 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1549 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1550 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1551 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1552 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1553 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1554 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1555 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1556 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1557 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1558 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1559 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1560 DATA iden/36/
1561C ----------
1562C BEGIN CODE
1563C ----------
1564 DO ihel=1,nexternal
1565 jc(ihel) = +1
1566 ENDDO
1567 ans = 0d0
1568 DO ihel=1,ncomb
1569 t=matrix_dc_ud_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1570 ans=ans+t
1571 ENDDO
1572 ans=ans/dble(iden)
1573 END
1574
1575
1576 REAL*8 FUNCTION matrix_dc_ud_noh(P,H1,H2,NHEL,IC)
1577C
1578C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1579C By the MadGraph Development Team
1580C Please visit us at https://launchpad.net/madgraph5
1581C
1582C Returns amplitude squared summed/avg over colors
1583C for the point with external lines W(0:6,NEXTERNAL)
1584C
1585C Process: d c > u d ta+ ta- / h QED=4
1586C
1587 IMPLICIT NONE
1588C
1589C CONSTANTS
1590C
1591 INTEGER NGRAPHS
1592 parameter(ngraphs=11)
1593 INTEGER NEXTERNAL
1594 parameter(nexternal=6)
1595 INTEGER NWAVEFUNCS, NCOLOR
1596 parameter(nwavefuncs=10, ncolor=1)
1597 real*8 zero
1598 parameter(zero=0d0)
1599 COMPLEX*16 IMAG1
1600 parameter(imag1=(0d0,1d0))
1601C
1602C ARGUMENTS
1603C
1604 real*8 p(0:3,nexternal)
1605 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1606C
1607C LOCAL VARIABLES
1608C
1609 INTEGER I,J
1610 COMPLEX*16 ZTEMP
1611 real*8 denom(ncolor), cf(ncolor,ncolor)
1612 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1613 COMPLEX*16 W(18,NWAVEFUNCS)
1614 COMPLEX*16 DUM0,DUM1
1615 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1616C
1617C GLOBAL VARIABLES
1618C
1619 include 'coupl.inc'
1620C
1621C COLOR DATA
1622C
1623 DATA denom(1)/1/
1624 DATA (cf(i, 1),i= 1, 1) / 9/
1625C 1 T(3,1) T(4,2)
1626
1627 INTEGER H1,H2
1628 real*8 matrix
1629 matrix_dc_ud_noh=0d0
1630 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1631 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1632
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,3),gc_100,mw,ww,w(1,7))
1644 CALL ffv2_3(w(1,2),w(1,4),gc_44,mw,ww,w(1,8))
1645 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1646C Amplitude(s) for diagram number 1
1647 CALL ffv2_0(w(1,5),w(1,9),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,8),w(1,7),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,8),w(1,7),w(1,10),gc_53,amp(3))
1654 CALL ffv2_2(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
1655C Amplitude(s) for diagram number 4
1656 CALL ffv1_0(w(1,6),w(1,4),w(1,9),gc_1,amp(4))
1657 CALL ffv2_1(w(1,4),w(1,7),gc_44,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_2,amp(5))
1660C Amplitude(s) for diagram number 6
1661 CALL ffv2_3_0(w(1,6),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1662C Amplitude(s) for diagram number 7
1663 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(7))
1664 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1665C Amplitude(s) for diagram number 8
1666 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1667 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1668C Amplitude(s) for diagram number 9
1669 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
1670C Amplitude(s) for diagram number 10
1671 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1672 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1673C Amplitude(s) for diagram number 11
1674 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,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_dc_ud_noh=matrix
1687 ENDIF ! CLOSES H1 IF
1688 ENDIF ! CLOSES H2 IF
1689 END
1690
1691
1692 SUBROUTINE dc_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: d c > 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_dc_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_dc_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_dc_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: d c > 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,1) T(4,2)
1860
1861 INTEGER H1,H2
1862 real*8 matrix
1863 matrix_dc_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
1867C ----------
1868C BEGIN CODE
1869C ----------
1870 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1871 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1872 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1873 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1874 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1875 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1876 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
1877 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,8))
1878 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1879C Amplitude(s) for diagram number 1
1880 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1881 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1882C Amplitude(s) for diagram number 2
1883 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1884 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1885C Amplitude(s) for diagram number 3
1886 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1887 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1888C Amplitude(s) for diagram number 4
1889 CALL ffv1_0(w(1,6),w(1,4),w(1,9),gc_1,amp(4))
1890 CALL ffv2_1(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1891C Amplitude(s) for diagram number 5
1892 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(5))
1893C Amplitude(s) for diagram number 6
1894 CALL ffv2_3_0(w(1,6),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1895C Amplitude(s) for diagram number 7
1896 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(7))
1897 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1898C Amplitude(s) for diagram number 8
1899 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1900 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1901C Amplitude(s) for diagram number 9
1902 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
1903C Amplitude(s) for diagram number 10
1904 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1905 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1906C Amplitude(s) for diagram number 11
1907 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
1908 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1909 $ -amp(9)-amp(10)-amp(11)
1910
1911 matrix = 0.d0
1912 DO i = 1, ncolor
1913 ztemp = (0.d0,0.d0)
1914 DO j = 1, ncolor
1915 ztemp = ztemp + cf(j,i)*jamp(j)
1916 ENDDO
1917 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1918 ENDDO
1919 matrix_dc_us_noh=matrix
1920 ENDIF ! CLOSES H1 IF
1921 ENDIF ! CLOSES H2 IF
1922 END
1923
1924
1925C ------- end subprocesses DC->jjtautau, no Higgs
1926C DONE 27.10.2015
1927
1928