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