C++ Interface to Tauola
USX.f
1c --------------------- subprocesses initiated by USbar
2
3 SUBROUTINE usx(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
9C I3=0,2 OR 4 I4=0,-1 OR -3
10
11 ans=0.d0
12
13C FINAL CSX IS EQUAL TO FINAL UDX
14
15 IF(key.EQ.1) THEN
16 IF(i3.EQ.2 .AND. i4.EQ.-1) CALL usx_udx_h(p,h1,h2,ans)
17 IF(i3.EQ.2 .AND. i4.EQ.-3) CALL usx_usx_h(p,h1,h2,ans)
18 IF(i3.EQ.4 .AND. i4.EQ.-1) CAll usx_cdx_h(p,h1,h2,ans)
19 IF(i3.EQ.4 .AND. i4.EQ.-3) CALL usx_udx_h(p,h1,h2,ans)
20 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
21 CALL usx_udx_h(p,h1,h2,ans1)
22 CALL usx_usx_h(p,h1,h2,ans2)
23 CALL usx_cdx_h(p,h1,h2,ans3)
24 CALL usx_udx_h(p,h1,h2,ans4)
25 ans=ans1+ans2+ans3+ans4
26 ENDIF
27 ELSE IF(key.EQ.0) THEN
28 IF(i3.EQ.2 .AND. i4.EQ.-1) CALL usx_udx_noh(p,h1,h2,ans)
29 IF(i3.EQ.2 .AND. i4.EQ.-3) CALL usx_usx_noh(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.-1) CAll usx_cdx_noh(p,h1,h2,ans)
31 IF(i3.EQ.4 .AND. i4.EQ.-3) CALL usx_udx_noh(p,h1,h2,ans)
32 IF(i3.EQ.0 .AND. i4.EQ.0) THEN
33 CALL usx_udx_noh(p,h1,h2,ans1)
34 CALL usx_usx_noh(p,h1,h2,ans2)
35 CALL usx_cdx_noh(p,h1,h2,ans3)
36 CALL usx_udx_noh(p,h1,h2,ans4)
37 ans=ans1+ans2+ans3+ans4
38 ENDIF
39 ELSE
40 WRITE(*,*) 'NOT FINISHED'
41 stop
42 ENDIF
43 END ! SUBROUTINE USX
44
45
46
47C ----- begin subprocesses USX->UDX with Higgs->tautau
48
49
50 SUBROUTINE usx_udx_h(P,H1,H2,ANS)
51C
52C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
53C By the MadGraph Development Team
54C Please visit us at https://launchpad.net/madgraph5
55C
56C MadGraph StandAlone Version
57C
58C Returns amplitude squared summed/avg over colors
59C and helicities
60C for the point in phase space P(0:3,NEXTERNAL)
61C
62C Process: u s~ > u d~ h WEIGHTED=6
63C * Decay: h > ta+ ta- WEIGHTED=2
64C Process: u s~ > c s~ h WEIGHTED=6
65C * Decay: h > ta+ ta- WEIGHTED=2
66C
67 IMPLICIT NONE
68C
69C CONSTANTS
70C
71 INTEGER NEXTERNAL
72 parameter(nexternal=6)
73 INTEGER NCOMB
74 parameter( ncomb=64)
75C
76C ARGUMENTS
77C
78 real*8 p(0:3,nexternal),ans
79 INTEGER H1,H2
80C
81C LOCAL VARIABLES
82C
83 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
84 real*8 t
85 real*8 matrix_usx_udx_h
86 INTEGER IHEL,IDEN, I
87 INTEGER JC(NEXTERNAL)
88 LOGICAL GOODHEL(NCOMB)
89 DATA ntry/0/
90 DATA goodhel/ncomb*.false./
91 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
92 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
93 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
94 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
95 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
96 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
97 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
98 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
99 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
100 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
101 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
102 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
103 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
104 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
105 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
106 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
107 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
108 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
109 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
110 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
111 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
112 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
113 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
114 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
115 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
116 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
117 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
118 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
119 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
120 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
121 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
122 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
123 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
124 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
125 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
126 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
127 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
128 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
129 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
130 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
131 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
132 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
133 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
134 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
135 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
136 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
137 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
138 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
139 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
140 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
141 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
142 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
143 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
144 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
145 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
146 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
147 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
148 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
149 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
150 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
151 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
152 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
153 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
154 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
155 DATA iden/36/
156C ----------
157C BEGIN CODE
158C ----------
159 DO ihel=1,nexternal
160 jc(ihel) = +1
161 ENDDO
162 ans = 0d0
163 DO ihel=1,ncomb
164 t=matrix_usx_udx_h(p ,h1,h2,nhel(1,ihel),jc(1))
165 ans=ans+t
166 ENDDO
167 ans=ans/dble(iden)
168 END
169
170
171 REAL*8 FUNCTION matrix_usx_udx_h(P,H1,H2,NHEL,IC)
172C
173C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
174C By the MadGraph Development Team
175C Please visit us at https://launchpad.net/madgraph5
176C
177C Returns amplitude squared summed/avg over colors
178C for the point with external lines W(0:6,NEXTERNAL)
179C
180C Process: u s~ > u d~ h WEIGHTED=6
181C * Decay: h > ta+ ta- WEIGHTED=2
182C Process: u s~ > c s~ h WEIGHTED=6
183C * Decay: h > ta+ ta- WEIGHTED=2
184C
185 IMPLICIT NONE
186C
187C CONSTANTS
188C
189 INTEGER NGRAPHS
190 parameter(ngraphs=1)
191 INTEGER NEXTERNAL
192 parameter(nexternal=6)
193 INTEGER NWAVEFUNCS, NCOLOR
194 parameter(nwavefuncs=7, ncolor=1)
195 real*8 zero
196 parameter(zero=0d0)
197 COMPLEX*16 IMAG1
198 parameter(imag1=(0d0,1d0))
199C
200C ARGUMENTS
201C
202 real*8 p(0:3,nexternal)
203 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
204C
205C LOCAL VARIABLES
206C
207 INTEGER I,J
208 COMPLEX*16 ZTEMP
209 real*8 denom(ncolor), cf(ncolor,ncolor)
210 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
211 COMPLEX*16 W(18,NWAVEFUNCS)
212 COMPLEX*16 DUM0,DUM1
213 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
214C
215C GLOBAL VARIABLES
216C
217 include 'coupl.inc'
218C
219C COLOR DATA
220C
221 DATA denom(1)/1/
222 DATA (cf(i, 1),i= 1, 1) / 9/
223C 1 T(2,1) T(3,4)
224
225 INTEGER H1,H2
226 real*8 matrix
227 matrix_usx_udx_h=0.d0
228 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
229 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
230
231
232C ----------
233C BEGIN CODE
234C ----------
235 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
236 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
237 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
238 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
239 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
240 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
241 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
242 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
243 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
244C Amplitude(s) for diagram number 1
245 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
246 jamp(1)=-amp(1)
247
248 matrix = 0.d0
249 DO i = 1, ncolor
250 ztemp = (0.d0,0.d0)
251 DO j = 1, ncolor
252 ztemp = ztemp + cf(j,i)*jamp(j)
253 ENDDO
254 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
255 ENDDO
256
257 matrix_usx_udx_h=matrix
258 ENDIF
259 ENDIF
260 END
261
262C ----- begin subprocesses USX->USX with Higgs->tautau
263
264 SUBROUTINE usx_usx_h(P,H1,H2,ANS)
265C
266C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
267C By the MadGraph Development Team
268C Please visit us at https://launchpad.net/madgraph5
269C
270C MadGraph StandAlone Version
271C
272C Returns amplitude squared summed/avg over colors
273C and helicities
274C for the point in phase space P(0:3,NEXTERNAL)
275C
276C Process: u s~ > u s~ h WEIGHTED=6
277C * Decay: h > ta+ ta- WEIGHTED=2
278C
279 IMPLICIT NONE
280C
281C CONSTANTS
282C
283 INTEGER NEXTERNAL
284 parameter(nexternal=6)
285 INTEGER NCOMB
286 parameter( ncomb=64)
287C
288C ARGUMENTS
289C
290 real*8 p(0:3,nexternal),ans
291 INTEGER H1,H2
292C
293C LOCAL VARIABLES
294C
295 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
296 real*8 t
297 real*8 matrix_usx_usx_h
298 INTEGER IHEL,IDEN, I
299 INTEGER JC(NEXTERNAL)
300 LOGICAL GOODHEL(NCOMB)
301 DATA ntry/0/
302 DATA goodhel/ncomb*.false./
303 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
304 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
305 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
306 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
307 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
308 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
309 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
310 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
311 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
312 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
313 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
314 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
315 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
316 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
317 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
318 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
319 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
320 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
321 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
322 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
323 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
324 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
325 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
326 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
327 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
328 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
329 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
330 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
331 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
332 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
333 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
334 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
335 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
336 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
337 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
338 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
339 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
340 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
341 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
342 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
343 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
344 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
345 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
346 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
347 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
348 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
349 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
350 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
351 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
352 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
353 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
354 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
355 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
356 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
357 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
358 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
359 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
360 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
361 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
362 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
363 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
364 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
365 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
366 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
367 DATA iden/36/
368C ----------
369C BEGIN CODE
370C ----------
371 DO ihel=1,nexternal
372 jc(ihel) = +1
373 ENDDO
374 ans = 0d0
375 DO ihel=1,ncomb
376 t=matrix_usx_usx_h(p ,h1,h2,nhel(1,ihel),jc(1))
377 ans=ans+t
378 ENDDO
379 ans=ans/dble(iden)
380 END
381
382
383 REAL*8 FUNCTION matrix_usx_usx_h(P,H1,H2,NHEL,IC)
384C
385C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
386C By the MadGraph Development Team
387C Please visit us at https://launchpad.net/madgraph5
388C
389C Returns amplitude squared summed/avg over colors
390C for the point with external lines W(0:6,NEXTERNAL)
391C
392C Process: u s~ > u s~ h WEIGHTED=6
393C * Decay: h > ta+ ta- WEIGHTED=2
394C
395 IMPLICIT NONE
396C
397C CONSTANTS
398C
399 INTEGER NGRAPHS
400 parameter(ngraphs=2)
401 INTEGER NEXTERNAL
402 parameter(nexternal=6)
403 INTEGER NWAVEFUNCS, NCOLOR
404 parameter(nwavefuncs=7, ncolor=2)
405 real*8 zero
406 parameter(zero=0d0)
407 COMPLEX*16 IMAG1
408 parameter(imag1=(0d0,1d0))
409C
410C ARGUMENTS
411C
412 real*8 p(0:3,nexternal)
413 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
414C
415C LOCAL VARIABLES
416C
417 INTEGER I,J
418 COMPLEX*16 ZTEMP
419 real*8 denom(ncolor), cf(ncolor,ncolor)
420 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
421 COMPLEX*16 W(18,NWAVEFUNCS)
422 COMPLEX*16 DUM0,DUM1
423 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
424C
425C GLOBAL VARIABLES
426C
427 include 'coupl.inc'
428C
429C COLOR DATA
430C
431 DATA denom(1)/1/
432 DATA (cf(i, 1),i= 1, 2) / 9, 3/
433C 1 T(2,1) T(3,4)
434 DATA denom(2)/1/
435 DATA (cf(i, 2),i= 1, 2) / 3, 9/
436C 1 T(2,4) T(3,1)
437
438 INTEGER H1,H2
439 real*8 matrix
440 matrix_usx_usx_h=0.d0
441 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
442 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
443
444
445C ----------
446C BEGIN CODE
447C ----------
448 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
449 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
450 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
451 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
452 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
453 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
454 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
455 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
456 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,5))
457C Amplitude(s) for diagram number 1
458 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_72,amp(1))
459 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
460 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
461C Amplitude(s) for diagram number 2
462 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
463 jamp(1)=-amp(1)
464 jamp(2)=+amp(2)
465
466 matrix = 0.d0
467 DO i = 1, ncolor
468 ztemp = (0.d0,0.d0)
469 DO j = 1, ncolor
470 ztemp = ztemp + cf(j,i)*jamp(j)
471 ENDDO
472 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
473 ENDDO
474 matrix_usx_usx_h=matrix
475 ENDIF
476 ENDIF
477 END
478
479C ----- begin subprocesses USX->CDX with Higgs->tautau
480
481 SUBROUTINE usx_cdx_h(P,H1,H2,ANS)
482C
483C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
484C By the MadGraph Development Team
485C Please visit us at https://launchpad.net/madgraph5
486C
487C MadGraph StandAlone Version
488C
489C Returns amplitude squared summed/avg over colors
490C and helicities
491C for the point in phase space P(0:3,NEXTERNAL)
492C
493C Process: u s~ > c d~ h WEIGHTED=6
494C * Decay: h > ta+ ta- WEIGHTED=2
495C
496 IMPLICIT NONE
497C
498C CONSTANTS
499C
500 INTEGER NEXTERNAL
501 parameter(nexternal=6)
502 INTEGER NCOMB
503 parameter( ncomb=64)
504C
505C ARGUMENTS
506C
507 real*8 p(0:3,nexternal),ans
508 INTEGER H1,H2
509C
510C LOCAL VARIABLES
511C
512 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
513 real*8 t
514 real*8 matrix_usx_cdx_h
515 INTEGER IHEL,IDEN, I
516 INTEGER JC(NEXTERNAL)
517 LOGICAL GOODHEL(NCOMB)
518 DATA ntry/0/
519 DATA goodhel/ncomb*.false./
520 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
521 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
522 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
523 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
524 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
525 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
526 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
527 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
528 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
529 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
530 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
531 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
532 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
533 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
534 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
535 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
536 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
537 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
538 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
539 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
540 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
541 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
542 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
543 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
544 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
545 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
546 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
547 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
548 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
549 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
550 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
551 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
552 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
553 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
554 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
555 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
556 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
557 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
558 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
559 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
560 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
561 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
562 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
563 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
564 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
565 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
566 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
567 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
568 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
569 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
570 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
571 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
572 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
573 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
574 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
575 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
576 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
577 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
578 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
579 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
580 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
581 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
582 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
583 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
584 DATA iden/36/
585C ----------
586C BEGIN CODE
587C ----------
588 DO ihel=1,nexternal
589 jc(ihel) = +1
590 ENDDO
591 ans = 0d0
592 DO ihel=1,ncomb
593 t=matrix_usx_cdx_h(p ,h1,h2,nhel(1,ihel),jc(1))
594 ans=ans+t
595 ENDDO
596 ans=ans/dble(iden)
597 END
598
599
600 REAL*8 FUNCTION matrix_usx_cdx_h(P,H1,H2,NHEL,IC)
601C
602C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
603C By the MadGraph Development Team
604C Please visit us at https://launchpad.net/madgraph5
605C
606C Returns amplitude squared summed/avg over colors
607C for the point with external lines W(0:6,NEXTERNAL)
608C
609C Process: u s~ > c d~ h WEIGHTED=6
610C * Decay: h > ta+ ta- WEIGHTED=2
611C
612 IMPLICIT NONE
613C
614C CONSTANTS
615C
616 INTEGER NGRAPHS
617 parameter(ngraphs=1)
618 INTEGER NEXTERNAL
619 parameter(nexternal=6)
620 INTEGER NWAVEFUNCS, NCOLOR
621 parameter(nwavefuncs=7, ncolor=1)
622 real*8 zero
623 parameter(zero=0d0)
624 COMPLEX*16 IMAG1
625 parameter(imag1=(0d0,1d0))
626C
627C ARGUMENTS
628C
629 real*8 p(0:3,nexternal)
630 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
631C
632C LOCAL VARIABLES
633C
634 INTEGER I,J
635 COMPLEX*16 ZTEMP
636 real*8 denom(ncolor), cf(ncolor,ncolor)
637 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
638 COMPLEX*16 W(18,NWAVEFUNCS)
639 COMPLEX*16 DUM0,DUM1
640 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
641C
642C GLOBAL VARIABLES
643C
644 include 'coupl.inc'
645C
646C COLOR DATA
647C
648 DATA denom(1)/1/
649 DATA (cf(i, 1),i= 1, 1) / 9/
650C 1 T(2,1) T(3,4)
651
652 INTEGER H1,H2
653 real*8 matrix
654 matrix_usx_cdx_h=0.d0
655 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
656 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
657
658C ----------
659C BEGIN CODE
660C ----------
661 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
662 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
663 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
664 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
665 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
666 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
667 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
668 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,6))
669 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,2))
670C Amplitude(s) for diagram number 1
671 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
672 jamp(1)=-amp(1)
673
674 matrix = 0.d0
675 DO i = 1, ncolor
676 ztemp = (0.d0,0.d0)
677 DO j = 1, ncolor
678 ztemp = ztemp + cf(j,i)*jamp(j)
679 ENDDO
680 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
681 ENDDO
682 matrix_usx_cdx_h=matrix
683 ENDIF
684 ENDIF
685 END
686
687C ----- begin subprocesses USX->CSX with Higgs->tautau
688C ----- THE SAME AS USX->UDX WITH Higgs->tautau
689
690C -------------- NOW PROCESSIS WITHOUT HIGGS
691
692C ----- begin subprocesses USX->UDX NO Higgs
693
694 SUBROUTINE usx_udx_noh(P,H1,H2,ANS)
695C
696C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
697C By the MadGraph Development Team
698C Please visit us at https://launchpad.net/madgraph5
699C
700C MadGraph StandAlone Version
701C
702C Returns amplitude squared summed/avg over colors
703C and helicities
704C for the point in phase space P(0:3,NEXTERNAL)
705C
706C Process: u s~ > u d~ ta+ ta- / h QED=4
707C Process: u s~ > c s~ ta+ ta- / h QED=4
708C
709 IMPLICIT NONE
710C
711C CONSTANTS
712C
713 INTEGER NEXTERNAL
714 parameter(nexternal=6)
715 INTEGER NCOMB
716 parameter( ncomb=64)
717C
718C ARGUMENTS
719C
720 real*8 p(0:3,nexternal),ans
721 INTEGER H1,H2
722C
723C LOCAL VARIABLES
724C
725 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
726 real*8 t
727 real*8 matrix_usx_udx_noh
728 INTEGER IHEL,IDEN, I
729 INTEGER JC(NEXTERNAL)
730 LOGICAL GOODHEL(NCOMB)
731 DATA ntry/0/
732 DATA goodhel/ncomb*.false./
733 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
734 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
735 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
736 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
737 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
738 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
739 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
740 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
741 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
742 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
743 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
744 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
745 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
746 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
747 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
748 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
749 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
750 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
751 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
752 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
753 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
754 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
755 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
756 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
757 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
758 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
759 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
760 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
761 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
762 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
763 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
764 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
765 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
766 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
767 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
768 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
769 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
770 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
771 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
772 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
773 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
774 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
775 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
776 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
777 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
778 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
779 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
780 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
781 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
782 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
783 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
784 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
785 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
786 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
787 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
788 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
789 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
790 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
791 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
792 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
793 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
794 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
795 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
796 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
797 DATA iden/36/
798C ----------
799C BEGIN CODE
800C ----------
801 DO ihel=1,nexternal
802 jc(ihel) = +1
803 ENDDO
804 ans = 0d0
805 DO ihel=1,ncomb
806 t=matrix_usx_udx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
807 ans=ans+t
808 ENDDO
809 ans=ans/dble(iden)
810 END
811
812
813 REAL*8 FUNCTION matrix_usx_udx_noh(P,H1,H2,NHEL,IC)
814C
815C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
816C By the MadGraph Development Team
817C Please visit us at https://launchpad.net/madgraph5
818C
819C Returns amplitude squared summed/avg over colors
820C for the point with external lines W(0:6,NEXTERNAL)
821C
822C Process: u s~ > u d~ ta+ ta- / h QED=4
823C Process: u s~ > c s~ ta+ ta- / h QED=4
824C
825 IMPLICIT NONE
826C
827C CONSTANTS
828C
829 INTEGER NGRAPHS
830 parameter(ngraphs=11)
831 INTEGER NEXTERNAL
832 parameter(nexternal=6)
833 INTEGER NWAVEFUNCS, NCOLOR
834 parameter(nwavefuncs=10, ncolor=1)
835 real*8 zero
836 parameter(zero=0d0)
837 COMPLEX*16 IMAG1
838 parameter(imag1=(0d0,1d0))
839C
840C ARGUMENTS
841C
842 real*8 p(0:3,nexternal)
843 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
844C
845C LOCAL VARIABLES
846C
847 INTEGER I,J
848 COMPLEX*16 ZTEMP
849 real*8 denom(ncolor), cf(ncolor,ncolor)
850 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
851 COMPLEX*16 W(18,NWAVEFUNCS)
852 COMPLEX*16 DUM0,DUM1
853 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
854C
855C GLOBAL VARIABLES
856C
857 include 'coupl.inc'
858C
859C COLOR DATA
860C
861 DATA denom(1)/1/
862 DATA (cf(i, 1),i= 1, 1) / 9/
863C 1 T(2,1) T(3,4)
864
865 INTEGER H1,H2
866 real*8 matrix
867 matrix_usx_udx_noh=0.d0
868 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
869 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
870
871C ----------
872C BEGIN CODE
873C ----------
874 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
875 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
876 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
877 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
878 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
879 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
880 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
881 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
882 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
883C Amplitude(s) for diagram number 1
884 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
885 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
886C Amplitude(s) for diagram number 2
887 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
888 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
889C Amplitude(s) for diagram number 3
890 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
891 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
892C Amplitude(s) for diagram number 4
893 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
894 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
895C Amplitude(s) for diagram number 5
896 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
897C Amplitude(s) for diagram number 6
898 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
899C Amplitude(s) for diagram number 7
900 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
901 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
902C Amplitude(s) for diagram number 8
903 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
904 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
905C Amplitude(s) for diagram number 9
906 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
907C Amplitude(s) for diagram number 10
908 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
909 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
910C Amplitude(s) for diagram number 11
911 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
912 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
913 $ -amp(9)-amp(10)-amp(11)
914
915 matrix = 0.d0
916 DO i = 1, ncolor
917 ztemp = (0.d0,0.d0)
918 DO j = 1, ncolor
919 ztemp = ztemp + cf(j,i)*jamp(j)
920 ENDDO
921 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
922 ENDDO
923 matrix_usx_udx_noh=matrix
924 ENDIF
925 ENDIF
926 END
927
928C ----- begin subprocesses USX->USX NO Higgs
929
930 SUBROUTINE usx_usx_noh(P,H1,H2,ANS)
931C
932C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
933C By the MadGraph Development Team
934C Please visit us at https://launchpad.net/madgraph5
935C
936C MadGraph StandAlone Version
937C
938C Returns amplitude squared summed/avg over colors
939C and helicities
940C for the point in phase space P(0:3,NEXTERNAL)
941C
942C Process: u s~ > u s~ ta+ ta- / h QED=4
943C
944 IMPLICIT NONE
945C
946C CONSTANTS
947C
948 INTEGER NEXTERNAL
949 parameter(nexternal=6)
950 INTEGER NCOMB
951 parameter( ncomb=64)
952C
953C ARGUMENTS
954C
955 real*8 p(0:3,nexternal),ans
956 INTEGER H1,H2
957C
958C LOCAL VARIABLES
959C
960 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
961 real*8 t
962 real*8 matrix_usx_usx_noh
963 INTEGER IHEL,IDEN, I
964 INTEGER JC(NEXTERNAL)
965 LOGICAL GOODHEL(NCOMB)
966 DATA ntry/0/
967 DATA goodhel/ncomb*.false./
968 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
969 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
970 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
971 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
972 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
973 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
974 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
975 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
976 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
977 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
978 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
979 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
980 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
981 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
982 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
983 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
984 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
985 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
986 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
987 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
988 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
989 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
990 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
991 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
992 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
993 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
994 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
995 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
996 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
997 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
998 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
999 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1000 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1001 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1002 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1003 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1004 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1005 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1006 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1007 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1008 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1009 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1010 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1011 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1012 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1013 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1014 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1015 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1016 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1017 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1018 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1019 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1020 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1021 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1022 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1023 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1024 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1025 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1026 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1027 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1028 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1029 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1030 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1031 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1032 DATA iden/36/
1033C ----------
1034C BEGIN CODE
1035C ----------
1036 DO ihel=1,nexternal
1037 jc(ihel) = +1
1038 ENDDO
1039 ans = 0d0
1040 DO ihel=1,ncomb
1041 t=matrix_usx_usx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1042 ans=ans+t
1043 ENDDO
1044 ans=ans/dble(iden)
1045 END
1046
1047
1048 REAL*8 FUNCTION matrix_usx_usx_noh(P,H1,H2,NHEL,IC)
1049C
1050C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1051C By the MadGraph Development Team
1052C Please visit us at https://launchpad.net/madgraph5
1053C
1054C Returns amplitude squared summed/avg over colors
1055C for the point with external lines W(0:6,NEXTERNAL)
1056C
1057C Process: u s~ > u s~ ta+ ta- / h QED=4
1058C
1059 IMPLICIT NONE
1060C
1061C CONSTANTS
1062C
1063 INTEGER NGRAPHS
1064 parameter(ngraphs=43)
1065 INTEGER NEXTERNAL
1066 parameter(nexternal=6)
1067 INTEGER NWAVEFUNCS, NCOLOR
1068 parameter(nwavefuncs=15, ncolor=2)
1069 real*8 zero
1070 parameter(zero=0d0)
1071 COMPLEX*16 IMAG1
1072 parameter(imag1=(0d0,1d0))
1073C
1074C ARGUMENTS
1075C
1076 real*8 p(0:3,nexternal)
1077 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1078C
1079C LOCAL VARIABLES
1080C
1081 INTEGER I,J
1082 COMPLEX*16 ZTEMP
1083 real*8 denom(ncolor), cf(ncolor,ncolor)
1084 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1085 COMPLEX*16 W(18,NWAVEFUNCS)
1086 COMPLEX*16 DUM0,DUM1
1087 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1088C
1089C GLOBAL VARIABLES
1090C
1091 include 'coupl.inc'
1092C
1093C COLOR DATA
1094C
1095 DATA denom(1)/1/
1096 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1097C 1 T(2,1) T(3,4)
1098 DATA denom(2)/1/
1099 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1100C 1 T(2,4) T(3,1)
1101
1102 INTEGER H1,H2
1103 real*8 matrix
1104 matrix_usx_usx_noh=0.d0
1105 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1106 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1107
1108C ----------
1109C BEGIN CODE
1110C ----------
1111 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1112 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1113 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1114 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1115 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1116 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1117 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1118 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1119 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1120C Amplitude(s) for diagram number 1
1121 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1122 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1123C Amplitude(s) for diagram number 2
1124 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1125 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1126C Amplitude(s) for diagram number 3
1127 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1128 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,11))
1129C Amplitude(s) for diagram number 4
1130 CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_1,amp(4))
1131 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,12))
1132C Amplitude(s) for diagram number 5
1133 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(5))
1134C Amplitude(s) for diagram number 6
1135 CALL ffv2_3_0(w(1,4),w(1,11),w(1,10),gc_50,gc_58,amp(6))
1136C Amplitude(s) for diagram number 7
1137 CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1138 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,12))
1139 CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,11))
1140 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1141C Amplitude(s) for diagram number 8
1142 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1143 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1144C Amplitude(s) for diagram number 9
1145 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1146 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,14))
1147C Amplitude(s) for diagram number 10
1148 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1149C Amplitude(s) for diagram number 11
1150 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1151 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,13))
1152 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1153C Amplitude(s) for diagram number 12
1154 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1155 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1156C Amplitude(s) for diagram number 13
1157 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1158C Amplitude(s) for diagram number 14
1159 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1160C Amplitude(s) for diagram number 15
1161 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1162 CALL ffv1_1(w(1,2),w(1,12),gc_1,zero,zero,w(1,15))
1163C Amplitude(s) for diagram number 16
1164 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(16))
1165 CALL ffv1_2(w(1,4),w(1,12),gc_1,zero,zero,w(1,5))
1166C Amplitude(s) for diagram number 17
1167 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(17))
1168C Amplitude(s) for diagram number 18
1169 CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(18))
1170C Amplitude(s) for diagram number 19
1171 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(19))
1172 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1173 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1174C Amplitude(s) for diagram number 20
1175 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(20))
1176 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1177C Amplitude(s) for diagram number 21
1178 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_1,amp(21))
1179C Amplitude(s) for diagram number 22
1180 CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(22))
1181C Amplitude(s) for diagram number 23
1182 CALL ffv2_3_0(w(1,12),w(1,2),w(1,10),gc_50,gc_58,amp(23))
1183 CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,12))
1184C Amplitude(s) for diagram number 24
1185 CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_1,amp(24))
1186 CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,15))
1187C Amplitude(s) for diagram number 25
1188 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_1,amp(25))
1189C Amplitude(s) for diagram number 26
1190 CALL ffv2_3_0(w(1,4),w(1,12),w(1,10),gc_50,gc_58,amp(26))
1191C Amplitude(s) for diagram number 27
1192 CALL ffv2_3_0(w(1,15),w(1,2),w(1,10),gc_50,gc_58,amp(27))
1193 CALL ffv1_2(w(1,1),w(1,11),gc_2,zero,zero,w(1,15))
1194C Amplitude(s) for diagram number 28
1195 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(28))
1196 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,12))
1197C Amplitude(s) for diagram number 29
1198 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_2,amp(29))
1199C Amplitude(s) for diagram number 30
1200 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(30))
1201 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,15))
1202C Amplitude(s) for diagram number 31
1203 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_2,amp(31))
1204 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1205 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1206C Amplitude(s) for diagram number 32
1207 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_2,amp(32))
1208C Amplitude(s) for diagram number 33
1209 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1210C Amplitude(s) for diagram number 34
1211 CALL ffv2_5_0(w(1,4),w(1,3),w(1,10),gc_51,gc_58,amp(34))
1212C Amplitude(s) for diagram number 35
1213 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1214 CALL ffv2_5_2(w(1,1),w(1,14),gc_51,gc_58,zero,zero,w(1,11))
1215C Amplitude(s) for diagram number 36
1216 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_2,amp(36))
1217C Amplitude(s) for diagram number 37
1218 CALL ffv2_5_0(w(1,12),w(1,3),w(1,14),gc_51,gc_58,amp(37))
1219C Amplitude(s) for diagram number 38
1220 CALL ffv2_5_0(w(1,11),w(1,3),w(1,10),gc_51,gc_58,amp(38))
1221C Amplitude(s) for diagram number 39
1222 CALL ffv2_5_0(w(1,15),w(1,3),w(1,14),gc_51,gc_58,amp(39))
1223 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,14))
1224C Amplitude(s) for diagram number 40
1225 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_1,amp(40))
1226C Amplitude(s) for diagram number 41
1227 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_101,amp(41))
1228C Amplitude(s) for diagram number 42
1229 CALL ffv2_3_0(w(1,14),w(1,2),w(1,10),gc_50,gc_58,amp(42))
1230C Amplitude(s) for diagram number 43
1231 CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_101,amp(43))
1232 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1233 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1234 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1235 $ -amp(40)-amp(41)-amp(42)-amp(43)
1236 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1237 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1238 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1239 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1240 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1241 $ +amp(38)+amp(39)
1242
1243 matrix = 0.d0
1244 DO i = 1, ncolor
1245 ztemp = (0.d0,0.d0)
1246 DO j = 1, ncolor
1247 ztemp = ztemp + cf(j,i)*jamp(j)
1248 ENDDO
1249 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1250 ENDDO
1251 matrix_usx_usx_noh=matrix
1252 ENDIF
1253 ENDIF
1254 END
1255
1256C ----- begin subprocesses USX->CDX NO Higgs
1257
1258 SUBROUTINE usx_cdx_noh(P,H1,H2,ANS)
1259C
1260C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1261C By the MadGraph Development Team
1262C Please visit us at https://launchpad.net/madgraph5
1263C
1264C MadGraph StandAlone Version
1265C
1266C Returns amplitude squared summed/avg over colors
1267C and helicities
1268C for the point in phase space P(0:3,NEXTERNAL)
1269C
1270C Process: u s~ > c d~ ta+ ta- / h QED=4
1271C
1272 IMPLICIT NONE
1273C
1274C CONSTANTS
1275C
1276 INTEGER NEXTERNAL
1277 parameter(nexternal=6)
1278 INTEGER NCOMB
1279 parameter( ncomb=64)
1280C
1281C ARGUMENTS
1282C
1283 real*8 p(0:3,nexternal),ans
1284 INTEGER H1,H2
1285C
1286C LOCAL VARIABLES
1287C
1288 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1289 real*8 t
1290 real*8 matrix_usx_cdx_noh
1291 INTEGER IHEL,IDEN, I
1292 INTEGER JC(NEXTERNAL)
1293 LOGICAL GOODHEL(NCOMB)
1294 DATA ntry/0/
1295 DATA goodhel/ncomb*.false./
1296 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1297 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1298 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1299 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1300 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1301 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1302 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1303 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1304 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1305 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1306 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1307 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1308 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1309 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1310 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1311 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1312 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1313 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1314 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1315 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1316 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1317 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1318 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1319 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1320 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1321 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1322 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1323 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1324 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1325 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1326 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1327 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1328 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1329 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1330 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1331 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1332 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1333 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1334 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1335 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1336 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1337 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1338 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1339 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1340 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1341 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1342 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1343 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1344 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1345 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1346 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1347 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1348 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1349 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1350 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1351 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1352 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1353 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1354 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1355 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1356 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1357 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1358 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1359 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1360 DATA iden/36/
1361C ----------
1362C BEGIN CODE
1363C ----------
1364 DO ihel=1,nexternal
1365 jc(ihel) = +1
1366 ENDDO
1367 ans = 0d0
1368 DO ihel=1,ncomb
1369 t=matrix_usx_cdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1370 ans=ans+t
1371 ENDDO
1372
1373 ans=ans/dble(iden)
1374 END
1375
1376
1377 REAL*8 FUNCTION matrix_usx_cdx_noh(P,H1,H2,NHEL,IC)
1378C
1379C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1380C By the MadGraph Development Team
1381C Please visit us at https://launchpad.net/madgraph5
1382C
1383C Returns amplitude squared summed/avg over colors
1384C for the point with external lines W(0:6,NEXTERNAL)
1385C
1386C Process: u s~ > c d~ ta+ ta- / h QED=4
1387C
1388 IMPLICIT NONE
1389C
1390C CONSTANTS
1391C
1392 INTEGER NGRAPHS
1393 parameter(ngraphs=11)
1394 INTEGER NEXTERNAL
1395 parameter(nexternal=6)
1396 INTEGER NWAVEFUNCS, NCOLOR
1397 parameter(nwavefuncs=10, ncolor=1)
1398 real*8 zero
1399 parameter(zero=0d0)
1400 COMPLEX*16 IMAG1
1401 parameter(imag1=(0d0,1d0))
1402C
1403C ARGUMENTS
1404C
1405 real*8 p(0:3,nexternal)
1406 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1407C
1408C LOCAL VARIABLES
1409C
1410 INTEGER I,J
1411 COMPLEX*16 ZTEMP
1412 real*8 denom(ncolor), cf(ncolor,ncolor)
1413 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1414 COMPLEX*16 W(18,NWAVEFUNCS)
1415 COMPLEX*16 DUM0,DUM1
1416 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1417C
1418C GLOBAL VARIABLES
1419C
1420 include 'coupl.inc'
1421C
1422C COLOR DATA
1423C
1424 DATA denom(1)/1/
1425 DATA (cf(i, 1),i= 1, 1) / 9/
1426C 1 T(2,1) T(3,4)
1427
1428 INTEGER H1,H2
1429 real*8 matrix
1430 matrix_usx_cdx_noh=0.d0
1431 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1432 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1433
1434C ----------
1435C BEGIN CODE
1436C ----------
1437 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1438 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1439 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1440 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1441 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1442 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1443 CALL ffv2_3(w(1,1),w(1,2),gc_101,mw,ww,w(1,7))
1444 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1445 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1446C Amplitude(s) for diagram number 1
1447 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1448 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1449C Amplitude(s) for diagram number 2
1450 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1451 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1452C Amplitude(s) for diagram number 3
1453 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1454 CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,6))
1455C Amplitude(s) for diagram number 4
1456 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1457 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
1458C Amplitude(s) for diagram number 5
1459 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1460C Amplitude(s) for diagram number 6
1461 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1462C Amplitude(s) for diagram number 7
1463 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1464 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
1465C Amplitude(s) for diagram number 8
1466 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1467 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1468C Amplitude(s) for diagram number 9
1469 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_101,amp(9))
1470C Amplitude(s) for diagram number 10
1471 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1472 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1473C Amplitude(s) for diagram number 11
1474 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_101,amp(11))
1475 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1476 $ -amp(9)-amp(10)-amp(11)
1477
1478 matrix = 0.d0
1479 DO i = 1, ncolor
1480 ztemp = (0.d0,0.d0)
1481 DO j = 1, ncolor
1482 ztemp = ztemp + cf(j,i)*jamp(j)
1483 ENDDO
1484 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1485 ENDDO
1486 matrix_usx_cdx_noh=matrix
1487 ENDIF
1488 ENDIF
1489 END
1490
1491C ----- begin subprocesses USX->CSX NO Higgs
1492C IS THE SAME AS USX->CSX NO Higgs
1493
1494
1495
1496C DONE 25.10.2015 STILL S=2 MISSING