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