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