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