C++ Interface to Tauola
SDX.f
1
2c --------------------- subprocesses initiated by s dbar: SDX 21.10.2015
3
4 SUBROUTINE sdx(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 sdx_uux_h(p,h1,h2,ans)
16 IF(i3.EQ.2 .AND. i4.EQ.-4) CALL sdx_ucx_h(p,h1,h2,ans)
17 IF(i3.EQ.3 .AND. i4.EQ.-1) CAll sdx_sdx_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.-2) CALL sdx_cux_h(p,h1,h2,ans)
19 IF(i3.EQ.4 .AND. i4.EQ.-4) CALL sdx_ccx_h(p,h1,h2,ans)
20 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
21 CALL sdx_uux_h(p,h1,h2,ans1)
22 CALL sdx_ucx_h(p,h1,h2,ans2)
23 CALL sdx_sdx_h(p,h1,h2,ans3)
24 CALL sdx_cux_h(p,h1,h2,ans4)
25 CALL sdx_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 sdx_uux_noh(p,h1,h2,ans)
30 IF(i3.EQ.2 .AND. i4.EQ.-4) CALL sdx_ucx_noh(p,h1,h2,ans)
31 IF(i3.EQ.3 .AND. i4.EQ.-1) CAll sdx_sdx_noh(p,h1,h2,ans)
32 IF(i3.EQ.4 .AND. i4.EQ.-2) CALL sdx_cux_noh(p,h1,h2,ans)
33 IF(i3.EQ.4 .AND. i4.EQ.-4) CALL sdx_ccx_noh(p,h1,h2,ans)
34 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
35 CALL sdx_uux_noh(p,h1,h2,ans1)
36 CALL sdx_ucx_noh(p,h1,h2,ans2)
37 CALL sdx_sdx_noh(p,h1,h2,ans3)
38 CALL sdx_cux_noh(p,h1,h2,ans4)
39 CALL sdx_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 SDX
48
49
50
51C ----- begin subprocesses SDX->jjH with Higgs->tautau
52
53 SUBROUTINE sdx_uux_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: s d~ > u u~ h WEIGHTED=6
66C * Decay: h > ta+ ta- WEIGHTED=2
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_sdx_uux_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_sdx_uux_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_sdx_uux_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: s d~ > u u~ h WEIGHTED=6
182C * Decay: h > ta+ ta- WEIGHTED=2
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_sdx_uux_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_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,6))
241 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
242C Amplitude(s) for diagram number 1
243 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,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_sdx_uux_h=matrix
255 ENDIF ! CLOSES H1 IF
256 ENDIF ! CLOSES H2 IF
257 END
258
259c==========================================
260
261 SUBROUTINE sdx_ucx_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: s d~ > u c~ h WEIGHTED=6
274C * Decay: h > ta+ ta- WEIGHTED=2
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_sdx_ucx_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_sdx_ucx_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_sdx_ucx_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: s d~ > u c~ h WEIGHTED=6
390C * Decay: h > ta+ ta- WEIGHTED=2
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_sdx_ucx_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
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_101,mw,ww,w(1,6))
450 CALL ffv2_3(w(1,4),w(1,2),gc_44,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_sdx_ucx_h=matrix
464 ENDIF
465 ENDIF
466 END
467
468
469
470C===============================================
471
472 SUBROUTINE sdx_sdx_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: s d~ > s d~ 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_sdx_sdx_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_sdx_sdx_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_sdx_sdx_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: s d~ > s d~ 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_sdx_sdx_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 matrix_sdx_sdx_h=matrix
674 ENDIF
675 ENDIF
676 END
677
678
679
680C========================================
681
682 SUBROUTINE sdx_cux_h(P,H1,H2,ANS)
683C
684C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
685C By the MadGraph Development Team
686C Please visit us at https://launchpad.net/madgraph5
687C
688C MadGraph StandAlone Version
689C
690C Returns amplitude squared summed/avg over colors
691C and helicities
692C for the point in phase space P(0:3,NEXTERNAL)
693C
694C Process: s d~ > c u~ h WEIGHTED=6
695C * Decay: h > ta+ ta- WEIGHTED=2
696C
697 IMPLICIT NONE
698C
699C CONSTANTS
700C
701 INTEGER NEXTERNAL
702 parameter(nexternal=6)
703 INTEGER NCOMB
704 parameter( ncomb=64)
705C
706C ARGUMENTS
707C
708 real*8 p(0:3,nexternal),ans
709 INTEGER H1,H2
710C
711C LOCAL VARIABLES
712C
713 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
714 real*8 t
715 real*8 matrix_sdx_cux_h
716 INTEGER IHEL,IDEN, I
717 INTEGER JC(NEXTERNAL)
718 LOGICAL GOODHEL(NCOMB)
719 DATA ntry/0/
720 DATA goodhel/ncomb*.false./
721 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
722 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
723 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
724 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
725 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
726 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
727 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
728 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
729 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
730 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
731 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
732 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
733 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
734 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
735 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
736 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
737 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
738 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
739 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
740 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
741 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
742 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
743 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
744 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
745 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
746 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
747 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
748 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
749 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
750 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
751 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
752 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
753 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
754 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
755 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
756 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
757 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
758 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
759 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
760 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
761 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
762 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
763 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
764 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
765 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
766 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
767 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
768 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
769 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
770 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
771 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
772 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
773 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
774 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
775 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
776 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
777 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
778 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
779 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
780 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
781 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
782 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
783 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
784 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
785 DATA iden/36/
786C ----------
787C BEGIN CODE
788C ----------
789 DO ihel=1,nexternal
790 jc(ihel) = +1
791 ENDDO
792 ans = 0d0
793 DO ihel=1,ncomb
794 t=matrix_sdx_cux_h(p ,h1,h2,nhel(1,ihel),jc(1))
795 ans=ans+t
796 ENDDO
797 ans=ans/dble(iden)
798 END
799
800
801 REAL*8 FUNCTION matrix_sdx_cux_h(P,H1,H2,NHEL,IC)
802C
803C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
804C By the MadGraph Development Team
805C Please visit us at https://launchpad.net/madgraph5
806C
807C Returns amplitude squared summed/avg over colors
808C for the point with external lines W(0:6,NEXTERNAL)
809C
810C Process: s d~ > c u~ h WEIGHTED=6
811C * Decay: h > ta+ ta- WEIGHTED=2
812C
813 IMPLICIT NONE
814C
815C CONSTANTS
816C
817 INTEGER NGRAPHS
818 parameter(ngraphs=1)
819 INTEGER NEXTERNAL
820 parameter(nexternal=6)
821 INTEGER NWAVEFUNCS, NCOLOR
822 parameter(nwavefuncs=7, ncolor=1)
823 real*8 zero
824 parameter(zero=0d0)
825 COMPLEX*16 IMAG1
826 parameter(imag1=(0d0,1d0))
827C
828C ARGUMENTS
829C
830 real*8 p(0:3,nexternal)
831 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
832C
833C LOCAL VARIABLES
834C
835 INTEGER I,J
836 COMPLEX*16 ZTEMP
837 real*8 denom(ncolor), cf(ncolor,ncolor)
838 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
839 COMPLEX*16 W(18,NWAVEFUNCS)
840 COMPLEX*16 DUM0,DUM1
841 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
842C
843C GLOBAL VARIABLES
844C
845 include 'coupl.inc'
846C
847C COLOR DATA
848C
849 DATA denom(1)/1/
850 DATA (cf(i, 1),i= 1, 1) / 9/
851C 1 T(2,4) T(3,1)
852
853 INTEGER H1,H2
854 real*8 matrix
855 matrix_sdx_cux_h=0.d0
856 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
857 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
858
859C ----------
860C BEGIN CODE
861C ----------
862 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
863 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
864 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
865 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
866 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
867 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
868 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
869 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
870 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
871C Amplitude(s) for diagram number 1
872 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
873 jamp(1)=+amp(1)
874
875 matrix = 0.d0
876 DO i = 1, ncolor
877 ztemp = (0.d0,0.d0)
878 DO j = 1, ncolor
879 ztemp = ztemp + cf(j,i)*jamp(j)
880 ENDDO
881 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
882 ENDDO
883 matrix_sdx_cux_h=matrix
884 ENDIF
885 ENDIF
886 END
887
888
889 SUBROUTINE sdx_ccx_h(P,H1,H2,ANS)
890C
891C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
892C By the MadGraph Development Team
893C Please visit us at https://launchpad.net/madgraph5
894C
895C MadGraph StandAlone Version
896C
897C Returns amplitude squared summed/avg over colors
898C and helicities
899C for the point in phase space P(0:3,NEXTERNAL)
900C
901C Process: s d~ > c c~ h WEIGHTED=6
902C * Decay: h > ta+ ta- WEIGHTED=2
903C
904 IMPLICIT NONE
905C
906C CONSTANTS
907C
908 INTEGER NEXTERNAL
909 parameter(nexternal=6)
910 INTEGER NCOMB
911 parameter( ncomb=64)
912C
913C ARGUMENTS
914C
915 real*8 p(0:3,nexternal),ans
916 INTEGER H1,H2
917C
918C LOCAL VARIABLES
919C
920 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
921 real*8 t
922 real*8 matrix_sdx_ccx_h
923 INTEGER IHEL,IDEN, I
924 INTEGER JC(NEXTERNAL)
925 LOGICAL GOODHEL(NCOMB)
926 DATA ntry/0/
927 DATA goodhel/ncomb*.false./
928 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
929 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
930 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
931 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
932 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
933 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
934 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
935 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
936 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
937 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
938 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
939 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
940 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
941 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
942 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
943 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
944 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
945 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
946 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
947 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
948 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
949 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
950 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
951 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
952 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
953 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
954 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
955 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
956 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
957 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
958 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
959 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
960 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
961 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
962 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
963 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
964 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
965 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
966 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
967 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
968 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
969 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
970 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
971 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
972 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
973 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
974 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
975 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
976 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
977 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
978 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
979 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
980 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
981 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
982 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
983 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
984 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
985 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
986 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
987 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
988 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
989 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
990 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
991 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
992 DATA iden/36/
993C ----------
994C BEGIN CODE
995C ----------
996 DO ihel=1,nexternal
997 jc(ihel) = +1
998 ENDDO
999 ans = 0d0
1000 DO ihel=1,ncomb
1001 t=matrix_sdx_ccx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1002 ans=ans+t
1003 ENDDO
1004 ans=ans/dble(iden)
1005 END
1006
1007
1008 REAL*8 FUNCTION matrix_sdx_ccx_h(P,H1,H2,NHEL,IC)
1009C
1010C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1011C By the MadGraph Development Team
1012C Please visit us at https://launchpad.net/madgraph5
1013C
1014C Returns amplitude squared summed/avg over colors
1015C for the point with external lines W(0:6,NEXTERNAL)
1016C
1017C Process: s d~ > c c~ h WEIGHTED=6
1018C * Decay: h > ta+ ta- WEIGHTED=2
1019C
1020 IMPLICIT NONE
1021C
1022C CONSTANTS
1023C
1024 INTEGER NGRAPHS
1025 parameter(ngraphs=1)
1026 INTEGER NEXTERNAL
1027 parameter(nexternal=6)
1028 INTEGER NWAVEFUNCS, NCOLOR
1029 parameter(nwavefuncs=7, ncolor=1)
1030 real*8 zero
1031 parameter(zero=0d0)
1032 COMPLEX*16 IMAG1
1033 parameter(imag1=(0d0,1d0))
1034C
1035C ARGUMENTS
1036C
1037 real*8 p(0:3,nexternal)
1038 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1039C
1040C LOCAL VARIABLES
1041C
1042 INTEGER I,J
1043 COMPLEX*16 ZTEMP
1044 real*8 denom(ncolor), cf(ncolor,ncolor)
1045 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1046 COMPLEX*16 W(18,NWAVEFUNCS)
1047 COMPLEX*16 DUM0,DUM1
1048 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1049C
1050C GLOBAL VARIABLES
1051C
1052 include 'coupl.inc'
1053C
1054C COLOR DATA
1055C
1056 DATA denom(1)/1/
1057 DATA (cf(i, 1),i= 1, 1) / 9/
1058C 1 T(2,4) T(3,1)
1059
1060 INTEGER H1,H2
1061 real*8 matrix
1062 matrix_sdx_ccx_h=0.d0
1063 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1064 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1065
1066C ----------
1067C BEGIN CODE
1068C ----------
1069 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1070 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1071 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1072 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1073 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1074 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1075 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1076 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
1077 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,3))
1078C Amplitude(s) for diagram number 1
1079 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
1080 jamp(1)=+amp(1)
1081
1082 matrix = 0.d0
1083 DO i = 1, ncolor
1084 ztemp = (0.d0,0.d0)
1085 DO j = 1, ncolor
1086 ztemp = ztemp + cf(j,i)*jamp(j)
1087 ENDDO
1088 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1089 ENDDO
1090
1091 matrix_sdx_ccx_h=matrix
1092 ENDIF
1093 ENDIF
1094 END
1095
1096C --- end of subprocess SDX->jjH, H->tautau
1097
1098C ----- begin subprocesses SDX->jjtautau, no Higgs
1099
1100 SUBROUTINE sdx_uux_noh(P,H1,H2,ANS)
1101C
1102C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1103C By the MadGraph Development Team
1104C Please visit us at https://launchpad.net/madgraph5
1105C
1106C MadGraph StandAlone Version
1107C
1108C Returns amplitude squared summed/avg over colors
1109C and helicities
1110C for the point in phase space P(0:3,NEXTERNAL)
1111C
1112C Process: s d~ > u u~ ta+ ta- / h Qed=4
1113C
1114 IMPLICIT NONE
1115C
1116C CONSTANTS
1117C
1118 INTEGER NEXTERNAL
1119 parameter(nexternal=6)
1120 INTEGER NCOMB
1121 parameter( ncomb=64)
1122C
1123C ARGUMENTS
1124C
1125 real*8 p(0:3,nexternal),ans
1126 INTEGER H1,H2
1127C
1128C LOCAL VARIABLES
1129C
1130 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1131 real*8 t
1132 real*8 matrix_sdx_uux_noh
1133 INTEGER IHEL,IDEN, I
1134 INTEGER JC(NEXTERNAL)
1135 LOGICAL GOODHEL(NCOMB)
1136 DATA ntry/0/
1137 DATA goodhel/ncomb*.false./
1138 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1139 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1140 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1141 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1142 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1143 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1144 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1145 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1146 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1147 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1148 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1149 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1150 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1151 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1152 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1153 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1154 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1155 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1156 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1157 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1158 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1159 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1160 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1161 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1162 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1163 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1164 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1165 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1166 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1167 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1168 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1169 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1170 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1171 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1172 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1173 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1174 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1175 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1176 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1177 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1178 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1179 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1180 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1181 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1182 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1183 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1184 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1185 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1186 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1187 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1188 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1189 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1190 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1191 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1192 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1193 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1194 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1195 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1196 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1197 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1198 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1199 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1200 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1201 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1202 DATA iden/36/
1203C ----------
1204C BEGIN CODE
1205C ----------
1206 DO ihel=1,nexternal
1207 jc(ihel) = +1
1208 ENDDO
1209 ans = 0d0
1210 DO ihel=1,ncomb
1211 t=matrix_sdx_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1212 ans=ans+t
1213 ENDDO
1214 ans=ans/dble(iden)
1215 END
1216
1217
1218 REAL*8 FUNCTION matrix_sdx_uux_noh(P,H1,H2,NHEL,IC)
1219C
1220C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1221C By the MadGraph Development Team
1222C Please visit us at https://launchpad.net/madgraph5
1223C
1224C Returns amplitude squared summed/avg over colors
1225C for the point with external lines W(0:6,NEXTERNAL)
1226C
1227C Process: s d~ > u u~ ta+ ta- / h Qed=4
1228C
1229 IMPLICIT NONE
1230C
1231C CONSTANTS
1232C
1233 INTEGER NGRAPHS
1234 parameter(ngraphs=11)
1235 INTEGER NEXTERNAL
1236 parameter(nexternal=6)
1237 INTEGER NWAVEFUNCS, NCOLOR
1238 parameter(nwavefuncs=10, ncolor=1)
1239 real*8 zero
1240 parameter(zero=0d0)
1241 COMPLEX*16 IMAG1
1242 parameter(imag1=(0d0,1d0))
1243C
1244C ARGUMENTS
1245C
1246 real*8 p(0:3,nexternal)
1247 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1248C
1249C LOCAL VARIABLES
1250C
1251 INTEGER I,J
1252 COMPLEX*16 ZTEMP
1253 real*8 denom(ncolor), cf(ncolor,ncolor)
1254 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1255 COMPLEX*16 W(18,NWAVEFUNCS)
1256 COMPLEX*16 DUM0,DUM1
1257 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1258C
1259C GLOBAL VARIABLES
1260C
1261 include 'coupl.inc'
1262C
1263C COLOR DATA
1264C
1265 DATA denom(1)/1/
1266 DATA (cf(i, 1),i= 1, 1) / 9/
1267C 1 T(2,4) T(3,1)
1268
1269 INTEGER H1,H2
1270 real*8 matrix
1271 matrix_sdx_uux_noh=0.d0
1272 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1273 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1274
1275C ----------
1276C BEGIN CODE
1277C ----------
1278 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1279 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1280 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1281 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1282 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1283 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1284 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,7))
1285 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
1286 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1287C Amplitude(s) for diagram number 1
1288 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1289 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1290C Amplitude(s) for diagram number 2
1291 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1292 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1293C Amplitude(s) for diagram number 3
1294 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1295 CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1296C Amplitude(s) for diagram number 4
1297 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1298 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1299C Amplitude(s) for diagram number 5
1300 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
1301C Amplitude(s) for diagram number 6
1302 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1303C Amplitude(s) for diagram number 7
1304 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
1305 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1306C Amplitude(s) for diagram number 8
1307 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1308 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1309C Amplitude(s) for diagram number 9
1310 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_101,amp(9))
1311C Amplitude(s) for diagram number 10
1312 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1313 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1314C Amplitude(s) for diagram number 11
1315 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_101,amp(11))
1316 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1317 $ +amp(9)+amp(10)+amp(11)
1318
1319 matrix = 0.d0
1320 DO i = 1, ncolor
1321 ztemp = (0.d0,0.d0)
1322 DO j = 1, ncolor
1323 ztemp = ztemp + cf(j,i)*jamp(j)
1324 ENDDO
1325 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1326 ENDDO
1327 matrix_sdx_uux_noh=matrix
1328 ENDIF ! CLOSES H1 IF
1329 ENDIF ! CLOSES H2 IF
1330 END
1331
1332
1333 SUBROUTINE sdx_ucx_noh(P,H1,H2,ANS)
1334C
1335C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1336C By the MadGraph Development Team
1337C Please visit us at https://launchpad.net/madgraph5
1338C
1339C MadGraph StandAlone Version
1340C
1341C Returns amplitude squared summed/avg over colors
1342C and helicities
1343C for the point in phase space P(0:3,NEXTERNAL)
1344C
1345C Process: s d~ > u c~ ta+ ta- / h Qed=4
1346C
1347 IMPLICIT NONE
1348C
1349C CONSTANTS
1350C
1351 INTEGER NEXTERNAL
1352 parameter(nexternal=6)
1353 INTEGER NCOMB
1354 parameter( ncomb=64)
1355C
1356C ARGUMENTS
1357C
1358 real*8 p(0:3,nexternal),ans
1359 INTEGER H1,H2
1360C
1361C LOCAL VARIABLES
1362C
1363 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1364 real*8 t
1365 real*8 matrix_sdx_ucx_noh
1366 INTEGER IHEL,IDEN, I
1367 INTEGER JC(NEXTERNAL)
1368 LOGICAL GOODHEL(NCOMB)
1369 DATA ntry/0/
1370 DATA goodhel/ncomb*.false./
1371 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1372 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1373 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1374 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1375 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1376 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1377 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1378 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1379 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1380 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1381 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1382 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1383 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1384 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1385 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1386 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1387 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1388 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1389 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1390 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1391 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1392 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1393 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1394 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1395 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1396 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1397 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1398 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1399 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1400 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1401 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1402 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1403 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1404 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1405 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1406 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1407 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1408 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1409 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1410 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1411 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1412 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1413 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1414 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1415 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1416 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1417 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1418 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1419 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1420 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1421 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1422 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1423 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1424 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1425 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1426 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1427 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1428 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1429 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1430 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1431 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1432 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1433 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1434 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1435 DATA iden/36/
1436C ----------
1437C BEGIN CODE
1438C ----------
1439 DO ihel=1,nexternal
1440 jc(ihel) = +1
1441 ENDDO
1442 ans = 0d0
1443 DO ihel=1,ncomb
1444 t=matrix_sdx_ucx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1445 ans=ans+t
1446 ENDDO
1447 ans=ans/dble(iden)
1448 END
1449
1450
1451 REAL*8 FUNCTION matrix_sdx_ucx_noh(P,H1,H2,NHEL,IC)
1452C
1453C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1454C By the MadGraph Development Team
1455C Please visit us at https://launchpad.net/madgraph5
1456C
1457C Returns amplitude squared summed/avg over colors
1458C for the point with external lines W(0:6,NEXTERNAL)
1459C
1460C Process: s d~ > u c~ ta+ ta- / h Qed=4
1461C
1462 IMPLICIT NONE
1463C
1464C CONSTANTS
1465C
1466 INTEGER NGRAPHS
1467 parameter(ngraphs=11)
1468 INTEGER NEXTERNAL
1469 parameter(nexternal=6)
1470 INTEGER NWAVEFUNCS, NCOLOR
1471 parameter(nwavefuncs=10, ncolor=1)
1472 real*8 zero
1473 parameter(zero=0d0)
1474 COMPLEX*16 IMAG1
1475 parameter(imag1=(0d0,1d0))
1476C
1477C ARGUMENTS
1478C
1479 real*8 p(0:3,nexternal)
1480 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1481C
1482C LOCAL VARIABLES
1483C
1484 INTEGER I,J
1485 COMPLEX*16 ZTEMP
1486 real*8 denom(ncolor), cf(ncolor,ncolor)
1487 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1488 COMPLEX*16 W(18,NWAVEFUNCS)
1489 COMPLEX*16 DUM0,DUM1
1490 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1491C
1492C GLOBAL VARIABLES
1493C
1494 include 'coupl.inc'
1495C
1496C COLOR DATA
1497C
1498 DATA denom(1)/1/
1499 DATA (cf(i, 1),i= 1, 1) / 9/
1500C 1 T(2,4) T(3,1)
1501
1502 INTEGER H1,H2
1503 real*8 matrix
1504 matrix_sdx_ucx_noh=0.d0
1505 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1506 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1507
1508C ----------
1509C BEGIN CODE
1510C ----------
1511 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1512 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1513 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1514 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1515 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1516 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1517 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,7))
1518 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,8))
1519 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1520C Amplitude(s) for diagram number 1
1521 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1522 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1523C Amplitude(s) for diagram number 2
1524 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1525 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1526C Amplitude(s) for diagram number 3
1527 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1528 CALL ffv2_1(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
1529C Amplitude(s) for diagram number 4
1530 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1531 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
1532C Amplitude(s) for diagram number 5
1533 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
1534C Amplitude(s) for diagram number 6
1535 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1536C Amplitude(s) for diagram number 7
1537 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
1538 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1539C Amplitude(s) for diagram number 8
1540 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1541 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1542C Amplitude(s) for diagram number 9
1543 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_101,amp(9))
1544C Amplitude(s) for diagram number 10
1545 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1546 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1547C Amplitude(s) for diagram number 11
1548 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_101,amp(11))
1549 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1550 $ +amp(9)+amp(10)+amp(11)
1551
1552 matrix = 0.d0
1553 DO i = 1, ncolor
1554 ztemp = (0.d0,0.d0)
1555 DO j = 1, ncolor
1556 ztemp = ztemp + cf(j,i)*jamp(j)
1557 ENDDO
1558 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1559 ENDDO
1560 matrix_sdx_ucx_noh=matrix
1561 ENDIF ! CLOSES H1 IF
1562 ENDIF ! CLOSES H2 IF
1563 END
1564
1565 SUBROUTINE sdx_sdx_noh(P,H1,H2,ANS)
1566C
1567C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1568C By the MadGraph Development Team
1569C Please visit us at https://launchpad.net/madgraph5
1570C
1571C MadGraph StandAlone Version
1572C
1573C Returns amplitude squared summed/avg over colors
1574C and helicities
1575C for the point in phase space P(0:3,NEXTERNAL)
1576C
1577C Process: s d~ > s d~ ta+ ta- / h Qed=4
1578C
1579 IMPLICIT NONE
1580C
1581C CONSTANTS
1582C
1583 INTEGER NEXTERNAL
1584 parameter(nexternal=6)
1585 INTEGER NCOMB
1586 parameter( ncomb=64)
1587C
1588C ARGUMENTS
1589C
1590 real*8 p(0:3,nexternal),ans
1591 INTEGER H1,H2
1592C
1593C LOCAL VARIABLES
1594C
1595 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1596 real*8 t
1597 real*8 matrix_sdx_sdx_noh
1598 INTEGER IHEL,IDEN, I
1599 INTEGER JC(NEXTERNAL)
1600 LOGICAL GOODHEL(NCOMB)
1601 DATA ntry/0/
1602 DATA goodhel/ncomb*.false./
1603 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1604 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1605 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1606 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1607 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1608 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1609 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1610 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1611 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1612 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1613 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1614 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1615 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1616 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1617 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1618 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1619 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1620 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1621 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1622 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1623 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1624 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1625 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1626 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1627 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1628 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1629 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1630 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1631 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1632 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1633 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1634 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1635 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1636 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1637 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1638 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1639 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1640 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1641 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1642 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1643 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1644 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1645 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1646 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1647 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1648 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1649 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1650 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1651 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1652 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1653 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1654 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1655 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1656 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1657 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1658 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1659 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1660 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1661 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1662 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1663 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1664 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1665 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1666 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1667 DATA iden/36/
1668C ----------
1669C BEGIN CODE
1670C ----------
1671 DO ihel=1,nexternal
1672 jc(ihel) = +1
1673 ENDDO
1674 ans = 0d0
1675 DO ihel=1,ncomb
1676 t=matrix_sdx_sdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1677 ans=ans+t
1678 ENDDO
1679 ans=ans/dble(iden)
1680 END
1681
1682
1683 REAL*8 FUNCTION matrix_sdx_sdx_noh(P,H1,H2,NHEL,IC)
1684C
1685C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1686C By the MadGraph Development Team
1687C Please visit us at https://launchpad.net/madgraph5
1688C
1689C Returns amplitude squared summed/avg over colors
1690C for the point with external lines W(0:6,NEXTERNAL)
1691C
1692C Process: s d~ > s d~ ta+ ta- / h Qed=4
1693C
1694 IMPLICIT NONE
1695C
1696C CONSTANTS
1697C
1698 INTEGER NGRAPHS
1699 parameter(ngraphs=32)
1700 INTEGER NEXTERNAL
1701 parameter(nexternal=6)
1702 INTEGER NWAVEFUNCS, NCOLOR
1703 parameter(nwavefuncs=13, ncolor=2)
1704 real*8 zero
1705 parameter(zero=0d0)
1706 COMPLEX*16 IMAG1
1707 parameter(imag1=(0d0,1d0))
1708C
1709C ARGUMENTS
1710C
1711 real*8 p(0:3,nexternal)
1712 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1713C
1714C LOCAL VARIABLES
1715C
1716 INTEGER I,J
1717 COMPLEX*16 ZTEMP
1718 real*8 denom(ncolor), cf(ncolor,ncolor)
1719 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1720 COMPLEX*16 W(18,NWAVEFUNCS)
1721 COMPLEX*16 DUM0,DUM1
1722 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1723C
1724C GLOBAL VARIABLES
1725C
1726 include 'coupl.inc'
1727C
1728C COLOR DATA
1729C
1730 DATA denom(1)/1/
1731 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1732C 1 T(2,1) T(3,4)
1733 DATA denom(2)/1/
1734 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1735C 1 T(2,4) T(3,1)
1736 INTEGER H1,H2
1737 real*8 matrix
1738 matrix_sdx_sdx_noh=0.d0
1739 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1740 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1741
1742C ----------
1743C BEGIN CODE
1744C ----------
1745 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1746 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1747 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1748 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1749 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1750 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1751 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,7))
1752 CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,8))
1753 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1754C Amplitude(s) for diagram number 1
1755 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1756 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1757C Amplitude(s) for diagram number 2
1758 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1759 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,11))
1760C Amplitude(s) for diagram number 3
1761 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1762C Amplitude(s) for diagram number 4
1763 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1764 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,10))
1765 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1766C Amplitude(s) for diagram number 5
1767 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1768 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1769C Amplitude(s) for diagram number 6
1770 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1771C Amplitude(s) for diagram number 7
1772 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1773C Amplitude(s) for diagram number 8
1774 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1775 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1776 CALL ffv1_1(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
1777C Amplitude(s) for diagram number 9
1778 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
1779 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1780C Amplitude(s) for diagram number 10
1781 CALL ffv1_0(w(1,13),w(1,2),w(1,12),gc_1,amp(10))
1782 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1783C Amplitude(s) for diagram number 11
1784 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
1785C Amplitude(s) for diagram number 12
1786 CALL ffv2_3_0(w(1,13),w(1,2),w(1,7),gc_50,gc_58,amp(12))
1787 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
1788 CALL ffv1_1(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
1789C Amplitude(s) for diagram number 13
1790 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
1791 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
1792C Amplitude(s) for diagram number 14
1793 CALL ffv1_0(w(1,6),w(1,2),w(1,12),gc_1,amp(14))
1794C Amplitude(s) for diagram number 15
1795 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
1796C Amplitude(s) for diagram number 16
1797 CALL ffv2_3_0(w(1,6),w(1,2),w(1,7),gc_50,gc_58,amp(16))
1798 CALL ffv2_3_1(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
1799C Amplitude(s) for diagram number 17
1800 CALL ffv1_0(w(1,4),w(1,6),w(1,12),gc_1,amp(17))
1801 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1802C Amplitude(s) for diagram number 18
1803 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_1,amp(18))
1804C Amplitude(s) for diagram number 19
1805 CALL ffv2_3_0(w(1,4),w(1,6),w(1,7),gc_50,gc_58,amp(19))
1806C Amplitude(s) for diagram number 20
1807 CALL ffv2_3_0(w(1,9),w(1,2),w(1,7),gc_50,gc_58,amp(20))
1808 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
1809C Amplitude(s) for diagram number 21
1810 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(21))
1811 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
1812C Amplitude(s) for diagram number 22
1813 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_1,amp(22))
1814C Amplitude(s) for diagram number 23
1815 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(23))
1816 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
1817C Amplitude(s) for diagram number 24
1818 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_1,amp(24))
1819 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,8))
1820 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1821C Amplitude(s) for diagram number 25
1822 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_1,amp(25))
1823C Amplitude(s) for diagram number 26
1824 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(26))
1825C Amplitude(s) for diagram number 27
1826 CALL ffv2_3_0(w(1,4),w(1,3),w(1,7),gc_50,gc_58,amp(27))
1827C Amplitude(s) for diagram number 28
1828 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(28))
1829 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
1830C Amplitude(s) for diagram number 29
1831 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_1,amp(29))
1832C Amplitude(s) for diagram number 30
1833 CALL ffv2_3_0(w(1,6),w(1,3),w(1,11),gc_50,gc_58,amp(30))
1834C Amplitude(s) for diagram number 31
1835 CALL ffv2_3_0(w(1,8),w(1,3),w(1,7),gc_50,gc_58,amp(31))
1836C Amplitude(s) for diagram number 32
1837 CALL ffv2_3_0(w(1,9),w(1,3),w(1,11),gc_50,gc_58,amp(32))
1838 jamp(1)=+1./2.*(+amp(13)+amp(14)+amp(15)+amp(16)+amp(25)+amp(26)
1839 $ +amp(27)+amp(28))
1840 jamp(2)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
1841 $ +amp(9)+amp(10)+amp(11)+amp(12)-1./6.*amp(13)-1./6.*amp(14)
1842 $ -1./6.*amp(15)-1./6.*amp(16)+amp(17)+amp(18)+amp(19)+amp(20)
1843 $ +amp(21)+amp(22)+amp(23)+amp(24)-1./6.*amp(25)-1./6.*amp(26)
1844 $ -1./6.*amp(27)-1./6.*amp(28)+amp(29)+amp(30)+amp(31)+amp(32)
1845
1846 matrix = 0.d0
1847 DO i = 1, ncolor
1848 ztemp = (0.d0,0.d0)
1849 DO j = 1, ncolor
1850 ztemp = ztemp + cf(j,i)*jamp(j)
1851 ENDDO
1852 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1853 ENDDO
1854 matrix_sdx_sdx_noh=matrix
1855 ENDIF ! CLOSES H1 IF
1856 ENDIF ! CLOSES H2 IF
1857 END
1858
1859
1860 SUBROUTINE sdx_cux_noh(P,H1,H2,ANS)
1861C
1862C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1863C By the MadGraph Development Team
1864C Please visit us at https://launchpad.net/madgraph5
1865C
1866C MadGraph StandAlone Version
1867C
1868C Returns amplitude squared summed/avg over colors
1869C and helicities
1870C for the point in phase space P(0:3,NEXTERNAL)
1871C
1872C Process: s d~ > c u~ ta+ ta- / h Qed=4
1873C
1874 IMPLICIT NONE
1875C
1876C CONSTANTS
1877C
1878 INTEGER NEXTERNAL
1879 parameter(nexternal=6)
1880 INTEGER NCOMB
1881 parameter( ncomb=64)
1882C
1883C ARGUMENTS
1884C
1885 real*8 p(0:3,nexternal),ans
1886 INTEGER H1,H2
1887C
1888C LOCAL VARIABLES
1889C
1890 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1891 real*8 t
1892 real*8 matrix_sdx_cux_noh
1893 INTEGER IHEL,IDEN, I
1894 INTEGER JC(NEXTERNAL)
1895 LOGICAL GOODHEL(NCOMB)
1896 DATA ntry/0/
1897 DATA goodhel/ncomb*.false./
1898 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1899 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1900 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1901 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1902 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1903 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1904 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1905 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1906 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1907 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1908 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1909 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1910 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1911 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1912 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1913 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1914 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1915 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1916 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1917 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1918 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1919 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1920 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1921 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1922 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1923 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1924 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1925 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1926 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1927 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1928 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1929 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1930 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1931 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1932 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1933 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1934 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1935 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1936 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1937 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1938 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1939 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1940 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1941 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1942 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1943 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1944 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1945 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1946 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1947 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1948 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1949 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1950 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1951 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1952 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1953 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1954 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1955 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1956 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1957 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1958 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1959 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1960 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1961 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1962 DATA iden/36/
1963C ----------
1964C BEGIN CODE
1965C ----------
1966 DO ihel=1,nexternal
1967 jc(ihel) = +1
1968 ENDDO
1969 ans = 0d0
1970 DO ihel=1,ncomb
1971 t=matrix_sdx_cux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1972 ans=ans+t
1973 ENDDO
1974 ans=ans/dble(iden)
1975 END
1976
1977
1978 REAL*8 FUNCTION matrix_sdx_cux_noh(P,H1,H2,NHEL,IC)
1979C
1980C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1981C By the MadGraph Development Team
1982C Please visit us at https://launchpad.net/madgraph5
1983C
1984C Returns amplitude squared summed/avg over colors
1985C for the point with external lines W(0:6,NEXTERNAL)
1986C
1987C Process: s d~ > c u~ ta+ ta- / h Qed=4
1988C
1989 IMPLICIT NONE
1990C
1991C CONSTANTS
1992C
1993 INTEGER NGRAPHS
1994 parameter(ngraphs=11)
1995 INTEGER NEXTERNAL
1996 parameter(nexternal=6)
1997 INTEGER NWAVEFUNCS, NCOLOR
1998 parameter(nwavefuncs=10, ncolor=1)
1999 real*8 zero
2000 parameter(zero=0d0)
2001 COMPLEX*16 IMAG1
2002 parameter(imag1=(0d0,1d0))
2003C
2004C ARGUMENTS
2005C
2006 real*8 p(0:3,nexternal)
2007 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2008C
2009C LOCAL VARIABLES
2010C
2011 INTEGER I,J
2012 COMPLEX*16 ZTEMP
2013 real*8 denom(ncolor), cf(ncolor,ncolor)
2014 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2015 COMPLEX*16 W(18,NWAVEFUNCS)
2016 COMPLEX*16 DUM0,DUM1
2017 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2018C
2019C GLOBAL VARIABLES
2020C
2021 include 'coupl.inc'
2022C
2023C COLOR DATA
2024C
2025 DATA denom(1)/1/
2026 DATA (cf(i, 1),i= 1, 1) / 9/
2027C 1 T(2,4) T(3,1)
2028
2029 INTEGER H1,H2
2030 real*8 matrix
2031 matrix_sdx_cux_noh=0.d0
2032 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2033 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2034
2035
2036C ----------
2037C BEGIN CODE
2038C ----------
2039 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2040 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2041 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2042 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2043 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2044 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2045 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
2046 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
2047 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2048C Amplitude(s) for diagram number 1
2049 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2050 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2051C Amplitude(s) for diagram number 2
2052 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2053 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2054C Amplitude(s) for diagram number 3
2055 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2056 CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
2057C Amplitude(s) for diagram number 4
2058 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2059 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
2060C Amplitude(s) for diagram number 5
2061 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2062C Amplitude(s) for diagram number 6
2063 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2064C Amplitude(s) for diagram number 7
2065 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2066 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
2067C Amplitude(s) for diagram number 8
2068 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2069 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2070C Amplitude(s) for diagram number 9
2071 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
2072C Amplitude(s) for diagram number 10
2073 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2074 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2075C Amplitude(s) for diagram number 11
2076 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
2077 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2078 $ +amp(9)+amp(10)+amp(11)
2079
2080 matrix = 0.d0
2081 DO i = 1, ncolor
2082 ztemp = (0.d0,0.d0)
2083 DO j = 1, ncolor
2084 ztemp = ztemp + cf(j,i)*jamp(j)
2085 ENDDO
2086 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2087 ENDDO
2088 matrix_sdx_cux_noh=matrix
2089 ENDIF ! CLOSES H1 IF
2090 ENDIF ! CLOSES H2 IF
2091 END
2092
2093
2094
2095 SUBROUTINE sdx_ccx_noh(P,H1,H2,ANS)
2096C
2097C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
2098C By the MadGraph Development Team
2099C Please visit us at https://launchpad.net/madgraph5
2100C
2101C MadGraph StandAlone Version
2102C
2103C Returns amplitude squared summed/avg over colors
2104C and helicities
2105C for the point in phase space P(0:3,NEXTERNAL)
2106C
2107C Process: s d~ > c c~ ta+ ta- / h Qed=4
2108C
2109 IMPLICIT NONE
2110C
2111C CONSTANTS
2112C
2113 INTEGER NEXTERNAL
2114 parameter(nexternal=6)
2115 INTEGER NCOMB
2116 parameter( ncomb=64)
2117C
2118C ARGUMENTS
2119C
2120 real*8 p(0:3,nexternal),ans
2121 INTEGER H1,H2
2122C
2123C LOCAL VARIABLES
2124C
2125 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2126 real*8 t
2127 real*8 matrix_sdx_ccx_noh
2128 INTEGER IHEL,IDEN, I
2129 INTEGER JC(NEXTERNAL)
2130 LOGICAL GOODHEL(NCOMB)
2131 DATA ntry/0/
2132 DATA goodhel/ncomb*.false./
2133 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2134 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2135 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2136 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2137 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2138 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2139 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2140 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2141 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2142 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2143 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2144 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2145 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2146 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2147 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2148 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2149 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2150 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2151 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2152 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2153 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2154 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2155 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2156 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2157 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2158 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2159 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2160 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2161 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2162 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2163 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2164 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2165 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2166 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2167 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2168 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2169 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2170 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2171 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2172 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2173 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2174 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2175 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2176 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2177 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2178 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2179 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2180 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2181 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2182 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2183 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2184 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2185 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2186 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2187 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2188 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2189 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2190 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2191 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2192 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2193 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2194 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2195 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2196 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2197 DATA iden/36/
2198C ----------
2199C BEGIN CODE
2200C ----------
2201 DO ihel=1,nexternal
2202 jc(ihel) = +1
2203 ENDDO
2204 ans = 0d0
2205 DO ihel=1,ncomb
2206 t=matrix_sdx_ccx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2207 ans=ans+t
2208 ENDDO
2209 ans=ans/dble(iden)
2210 END
2211
2212
2213 REAL*8 FUNCTION matrix_sdx_ccx_noh(P,H1,H2,NHEL,IC)
2214C
2215C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
2216C By the MadGraph Development Team
2217C Please visit us at https://launchpad.net/madgraph5
2218C
2219C Returns amplitude squared summed/avg over colors
2220C for the point with external lines W(0:6,NEXTERNAL)
2221C
2222C Process: s d~ > c c~ ta+ ta- / h Qed=4
2223C
2224 IMPLICIT NONE
2225C
2226C CONSTANTS
2227C
2228 INTEGER NGRAPHS
2229 parameter(ngraphs=11)
2230 INTEGER NEXTERNAL
2231 parameter(nexternal=6)
2232 INTEGER NWAVEFUNCS, NCOLOR
2233 parameter(nwavefuncs=10, ncolor=1)
2234 real*8 zero
2235 parameter(zero=0d0)
2236 COMPLEX*16 IMAG1
2237 parameter(imag1=(0d0,1d0))
2238C
2239C ARGUMENTS
2240C
2241 real*8 p(0:3,nexternal)
2242 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2243C
2244C LOCAL VARIABLES
2245C
2246 INTEGER I,J
2247 COMPLEX*16 ZTEMP
2248 real*8 denom(ncolor), cf(ncolor,ncolor)
2249 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2250 COMPLEX*16 W(18,NWAVEFUNCS)
2251 COMPLEX*16 DUM0,DUM1
2252 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2253C
2254C GLOBAL VARIABLES
2255C
2256 include 'coupl.inc'
2257C
2258C COLOR DATA
2259C
2260 DATA denom(1)/1/
2261 DATA (cf(i, 1),i= 1, 1) / 9/
2262C 1 T(2,4) T(3,1)
2263
2264 INTEGER H1,H2
2265 real*8 matrix
2266 matrix_sdx_ccx_noh=0.d0
2267 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
2268 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
2269
2270C ----------
2271C BEGIN CODE
2272C ----------
2273 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2274 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2275 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2276 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2277 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2278 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2279 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
2280 CALL ffv2_3(w(1,4),w(1,2),gc_44,mw,ww,w(1,8))
2281 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
2282C Amplitude(s) for diagram number 1
2283 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
2284 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2285C Amplitude(s) for diagram number 2
2286 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
2287 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2288C Amplitude(s) for diagram number 3
2289 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
2290 CALL ffv2_1(w(1,2),w(1,7),gc_44,zero,zero,w(1,6))
2291C Amplitude(s) for diagram number 4
2292 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
2293 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
2294C Amplitude(s) for diagram number 5
2295 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(5))
2296C Amplitude(s) for diagram number 6
2297 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
2298C Amplitude(s) for diagram number 7
2299 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(7))
2300 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
2301C Amplitude(s) for diagram number 8
2302 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
2303 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
2304C Amplitude(s) for diagram number 9
2305 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
2306C Amplitude(s) for diagram number 10
2307 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
2308 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
2309C Amplitude(s) for diagram number 11
2310 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
2311 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2312 $ +amp(9)+amp(10)+amp(11)
2313
2314 matrix = 0.d0
2315 DO i = 1, ncolor
2316 ztemp = (0.d0,0.d0)
2317 DO j = 1, ncolor
2318 ztemp = ztemp + cf(j,i)*jamp(j)
2319 ENDDO
2320 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2321 ENDDO
2322 matrix_sdx_ccx_noh=matrix
2323 ENDIF ! CLOSES H1 IF
2324 ENDIF ! CLOSES H2 IF
2325 END
2326
2327C ------- end subprocesses SDX->jjtautau, no Higgs
2328
2329