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