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