C++ Interface to Tauola
GU.f
1 C JK. done 20.11.2014
2 c------------------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 
11 C 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 
40 c ------- no subprocesses for gluon-U->jjH
41 
42 c --------------- begin subprocesses initiated by gloun-U, no Higgs
43 
44  SUBROUTINE gu_gu_noh(P,H1,H2,ANS)
45 C
46 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
47 C By the MadGraph Development Team
48 C Please visit us at https://launchpad.net/madgraph5
49 C
50 C MadGraph StandAlone Version
51 C
52 C Returns amplitude squared summed/avg over colors
53 C and helicities
54 C for the point in phase space P(0:3,NEXTERNAL)
55 C
56 C Process: g u > g u ta+ ta- / h QED=4
57 C Process: g c > g c ta+ ta- / h QED=4
58 C
59  IMPLICIT NONE
60 C
61 C CONSTANTS
62 C
63  INTEGER NEXTERNAL
64  parameter(nexternal=6)
65  INTEGER NCOMB
66  parameter( ncomb=64)
67 C
68 C ARGUMENTS
69 C
70  real*8 p(0:3,nexternal),ans
71  INTEGER H1,H2
72 C
73 C LOCAL VARIABLES
74 C
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/
148 C ----------
149 C BEGIN CODE
150 C ----------
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)
164 C
165 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
166 C By the MadGraph Development Team
167 C Please visit us at https://launchpad.net/madgraph5
168 C
169 C Returns amplitude squared summed/avg over colors
170 C for the point with external lines W(0:6,NEXTERNAL)
171 C
172 C Process: g u > g u ta+ ta- / h QED=4
173 C Process: g c > g c ta+ ta- / h QED=4
174 C
175  IMPLICIT NONE
176 C
177 C CONSTANTS
178 C
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))
189 C
190 C ARGUMENTS
191 C
192  real*8 p(0:3,nexternal)
193  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
194 C
195 C LOCAL VARIABLES
196 C
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)/
204 C
205 C GLOBAL VARIABLES
206 C
207  include 'coupl.inc'
208 C
209 C COLOR DATA
210 C
211  DATA denom(1)/3/
212  DATA (cf(i, 1),i= 1, 2) / 16, -2/
213 C 1 T(1,3,4,2)
214  DATA denom(2)/3/
215  DATA (cf(i, 2),i= 1, 2) / -2, 16/
216 C 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 
224 C ----------
225 C BEGIN CODE
226 C ----------
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))
236 C 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))
239 C 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))
242 C 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))
244 C 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))
248 C Amplitude(s) for diagram number 5
249  CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_2,amp(5))
250 C 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))
253 C Amplitude(s) for diagram number 7
254  CALL ffv1_0(w(1,2),w(1,6),w(1,8),gc_2,amp(7))
255 C 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))
259 C Amplitude(s) for diagram number 9
260  CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_2,amp(9))
261 C 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))
264 C Amplitude(s) for diagram number 11
265  CALL ffv1_0(w(1,5),w(1,4),w(1,8),gc_2,amp(11))
266 C 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))
269 C Amplitude(s) for diagram number 13
270  CALL ffv1_0(w(1,2),w(1,5),w(1,8),gc_2,amp(13))
271 C 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))
274 C Amplitude(s) for diagram number 15
275  CALL ffv1_0(w(1,5),w(1,4),w(1,8),gc_2,amp(15))
276 C 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 
Generated by doxygen 1.9.1