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