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