C++ Interface to Tauola
SU.f
1
2c --------------------- subprocesses initiated by SU: 27.10.2015
3
4 SUBROUTINE su(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,ans3,ans4
9
10C I3=0,2 OR 4 I4=0,1 OR 3
11
12 ans=0.d0
13
14 IF(key.EQ.1) THEN
15 IF(i3.EQ.2 .AND. i4.EQ.1) CALL su_ud_h(p,h1,h2,ans)
16 IF(i3.EQ.2 .AND. i4.EQ.3) CALL su_us_h(p,h1,h2,ans)
17 IF(i3.EQ.4 .AND. i4.EQ.1) CAll su_cd_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.3) CALL su_cs_h(p,h1,h2,ans)
19 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
20 CALL su_ud_h(p,h1,h2,ans1)
21 CALL su_us_h(p,h1,h2,ans2)
22 CALL su_cd_h(p,h1,h2,ans3)
23 CALL su_cs_h(p,h1,h2,ans4)
24 ans=ans1+ans2+ans3+ans4
25 ENDIF
26 ELSE IF(key.EQ.0) THEN
27 IF(i3.EQ.2 .AND. i4.EQ.1) CALL su_ud_noh(p,h1,h2,ans)
28 IF(i3.EQ.2 .AND. i4.EQ.3) CALL su_us_noh(p,h1,h2,ans)
29 IF(i3.EQ.4 .AND. i4.EQ.1) CAll su_cd_noh(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.3) CALL su_cs_noh(p,h1,h2,ans)
31 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
32 CALL su_ud_noh(p,h1,h2,ans1)
33 CALL su_us_noh(p,h1,h2,ans2)
34 CALL su_cd_noh(p,h1,h2,ans3)
35 CALL su_cs_noh(p,h1,h2,ans4)
36 ans=ans1+ans2+ans3+ans4
37 ENDIF
38 ELSE
39 WRITE(*,*) 'NOT FINISHED'
40 stop
41 ENDIF
42 END ! SUBROUTINE SU
43
44
45
46C ----- begin subprocesses US->jjH with Higgs->tautau
47
48 SUBROUTINE su_cd_h(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: s u > c d h WEIGHTED=6
61C * Decay: h > ta+ ta- WEIGHTED=2
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_su_cd_h
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/36/
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_su_cd_h(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_su_cd_h(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: s u > c d h WEIGHTED=6
177C * Decay: h > ta+ ta- WEIGHTED=2
178C
179 IMPLICIT NONE
180C
181C CONSTANTS
182C
183 INTEGER NGRAPHS
184 parameter(ngraphs=1)
185 INTEGER NEXTERNAL
186 parameter(nexternal=6)
187 INTEGER NWAVEFUNCS, NCOLOR
188 parameter(nwavefuncs=7, ncolor=1)
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)/1/
216 DATA (cf(i, 1),i= 1, 1) / 9/
217C 1 T(3,1) T(4,2)
218
219 INTEGER H1,H2
220 real*8 matrix
221 matrix_su_cd_h=0.d0
222 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
223 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
224
225C ----------
226C BEGIN CODE
227C ----------
228 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
229 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
230 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
231 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
232 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
233 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
234 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
235 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
236 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,3))
237C Amplitude(s) for diagram number 1
238 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
239 jamp(1)=-amp(1)
240
241 matrix = 0.d0
242 DO i = 1, ncolor
243 ztemp = (0.d0,0.d0)
244 DO j = 1, ncolor
245 ztemp = ztemp + cf(j,i)*jamp(j)
246 ENDDO
247 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
248 ENDDO
249 matrix_su_cd_h=matrix
250 ENDIF ! CLOSES H1 IF
251 ENDIF ! CLOSES H2 IF
252 END
253
254c==========================================
255
256 SUBROUTINE su_cs_h(P,H1,H2,ANS)
257C
258C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
259C By the MadGraph Development Team
260C Please visit us at https://launchpad.net/madgraph5
261C
262C MadGraph StandAlone Version
263C
264C Returns amplitude squared summed/avg over colors
265C and helicities
266C for the point in phase space P(0:3,NEXTERNAL)
267C
268C Process: s u > c s h WEIGHTED=6
269C * Decay: h > ta+ ta- WEIGHTED=2
270C
271 IMPLICIT NONE
272C
273C CONSTANTS
274C
275 INTEGER NEXTERNAL
276 parameter(nexternal=6)
277 INTEGER NCOMB
278 parameter( ncomb=64)
279C
280C ARGUMENTS
281C
282 real*8 p(0:3,nexternal),ans
283 INTEGER H1,H2
284C
285C LOCAL VARIABLES
286C
287 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
288 real*8 t
289 real*8 matrix_su_cs_h
290 INTEGER IHEL,IDEN, I
291 INTEGER JC(NEXTERNAL)
292 LOGICAL GOODHEL(NCOMB)
293 DATA ntry/0/
294 DATA goodhel/ncomb*.false./
295 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
296 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
297 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
298 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
299 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
300 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
301 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
302 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
303 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
304 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
305 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
306 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
307 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
308 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
309 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
310 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
311 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
312 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
313 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
314 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
315 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
316 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
317 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
318 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
319 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
320 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
321 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
322 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
323 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
324 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
325 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
326 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
327 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
328 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
329 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
330 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
331 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
332 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
333 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
334 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
335 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
336 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
337 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
338 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
339 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
340 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
341 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
342 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
343 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
344 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
345 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
346 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
347 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
348 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
349 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
350 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
351 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
352 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
353 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
354 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
355 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
356 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
357 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
358 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
359 DATA iden/36/
360C ----------
361C BEGIN CODE
362C ----------
363 DO ihel=1,nexternal
364 jc(ihel) = +1
365 ENDDO
366 ans = 0d0
367 DO ihel=1,ncomb
368 t=matrix_su_cs_h(p ,h1,h2,nhel(1,ihel),jc(1))
369 ans=ans+t
370 ENDDO
371 ans=ans/dble(iden)
372 END
373
374
375 REAL*8 FUNCTION matrix_su_cs_h(P,H1,H2,NHEL,IC)
376C
377C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
378C By the MadGraph Development Team
379C Please visit us at https://launchpad.net/madgraph5
380C
381C Returns amplitude squared summed/avg over colors
382C for the point with external lines W(0:6,NEXTERNAL)
383C
384C Process: s u > c s h WEIGHTED=6
385C * Decay: h > ta+ ta- WEIGHTED=2
386C
387 IMPLICIT NONE
388C
389C CONSTANTS
390C
391 INTEGER NGRAPHS
392 parameter(ngraphs=1)
393 INTEGER NEXTERNAL
394 parameter(nexternal=6)
395 INTEGER NWAVEFUNCS, NCOLOR
396 parameter(nwavefuncs=7, ncolor=1)
397 real*8 zero
398 parameter(zero=0d0)
399 COMPLEX*16 IMAG1
400 parameter(imag1=(0d0,1d0))
401C
402C ARGUMENTS
403C
404 real*8 p(0:3,nexternal)
405 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
406C
407C LOCAL VARIABLES
408C
409 INTEGER I,J
410 COMPLEX*16 ZTEMP
411 real*8 denom(ncolor), cf(ncolor,ncolor)
412 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
413 COMPLEX*16 W(18,NWAVEFUNCS)
414 COMPLEX*16 DUM0,DUM1
415 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
416C
417C GLOBAL VARIABLES
418C
419 include 'coupl.inc'
420C
421C COLOR DATA
422C
423 DATA denom(1)/1/
424 DATA (cf(i, 1),i= 1, 1) / 9/
425C 1 T(3,1) T(4,2)
426
427 INTEGER H1,H2
428 real*8 matrix
429 matrix_su_cs_h=0.d0
430 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
431 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
432
433C ----------
434C BEGIN CODE
435C ----------
436 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
437 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
438 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
439 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
440 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
441 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
442 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
443 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
444 CALL ffv2_3(w(1,2),w(1,4),gc_101,mw,ww,w(1,3))
445C Amplitude(s) for diagram number 1
446 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
447 jamp(1)=-amp(1)
448
449 matrix = 0.d0
450 DO i = 1, ncolor
451 ztemp = (0.d0,0.d0)
452 DO j = 1, ncolor
453 ztemp = ztemp + cf(j,i)*jamp(j)
454 ENDDO
455 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
456 ENDDO
457 matrix_su_cs_h=matrix
458 ENDIF
459 ENDIF
460 END
461
462
463
464C===============================================
465
466 SUBROUTINE su_ud_h(P,H1,H2,ANS)
467C
468C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
469C By the MadGraph Development Team
470C Please visit us at https://launchpad.net/madgraph5
471C
472C MadGraph StandAlone Version
473C
474C Returns amplitude squared summed/avg over colors
475C and helicities
476C for the point in phase space P(0:3,NEXTERNAL)
477C
478C Process: s u > u d h WEIGHTED=6
479C * Decay: h > ta+ ta- WEIGHTED=2
480C
481 IMPLICIT NONE
482C
483C CONSTANTS
484C
485 INTEGER NEXTERNAL
486 parameter(nexternal=6)
487 INTEGER NCOMB
488 parameter( ncomb=64)
489C
490C ARGUMENTS
491C
492 real*8 p(0:3,nexternal),ans
493 INTEGER H1,H2
494C
495C LOCAL VARIABLES
496C
497 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
498 real*8 t
499 real*8 matrix_su_ud_h
500 INTEGER IHEL,IDEN, I
501 INTEGER JC(NEXTERNAL)
502 LOGICAL GOODHEL(NCOMB)
503 DATA ntry/0/
504 DATA goodhel/ncomb*.false./
505 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
506 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
507 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
508 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
509 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
510 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
511 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
512 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
513 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
514 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
515 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
516 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
517 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
518 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
519 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
520 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
521 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
522 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
523 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
524 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
525 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
526 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
527 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
528 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
529 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
530 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
531 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
532 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
533 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
534 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
535 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
536 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
537 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
538 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
539 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
540 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
541 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
542 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
543 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
544 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
545 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
546 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
547 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
548 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
549 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
550 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
551 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
552 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
553 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
554 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
555 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
556 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
557 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
558 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
559 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
560 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
561 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
562 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
563 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
564 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
565 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
566 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
567 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
568 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
569 DATA iden/36/
570C ----------
571C BEGIN CODE
572C ----------
573 DO ihel=1,nexternal
574 jc(ihel) = +1
575 ENDDO
576 ans = 0d0
577 DO ihel=1,ncomb
578 t=matrix_su_ud_h(p ,h1,h2,nhel(1,ihel),jc(1))
579 ans=ans+t
580 ENDDO
581 ans=ans/dble(iden)
582 END
583
584
585 REAL*8 FUNCTION matrix_su_ud_h(P,H1,H2,NHEL,IC)
586C
587C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
588C By the MadGraph Development Team
589C Please visit us at https://launchpad.net/madgraph5
590C
591C Returns amplitude squared summed/avg over colors
592C for the point with external lines W(0:6,NEXTERNAL)
593C
594C Process: s u > u d h WEIGHTED=6
595C * Decay: h > ta+ ta- WEIGHTED=2
596C
597 IMPLICIT NONE
598C
599C CONSTANTS
600C
601 INTEGER NGRAPHS
602 parameter(ngraphs=1)
603 INTEGER NEXTERNAL
604 parameter(nexternal=6)
605 INTEGER NWAVEFUNCS, NCOLOR
606 parameter(nwavefuncs=7, ncolor=1)
607 real*8 zero
608 parameter(zero=0d0)
609 COMPLEX*16 IMAG1
610 parameter(imag1=(0d0,1d0))
611C
612C ARGUMENTS
613C
614 real*8 p(0:3,nexternal)
615 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
616C
617C LOCAL VARIABLES
618C
619 INTEGER I,J
620 COMPLEX*16 ZTEMP
621 real*8 denom(ncolor), cf(ncolor,ncolor)
622 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
623 COMPLEX*16 W(18,NWAVEFUNCS)
624 COMPLEX*16 DUM0,DUM1
625 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
626C
627C GLOBAL VARIABLES
628C
629 include 'coupl.inc'
630C
631C COLOR DATA
632C
633 DATA denom(1)/1/
634 DATA (cf(i, 1),i= 1, 1) / 9/
635C 1 T(3,1) T(4,2)
636
637 INTEGER H1,H2
638 real*8 matrix
639 matrix_su_ud_h=0.d0
640 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
641 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
642
643C ----------
644C BEGIN CODE
645C ----------
646 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
647 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
648 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
649 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
650 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
651 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
652 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
653 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,6))
654 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,3))
655C Amplitude(s) for diagram number 1
656 CALL vvs1_0(w(1,3),w(1,6),w(1,7),gc_72,amp(1))
657 jamp(1)=-amp(1)
658
659 matrix = 0.d0
660 DO i = 1, ncolor
661 ztemp = (0.d0,0.d0)
662 DO j = 1, ncolor
663 ztemp = ztemp + cf(j,i)*jamp(j)
664 ENDDO
665 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
666 ENDDO
667 matrix_su_ud_h=matrix
668 ENDIF
669 ENDIF
670 END
671
672
673
674C========================================
675
676 SUBROUTINE su_us_h(P,H1,H2,ANS)
677C
678C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
679C By the MadGraph Development Team
680C Please visit us at https://launchpad.net/madgraph5
681C
682C MadGraph StandAlone Version
683C
684C Returns amplitude squared summed/avg over colors
685C and helicities
686C for the point in phase space P(0:3,NEXTERNAL)
687C
688C Process: s u > u s h WEIGHTED=6
689C * Decay: h > ta+ ta- WEIGHTED=2
690C
691 IMPLICIT NONE
692C
693C CONSTANTS
694C
695 INTEGER NEXTERNAL
696 parameter(nexternal=6)
697 INTEGER NCOMB
698 parameter( ncomb=64)
699C
700C ARGUMENTS
701C
702 real*8 p(0:3,nexternal),ans
703 INTEGER H1,H2
704C
705C LOCAL VARIABLES
706C
707 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
708 real*8 t
709 real*8 matrix_su_us_h
710 INTEGER IHEL,IDEN, I
711 INTEGER JC(NEXTERNAL)
712 LOGICAL GOODHEL(NCOMB)
713 DATA ntry/0/
714 DATA goodhel/ncomb*.false./
715 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
716 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
717 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
718 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
719 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
720 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
721 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
722 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
723 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
724 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
725 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
726 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
727 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
728 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
729 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
730 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
731 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
732 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
733 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
734 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
735 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
736 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
737 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
738 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
739 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
740 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
741 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
742 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
743 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
744 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
745 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
746 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
747 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
748 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
749 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
750 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
751 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
752 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
753 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
754 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
755 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
756 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
757 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
758 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
759 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
760 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
761 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
762 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
763 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
764 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
765 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
766 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
767 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
768 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
769 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
770 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
771 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
772 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
773 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
774 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
775 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
776 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
777 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
778 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
779 DATA iden/36/
780C ----------
781C BEGIN CODE
782C ----------
783 DO ihel=1,nexternal
784 jc(ihel) = +1
785 ENDDO
786 ans = 0d0
787 DO ihel=1,ncomb
788 t=matrix_su_us_h(p ,h1,h2,nhel(1,ihel),jc(1))
789 ans=ans+t
790 ENDDO
791 ans=ans/dble(iden)
792 END
793
794
795 REAL*8 FUNCTION matrix_su_us_h(P,H1,H2,NHEL,IC)
796C
797C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
798C By the MadGraph Development Team
799C Please visit us at https://launchpad.net/madgraph5
800C
801C Returns amplitude squared summed/avg over colors
802C for the point with external lines W(0:6,NEXTERNAL)
803C
804C Process: s u > u s h WEIGHTED=6
805C * Decay: h > ta+ ta- WEIGHTED=2
806C
807 IMPLICIT NONE
808C
809C CONSTANTS
810C
811 INTEGER NGRAPHS
812 parameter(ngraphs=2)
813 INTEGER NEXTERNAL
814 parameter(nexternal=6)
815 INTEGER NWAVEFUNCS, NCOLOR
816 parameter(nwavefuncs=7, ncolor=2)
817 real*8 zero
818 parameter(zero=0d0)
819 COMPLEX*16 IMAG1
820 parameter(imag1=(0d0,1d0))
821C
822C ARGUMENTS
823C
824 real*8 p(0:3,nexternal)
825 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
826C
827C LOCAL VARIABLES
828C
829 INTEGER I,J
830 COMPLEX*16 ZTEMP
831 real*8 denom(ncolor), cf(ncolor,ncolor)
832 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
833 COMPLEX*16 W(18,NWAVEFUNCS)
834 COMPLEX*16 DUM0,DUM1
835 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
836C
837C GLOBAL VARIABLES
838C
839 include 'coupl.inc'
840C
841C COLOR DATA
842C
843 DATA denom(1)/1/
844 DATA (cf(i, 1),i= 1, 2) / 9, 3/
845C 1 T(3,1) T(4,2)
846 DATA denom(2)/1/
847 DATA (cf(i, 2),i= 1, 2) / 3, 9/
848C 1 T(3,2) T(4,1)
849
850 INTEGER H1,H2
851 real*8 matrix
852 matrix_su_us_h=0.d0
853 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
854 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
855
856C ----------
857C BEGIN CODE
858C ----------
859 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
860 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
861 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
862 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
863 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
864 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
865 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
866 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,6))
867 CALL ffv2_3(w(1,2),w(1,4),gc_101,mw,ww,w(1,5))
868C Amplitude(s) for diagram number 1
869 CALL vvs1_0(w(1,5),w(1,6),w(1,7),gc_72,amp(1))
870 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
871 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,4))
872C Amplitude(s) for diagram number 2
873 CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_81,amp(2))
874 jamp(1)=-amp(1)
875 jamp(2)=+amp(2)
876
877 matrix = 0.d0
878 DO i = 1, ncolor
879 ztemp = (0.d0,0.d0)
880 DO j = 1, ncolor
881 ztemp = ztemp + cf(j,i)*jamp(j)
882 ENDDO
883 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
884 ENDDO
885 matrix_su_us_h=matrix
886 ENDIF
887 ENDIF
888 END
889
890
891
892
893C --- end of subprocess SU->jjH, H->tautau
894
895C ----- begin subprocesses SU->jjtautau, no Higgs
896
897 SUBROUTINE su_cd_noh(P,H1,H2,ANS)
898C
899C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
900C By the MadGraph Development Team
901C Please visit us at https://launchpad.net/madgraph5
902C
903C MadGraph StandAlone Version
904C
905C Returns amplitude squared summed/avg over colors
906C and helicities
907C for the point in phase space P(0:3,NEXTERNAL)
908C
909C Process: s u > c d ta+ ta- / h QED=4
910C
911 IMPLICIT NONE
912C
913C CONSTANTS
914C
915 INTEGER NEXTERNAL
916 parameter(nexternal=6)
917 INTEGER NCOMB
918 parameter( ncomb=64)
919C
920C ARGUMENTS
921C
922 real*8 p(0:3,nexternal),ans
923 INTEGER H1,H2
924C
925C LOCAL VARIABLES
926C
927 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
928 real*8 t
929 real*8 matrix_su_cd_noh
930 INTEGER IHEL,IDEN, I
931 INTEGER JC(NEXTERNAL)
932 LOGICAL GOODHEL(NCOMB)
933 DATA ntry/0/
934 DATA goodhel/ncomb*.false./
935 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
936 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
937 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
938 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
939 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
940 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
941 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
942 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
943 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
944 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
945 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
946 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
947 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
948 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
949 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
950 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
951 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
952 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
953 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
954 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
955 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
956 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
957 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
958 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
959 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
960 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
961 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
962 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
963 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
964 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
965 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
966 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
967 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
968 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
969 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
970 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
971 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
972 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
973 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
974 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
975 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
976 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
977 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
978 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
979 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
980 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
981 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
982 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
983 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
984 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
985 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
986 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
987 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
988 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
989 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
990 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
991 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
992 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
993 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
994 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
995 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
996 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
997 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
998 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
999 DATA iden/36/
1000C ----------
1001C BEGIN CODE
1002C ----------
1003 DO ihel=1,nexternal
1004 jc(ihel) = +1
1005 ENDDO
1006 ans = 0d0
1007 DO ihel=1,ncomb
1008 t=matrix_su_cd_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1009 ans=ans+t
1010 ENDDO
1011 ans=ans/dble(iden)
1012 END
1013
1014
1015 REAL*8 FUNCTION matrix_su_cd_noh(P,H1,H2,NHEL,IC)
1016C
1017C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1018C By the MadGraph Development Team
1019C Please visit us at https://launchpad.net/madgraph5
1020C
1021C Returns amplitude squared summed/avg over colors
1022C for the point with external lines W(0:6,NEXTERNAL)
1023C
1024C Process: s u > c d ta+ ta- / h QED=4
1025C
1026 IMPLICIT NONE
1027C
1028C CONSTANTS
1029C
1030 INTEGER NGRAPHS
1031 parameter(ngraphs=11)
1032 INTEGER NEXTERNAL
1033 parameter(nexternal=6)
1034 INTEGER NWAVEFUNCS, NCOLOR
1035 parameter(nwavefuncs=10, ncolor=1)
1036 real*8 zero
1037 parameter(zero=0d0)
1038 COMPLEX*16 IMAG1
1039 parameter(imag1=(0d0,1d0))
1040C
1041C ARGUMENTS
1042C
1043 real*8 p(0:3,nexternal)
1044 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1045C
1046C LOCAL VARIABLES
1047C
1048 INTEGER I,J
1049 COMPLEX*16 ZTEMP
1050 real*8 denom(ncolor), cf(ncolor,ncolor)
1051 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1052 COMPLEX*16 W(18,NWAVEFUNCS)
1053 COMPLEX*16 DUM0,DUM1
1054 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1055C
1056C GLOBAL VARIABLES
1057C
1058 include 'coupl.inc'
1059C
1060C COLOR DATA
1061C
1062 DATA denom(1)/1/
1063 DATA (cf(i, 1),i= 1, 1) / 9/
1064C 1 T(3,1) T(4,2)
1065
1066 INTEGER H1,H2
1067 real*8 matrix
1068 matrix_su_cd_noh=0.d0
1069 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1070 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1071
1072C ----------
1073C BEGIN CODE
1074C ----------
1075 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1076 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1077 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1078 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1079 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1080 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1081 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
1082 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,8))
1083 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1084C Amplitude(s) for diagram number 1
1085 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1086 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1087C Amplitude(s) for diagram number 2
1088 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1089 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1090C Amplitude(s) for diagram number 3
1091 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1092 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1093C Amplitude(s) for diagram number 4
1094 CALL ffv1_0(w(1,6),w(1,4),w(1,9),gc_1,amp(4))
1095 CALL ffv2_1(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1096C Amplitude(s) for diagram number 5
1097 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(5))
1098C Amplitude(s) for diagram number 6
1099 CALL ffv2_3_0(w(1,6),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1100C Amplitude(s) for diagram number 7
1101 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(7))
1102 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1103C Amplitude(s) for diagram number 8
1104 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1105 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1106C Amplitude(s) for diagram number 9
1107 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
1108C Amplitude(s) for diagram number 10
1109 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1110 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1111C Amplitude(s) for diagram number 11
1112 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
1113 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1114 $ -amp(9)-amp(10)-amp(11)
1115
1116 matrix = 0.d0
1117 DO i = 1, ncolor
1118 ztemp = (0.d0,0.d0)
1119 DO j = 1, ncolor
1120 ztemp = ztemp + cf(j,i)*jamp(j)
1121 ENDDO
1122 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1123 ENDDO
1124 matrix_su_cd_noh=matrix
1125 ENDIF ! CLOSES H1 IF
1126 ENDIF ! CLOSES H2 IF
1127 END
1128
1129
1130 SUBROUTINE su_cs_noh(P,H1,H2,ANS)
1131C
1132C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1133C By the MadGraph Development Team
1134C Please visit us at https://launchpad.net/madgraph5
1135C
1136C MadGraph StandAlone Version
1137C
1138C Returns amplitude squared summed/avg over colors
1139C and helicities
1140C for the point in phase space P(0:3,NEXTERNAL)
1141C
1142C Process: s u > c s ta+ ta- / h QED=4
1143C
1144 IMPLICIT NONE
1145C
1146C CONSTANTS
1147C
1148 INTEGER NEXTERNAL
1149 parameter(nexternal=6)
1150 INTEGER NCOMB
1151 parameter( ncomb=64)
1152C
1153C ARGUMENTS
1154C
1155 real*8 p(0:3,nexternal),ans
1156 INTEGER H1,H2
1157C
1158C LOCAL VARIABLES
1159C
1160 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1161 real*8 t
1162 real*8 matrix_su_cs_noh
1163 INTEGER IHEL,IDEN, I
1164 INTEGER JC(NEXTERNAL)
1165 LOGICAL GOODHEL(NCOMB)
1166 DATA ntry/0/
1167 DATA goodhel/ncomb*.false./
1168 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1169 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1170 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1171 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1172 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1173 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1174 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1175 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1176 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1177 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1178 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1179 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1180 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1181 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1182 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1183 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1184 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1185 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1186 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1187 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1188 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1189 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1190 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1191 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1192 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1193 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1194 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1195 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1196 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1197 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1198 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1199 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1200 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1201 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1202 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1203 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1204 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1205 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1206 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1207 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1208 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1209 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1210 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1211 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1212 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1213 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1214 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1215 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1216 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1217 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1218 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1219 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1220 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1221 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1222 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1223 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1224 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1225 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1226 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1227 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1228 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1229 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1230 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1231 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1232 DATA iden/36/
1233C ----------
1234C BEGIN CODE
1235C ----------
1236 DO ihel=1,nexternal
1237 jc(ihel) = +1
1238 ENDDO
1239 ans = 0d0
1240 DO ihel=1,ncomb
1241 t=matrix_su_cs_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1242 ans=ans+t
1243 ENDDO
1244 ans=ans/dble(iden)
1245 END
1246
1247
1248 REAL*8 FUNCTION matrix_su_cs_noh(P,H1,H2,NHEL,IC)
1249C
1250C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1251C By the MadGraph Development Team
1252C Please visit us at https://launchpad.net/madgraph5
1253C
1254C Returns amplitude squared summed/avg over colors
1255C for the point with external lines W(0:6,NEXTERNAL)
1256C
1257C Process: s u > c s ta+ ta- / h QED=4
1258C
1259 IMPLICIT NONE
1260C
1261C CONSTANTS
1262C
1263 INTEGER NGRAPHS
1264 parameter(ngraphs=11)
1265 INTEGER NEXTERNAL
1266 parameter(nexternal=6)
1267 INTEGER NWAVEFUNCS, NCOLOR
1268 parameter(nwavefuncs=10, ncolor=1)
1269 real*8 zero
1270 parameter(zero=0d0)
1271 COMPLEX*16 IMAG1
1272 parameter(imag1=(0d0,1d0))
1273C
1274C ARGUMENTS
1275C
1276 real*8 p(0:3,nexternal)
1277 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1278C
1279C LOCAL VARIABLES
1280C
1281 INTEGER I,J
1282 COMPLEX*16 ZTEMP
1283 real*8 denom(ncolor), cf(ncolor,ncolor)
1284 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1285 COMPLEX*16 W(18,NWAVEFUNCS)
1286 COMPLEX*16 DUM0,DUM1
1287 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1288C
1289C GLOBAL VARIABLES
1290C
1291 include 'coupl.inc'
1292C
1293C COLOR DATA
1294C
1295 DATA denom(1)/1/
1296 DATA (cf(i, 1),i= 1, 1) / 9/
1297C 1 T(3,1) T(4,2)
1298
1299 INTEGER H1,H2
1300 real*8 matrix
1301 matrix_su_cs_noh=0.d0
1302 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1303 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1304
1305C ----------
1306C BEGIN CODE
1307C ----------
1308 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1309 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1310 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1311 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1312 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1313 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1314 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
1315 CALL ffv2_3(w(1,2),w(1,4),gc_101,mw,ww,w(1,8))
1316 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1317C Amplitude(s) for diagram number 1
1318 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1319 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1320C Amplitude(s) for diagram number 2
1321 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1322 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1323C Amplitude(s) for diagram number 3
1324 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1325 CALL ffv2_2(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
1326C Amplitude(s) for diagram number 4
1327 CALL ffv1_0(w(1,6),w(1,4),w(1,9),gc_1,amp(4))
1328 CALL ffv2_1(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
1329C Amplitude(s) for diagram number 5
1330 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(5))
1331C Amplitude(s) for diagram number 6
1332 CALL ffv2_3_0(w(1,6),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1333C Amplitude(s) for diagram number 7
1334 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(7))
1335 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1336C Amplitude(s) for diagram number 8
1337 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1338 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1339C Amplitude(s) for diagram number 9
1340 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
1341C Amplitude(s) for diagram number 10
1342 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1343 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1344C Amplitude(s) for diagram number 11
1345 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
1346 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1347 $ -amp(9)-amp(10)-amp(11)
1348
1349 matrix = 0.d0
1350 DO i = 1, ncolor
1351 ztemp = (0.d0,0.d0)
1352 DO j = 1, ncolor
1353 ztemp = ztemp + cf(j,i)*jamp(j)
1354 ENDDO
1355 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1356 ENDDO
1357 matrix_su_cs_noh=matrix
1358 ENDIF ! CLOSES H1 IF
1359 ENDIF ! CLOSES H2 IF
1360 END
1361
1362 SUBROUTINE su_ud_noh(P,H1,H2,ANS)
1363C
1364C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1365C By the MadGraph Development Team
1366C Please visit us at https://launchpad.net/madgraph5
1367C
1368C MadGraph StandAlone Version
1369C
1370C Returns amplitude squared summed/avg over colors
1371C and helicities
1372C for the point in phase space P(0:3,NEXTERNAL)
1373C
1374C Process: s u > u d ta+ ta- / h QED=4
1375C
1376 IMPLICIT NONE
1377C
1378C CONSTANTS
1379C
1380 INTEGER NEXTERNAL
1381 parameter(nexternal=6)
1382 INTEGER NCOMB
1383 parameter( ncomb=64)
1384C
1385C ARGUMENTS
1386C
1387 real*8 p(0:3,nexternal),ans
1388 INTEGER H1,H2
1389C
1390C LOCAL VARIABLES
1391C
1392 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1393 real*8 t
1394 real*8 matrix_su_ud_noh
1395 INTEGER IHEL,IDEN, I
1396 INTEGER JC(NEXTERNAL)
1397 LOGICAL GOODHEL(NCOMB)
1398 DATA ntry/0/
1399 DATA goodhel/ncomb*.false./
1400 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1401 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1402 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1403 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1404 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1405 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1406 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1407 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1408 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1409 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1410 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1411 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1412 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1413 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1414 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1415 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1416 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1417 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1418 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1419 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1420 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1421 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1422 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1423 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1424 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1425 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1426 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1427 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1428 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1429 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1430 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1431 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1432 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1433 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1434 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1435 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1436 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1437 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1438 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1439 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1440 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1441 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1442 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1443 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1444 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1445 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1446 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1447 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1448 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1449 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1450 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1451 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1452 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1453 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1454 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1455 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1456 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1457 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1458 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1459 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1460 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1461 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1462 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1463 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1464 DATA iden/36/
1465C ----------
1466C BEGIN CODE
1467C ----------
1468 DO ihel=1,nexternal
1469 jc(ihel) = +1
1470 ENDDO
1471 ans = 0d0
1472 DO ihel=1,ncomb
1473 t=matrix_su_ud_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1474 ans=ans+t
1475 ENDDO
1476 ans=ans/dble(iden)
1477 END
1478
1479
1480 REAL*8 FUNCTION matrix_su_ud_noh(P,H1,H2,NHEL,IC)
1481C
1482C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1483C By the MadGraph Development Team
1484C Please visit us at https://launchpad.net/madgraph5
1485C
1486C Returns amplitude squared summed/avg over colors
1487C for the point with external lines W(0:6,NEXTERNAL)
1488C
1489C Process: s u > u d ta+ ta- / h QED=4
1490C
1491 IMPLICIT NONE
1492C
1493C CONSTANTS
1494C
1495 INTEGER NGRAPHS
1496 parameter(ngraphs=11)
1497 INTEGER NEXTERNAL
1498 parameter(nexternal=6)
1499 INTEGER NWAVEFUNCS, NCOLOR
1500 parameter(nwavefuncs=10, ncolor=1)
1501 real*8 zero
1502 parameter(zero=0d0)
1503 COMPLEX*16 IMAG1
1504 parameter(imag1=(0d0,1d0))
1505C
1506C ARGUMENTS
1507C
1508 real*8 p(0:3,nexternal)
1509 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1510C
1511C LOCAL VARIABLES
1512C
1513 INTEGER I,J
1514 COMPLEX*16 ZTEMP
1515 real*8 denom(ncolor), cf(ncolor,ncolor)
1516 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1517 COMPLEX*16 W(18,NWAVEFUNCS)
1518 COMPLEX*16 DUM0,DUM1
1519 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1520C
1521C GLOBAL VARIABLES
1522C
1523 include 'coupl.inc'
1524C
1525C COLOR DATA
1526C
1527 DATA denom(1)/1/
1528 DATA (cf(i, 1),i= 1, 1) / 9/
1529C 1 T(3,1) T(4,2)
1530
1531 INTEGER H1,H2
1532 real*8 matrix
1533 matrix_su_ud_noh=0.d0
1534 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1535 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1536
1537C ----------
1538C BEGIN CODE
1539C ----------
1540 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1541 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1542 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1543 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1544 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1545 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1546 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,7))
1547 CALL ffv2_3(w(1,2),w(1,4),gc_100,mw,ww,w(1,8))
1548 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1549C Amplitude(s) for diagram number 1
1550 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1551 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1552C Amplitude(s) for diagram number 2
1553 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1554 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1555C Amplitude(s) for diagram number 3
1556 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1557 CALL ffv2_2(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
1558C Amplitude(s) for diagram number 4
1559 CALL ffv1_0(w(1,6),w(1,4),w(1,9),gc_1,amp(4))
1560 CALL ffv2_1(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1561C Amplitude(s) for diagram number 5
1562 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(5))
1563C Amplitude(s) for diagram number 6
1564 CALL ffv2_3_0(w(1,6),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1565C Amplitude(s) for diagram number 7
1566 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(7))
1567 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1568C Amplitude(s) for diagram number 8
1569 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(8))
1570 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,2))
1571C Amplitude(s) for diagram number 9
1572 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_101,amp(9))
1573C Amplitude(s) for diagram number 10
1574 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(10))
1575 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1576C Amplitude(s) for diagram number 11
1577 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_101,amp(11))
1578 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1579 $ -amp(9)-amp(10)-amp(11)
1580
1581 matrix = 0.d0
1582 DO i = 1, ncolor
1583 ztemp = (0.d0,0.d0)
1584 DO j = 1, ncolor
1585 ztemp = ztemp + cf(j,i)*jamp(j)
1586 ENDDO
1587 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1588 ENDDO
1589 matrix_su_ud_noh=matrix
1590 ENDIF ! CLOSES H1 IF
1591 ENDIF ! CLOSES H2 IF
1592 END
1593
1594
1595 SUBROUTINE su_us_noh(P,H1,H2,ANS)
1596C
1597C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1598C By the MadGraph Development Team
1599C Please visit us at https://launchpad.net/madgraph5
1600C
1601C MadGraph StandAlone Version
1602C
1603C Returns amplitude squared summed/avg over colors
1604C and helicities
1605C for the point in phase space P(0:3,NEXTERNAL)
1606C
1607C Process: s u > u s ta+ ta- / h QED=4
1608C
1609 IMPLICIT NONE
1610C
1611C CONSTANTS
1612C
1613 INTEGER NEXTERNAL
1614 parameter(nexternal=6)
1615 INTEGER NCOMB
1616 parameter( ncomb=64)
1617C
1618C ARGUMENTS
1619C
1620 real*8 p(0:3,nexternal),ans
1621 INTEGER H1,H2
1622C
1623C LOCAL VARIABLES
1624C
1625 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1626 real*8 t
1627 real*8 matrix_su_us_noh
1628 INTEGER IHEL,IDEN, I
1629 INTEGER JC(NEXTERNAL)
1630 LOGICAL GOODHEL(NCOMB)
1631 DATA ntry/0/
1632 DATA goodhel/ncomb*.false./
1633 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1634 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1635 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1636 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1637 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1638 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1639 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1640 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1641 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1642 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1643 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1644 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1645 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1646 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1647 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1648 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1649 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1650 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1651 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1652 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1653 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1654 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1655 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1656 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1657 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1658 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1659 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1660 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1661 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1662 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1663 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1664 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1665 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1666 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1667 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1668 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1669 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1670 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1671 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1672 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1673 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1674 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1675 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1676 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1677 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1678 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1679 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1680 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1681 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1682 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1683 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1684 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1685 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1686 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1687 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1688 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1689 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1690 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1691 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1692 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1693 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1694 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1695 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1696 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1697 DATA iden/36/
1698C ----------
1699C BEGIN CODE
1700C ----------
1701 DO ihel=1,nexternal
1702 jc(ihel) = +1
1703 ENDDO
1704 ans = 0d0
1705 DO ihel=1,ncomb
1706 t=matrix_su_us_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1707 ans=ans+t
1708 ENDDO
1709 ans=ans/dble(iden)
1710 END
1711
1712
1713 REAL*8 FUNCTION matrix_su_us_noh(P,H1,H2,NHEL,IC)
1714C
1715C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1716C By the MadGraph Development Team
1717C Please visit us at https://launchpad.net/madgraph5
1718C
1719C Returns amplitude squared summed/avg over colors
1720C for the point with external lines W(0:6,NEXTERNAL)
1721C
1722C Process: s u > u s ta+ ta- / h QED=4
1723C
1724 IMPLICIT NONE
1725C
1726C CONSTANTS
1727C
1728 INTEGER NGRAPHS
1729 parameter(ngraphs=43)
1730 INTEGER NEXTERNAL
1731 parameter(nexternal=6)
1732 INTEGER NWAVEFUNCS, NCOLOR
1733 parameter(nwavefuncs=15, ncolor=2)
1734 real*8 zero
1735 parameter(zero=0d0)
1736 COMPLEX*16 IMAG1
1737 parameter(imag1=(0d0,1d0))
1738C
1739C ARGUMENTS
1740C
1741 real*8 p(0:3,nexternal)
1742 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1743C
1744C LOCAL VARIABLES
1745C
1746 INTEGER I,J
1747 COMPLEX*16 ZTEMP
1748 real*8 denom(ncolor), cf(ncolor,ncolor)
1749 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1750 COMPLEX*16 W(18,NWAVEFUNCS)
1751 COMPLEX*16 DUM0,DUM1
1752 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1753C
1754C GLOBAL VARIABLES
1755C
1756 include 'coupl.inc'
1757C
1758C COLOR DATA
1759C
1760 DATA denom(1)/1/
1761 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1762C 1 T(3,1) T(4,2)
1763 DATA denom(2)/1/
1764 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1765C 1 T(3,2) T(4,1)
1766
1767 INTEGER H1,H2
1768 real*8 matrix
1769 matrix_su_us_noh=0.d0
1770 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1771 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1772
1773
1774C ----------
1775C BEGIN CODE
1776C ----------
1777 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1778 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
1779 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1780 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
1781 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1782 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1783 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,7))
1784 CALL ffv2_3(w(1,2),w(1,4),gc_101,mw,ww,w(1,8))
1785 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1786C Amplitude(s) for diagram number 1
1787 CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1788 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1789C Amplitude(s) for diagram number 2
1790 CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1791 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1792C Amplitude(s) for diagram number 3
1793 CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1794 CALL ffv2_2(w(1,2),w(1,7),gc_101,zero,zero,w(1,11))
1795C Amplitude(s) for diagram number 4
1796 CALL ffv1_0(w(1,11),w(1,4),w(1,9),gc_1,amp(4))
1797 CALL ffv2_1(w(1,4),w(1,7),gc_101,zero,zero,w(1,12))
1798C Amplitude(s) for diagram number 5
1799 CALL ffv1_0(w(1,2),w(1,12),w(1,9),gc_2,amp(5))
1800C Amplitude(s) for diagram number 6
1801 CALL ffv2_3_0(w(1,11),w(1,4),w(1,10),gc_50,gc_58,amp(6))
1802C Amplitude(s) for diagram number 7
1803 CALL ffv2_5_0(w(1,2),w(1,12),w(1,10),gc_51,gc_58,amp(7))
1804 CALL ffv1p0_3(w(1,1),w(1,4),gc_1,zero,zero,w(1,12))
1805 CALL ffv1p0_3(w(1,2),w(1,3),gc_2,zero,zero,w(1,11))
1806 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1807C Amplitude(s) for diagram number 8
1808 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1809 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1810C Amplitude(s) for diagram number 9
1811 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1812 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,14))
1813C Amplitude(s) for diagram number 10
1814 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1815C Amplitude(s) for diagram number 11
1816 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1817 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,13))
1818 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1819C Amplitude(s) for diagram number 12
1820 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1821 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1822C Amplitude(s) for diagram number 13
1823 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1824C Amplitude(s) for diagram number 14
1825 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1826C Amplitude(s) for diagram number 15
1827 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1828 CALL ffv1_2(w(1,2),w(1,12),gc_2,zero,zero,w(1,15))
1829C Amplitude(s) for diagram number 16
1830 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(16))
1831 CALL ffv1_1(w(1,3),w(1,12),gc_2,zero,zero,w(1,5))
1832C Amplitude(s) for diagram number 17
1833 CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_2,amp(17))
1834C Amplitude(s) for diagram number 18
1835 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(18))
1836C Amplitude(s) for diagram number 19
1837 CALL ffv2_5_0(w(1,2),w(1,5),w(1,10),gc_51,gc_58,amp(19))
1838 CALL ffv1p0_3(w(1,1),w(1,4),gc_11,zero,zero,w(1,5))
1839 CALL ffv1_2(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1840C Amplitude(s) for diagram number 20
1841 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(20))
1842 CALL ffv1_1(w(1,3),w(1,5),gc_11,zero,zero,w(1,12))
1843C Amplitude(s) for diagram number 21
1844 CALL ffv1_0(w(1,2),w(1,12),w(1,9),gc_2,amp(21))
1845C Amplitude(s) for diagram number 22
1846 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(22))
1847C Amplitude(s) for diagram number 23
1848 CALL ffv2_5_0(w(1,2),w(1,12),w(1,10),gc_51,gc_58,amp(23))
1849 CALL ffv2_5_2(w(1,2),w(1,13),gc_51,gc_58,zero,zero,w(1,12))
1850C Amplitude(s) for diagram number 24
1851 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(24))
1852 CALL ffv2_5_1(w(1,3),w(1,13),gc_51,gc_58,zero,zero,w(1,15))
1853C Amplitude(s) for diagram number 25
1854 CALL ffv1_0(w(1,2),w(1,15),w(1,9),gc_2,amp(25))
1855C Amplitude(s) for diagram number 26
1856 CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(26))
1857C Amplitude(s) for diagram number 27
1858 CALL ffv2_5_0(w(1,2),w(1,15),w(1,10),gc_51,gc_58,amp(27))
1859 CALL ffv1_2(w(1,1),w(1,11),gc_1,zero,zero,w(1,15))
1860C Amplitude(s) for diagram number 28
1861 CALL ffv1_0(w(1,15),w(1,4),w(1,9),gc_1,amp(28))
1862 CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,12))
1863C Amplitude(s) for diagram number 29
1864 CALL ffv1_0(w(1,12),w(1,4),w(1,11),gc_1,amp(29))
1865C Amplitude(s) for diagram number 30
1866 CALL ffv2_3_0(w(1,15),w(1,4),w(1,10),gc_50,gc_58,amp(30))
1867 CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,15))
1868C Amplitude(s) for diagram number 31
1869 CALL ffv1_0(w(1,15),w(1,4),w(1,11),gc_1,amp(31))
1870 CALL ffv1p0_3(w(1,2),w(1,3),gc_11,zero,zero,w(1,11))
1871 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,2))
1872C Amplitude(s) for diagram number 32
1873 CALL ffv1_0(w(1,2),w(1,4),w(1,9),gc_1,amp(32))
1874C Amplitude(s) for diagram number 33
1875 CALL ffv1_0(w(1,12),w(1,4),w(1,11),gc_11,amp(33))
1876C Amplitude(s) for diagram number 34
1877 CALL ffv2_3_0(w(1,2),w(1,4),w(1,10),gc_50,gc_58,amp(34))
1878C Amplitude(s) for diagram number 35
1879 CALL ffv1_0(w(1,15),w(1,4),w(1,11),gc_11,amp(35))
1880 CALL ffv2_3_2(w(1,1),w(1,14),gc_50,gc_58,zero,zero,w(1,11))
1881C Amplitude(s) for diagram number 36
1882 CALL ffv1_0(w(1,11),w(1,4),w(1,9),gc_1,amp(36))
1883C Amplitude(s) for diagram number 37
1884 CALL ffv2_3_0(w(1,12),w(1,4),w(1,14),gc_50,gc_58,amp(37))
1885C Amplitude(s) for diagram number 38
1886 CALL ffv2_3_0(w(1,11),w(1,4),w(1,10),gc_50,gc_58,amp(38))
1887C Amplitude(s) for diagram number 39
1888 CALL ffv2_3_0(w(1,15),w(1,4),w(1,14),gc_50,gc_58,amp(39))
1889 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,14))
1890C Amplitude(s) for diagram number 40
1891 CALL ffv1_0(w(1,14),w(1,3),w(1,9),gc_2,amp(40))
1892C Amplitude(s) for diagram number 41
1893 CALL ffv2_0(w(1,12),w(1,3),w(1,8),gc_101,amp(41))
1894C Amplitude(s) for diagram number 42
1895 CALL ffv2_5_0(w(1,14),w(1,3),w(1,10),gc_51,gc_58,amp(42))
1896C Amplitude(s) for diagram number 43
1897 CALL ffv2_0(w(1,15),w(1,3),w(1,8),gc_101,amp(43))
1898 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1899 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1900 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1901 $ -amp(40)-amp(41)-amp(42)-amp(43)
1902 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1903 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1904 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1905 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1906 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1907 $ +amp(38)+amp(39)
1908
1909 matrix = 0.d0
1910 DO i = 1, ncolor
1911 ztemp = (0.d0,0.d0)
1912 DO j = 1, ncolor
1913 ztemp = ztemp + cf(j,i)*jamp(j)
1914 ENDDO
1915 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1916 ENDDO
1917 matrix_su_us_noh=matrix
1918 ENDIF ! CLOSES H1 IF
1919 ENDIF ! CLOSES H2 IF
1920 END
1921
1922
1923C ------- end subprocesses SU->jjtautau, no Higgs
1924C DONE 27.10.2015
1925