C++ Interface to Tauola
GD.f
1 C JK. done 20.11.2014
2 c------------------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 
11 C 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 
44 c ------- no subprocesses for gluon-d->jjH
45 
46 c --------------- begin subprocesses initiated by gloun-d, no Higgs
47 
48  SUBROUTINE gd_gd_noh(P,H1,H2,ANS)
49 C
50 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
51 C By the MadGraph Development Team
52 C Please visit us at https://launchpad.net/madgraph5
53 C
54 C MadGraph StandAlone Version
55 C
56 C Returns amplitude squared summed/avg over colors
57 C and helicities
58 C for the point in phase space P(0:3,NEXTERNAL)
59 C
60 C Process: g d > g d ta+ ta- / h QED=4
61 C Process: g s > g s ta+ ta- / h QED=4
62 C
63  IMPLICIT NONE
64 C
65 C CONSTANTS
66 C
67  INTEGER NEXTERNAL
68  parameter(nexternal=6)
69  INTEGER NCOMB
70  parameter( ncomb=64)
71 C
72 C ARGUMENTS
73 C
74  real*8 p(0:3,nexternal),ans
75  INTEGER H1,H2
76 C
77 C LOCAL VARIABLES
78 C
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/
152 C ----------
153 C BEGIN CODE
154 C ----------
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)
168 C
169 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
170 C By the MadGraph Development Team
171 C Please visit us at https://launchpad.net/madgraph5
172 C
173 C Returns amplitude squared summed/avg over colors
174 C for the point with external lines W(0:6,NEXTERNAL)
175 C
176 C Process: g d > g d ta+ ta- / h QED=4
177 C Process: g s > g s ta+ ta- / h QED=4
178 C
179  IMPLICIT NONE
180 C
181 C CONSTANTS
182 C
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))
193 C
194 C ARGUMENTS
195 C
196  real*8 p(0:3,nexternal)
197  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
198 C
199 C LOCAL VARIABLES
200 C
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)/
208 C
209 C GLOBAL VARIABLES
210 C
211  include 'coupl.inc'
212 C
213 C COLOR DATA
214 C
215  DATA denom(1)/3/
216  DATA (cf(i, 1),i= 1, 2) / 16, -2/
217 C 1 T(1,3,4,2)
218  DATA denom(2)/3/
219  DATA (cf(i, 2),i= 1, 2) / -2, 16/
220 C 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 
228 C ----------
229 C BEGIN CODE
230 C ----------
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))
240 C 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))
243 C 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))
246 C 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))
248 C 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))
252 C Amplitude(s) for diagram number 5
253  CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_1,amp(5))
254 C 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))
257 C Amplitude(s) for diagram number 7
258  CALL ffv1_0(w(1,2),w(1,6),w(1,8),gc_1,amp(7))
259 C 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))
263 C Amplitude(s) for diagram number 9
264  CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_1,amp(9))
265 C 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))
268 C Amplitude(s) for diagram number 11
269  CALL ffv1_0(w(1,5),w(1,4),w(1,8),gc_1,amp(11))
270 C 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))
273 C Amplitude(s) for diagram number 13
274  CALL ffv1_0(w(1,2),w(1,5),w(1,8),gc_1,amp(13))
275 C 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))
278 C Amplitude(s) for diagram number 15
279  CALL ffv1_0(w(1,5),w(1,4),w(1,8),gc_1,amp(15))
280 C 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 
Generated by doxygen 1.9.1