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