C++ Interface to Tauola
DDX.f
1C --------- begin processes initiated by DDbar
2
3 SUBROUTINE ddx(P,I3,I4,H1,H2,KEY,ANS)
4 IMPLICIT NONE
5
6 INTEGER I3,I4, H1,H2, KEY
7 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4, ans5,ans6
8
9 real*8 pp(0:3,6)
10
11C I3=0,1,2,3 OR 4 I4=0,-1,-2,-3, OR -4
12C or I3=21 AND I4=21 -> GLUONS AS FINAL JETS
13
14 ans=0.d0
15
16 ! switch 4-vectors for I3, I4 and provide/(add contribution from) new configuration
17 pp(0:3,1) = p(0:3,1)
18 pp(0:3,2) = p(0:3,2)
19 pp(0:3,5) = p(0:3,5)
20 pp(0:3,6) = p(0:3,6)
21
22 pp(0:3,3) = p(0:3,4)
23 pp(0:3,4) = p(0:3,3)
24
25 IF(key.EQ.1) THEN
26 IF(i3.EQ.1 .AND. i4.EQ.-1) CALL ddx_ddx_h(p,h1,h2,ans)
27 IF(i3.EQ.2 .AND. i4.EQ.-2) CALL ddx_uux_h(p,h1,h2,ans)
28 IF(i3.EQ.3 .AND. i4.EQ.-3) CAll ddx_ssx_h(p,h1,h2,ans)
29 IF(i3.EQ.4 .AND. i4.EQ.-4) CALL ddx_ccx_h(p,h1,h2,ans)
30 IF(i3.EQ.2 .AND. i4.EQ.-4) CALL ddx_ucx_h(p,h1,h2,ans)
31 IF(i3.EQ.4 .AND. i4.EQ.-2) CALL ddx_cux_h(p,h1,h2,ans)
32 IF(i3.EQ.-1 .AND. i4.EQ.1) CALL ddx_ddx_h(pp,h1,h2,ans)
33 IF(i3.EQ.-2 .AND. i4.EQ.2) CALL ddx_uux_h(pp,h1,h2,ans)
34 IF(i3.EQ.-3 .AND. i4.EQ.3) CAll ddx_ssx_h(pp,h1,h2,ans)
35 IF(i3.EQ.-4 .AND. i4.EQ.4) CALL ddx_ccx_h(pp,h1,h2,ans)
36 IF(i3.EQ.-2 .AND. i4.EQ.4) CALL ddx_cux_h(pp,h1,h2,ans)
37 IF(i3.EQ.-4 .AND. i4.EQ.2) CALL ddx_ucx_h(pp,h1,h2,ans)
38 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
39 CALL ddx_ddx_h(p,h1,h2,ans1)
40 CALL ddx_uux_h(p,h1,h2,ans2)
41 CALL ddx_ssx_h(p,h1,h2,ans3)
42 CALL ddx_ccx_h(p,h1,h2,ans4)
43 CALL ddx_ucx_h(p,h1,h2,ans5)
44 CALL ddx_cux_h(p,h1,h2,ans6)
45 ans=ans1+ans2+ans3+ans4+ans5+ans6
46 CALL ddx_ddx_h(pp,h1,h2,ans1)
47 CALL ddx_uux_h(pp,h1,h2,ans2)
48 CALL ddx_ssx_h(pp,h1,h2,ans3)
49 CALL ddx_ccx_h(pp,h1,h2,ans4)
50 CALL ddx_ucx_h(pp,h1,h2,ans5)
51 CALL ddx_cux_h(pp,h1,h2,ans6)
52 ans=ans1+ans2+ans3+ans4+ans5+ans6+ans
53 ENDIF
54 ELSE IF(key.EQ.0) THEN
55 IF(i3.EQ.1 .AND. i4.EQ.-1) CALL ddx_ddx_noh(p,h1,h2,ans)
56 IF(i3.EQ.2 .AND. i4.EQ.-2) CALL ddx_uux_noh(p,h1,h2,ans)
57 IF(i3.EQ.3 .AND. i4.EQ.-3) CAll ddx_ssx_noh(p,h1,h2,ans)
58 IF(i3.EQ.4 .AND. i4.EQ.-4) CALL ddx_ccx_noh(p,h1,h2,ans)
59 IF(i3.EQ.2 .AND. i4.EQ.-4) CALL ddx_ucx_noh(p,h1,h2,ans)
60 IF(i3.EQ.4 .AND. i4.EQ.-2) CALL ddx_cux_noh(p,h1,h2,ans)
61
62 IF(i3.EQ.21 .AND. i4.EQ.21) CALL ddx_gg_noh(p,h1,h2,ans)
63 IF(i3.EQ.-1 .AND. i4.EQ.1) CALL ddx_ddx_noh(pp,h1,h2,ans)
64 IF(i3.EQ.-2 .AND. i4.EQ.2) CALL ddx_uux_noh(pp,h1,h2,ans)
65 IF(i3.EQ.-3 .AND. i4.EQ.3) CAll ddx_ssx_noh(pp,h1,h2,ans)
66 IF(i3.EQ.-4 .AND. i4.EQ.4) CALL ddx_ccx_noh(pp,h1,h2,ans)
67 IF(i3.EQ.-2 .AND. i4.EQ.4) CALL ddx_cux_noh(pp,h1,h2,ans)
68 IF(i3.EQ.-4 .AND. i4.EQ.2) CALL ddx_ucx_noh(pp,h1,h2,ans)
69 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
70 CALL ddx_ddx_noh(p,h1,h2,ans1)
71 CALL ddx_uux_noh(p,h1,h2,ans2)
72 CALL ddx_ssx_noh(p,h1,h2,ans3)
73 CALL ddx_ccx_noh(p,h1,h2,ans4)
74 CALL ddx_ucx_noh(p,h1,h2,ans5)
75 CALL ddx_cux_noh(p,h1,h2,ans6)
76 ans=ans1+ans2+ans3+ans4+ans5+ans6
77 CALL ddx_ddx_noh(pp,h1,h2,ans1)
78 CALL ddx_uux_noh(pp,h1,h2,ans2)
79 CALL ddx_ssx_noh(pp,h1,h2,ans3)
80 CALL ddx_ccx_noh(pp,h1,h2,ans4)
81 CALL ddx_ucx_noh(pp,h1,h2,ans5)
82 CALL ddx_cux_noh(pp,h1,h2,ans6)
83 ans=ans1+ans2+ans3+ans4+ans5+ans6+ans
84 CALL ddx_gg_noh(p,h1,h2,ans1)
85 ans=ans1+ans
86 ENDIF
87 ELSE
88 WRITE(*,*) 'NOT FINISHED'
89 stop
90 ENDIF
91 END SUBROUTINE ddx
92
93
94
95C ----- begin subprocesses DDX->DDX with Higgs->tautau
96
97 SUBROUTINE ddx_ddx_h(P,H1,H2,ANS)
98C
99C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
100C By the MadGraph Development Team
101C Please visit us at https://launchpad.net/madgraph5
102C
103C MadGraph StandAlone Version
104C
105C Returns amplitude squared summed/avg over colors
106C and helicities
107C for the point in phase space P(0:3,NEXTERNAL)
108C
109C Process: d d~ > d d~ h WEIGHTED=6
110C * Decay: h > ta+ ta- WEIGHTED=2
111C
112 IMPLICIT NONE
113C
114C CONSTANTS
115C
116 INTEGER NEXTERNAL
117 parameter(nexternal=6)
118 INTEGER NCOMB
119 parameter( ncomb=64)
120C
121C ARGUMENTS
122C
123 real*8 p(0:3,nexternal),ans
124 INTEGER H1,H2
125C
126C LOCAL VARIABLES
127C
128 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
129 real*8 t
130 real*8 matrix_ddx_ddx_h
131 INTEGER IHEL,IDEN, I
132 INTEGER JC(NEXTERNAL)
133 LOGICAL GOODHEL(NCOMB)
134 DATA ntry/0/
135 DATA goodhel/ncomb*.false./
136 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
137 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
138 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
139 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
140 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
141 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
142 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
143 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
144 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
145 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
146 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
147 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
148 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
149 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
150 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
151 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
152 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
153 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
154 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
155 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
156 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
157 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
158 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
159 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
160 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
161 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
162 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
163 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
164 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
165 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
166 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
167 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
168 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
169 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
170 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
171 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
172 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
173 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
174 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
175 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
176 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
177 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
178 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
179 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
180 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
181 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
182 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
183 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
184 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
185 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
186 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
187 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
188 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
189 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
190 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
191 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
192 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
193 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
194 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
195 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
196 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
197 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
198 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
199 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
200 DATA iden/36/
201C ----------
202C BEGIN CODE
203C ----------
204 DO ihel=1,nexternal
205 jc(ihel) = +1
206 ENDDO
207 ans = 0d0
208 DO ihel=1,ncomb
209 t=matrix_ddx_ddx_h(p ,h1,h2,nhel(1,ihel),jc(1))
210 ans=ans+t
211 ENDDO
212 ans=ans/dble(iden)
213 END
214
215
216 REAL*8 FUNCTION matrix_ddx_ddx_h(P,H1,H2,NHEL,IC)
217C
218C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
219C By the MadGraph Development Team
220C Please visit us at https://launchpad.net/madgraph5
221C
222C Returns amplitude squared summed/avg over colors
223C for the point with external lines W(0:6,NEXTERNAL)
224C
225C Process: d d~ > d d~ h WEIGHTED=6
226C * Decay: h > ta+ ta- WEIGHTED=2
227C
228 IMPLICIT NONE
229C
230C CONSTANTS
231C
232 INTEGER NGRAPHS
233 parameter(ngraphs=2)
234 INTEGER NEXTERNAL
235 parameter(nexternal=6)
236 INTEGER NWAVEFUNCS, NCOLOR
237 parameter(nwavefuncs=7, ncolor=2)
238 real*8 zero
239 parameter(zero=0d0)
240 COMPLEX*16 IMAG1
241 parameter(imag1=(0d0,1d0))
242C
243C ARGUMENTS
244C
245 real*8 p(0:3,nexternal)
246 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
247C
248C LOCAL VARIABLES
249C
250 INTEGER I,J
251 COMPLEX*16 ZTEMP
252 real*8 denom(ncolor), cf(ncolor,ncolor)
253 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
254 COMPLEX*16 W(18,NWAVEFUNCS)
255 COMPLEX*16 DUM0,DUM1
256 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
257C
258C GLOBAL VARIABLES
259C
260 include 'coupl.inc'
261C
262C COLOR DATA
263C
264 DATA denom(1)/1/
265 DATA (cf(i, 1),i= 1, 2) / 9, 3/
266C 1 T(2,1) T(3,4)
267 DATA denom(2)/1/
268 DATA (cf(i, 2),i= 1, 2) / 3, 9/
269C 1 T(2,4) T(3,1)
270
271 INTEGER H1,H2
272 real*8 matrix
273 matrix_ddx_ddx_h=0.d0
274 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
275 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
276
277
278C ----------
279C BEGIN CODE
280C ----------
281 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
282 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
283 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
284 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
285 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
286 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
287 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
288 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
289 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
290C Amplitude(s) for diagram number 1
291 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
292 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
293 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
294C Amplitude(s) for diagram number 2
295 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
296 jamp(1)=-amp(1)
297 jamp(2)=+amp(2)
298
299 matrix = 0.d0
300 DO i = 1, ncolor
301 ztemp = (0.d0,0.d0)
302 DO j = 1, ncolor
303 ztemp = ztemp + cf(j,i)*jamp(j)
304 ENDDO
305 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
306 ENDDO
307 matrix_ddx_ddx_h=matrix
308 ENDIF
309 ENDIF
310 END
311
312C ----- begin subprocesses DDX->UUX with Higgs->tautau
313
314 SUBROUTINE ddx_uux_h(P,H1,H2,ANS)
315C
316C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
317C By the MadGraph Development Team
318C Please visit us at https://launchpad.net/madgraph5
319C
320C MadGraph StandAlone Version
321C
322C Returns amplitude squared summed/avg over colors
323C and helicities
324C for the point in phase space P(0:3,NEXTERNAL)
325C
326C Process: d d~ > u u~ h WEIGHTED=6
327C * Decay: h > ta+ ta- WEIGHTED=2
328C
329 IMPLICIT NONE
330C
331C CONSTANTS
332C
333 INTEGER NEXTERNAL
334 parameter(nexternal=6)
335 INTEGER NCOMB
336 parameter( ncomb=64)
337C
338C ARGUMENTS
339C
340 real*8 p(0:3,nexternal),ans
341 INTEGER H1,H2
342C
343C LOCAL VARIABLES
344C
345 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
346 real*8 t
347 real*8 matrix_ddx_uux_h
348 INTEGER IHEL,IDEN, I
349 INTEGER JC(NEXTERNAL)
350 LOGICAL GOODHEL(NCOMB)
351 DATA ntry/0/
352 DATA goodhel/ncomb*.false./
353 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
354 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
355 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
356 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
357 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
358 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
359 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
360 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
361 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
362 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
363 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
364 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
365 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
366 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
367 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
368 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
369 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
370 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
371 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
372 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
373 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
374 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
375 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
376 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
377 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
378 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
379 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
380 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
381 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
382 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
383 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
384 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
385 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
386 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
387 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
388 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
389 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
390 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
391 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
392 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
393 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
394 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
395 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
396 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
397 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
398 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
399 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
400 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
401 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
402 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
403 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
404 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
405 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
406 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
407 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
408 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
409 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
410 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
411 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
412 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
413 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
414 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
415 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
416 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
417 DATA iden/36/
418C ----------
419C BEGIN CODE
420C ----------
421 DO ihel=1,nexternal
422 jc(ihel) = +1
423 ENDDO
424 ans = 0d0
425 DO ihel=1,ncomb
426 t=matrix_ddx_uux_h(p ,h1,h2,nhel(1,ihel),jc(1))
427 ans=ans+t
428 ENDDO
429 ans=ans/dble(iden)
430 END
431
432
433 REAL*8 FUNCTION matrix_ddx_uux_h(P,H1,H2,NHEL,IC)
434C
435C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
436C By the MadGraph Development Team
437C Please visit us at https://launchpad.net/madgraph5
438C
439C Returns amplitude squared summed/avg over colors
440C for the point with external lines W(0:6,NEXTERNAL)
441C
442C Process: d d~ > u u~ h WEIGHTED=6
443C * Decay: h > ta+ ta- WEIGHTED=2
444C
445 IMPLICIT NONE
446C
447C CONSTANTS
448C
449 INTEGER NGRAPHS
450 parameter(ngraphs=2)
451 INTEGER NEXTERNAL
452 parameter(nexternal=6)
453 INTEGER NWAVEFUNCS, NCOLOR
454 parameter(nwavefuncs=7, ncolor=2)
455 real*8 zero
456 parameter(zero=0d0)
457 COMPLEX*16 IMAG1
458 parameter(imag1=(0d0,1d0))
459C
460C ARGUMENTS
461C
462 real*8 p(0:3,nexternal)
463 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
464C
465C LOCAL VARIABLES
466C
467 INTEGER I,J
468 COMPLEX*16 ZTEMP
469 real*8 denom(ncolor), cf(ncolor,ncolor)
470 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
471 COMPLEX*16 W(18,NWAVEFUNCS)
472 COMPLEX*16 DUM0,DUM1
473 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
474C
475C GLOBAL VARIABLES
476C
477 include 'coupl.inc'
478C
479C COLOR DATA
480C
481 DATA denom(1)/1/
482 DATA (cf(i, 1),i= 1, 2) / 9, 3/
483C 1 T(2,1) T(3,4)
484 DATA denom(2)/1/
485 DATA (cf(i, 2),i= 1, 2) / 3, 9/
486C 1 T(2,4) T(3,1)
487
488 INTEGER H1,H2
489 real*8 matrix
490 matrix_ddx_uux_h=0.d0
491 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
492 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
493
494C ----------
495C BEGIN CODE
496C ----------
497 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
498 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
499 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
500 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
501 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
502 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
503 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
504 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
505 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
506C Amplitude(s) for diagram number 1
507 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
508 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,5))
509 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
510C Amplitude(s) for diagram number 2
511 CALL vvs1_0(w(1,3),w(1,5),w(1,7),gc_72,amp(2))
512 jamp(1)=-amp(1)
513 jamp(2)=+amp(2)
514
515 matrix = 0.d0
516 DO i = 1, ncolor
517 ztemp = (0.d0,0.d0)
518 DO j = 1, ncolor
519 ztemp = ztemp + cf(j,i)*jamp(j)
520 ENDDO
521 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
522 ENDDO
523 matrix_ddx_uux_h=matrix
524 ENDIF
525 ENDIF
526 END
527
528C ----- begin subprocesses DDX->SSX with Higgs->tautau
529
530 SUBROUTINE ddx_ssx_h(P,H1,H2,ANS)
531C
532C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
533C By the MadGraph Development Team
534C Please visit us at https://launchpad.net/madgraph5
535C
536C MadGraph StandAlone Version
537C
538C Returns amplitude squared summed/avg over colors
539C and helicities
540C for the point in phase space P(0:3,NEXTERNAL)
541C
542C Process: d d~ > s s~ h WEIGHTED=6
543C * Decay: h > ta+ ta- WEIGHTED=2
544C
545 IMPLICIT NONE
546C
547C CONSTANTS
548C
549 INTEGER NEXTERNAL
550 parameter(nexternal=6)
551 INTEGER NCOMB
552 parameter( ncomb=64)
553C
554C ARGUMENTS
555C
556 real*8 p(0:3,nexternal),ans
557 INTEGER H1,H2
558C
559C LOCAL VARIABLES
560C
561 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
562 real*8 t
563 real*8 matrix_ddx_ssx_h
564 INTEGER IHEL,IDEN, I
565 INTEGER JC(NEXTERNAL)
566 LOGICAL GOODHEL(NCOMB)
567 DATA ntry/0/
568 DATA goodhel/ncomb*.false./
569 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
570 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
571 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
572 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
573 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
574 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
575 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
576 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
577 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
578 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
579 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
580 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
581 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
582 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
583 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
584 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
585 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
586 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
587 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
588 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
589 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
590 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
591 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
592 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
593 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
594 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
595 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
596 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
597 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
598 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
599 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
600 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
601 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
602 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
603 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
604 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
605 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
606 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
607 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
608 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
609 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
610 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
611 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
612 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
613 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
614 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
615 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
616 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
617 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
618 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
619 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
620 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
621 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
622 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
623 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
624 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
625 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
626 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
627 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
628 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
629 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
630 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
631 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
632 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
633 DATA iden/36/
634C ----------
635C BEGIN CODE
636C ----------
637 DO ihel=1,nexternal
638 jc(ihel) = +1
639 ENDDO
640 ans = 0d0
641 DO ihel=1,ncomb
642 t=matrix_ddx_ssx_h(p , h1, h2, nhel(1,ihel),jc(1))
643 ans=ans+t
644 ENDDO
645 ans=ans/dble(iden)
646 END
647
648
649 REAL*8 FUNCTION matrix_ddx_ssx_h(P,H1,H2,NHEL,IC)
650C
651C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
652C By the MadGraph Development Team
653C Please visit us at https://launchpad.net/madgraph5
654C
655C Returns amplitude squared summed/avg over colors
656C for the point with external lines W(0:6,NEXTERNAL)
657C
658C Process: d d~ > s s~ h WEIGHTED=6
659C * Decay: h > ta+ ta- WEIGHTED=2
660C
661 IMPLICIT NONE
662C
663C CONSTANTS
664C
665 INTEGER NGRAPHS
666 parameter(ngraphs=1)
667 INTEGER NEXTERNAL
668 parameter(nexternal=6)
669 INTEGER NWAVEFUNCS, NCOLOR
670 parameter(nwavefuncs=7, ncolor=1)
671 real*8 zero
672 parameter(zero=0d0)
673 COMPLEX*16 IMAG1
674 parameter(imag1=(0d0,1d0))
675C
676C ARGUMENTS
677C
678 real*8 p(0:3,nexternal)
679 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
680C
681C LOCAL VARIABLES
682C
683 INTEGER I,J
684 COMPLEX*16 ZTEMP
685 real*8 denom(ncolor), cf(ncolor,ncolor)
686 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
687 COMPLEX*16 W(18,NWAVEFUNCS)
688 COMPLEX*16 DUM0,DUM1
689 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
690C
691C GLOBAL VARIABLES
692C
693 include 'coupl.inc'
694C
695C COLOR DATA
696C
697 DATA denom(1)/1/
698 DATA (cf(i, 1),i= 1, 1) / 9/
699C 1 T(2,1) T(3,4)
700
701 INTEGER H1,H2
702 real*8 matrix
703 matrix_ddx_ssx_h=0.d0
704 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
705 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
706
707C ----------
708C BEGIN CODE
709C ----------
710 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
711 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
712 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
713 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
714 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
715 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
716 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
717 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
718 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,2))
719C Amplitude(s) for diagram number 1
720 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_81,amp(1))
721 jamp(1)=-amp(1)
722
723 matrix = 0.d0
724 DO i = 1, ncolor
725 ztemp = (0.d0,0.d0)
726 DO j = 1, ncolor
727 ztemp = ztemp + cf(j,i)*jamp(j)
728 ENDDO
729 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
730 ENDDO
731 matrix_ddx_ssx_h=matrix
732 ENDIF
733 ENDIF
734 END
735
736C ----- begin subprocesses DDX->CCX with Higgs->tautau
737
738 SUBROUTINE ddx_ccx_h(P,H1,H2,ANS)
739C
740C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
741C By the MadGraph Development Team
742C Please visit us at https://launchpad.net/madgraph5
743C
744C MadGraph StandAlone Version
745C
746C Returns amplitude squared summed/avg over colors
747C and helicities
748C for the point in phase space P(0:3,NEXTERNAL)
749C
750C Process: d d~ > c c~ h WEIGHTED=6
751C * Decay: h > ta+ ta- WEIGHTED=2
752C
753 IMPLICIT NONE
754C
755C CONSTANTS
756C
757 INTEGER NEXTERNAL
758 parameter(nexternal=6)
759 INTEGER NCOMB
760 parameter( ncomb=64)
761C
762C ARGUMENTS
763C
764 real*8 p(0:3,nexternal),ans
765 INTEGER H1,H2
766C
767C LOCAL VARIABLES
768C
769 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
770 real*8 t
771 real*8 matrix_ddx_ccx_h
772 INTEGER IHEL,IDEN, I
773 INTEGER JC(NEXTERNAL)
774 LOGICAL GOODHEL(NCOMB)
775 DATA ntry/0/
776 DATA goodhel/ncomb*.false./
777 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
778 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
779 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
780 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
781 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
782 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
783 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
784 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
785 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
786 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
787 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
788 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
789 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
790 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
791 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
792 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
793 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
794 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
795 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
796 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
797 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
798 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
799 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
800 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
801 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
802 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
803 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
804 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
805 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
806 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
807 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
808 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
809 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
810 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
811 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
812 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
813 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
814 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
815 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
816 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
817 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
818 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
819 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
820 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
821 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
822 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
823 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
824 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
825 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
826 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
827 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
828 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
829 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
830 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
831 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
832 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
833 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
834 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
835 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
836 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
837 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
838 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
839 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
840 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
841 DATA iden/36/
842C ----------
843C BEGIN CODE
844C ----------
845 DO ihel=1,nexternal
846 jc(ihel) = +1
847 ENDDO
848 ans = 0d0
849 DO ihel=1,ncomb
850 t=matrix_ddx_ccx_h(p ,h1,h2,nhel(1,ihel),jc(1))
851 ans=ans+t
852 ENDDO
853 ans=ans/dble(iden)
854 END
855
856
857 REAL*8 FUNCTION matrix_ddx_ccx_h(P,H1,H2,NHEL,IC)
858C
859C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
860C By the MadGraph Development Team
861C Please visit us at https://launchpad.net/madgraph5
862C
863C Returns amplitude squared summed/avg over colors
864C for the point with external lines W(0:6,NEXTERNAL)
865C
866C Process: d d~ > c c~ h WEIGHTED=6
867C * Decay: h > ta+ ta- WEIGHTED=2
868C
869 IMPLICIT NONE
870C
871C CONSTANTS
872C
873 INTEGER NGRAPHS
874 parameter(ngraphs=2)
875 INTEGER NEXTERNAL
876 parameter(nexternal=6)
877 INTEGER NWAVEFUNCS, NCOLOR
878 parameter(nwavefuncs=7, ncolor=2)
879 real*8 zero
880 parameter(zero=0d0)
881 COMPLEX*16 IMAG1
882 parameter(imag1=(0d0,1d0))
883C
884C ARGUMENTS
885C
886 real*8 p(0:3,nexternal)
887 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
888C
889C LOCAL VARIABLES
890C
891 INTEGER I,J
892 COMPLEX*16 ZTEMP
893 real*8 denom(ncolor), cf(ncolor,ncolor)
894 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
895 COMPLEX*16 W(18,NWAVEFUNCS)
896 COMPLEX*16 DUM0,DUM1
897 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
898C
899C GLOBAL VARIABLES
900C
901 include 'coupl.inc'
902C
903C COLOR DATA
904C
905 DATA denom(1)/1/
906 DATA (cf(i, 1),i= 1, 2) / 9, 3/
907C 1 T(2,1) T(3,4)
908 DATA denom(2)/1/
909 DATA (cf(i, 2),i= 1, 2) / 3, 9/
910C 1 T(2,4) T(3,1)
911
912 INTEGER H1,H2
913 real*8 matrix
914 matrix_ddx_ccx_h=0.d0
915 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
916 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
917
918C ----------
919C BEGIN CODE
920C ----------
921 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
922 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
923 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
924 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
925 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
926 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
927 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
928 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,6))
929 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
930C Amplitude(s) for diagram number 1
931 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
932 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,5))
933 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,3))
934C Amplitude(s) for diagram number 2
935 CALL vvs1_0(w(1,3),w(1,5),w(1,7),gc_72,amp(2))
936 jamp(1)=-amp(1)
937 jamp(2)=+amp(2)
938
939 matrix = 0.d0
940 DO i = 1, ncolor
941 ztemp = (0.d0,0.d0)
942 DO j = 1, ncolor
943 ztemp = ztemp + cf(j,i)*jamp(j)
944 ENDDO
945 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
946 ENDDO
947 matrix_ddx_ccx_h=matrix
948 ENDIF
949 ENDIF
950 END
951
952C ----- begin subprocesses DDX->UCX with Higgs->tautau
953
954 SUBROUTINE ddx_ucx_h(P,H1,H2,ANS)
955C
956C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
957C By the MadGraph Development Team
958C Please visit us at https://launchpad.net/madgraph5
959C
960C MadGraph StandAlone Version
961C
962C Returns amplitude squared summed/avg over colors
963C and helicities
964C for the point in phase space P(0:3,NEXTERNAL)
965C
966C Process: d d~ > u c~ h WEIGHTED=6
967C * Decay: h > ta+ ta- WEIGHTED=2
968C
969 IMPLICIT NONE
970C
971C CONSTANTS
972C
973 INTEGER NEXTERNAL
974 parameter(nexternal=6)
975 INTEGER NCOMB
976 parameter( ncomb=64)
977C
978C ARGUMENTS
979C
980 real*8 p(0:3,nexternal),ans
981 INTEGER H1,H2
982C
983C LOCAL VARIABLES
984C
985 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
986 real*8 t
987 real*8 matrix_ddx_ucx_h
988 INTEGER IHEL,IDEN, I
989 INTEGER JC(NEXTERNAL)
990 LOGICAL GOODHEL(NCOMB)
991 DATA ntry/0/
992 DATA goodhel/ncomb*.false./
993 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
994 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
995 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
996 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
997 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
998 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
999 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1000 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1001 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1002 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1003 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1004 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1005 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1006 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1007 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1008 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1009 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1010 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1011 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1012 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1013 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1014 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1015 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1016 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1017 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1018 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1019 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1020 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1021 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1022 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1023 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1024 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1025 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1026 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1027 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1028 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1029 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1030 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1031 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1032 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1033 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1034 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1035 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1036 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1037 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1038 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1039 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1040 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1041 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1042 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1043 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1044 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1045 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1046 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1047 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1048 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1049 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1050 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1051 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1052 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1053 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1054 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1055 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1056 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1057 DATA iden/36/
1058C ----------
1059C BEGIN CODE
1060C ----------
1061 DO ihel=1,nexternal
1062 jc(ihel) = +1
1063 ENDDO
1064 ans = 0d0
1065 DO ihel=1,ncomb
1066 t=matrix_ddx_ucx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1067 ans=ans+t
1068 ENDDO
1069 ans=ans/dble(iden)
1070 END
1071
1072
1073 REAL*8 FUNCTION matrix_ddx_ucx_h(P,H1,H2,NHEL,IC)
1074C
1075C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1076C By the MadGraph Development Team
1077C Please visit us at https://launchpad.net/madgraph5
1078C
1079C Returns amplitude squared summed/avg over colors
1080C for the point with external lines W(0:6,NEXTERNAL)
1081C
1082C Process: d d~ > u c~ h WEIGHTED=6
1083C * Decay: h > ta+ ta- WEIGHTED=2
1084C
1085 IMPLICIT NONE
1086C
1087C CONSTANTS
1088C
1089 INTEGER NGRAPHS
1090 parameter(ngraphs=1)
1091 INTEGER NEXTERNAL
1092 parameter(nexternal=6)
1093 INTEGER NWAVEFUNCS, NCOLOR
1094 parameter(nwavefuncs=7, ncolor=1)
1095 real*8 zero
1096 parameter(zero=0d0)
1097 COMPLEX*16 IMAG1
1098 parameter(imag1=(0d0,1d0))
1099C
1100C ARGUMENTS
1101C
1102 real*8 p(0:3,nexternal)
1103 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1104C
1105C LOCAL VARIABLES
1106C
1107 INTEGER I,J
1108 COMPLEX*16 ZTEMP
1109 real*8 denom(ncolor), cf(ncolor,ncolor)
1110 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1111 COMPLEX*16 W(18,NWAVEFUNCS)
1112 COMPLEX*16 DUM0,DUM1
1113 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1114C
1115C GLOBAL VARIABLES
1116C
1117 include 'coupl.inc'
1118C
1119C COLOR DATA
1120C
1121 DATA denom(1)/1/
1122 DATA (cf(i, 1),i= 1, 1) / 9/
1123C 1 T(2,4) T(3,1)
1124
1125 INTEGER H1,H2
1126 real*8 matrix
1127 matrix_ddx_ucx_h=0.d0
1128 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1129 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1130
1131C ----------
1132C BEGIN CODE
1133C ----------
1134 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1135 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1136 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1137 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1138 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1139 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1140 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1141 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
1142 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,3))
1143C Amplitude(s) for diagram number 1
1144 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1145 jamp(1)=+amp(1)
1146
1147 matrix = 0.d0
1148 DO i = 1, ncolor
1149 ztemp = (0.d0,0.d0)
1150 DO j = 1, ncolor
1151 ztemp = ztemp + cf(j,i)*jamp(j)
1152 ENDDO
1153 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1154 ENDDO
1155 matrix_ddx_ucx_h=matrix
1156 ENDIF
1157 ENDIF
1158 END
1159
1160C ------- BEFFINING SUBPROCESS DDX->CUX with H->tautau
1161
1162 SUBROUTINE ddx_cux_h(P,H1,H2,ANS)
1163C
1164C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1165C By the MadGraph Development Team
1166C Please visit us at https://launchpad.net/madgraph5
1167C
1168C MadGraph StandAlone Version
1169C
1170C Returns amplitude squared summed/avg over colors
1171C and helicities
1172C for the point in phase space P(0:3,NEXTERNAL)
1173C
1174C Process: d d~ > c u~ h WEIGHTED=6
1175C * Decay: h > ta+ ta- WEIGHTED=2
1176C
1177 IMPLICIT NONE
1178C
1179C CONSTANTS
1180C
1181 INTEGER NEXTERNAL
1182 parameter(nexternal=6)
1183 INTEGER NCOMB
1184 parameter( ncomb=64)
1185C
1186C ARGUMENTS
1187C
1188 real*8 p(0:3,nexternal),ans
1189 INTEGER H1,H2
1190C
1191C LOCAL VARIABLES
1192C
1193 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1194 real*8 t
1195 real*8 matrix_ddx_cux_h
1196 INTEGER IHEL,IDEN, I
1197 INTEGER JC(NEXTERNAL)
1198 LOGICAL GOODHEL(NCOMB)
1199 DATA ntry/0/
1200 DATA goodhel/ncomb*.false./
1201 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1202 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1203 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1204 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1205 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1206 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1207 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1208 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1209 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1210 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1211 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1212 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1213 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1214 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1215 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1216 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1217 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1218 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1219 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1220 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1221 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1222 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1223 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1224 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1225 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1226 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1227 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1228 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1229 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1230 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1231 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1232 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1233 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1234 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1235 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1236 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1237 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1238 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1239 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1240 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1241 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1242 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1243 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1244 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1245 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1246 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1247 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1248 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1249 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1250 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1251 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1252 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1253 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1254 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1255 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1256 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1257 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1258 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1259 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1260 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1261 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1262 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1263 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1264 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1265 DATA iden/36/
1266C ----------
1267C BEGIN CODE
1268C ----------
1269 DO ihel=1,nexternal
1270 jc(ihel) = +1
1271 ENDDO
1272 ans = 0d0
1273 DO ihel=1,ncomb
1274 t=matrix_ddx_cux_h(p ,h1,h2,nhel(1,ihel),jc(1))
1275 ans=ans+t
1276 ENDDO
1277 ans=ans/dble(iden)
1278 END
1279
1280
1281 REAL*8 FUNCTION matrix_ddx_cux_h(P,H1,H2,NHEL,IC)
1282C
1283C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1284C By the MadGraph Development Team
1285C Please visit us at https://launchpad.net/madgraph5
1286C
1287C Returns amplitude squared summed/avg over colors
1288C for the point with external lines W(0:6,NEXTERNAL)
1289C
1290C Process: d d~ > c u~ h WEIGHTED=6
1291C * Decay: h > ta+ ta- WEIGHTED=2
1292C
1293 IMPLICIT NONE
1294C
1295C CONSTANTS
1296C
1297 INTEGER NGRAPHS
1298 parameter(ngraphs=1)
1299 INTEGER NEXTERNAL
1300 parameter(nexternal=6)
1301 INTEGER NWAVEFUNCS, NCOLOR
1302 parameter(nwavefuncs=7, ncolor=1)
1303 real*8 zero
1304 parameter(zero=0d0)
1305 COMPLEX*16 IMAG1
1306 parameter(imag1=(0d0,1d0))
1307C
1308C ARGUMENTS
1309C
1310 real*8 p(0:3,nexternal)
1311 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1312C
1313C LOCAL VARIABLES
1314C
1315 INTEGER I,J
1316 COMPLEX*16 ZTEMP
1317 real*8 denom(ncolor), cf(ncolor,ncolor)
1318 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1319 COMPLEX*16 W(18,NWAVEFUNCS)
1320 COMPLEX*16 DUM0,DUM1
1321 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1322C
1323C GLOBAL VARIABLES
1324C
1325 include 'coupl.inc'
1326C
1327C COLOR DATA
1328C
1329 DATA denom(1)/1/
1330 DATA (cf(i, 1),i= 1, 1) / 9/
1331C 1 T(2,4) T(3,1)
1332
1333 INTEGER H1,H2
1334 real*8 matrix
1335 matrix_ddx_cux_h=0.d0
1336 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1337 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1338
1339C ----------
1340C BEGIN CODE
1341C ----------
1342 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1343 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1344 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1345 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1346 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1347 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1348 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1349 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,6))
1350 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
1351C Amplitude(s) for diagram number 1
1352 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1353 jamp(1)=+amp(1)
1354
1355 matrix = 0.d0
1356 DO i = 1, ncolor
1357 ztemp = (0.d0,0.d0)
1358 DO j = 1, ncolor
1359 ztemp = ztemp + cf(j,i)*jamp(j)
1360 ENDDO
1361 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1362 ENDDO
1363 matrix_ddx_cux_h=matrix
1364 ENDIF
1365 ENDIF
1366 END
1367
1368CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
1369C ----- begin subprocesses ud->jjtautau, no Higgs
1370
1371C-----------SUBPROCESS DDX->DDX NO HIGGS
1372
1373 SUBROUTINE ddx_ddx_noh(P,H1,H2,ANS)
1374C
1375C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1376C By the MadGraph Development Team
1377C Please visit us at https://launchpad.net/madgraph5
1378C
1379C MadGraph StandAlone Version
1380C
1381C Returns amplitude squared summed/avg over colors
1382C and helicities
1383C for the point in phase space P(0:3,NEXTERNAL)
1384C
1385C Process: d d~ > d d~ ta+ ta- / h QED=4
1386C
1387 IMPLICIT NONE
1388C
1389C CONSTANTS
1390C
1391 INTEGER NEXTERNAL
1392 parameter(nexternal=6)
1393 INTEGER NCOMB
1394 parameter( ncomb=64)
1395C
1396C ARGUMENTS
1397C
1398 real*8 p(0:3,nexternal),ans
1399 INTEGER H1,H2
1400C
1401C LOCAL VARIABLES
1402C
1403 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1404 real*8 t
1405 real*8 matrix_ddx_ddx_noh
1406 INTEGER IHEL,IDEN, I
1407 INTEGER JC(NEXTERNAL)
1408 LOGICAL GOODHEL(NCOMB)
1409 DATA ntry/0/
1410 DATA goodhel/ncomb*.false./
1411 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1412 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1413 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1414 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1415 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1416 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1417 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1418 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1419 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1420 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1421 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1422 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1423 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1424 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1425 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1426 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1427 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1428 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1429 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1430 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1431 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1432 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1433 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1434 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1435 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1436 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1437 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1438 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1439 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1440 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1441 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1442 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1443 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1444 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1445 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1446 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1447 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1448 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1449 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1450 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1451 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1452 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1453 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1454 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1455 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1456 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1457 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1458 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1459 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1460 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1461 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1462 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1463 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1464 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1465 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1466 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1467 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1468 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1469 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1470 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1471 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1472 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1473 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1474 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1475 DATA iden/36/
1476C ----------
1477C BEGIN CODE
1478C ----------
1479 DO ihel=1,nexternal
1480 jc(ihel) = +1
1481 ENDDO
1482 ans = 0d0
1483 DO ihel=1,ncomb
1484 t=matrix_ddx_ddx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1485 ans=ans+t
1486 ENDDO
1487 ans=ans/dble(iden)
1488 END
1489
1490
1491 REAL*8 FUNCTION matrix_ddx_ddx_noh(P,H1,H2,NHEL,IC)
1492C
1493C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1494C By the MadGraph Development Team
1495C Please visit us at https://launchpad.net/madgraph5
1496C
1497C Returns amplitude squared summed/avg over colors
1498C for the point with external lines W(0:6,NEXTERNAL)
1499C
1500C Process: d d~ > d d~ ta+ ta- / h QED=4
1501C
1502 IMPLICIT NONE
1503C
1504C CONSTANTS
1505C
1506 INTEGER NGRAPHS
1507 parameter(ngraphs=64)
1508 INTEGER NEXTERNAL
1509 parameter(nexternal=6)
1510 INTEGER NWAVEFUNCS, NCOLOR
1511 parameter(nwavefuncs=16, ncolor=2)
1512 real*8 zero
1513 parameter(zero=0d0)
1514 COMPLEX*16 IMAG1
1515 parameter(imag1=(0d0,1d0))
1516C
1517C ARGUMENTS
1518C
1519 real*8 p(0:3,nexternal)
1520 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1521C
1522C LOCAL VARIABLES
1523C
1524 INTEGER I,J
1525 COMPLEX*16 ZTEMP
1526 real*8 denom(ncolor), cf(ncolor,ncolor)
1527 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1528 COMPLEX*16 W(18,NWAVEFUNCS)
1529 COMPLEX*16 DUM0,DUM1
1530 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1531C
1532C GLOBAL VARIABLES
1533C
1534 include 'coupl.inc'
1535C
1536C COLOR DATA
1537C
1538 DATA denom(1)/1/
1539 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1540C 1 T(2,1) T(3,4)
1541 DATA denom(2)/1/
1542 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1543C 1 T(2,4) T(3,1)
1544
1545 INTEGER H1,H2
1546 real*8 matrix
1547 matrix_ddx_ddx_noh=0.d0
1548 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1549 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1550
1551C ----------
1552C BEGIN CODE
1553C ----------
1554 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1555 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1556 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1557 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1558 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1559 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1560 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
1561 CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
1562 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1563C Amplitude(s) for diagram number 1
1564 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1565 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1566C Amplitude(s) for diagram number 2
1567 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1568 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
1569C Amplitude(s) for diagram number 3
1570 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1571C Amplitude(s) for diagram number 4
1572 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1573 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
1574 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1575C Amplitude(s) for diagram number 5
1576 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1577 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1578C Amplitude(s) for diagram number 6
1579 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1580C Amplitude(s) for diagram number 7
1581 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1582C Amplitude(s) for diagram number 8
1583 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1584 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1585 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
1586C Amplitude(s) for diagram number 9
1587 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
1588 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1589C Amplitude(s) for diagram number 10
1590 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
1591 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1592C Amplitude(s) for diagram number 11
1593 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
1594C Amplitude(s) for diagram number 12
1595 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
1596 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1597 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1598C Amplitude(s) for diagram number 13
1599 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
1600 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1601C Amplitude(s) for diagram number 14
1602 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(14))
1603C Amplitude(s) for diagram number 15
1604 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
1605C Amplitude(s) for diagram number 16
1606 CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(16))
1607 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
1608C Amplitude(s) for diagram number 17
1609 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_1,amp(17))
1610 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1611C Amplitude(s) for diagram number 18
1612 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
1613C Amplitude(s) for diagram number 19
1614 CALL ffv2_3_0(w(1,4),w(1,14),w(1,7),gc_50,gc_58,amp(19))
1615C Amplitude(s) for diagram number 20
1616 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
1617 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,9))
1618 CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,14))
1619 CALL ffv1_2(w(1,5),w(1,9),gc_3,mta,zero,w(1,10))
1620C Amplitude(s) for diagram number 21
1621 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(21))
1622 CALL ffv1_1(w(1,6),w(1,9),gc_3,mta,zero,w(1,13))
1623C Amplitude(s) for diagram number 22
1624 CALL ffv1_0(w(1,5),w(1,13),w(1,14),gc_3,amp(22))
1625 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,15))
1626C Amplitude(s) for diagram number 23
1627 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(23))
1628C Amplitude(s) for diagram number 24
1629 CALL ffv2_4_0(w(1,5),w(1,13),w(1,15),gc_50,gc_59,amp(24))
1630 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,13))
1631 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,10))
1632C Amplitude(s) for diagram number 25
1633 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(25))
1634 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,16))
1635C Amplitude(s) for diagram number 26
1636 CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
1637C Amplitude(s) for diagram number 27
1638 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(27))
1639C Amplitude(s) for diagram number 28
1640 CALL ffv2_4_0(w(1,5),w(1,16),w(1,15),gc_50,gc_59,amp(28))
1641 CALL ffv1_1(w(1,2),w(1,9),gc_1,zero,zero,w(1,16))
1642C Amplitude(s) for diagram number 29
1643 CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_1,amp(29))
1644 CALL ffv1_2(w(1,4),w(1,9),gc_1,zero,zero,w(1,5))
1645C Amplitude(s) for diagram number 30
1646 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(30))
1647C Amplitude(s) for diagram number 31
1648 CALL ffv2_3_0(w(1,4),w(1,16),w(1,7),gc_50,gc_58,amp(31))
1649C Amplitude(s) for diagram number 32
1650 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(32))
1651 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1652 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,16))
1653C Amplitude(s) for diagram number 33
1654 CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_1,amp(33))
1655 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,9))
1656C Amplitude(s) for diagram number 34
1657 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(34))
1658C Amplitude(s) for diagram number 35
1659 CALL ffv2_3_0(w(1,4),w(1,16),w(1,7),gc_50,gc_58,amp(35))
1660C Amplitude(s) for diagram number 36
1661 CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(36))
1662 CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,9))
1663C Amplitude(s) for diagram number 37
1664 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(37))
1665 CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,16))
1666C Amplitude(s) for diagram number 38
1667 CALL ffv1_0(w(1,16),w(1,2),w(1,12),gc_1,amp(38))
1668C Amplitude(s) for diagram number 39
1669 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(39))
1670C Amplitude(s) for diagram number 40
1671 CALL ffv2_3_0(w(1,16),w(1,2),w(1,7),gc_50,gc_58,amp(40))
1672 CALL ffv1_2(w(1,1),w(1,14),gc_1,zero,zero,w(1,16))
1673C Amplitude(s) for diagram number 41
1674 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_1,amp(41))
1675 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,9))
1676C Amplitude(s) for diagram number 42
1677 CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_1,amp(42))
1678C Amplitude(s) for diagram number 43
1679 CALL ffv2_3_0(w(1,16),w(1,3),w(1,7),gc_50,gc_58,amp(43))
1680 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,16))
1681C Amplitude(s) for diagram number 44
1682 CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_1,amp(44))
1683 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,14))
1684 CALL ffv1_2(w(1,1),w(1,14),gc_11,zero,zero,w(1,13))
1685C Amplitude(s) for diagram number 45
1686 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(45))
1687C Amplitude(s) for diagram number 46
1688 CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_11,amp(46))
1689C Amplitude(s) for diagram number 47
1690 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(47))
1691C Amplitude(s) for diagram number 48
1692 CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_11,amp(48))
1693 CALL ffv2_3_2(w(1,1),w(1,15),gc_50,gc_58,zero,zero,w(1,14))
1694C Amplitude(s) for diagram number 49
1695 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(49))
1696C Amplitude(s) for diagram number 50
1697 CALL ffv2_3_0(w(1,9),w(1,3),w(1,15),gc_50,gc_58,amp(50))
1698C Amplitude(s) for diagram number 51
1699 CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(51))
1700C Amplitude(s) for diagram number 52
1701 CALL ffv2_3_0(w(1,16),w(1,3),w(1,15),gc_50,gc_58,amp(52))
1702 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,15))
1703C Amplitude(s) for diagram number 53
1704 CALL ffv1_0(w(1,15),w(1,2),w(1,12),gc_1,amp(53))
1705C Amplitude(s) for diagram number 54
1706 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_1,amp(54))
1707C Amplitude(s) for diagram number 55
1708 CALL ffv2_3_0(w(1,15),w(1,2),w(1,7),gc_50,gc_58,amp(55))
1709C Amplitude(s) for diagram number 56
1710 CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_1,amp(56))
1711 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
1712 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1713C Amplitude(s) for diagram number 57
1714 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(57))
1715C Amplitude(s) for diagram number 58
1716 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(58))
1717C Amplitude(s) for diagram number 59
1718 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(59))
1719C Amplitude(s) for diagram number 60
1720 CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_11,amp(60))
1721 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
1722C Amplitude(s) for diagram number 61
1723 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(61))
1724C Amplitude(s) for diagram number 62
1725 CALL ffv2_3_0(w(1,9),w(1,2),w(1,11),gc_50,gc_58,amp(62))
1726C Amplitude(s) for diagram number 63
1727 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(63))
1728C Amplitude(s) for diagram number 64
1729 CALL ffv2_3_0(w(1,16),w(1,2),w(1,11),gc_50,gc_58,amp(64))
1730 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1731 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1732 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1733 $ +1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)+1./2.*amp(36)
1734 $ +1./2.*amp(45)+1./2.*amp(46)+1./2.*amp(47)+1./2.*amp(48)
1735 $ -amp(53)-amp(54)-amp(55)-amp(56)+1./6.*amp(57)+1./6.*amp(58)
1736 $ +1./6.*amp(59)+1./6.*amp(60)-amp(61)-amp(62)-amp(63)-amp(64)
1737 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1738 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1739 $ +amp(28)+amp(29)+amp(30)+amp(31)+amp(32)-1./6.*amp(33)
1740 $ -1./6.*amp(34)-1./6.*amp(35)-1./6.*amp(36)+amp(37)+amp(38)
1741 $ +amp(39)+amp(40)+amp(41)+amp(42)+amp(43)+amp(44)-1./6.*amp(45)
1742 $ -1./6.*amp(46)-1./6.*amp(47)-1./6.*amp(48)+amp(49)+amp(50)
1743 $ +amp(51)+amp(52)-1./2.*amp(57)-1./2.*amp(58)-1./2.*amp(59)
1744 $ -1./2.*amp(60)
1745
1746 matrix = 0.d0
1747 DO i = 1, ncolor
1748 ztemp = (0.d0,0.d0)
1749 DO j = 1, ncolor
1750 ztemp = ztemp + cf(j,i)*jamp(j)
1751 ENDDO
1752 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1753 ENDDO
1754 matrix_ddx_ddx_noh=matrix
1755 ENDIF
1756 ENDIF
1757 END
1758
1759C-----------SUBPROCESS DDX->UUX NO HIGGS
1760
1761 SUBROUTINE ddx_uux_noh(P,H1,H2,ANS)
1762C
1763C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1764C By the MadGraph Development Team
1765C Please visit us at https://launchpad.net/madgraph5
1766C
1767C MadGraph StandAlone Version
1768C
1769C Returns amplitude squared summed/avg over colors
1770C and helicities
1771C for the point in phase space P(0:3,NEXTERNAL)
1772C
1773C Process: d d~ > u u~ ta+ ta- / h QED=4
1774C
1775 IMPLICIT NONE
1776C
1777C CONSTANTS
1778C
1779 INTEGER NEXTERNAL
1780 parameter(nexternal=6)
1781 INTEGER NCOMB
1782 parameter( ncomb=64)
1783C
1784C ARGUMENTS
1785C
1786 real*8 p(0:3,nexternal),ans
1787 INTEGER H1,H2
1788C
1789C LOCAL VARIABLES
1790C
1791 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1792 real*8 t
1793 real*8 matrix_ddx_uux_noh
1794 INTEGER IHEL,IDEN, I
1795 INTEGER JC(NEXTERNAL)
1796 LOGICAL GOODHEL(NCOMB)
1797 DATA ntry/0/
1798 DATA goodhel/ncomb*.false./
1799 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1800 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1801 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1802 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1803 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1804 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1805 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1806 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1807 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1808 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1809 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1810 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1811 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1812 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1813 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1814 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1815 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1816 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1817 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1818 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1819 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1820 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1821 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1822 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1823 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1824 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1825 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1826 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1827 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1828 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1829 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1830 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1831 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1832 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1833 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1834 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1835 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1836 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1837 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1838 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1839 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1840 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1841 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1842 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1843 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1844 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1845 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1846 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1847 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1848 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1849 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1850 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1851 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1852 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1853 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1854 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1855 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1856 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1857 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1858 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1859 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1860 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1861 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1862 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1863 DATA iden/36/
1864C ----------
1865C BEGIN CODE
1866C ----------
1867 DO ihel=1,nexternal
1868 jc(ihel) = +1
1869 ENDDO
1870 ans = 0d0
1871 DO ihel=1,ncomb
1872 t=matrix_ddx_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1873 ans=ans+t
1874 ENDDO
1875 ans=ans/dble(iden)
1876 END
1877
1878
1879 REAL*8 FUNCTION matrix_ddx_uux_noh(P,H1,H2,NHEL,IC)
1880C
1881C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
1882C By the MadGraph Development Team
1883C Please visit us at https://launchpad.net/madgraph5
1884C
1885C Returns amplitude squared summed/avg over colors
1886C for the point with external lines W(0:6,NEXTERNAL)
1887C
1888C Process: d d~ > u u~ ta+ ta- / h QED=4
1889C
1890 IMPLICIT NONE
1891C
1892C CONSTANTS
1893C
1894 INTEGER NGRAPHS
1895 parameter(ngraphs=43)
1896 INTEGER NEXTERNAL
1897 parameter(nexternal=6)
1898 INTEGER NWAVEFUNCS, NCOLOR
1899 parameter(nwavefuncs=14, ncolor=2)
1900 real*8 zero
1901 parameter(zero=0d0)
1902 COMPLEX*16 IMAG1
1903 parameter(imag1=(0d0,1d0))
1904C
1905C ARGUMENTS
1906C
1907 real*8 p(0:3,nexternal)
1908 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1909C
1910C LOCAL VARIABLES
1911C
1912 INTEGER I,J
1913 COMPLEX*16 ZTEMP
1914 real*8 denom(ncolor), cf(ncolor,ncolor)
1915 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1916 COMPLEX*16 W(18,NWAVEFUNCS)
1917 COMPLEX*16 DUM0,DUM1
1918 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1919C
1920C GLOBAL VARIABLES
1921C
1922 include 'coupl.inc'
1923C
1924C COLOR DATA
1925C
1926 DATA denom(1)/1/
1927 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1928C 1 T(2,1) T(3,4)
1929 DATA denom(2)/1/
1930 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1931C 1 T(2,4) T(3,1)
1932
1933 INTEGER H1,H2
1934 real*8 matrix
1935 matrix_ddx_uux_noh=0.d0
1936 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1937 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1938
1939C ----------
1940C BEGIN CODE
1941C ----------
1942 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1943 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1944 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1945 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1946 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1947 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1948 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
1949 CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
1950 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1951C Amplitude(s) for diagram number 1
1952 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1953 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1954C Amplitude(s) for diagram number 2
1955 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1956 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
1957C Amplitude(s) for diagram number 3
1958 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1959C Amplitude(s) for diagram number 4
1960 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1961 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
1962 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1963C Amplitude(s) for diagram number 5
1964 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1965 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1966C Amplitude(s) for diagram number 6
1967 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1968C Amplitude(s) for diagram number 7
1969 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1970C Amplitude(s) for diagram number 8
1971 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1972 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1973 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
1974C Amplitude(s) for diagram number 9
1975 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
1976 CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
1977C Amplitude(s) for diagram number 10
1978 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
1979 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1980C Amplitude(s) for diagram number 11
1981 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
1982C Amplitude(s) for diagram number 12
1983 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
1984 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1985 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1986C Amplitude(s) for diagram number 13
1987 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
1988 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1989C Amplitude(s) for diagram number 14
1990 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
1991C Amplitude(s) for diagram number 15
1992 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
1993C Amplitude(s) for diagram number 16
1994 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
1995 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
1996C Amplitude(s) for diagram number 17
1997 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
1998 CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
1999C Amplitude(s) for diagram number 18
2000 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
2001C Amplitude(s) for diagram number 19
2002 CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
2003C Amplitude(s) for diagram number 20
2004 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
2005 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,9))
2006 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,14))
2007 CALL ffv2_1(w(1,6),w(1,9),gc_108,zero,zero,w(1,10))
2008C Amplitude(s) for diagram number 21
2009 CALL ffv2_0(w(1,5),w(1,10),w(1,14),gc_108,amp(21))
2010C Amplitude(s) for diagram number 22
2011 CALL vvv1_0(w(1,12),w(1,14),w(1,9),gc_4,amp(22))
2012C Amplitude(s) for diagram number 23
2013 CALL vvv1_0(w(1,14),w(1,9),w(1,7),gc_53,amp(23))
2014 CALL ffv2_1(w(1,2),w(1,9),gc_100,zero,zero,w(1,10))
2015C Amplitude(s) for diagram number 24
2016 CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_2,amp(24))
2017 CALL ffv2_2(w(1,4),w(1,9),gc_100,zero,zero,w(1,5))
2018C Amplitude(s) for diagram number 25
2019 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(25))
2020C Amplitude(s) for diagram number 26
2021 CALL ffv2_5_0(w(1,4),w(1,10),w(1,7),gc_51,gc_58,amp(26))
2022C Amplitude(s) for diagram number 27
2023 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2024 CALL ffv2_2(w(1,1),w(1,14),gc_100,zero,zero,w(1,5))
2025C Amplitude(s) for diagram number 28
2026 CALL ffv1_0(w(1,5),w(1,3),w(1,12),gc_2,amp(28))
2027 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,10))
2028C Amplitude(s) for diagram number 29
2029 CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_100,amp(29))
2030C Amplitude(s) for diagram number 30
2031 CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(30))
2032 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,5))
2033C Amplitude(s) for diagram number 31
2034 CALL ffv2_0(w(1,5),w(1,3),w(1,14),gc_100,amp(31))
2035 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,14))
2036C Amplitude(s) for diagram number 32
2037 CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_1,amp(32))
2038C Amplitude(s) for diagram number 33
2039 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(33))
2040C Amplitude(s) for diagram number 34
2041 CALL ffv2_3_0(w(1,14),w(1,2),w(1,7),gc_50,gc_58,amp(34))
2042C Amplitude(s) for diagram number 35
2043 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(35))
2044 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2045 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2046C Amplitude(s) for diagram number 36
2047 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(36))
2048C Amplitude(s) for diagram number 37
2049 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
2050C Amplitude(s) for diagram number 38
2051 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(38))
2052C Amplitude(s) for diagram number 39
2053 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(39))
2054 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2055C Amplitude(s) for diagram number 40
2056 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(40))
2057C Amplitude(s) for diagram number 41
2058 CALL ffv2_3_0(w(1,10),w(1,2),w(1,11),gc_50,gc_58,amp(41))
2059C Amplitude(s) for diagram number 42
2060 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(42))
2061C Amplitude(s) for diagram number 43
2062 CALL ffv2_3_0(w(1,5),w(1,2),w(1,11),gc_50,gc_58,amp(43))
2063 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2064 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2065 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2066 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
2067 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
2068 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2069 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2070 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
2071 $ -1./2.*amp(38)-1./2.*amp(39)
2072
2073 matrix = 0.d0
2074 DO i = 1, ncolor
2075 ztemp = (0.d0,0.d0)
2076 DO j = 1, ncolor
2077 ztemp = ztemp + cf(j,i)*jamp(j)
2078 ENDDO
2079 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2080 ENDDO
2081 matrix_ddx_uux_noh=matrix
2082 ENDIF
2083 ENDIF
2084 END
2085
2086C-----------SUBPROCESS DDX->SSX NO HIGGS
2087
2088 SUBROUTINE ddx_ssx_noh(P,H1,H2,ANS)
2089C
2090C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2091C By the MadGraph Development Team
2092C Please visit us at https://launchpad.net/madgraph5
2093C
2094C MadGraph StandAlone Version
2095C
2096C Returns amplitude squared summed/avg over colors
2097C and helicities
2098C for the point in phase space P(0:3,NEXTERNAL)
2099C
2100C Process: d d~ > s s~ ta+ ta- / h QED=4
2101C
2102 IMPLICIT NONE
2103C
2104C CONSTANTS
2105C
2106 INTEGER NEXTERNAL
2107 parameter(nexternal=6)
2108 INTEGER NCOMB
2109 parameter( ncomb=64)
2110C
2111C ARGUMENTS
2112C
2113 real*8 p(0:3,nexternal),ans
2114 INTEGER H1,H2
2115C
2116C LOCAL VARIABLES
2117C
2118 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2119 real*8 t
2120 real*8 matrix_ddx_ssx_noh
2121 INTEGER IHEL,IDEN, I
2122 INTEGER JC(NEXTERNAL)
2123 LOGICAL GOODHEL(NCOMB)
2124 DATA ntry/0/
2125 DATA goodhel/ncomb*.false./
2126 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2127 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2128 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2129 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2130 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2131 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2132 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2133 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2134 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2135 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2136 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2137 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2138 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2139 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2140 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2141 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2142 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2143 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2144 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2145 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2146 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2147 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2148 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2149 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2150 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2151 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2152 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2153 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2154 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2155 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2156 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2157 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2158 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2159 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2160 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2161 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2162 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2163 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2164 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2165 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2166 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2167 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2168 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2169 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2170 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2171 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2172 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2173 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2174 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2175 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2176 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2177 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2178 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2179 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2180 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2181 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2182 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2183 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2184 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2185 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2186 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2187 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2188 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2189 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2190 DATA iden/36/
2191C ----------
2192C BEGIN CODE
2193C ----------
2194 DO ihel=1,nexternal
2195 jc(ihel) = +1
2196 ENDDO
2197 ans = 0d0
2198 DO ihel=1,ncomb
2199 t=matrix_ddx_ssx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2200 ans=ans+t
2201 ENDDO
2202 ans=ans/dble(iden)
2203 END
2204
2205
2206 REAL*8 FUNCTION matrix_ddx_ssx_noh(P,H1,H2,NHEL,IC)
2207C
2208C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2209C By the MadGraph Development Team
2210C Please visit us at https://launchpad.net/madgraph5
2211C
2212C Returns amplitude squared summed/avg over colors
2213C for the point with external lines W(0:6,NEXTERNAL)
2214C
2215C Process: d d~ > s s~ ta+ ta- / h QED=4
2216C
2217 IMPLICIT NONE
2218C
2219C CONSTANTS
2220C
2221 INTEGER NGRAPHS
2222 parameter(ngraphs=32)
2223 INTEGER NEXTERNAL
2224 parameter(nexternal=6)
2225 INTEGER NWAVEFUNCS, NCOLOR
2226 parameter(nwavefuncs=13, ncolor=2)
2227 real*8 zero
2228 parameter(zero=0d0)
2229 COMPLEX*16 IMAG1
2230 parameter(imag1=(0d0,1d0))
2231C
2232C ARGUMENTS
2233C
2234 real*8 p(0:3,nexternal)
2235 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2236C
2237C LOCAL VARIABLES
2238C
2239 INTEGER I,J
2240 COMPLEX*16 ZTEMP
2241 real*8 denom(ncolor), cf(ncolor,ncolor)
2242 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2243 COMPLEX*16 W(18,NWAVEFUNCS)
2244 COMPLEX*16 DUM0,DUM1
2245 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2246C
2247C GLOBAL VARIABLES
2248C
2249 include 'coupl.inc'
2250C
2251C COLOR DATA
2252C
2253 DATA denom(1)/1/
2254 DATA (cf(i, 1),i= 1, 2) / 9, 3/
2255C 1 T(2,1) T(3,4)
2256 DATA denom(2)/1/
2257 DATA (cf(i, 2),i= 1, 2) / 3, 9/
2258C 1 T(2,4) T(3,1)
2259
2260 INTEGER H1,H2
2261 real*8 matrix
2262 matrix_ddx_ssx_noh=0.d0
2263 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2264 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2265
2266C ----------
2267C BEGIN CODE
2268C ----------
2269 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2270 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2271 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2272 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2273 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2274 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2275 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
2276 CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
2277 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2278C Amplitude(s) for diagram number 1
2279 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2280 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2281C Amplitude(s) for diagram number 2
2282 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2283 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
2284C Amplitude(s) for diagram number 3
2285 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2286C Amplitude(s) for diagram number 4
2287 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2288 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
2289 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2290C Amplitude(s) for diagram number 5
2291 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2292 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2293C Amplitude(s) for diagram number 6
2294 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2295C Amplitude(s) for diagram number 7
2296 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2297C Amplitude(s) for diagram number 8
2298 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2299 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2300 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
2301C Amplitude(s) for diagram number 9
2302 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
2303 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
2304C Amplitude(s) for diagram number 10
2305 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
2306 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2307C Amplitude(s) for diagram number 11
2308 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
2309C Amplitude(s) for diagram number 12
2310 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
2311 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2312 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2313C Amplitude(s) for diagram number 13
2314 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
2315 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
2316C Amplitude(s) for diagram number 14
2317 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_1,amp(14))
2318C Amplitude(s) for diagram number 15
2319 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
2320C Amplitude(s) for diagram number 16
2321 CALL ffv2_3_0(w(1,6),w(1,3),w(1,7),gc_50,gc_58,amp(16))
2322 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
2323C Amplitude(s) for diagram number 17
2324 CALL ffv1_0(w(1,4),w(1,6),w(1,12),gc_1,amp(17))
2325 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
2326C Amplitude(s) for diagram number 18
2327 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
2328C Amplitude(s) for diagram number 19
2329 CALL ffv2_3_0(w(1,4),w(1,6),w(1,7),gc_50,gc_58,amp(19))
2330C Amplitude(s) for diagram number 20
2331 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
2332 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
2333C Amplitude(s) for diagram number 21
2334 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(21))
2335 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
2336C Amplitude(s) for diagram number 22
2337 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_1,amp(22))
2338C Amplitude(s) for diagram number 23
2339 CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(23))
2340 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
2341C Amplitude(s) for diagram number 24
2342 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_1,amp(24))
2343 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2344 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2345C Amplitude(s) for diagram number 25
2346 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(25))
2347C Amplitude(s) for diagram number 26
2348 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(26))
2349C Amplitude(s) for diagram number 27
2350 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2351C Amplitude(s) for diagram number 28
2352 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(28))
2353 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2354C Amplitude(s) for diagram number 29
2355 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(29))
2356C Amplitude(s) for diagram number 30
2357 CALL ffv2_3_0(w(1,6),w(1,2),w(1,11),gc_50,gc_58,amp(30))
2358C Amplitude(s) for diagram number 31
2359 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(31))
2360C Amplitude(s) for diagram number 32
2361 CALL ffv2_3_0(w(1,9),w(1,2),w(1,11),gc_50,gc_58,amp(32))
2362 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2363 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2364 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2365 $ -amp(21)-amp(22)-amp(23)-amp(24)+1./6.*amp(25)+1./6.*amp(26)
2366 $ +1./6.*amp(27)+1./6.*amp(28)-amp(29)-amp(30)-amp(31)-amp(32)
2367 jamp(2)=+1./2.*(-amp(13)-amp(14)-amp(15)-amp(16)-amp(25)-amp(26)
2368 $ -amp(27)-amp(28))
2369
2370 matrix = 0.d0
2371 DO i = 1, ncolor
2372 ztemp = (0.d0,0.d0)
2373 DO j = 1, ncolor
2374 ztemp = ztemp + cf(j,i)*jamp(j)
2375 ENDDO
2376 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2377 ENDDO
2378 matrix_ddx_ssx_noh=matrix
2379 ENDIF
2380 ENDIF
2381 END
2382
2383C-----------SUBPROCESS DDX->CCX NO HIGGS
2384
2385 SUBROUTINE ddx_ccx_noh(P,H1,H2,ANS)
2386C
2387C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2388C By the MadGraph Development Team
2389C Please visit us at https://launchpad.net/madgraph5
2390C
2391C MadGraph StandAlone Version
2392C
2393C Returns amplitude squared summed/avg over colors
2394C and helicities
2395C for the point in phase space P(0:3,NEXTERNAL)
2396C
2397C Process: d d~ > c c~ ta+ ta- / h QED=4
2398C
2399 IMPLICIT NONE
2400C
2401C CONSTANTS
2402C
2403 INTEGER NEXTERNAL
2404 parameter(nexternal=6)
2405 INTEGER NCOMB
2406 parameter( ncomb=64)
2407C
2408C ARGUMENTS
2409C
2410 real*8 p(0:3,nexternal),ans
2411 INTEGER H1,H2
2412C
2413C LOCAL VARIABLES
2414C
2415 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2416 real*8 t
2417 real*8 matrix_ddx_ccx_noh
2418 INTEGER IHEL,IDEN, I
2419 INTEGER JC(NEXTERNAL)
2420 LOGICAL GOODHEL(NCOMB)
2421 DATA ntry/0/
2422 DATA goodhel/ncomb*.false./
2423 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2424 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2425 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2426 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2427 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2428 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2429 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2430 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2431 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2432 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2433 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2434 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2435 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2436 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2437 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2438 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2439 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2440 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2441 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2442 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2443 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2444 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2445 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2446 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2447 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2448 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2449 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2450 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2451 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2452 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2453 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2454 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2455 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2456 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2457 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2458 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2459 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2460 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2461 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2462 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2463 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2464 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2465 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2466 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2467 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2468 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2469 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2470 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2471 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2472 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2473 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2474 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2475 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2476 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2477 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2478 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2479 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2480 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2481 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2482 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2483 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2484 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2485 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2486 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2487 DATA iden/36/
2488C ----------
2489C BEGIN CODE
2490C ----------
2491 DO ihel=1,nexternal
2492 jc(ihel) = +1
2493 ENDDO
2494 ans = 0d0
2495 DO ihel=1,ncomb
2496 t=matrix_ddx_ccx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2497 ans=ans+t
2498 ENDDO
2499 ans=ans/dble(iden)
2500 END
2501
2502
2503 REAL*8 FUNCTION matrix_ddx_ccx_noh(P,H1,H2,NHEL,IC)
2504C
2505C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2506C By the MadGraph Development Team
2507C Please visit us at https://launchpad.net/madgraph5
2508C
2509C Returns amplitude squared summed/avg over colors
2510C for the point with external lines W(0:6,NEXTERNAL)
2511C
2512C Process: d d~ > c c~ ta+ ta- / h QED=4
2513C
2514 IMPLICIT NONE
2515C
2516C CONSTANTS
2517C
2518 INTEGER NGRAPHS
2519 parameter(ngraphs=43)
2520 INTEGER NEXTERNAL
2521 parameter(nexternal=6)
2522 INTEGER NWAVEFUNCS, NCOLOR
2523 parameter(nwavefuncs=14, ncolor=2)
2524 real*8 zero
2525 parameter(zero=0d0)
2526 COMPLEX*16 IMAG1
2527 parameter(imag1=(0d0,1d0))
2528C
2529C ARGUMENTS
2530C
2531 real*8 p(0:3,nexternal)
2532 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2533C
2534C LOCAL VARIABLES
2535C
2536 INTEGER I,J
2537 COMPLEX*16 ZTEMP
2538 real*8 denom(ncolor), cf(ncolor,ncolor)
2539 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2540 COMPLEX*16 W(18,NWAVEFUNCS)
2541 COMPLEX*16 DUM0,DUM1
2542 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2543C
2544C GLOBAL VARIABLES
2545C
2546 include 'coupl.inc'
2547C
2548C COLOR DATA
2549C
2550 DATA denom(1)/1/
2551 DATA (cf(i, 1),i= 1, 2) / 9, 3/
2552C 1 T(2,1) T(3,4)
2553 DATA denom(2)/1/
2554 DATA (cf(i, 2),i= 1, 2) / 3, 9/
2555C 1 T(2,4) T(3,1)
2556
2557 INTEGER H1,H2
2558 real*8 matrix
2559 matrix_ddx_ccx_noh=0.d0
2560 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2561 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2562
2563
2564C ----------
2565C BEGIN CODE
2566C ----------
2567 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2568 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2569 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2570 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2571 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2572 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2573 CALL ffv1p0_3(w(1,1),w(1,2),gc_1,zero,zero,w(1,7))
2574 CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
2575 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2576C Amplitude(s) for diagram number 1
2577 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2578 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2579C Amplitude(s) for diagram number 2
2580 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2581 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
2582C Amplitude(s) for diagram number 3
2583 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2584C Amplitude(s) for diagram number 4
2585 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2586 CALL ffv2_3_3(w(1,1),w(1,2),gc_50,gc_58,mz,wz,w(1,10))
2587 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2588C Amplitude(s) for diagram number 5
2589 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2590 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2591C Amplitude(s) for diagram number 6
2592 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2593C Amplitude(s) for diagram number 7
2594 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2595C Amplitude(s) for diagram number 8
2596 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2597 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2598 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
2599C Amplitude(s) for diagram number 9
2600 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
2601 CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
2602C Amplitude(s) for diagram number 10
2603 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
2604 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2605C Amplitude(s) for diagram number 11
2606 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
2607C Amplitude(s) for diagram number 12
2608 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
2609 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2610 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2611C Amplitude(s) for diagram number 13
2612 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
2613 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
2614C Amplitude(s) for diagram number 14
2615 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
2616C Amplitude(s) for diagram number 15
2617 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
2618C Amplitude(s) for diagram number 16
2619 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
2620 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
2621C Amplitude(s) for diagram number 17
2622 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
2623 CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
2624C Amplitude(s) for diagram number 18
2625 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
2626C Amplitude(s) for diagram number 19
2627 CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
2628C Amplitude(s) for diagram number 20
2629 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
2630 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,9))
2631 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,14))
2632 CALL ffv2_1(w(1,6),w(1,9),gc_108,zero,zero,w(1,10))
2633C Amplitude(s) for diagram number 21
2634 CALL ffv2_0(w(1,5),w(1,10),w(1,14),gc_108,amp(21))
2635C Amplitude(s) for diagram number 22
2636 CALL vvv1_0(w(1,12),w(1,14),w(1,9),gc_4,amp(22))
2637C Amplitude(s) for diagram number 23
2638 CALL vvv1_0(w(1,14),w(1,9),w(1,7),gc_53,amp(23))
2639 CALL ffv2_1(w(1,2),w(1,9),gc_44,zero,zero,w(1,10))
2640C Amplitude(s) for diagram number 24
2641 CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_2,amp(24))
2642 CALL ffv2_2(w(1,4),w(1,9),gc_44,zero,zero,w(1,5))
2643C Amplitude(s) for diagram number 25
2644 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_1,amp(25))
2645C Amplitude(s) for diagram number 26
2646 CALL ffv2_5_0(w(1,4),w(1,10),w(1,7),gc_51,gc_58,amp(26))
2647C Amplitude(s) for diagram number 27
2648 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(27))
2649 CALL ffv2_2(w(1,1),w(1,14),gc_44,zero,zero,w(1,5))
2650C Amplitude(s) for diagram number 28
2651 CALL ffv1_0(w(1,5),w(1,3),w(1,12),gc_2,amp(28))
2652 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,10))
2653C Amplitude(s) for diagram number 29
2654 CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_44,amp(29))
2655C Amplitude(s) for diagram number 30
2656 CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(30))
2657 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,5))
2658C Amplitude(s) for diagram number 31
2659 CALL ffv2_0(w(1,5),w(1,3),w(1,14),gc_44,amp(31))
2660 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,14))
2661C Amplitude(s) for diagram number 32
2662 CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_1,amp(32))
2663C Amplitude(s) for diagram number 33
2664 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(33))
2665C Amplitude(s) for diagram number 34
2666 CALL ffv2_3_0(w(1,14),w(1,2),w(1,7),gc_50,gc_58,amp(34))
2667C Amplitude(s) for diagram number 35
2668 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(35))
2669 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2670 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2671C Amplitude(s) for diagram number 36
2672 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_1,amp(36))
2673C Amplitude(s) for diagram number 37
2674 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
2675C Amplitude(s) for diagram number 38
2676 CALL ffv2_3_0(w(1,4),w(1,2),w(1,7),gc_50,gc_58,amp(38))
2677C Amplitude(s) for diagram number 39
2678 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_11,amp(39))
2679 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
2680C Amplitude(s) for diagram number 40
2681 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_1,amp(40))
2682C Amplitude(s) for diagram number 41
2683 CALL ffv2_3_0(w(1,10),w(1,2),w(1,11),gc_50,gc_58,amp(41))
2684C Amplitude(s) for diagram number 42
2685 CALL ffv2_3_0(w(1,8),w(1,2),w(1,7),gc_50,gc_58,amp(42))
2686C Amplitude(s) for diagram number 43
2687 CALL ffv2_3_0(w(1,5),w(1,2),w(1,11),gc_50,gc_58,amp(43))
2688 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2689 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2690 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2691 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
2692 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
2693 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2694 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2695 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
2696 $ -1./2.*amp(38)-1./2.*amp(39)
2697
2698 matrix = 0.d0
2699 DO i = 1, ncolor
2700 ztemp = (0.d0,0.d0)
2701 DO j = 1, ncolor
2702 ztemp = ztemp + cf(j,i)*jamp(j)
2703 ENDDO
2704 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2705 ENDDO
2706 matrix_ddx_ccx_noh=matrix
2707 ENDIF
2708 ENDIF
2709 END
2710
2711C -----------------SUBPROCESS DDX->UCX, NO HIGGS
2712
2713 SUBROUTINE ddx_ucx_noh(P,H1,H2,ANS)
2714C
2715C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2716C By the MadGraph Development Team
2717C Please visit us at https://launchpad.net/madgraph5
2718C
2719C MadGraph StandAlone Version
2720C
2721C Returns amplitude squared summed/avg over colors
2722C and helicities
2723C for the point in phase space P(0:3,NEXTERNAL)
2724C
2725C Process: d d~ > u c~ ta+ ta- / h QED=4
2726C
2727 IMPLICIT NONE
2728C
2729C CONSTANTS
2730C
2731 INTEGER NEXTERNAL
2732 parameter(nexternal=6)
2733 INTEGER NCOMB
2734 parameter( ncomb=64)
2735C
2736C ARGUMENTS
2737C
2738 real*8 p(0:3,nexternal),ans
2739 INTEGER H1,H2
2740C
2741C LOCAL VARIABLES
2742C
2743 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2744 real*8 t
2745 real*8 matrix_ddx_ucx_noh
2746 INTEGER IHEL,IDEN, I
2747 INTEGER JC(NEXTERNAL)
2748 LOGICAL GOODHEL(NCOMB)
2749 DATA ntry/0/
2750 DATA goodhel/ncomb*.false./
2751 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2752 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2753 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2754 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2755 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2756 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2757 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2758 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2759 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2760 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2761 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2762 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2763 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2764 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2765 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2766 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2767 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2768 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2769 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2770 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2771 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2772 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2773 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2774 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2775 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2776 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2777 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2778 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2779 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2780 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2781 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2782 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2783 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2784 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2785 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2786 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2787 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2788 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2789 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2790 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2791 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2792 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2793 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2794 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2795 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2796 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2797 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2798 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2799 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2800 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2801 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2802 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2803 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2804 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2805 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2806 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2807 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2808 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2809 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2810 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2811 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2812 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2813 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2814 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2815 DATA iden/36/
2816C ----------
2817C BEGIN CODE
2818C ----------
2819 DO ihel=1,nexternal
2820 jc(ihel) = +1
2821 ENDDO
2822 ans = 0d0
2823 DO ihel=1,ncomb
2824 t=matrix_ddx_ucx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2825 ans=ans+t
2826 ENDDO
2827 ans=ans/dble(iden)
2828 END
2829
2830
2831 REAL*8 FUNCTION matrix_ddx_ucx_noh(P,H1,H2,NHEL,IC)
2832C
2833C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2834C By the MadGraph Development Team
2835C Please visit us at https://launchpad.net/madgraph5
2836C
2837C Returns amplitude squared summed/avg over colors
2838C for the point with external lines W(0:6,NEXTERNAL)
2839C
2840C Process: d d~ > u c~ ta+ ta- / h QED=4
2841C
2842 IMPLICIT NONE
2843C
2844C CONSTANTS
2845C
2846 INTEGER NGRAPHS
2847 parameter(ngraphs=11)
2848 INTEGER NEXTERNAL
2849 parameter(nexternal=6)
2850 INTEGER NWAVEFUNCS, NCOLOR
2851 parameter(nwavefuncs=10, ncolor=1)
2852 real*8 zero
2853 parameter(zero=0d0)
2854 COMPLEX*16 IMAG1
2855 parameter(imag1=(0d0,1d0))
2856C
2857C ARGUMENTS
2858C
2859 real*8 p(0:3,nexternal)
2860 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2861C
2862C LOCAL VARIABLES
2863C
2864 INTEGER I,J
2865 COMPLEX*16 ZTEMP
2866 real*8 denom(ncolor), cf(ncolor,ncolor)
2867 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2868 COMPLEX*16 W(18,NWAVEFUNCS)
2869 COMPLEX*16 DUM0,DUM1
2870 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2871C
2872C GLOBAL VARIABLES
2873C
2874 include 'coupl.inc'
2875C
2876C COLOR DATA
2877C
2878 DATA denom(1)/1/
2879 DATA (cf(i, 1),i= 1, 1) / 9/
2880C 1 T(2,4) T(3,1)
2881
2882 INTEGER H1,H2
2883 real*8 matrix
2884 matrix_ddx_ucx_noh=0.d0
2885 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2886 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2887
2888C ----------
2889C BEGIN CODE
2890C ----------
2891 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2892 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2893 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2894 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2895 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2896 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2897 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
2898 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,8))
2899 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2900C Amplitude(s) for diagram number 1
2901 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2902 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2903C Amplitude(s) for diagram number 2
2904 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2905 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2906C Amplitude(s) for diagram number 3
2907 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2908 CALL ffv2_1(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
2909C Amplitude(s) for diagram number 4
2910 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2911 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
2912C Amplitude(s) for diagram number 5
2913 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2914C Amplitude(s) for diagram number 6
2915 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2916C Amplitude(s) for diagram number 7
2917 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2918 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
2919C Amplitude(s) for diagram number 8
2920 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2921 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2922C Amplitude(s) for diagram number 9
2923 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
2924C Amplitude(s) for diagram number 10
2925 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2926 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2927C Amplitude(s) for diagram number 11
2928 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
2929 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2930 $ +amp(9)+amp(10)+amp(11)
2931
2932 matrix = 0.d0
2933 DO i = 1, ncolor
2934 ztemp = (0.d0,0.d0)
2935 DO j = 1, ncolor
2936 ztemp = ztemp + cf(j,i)*jamp(j)
2937 ENDDO
2938 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2939 ENDDO
2940 matrix_ddx_ucx_noh=matrix
2941 ENDIF
2942 ENDIF
2943 END
2944
2945C ----------------SUBPROCESS DDX->CUX, NO HIGGS
2946
2947 SUBROUTINE ddx_cux_noh(P,H1,H2,ANS)
2948C
2949C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
2950C By the MadGraph Development Team
2951C Please visit us at https://launchpad.net/madgraph5
2952C
2953C MadGraph StandAlone Version
2954C
2955C Returns amplitude squared summed/avg over colors
2956C and helicities
2957C for the point in phase space P(0:3,NEXTERNAL)
2958C
2959C Process: d d~ > c u~ ta+ ta- / h QED=4
2960C
2961 IMPLICIT NONE
2962C
2963C CONSTANTS
2964C
2965 INTEGER NEXTERNAL
2966 parameter(nexternal=6)
2967 INTEGER NCOMB
2968 parameter( ncomb=64)
2969C
2970C ARGUMENTS
2971C
2972 real*8 p(0:3,nexternal),ans
2973 INTEGER H1,H2
2974C
2975C LOCAL VARIABLES
2976C
2977 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2978 real*8 t
2979 real*8 matrix_ddx_cux_noh
2980 INTEGER IHEL,IDEN, I
2981 INTEGER JC(NEXTERNAL)
2982 LOGICAL GOODHEL(NCOMB)
2983 DATA ntry/0/
2984 DATA goodhel/ncomb*.false./
2985 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2986 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2987 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2988 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2989 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2990 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2991 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2992 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2993 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2994 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2995 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2996 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2997 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2998 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2999 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3000 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3001 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3002 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3003 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3004 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3005 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3006 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3007 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3008 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3009 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3010 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3011 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3012 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3013 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3014 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3015 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3016 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3017 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3018 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3019 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3020 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3021 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3022 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3023 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3024 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3025 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3026 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3027 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3028 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3029 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3030 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3031 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3032 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3033 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3034 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3035 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3036 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3037 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3038 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3039 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3040 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3041 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3042 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3043 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3044 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3045 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3046 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3047 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3048 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3049 DATA iden/36/
3050C ----------
3051C BEGIN CODE
3052C ----------
3053 DO ihel=1,nexternal
3054 jc(ihel) = +1
3055 ENDDO
3056 ans = 0d0
3057 DO ihel=1,ncomb
3058 t=matrix_ddx_cux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3059 ans=ans+t
3060 ENDDO
3061 ans=ans/dble(iden)
3062 END
3063
3064
3065 REAL*8 FUNCTION matrix_ddx_cux_noh(P,H1,H2,NHEL,IC)
3066C
3067C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
3068C By the MadGraph Development Team
3069C Please visit us at https://launchpad.net/madgraph5
3070C
3071C Returns amplitude squared summed/avg over colors
3072C for the point with external lines W(0:6,NEXTERNAL)
3073C
3074C Process: d d~ > c u~ ta+ ta- / h QED=4
3075C
3076 IMPLICIT NONE
3077C
3078C CONSTANTS
3079C
3080 INTEGER NGRAPHS
3081 parameter(ngraphs=11)
3082 INTEGER NEXTERNAL
3083 parameter(nexternal=6)
3084 INTEGER NWAVEFUNCS, NCOLOR
3085 parameter(nwavefuncs=10, ncolor=1)
3086 real*8 zero
3087 parameter(zero=0d0)
3088 COMPLEX*16 IMAG1
3089 parameter(imag1=(0d0,1d0))
3090C
3091C ARGUMENTS
3092C
3093 real*8 p(0:3,nexternal)
3094 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
3095C
3096C LOCAL VARIABLES
3097C
3098 INTEGER I,J
3099 COMPLEX*16 ZTEMP
3100 real*8 denom(ncolor), cf(ncolor,ncolor)
3101 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
3102 COMPLEX*16 W(18,NWAVEFUNCS)
3103 COMPLEX*16 DUM0,DUM1
3104 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3105C
3106C GLOBAL VARIABLES
3107C
3108 include 'coupl.inc'
3109C
3110C COLOR DATA
3111C
3112 DATA denom(1)/1/
3113 DATA (cf(i, 1),i= 1, 1) / 9/
3114C 1 T(2,4) T(3,1)
3115
3116 INTEGER H1,H2
3117 real*8 matrix
3118 matrix_ddx_cux_noh=0.d0
3119 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
3120 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
3121
3122C ----------
3123C BEGIN CODE
3124C ----------
3125 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3126 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3127 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3128 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
3129 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3130 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3131 CALL ffv2_3(w(1,1),w(1,3),gc_44,mw,ww,w(1,7))
3132 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
3133 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
3134C Amplitude(s) for diagram number 1
3135 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
3136 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
3137C Amplitude(s) for diagram number 2
3138 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
3139 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
3140C Amplitude(s) for diagram number 3
3141 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
3142 CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
3143C Amplitude(s) for diagram number 4
3144 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
3145 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
3146C Amplitude(s) for diagram number 5
3147 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
3148C Amplitude(s) for diagram number 6
3149 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
3150C Amplitude(s) for diagram number 7
3151 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
3152 CALL ffv2_2(w(1,1),w(1,8),gc_44,zero,zero,w(1,5))
3153C Amplitude(s) for diagram number 8
3154 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
3155 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
3156C Amplitude(s) for diagram number 9
3157 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_44,amp(9))
3158C Amplitude(s) for diagram number 10
3159 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
3160 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
3161C Amplitude(s) for diagram number 11
3162 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_44,amp(11))
3163 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
3164 $ +amp(9)+amp(10)+amp(11)
3165
3166 matrix = 0.d0
3167 DO i = 1, ncolor
3168 ztemp = (0.d0,0.d0)
3169 DO j = 1, ncolor
3170 ztemp = ztemp + cf(j,i)*jamp(j)
3171 ENDDO
3172 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3173 ENDDO
3174 matrix_ddx_cux_noh=matrix
3175 ENDIF
3176 ENDIF
3177 END
3178
3179C DONE 19.11.2014
3180
3181C ------------------------------------------------------
3182C WITH GLUONS AS JETS ADDED 9.07.2015
3183
3184C --- BEGIN SUBROUTINE DDX->GG NO HIGGS
3185
3186
3187 SUBROUTINE ddx_gg_noh(P,H1,H2,ANS)
3188C
3189C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
3190C By the MadGraph Development Team
3191C Please visit us at https://launchpad.net/madgraph5
3192C
3193C MadGraph StandAlone Version
3194C
3195C Returns amplitude squared summed/avg over colors
3196C and helicities
3197C for the point in phase space P(0:3,NEXTERNAL)
3198C
3199C Process: d d~ > g g ta+ ta- / h QED=4
3200C Process: s s~ > g g ta+ ta- / h QED=4
3201C
3202 IMPLICIT NONE
3203C
3204C CONSTANTS
3205C
3206 INTEGER NEXTERNAL
3207 parameter(nexternal=6)
3208 INTEGER NCOMB
3209 parameter( ncomb=64)
3210C
3211C ARGUMENTS
3212C
3213 real*8 p(0:3,nexternal),ans
3214 INTEGER H1,H2
3215C
3216C LOCAL VARIABLES
3217C
3218 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
3219 real*8 t
3220 real*8 matrix_ddx_gg_noh
3221 INTEGER IHEL,IDEN, I
3222 INTEGER JC(NEXTERNAL)
3223 LOGICAL GOODHEL(NCOMB)
3224 DATA ntry/0/
3225 DATA goodhel/ncomb*.false./
3226 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
3227 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
3228 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
3229 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
3230 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
3231 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
3232 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
3233 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
3234 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
3235 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
3236 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
3237 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
3238 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
3239 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
3240 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3241 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3242 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3243 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3244 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3245 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3246 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3247 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3248 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3249 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3250 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3251 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3252 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3253 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3254 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3255 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3256 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3257 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3258 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3259 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3260 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3261 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3262 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3263 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3264 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3265 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3266 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3267 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3268 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3269 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3270 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3271 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3272 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3273 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3274 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3275 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3276 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3277 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3278 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3279 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3280 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3281 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3282 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3283 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3284 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3285 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3286 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3287 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3288 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3289 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3290 DATA iden/72/
3291C ----------
3292C BEGIN CODE
3293C ----------
3294 DO ihel=1,nexternal
3295 jc(ihel) = +1
3296 ENDDO
3297 ans = 0d0
3298 DO ihel=1,ncomb
3299 t=matrix_ddx_gg_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3300 ans=ans+t
3301 ENDDO
3302 ans=ans/dble(iden)
3303 END
3304
3305
3306 REAL*8 FUNCTION matrix_ddx_gg_noh(P,H1,H2,NHEL,IC)
3307C
3308C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
3309C By the MadGraph Development Team
3310C Please visit us at https://launchpad.net/madgraph5
3311C
3312C Returns amplitude squared summed/avg over colors
3313C for the point with external lines W(0:6,NEXTERNAL)
3314C
3315C Process: d d~ > g g ta+ ta- / h QED=4
3316C Process: s s~ > g g ta+ ta- / h QED=4
3317C
3318 IMPLICIT NONE
3319C
3320C CONSTANTS
3321C
3322 INTEGER NGRAPHS
3323 parameter(ngraphs=16)
3324 INTEGER NEXTERNAL
3325 parameter(nexternal=6)
3326 INTEGER NWAVEFUNCS, NCOLOR
3327 parameter(nwavefuncs=10, ncolor=2)
3328 real*8 zero
3329 parameter(zero=0d0)
3330 COMPLEX*16 IMAG1
3331 parameter(imag1=(0d0,1d0))
3332C
3333C ARGUMENTS
3334C
3335 real*8 p(0:3,nexternal)
3336 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
3337C
3338C LOCAL VARIABLES
3339C
3340 INTEGER I,J
3341 COMPLEX*16 ZTEMP
3342 real*8 denom(ncolor), cf(ncolor,ncolor)
3343 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
3344 COMPLEX*16 W(18,NWAVEFUNCS)
3345 COMPLEX*16 DUM0,DUM1
3346 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3347C
3348C GLOBAL VARIABLES
3349C
3350 include 'coupl.inc'
3351C
3352C COLOR DATA
3353C
3354 DATA denom(1)/3/
3355 DATA (cf(i, 1),i= 1, 2) / 16, -2/
3356C 1 T(3,4,2,1)
3357 DATA denom(2)/3/
3358 DATA (cf(i, 2),i= 1, 2) / -2, 16/
3359C 1 T(4,3,2,1)
3360
3361 INTEGER H1,H2
3362 real*8 matrix
3363 matrix_ddx_gg_noh=0.d0
3364 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
3365 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
3366
3367C ----------
3368C BEGIN CODE
3369C ----------
3370 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3371 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3372 CALL vxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3373 CALL vxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
3374 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3375 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3376 CALL vvv1p0_1(w(1,3),w(1,4),gc_10,zero,zero,w(1,7))
3377 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,8))
3378 CALL ffv1_1(w(1,2),w(1,7),gc_11,zero,zero,w(1,9))
3379C Amplitude(s) for diagram number 1
3380 CALL ffv1_0(w(1,1),w(1,9),w(1,8),gc_1,amp(1))
3381 CALL ffv1_2(w(1,1),w(1,7),gc_11,zero,zero,w(1,10))
3382C Amplitude(s) for diagram number 2
3383 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_1,amp(2))
3384 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
3385C Amplitude(s) for diagram number 3
3386 CALL ffv2_3_0(w(1,1),w(1,9),w(1,7),gc_50,gc_58,amp(3))
3387C Amplitude(s) for diagram number 4
3388 CALL ffv2_3_0(w(1,10),w(1,2),w(1,7),gc_50,gc_58,amp(4))
3389 CALL ffv1_1(w(1,2),w(1,3),gc_11,zero,zero,w(1,10))
3390 CALL ffv1_2(w(1,1),w(1,4),gc_11,zero,zero,w(1,9))
3391C Amplitude(s) for diagram number 5
3392 CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_1,amp(5))
3393C Amplitude(s) for diagram number 6
3394 CALL ffv2_3_0(w(1,9),w(1,10),w(1,7),gc_50,gc_58,amp(6))
3395 CALL ffv1_1(w(1,10),w(1,4),gc_11,zero,zero,w(1,6))
3396C Amplitude(s) for diagram number 7
3397 CALL ffv1_0(w(1,1),w(1,6),w(1,8),gc_1,amp(7))
3398C Amplitude(s) for diagram number 8
3399 CALL ffv2_3_0(w(1,1),w(1,6),w(1,7),gc_50,gc_58,amp(8))
3400 CALL ffv1_2(w(1,1),w(1,3),gc_11,zero,zero,w(1,6))
3401 CALL ffv1_1(w(1,2),w(1,4),gc_11,zero,zero,w(1,10))
3402C Amplitude(s) for diagram number 9
3403 CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_1,amp(9))
3404C Amplitude(s) for diagram number 10
3405 CALL ffv2_3_0(w(1,6),w(1,10),w(1,7),gc_50,gc_58,amp(10))
3406 CALL ffv1_2(w(1,6),w(1,4),gc_11,zero,zero,w(1,5))
3407C Amplitude(s) for diagram number 11
3408 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(11))
3409C Amplitude(s) for diagram number 12
3410 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(12))
3411 CALL ffv1_1(w(1,10),w(1,3),gc_11,zero,zero,w(1,5))
3412C Amplitude(s) for diagram number 13
3413 CALL ffv1_0(w(1,1),w(1,5),w(1,8),gc_1,amp(13))
3414C Amplitude(s) for diagram number 14
3415 CALL ffv2_3_0(w(1,1),w(1,5),w(1,7),gc_50,gc_58,amp(14))
3416 CALL ffv1_2(w(1,9),w(1,3),gc_11,zero,zero,w(1,5))
3417C Amplitude(s) for diagram number 15
3418 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_1,amp(15))
3419C Amplitude(s) for diagram number 16
3420 CALL ffv2_3_0(w(1,5),w(1,2),w(1,7),gc_50,gc_58,amp(16))
3421 jamp(1)=-imag1*amp(1)-imag1*amp(2)-imag1*amp(3)-imag1*amp(4)
3422 $ +amp(5)+amp(6)+amp(7)+amp(8)+amp(15)+amp(16)
3423 jamp(2)=+imag1*amp(1)+imag1*amp(2)+imag1*amp(3)+imag1*amp(4)
3424 $ +amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
3425
3426 matrix = 0.d0
3427 DO i = 1, ncolor
3428 ztemp = (0.d0,0.d0)
3429 DO j = 1, ncolor
3430 ztemp = ztemp + cf(j,i)*jamp(j)
3431 ENDDO
3432 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3433 ENDDO
3434 matrix_ddx_gg_noh=matrix
3435 ENDIF
3436 ENDIF
3437 END
3438
3439