C++InterfacetoTauola
GG.f
1 C JK. done 20.11.2014
2 c------------------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 
11 C I3=1 I4=-1 ! gg->ddbar
12 C I3=2 I4=-2 ! gg->uubar
13 c 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 
54 c ------- no subprocesses for gluon-gluon->jjH
55 
56 c --------------- begin subprocesses initiated by gloun-gluon, no Higgs
57 
58 c------- subprocess gg->ddbar, no Higgs
59 
60  SUBROUTINE gg_ddx_noh(P,H1,H2,ANS)
61 C
62 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
63 C By the MadGraph Development Team
64 C Please visit us at https://launchpad.net/madgraph5
65 C
66 C MadGraph StandAlone Version
67 C
68 C Returns amplitude squared summed/avg over colors
69 C and helicities
70 C for the point in phase space P(0:3,NEXTERNAL)
71 C
72 C Process: g g > d d~ ta+ ta- / h QED=4
73 C Process: g g > s s~ ta+ ta- / h QED=4
74 C
75  IMPLICIT NONE
76 C
77 C CONSTANTS
78 C
79  INTEGER nexternal
80  parameter(nexternal=6)
81  INTEGER ncomb
82  parameter( ncomb=64)
83 C
84 C ARGUMENTS
85 C
86  REAL*8 p(0:3,nexternal),ans
87  INTEGER h1,h2
88 C
89 C LOCAL VARIABLES
90 C
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/
164 C ----------
165 C BEGIN CODE
166 C ----------
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)
180 C
181 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
182 C By the MadGraph Development Team
183 C Please visit us at https://launchpad.net/madgraph5
184 C
185 C Returns amplitude squared summed/avg over colors
186 C for the point with external lines W(0:6,NEXTERNAL)
187 C
188 C Process: g g > d d~ ta+ ta- / h QED=4
189 C Process: g g > s s~ ta+ ta- / h QED=4
190 C
191  IMPLICIT NONE
192 C
193 C CONSTANTS
194 C
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))
205 C
206 C ARGUMENTS
207 C
208  REAL*8 p(0:3,nexternal)
209  INTEGER nhel(nexternal), ic(nexternal)
210 C
211 C LOCAL VARIABLES
212 C
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)/
220 C
221 C GLOBAL VARIABLES
222 C
223  include 'coupl.inc'
224 C
225 C COLOR DATA
226 C
227  DATA denom(1)/3/
228  DATA (cf(i, 1),i= 1, 2) / 16, -2/
229 C 1 T(1,2,3,4)
230  DATA denom(2)/3/
231  DATA (cf(i, 2),i= 1, 2) / -2, 16/
232 C 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 
240 C ----------
241 C BEGIN CODE
242 C ----------
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))
252 C 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))
255 C 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))
258 C 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))
260 C 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))
264 C Amplitude(s) for diagram number 5
265  CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_1,amp(5))
266 C 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))
269 C Amplitude(s) for diagram number 7
270  CALL ffv1_0(w(1,4),w(1,6),w(1,8),gc_1,amp(7))
271 C 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))
275 C Amplitude(s) for diagram number 9
276  CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_1,amp(9))
277 C 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))
280 C Amplitude(s) for diagram number 11
281  CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_1,amp(11))
282 C 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))
285 C Amplitude(s) for diagram number 13
286  CALL ffv1_0(w(1,4),w(1,5),w(1,8),gc_1,amp(13))
287 C 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))
290 C Amplitude(s) for diagram number 15
291  CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_1,amp(15))
292 C 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 
313 c--------subprocess gg->uubar, no Higgs
314 
315  SUBROUTINE gg_uux_noh(P,H1,H2,ANS)
316 C
317 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
318 C By the MadGraph Development Team
319 C Please visit us at https://launchpad.net/madgraph5
320 C
321 C MadGraph StandAlone Version
322 C
323 C Returns amplitude squared summed/avg over colors
324 C and helicities
325 C for the point in phase space P(0:3,NEXTERNAL)
326 C
327 C Process: g g > u u~ ta+ ta- / h QED=4
328 C Process: g g > c c~ ta+ ta- / h QED=4
329 C
330  IMPLICIT NONE
331 C
332 C CONSTANTS
333 C
334  INTEGER nexternal
335  parameter(nexternal=6)
336  INTEGER ncomb
337  parameter( ncomb=64)
338 C
339 C ARGUMENTS
340 C
341  REAL*8 p(0:3,nexternal),ans
342  INTEGER h1,h2
343 C
344 C LOCAL VARIABLES
345 C
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/
419 C ----------
420 C BEGIN CODE
421 C ----------
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)
435 C
436 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
437 C By the MadGraph Development Team
438 C Please visit us at https://launchpad.net/madgraph5
439 C
440 C Returns amplitude squared summed/avg over colors
441 C for the point with external lines W(0:6,NEXTERNAL)
442 C
443 C Process: g g > u u~ ta+ ta- / h QED=4
444 C Process: g g > c c~ ta+ ta- / h QED=4
445 C
446  IMPLICIT NONE
447 C
448 C CONSTANTS
449 C
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))
460 C
461 C ARGUMENTS
462 C
463  REAL*8 p(0:3,nexternal)
464  INTEGER nhel(nexternal), ic(nexternal)
465 C
466 C LOCAL VARIABLES
467 C
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)/
475 C
476 C GLOBAL VARIABLES
477 C
478  include 'coupl.inc'
479 C
480 C COLOR DATA
481 C
482  DATA denom(1)/3/
483  DATA (cf(i, 1),i= 1, 2) / 16, -2/
484 C 1 T(1,2,3,4)
485  DATA denom(2)/3/
486  DATA (cf(i, 2),i= 1, 2) / -2, 16/
487 C 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 
495 C ----------
496 C BEGIN CODE
497 C ----------
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))
507 C 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))
510 C 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))
513 C 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))
515 C 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))
519 C Amplitude(s) for diagram number 5
520  CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_2,amp(5))
521 C 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))
524 C Amplitude(s) for diagram number 7
525  CALL ffv1_0(w(1,4),w(1,6),w(1,8),gc_2,amp(7))
526 C 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))
530 C Amplitude(s) for diagram number 9
531  CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_2,amp(9))
532 C 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))
535 C Amplitude(s) for diagram number 11
536  CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_2,amp(11))
537 C 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))
540 C Amplitude(s) for diagram number 13
541  CALL ffv1_0(w(1,4),w(1,5),w(1,8),gc_2,amp(13))
542 C 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))
545 C Amplitude(s) for diagram number 15
546  CALL ffv1_0(w(1,5),w(1,3),w(1,8),gc_2,amp(15))
547 C 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
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7