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