C++ Interface to Tauola
GG.f
1C JK. done 20.11.2014
2c------------------begin subprocess initiated by GG----------
3
4 SUBROUTINE gg(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
9 real*8 pp(0:3,6)
10
11C I3=1 I4=-1 ! gg->ddbar
12C I3=2 I4=-2 ! gg->uubar
13c I3=I4=0 ! gg->ddbar+uubar
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 ans=0.d0
28 ELSE IF(key.EQ.0) THEN
29 IF(i3.EQ.1 .AND. i4.EQ.-1) CALL gg_ddx_noh(p,h1,h2,ans)
30 IF(i3.EQ.-1 .AND. i4.EQ.1) CALL gg_ddx_noh(pp,h1,h2,ans)
31 IF(i3.EQ.2 .AND. i4.EQ.-2) CALL gg_uux_noh(p,h1,h2,ans)
32 IF(i3.EQ.-2 .AND. i4.EQ.2) CALL gg_uux_noh(pp,h1,h2,ans)
33 IF(i3.EQ.3 .AND. i4.EQ.-3) CALL gg_ddx_noh(p,h1,h2,ans)
34 IF(i3.EQ.-3 .AND. i4.EQ.3) CALL gg_ddx_noh(pp,h1,h2,ans)
35 IF(i3.EQ.4 .AND. i4.EQ.-4) CALL gg_uux_noh(p,h1,h2,ans)
36 IF(i3.EQ.-4 .AND. i4.EQ.4) CALL gg_uux_noh(pp,h1,h2,ans)
37 IF(i3.EQ.5 .AND. i4.EQ.-5) CALL gg_ddx_noh(p,h1,h2,ans)
38 IF(i3.EQ.-5 .AND. i4.EQ.5) CALL gg_ddx_noh(pp,h1,h2,ans)
39
40 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
41 CALL gg_ddx_noh(p,h1,h2,ans1)
42 CALL gg_uux_noh(p,h1,h2,ans2)
43 ans=ans1+ans2
44 CALL gg_ddx_noh(pp,h1,h2,ans1)
45 CALL gg_uux_noh(pp,h1,h2,ans2)
46 ans=ans1+ans2+ans
47 ENDIF
48 ELSE
49 WRITE(*,*) 'SPIN=2 NOT FINISHED'
50 stop
51 ENDIF
52 END SUBROUTINE gg
53
54c ------- no subprocesses for gluon-gluon->jjH
55
56c --------------- begin subprocesses initiated by gloun-gluon, no Higgs
57
58c------- subprocess gg->ddbar, no Higgs
59
60 SUBROUTINE gg_ddx_noh(P,H1,H2,ANS)
61C
62C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
63C By the MadGraph Development Team
64C Please visit us at https://launchpad.net/madgraph5
65C
66C MadGraph StandAlone Version
67C
68C Returns amplitude squared summed/avg over colors
69C and helicities
70C for the point in phase space P(0:3,NEXTERNAL)
71C
72C Process: g g > d d~ ta+ ta- / h QED=4
73C Process: g g > s s~ ta+ ta- / h QED=4
74C
75 IMPLICIT NONE
76C
77C CONSTANTS
78C
79 INTEGER NEXTERNAL
80 parameter(nexternal=6)
81 INTEGER NCOMB
82 parameter( ncomb=64)
83C
84C ARGUMENTS
85C
86 real*8 p(0:3,nexternal),ans
87 INTEGER H1,H2
88C
89C LOCAL VARIABLES
90C
91 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
92 real*8 t
93 real*8 matrix_gg_ddx_noh
94 INTEGER IHEL,IDEN, I
95 INTEGER JC(NEXTERNAL)
96 LOGICAL GOODHEL(NCOMB)
97 DATA ntry/0/
98 DATA goodhel/ncomb*.false./
99 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
100 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
101 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
102 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
103 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
104 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
105 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
106 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
107 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
108 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
109 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
110 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
111 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
112 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
113 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
114 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
115 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
116 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
117 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
118 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
119 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
120 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
121 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
122 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
123 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
124 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
125 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
126 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
127 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
128 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
129 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
130 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
131 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
132 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
133 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
134 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
135 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
136 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
137 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
138 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
139 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
140 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
141 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
142 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
143 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
144 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
145 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
146 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
147 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
148 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
149 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
150 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
151 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
152 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
153 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
154 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
155 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
156 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
157 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
158 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
159 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
160 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
161 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
162 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
163 DATA iden/256/
164C ----------
165C BEGIN CODE
166C ----------
167 DO ihel=1,nexternal
168 jc(ihel) = +1
169 ENDDO
170 ans = 0d0
171 DO ihel=1,ncomb
172 t=matrix_gg_ddx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
173 ans=ans+t
174 ENDDO
175 ans=ans/dble(iden)
176 END
177
178
179 REAL*8 FUNCTION matrix_gg_ddx_noh(P,H1,H2,NHEL,IC)
180C
181C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
182C By the MadGraph Development Team
183C Please visit us at https://launchpad.net/madgraph5
184C
185C Returns amplitude squared summed/avg over colors
186C for the point with external lines W(0:6,NEXTERNAL)
187C
188C Process: g g > d d~ ta+ ta- / h QED=4
189C Process: g g > s s~ ta+ ta- / h QED=4
190C
191 IMPLICIT NONE
192C
193C CONSTANTS
194C
195 INTEGER NGRAPHS
196 parameter(ngraphs=16)
197 INTEGER NEXTERNAL
198 parameter(nexternal=6)
199 INTEGER NWAVEFUNCS, NCOLOR
200 parameter(nwavefuncs=10, ncolor=2)
201 real*8 zero
202 parameter(zero=0d0)
203 COMPLEX*16 IMAG1
204 parameter(imag1=(0d0,1d0))
205C
206C ARGUMENTS
207C
208 real*8 p(0:3,nexternal)
209 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
210C
211C LOCAL VARIABLES
212C
213 INTEGER I,J
214 COMPLEX*16 ZTEMP
215 real*8 denom(ncolor), cf(ncolor,ncolor)
216 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
217 COMPLEX*16 W(18,NWAVEFUNCS)
218 COMPLEX*16 DUM0,DUM1
219 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
220C
221C GLOBAL VARIABLES
222C
223 include 'coupl.inc'
224C
225C COLOR DATA
226C
227 DATA denom(1)/3/
228 DATA (cf(i, 1),i= 1, 2) / 16, -2/
229C 1 T(1,2,3,4)
230 DATA denom(2)/3/
231 DATA (cf(i, 2),i= 1, 2) / -2, 16/
232C 1 T(2,1,3,4)
233
234 INTEGER H1,H2
235 real*8 matrix
236 matrix_gg_ddx_noh=0.d0
237 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
238 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
239
240C ----------
241C BEGIN CODE
242C ----------
243 CALL vxxxxx(p(0,1),zero,nhel(1),-1*ic(1),w(1,1))
244 CALL vxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
245 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
246 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
247 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
248 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
249 CALL vvv1p0_1(w(1,1),w(1,2),gc_10,zero,zero,w(1,7))
250 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,8))
251 CALL ffv1_1(w(1,3),w(1,7),gc_11,zero,zero,w(1,9))
252C Amplitude(s) for diagram number 1
253 CALL ffv1_0(w(1,4),w(1,9),w(1,8),gc_1,amp(1))
254 CALL ffv1_2(w(1,4),w(1,7),gc_11,zero,zero,w(1,10))
255C Amplitude(s) for diagram number 2
256 CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_1,amp(2))
257 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
258C Amplitude(s) for diagram number 3
259 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(3))
260C Amplitude(s) for diagram number 4
261 CALL ffv2_3_0(w(1,10),w(1,3),w(1,7),gc_50,gc_58,amp(4))
262 CALL ffv1_1(w(1,3),w(1,1),gc_11,zero,zero,w(1,10))
263 CALL ffv1_2(w(1,4),w(1,2),gc_11,zero,zero,w(1,9))
264C Amplitude(s) for diagram number 5
265 CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_1,amp(5))
266C Amplitude(s) for diagram number 6
267 CALL ffv2_3_0(w(1,9),w(1,10),w(1,7),gc_50,gc_58,amp(6))
268 CALL ffv1_1(w(1,10),w(1,2),gc_11,zero,zero,w(1,6))
269C Amplitude(s) for diagram number 7
270 CALL ffv1_0(w(1,4),w(1,6),w(1,8),gc_1,amp(7))
271C Amplitude(s) for diagram number 8
272 CALL ffv2_3_0(w(1,4),w(1,6),w(1,7),gc_50,gc_58,amp(8))
273 CALL ffv1_2(w(1,4),w(1,1),gc_11,zero,zero,w(1,6))
274 CALL ffv1_1(w(1,3),w(1,2),gc_11,zero,zero,w(1,10))
275C Amplitude(s) for diagram number 9
276 CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_1,amp(9))
277C Amplitude(s) for diagram number 10
278 CALL ffv2_3_0(w(1,6),w(1,10),w(1,7),gc_50,gc_58,amp(10))
279 CALL ffv1_2(w(1,6),w(1,2),gc_11,zero,zero,w(1,5))
280C Amplitude(s) for diagram number 11
281 CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_1,amp(11))
282C Amplitude(s) for diagram number 12
283 CALL ffv2_3_0(w(1,5),w(1,3),w(1,7),gc_50,gc_58,amp(12))
284 CALL ffv1_1(w(1,10),w(1,1),gc_11,zero,zero,w(1,5))
285C Amplitude(s) for diagram number 13
286 CALL ffv1_0(w(1,4),w(1,5),w(1,8),gc_1,amp(13))
287C Amplitude(s) for diagram number 14
288 CALL ffv2_3_0(w(1,4),w(1,5),w(1,7),gc_50,gc_58,amp(14))
289 CALL ffv1_2(w(1,9),w(1,1),gc_11,zero,zero,w(1,5))
290C Amplitude(s) for diagram number 15
291 CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_1,amp(15))
292C Amplitude(s) for diagram number 16
293 CALL ffv2_3_0(w(1,5),w(1,3),w(1,7),gc_50,gc_58,amp(16))
294 jamp(1)=+imag1*amp(1)+imag1*amp(2)+imag1*amp(3)+imag1*amp(4)
295 $ -amp(5)-amp(6)-amp(7)-amp(8)-amp(15)-amp(16)
296 jamp(2)=-imag1*amp(1)-imag1*amp(2)-imag1*amp(3)-imag1*amp(4)
297 $ -amp(9)-amp(10)-amp(11)-amp(12)-amp(13)-amp(14)
298
299 matrix = 0.d0
300 DO i = 1, ncolor
301 ztemp = (0.d0,0.d0)
302 DO j = 1, ncolor
303 ztemp = ztemp + cf(j,i)*jamp(j)
304 ENDDO
305 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
306 ENDDO
307 matrix_gg_ddx_noh=matrix
308 ENDIF
309 ENDIF
310 END
311
312
313c--------subprocess gg->uubar, no Higgs
314
315 SUBROUTINE gg_uux_noh(P,H1,H2,ANS)
316C
317C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
318C By the MadGraph Development Team
319C Please visit us at https://launchpad.net/madgraph5
320C
321C MadGraph StandAlone Version
322C
323C Returns amplitude squared summed/avg over colors
324C and helicities
325C for the point in phase space P(0:3,NEXTERNAL)
326C
327C Process: g g > u u~ ta+ ta- / h QED=4
328C Process: g g > c c~ ta+ ta- / h QED=4
329C
330 IMPLICIT NONE
331C
332C CONSTANTS
333C
334 INTEGER NEXTERNAL
335 parameter(nexternal=6)
336 INTEGER NCOMB
337 parameter( ncomb=64)
338C
339C ARGUMENTS
340C
341 real*8 p(0:3,nexternal),ans
342 INTEGER H1,H2
343C
344C LOCAL VARIABLES
345C
346 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
347 real*8 t
348 real*8 matrix_gg_uux_noh
349 INTEGER IHEL,IDEN, I
350 INTEGER JC(NEXTERNAL)
351 LOGICAL GOODHEL(NCOMB)
352 DATA ntry/0/
353 DATA goodhel/ncomb*.false./
354 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
355 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
356 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
357 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
358 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
359 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
360 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
361 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
362 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
363 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
364 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
365 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
366 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
367 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
368 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
369 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
370 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
371 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
372 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
373 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
374 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
375 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
376 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
377 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
378 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
379 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
380 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
381 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
382 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
383 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
384 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
385 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
386 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
387 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
388 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
389 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
390 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
391 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
392 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
393 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
394 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
395 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
396 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
397 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
398 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
399 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
400 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
401 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
402 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
403 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
404 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
405 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
406 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
407 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
408 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
409 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
410 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
411 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
412 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
413 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
414 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
415 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
416 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
417 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
418 DATA iden/256/
419C ----------
420C BEGIN CODE
421C ----------
422 DO ihel=1,nexternal
423 jc(ihel) = +1
424 ENDDO
425 ans = 0d0
426 DO ihel=1,ncomb
427 t=matrix_gg_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
428 ans=ans+t
429 ENDDO
430 ans=ans/dble(iden)
431 END
432
433
434 REAL*8 FUNCTION matrix_gg_uux_noh(P,H1,H2,NHEL,IC)
435C
436C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
437C By the MadGraph Development Team
438C Please visit us at https://launchpad.net/madgraph5
439C
440C Returns amplitude squared summed/avg over colors
441C for the point with external lines W(0:6,NEXTERNAL)
442C
443C Process: g g > u u~ ta+ ta- / h QED=4
444C Process: g g > c c~ ta+ ta- / h QED=4
445C
446 IMPLICIT NONE
447C
448C CONSTANTS
449C
450 INTEGER NGRAPHS
451 parameter(ngraphs=16)
452 INTEGER NEXTERNAL
453 parameter(nexternal=6)
454 INTEGER NWAVEFUNCS, NCOLOR
455 parameter(nwavefuncs=10, ncolor=2)
456 real*8 zero
457 parameter(zero=0d0)
458 COMPLEX*16 IMAG1
459 parameter(imag1=(0d0,1d0))
460C
461C ARGUMENTS
462C
463 real*8 p(0:3,nexternal)
464 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
465C
466C LOCAL VARIABLES
467C
468 INTEGER I,J
469 COMPLEX*16 ZTEMP
470 real*8 denom(ncolor), cf(ncolor,ncolor)
471 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
472 COMPLEX*16 W(18,NWAVEFUNCS)
473 COMPLEX*16 DUM0,DUM1
474 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
475C
476C GLOBAL VARIABLES
477C
478 include 'coupl.inc'
479C
480C COLOR DATA
481C
482 DATA denom(1)/3/
483 DATA (cf(i, 1),i= 1, 2) / 16, -2/
484C 1 T(1,2,3,4)
485 DATA denom(2)/3/
486 DATA (cf(i, 2),i= 1, 2) / -2, 16/
487C 1 T(2,1,3,4)
488
489 INTEGER H1,H2
490 real*8 matrix
491 matrix_gg_uux_noh=0.d0
492 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
493 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
494
495C ----------
496C BEGIN CODE
497C ----------
498 CALL vxxxxx(p(0,1),zero,nhel(1),-1*ic(1),w(1,1))
499 CALL vxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
500 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
501 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
502 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
503 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
504 CALL vvv1p0_1(w(1,1),w(1,2),gc_10,zero,zero,w(1,7))
505 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,8))
506 CALL ffv1_1(w(1,3),w(1,7),gc_11,zero,zero,w(1,9))
507C Amplitude(s) for diagram number 1
508 CALL ffv1_0(w(1,4),w(1,9),w(1,8),gc_2,amp(1))
509 CALL ffv1_2(w(1,4),w(1,7),gc_11,zero,zero,w(1,10))
510C Amplitude(s) for diagram number 2
511 CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_2,amp(2))
512 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
513C Amplitude(s) for diagram number 3
514 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(3))
515C Amplitude(s) for diagram number 4
516 CALL ffv2_5_0(w(1,10),w(1,3),w(1,7),gc_51,gc_58,amp(4))
517 CALL ffv1_1(w(1,3),w(1,1),gc_11,zero,zero,w(1,10))
518 CALL ffv1_2(w(1,4),w(1,2),gc_11,zero,zero,w(1,9))
519C Amplitude(s) for diagram number 5
520 CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_2,amp(5))
521C Amplitude(s) for diagram number 6
522 CALL ffv2_5_0(w(1,9),w(1,10),w(1,7),gc_51,gc_58,amp(6))
523 CALL ffv1_1(w(1,10),w(1,2),gc_11,zero,zero,w(1,6))
524C Amplitude(s) for diagram number 7
525 CALL ffv1_0(w(1,4),w(1,6),w(1,8),gc_2,amp(7))
526C Amplitude(s) for diagram number 8
527 CALL ffv2_5_0(w(1,4),w(1,6),w(1,7),gc_51,gc_58,amp(8))
528 CALL ffv1_2(w(1,4),w(1,1),gc_11,zero,zero,w(1,6))
529 CALL ffv1_1(w(1,3),w(1,2),gc_11,zero,zero,w(1,10))
530C Amplitude(s) for diagram number 9
531 CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_2,amp(9))
532C Amplitude(s) for diagram number 10
533 CALL ffv2_5_0(w(1,6),w(1,10),w(1,7),gc_51,gc_58,amp(10))
534 CALL ffv1_2(w(1,6),w(1,2),gc_11,zero,zero,w(1,5))
535C Amplitude(s) for diagram number 11
536 CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_2,amp(11))
537C Amplitude(s) for diagram number 12
538 CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(12))
539 CALL ffv1_1(w(1,10),w(1,1),gc_11,zero,zero,w(1,5))
540C Amplitude(s) for diagram number 13
541 CALL ffv1_0(w(1,4),w(1,5),w(1,8),gc_2,amp(13))
542C Amplitude(s) for diagram number 14
543 CALL ffv2_5_0(w(1,4),w(1,5),w(1,7),gc_51,gc_58,amp(14))
544 CALL ffv1_2(w(1,9),w(1,1),gc_11,zero,zero,w(1,5))
545C Amplitude(s) for diagram number 15
546 CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_2,amp(15))
547C Amplitude(s) for diagram number 16
548 CALL ffv2_5_0(w(1,5),w(1,3),w(1,7),gc_51,gc_58,amp(16))
549 jamp(1)=+imag1*amp(1)+imag1*amp(2)+imag1*amp(3)+imag1*amp(4)
550 $ -amp(5)-amp(6)-amp(7)-amp(8)-amp(15)-amp(16)
551 jamp(2)=-imag1*amp(1)-imag1*amp(2)-imag1*amp(3)-imag1*amp(4)
552 $ -amp(9)-amp(10)-amp(11)-amp(12)-amp(13)-amp(14)
553
554 matrix = 0.d0
555 DO i = 1, ncolor
556 ztemp = (0.d0,0.d0)
557 DO j = 1, ncolor
558 ztemp = ztemp + cf(j,i)*jamp(j)
559 ENDDO
560 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
561 ENDDO
562 matrix_gg_uux_noh=matrix
563 ENDIF
564 ENDIF
565 END