C++ Interface to Tauola
SUX.f
1c --------------------- subprocesses initiated by SUbar
2
3 SUBROUTINE sux(P,I3,I4,H1,H2,KEY,ANS)
4 IMPLICIT NONE
5
6 INTEGER I3,I4, H1,H2, KEY
7 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
8
9C I3=0,1 OR 3 I4=0,-2 OR -4
10
11 ans=0.d0
12
13C FINAL SCX IS EQUAL TO FINAL DUX
14
15 IF(key.EQ.1) THEN
16 IF(i3.EQ.1 .AND. i4.EQ.-2) CALL sux_dux_h(p,h1,h2,ans)
17 IF(i3.EQ.3 .AND. i4.EQ.-2) CALL sux_sux_h(p,h1,h2,ans)
18 IF(i3.EQ.1 .AND. i4.EQ.-4) CAll sux_dcx_h(p,h1,h2,ans)
19 IF(i3.EQ.3 .AND. i4.EQ.-4) CALL sux_dux_h(p,h1,h2,ans) !FINAL SCX EQUAL TO DUX
20 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
21 CALL sux_dux_h(p,h1,h2,ans1)
22 CALL sux_sux_h(p,h1,h2,ans2)
23 CALL sux_dcx_h(p,h1,h2,ans3)
24 CALL sux_dux_h(p,h1,h2,ans4)
25 ans=ans1+ans2+ans3+ans4
26 ENDIF
27 ELSE IF(key.EQ.0) THEN
28 IF(i3.EQ.1 .AND. i4.EQ.-2) CALL sux_dux_noh(p,h1,h2,ans)
29 IF(i3.EQ.3 .AND. i4.EQ.-2) CALL sux_sux_noh(p,h1,h2,ans)
30 IF(i3.EQ.1 .AND. i4.EQ.-4) CAll sux_dcx_noh(p,h1,h2,ans)
31 IF(i3.EQ.3 .AND. i4.EQ.-4) CALL sux_dux_noh(p,h1,h2,ans)
32 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
33 CALL sux_dux_noh(p,h1,h2,ans1)
34 CALL sux_sux_noh(p,h1,h2,ans2)
35 CALL sux_dcx_noh(p,h1,h2,ans3)
36 CALL sux_dux_noh(p,h1,h2,ans4)
37 ans=ans1+ans2+ans3+ans4
38 ENDIF
39 ELSE
40 WRITE(*,*) 'NOT FINISHED'
41 stop
42 ENDIF
43 END ! SUBROUTINE SUX
44
45
46
47C ----- begin subprocesses SUX->DUX with Higgs->tautau
48
49
50 SUBROUTINE sux_dux_h(P,H1,H2,ANS)
51C
52C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
53C By the MadGraph Development Team
54C Please visit us at https://launchpad.net/madgraph5
55C
56C MadGraph StandAlone Version
57C
58C Returns amplitude squared summed/avg over colors
59C and helicities
60C for the point in phase space P(0:3,NEXTERNAL)
61C
62C Process: s u~ > d u~ h WEIGHTED=6
63C * Decay: h > ta+ ta- WEIGHTED=2
64C Process: s u~ > s c~ h WEIGHTED=6
65C * Decay: h > ta+ ta- WEIGHTED=2
66C
67 IMPLICIT NONE
68C
69C CONSTANTS
70C
71 INTEGER NEXTERNAL
72 parameter(nexternal=6)
73 INTEGER NCOMB
74 parameter( ncomb=64)
75C
76C ARGUMENTS
77C
78 real*8 p(0:3,nexternal),ans
79 INTEGER H1,H2
80C
81C LOCAL VARIABLES
82C
83 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
84 real*8 t
85 real*8 matrix_sux_dux_h
86 INTEGER IHEL,IDEN, I
87 INTEGER JC(NEXTERNAL)
88 LOGICAL GOODHEL(NCOMB)
89 DATA ntry/0/
90 DATA goodhel/ncomb*.false./
91 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
92 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
93 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
94 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
95 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
96 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
97 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
98 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
99 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
100 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
101 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
102 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
103 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
104 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
105 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
106 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
107 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
108 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
109 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
110 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
111 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
112 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
113 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
114 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
115 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
116 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
117 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
118 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
119 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
120 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
121 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
122 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
123 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
124 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
125 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
126 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
127 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
128 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
129 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
130 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
131 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
132 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
133 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
134 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
135 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
136 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
137 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
138 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
139 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
140 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
141 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
142 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
143 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
144 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
145 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
146 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
147 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
148 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
149 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
150 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
151 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
152 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
153 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
154 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
155 DATA iden/36/
156C ----------
157C BEGIN CODE
158C ----------
159 DO ihel=1,nexternal
160 jc(ihel) = +1
161 ENDDO
162 ans = 0d0
163 DO ihel=1,ncomb
164 t=matrix_sux_dux_h(p ,h1,h2,nhel(1,ihel),jc(1))
165 ans=ans+t
166 ENDDO
167 ans=ans/dble(iden)
168 END
169
170
171 REAL*8 FUNCTION matrix_sux_dux_h(P,H1,H2,NHEL,IC)
172C
173C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
174C By the MadGraph Development Team
175C Please visit us at https://launchpad.net/madgraph5
176C
177C Returns amplitude squared summed/avg over colors
178C for the point with external lines W(0:6,NEXTERNAL)
179C
180C Process: s u~ > d u~ h WEIGHTED=6
181C * Decay: h > ta+ ta- WEIGHTED=2
182C Process: s u~ > s c~ h WEIGHTED=6
183C * Decay: h > ta+ ta- WEIGHTED=2
184C
185 IMPLICIT NONE
186C
187C CONSTANTS
188C
189 INTEGER NGRAPHS
190 parameter(ngraphs=1)
191 INTEGER NEXTERNAL
192 parameter(nexternal=6)
193 INTEGER NWAVEFUNCS, NCOLOR
194 parameter(nwavefuncs=7, ncolor=1)
195 real*8 zero
196 parameter(zero=0d0)
197 COMPLEX*16 IMAG1
198 parameter(imag1=(0d0,1d0))
199C
200C ARGUMENTS
201C
202 real*8 p(0:3,nexternal)
203 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
204C
205C LOCAL VARIABLES
206C
207 INTEGER I,J
208 COMPLEX*16 ZTEMP
209 real*8 denom(ncolor), cf(ncolor,ncolor)
210 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
211 COMPLEX*16 W(18,NWAVEFUNCS)
212 COMPLEX*16 DUM0,DUM1
213 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
214C
215C GLOBAL VARIABLES
216C
217 include 'coupl.inc'
218C
219C COLOR DATA
220C
221 DATA denom(1)/1/
222 DATA (cf(i, 1),i= 1, 1) / 9/
223C 1 T(2,1) T(3,4)
224
225 INTEGER H1,H2
226 real*8 matrix
227 matrix_sux_dux_h=0.d0
228 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
229 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
230
231C ----------
232C BEGIN CODE
233C ----------
234 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
235 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
236 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
237 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
238 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
239 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
240 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
241 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
242 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
243C Amplitude(s) for diagram number 1
244 CALL vvs1_0(w(1,2),w(1,6),w(1,7),gc_72,amp(1))
245 jamp(1)=-amp(1)
246
247 matrix = 0.d0
248 DO i = 1, ncolor
249 ztemp = (0.d0,0.d0)
250 DO j = 1, ncolor
251 ztemp = ztemp + cf(j,i)*jamp(j)
252 ENDDO
253 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
254 ENDDO
255
256 matrix_sux_dux_h=matrix
257 ENDIF
258 ENDIF
259 END
260
261C ----- begin subprocesses SUX->SUX with Higgs->tautau
262
263 SUBROUTINE sux_sux_h(P,H1,H2,ANS)
264C
265C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
266C By the MadGraph Development Team
267C Please visit us at https://launchpad.net/madgraph5
268C
269C MadGraph StandAlone Version
270C
271C Returns amplitude squared summed/avg over colors
272C and helicities
273C for the point in phase space P(0:3,NEXTERNAL)
274C
275C Process: s u~ > s u~ h WEIGHTED=6
276C * Decay: h > ta+ ta- WEIGHTED=2
277C
278 IMPLICIT NONE
279C
280C CONSTANTS
281C
282 INTEGER NEXTERNAL
283 parameter(nexternal=6)
284 INTEGER NCOMB
285 parameter( ncomb=64)
286C
287C ARGUMENTS
288C
289 real*8 p(0:3,nexternal),ans
290 INTEGER H1,H2
291C
292C LOCAL VARIABLES
293C
294 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
295 real*8 t
296 real*8 matrix_sux_sux_h
297 INTEGER IHEL,IDEN, I
298 INTEGER JC(NEXTERNAL)
299 LOGICAL GOODHEL(NCOMB)
300 DATA ntry/0/
301 DATA goodhel/ncomb*.false./
302 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
303 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
304 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
305 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
306 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
307 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
308 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
309 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
310 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
311 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
312 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
313 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
314 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
315 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
316 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
317 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
318 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
319 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
320 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
321 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
322 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
323 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
324 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
325 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
326 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
327 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
328 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
329 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
330 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
331 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
332 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
333 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
334 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
335 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
336 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
337 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
338 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
339 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
340 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
341 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
342 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
343 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
344 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
345 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
346 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
347 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
348 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
349 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
350 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
351 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
352 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
353 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
354 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
355 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
356 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
357 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
358 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
359 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
360 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
361 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
362 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
363 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
364 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
365 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
366 DATA iden/36/
367C ----------
368C BEGIN CODE
369C ----------
370 DO ihel=1,nexternal
371 jc(ihel) = +1
372 ENDDO
373 ans = 0d0
374 DO ihel=1,ncomb
375 t=matrix_sux_sux_h(p ,h1,h2,nhel(1,ihel),jc(1))
376 ans=ans+t
377 ENDDO
378 ans=ans/dble(iden)
379 END
380
381
382 REAL*8 FUNCTION matrix_sux_sux_h(P,H1,H2,NHEL,IC)
383C
384C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
385C By the MadGraph Development Team
386C Please visit us at https://launchpad.net/madgraph5
387C
388C Returns amplitude squared summed/avg over colors
389C for the point with external lines W(0:6,NEXTERNAL)
390C
391C Process: s u~ > s u~ h WEIGHTED=6
392C * Decay: h > ta+ ta- WEIGHTED=2
393C
394 IMPLICIT NONE
395C
396C CONSTANTS
397C
398 INTEGER NGRAPHS
399 parameter(ngraphs=2)
400 INTEGER NEXTERNAL
401 parameter(nexternal=6)
402 INTEGER NWAVEFUNCS, NCOLOR
403 parameter(nwavefuncs=7, ncolor=2)
404 real*8 zero
405 parameter(zero=0d0)
406 COMPLEX*16 IMAG1
407 parameter(imag1=(0d0,1d0))
408C
409C ARGUMENTS
410C
411 real*8 p(0:3,nexternal)
412 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
413C
414C LOCAL VARIABLES
415C
416 INTEGER I,J
417 COMPLEX*16 ZTEMP
418 real*8 denom(ncolor), cf(ncolor,ncolor)
419 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
420 COMPLEX*16 W(18,NWAVEFUNCS)
421 COMPLEX*16 DUM0,DUM1
422 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
423C
424C GLOBAL VARIABLES
425C
426 include 'coupl.inc'
427C
428C COLOR DATA
429C
430 DATA denom(1)/1/
431 DATA (cf(i, 1),i= 1, 2) / 9, 3/
432C 1 T(2,1) T(3,4)
433 DATA denom(2)/1/
434 DATA (cf(i, 2),i= 1, 2) / 3, 9/
435C 1 T(2,4) T(3,1)
436
437 INTEGER H1,H2
438 real*8 matrix
439 matrix_sux_sux_h=0.d0
440 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
441 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
442
443C ----------
444C BEGIN CODE
445C ----------
446 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
447 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
448 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
449 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
450 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
451 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
452 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
453 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
454 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,5))
455C Amplitude(s) for diagram number 1
456 CALL vvs1_0(w(1,5),w(1,6),w(1,7),gc_72,amp(1))
457 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
458 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,3))
459C Amplitude(s) for diagram number 2
460 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
461 jamp(1)=-amp(1)
462 jamp(2)=+amp(2)
463
464 matrix = 0.d0
465 DO i = 1, ncolor
466 ztemp = (0.d0,0.d0)
467 DO j = 1, ncolor
468 ztemp = ztemp + cf(j,i)*jamp(j)
469 ENDDO
470 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
471 ENDDO
472 matrix_sux_sux_h=matrix
473 ENDIF
474 ENDIF
475 END
476
477C ----- begin subprocesses SUX->DCX with Higgs->tautau
478
479 SUBROUTINE sux_dcx_h(P,H1,H2,ANS)
480C
481C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
482C By the MadGraph Development Team
483C Please visit us at https://launchpad.net/madgraph5
484C
485C MadGraph StandAlone Version
486C
487C Returns amplitude squared summed/avg over colors
488C and helicities
489C for the point in phase space P(0:3,NEXTERNAL)
490C
491C Process: s u~ > d c~ h WEIGHTED=6
492C * Decay: h > ta+ ta- WEIGHTED=2
493C
494 IMPLICIT NONE
495C
496C CONSTANTS
497C
498 INTEGER NEXTERNAL
499 parameter(nexternal=6)
500 INTEGER NCOMB
501 parameter( ncomb=64)
502C
503C ARGUMENTS
504C
505 real*8 p(0:3,nexternal),ans
506 INTEGER H1,H2
507C
508C LOCAL VARIABLES
509C
510 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
511 real*8 t
512 real*8 matrix_sux_dcx_h
513 INTEGER IHEL,IDEN, I
514 INTEGER JC(NEXTERNAL)
515 LOGICAL GOODHEL(NCOMB)
516 DATA ntry/0/
517 DATA goodhel/ncomb*.false./
518 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
519 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
520 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
521 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
522 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
523 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
524 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
525 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
526 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
527 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
528 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
529 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
530 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
531 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
532 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
533 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
534 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
535 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
536 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
537 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
538 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
539 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
540 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
541 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
542 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
543 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
544 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
545 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
546 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
547 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
548 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
549 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
550 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
551 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
552 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
553 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
554 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
555 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
556 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
557 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
558 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
559 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
560 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
561 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
562 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
563 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
564 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
565 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
566 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
567 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
568 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
569 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
570 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
571 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
572 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
573 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
574 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
575 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
576 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
577 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
578 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
579 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
580 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
581 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
582 DATA iden/36/
583C ----------
584C BEGIN CODE
585C ----------
586 DO ihel=1,nexternal
587 jc(ihel) = +1
588 ENDDO
589 ans = 0d0
590 DO ihel=1,ncomb
591 t=matrix_sux_dcx_h(p ,h1,h2,nhel(1,ihel),jc(1))
592 ans=ans+t
593 ENDDO
594 ans=ans/dble(iden)
595 END
596
597
598 REAL*8 FUNCTION matrix_sux_dcx_h(P,H1,H2,NHEL,IC)
599C
600C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
601C By the MadGraph Development Team
602C Please visit us at https://launchpad.net/madgraph5
603C
604C Returns amplitude squared summed/avg over colors
605C for the point with external lines W(0:6,NEXTERNAL)
606C
607C Process: s u~ > d c~ h WEIGHTED=6
608C * Decay: h > ta+ ta- WEIGHTED=2
609C
610 IMPLICIT NONE
611C
612C CONSTANTS
613C
614 INTEGER NGRAPHS
615 parameter(ngraphs=1)
616 INTEGER NEXTERNAL
617 parameter(nexternal=6)
618 INTEGER NWAVEFUNCS, NCOLOR
619 parameter(nwavefuncs=7, ncolor=1)
620 real*8 zero
621 parameter(zero=0d0)
622 COMPLEX*16 IMAG1
623 parameter(imag1=(0d0,1d0))
624C
625C ARGUMENTS
626C
627 real*8 p(0:3,nexternal)
628 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
629C
630C LOCAL VARIABLES
631C
632 INTEGER I,J
633 COMPLEX*16 ZTEMP
634 real*8 denom(ncolor), cf(ncolor,ncolor)
635 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
636 COMPLEX*16 W(18,NWAVEFUNCS)
637 COMPLEX*16 DUM0,DUM1
638 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
639C
640C GLOBAL VARIABLES
641C
642 include 'coupl.inc'
643C
644C COLOR DATA
645C
646 DATA denom(1)/1/
647 DATA (cf(i, 1),i= 1, 1) / 9/
648C 1 T(2,1) T(3,4)
649
650 INTEGER H1,H2
651 real*8 matrix
652 matrix_sux_dcx_h=0.d0
653 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
654 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
655
656C ----------
657C BEGIN CODE
658C ----------
659 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
660 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
661 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
662 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
663 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
664 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
665 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
666 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
667 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,2))
668C Amplitude(s) for diagram number 1
669 CALL vvs1_0(w(1,2),w(1,6),w(1,7),gc_72,amp(1))
670 jamp(1)=-amp(1)
671
672 matrix = 0.d0
673 DO i = 1, ncolor
674 ztemp = (0.d0,0.d0)
675 DO j = 1, ncolor
676 ztemp = ztemp + cf(j,i)*jamp(j)
677 ENDDO
678 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
679 ENDDO
680 matrix_sux_dcx_h=matrix
681 ENDIF
682 ENDIF
683 END
684
685C ----- begin subprocesses SUX->SCX with Higgs->tautau
686C ----- THE SAME AS SUX->DCX WITH Higgs->tautau
687
688C -------------- NOW PROCESSIS WITHOUT HIGGS
689
690C ----- begin subprocesses SUX->DUX NO Higgs
691
692 SUBROUTINE sux_dux_noh(P,H1,H2,ANS)
693C
694C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
695C By the MadGraph Development Team
696C Please visit us at https://launchpad.net/madgraph5
697C
698C MadGraph StandAlone Version
699C
700C Returns amplitude squared summed/avg over colors
701C and helicities
702C for the point in phase space P(0:3,NEXTERNAL)
703C
704C Process: s u~ > d u~ ta+ ta- / h QED=4
705C Process: s u~ > s c~ ta+ ta- / h QED=4
706C
707 IMPLICIT NONE
708C
709C CONSTANTS
710C
711 INTEGER NEXTERNAL
712 parameter(nexternal=6)
713 INTEGER NCOMB
714 parameter( ncomb=64)
715C
716C ARGUMENTS
717C
718 real*8 p(0:3,nexternal),ans
719 INTEGER H1,H2
720C
721C LOCAL VARIABLES
722C
723 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
724 real*8 t
725 real*8 matrix_sux_dux_noh
726 INTEGER IHEL,IDEN, I
727 INTEGER JC(NEXTERNAL)
728 LOGICAL GOODHEL(NCOMB)
729 DATA ntry/0/
730 DATA goodhel/ncomb*.false./
731 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
732 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
733 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
734 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
735 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
736 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
737 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
738 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
739 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
740 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
741 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
742 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
743 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
744 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
745 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
746 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
747 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
748 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
749 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
750 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
751 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
752 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
753 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
754 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
755 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
756 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
757 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
758 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
759 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
760 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
761 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
762 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
763 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
764 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
765 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
766 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
767 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
768 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
769 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
770 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
771 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
772 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
773 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
774 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
775 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
776 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
777 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
778 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
779 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
780 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
781 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
782 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
783 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
784 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
785 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
786 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
787 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
788 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
789 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
790 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
791 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
792 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
793 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
794 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
795 DATA iden/36/
796C ----------
797C BEGIN CODE
798C ----------
799 DO ihel=1,nexternal
800 jc(ihel) = +1
801 ENDDO
802 ans = 0d0
803 DO ihel=1,ncomb
804 t=matrix_sux_dux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
805 ans=ans+t
806 ENDDO
807 ans=ans/dble(iden)
808 END
809
810
811 REAL*8 FUNCTION matrix_sux_dux_noh(P,H1,H2,NHEL,IC)
812C
813C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
814C By the MadGraph Development Team
815C Please visit us at https://launchpad.net/madgraph5
816C
817C Returns amplitude squared summed/avg over colors
818C for the point with external lines W(0:6,NEXTERNAL)
819C
820C Process: s u~ > d u~ ta+ ta- / h QED=4
821C Process: s u~ > s c~ ta+ ta- / h QED=4
822C
823 IMPLICIT NONE
824C
825C CONSTANTS
826C
827 INTEGER NGRAPHS
828 parameter(ngraphs=11)
829 INTEGER NEXTERNAL
830 parameter(nexternal=6)
831 INTEGER NWAVEFUNCS, NCOLOR
832 parameter(nwavefuncs=10, ncolor=1)
833 real*8 zero
834 parameter(zero=0d0)
835 COMPLEX*16 IMAG1
836 parameter(imag1=(0d0,1d0))
837C
838C ARGUMENTS
839C
840 real*8 p(0:3,nexternal)
841 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
842C
843C LOCAL VARIABLES
844C
845 INTEGER I,J
846 COMPLEX*16 ZTEMP
847 real*8 denom(ncolor), cf(ncolor,ncolor)
848 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
849 COMPLEX*16 W(18,NWAVEFUNCS)
850 COMPLEX*16 DUM0,DUM1
851 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
852C
853C GLOBAL VARIABLES
854C
855 include 'coupl.inc'
856C
857C COLOR DATA
858C
859 DATA denom(1)/1/
860 DATA (cf(i, 1),i= 1, 1) / 9/
861C 1 T(2,1) T(3,4)
862
863 INTEGER H1,H2
864 real*8 matrix
865 matrix_sux_dux_noh=0.d0
866 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
867 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
868
869C ----------
870C BEGIN CODE
871C ----------
872 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
873 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
874 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
875 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
876 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
877 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
878 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
879 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
880 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
881C Amplitude(s) for diagram number 1
882 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
883 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
884C Amplitude(s) for diagram number 2
885 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
886 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
887C Amplitude(s) for diagram number 3
888 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
889 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
890C Amplitude(s) for diagram number 4
891 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
892 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
893C Amplitude(s) for diagram number 5
894 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
895C Amplitude(s) for diagram number 6
896 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
897C Amplitude(s) for diagram number 7
898 CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(7))
899 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
900C Amplitude(s) for diagram number 8
901 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(8))
902 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,3))
903C Amplitude(s) for diagram number 9
904 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
905C Amplitude(s) for diagram number 10
906 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
907 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
908C Amplitude(s) for diagram number 11
909 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
910 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
911 $ -amp(9)-amp(10)-amp(11)
912
913 matrix = 0.d0
914 DO i = 1, ncolor
915 ztemp = (0.d0,0.d0)
916 DO j = 1, ncolor
917 ztemp = ztemp + cf(j,i)*jamp(j)
918 ENDDO
919 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
920 ENDDO
921 matrix_sux_dux_noh=matrix
922 ENDIF
923 ENDIF
924 END
925
926C ----- begin subprocesses SUX->SUX NO Higgs
927
928 SUBROUTINE sux_sux_noh(P,H1,H2,ANS)
929C
930C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
931C By the MadGraph Development Team
932C Please visit us at https://launchpad.net/madgraph5
933C
934C MadGraph StandAlone Version
935C
936C Returns amplitude squared summed/avg over colors
937C and helicities
938C for the point in phase space P(0:3,NEXTERNAL)
939C
940C Process: s u~ > s u~ ta+ ta- / h QED=4
941C
942 IMPLICIT NONE
943C
944C CONSTANTS
945C
946 INTEGER NEXTERNAL
947 parameter(nexternal=6)
948 INTEGER NCOMB
949 parameter( ncomb=64)
950C
951C ARGUMENTS
952C
953 real*8 p(0:3,nexternal),ans
954 INTEGER H1,H2
955C
956C LOCAL VARIABLES
957C
958 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
959 real*8 t
960 real*8 matrix_sux_sux_noh
961 INTEGER IHEL,IDEN, I
962 INTEGER JC(NEXTERNAL)
963 LOGICAL GOODHEL(NCOMB)
964 DATA ntry/0/
965 DATA goodhel/ncomb*.false./
966 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
967 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
968 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
969 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
970 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
971 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
972 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
973 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
974 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
975 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
976 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
977 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
978 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
979 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
980 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
981 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
982 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
983 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
984 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
985 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
986 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
987 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
988 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
989 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
990 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
991 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
992 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
993 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
994 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
995 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
996 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
997 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
998 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
999 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1000 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1001 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1002 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1003 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1004 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1005 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1006 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1007 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1008 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1009 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1010 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1011 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1012 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1013 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1014 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1015 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1016 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1017 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1018 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1019 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1020 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1021 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1022 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1023 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1024 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1025 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1026 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1027 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1028 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1029 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1030 DATA iden/36/
1031C ----------
1032C BEGIN CODE
1033C ----------
1034 DO ihel=1,nexternal
1035 jc(ihel) = +1
1036 ENDDO
1037 ans = 0d0
1038 DO ihel=1,ncomb
1039 t=matrix_sux_sux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1040 ans=ans+t
1041 ENDDO
1042 ans=ans/dble(iden)
1043 END
1044
1045
1046 REAL*8 FUNCTION matrix_sux_sux_noh(P,H1,H2,NHEL,IC)
1047C
1048C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1049C By the MadGraph Development Team
1050C Please visit us at https://launchpad.net/madgraph5
1051C
1052C Returns amplitude squared summed/avg over colors
1053C for the point with external lines W(0:6,NEXTERNAL)
1054C
1055C Process: s u~ > s u~ ta+ ta- / h QED=4
1056C
1057 IMPLICIT NONE
1058C
1059C CONSTANTS
1060C
1061 INTEGER NGRAPHS
1062 parameter(ngraphs=43)
1063 INTEGER NEXTERNAL
1064 parameter(nexternal=6)
1065 INTEGER NWAVEFUNCS, NCOLOR
1066 parameter(nwavefuncs=15, ncolor=2)
1067 real*8 zero
1068 parameter(zero=0d0)
1069 COMPLEX*16 IMAG1
1070 parameter(imag1=(0d0,1d0))
1071C
1072C ARGUMENTS
1073C
1074 real*8 p(0:3,nexternal)
1075 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1076C
1077C LOCAL VARIABLES
1078C
1079 INTEGER I,J
1080 COMPLEX*16 ZTEMP
1081 real*8 denom(ncolor), cf(ncolor,ncolor)
1082 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1083 COMPLEX*16 W(18,NWAVEFUNCS)
1084 COMPLEX*16 DUM0,DUM1
1085 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1086C
1087C GLOBAL VARIABLES
1088C
1089 include 'coupl.inc'
1090C
1091C COLOR DATA
1092C
1093 DATA denom(1)/1/
1094 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1095C 1 T(2,1) T(3,4)
1096 DATA denom(2)/1/
1097 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1098C 1 T(2,4) T(3,1)
1099
1100 INTEGER H1,H2
1101 real*8 matrix
1102 matrix_sux_sux_noh=0.d0
1103 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1104 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1105
1106C ----------
1107C BEGIN CODE
1108C ----------
1109 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1110 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1111 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1112 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1113 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1114 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1115 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1116 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1117 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1118C Amplitude(s) for diagram number 1
1119 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1120 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1121C Amplitude(s) for diagram number 2
1122 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1123 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1124C Amplitude(s) for diagram number 3
1125 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1126 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,11))
1127C Amplitude(s) for diagram number 4
1128 CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_2,amp(4))
1129 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,12))
1130C Amplitude(s) for diagram number 5
1131 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_1,amp(5))
1132C Amplitude(s) for diagram number 6
1133 CALL ffv2_5_0(w(1,4),w(1,11),w(1,10),gc_51,gc_58,amp(6))
1134C Amplitude(s) for diagram number 7
1135 CALL ffv2_3_0(w(1,12),w(1,3),w(1,10),gc_50,gc_58,amp(7))
1136 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,12))
1137 CALL ffv1p0_3(w(1,4),w(1,2),gc_2,zero,zero,w(1,11))
1138 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1139C Amplitude(s) for diagram number 8
1140 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1141 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1142C Amplitude(s) for diagram number 9
1143 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1144 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,14))
1145C Amplitude(s) for diagram number 10
1146 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1147C Amplitude(s) for diagram number 11
1148 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1149 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,13))
1150 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1151C Amplitude(s) for diagram number 12
1152 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1153 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1154C Amplitude(s) for diagram number 13
1155 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1156C Amplitude(s) for diagram number 14
1157 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1158C Amplitude(s) for diagram number 15
1159 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1160 CALL ffv1_1(w(1,2),w(1,12),gc_2,zero,zero,w(1,15))
1161C Amplitude(s) for diagram number 16
1162 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_2,amp(16))
1163 CALL ffv1_2(w(1,4),w(1,12),gc_2,zero,zero,w(1,5))
1164C Amplitude(s) for diagram number 17
1165 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(17))
1166C Amplitude(s) for diagram number 18
1167 CALL ffv2_5_0(w(1,4),w(1,15),w(1,10),gc_51,gc_58,amp(18))
1168C Amplitude(s) for diagram number 19
1169 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(19))
1170 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1171 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1172C Amplitude(s) for diagram number 20
1173 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_2,amp(20))
1174 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1175C Amplitude(s) for diagram number 21
1176 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_2,amp(21))
1177C Amplitude(s) for diagram number 22
1178 CALL ffv2_5_0(w(1,4),w(1,15),w(1,10),gc_51,gc_58,amp(22))
1179C Amplitude(s) for diagram number 23
1180 CALL ffv2_5_0(w(1,12),w(1,2),w(1,10),gc_51,gc_58,amp(23))
1181 CALL ffv2_5_1(w(1,2),w(1,13),gc_51,gc_58,zero,zero,w(1,12))
1182C Amplitude(s) for diagram number 24
1183 CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_2,amp(24))
1184 CALL ffv2_5_2(w(1,4),w(1,13),gc_51,gc_58,zero,zero,w(1,15))
1185C Amplitude(s) for diagram number 25
1186 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_2,amp(25))
1187C Amplitude(s) for diagram number 26
1188 CALL ffv2_5_0(w(1,4),w(1,12),w(1,10),gc_51,gc_58,amp(26))
1189C Amplitude(s) for diagram number 27
1190 CALL ffv2_5_0(w(1,15),w(1,2),w(1,10),gc_51,gc_58,amp(27))
1191 CALL ffv1_2(w(1,1),w(1,11),gc_1,zero,zero,w(1,15))
1192C Amplitude(s) for diagram number 28
1193 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_1,amp(28))
1194 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,12))
1195C Amplitude(s) for diagram number 29
1196 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_1,amp(29))
1197C Amplitude(s) for diagram number 30
1198 CALL ffv2_3_0(w(1,15),w(1,3),w(1,10),gc_50,gc_58,amp(30))
1199 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,15))
1200C Amplitude(s) for diagram number 31
1201 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_1,amp(31))
1202 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1203 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1204C Amplitude(s) for diagram number 32
1205 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_1,amp(32))
1206C Amplitude(s) for diagram number 33
1207 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1208C Amplitude(s) for diagram number 34
1209 CALL ffv2_3_0(w(1,4),w(1,3),w(1,10),gc_50,gc_58,amp(34))
1210C Amplitude(s) for diagram number 35
1211 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1212 CALL ffv2_3_2(w(1,1),w(1,14),gc_50,gc_58,zero,zero,w(1,11))
1213C Amplitude(s) for diagram number 36
1214 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_1,amp(36))
1215C Amplitude(s) for diagram number 37
1216 CALL ffv2_3_0(w(1,12),w(1,3),w(1,14),gc_50,gc_58,amp(37))
1217C Amplitude(s) for diagram number 38
1218 CALL ffv2_3_0(w(1,11),w(1,3),w(1,10),gc_50,gc_58,amp(38))
1219C Amplitude(s) for diagram number 39
1220 CALL ffv2_3_0(w(1,15),w(1,3),w(1,14),gc_50,gc_58,amp(39))
1221 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,14))
1222C Amplitude(s) for diagram number 40
1223 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_2,amp(40))
1224C Amplitude(s) for diagram number 41
1225 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_101,amp(41))
1226C Amplitude(s) for diagram number 42
1227 CALL ffv2_5_0(w(1,14),w(1,2),w(1,10),gc_51,gc_58,amp(42))
1228C Amplitude(s) for diagram number 43
1229 CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_101,amp(43))
1230 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1231 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1232 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1233 $ -amp(40)-amp(41)-amp(42)-amp(43)
1234 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1235 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1236 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1237 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1238 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1239 $ +amp(38)+amp(39)
1240
1241 matrix = 0.d0
1242 DO i = 1, ncolor
1243 ztemp = (0.d0,0.d0)
1244 DO j = 1, ncolor
1245 ztemp = ztemp + cf(j,i)*jamp(j)
1246 ENDDO
1247 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1248 ENDDO
1249 matrix_sux_sux_noh=matrix
1250 ENDIF
1251 ENDIF
1252 END
1253
1254C ----- begin subprocesses SUX->DCX NO Higgs
1255
1256 SUBROUTINE sux_dcx_noh(P,H1,H2,ANS)
1257C
1258C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1259C By the MadGraph Development Team
1260C Please visit us at https://launchpad.net/madgraph5
1261C
1262C MadGraph StandAlone Version
1263C
1264C Returns amplitude squared summed/avg over colors
1265C and helicities
1266C for the point in phase space P(0:3,NEXTERNAL)
1267C
1268C Process: s u~ > d c~ ta+ ta- / h QED=4
1269C
1270 IMPLICIT NONE
1271C
1272C CONSTANTS
1273C
1274 INTEGER NEXTERNAL
1275 parameter(nexternal=6)
1276 INTEGER NCOMB
1277 parameter( ncomb=64)
1278C
1279C ARGUMENTS
1280C
1281 real*8 p(0:3,nexternal),ans
1282 INTEGER H1,H2
1283C
1284C LOCAL VARIABLES
1285C
1286 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1287 real*8 t
1288 real*8 matrix_sux_dcx_noh
1289 INTEGER IHEL,IDEN, I
1290 INTEGER JC(NEXTERNAL)
1291 LOGICAL GOODHEL(NCOMB)
1292 DATA ntry/0/
1293 DATA goodhel/ncomb*.false./
1294 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1295 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1296 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1297 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1298 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1299 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1300 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1301 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1302 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1303 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1304 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1305 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1306 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1307 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1308 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1309 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1310 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1311 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1312 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1313 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1314 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1315 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1316 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1317 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1318 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1319 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1320 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1321 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1322 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1323 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1324 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1325 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1326 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1327 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1328 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1329 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1330 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1331 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1332 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1333 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1334 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1335 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1336 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1337 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1338 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1339 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1340 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1341 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1342 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1343 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1344 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1345 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1346 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1347 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1348 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1349 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1350 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1351 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1352 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1353 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1354 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1355 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1356 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1357 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1358 DATA iden/36/
1359C ----------
1360C BEGIN CODE
1361C ----------
1362 DO ihel=1,nexternal
1363 jc(ihel) = +1
1364 ENDDO
1365 ans = 0d0
1366 DO ihel=1,ncomb
1367 t=matrix_sux_dcx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1368 ans=ans+t
1369 ENDDO
1370 ans=ans/dble(iden)
1371 END
1372
1373
1374 REAL*8 FUNCTION matrix_sux_dcx_noh(P,H1,H2,NHEL,IC)
1375C
1376C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1377C By the MadGraph Development Team
1378C Please visit us at https://launchpad.net/madgraph5
1379C
1380C Returns amplitude squared summed/avg over colors
1381C for the point with external lines W(0:6,NEXTERNAL)
1382C
1383C Process: s u~ > d c~ ta+ ta- / h QED=4
1384C
1385 IMPLICIT NONE
1386C
1387C CONSTANTS
1388C
1389 INTEGER NGRAPHS
1390 parameter(ngraphs=11)
1391 INTEGER NEXTERNAL
1392 parameter(nexternal=6)
1393 INTEGER NWAVEFUNCS, NCOLOR
1394 parameter(nwavefuncs=10, ncolor=1)
1395 real*8 zero
1396 parameter(zero=0d0)
1397 COMPLEX*16 IMAG1
1398 parameter(imag1=(0d0,1d0))
1399C
1400C ARGUMENTS
1401C
1402 real*8 p(0:3,nexternal)
1403 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1404C
1405C LOCAL VARIABLES
1406C
1407 INTEGER I,J
1408 COMPLEX*16 ZTEMP
1409 real*8 denom(ncolor), cf(ncolor,ncolor)
1410 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1411 COMPLEX*16 W(18,NWAVEFUNCS)
1412 COMPLEX*16 DUM0,DUM1
1413 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1414C
1415C GLOBAL VARIABLES
1416C
1417 include 'coupl.inc'
1418C
1419C COLOR DATA
1420C
1421 DATA denom(1)/1/
1422 DATA (cf(i, 1),i= 1, 1) / 9/
1423C 1 T(2,1) T(3,4)
1424
1425 INTEGER H1,H2
1426 real*8 matrix
1427 matrix_sux_dcx_noh=0.d0
1428 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1429 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1430
1431C ----------
1432C BEGIN CODE
1433C ----------
1434 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1435 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1436 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1437 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1438 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1439 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1440 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1441 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1442 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1443C Amplitude(s) for diagram number 1
1444 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1445 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1446C Amplitude(s) for diagram number 2
1447 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1448 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1449C Amplitude(s) for diagram number 3
1450 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1451 CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,6))
1452C Amplitude(s) for diagram number 4
1453 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1454 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
1455C Amplitude(s) for diagram number 5
1456 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1457C Amplitude(s) for diagram number 6
1458 CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1459C Amplitude(s) for diagram number 7
1460 CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(7))
1461 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1462C Amplitude(s) for diagram number 8
1463 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(8))
1464 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,3))
1465C Amplitude(s) for diagram number 9
1466 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
1467C Amplitude(s) for diagram number 10
1468 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
1469 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1470C Amplitude(s) for diagram number 11
1471 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
1472 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1473 $ -amp(9)-amp(10)-amp(11)
1474
1475 matrix = 0.d0
1476 DO i = 1, ncolor
1477 ztemp = (0.d0,0.d0)
1478 DO j = 1, ncolor
1479 ztemp = ztemp + cf(j,i)*jamp(j)
1480 ENDDO
1481 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1482 ENDDO
1483 matrix_sux_dcx_noh=matrix
1484 ENDIF
1485 ENDIF
1486 END
1487
1488C ----- begin subprocesses SUX->SCX NO Higgs
1489C IS THE SAME AS SUX->SCX NO Higgs
1490
1491
1492
1493C DONE 26.10.2015 STILL S=2 MISSING