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