C++ Interface to Tauola
GU.f
1C JK. done 20.11.2014
2c------------------begin subprocess initiated by GU----------
3
4 SUBROUTINE gu(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=0,21 GLOUN I4=0,2 ! only gu->gu
12
13 ans=0.d0
14
15 IF(key.EQ.1) THEN
16 ans=0.d0
17 ELSE IF(key.EQ.0) THEN
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 IF(i3.EQ.21 .AND. i4.EQ.2) CALL gu_gu_noh(p,h1,h2,ans)
26 IF(i3.EQ.2 .AND. i4.EQ.21) CALL gu_gu_noh(pp,h1,h2,ans)
27 IF(i3.EQ.21 .AND. i4.EQ.4) CALL gu_gu_noh(p,h1,h2,ans)
28 IF(i3.EQ.4 .AND. i4.EQ.21) CALL gu_gu_noh(pp,h1,h2,ans)
29 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
30 CALL gu_gu_noh(p,h1,h2,ans1)
31 CALL gu_gu_noh(pp,h1,h2,ans2)
32 ans=ans1+ans2
33 ENDIF
34 ELSE
35 WRITE(*,*) 'NOT FINISHED'
36 stop
37 ENDIF
38 END SUBROUTINE gu
39
40c ------- no subprocesses for gluon-U->jjH
41
42c --------------- begin subprocesses initiated by gloun-U, no Higgs
43
44 SUBROUTINE gu_gu_noh(P,H1,H2,ANS)
45C
46C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
47C By the MadGraph Development Team
48C Please visit us at https://launchpad.net/madgraph5
49C
50C MadGraph StandAlone Version
51C
52C Returns amplitude squared summed/avg over colors
53C and helicities
54C for the point in phase space P(0:3,NEXTERNAL)
55C
56C Process: g u > g u ta+ ta- / h QED=4
57C Process: g c > g c ta+ ta- / h QED=4
58C
59 IMPLICIT NONE
60C
61C CONSTANTS
62C
63 INTEGER NEXTERNAL
64 parameter(nexternal=6)
65 INTEGER NCOMB
66 parameter( ncomb=64)
67C
68C ARGUMENTS
69C
70 real*8 p(0:3,nexternal),ans
71 INTEGER H1,H2
72C
73C LOCAL VARIABLES
74C
75 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
76 real*8 t
77 real*8 matrix_gu_gu_noh
78 INTEGER IHEL,IDEN, I
79 INTEGER JC(NEXTERNAL)
80 LOGICAL GOODHEL(NCOMB)
81 DATA ntry/0/
82 DATA goodhel/ncomb*.false./
83 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
84 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
85 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
86 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
87 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
88 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
89 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
90 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
91 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
92 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
93 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
94 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
95 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
96 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
97 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
98 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
99 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
100 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
101 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
102 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
103 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
104 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
105 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
106 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
107 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
108 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
109 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
110 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
111 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
112 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
113 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
114 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
115 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
116 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
117 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
118 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
119 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
120 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
121 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
122 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
123 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
124 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
125 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
126 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
127 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
128 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
129 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
130 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
131 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
132 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
133 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
134 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
135 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
136 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
137 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
138 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
139 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
140 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
141 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
142 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
143 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
144 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
145 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
146 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
147 DATA iden/96/
148C ----------
149C BEGIN CODE
150C ----------
151 DO ihel=1,nexternal
152 jc(ihel) = +1
153 ENDDO
154 ans = 0d0
155 DO ihel=1,ncomb
156 t=matrix_gu_gu_noh(p ,h1,h2,nhel(1,ihel),jc(1))
157 ans=ans+t
158 ENDDO
159 ans=ans/dble(iden)
160 END
161
162
163 REAL*8 FUNCTION matrix_gu_gu_noh(P,H1,H2,NHEL,IC)
164C
165C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
166C By the MadGraph Development Team
167C Please visit us at https://launchpad.net/madgraph5
168C
169C Returns amplitude squared summed/avg over colors
170C for the point with external lines W(0:6,NEXTERNAL)
171C
172C Process: g u > g u ta+ ta- / h QED=4
173C Process: g c > g c ta+ ta- / h QED=4
174C
175 IMPLICIT NONE
176C
177C CONSTANTS
178C
179 INTEGER NGRAPHS
180 parameter(ngraphs=16)
181 INTEGER NEXTERNAL
182 parameter(nexternal=6)
183 INTEGER NWAVEFUNCS, NCOLOR
184 parameter(nwavefuncs=10, ncolor=2)
185 real*8 zero
186 parameter(zero=0d0)
187 COMPLEX*16 IMAG1
188 parameter(imag1=(0d0,1d0))
189C
190C ARGUMENTS
191C
192 real*8 p(0:3,nexternal)
193 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
194C
195C LOCAL VARIABLES
196C
197 INTEGER I,J
198 COMPLEX*16 ZTEMP
199 real*8 denom(ncolor), cf(ncolor,ncolor)
200 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
201 COMPLEX*16 W(18,NWAVEFUNCS)
202 COMPLEX*16 DUM0,DUM1
203 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
204C
205C GLOBAL VARIABLES
206C
207 include 'coupl.inc'
208C
209C COLOR DATA
210C
211 DATA denom(1)/3/
212 DATA (cf(i, 1),i= 1, 2) / 16, -2/
213C 1 T(1,3,4,2)
214 DATA denom(2)/3/
215 DATA (cf(i, 2),i= 1, 2) / -2, 16/
216C 1 T(3,1,4,2)
217
218 INTEGER H1,H2
219 real*8 matrix
220 matrix_gu_gu_noh=0.d0
221 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
222 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
223
224C ----------
225C BEGIN CODE
226C ----------
227 CALL vxxxxx(p(0,1),zero,nhel(1),-1*ic(1),w(1,1))
228 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
229 CALL vxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
230 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
231 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
232 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
233 CALL vvv1p0_1(w(1,1),w(1,3),gc_10,zero,zero,w(1,7))
234 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,8))
235 CALL ffv1_1(w(1,4),w(1,7),gc_11,zero,zero,w(1,9))
236C Amplitude(s) for diagram number 1
237 CALL ffv1_0(w(1,2),w(1,9),w(1,8),gc_2,amp(1))
238 CALL ffv1_2(w(1,2),w(1,7),gc_11,zero,zero,w(1,10))
239C Amplitude(s) for diagram number 2
240 CALL ffv1_0(w(1,10),w(1,4),w(1,8),gc_2,amp(2))
241 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
242C Amplitude(s) for diagram number 3
243 CALL ffv2_5_0(w(1,2),w(1,9),w(1,7),gc_51,gc_58,amp(3))
244C Amplitude(s) for diagram number 4
245 CALL ffv2_5_0(w(1,10),w(1,4),w(1,7),gc_51,gc_58,amp(4))
246 CALL ffv1_1(w(1,4),w(1,1),gc_11,zero,zero,w(1,10))
247 CALL ffv1_2(w(1,2),w(1,3),gc_11,zero,zero,w(1,9))
248C Amplitude(s) for diagram number 5
249 CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_2,amp(5))
250C Amplitude(s) for diagram number 6
251 CALL ffv2_5_0(w(1,9),w(1,10),w(1,7),gc_51,gc_58,amp(6))
252 CALL ffv1_1(w(1,10),w(1,3),gc_11,zero,zero,w(1,6))
253C Amplitude(s) for diagram number 7
254 CALL ffv1_0(w(1,2),w(1,6),w(1,8),gc_2,amp(7))
255C Amplitude(s) for diagram number 8
256 CALL ffv2_5_0(w(1,2),w(1,6),w(1,7),gc_51,gc_58,amp(8))
257 CALL ffv1_2(w(1,2),w(1,1),gc_11,zero,zero,w(1,6))
258 CALL ffv1_1(w(1,4),w(1,3),gc_11,zero,zero,w(1,10))
259C Amplitude(s) for diagram number 9
260 CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_2,amp(9))
261C Amplitude(s) for diagram number 10
262 CALL ffv2_5_0(w(1,6),w(1,10),w(1,7),gc_51,gc_58,amp(10))
263 CALL ffv1_2(w(1,6),w(1,3),gc_11,zero,zero,w(1,5))
264C Amplitude(s) for diagram number 11
265 CALL ffv1_0(w(1,5),w(1,4),w(1,8),gc_2,amp(11))
266C Amplitude(s) for diagram number 12
267 CALL ffv2_5_0(w(1,5),w(1,4),w(1,7),gc_51,gc_58,amp(12))
268 CALL ffv1_1(w(1,10),w(1,1),gc_11,zero,zero,w(1,5))
269C Amplitude(s) for diagram number 13
270 CALL ffv1_0(w(1,2),w(1,5),w(1,8),gc_2,amp(13))
271C Amplitude(s) for diagram number 14
272 CALL ffv2_5_0(w(1,2),w(1,5),w(1,7),gc_51,gc_58,amp(14))
273 CALL ffv1_2(w(1,9),w(1,1),gc_11,zero,zero,w(1,5))
274C Amplitude(s) for diagram number 15
275 CALL ffv1_0(w(1,5),w(1,4),w(1,8),gc_2,amp(15))
276C Amplitude(s) for diagram number 16
277 CALL ffv2_5_0(w(1,5),w(1,4),w(1,7),gc_51,gc_58,amp(16))
278 jamp(1)=-imag1*amp(1)-imag1*amp(2)-imag1*amp(3)-imag1*amp(4)
279 $ +amp(5)+amp(6)+amp(7)+amp(8)+amp(15)+amp(16)
280 jamp(2)=+imag1*amp(1)+imag1*amp(2)+imag1*amp(3)+imag1*amp(4)
281 $ +amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
282
283 matrix = 0.d0
284 DO i = 1, ncolor
285 ztemp = (0.d0,0.d0)
286 DO j = 1, ncolor
287 ztemp = ztemp + cf(j,i)*jamp(j)
288 ENDDO
289 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
290 ENDDO
291 matrix_gu_gu_noh=matrix
292 ENDIF
293 ENDIF
294 END
295
296