C++InterfacetoTauola
CS.f
1 
2 c --------------------- 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 
10 C 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 
46 C ----- begin subprocesses CS->jjH with Higgs->tautau
47 
48  SUBROUTINE cs_cd_h(P,H1,H2,ANS)
49 C
50 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
51 C By the MadGraph Development Team
52 C Please visit us at https://launchpad.net/madgraph5
53 C
54 C MadGraph StandAlone Version
55 C
56 C Returns amplitude squared summed/avg over colors
57 C and helicities
58 C for the point in phase space P(0:3,NEXTERNAL)
59 C
60 C Process: c s > c d h WEIGHTED=6
61 C * Decay: h > ta+ ta- WEIGHTED=2
62 C
63  IMPLICIT NONE
64 C
65 C CONSTANTS
66 C
67  INTEGER nexternal
68  parameter(nexternal=6)
69  INTEGER ncomb
70  parameter( ncomb=64)
71 C
72 C ARGUMENTS
73 C
74  REAL*8 p(0:3,nexternal),ans
75  INTEGER h1,h2
76 C
77 C LOCAL VARIABLES
78 C
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/
152 C ----------
153 C BEGIN CODE
154 C ----------
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)
168 C
169 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
170 C By the MadGraph Development Team
171 C Please visit us at https://launchpad.net/madgraph5
172 C
173 C Returns amplitude squared summed/avg over colors
174 C for the point with external lines W(0:6,NEXTERNAL)
175 C
176 C Process: c s > c d h WEIGHTED=6
177 C * Decay: h > ta+ ta- WEIGHTED=2
178 C
179  IMPLICIT NONE
180 C
181 C CONSTANTS
182 C
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))
193 C
194 C ARGUMENTS
195 C
196  REAL*8 p(0:3,nexternal)
197  INTEGER nhel(nexternal), ic(nexternal)
198 C
199 C LOCAL VARIABLES
200 C
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)/
208 C
209 C GLOBAL VARIABLES
210 C
211  include 'coupl.inc'
212 C
213 C COLOR DATA
214 C
215  DATA denom(1)/1/
216  DATA (cf(i, 1),i= 1, 1) / 9/
217 C 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 
225 C ----------
226 C BEGIN CODE
227 C ----------
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))
237 C 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 
254 c==========================================
255 
256  SUBROUTINE cs_cs_h(P,H1,H2,ANS)
257 C
258 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
259 C By the MadGraph Development Team
260 C Please visit us at https://launchpad.net/madgraph5
261 C
262 C MadGraph StandAlone Version
263 C
264 C Returns amplitude squared summed/avg over colors
265 C and helicities
266 C for the point in phase space P(0:3,NEXTERNAL)
267 C
268 C Process: c s > c s h WEIGHTED=6
269 C * Decay: h > ta+ ta- WEIGHTED=2
270 C
271  IMPLICIT NONE
272 C
273 C CONSTANTS
274 C
275  INTEGER nexternal
276  parameter(nexternal=6)
277  INTEGER ncomb
278  parameter( ncomb=64)
279 C
280 C ARGUMENTS
281 C
282  REAL*8 p(0:3,nexternal),ans
283  INTEGER h1,h2
284 C
285 C LOCAL VARIABLES
286 C
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/
360 C ----------
361 C BEGIN CODE
362 C ----------
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)
376 C
377 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
378 C By the MadGraph Development Team
379 C Please visit us at https://launchpad.net/madgraph5
380 C
381 C Returns amplitude squared summed/avg over colors
382 C for the point with external lines W(0:6,NEXTERNAL)
383 C
384 C Process: c s > c s h WEIGHTED=6
385 C * Decay: h > ta+ ta- WEIGHTED=2
386 C
387  IMPLICIT NONE
388 C
389 C CONSTANTS
390 C
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))
401 C
402 C ARGUMENTS
403 C
404  REAL*8 p(0:3,nexternal)
405  INTEGER nhel(nexternal), ic(nexternal)
406 C
407 C LOCAL VARIABLES
408 C
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)/
416 C
417 C GLOBAL VARIABLES
418 C
419  include 'coupl.inc'
420 C
421 C COLOR DATA
422 C
423  DATA denom(1)/1/
424  DATA (cf(i, 1),i= 1, 2) / 9, 3/
425 C 1 T(3,1) T(4,2)
426  DATA denom(2)/1/
427  DATA (cf(i, 2),i= 1, 2) / 3, 9/
428 C 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 
438 C ----------
439 C BEGIN CODE
440 C ----------
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))
450 C 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))
454 C 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 
474 C===============================================
475 
476  SUBROUTINE cs_ud_h(P,H1,H2,ANS)
477 C
478 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
479 C By the MadGraph Development Team
480 C Please visit us at https://launchpad.net/madgraph5
481 C
482 C MadGraph StandAlone Version
483 C
484 C Returns amplitude squared summed/avg over colors
485 C and helicities
486 C for the point in phase space P(0:3,NEXTERNAL)
487 C
488 C Process: c s > u d h WEIGHTED=6
489 C * Decay: h > ta+ ta- WEIGHTED=2
490 C
491  IMPLICIT NONE
492 C
493 C CONSTANTS
494 C
495  INTEGER nexternal
496  parameter(nexternal=6)
497  INTEGER ncomb
498  parameter( ncomb=64)
499 C
500 C ARGUMENTS
501 C
502  REAL*8 p(0:3,nexternal),ans
503  INTEGER h1,h2
504 C
505 C LOCAL VARIABLES
506 C
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/
580 C ----------
581 C BEGIN CODE
582 C ----------
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)
596 C
597 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
598 C By the MadGraph Development Team
599 C Please visit us at https://launchpad.net/madgraph5
600 C
601 C Returns amplitude squared summed/avg over colors
602 C for the point with external lines W(0:6,NEXTERNAL)
603 C
604 C Process: c s > u d h WEIGHTED=6
605 C * Decay: h > ta+ ta- WEIGHTED=2
606 C
607  IMPLICIT NONE
608 C
609 C CONSTANTS
610 C
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))
621 C
622 C ARGUMENTS
623 C
624  REAL*8 p(0:3,nexternal)
625  INTEGER nhel(nexternal), ic(nexternal)
626 C
627 C LOCAL VARIABLES
628 C
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)/
636 C
637 C GLOBAL VARIABLES
638 C
639  include 'coupl.inc'
640 C
641 C COLOR DATA
642 C
643  DATA denom(1)/1/
644  DATA (cf(i, 1),i= 1, 1) / 9/
645 C 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 
654 C ----------
655 C BEGIN CODE
656 C ----------
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))
666 C 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 
685 C========================================
686 
687  SUBROUTINE cs_us_h(P,H1,H2,ANS)
688 C
689 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
690 C By the MadGraph Development Team
691 C Please visit us at https://launchpad.net/madgraph5
692 C
693 C MadGraph StandAlone Version
694 C
695 C Returns amplitude squared summed/avg over colors
696 C and helicities
697 C for the point in phase space P(0:3,NEXTERNAL)
698 C
699 C Process: c s > u s h WEIGHTED=6
700 C * Decay: h > ta+ ta- WEIGHTED=2
701 C
702  IMPLICIT NONE
703 C
704 C CONSTANTS
705 C
706  INTEGER nexternal
707  parameter(nexternal=6)
708  INTEGER ncomb
709  parameter( ncomb=64)
710 C
711 C ARGUMENTS
712 C
713  REAL*8 p(0:3,nexternal),ans
714  INTEGER h1,h2
715 C
716 C LOCAL VARIABLES
717 C
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/
791 C ----------
792 C BEGIN CODE
793 C ----------
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)
807 C
808 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
809 C By the MadGraph Development Team
810 C Please visit us at https://launchpad.net/madgraph5
811 C
812 C Returns amplitude squared summed/avg over colors
813 C for the point with external lines W(0:6,NEXTERNAL)
814 C
815 C Process: c s > u s h WEIGHTED=6
816 C * Decay: h > ta+ ta- WEIGHTED=2
817 C
818  IMPLICIT NONE
819 C
820 C CONSTANTS
821 C
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))
832 C
833 C ARGUMENTS
834 C
835  REAL*8 p(0:3,nexternal)
836  INTEGER nhel(nexternal), ic(nexternal)
837 C
838 C LOCAL VARIABLES
839 C
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)/
847 C
848 C GLOBAL VARIABLES
849 C
850  include 'coupl.inc'
851 C
852 C COLOR DATA
853 C
854  DATA denom(1)/1/
855  DATA (cf(i, 1),i= 1, 1) / 9/
856 C 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 
865 C ----------
866 C BEGIN CODE
867 C ----------
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))
877 C 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 
897 C --- end of subprocess CS->jjH, H->tautau
898 
899 C ----- begin subprocesses CS->jjtautau, no Higgs
900 
901  SUBROUTINE cs_cd_noh(P,H1,H2,ANS)
902 C
903 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
904 C By the MadGraph Development Team
905 C Please visit us at https://launchpad.net/madgraph5
906 C
907 C MadGraph StandAlone Version
908 C
909 C Returns amplitude squared summed/avg over colors
910 C and helicities
911 C for the point in phase space P(0:3,NEXTERNAL)
912 C
913 C Process: c s > c d ta+ ta- / h QED=4
914 C
915  IMPLICIT NONE
916 C
917 C CONSTANTS
918 C
919  INTEGER nexternal
920  parameter(nexternal=6)
921  INTEGER ncomb
922  parameter( ncomb=64)
923 C
924 C ARGUMENTS
925 C
926  REAL*8 p(0:3,nexternal),ans
927  INTEGER h1,h2
928 C
929 C LOCAL VARIABLES
930 C
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/
1004 C ----------
1005 C BEGIN CODE
1006 C ----------
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)
1020 C
1021 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1022 C By the MadGraph Development Team
1023 C Please visit us at https://launchpad.net/madgraph5
1024 C
1025 C Returns amplitude squared summed/avg over colors
1026 C for the point with external lines W(0:6,NEXTERNAL)
1027 C
1028 C Process: c s > c d ta+ ta- / h QED=4
1029 C
1030  IMPLICIT NONE
1031 C
1032 C CONSTANTS
1033 C
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))
1044 C
1045 C ARGUMENTS
1046 C
1047  REAL*8 p(0:3,nexternal)
1048  INTEGER nhel(nexternal), ic(nexternal)
1049 C
1050 C LOCAL VARIABLES
1051 C
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)/
1059 C
1060 C GLOBAL VARIABLES
1061 C
1062  include 'coupl.inc'
1063 C
1064 C COLOR DATA
1065 C
1066  DATA denom(1)/1/
1067  DATA (cf(i, 1),i= 1, 1) / 9/
1068 C 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 
1076 C ----------
1077 C BEGIN CODE
1078 C ----------
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))
1088 C 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))
1091 C 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))
1094 C 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))
1097 C 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))
1100 C Amplitude(s) for diagram number 5
1101  CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1102 C 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))
1104 C 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))
1107 C 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))
1110 C Amplitude(s) for diagram number 9
1111  CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_44,amp(9))
1112 C 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))
1115 C 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)
1135 C
1136 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1137 C By the MadGraph Development Team
1138 C Please visit us at https://launchpad.net/madgraph5
1139 C
1140 C MadGraph StandAlone Version
1141 C
1142 C Returns amplitude squared summed/avg over colors
1143 C and helicities
1144 C for the point in phase space P(0:3,NEXTERNAL)
1145 C
1146 C Process: c s > c s ta+ ta- / h QED=4
1147 C
1148  IMPLICIT NONE
1149 C
1150 C CONSTANTS
1151 C
1152  INTEGER nexternal
1153  parameter(nexternal=6)
1154  INTEGER ncomb
1155  parameter( ncomb=64)
1156 C
1157 C ARGUMENTS
1158 C
1159  REAL*8 p(0:3,nexternal),ans
1160  INTEGER h1,h2
1161 C
1162 C LOCAL VARIABLES
1163 C
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/
1237 C ----------
1238 C BEGIN CODE
1239 C ----------
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)
1253 C
1254 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1255 C By the MadGraph Development Team
1256 C Please visit us at https://launchpad.net/madgraph5
1257 C
1258 C Returns amplitude squared summed/avg over colors
1259 C for the point with external lines W(0:6,NEXTERNAL)
1260 C
1261 C Process: c s > c s ta+ ta- / h QED=4
1262 C
1263  IMPLICIT NONE
1264 C
1265 C CONSTANTS
1266 C
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))
1277 C
1278 C ARGUMENTS
1279 C
1280  REAL*8 p(0:3,nexternal)
1281  INTEGER nhel(nexternal), ic(nexternal)
1282 C
1283 C LOCAL VARIABLES
1284 C
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)/
1292 C
1293 C GLOBAL VARIABLES
1294 C
1295  include 'coupl.inc'
1296 C
1297 C COLOR DATA
1298 C
1299  DATA denom(1)/1/
1300  DATA (cf(i, 1),i= 1, 2) / 9, 3/
1301 C 1 T(3,1) T(4,2)
1302  DATA denom(2)/1/
1303  DATA (cf(i, 2),i= 1, 2) / 3, 9/
1304 C 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 
1313 C ----------
1314 C BEGIN CODE
1315 C ----------
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))
1325 C 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))
1328 C 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))
1331 C 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))
1333 C 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))
1337 C 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))
1340 C Amplitude(s) for diagram number 6
1341  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1342 C 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))
1344 C 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))
1348 C 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))
1351 C 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))
1354 C 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))
1356 C 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))
1360 C 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))
1363 C Amplitude(s) for diagram number 14
1364  CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_1,amp(14))
1365 C 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))
1367 C 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))
1370 C 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))
1373 C Amplitude(s) for diagram number 18
1374  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
1375 C 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))
1377 C 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))
1382 C Amplitude(s) for diagram number 21
1383  CALL ffv2_0(w(1,10),w(1,6),w(1,14),gc_108,amp(21))
1384 C Amplitude(s) for diagram number 22
1385  CALL vvv1_0(w(1,12),w(1,9),w(1,14),gc_4,amp(22))
1386 C 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))
1389 C 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))
1392 C Amplitude(s) for diagram number 25
1393  CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(25))
1394 C 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))
1396 C 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))
1399 C 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))
1402 C Amplitude(s) for diagram number 29
1403  CALL ffv2_0(w(1,10),w(1,4),w(1,14),gc_100,amp(29))
1404 C 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))
1407 C 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))
1410 C Amplitude(s) for diagram number 32
1411  CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(32))
1412 C Amplitude(s) for diagram number 33
1413  CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_2,amp(33))
1414 C 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))
1416 C 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))
1420 C Amplitude(s) for diagram number 36
1421  CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(36))
1422 C Amplitude(s) for diagram number 37
1423  CALL ffv1_0(w(1,10),w(1,3),w(1,8),gc_11,amp(37))
1424 C 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))
1426 C 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))
1429 C Amplitude(s) for diagram number 40
1430  CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(40))
1431 C 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))
1433 C 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))
1435 C 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)
1461 C
1462 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1463 C By the MadGraph Development Team
1464 C Please visit us at https://launchpad.net/madgraph5
1465 C
1466 C MadGraph StandAlone Version
1467 C
1468 C Returns amplitude squared summed/avg over colors
1469 C and helicities
1470 C for the point in phase space P(0:3,NEXTERNAL)
1471 C
1472 C Process: c s > u d ta+ ta- / h QED=4
1473 C
1474  IMPLICIT NONE
1475 C
1476 C CONSTANTS
1477 C
1478  INTEGER nexternal
1479  parameter(nexternal=6)
1480  INTEGER ncomb
1481  parameter( ncomb=64)
1482 C
1483 C ARGUMENTS
1484 C
1485  REAL*8 p(0:3,nexternal),ans
1486  INTEGER h1,h2
1487 C
1488 C LOCAL VARIABLES
1489 C
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/
1563 C ----------
1564 C BEGIN CODE
1565 C ----------
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)
1579 C
1580 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1581 C By the MadGraph Development Team
1582 C Please visit us at https://launchpad.net/madgraph5
1583 C
1584 C Returns amplitude squared summed/avg over colors
1585 C for the point with external lines W(0:6,NEXTERNAL)
1586 C
1587 C Process: c s > u d ta+ ta- / h QED=4
1588 C
1589  IMPLICIT NONE
1590 C
1591 C CONSTANTS
1592 C
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))
1603 C
1604 C ARGUMENTS
1605 C
1606  REAL*8 p(0:3,nexternal)
1607  INTEGER nhel(nexternal), ic(nexternal)
1608 C
1609 C LOCAL VARIABLES
1610 C
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)/
1618 C
1619 C GLOBAL VARIABLES
1620 C
1621  include 'coupl.inc'
1622 C
1623 C COLOR DATA
1624 C
1625  DATA denom(1)/1/
1626  DATA (cf(i, 1),i= 1, 1) / 9/
1627 C 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 
1636 C ----------
1637 C BEGIN CODE
1638 C ----------
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))
1648 C 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))
1651 C 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))
1654 C 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))
1657 C 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))
1660 C Amplitude(s) for diagram number 5
1661  CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1662 C 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))
1664 C 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))
1667 C 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))
1670 C Amplitude(s) for diagram number 9
1671  CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_44,amp(9))
1672 C 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))
1675 C 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)
1695 C
1696 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1697 C By the MadGraph Development Team
1698 C Please visit us at https://launchpad.net/madgraph5
1699 C
1700 C MadGraph StandAlone Version
1701 C
1702 C Returns amplitude squared summed/avg over colors
1703 C and helicities
1704 C for the point in phase space P(0:3,NEXTERNAL)
1705 C
1706 C Process: c s > u s ta+ ta- / h QED=4
1707 C
1708  IMPLICIT NONE
1709 C
1710 C CONSTANTS
1711 C
1712  INTEGER nexternal
1713  parameter(nexternal=6)
1714  INTEGER ncomb
1715  parameter( ncomb=64)
1716 C
1717 C ARGUMENTS
1718 C
1719  REAL*8 p(0:3,nexternal),ans
1720  INTEGER h1,h2
1721 C
1722 C LOCAL VARIABLES
1723 C
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/
1797 C ----------
1798 C BEGIN CODE
1799 C ----------
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)
1813 C
1814 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1815 C By the MadGraph Development Team
1816 C Please visit us at https://launchpad.net/madgraph5
1817 C
1818 C Returns amplitude squared summed/avg over colors
1819 C for the point with external lines W(0:6,NEXTERNAL)
1820 C
1821 C Process: c s > u s ta+ ta- / h QED=4
1822 C
1823  IMPLICIT NONE
1824 C
1825 C CONSTANTS
1826 C
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))
1837 C
1838 C ARGUMENTS
1839 C
1840  REAL*8 p(0:3,nexternal)
1841  INTEGER nhel(nexternal), ic(nexternal)
1842 C
1843 C LOCAL VARIABLES
1844 C
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)/
1852 C
1853 C GLOBAL VARIABLES
1854 C
1855  include 'coupl.inc'
1856 C
1857 C COLOR DATA
1858 C
1859  DATA denom(1)/1/
1860  DATA (cf(i, 1),i= 1, 1) / 9/
1861 C 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 
1870 C ----------
1871 C BEGIN CODE
1872 C ----------
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))
1882 C 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))
1885 C 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))
1888 C 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))
1891 C 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))
1894 C Amplitude(s) for diagram number 5
1895  CALL ffv1_0(w(1,2),w(1,5),w(1,9),gc_1,amp(5))
1896 C 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))
1898 C 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))
1901 C 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))
1904 C Amplitude(s) for diagram number 9
1905  CALL ffv2_0(w(1,2),w(1,4),w(1,8),gc_100,amp(9))
1906 C 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))
1909 C 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 
1928 C ------- end subprocesses CS->jjtautau, no Higgs
1929 
1930 
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7