C++InterfacetoTauola
SCX.f
1 c --------------------- subprocesses initiated by SCbar
2 
3  SUBROUTINE scx(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,1 OR 3 I4=0,-2 OR -4
10 
11  ans=0.d0
12 
13 
14  IF(key.EQ.1) THEN
15  IF(i3.EQ.1 .AND. i4.EQ.-2) CALL scx_dux_h(p,h1,h2,ans)
16  IF(i3.EQ.3 .AND. i4.EQ.-2) CALL scx_sux_h(p,h1,h2,ans)
17  IF(i3.EQ.1 .AND. i4.EQ.-4) CAll scx_dcx_h(p,h1,h2,ans)
18  IF(i3.EQ.3 .AND. i4.EQ.-4) CALL scx_scx_h(p,h1,h2,ans)
19  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
20  CALL scx_dux_h(p,h1,h2,ans1)
21  CALL scx_sux_h(p,h1,h2,ans2)
22  CALL scx_dcx_h(p,h1,h2,ans3)
23  CALL scx_scx_h(p,h1,h2,ans4)
24  ans=ans1+ans2+ans3+ans4
25  ENDIF
26  ELSE IF(key.EQ.0) THEN
27  IF(i3.EQ.1 .AND. i4.EQ.-2) CALL scx_dux_noh(p,h1,h2,ans)
28  IF(i3.EQ.3 .AND. i4.EQ.-2) CALL scx_sux_noh(p,h1,h2,ans)
29  IF(i3.EQ.1 .AND. i4.EQ.-4) CAll scx_dcx_noh(p,h1,h2,ans)
30  IF(i3.EQ.3 .AND. i4.EQ.-4) CALL scx_scx_noh(p,h1,h2,ans)
31  IF(i3.EQ.0 .AND. i4.EQ.0) THEN
32  CALL scx_dux_noh(p,h1,h2,ans1)
33  CALL scx_sux_noh(p,h1,h2,ans2)
34  CALL scx_dcx_noh(p,h1,h2,ans3)
35  CALL scx_scx_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 SCX
43 
44 
45 
46 C ----- begin subprocesses SCX->DUX with Higgs->tautau
47 
48 
49  SUBROUTINE scx_dux_h(P,H1,H2,ANS)
50 C
51 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
52 C By the MadGraph Development Team
53 C Please visit us at https://launchpad.net/madgraph5
54 C
55 C MadGraph StandAlone Version
56 C
57 C Returns amplitude squared summed/avg over colors
58 C and helicities
59 C for the point in phase space P(0:3,NEXTERNAL)
60 C
61 C Process: s c~ > d u~ h WEIGHTED=6
62 C * Decay: h > ta+ ta- WEIGHTED=2
63 C
64  IMPLICIT NONE
65 C
66 C CONSTANTS
67 C
68  INTEGER nexternal
69  parameter(nexternal=6)
70  INTEGER ncomb
71  parameter( ncomb=64)
72 C
73 C ARGUMENTS
74 C
75  REAL*8 p(0:3,nexternal),ans
76  INTEGER h1,h2
77 C
78 C LOCAL VARIABLES
79 C
80  INTEGER nhel(nexternal,ncomb),ntry
81  REAL*8 t
82  REAL*8 matrix_scx_dux_h
83  INTEGER ihel,iden, i
84  INTEGER jc(nexternal)
85  LOGICAL goodhel(ncomb)
86  DATA ntry/0/
87  DATA goodhel/ncomb*.false./
88  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
89  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
90  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
91  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
92  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
93  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
94  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
95  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
96  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
97  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
98  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
99  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
100  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
101  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
102  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
103  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
104  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
105  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
106  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
107  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
108  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
109  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
110  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
111  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
112  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
113  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
114  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
115  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
116  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
117  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
118  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
119  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
120  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
121  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
122  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
123  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
124  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
125  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
126  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
127  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
128  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
129  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
130  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
131  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
132  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
133  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
134  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
135  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
136  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
137  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
138  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
139  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
140  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
141  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
142  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
143  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
144  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
145  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
146  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
147  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
148  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
149  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
150  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
151  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
152  DATA iden/36/
153 C ----------
154 C BEGIN CODE
155 C ----------
156  DO ihel=1,nexternal
157  jc(ihel) = +1
158  ENDDO
159  ans = 0d0
160  DO ihel=1,ncomb
161  t=matrix_scx_dux_h(p ,h1,h2,nhel(1,ihel),jc(1))
162  ans=ans+t
163  ENDDO
164  ans=ans/dble(iden)
165  END
166 
167 
168  REAL*8 FUNCTION matrix_scx_dux_h(P,H1,H2,NHEL,IC)
169 C
170 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
171 C By the MadGraph Development Team
172 C Please visit us at https://launchpad.net/madgraph5
173 C
174 C Returns amplitude squared summed/avg over colors
175 C for the point with external lines W(0:6,NEXTERNAL)
176 C
177 C Process: s c~ > d u~ h WEIGHTED=6
178 C * Decay: h > ta+ ta- WEIGHTED=2
179 C
180  IMPLICIT NONE
181 C
182 C CONSTANTS
183 C
184  INTEGER ngraphs
185  parameter(ngraphs=1)
186  INTEGER nexternal
187  parameter(nexternal=6)
188  INTEGER nwavefuncs, ncolor
189  parameter(nwavefuncs=7, ncolor=1)
190  REAL*8 zero
191  parameter(zero=0d0)
192  COMPLEX*16 imag1
193  parameter(imag1=(0d0,1d0))
194 C
195 C ARGUMENTS
196 C
197  REAL*8 p(0:3,nexternal)
198  INTEGER nhel(nexternal), ic(nexternal)
199 C
200 C LOCAL VARIABLES
201 C
202  INTEGER i,j
203  COMPLEX*16 ztemp
204  REAL*8 denom(ncolor), cf(ncolor,ncolor)
205  COMPLEX*16 amp(ngraphs), jamp(ncolor)
206  COMPLEX*16 w(18,nwavefuncs)
207  COMPLEX*16 dum0,dum1
208  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
209 C
210 C GLOBAL VARIABLES
211 C
212  include 'coupl.inc'
213 C
214 C COLOR DATA
215 C
216  DATA denom(1)/1/
217  DATA (cf(i, 1),i= 1, 1) / 9/
218 C 1 T(2,1) T(3,4)
219 
220  INTEGER h1,h2
221  REAL*8 matrix
222  matrix_scx_dux_h=0.d0
223  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
224  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
225 
226 C ----------
227 C BEGIN CODE
228 C ----------
229  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
230  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
231  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
232  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
233  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
234  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
235  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
236  CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
237  CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
238 C Amplitude(s) for diagram number 1
239  CALL vvs1_0(w(1,2),w(1,6),w(1,7),gc_72,amp(1))
240  jamp(1)=-amp(1)
241 
242  matrix = 0.d0
243  DO i = 1, ncolor
244  ztemp = (0.d0,0.d0)
245  DO j = 1, ncolor
246  ztemp = ztemp + cf(j,i)*jamp(j)
247  ENDDO
248  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
249  ENDDO
250  matrix_scx_dux_h=matrix
251  ENDIF
252  ENDIF
253  END
254 
255 C ----- begin subprocesses SCX->SCX with Higgs->tautau
256 
257  SUBROUTINE scx_scx_h(P,H1,H2,ANS)
258 C
259 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
260 C By the MadGraph Development Team
261 C Please visit us at https://launchpad.net/madgraph5
262 C
263 C MadGraph StandAlone Version
264 C
265 C Returns amplitude squared summed/avg over colors
266 C and helicities
267 C for the point in phase space P(0:3,NEXTERNAL)
268 C
269 C Process: s c~ > s c~ h WEIGHTED=6
270 C * Decay: h > ta+ ta- WEIGHTED=2
271 C
272  IMPLICIT NONE
273 C
274 C CONSTANTS
275 C
276  INTEGER nexternal
277  parameter(nexternal=6)
278  INTEGER ncomb
279  parameter( ncomb=64)
280 C
281 C ARGUMENTS
282 C
283  REAL*8 p(0:3,nexternal),ans
284  INTEGER h1,h2
285 C
286 C LOCAL VARIABLES
287 C
288  INTEGER nhel(nexternal,ncomb),ntry
289  REAL*8 t
290  REAL*8 matrix_scx_scx_h
291  INTEGER ihel,iden, i
292  INTEGER jc(nexternal)
293  LOGICAL goodhel(ncomb)
294  DATA ntry/0/
295  DATA goodhel/ncomb*.false./
296  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
297  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
298  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
299  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
300  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
301  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
302  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
303  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
304  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
305  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
306  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
307  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
308  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
309  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
310  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
311  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
312  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
313  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
314  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
315  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
316  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
317  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
318  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
319  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
320  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
321  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
322  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
323  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
324  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
325  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
326  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
327  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
328  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
329  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
330  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
331  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
332  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
333  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
334  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
335  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
336  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
337  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
338  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
339  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
340  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
341  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
342  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
343  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
344  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
345  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
346  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
347  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
348  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
349  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
350  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
351  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
352  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
353  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
354  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
355  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
356  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
357  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
358  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
359  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
360  DATA iden/36/
361 C ----------
362 C BEGIN CODE
363 C ----------
364  DO ihel=1,nexternal
365  jc(ihel) = +1
366  ENDDO
367  ans = 0d0
368  DO ihel=1,ncomb
369  t=matrix_scx_scx_h(p ,h1,h2,nhel(1,ihel),jc(1))
370  ans=ans+t
371  ENDDO
372  ans=ans/dble(iden)
373  END
374 
375 
376  REAL*8 FUNCTION matrix_scx_scx_h(P,H1,H2,NHEL,IC)
377 C
378 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
379 C By the MadGraph Development Team
380 C Please visit us at https://launchpad.net/madgraph5
381 C
382 C Returns amplitude squared summed/avg over colors
383 C for the point with external lines W(0:6,NEXTERNAL)
384 C
385 C Process: s c~ > s c~ h WEIGHTED=6
386 C * Decay: h > ta+ ta- WEIGHTED=2
387 C
388  IMPLICIT NONE
389 C
390 C CONSTANTS
391 C
392  INTEGER ngraphs
393  parameter(ngraphs=2)
394  INTEGER nexternal
395  parameter(nexternal=6)
396  INTEGER nwavefuncs, ncolor
397  parameter(nwavefuncs=7, ncolor=2)
398  REAL*8 zero
399  parameter(zero=0d0)
400  COMPLEX*16 imag1
401  parameter(imag1=(0d0,1d0))
402 C
403 C ARGUMENTS
404 C
405  REAL*8 p(0:3,nexternal)
406  INTEGER nhel(nexternal), ic(nexternal)
407 C
408 C LOCAL VARIABLES
409 C
410  INTEGER i,j
411  COMPLEX*16 ztemp
412  REAL*8 denom(ncolor), cf(ncolor,ncolor)
413  COMPLEX*16 amp(ngraphs), jamp(ncolor)
414  COMPLEX*16 w(18,nwavefuncs)
415  COMPLEX*16 dum0,dum1
416  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
417 C
418 C GLOBAL VARIABLES
419 C
420  include 'coupl.inc'
421 C
422 C COLOR DATA
423 C
424  DATA denom(1)/1/
425  DATA (cf(i, 1),i= 1, 2) / 9, 3/
426 C 1 T(2,1) T(3,4)
427  DATA denom(2)/1/
428  DATA (cf(i, 2),i= 1, 2) / 3, 9/
429 C 1 T(2,4) T(3,1)
430 
431  INTEGER h1,h2
432  REAL*8 matrix
433  matrix_scx_scx_h=0.d0
434  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
435  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
436 
437 C ----------
438 C BEGIN CODE
439 C ----------
440  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
441  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
442  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
443  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
444  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
445  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
446  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
447  CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
448  CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,5))
449 C Amplitude(s) for diagram number 1
450  CALL vvs1_0(w(1,5),w(1,6),w(1,7),gc_72,amp(1))
451  CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
452  CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,3))
453 C Amplitude(s) for diagram number 2
454  CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
455  jamp(1)=-amp(1)
456  jamp(2)=+amp(2)
457 
458  matrix = 0.d0
459  DO i = 1, ncolor
460  ztemp = (0.d0,0.d0)
461  DO j = 1, ncolor
462  ztemp = ztemp + cf(j,i)*jamp(j)
463  ENDDO
464  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
465  ENDDO
466  matrix_scx_scx_h=matrix
467  ENDIF
468  ENDIF
469  END
470 
471 C ----- begin subprocesses SCX->DCX with Higgs->tautau
472 
473  SUBROUTINE scx_dcx_h(P,H1,H2,ANS)
474 C
475 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
476 C By the MadGraph Development Team
477 C Please visit us at https://launchpad.net/madgraph5
478 C
479 C MadGraph StandAlone Version
480 C
481 C Returns amplitude squared summed/avg over colors
482 C and helicities
483 C for the point in phase space P(0:3,NEXTERNAL)
484 C
485 C Process: s c~ > d c~ h WEIGHTED=6
486 C * Decay: h > ta+ ta- WEIGHTED=2
487 C
488  IMPLICIT NONE
489 C
490 C CONSTANTS
491 C
492  INTEGER nexternal
493  parameter(nexternal=6)
494  INTEGER ncomb
495  parameter( ncomb=64)
496 C
497 C ARGUMENTS
498 C
499  REAL*8 p(0:3,nexternal),ans
500  INTEGER h1,h2
501 C
502 C LOCAL VARIABLES
503 C
504  INTEGER nhel(nexternal,ncomb),ntry
505  REAL*8 t
506  REAL*8 matrix_scx_dcx_h
507  INTEGER ihel,iden, i
508  INTEGER jc(nexternal)
509  LOGICAL goodhel(ncomb)
510  DATA ntry/0/
511  DATA goodhel/ncomb*.false./
512  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
513  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
514  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
515  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
516  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
517  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
518  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
519  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
520  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
521  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
522  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
523  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
524  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
525  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
526  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
527  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
528  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
529  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
530  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
531  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
532  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
533  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
534  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
535  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
536  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
537  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
538  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
539  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
540  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
541  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
542  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
543  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
544  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
545  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
546  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
547  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
548  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
549  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
550  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
551  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
552  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
553  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
554  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
555  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
556  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
557  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
558  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
559  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
560  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
561  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
562  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
563  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
564  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
565  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
566  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
567  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
568  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
569  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
570  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
571  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
572  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
573  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
574  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
575  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
576  DATA iden/36/
577 C ----------
578 C BEGIN CODE
579 C ----------
580  DO ihel=1,nexternal
581  jc(ihel) = +1
582  ENDDO
583  ans = 0d0
584  DO ihel=1,ncomb
585  t=matrix_scx_dcx_h(p ,h1,h2,nhel(1,ihel),jc(1))
586  ans=ans+t
587  ENDDO
588  ans=ans/dble(iden)
589  END
590 
591 
592  REAL*8 FUNCTION matrix_scx_dcx_h(P,H1,H2,NHEL,IC)
593 C
594 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
595 C By the MadGraph Development Team
596 C Please visit us at https://launchpad.net/madgraph5
597 C
598 C Returns amplitude squared summed/avg over colors
599 C for the point with external lines W(0:6,NEXTERNAL)
600 C
601 C Process: s c~ > d c~ h WEIGHTED=6
602 C * Decay: h > ta+ ta- WEIGHTED=2
603 C
604  IMPLICIT NONE
605 C
606 C CONSTANTS
607 C
608  INTEGER ngraphs
609  parameter(ngraphs=1)
610  INTEGER nexternal
611  parameter(nexternal=6)
612  INTEGER nwavefuncs, ncolor
613  parameter(nwavefuncs=7, ncolor=1)
614  REAL*8 zero
615  parameter(zero=0d0)
616  COMPLEX*16 imag1
617  parameter(imag1=(0d0,1d0))
618 C
619 C ARGUMENTS
620 C
621  REAL*8 p(0:3,nexternal)
622  INTEGER nhel(nexternal), ic(nexternal)
623 C
624 C LOCAL VARIABLES
625 C
626  INTEGER i,j
627  COMPLEX*16 ztemp
628  REAL*8 denom(ncolor), cf(ncolor,ncolor)
629  COMPLEX*16 amp(ngraphs), jamp(ncolor)
630  COMPLEX*16 w(18,nwavefuncs)
631  COMPLEX*16 dum0,dum1
632  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
633 C
634 C GLOBAL VARIABLES
635 C
636  include 'coupl.inc'
637 C
638 C COLOR DATA
639 C
640  DATA denom(1)/1/
641  DATA (cf(i, 1),i= 1, 1) / 9/
642 C 1 T(2,1) T(3,4)
643 
644  INTEGER h1,h2
645  REAL*8 matrix
646  matrix_scx_dcx_h=0.d0
647  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
648  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
649 
650 C ----------
651 C BEGIN CODE
652 C ----------
653  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
654  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
655  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
656  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
657  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
658  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
659  CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
660  CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
661  CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,2))
662 C Amplitude(s) for diagram number 1
663  CALL vvs1_0(w(1,2),w(1,6),w(1,7),gc_72,amp(1))
664  jamp(1)=-amp(1)
665 
666  matrix = 0.d0
667  DO i = 1, ncolor
668  ztemp = (0.d0,0.d0)
669  DO j = 1, ncolor
670  ztemp = ztemp + cf(j,i)*jamp(j)
671  ENDDO
672  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
673  ENDDO
674  matrix_scx_dcx_h=matrix
675  ENDIF
676  ENDIF
677  END
678 
679 C ----- begin subprocesses SCX->SUX with Higgs->tautau
680 
681  SUBROUTINE scx_sux_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: s c~ > s u~ 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_scx_sux_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_scx_sux_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_scx_sux_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: s c~ > s u~ 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_scx_sux_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_101,mw,ww,w(1,2))
870 C Amplitude(s) for diagram number 1
871  CALL vvs1_0(w(1,2),w(1,6),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 
883  matrix_scx_sux_h=matrix
884  ENDIF
885  ENDIF
886  END
887 
888 C -------------- NOW PROCESSIS WITHOUT HIGGS
889 
890 C ----- begin subprocesses SCX->DUX NO Higgs
891 
892  SUBROUTINE scx_dux_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: s c~ > d u~ 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_scx_dux_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_scx_dux_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_scx_dux_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: s c~ > d u~ 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_scx_dux_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_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1079 C Amplitude(s) for diagram number 1
1080  CALL ffv2_0(w(1,5),w(1,9),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,8),w(1,7),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,8),w(1,7),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_2,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_1,amp(5))
1093 C Amplitude(s) for diagram number 6
1094  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1095 C Amplitude(s) for diagram number 7
1096  CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,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_2,amp(8))
1100  CALL ffv1_2(w(1,1),w(1,9),gc_1,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_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
1105  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,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_scx_dux_noh=matrix
1120  ENDIF
1121  ENDIF
1122  END
1123 
1124 C ----- begin subprocesses SCX->SCX NO Higgs
1125 
1126  SUBROUTINE scx_scx_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: s c~ > s c~ 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_scx_scx_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_scx_scx_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_scx_scx_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: s c~ > s c~ 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_scx_scx_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_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1316 C Amplitude(s) for diagram number 1
1317  CALL ffv2_0(w(1,5),w(1,9),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,8),w(1,7),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,8),w(1,7),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_2,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_1,amp(5))
1330 C Amplitude(s) for diagram number 6
1331  CALL ffv2_5_0(w(1,4),w(1,11),w(1,10),gc_51,gc_58,amp(6))
1332 C Amplitude(s) for diagram number 7
1333  CALL ffv2_3_0(w(1,12),w(1,3),w(1,10),gc_50,gc_58,amp(7))
1334  CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,12))
1335  CALL ffv1p0_3(w(1,4),w(1,2),gc_2,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_5_3(w(1,4),w(1,2),gc_51,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_3_3(w(1,1),w(1,3),gc_50,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_2,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_2,amp(16))
1361  CALL ffv1_2(w(1,4),w(1,12),gc_2,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_2,amp(17))
1364 C Amplitude(s) for diagram number 18
1365  CALL ffv2_5_0(w(1,4),w(1,15),w(1,10),gc_51,gc_58,amp(18))
1366 C Amplitude(s) for diagram number 19
1367  CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,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_2,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_2,amp(21))
1375 C Amplitude(s) for diagram number 22
1376  CALL ffv2_5_0(w(1,4),w(1,15),w(1,10),gc_51,gc_58,amp(22))
1377 C Amplitude(s) for diagram number 23
1378  CALL ffv2_5_0(w(1,12),w(1,2),w(1,10),gc_51,gc_58,amp(23))
1379  CALL ffv2_5_1(w(1,2),w(1,13),gc_51,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_2,amp(24))
1382  CALL ffv2_5_2(w(1,4),w(1,13),gc_51,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_2,amp(25))
1385 C Amplitude(s) for diagram number 26
1386  CALL ffv2_5_0(w(1,4),w(1,12),w(1,10),gc_51,gc_58,amp(26))
1387 C Amplitude(s) for diagram number 27
1388  CALL ffv2_5_0(w(1,15),w(1,2),w(1,10),gc_51,gc_58,amp(27))
1389  CALL ffv1_2(w(1,1),w(1,11),gc_1,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_1,amp(28))
1392  CALL ffv1_2(w(1,1),w(1,9),gc_1,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_1,amp(29))
1395 C Amplitude(s) for diagram number 30
1396  CALL ffv2_3_0(w(1,15),w(1,3),w(1,10),gc_50,gc_58,amp(30))
1397  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,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_1,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_1,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_3_0(w(1,4),w(1,3),w(1,10),gc_50,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_3_2(w(1,1),w(1,14),gc_50,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_1,amp(36))
1413 C Amplitude(s) for diagram number 37
1414  CALL ffv2_3_0(w(1,12),w(1,3),w(1,14),gc_50,gc_58,amp(37))
1415 C Amplitude(s) for diagram number 38
1416  CALL ffv2_3_0(w(1,11),w(1,3),w(1,10),gc_50,gc_58,amp(38))
1417 C Amplitude(s) for diagram number 39
1418  CALL ffv2_3_0(w(1,15),w(1,3),w(1,14),gc_50,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_2,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_5_0(w(1,14),w(1,2),w(1,10),gc_51,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_scx_scx_noh=matrix
1448  ENDIF
1449  ENDIF
1450  END
1451 
1452 C ----- begin subprocesses SCX->DCX NO Higgs
1453 
1454  SUBROUTINE scx_dcx_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: s c~ > d c~ 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_scx_dcx_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_scx_dcx_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_scx_dcx_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: s c~ > d c~ 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_scx_dcx_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_44,mw,ww,w(1,8))
1640  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1641 C Amplitude(s) for diagram number 1
1642  CALL ffv2_0(w(1,5),w(1,9),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,8),w(1,7),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,8),w(1,7),w(1,10),gc_53,amp(3))
1649  CALL ffv2_1(w(1,3),w(1,7),gc_44,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_2,amp(4))
1652  CALL ffv2_2(w(1,4),w(1,7),gc_44,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_1,amp(5))
1655 C Amplitude(s) for diagram number 6
1656  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1657 C Amplitude(s) for diagram number 7
1658  CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,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_2,amp(8))
1662  CALL ffv1_2(w(1,1),w(1,9),gc_1,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_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
1667  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,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_scx_dcx_noh=matrix
1682  ENDIF
1683  ENDIF
1684  END
1685 
1686 C ----- begin subprocesses SCX->SUX NO Higgs
1687 
1688  SUBROUTINE scx_sux_noh(P,H1,H2,ANS)
1689 C
1690 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1691 C By the MadGraph Development Team
1692 C Please visit us at https://launchpad.net/madgraph5
1693 C
1694 C MadGraph StandAlone Version
1695 C
1696 C Returns amplitude squared summed/avg over colors
1697 C and helicities
1698 C for the point in phase space P(0:3,NEXTERNAL)
1699 C
1700 C Process: s c~ > s u~ ta+ ta- / h QED=4
1701 C
1702  IMPLICIT NONE
1703 C
1704 C CONSTANTS
1705 C
1706  INTEGER nexternal
1707  parameter(nexternal=6)
1708  INTEGER ncomb
1709  parameter( ncomb=64)
1710 C
1711 C ARGUMENTS
1712 C
1713  REAL*8 p(0:3,nexternal),ans
1714  INTEGER h1,h2
1715 C
1716 C LOCAL VARIABLES
1717 C
1718  INTEGER nhel(nexternal,ncomb),ntry
1719  REAL*8 t
1720  REAL*8 matrix_scx_sux_noh
1721  INTEGER ihel,iden, i
1722  INTEGER jc(nexternal)
1723  LOGICAL goodhel(ncomb)
1724  DATA ntry/0/
1725  DATA goodhel/ncomb*.false./
1726  DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1727  DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1728  DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1729  DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1730  DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1731  DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1732  DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1733  DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1734  DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1735  DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1736  DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1737  DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1738  DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1739  DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1740  DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1741  DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1742  DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1743  DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1744  DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1745  DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1746  DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1747  DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1748  DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1749  DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1750  DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1751  DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1752  DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1753  DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1754  DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1755  DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1756  DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1757  DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1758  DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1759  DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1760  DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1761  DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1762  DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1763  DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1764  DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1765  DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1766  DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1767  DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1768  DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1769  DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1770  DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1771  DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1772  DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1773  DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1774  DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1775  DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1776  DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1777  DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1778  DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1779  DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1780  DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1781  DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1782  DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1783  DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1784  DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1785  DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1786  DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1787  DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1788  DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1789  DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1790  DATA iden/36/
1791 C ----------
1792 C BEGIN CODE
1793 C ----------
1794  DO ihel=1,nexternal
1795  jc(ihel) = +1
1796  ENDDO
1797  ans = 0d0
1798  DO ihel=1,ncomb
1799  t=matrix_scx_sux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1800  ans=ans+t
1801  ENDDO
1802  ans=ans/dble(iden)
1803  END
1804 
1805 
1806  REAL*8 FUNCTION matrix_scx_sux_noh(P,H1,H2,NHEL,IC)
1807 C
1808 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1809 C By the MadGraph Development Team
1810 C Please visit us at https://launchpad.net/madgraph5
1811 C
1812 C Returns amplitude squared summed/avg over colors
1813 C for the point with external lines W(0:6,NEXTERNAL)
1814 C
1815 C Process: s c~ > s u~ ta+ ta- / h QED=4
1816 C
1817  IMPLICIT NONE
1818 C
1819 C CONSTANTS
1820 C
1821  INTEGER ngraphs
1822  parameter(ngraphs=11)
1823  INTEGER nexternal
1824  parameter(nexternal=6)
1825  INTEGER nwavefuncs, ncolor
1826  parameter(nwavefuncs=10, ncolor=1)
1827  REAL*8 zero
1828  parameter(zero=0d0)
1829  COMPLEX*16 imag1
1830  parameter(imag1=(0d0,1d0))
1831 C
1832 C ARGUMENTS
1833 C
1834  REAL*8 p(0:3,nexternal)
1835  INTEGER nhel(nexternal), ic(nexternal)
1836 C
1837 C LOCAL VARIABLES
1838 C
1839  INTEGER i,j
1840  COMPLEX*16 ztemp
1841  REAL*8 denom(ncolor), cf(ncolor,ncolor)
1842  COMPLEX*16 amp(ngraphs), jamp(ncolor)
1843  COMPLEX*16 w(18,nwavefuncs)
1844  COMPLEX*16 dum0,dum1
1845  DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1846 C
1847 C GLOBAL VARIABLES
1848 C
1849  include 'coupl.inc'
1850 C
1851 C COLOR DATA
1852 C
1853  DATA denom(1)/1/
1854  DATA (cf(i, 1),i= 1, 1) / 9/
1855 C 1 T(2,1) T(3,4)
1856 
1857  INTEGER h1,h2
1858  REAL*8 matrix
1859  matrix_scx_sux_noh=0.d0
1860  IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
1861  IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
1862 
1863 C ----------
1864 C BEGIN CODE
1865 C ----------
1866  CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1867  CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1868  CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1869  CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1870  CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1871  CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1872  CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1873  CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1874  CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1875 C Amplitude(s) for diagram number 1
1876  CALL ffv2_0(w(1,5),w(1,9),w(1,8),gc_108,amp(1))
1877  CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1878 C Amplitude(s) for diagram number 2
1879  CALL vvv1_0(w(1,9),w(1,8),w(1,7),gc_4,amp(2))
1880  CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1881 C Amplitude(s) for diagram number 3
1882  CALL vvv1_0(w(1,8),w(1,7),w(1,10),gc_53,amp(3))
1883  CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,6))
1884 C Amplitude(s) for diagram number 4
1885  CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_2,amp(4))
1886  CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
1887 C Amplitude(s) for diagram number 5
1888  CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1889 C Amplitude(s) for diagram number 6
1890  CALL ffv2_5_0(w(1,4),w(1,6),w(1,10),gc_51,gc_58,amp(6))
1891 C Amplitude(s) for diagram number 7
1892  CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(7))
1893  CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1894 C Amplitude(s) for diagram number 8
1895  CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(8))
1896  CALL ffv1_2(w(1,1),w(1,9),gc_1,zero,zero,w(1,3))
1897 C Amplitude(s) for diagram number 9
1898  CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1899 C Amplitude(s) for diagram number 10
1900  CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(10))
1901  CALL ffv2_3_2(w(1,1),w(1,10),gc_50,gc_58,zero,zero,w(1,5))
1902 C Amplitude(s) for diagram number 11
1903  CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_100,amp(11))
1904  jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1905  $ -amp(9)-amp(10)-amp(11)
1906 
1907  matrix = 0.d0
1908  DO i = 1, ncolor
1909  ztemp = (0.d0,0.d0)
1910  DO j = 1, ncolor
1911  ztemp = ztemp + cf(j,i)*jamp(j)
1912  ENDDO
1913  matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1914  ENDDO
1915  matrix_scx_sux_noh=matrix
1916  ENDIF
1917  ENDIF
1918  END
1919 
1920 
1921 
1922 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