C++ Interface to Tauola
SCX.f
1c --------------------- 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
9C 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
46C ----- begin subprocesses SCX->DUX with Higgs->tautau
47
48
49 SUBROUTINE scx_dux_h(P,H1,H2,ANS)
50C
51C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
52C By the MadGraph Development Team
53C Please visit us at https://launchpad.net/madgraph5
54C
55C MadGraph StandAlone Version
56C
57C Returns amplitude squared summed/avg over colors
58C and helicities
59C for the point in phase space P(0:3,NEXTERNAL)
60C
61C Process: s c~ > d u~ h WEIGHTED=6
62C * Decay: h > ta+ ta- WEIGHTED=2
63C
64 IMPLICIT NONE
65C
66C CONSTANTS
67C
68 INTEGER NEXTERNAL
69 parameter(nexternal=6)
70 INTEGER NCOMB
71 parameter( ncomb=64)
72C
73C ARGUMENTS
74C
75 real*8 p(0:3,nexternal),ans
76 INTEGER H1,H2
77C
78C LOCAL VARIABLES
79C
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/
153C ----------
154C BEGIN CODE
155C ----------
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)
169C
170C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
171C By the MadGraph Development Team
172C Please visit us at https://launchpad.net/madgraph5
173C
174C Returns amplitude squared summed/avg over colors
175C for the point with external lines W(0:6,NEXTERNAL)
176C
177C Process: s c~ > d u~ h WEIGHTED=6
178C * Decay: h > ta+ ta- WEIGHTED=2
179C
180 IMPLICIT NONE
181C
182C CONSTANTS
183C
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))
194C
195C ARGUMENTS
196C
197 real*8 p(0:3,nexternal)
198 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
199C
200C LOCAL VARIABLES
201C
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)/
209C
210C GLOBAL VARIABLES
211C
212 include 'coupl.inc'
213C
214C COLOR DATA
215C
216 DATA denom(1)/1/
217 DATA (cf(i, 1),i= 1, 1) / 9/
218C 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
226C ----------
227C BEGIN CODE
228C ----------
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))
238C 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
255C ----- begin subprocesses SCX->SCX with Higgs->tautau
256
257 SUBROUTINE scx_scx_h(P,H1,H2,ANS)
258C
259C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
260C By the MadGraph Development Team
261C Please visit us at https://launchpad.net/madgraph5
262C
263C MadGraph StandAlone Version
264C
265C Returns amplitude squared summed/avg over colors
266C and helicities
267C for the point in phase space P(0:3,NEXTERNAL)
268C
269C Process: s c~ > s c~ h WEIGHTED=6
270C * Decay: h > ta+ ta- WEIGHTED=2
271C
272 IMPLICIT NONE
273C
274C CONSTANTS
275C
276 INTEGER NEXTERNAL
277 parameter(nexternal=6)
278 INTEGER NCOMB
279 parameter( ncomb=64)
280C
281C ARGUMENTS
282C
283 real*8 p(0:3,nexternal),ans
284 INTEGER H1,H2
285C
286C LOCAL VARIABLES
287C
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/
361C ----------
362C BEGIN CODE
363C ----------
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)
377C
378C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
379C By the MadGraph Development Team
380C Please visit us at https://launchpad.net/madgraph5
381C
382C Returns amplitude squared summed/avg over colors
383C for the point with external lines W(0:6,NEXTERNAL)
384C
385C Process: s c~ > s c~ h WEIGHTED=6
386C * Decay: h > ta+ ta- WEIGHTED=2
387C
388 IMPLICIT NONE
389C
390C CONSTANTS
391C
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))
402C
403C ARGUMENTS
404C
405 real*8 p(0:3,nexternal)
406 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
407C
408C LOCAL VARIABLES
409C
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)/
417C
418C GLOBAL VARIABLES
419C
420 include 'coupl.inc'
421C
422C COLOR DATA
423C
424 DATA denom(1)/1/
425 DATA (cf(i, 1),i= 1, 2) / 9, 3/
426C 1 T(2,1) T(3,4)
427 DATA denom(2)/1/
428 DATA (cf(i, 2),i= 1, 2) / 3, 9/
429C 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
437C ----------
438C BEGIN CODE
439C ----------
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))
449C 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))
453C 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
471C ----- begin subprocesses SCX->DCX with Higgs->tautau
472
473 SUBROUTINE scx_dcx_h(P,H1,H2,ANS)
474C
475C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
476C By the MadGraph Development Team
477C Please visit us at https://launchpad.net/madgraph5
478C
479C MadGraph StandAlone Version
480C
481C Returns amplitude squared summed/avg over colors
482C and helicities
483C for the point in phase space P(0:3,NEXTERNAL)
484C
485C Process: s c~ > d c~ h WEIGHTED=6
486C * Decay: h > ta+ ta- WEIGHTED=2
487C
488 IMPLICIT NONE
489C
490C CONSTANTS
491C
492 INTEGER NEXTERNAL
493 parameter(nexternal=6)
494 INTEGER NCOMB
495 parameter( ncomb=64)
496C
497C ARGUMENTS
498C
499 real*8 p(0:3,nexternal),ans
500 INTEGER H1,H2
501C
502C LOCAL VARIABLES
503C
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/
577C ----------
578C BEGIN CODE
579C ----------
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)
593C
594C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
595C By the MadGraph Development Team
596C Please visit us at https://launchpad.net/madgraph5
597C
598C Returns amplitude squared summed/avg over colors
599C for the point with external lines W(0:6,NEXTERNAL)
600C
601C Process: s c~ > d c~ h WEIGHTED=6
602C * Decay: h > ta+ ta- WEIGHTED=2
603C
604 IMPLICIT NONE
605C
606C CONSTANTS
607C
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))
618C
619C ARGUMENTS
620C
621 real*8 p(0:3,nexternal)
622 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
623C
624C LOCAL VARIABLES
625C
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)/
633C
634C GLOBAL VARIABLES
635C
636 include 'coupl.inc'
637C
638C COLOR DATA
639C
640 DATA denom(1)/1/
641 DATA (cf(i, 1),i= 1, 1) / 9/
642C 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
650C ----------
651C BEGIN CODE
652C ----------
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))
662C 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
679C ----- begin subprocesses SCX->SUX with Higgs->tautau
680
681 SUBROUTINE scx_sux_h(P,H1,H2,ANS)
682C
683C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
684C By the MadGraph Development Team
685C Please visit us at https://launchpad.net/madgraph5
686C
687C MadGraph StandAlone Version
688C
689C Returns amplitude squared summed/avg over colors
690C and helicities
691C for the point in phase space P(0:3,NEXTERNAL)
692C
693C Process: s c~ > s u~ h WEIGHTED=6
694C * Decay: h > ta+ ta- WEIGHTED=2
695C
696 IMPLICIT NONE
697C
698C CONSTANTS
699C
700 INTEGER NEXTERNAL
701 parameter(nexternal=6)
702 INTEGER NCOMB
703 parameter( ncomb=64)
704C
705C ARGUMENTS
706C
707 real*8 p(0:3,nexternal),ans
708 INTEGER H1,H2
709C
710C LOCAL VARIABLES
711C
712 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
713 real*8 t
714 real*8 matrix_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/
785C ----------
786C BEGIN CODE
787C ----------
788 DO ihel=1,nexternal
789 jc(ihel) = +1
790 ENDDO
791 ans = 0d0
792 DO ihel=1,ncomb
793 t=matrix_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)
801C
802C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
803C By the MadGraph Development Team
804C Please visit us at https://launchpad.net/madgraph5
805C
806C Returns amplitude squared summed/avg over colors
807C for the point with external lines W(0:6,NEXTERNAL)
808C
809C Process: s c~ > s u~ h WEIGHTED=6
810C * Decay: h > ta+ ta- WEIGHTED=2
811C
812 IMPLICIT NONE
813C
814C CONSTANTS
815C
816 INTEGER NGRAPHS
817 parameter(ngraphs=1)
818 INTEGER NEXTERNAL
819 parameter(nexternal=6)
820 INTEGER NWAVEFUNCS, NCOLOR
821 parameter(nwavefuncs=7, ncolor=1)
822 real*8 zero
823 parameter(zero=0d0)
824 COMPLEX*16 IMAG1
825 parameter(imag1=(0d0,1d0))
826C
827C ARGUMENTS
828C
829 real*8 p(0:3,nexternal)
830 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
831C
832C LOCAL VARIABLES
833C
834 INTEGER I,J
835 COMPLEX*16 ZTEMP
836 real*8 denom(ncolor), cf(ncolor,ncolor)
837 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
838 COMPLEX*16 W(18,NWAVEFUNCS)
839 COMPLEX*16 DUM0,DUM1
840 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
841C
842C GLOBAL VARIABLES
843C
844 include 'coupl.inc'
845C
846C COLOR DATA
847C
848 DATA denom(1)/1/
849 DATA (cf(i, 1),i= 1, 1) / 9/
850C 1 T(2,1) T(3,4)
851
852 INTEGER H1,H2
853 real*8 matrix
854 matrix_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
858C ----------
859C BEGIN CODE
860C ----------
861 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
862 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
863 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
864 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
865 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
866 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
867 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
868 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
869 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,2))
870C 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
888C -------------- NOW PROCESSIS WITHOUT HIGGS
889
890C ----- begin subprocesses SCX->DUX NO Higgs
891
892 SUBROUTINE scx_dux_noh(P,H1,H2,ANS)
893C
894C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
895C By the MadGraph Development Team
896C Please visit us at https://launchpad.net/madgraph5
897C
898C MadGraph StandAlone Version
899C
900C Returns amplitude squared summed/avg over colors
901C and helicities
902C for the point in phase space P(0:3,NEXTERNAL)
903C
904C Process: s c~ > d u~ ta+ ta- / h QED=4
905C
906 IMPLICIT NONE
907C
908C CONSTANTS
909C
910 INTEGER NEXTERNAL
911 parameter(nexternal=6)
912 INTEGER NCOMB
913 parameter( ncomb=64)
914C
915C ARGUMENTS
916C
917 real*8 p(0:3,nexternal),ans
918 INTEGER H1,H2
919C
920C LOCAL VARIABLES
921C
922 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
923 real*8 t
924 real*8 matrix_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/
995C ----------
996C BEGIN CODE
997C ----------
998 DO ihel=1,nexternal
999 jc(ihel) = +1
1000 ENDDO
1001 ans = 0d0
1002 DO ihel=1,ncomb
1003 t=matrix_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)
1011C
1012C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1013C By the MadGraph Development Team
1014C Please visit us at https://launchpad.net/madgraph5
1015C
1016C Returns amplitude squared summed/avg over colors
1017C for the point with external lines W(0:6,NEXTERNAL)
1018C
1019C Process: s c~ > d u~ ta+ ta- / h QED=4
1020C
1021 IMPLICIT NONE
1022C
1023C CONSTANTS
1024C
1025 INTEGER NGRAPHS
1026 parameter(ngraphs=11)
1027 INTEGER NEXTERNAL
1028 parameter(nexternal=6)
1029 INTEGER NWAVEFUNCS, NCOLOR
1030 parameter(nwavefuncs=10, ncolor=1)
1031 real*8 zero
1032 parameter(zero=0d0)
1033 COMPLEX*16 IMAG1
1034 parameter(imag1=(0d0,1d0))
1035C
1036C ARGUMENTS
1037C
1038 real*8 p(0:3,nexternal)
1039 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1040C
1041C LOCAL VARIABLES
1042C
1043 INTEGER I,J
1044 COMPLEX*16 ZTEMP
1045 real*8 denom(ncolor), cf(ncolor,ncolor)
1046 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1047 COMPLEX*16 W(18,NWAVEFUNCS)
1048 COMPLEX*16 DUM0,DUM1
1049 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1050C
1051C GLOBAL VARIABLES
1052C
1053 include 'coupl.inc'
1054C
1055C COLOR DATA
1056C
1057 DATA denom(1)/1/
1058 DATA (cf(i, 1),i= 1, 1) / 9/
1059C 1 T(2,1) T(3,4)
1060
1061 INTEGER H1,H2
1062 real*8 matrix
1063 matrix_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
1067C ----------
1068C BEGIN CODE
1069C ----------
1070 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1071 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1072 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1073 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1074 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1075 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1076 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1077 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
1078 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1079C 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))
1082C 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))
1085C 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))
1088C 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))
1091C Amplitude(s) for diagram number 5
1092 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1093C 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))
1095C 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))
1098C 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))
1101C Amplitude(s) for diagram number 9
1102 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1103C Amplitude(s) for diagram number 10
1104 CALL ffv2_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))
1106C Amplitude(s) for diagram number 11
1107 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_100,amp(11))
1108 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1109 $ -amp(9)-amp(10)-amp(11)
1110
1111 matrix = 0.d0
1112 DO i = 1, ncolor
1113 ztemp = (0.d0,0.d0)
1114 DO j = 1, ncolor
1115 ztemp = ztemp + cf(j,i)*jamp(j)
1116 ENDDO
1117 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1118 ENDDO
1119 matrix_scx_dux_noh=matrix
1120 ENDIF
1121 ENDIF
1122 END
1123
1124C ----- begin subprocesses SCX->SCX NO Higgs
1125
1126 SUBROUTINE scx_scx_noh(P,H1,H2,ANS)
1127C
1128C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1129C By the MadGraph Development Team
1130C Please visit us at https://launchpad.net/madgraph5
1131C
1132C MadGraph StandAlone Version
1133C
1134C Returns amplitude squared summed/avg over colors
1135C and helicities
1136C for the point in phase space P(0:3,NEXTERNAL)
1137C
1138C Process: s c~ > s c~ ta+ ta- / h QED=4
1139C
1140 IMPLICIT NONE
1141C
1142C CONSTANTS
1143C
1144 INTEGER NEXTERNAL
1145 parameter(nexternal=6)
1146 INTEGER NCOMB
1147 parameter( ncomb=64)
1148C
1149C ARGUMENTS
1150C
1151 real*8 p(0:3,nexternal),ans
1152 INTEGER H1,H2
1153C
1154C LOCAL VARIABLES
1155C
1156 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1157 real*8 t
1158 real*8 matrix_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/
1229C ----------
1230C BEGIN CODE
1231C ----------
1232 DO ihel=1,nexternal
1233 jc(ihel) = +1
1234 ENDDO
1235 ans = 0d0
1236 DO ihel=1,ncomb
1237 t=matrix_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)
1245C
1246C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1247C By the MadGraph Development Team
1248C Please visit us at https://launchpad.net/madgraph5
1249C
1250C Returns amplitude squared summed/avg over colors
1251C for the point with external lines W(0:6,NEXTERNAL)
1252C
1253C Process: s c~ > s c~ ta+ ta- / h QED=4
1254C
1255 IMPLICIT NONE
1256C
1257C CONSTANTS
1258C
1259 INTEGER NGRAPHS
1260 parameter(ngraphs=43)
1261 INTEGER NEXTERNAL
1262 parameter(nexternal=6)
1263 INTEGER NWAVEFUNCS, NCOLOR
1264 parameter(nwavefuncs=15, ncolor=2)
1265 real*8 zero
1266 parameter(zero=0d0)
1267 COMPLEX*16 IMAG1
1268 parameter(imag1=(0d0,1d0))
1269C
1270C ARGUMENTS
1271C
1272 real*8 p(0:3,nexternal)
1273 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1274C
1275C LOCAL VARIABLES
1276C
1277 INTEGER I,J
1278 COMPLEX*16 ZTEMP
1279 real*8 denom(ncolor), cf(ncolor,ncolor)
1280 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1281 COMPLEX*16 W(18,NWAVEFUNCS)
1282 COMPLEX*16 DUM0,DUM1
1283 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1284C
1285C GLOBAL VARIABLES
1286C
1287 include 'coupl.inc'
1288C
1289C COLOR DATA
1290C
1291 DATA denom(1)/1/
1292 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1293C 1 T(2,1) T(3,4)
1294 DATA denom(2)/1/
1295 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1296C 1 T(2,4) T(3,1)
1297
1298 INTEGER H1,H2
1299 real*8 matrix
1300 matrix_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
1304C ----------
1305C BEGIN CODE
1306C ----------
1307 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1308 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1309 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1310 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1311 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1312 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1313 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1314 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
1315 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1316C 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))
1319C 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))
1322C 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))
1325C 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))
1328C Amplitude(s) for diagram number 5
1329 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_1,amp(5))
1330C 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))
1332C 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))
1337C Amplitude(s) for diagram number 8
1338 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1339 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1340C Amplitude(s) for diagram number 9
1341 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1342 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,14))
1343C Amplitude(s) for diagram number 10
1344 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1345C Amplitude(s) for diagram number 11
1346 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1347 CALL ffv2_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))
1349C Amplitude(s) for diagram number 12
1350 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1351 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1352C Amplitude(s) for diagram number 13
1353 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1354C Amplitude(s) for diagram number 14
1355 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1356C Amplitude(s) for diagram number 15
1357 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1358 CALL ffv1_1(w(1,2),w(1,12),gc_2,zero,zero,w(1,15))
1359C Amplitude(s) for diagram number 16
1360 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_2,amp(16))
1361 CALL ffv1_2(w(1,4),w(1,12),gc_2,zero,zero,w(1,5))
1362C Amplitude(s) for diagram number 17
1363 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(17))
1364C 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))
1366C 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))
1370C 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))
1373C Amplitude(s) for diagram number 21
1374 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_2,amp(21))
1375C 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))
1377C 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))
1380C 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))
1383C Amplitude(s) for diagram number 25
1384 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_2,amp(25))
1385C 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))
1387C 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))
1390C 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))
1393C Amplitude(s) for diagram number 29
1394 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_1,amp(29))
1395C 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))
1398C 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))
1402C Amplitude(s) for diagram number 32
1403 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_1,amp(32))
1404C Amplitude(s) for diagram number 33
1405 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1406C Amplitude(s) for diagram number 34
1407 CALL ffv2_3_0(w(1,4),w(1,3),w(1,10),gc_50,gc_58,amp(34))
1408C Amplitude(s) for diagram number 35
1409 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1410 CALL ffv2_3_2(w(1,1),w(1,14),gc_50,gc_58,zero,zero,w(1,11))
1411C Amplitude(s) for diagram number 36
1412 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_1,amp(36))
1413C 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))
1415C 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))
1417C 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))
1420C Amplitude(s) for diagram number 40
1421 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_2,amp(40))
1422C Amplitude(s) for diagram number 41
1423 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_100,amp(41))
1424C Amplitude(s) for diagram number 42
1425 CALL ffv2_5_0(w(1,14),w(1,2),w(1,10),gc_51,gc_58,amp(42))
1426C Amplitude(s) for diagram number 43
1427 CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_100,amp(43))
1428 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1429 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1430 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1431 $ -amp(40)-amp(41)-amp(42)-amp(43)
1432 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1433 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1434 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1435 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1436 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1437 $ +amp(38)+amp(39)
1438
1439 matrix = 0.d0
1440 DO i = 1, ncolor
1441 ztemp = (0.d0,0.d0)
1442 DO j = 1, ncolor
1443 ztemp = ztemp + cf(j,i)*jamp(j)
1444 ENDDO
1445 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1446 ENDDO
1447 matrix_scx_scx_noh=matrix
1448 ENDIF
1449 ENDIF
1450 END
1451
1452C ----- begin subprocesses SCX->DCX NO Higgs
1453
1454 SUBROUTINE scx_dcx_noh(P,H1,H2,ANS)
1455C
1456C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1457C By the MadGraph Development Team
1458C Please visit us at https://launchpad.net/madgraph5
1459C
1460C MadGraph StandAlone Version
1461C
1462C Returns amplitude squared summed/avg over colors
1463C and helicities
1464C for the point in phase space P(0:3,NEXTERNAL)
1465C
1466C Process: s c~ > d c~ ta+ ta- / h QED=4
1467C
1468 IMPLICIT NONE
1469C
1470C CONSTANTS
1471C
1472 INTEGER NEXTERNAL
1473 parameter(nexternal=6)
1474 INTEGER NCOMB
1475 parameter( ncomb=64)
1476C
1477C ARGUMENTS
1478C
1479 real*8 p(0:3,nexternal),ans
1480 INTEGER H1,H2
1481C
1482C LOCAL VARIABLES
1483C
1484 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1485 real*8 t
1486 real*8 matrix_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/
1557C ----------
1558C BEGIN CODE
1559C ----------
1560 DO ihel=1,nexternal
1561 jc(ihel) = +1
1562 ENDDO
1563 ans = 0d0
1564 DO ihel=1,ncomb
1565 t=matrix_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)
1573C
1574C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1575C By the MadGraph Development Team
1576C Please visit us at https://launchpad.net/madgraph5
1577C
1578C Returns amplitude squared summed/avg over colors
1579C for the point with external lines W(0:6,NEXTERNAL)
1580C
1581C Process: s c~ > d c~ ta+ ta- / h QED=4
1582C
1583 IMPLICIT NONE
1584C
1585C CONSTANTS
1586C
1587 INTEGER NGRAPHS
1588 parameter(ngraphs=11)
1589 INTEGER NEXTERNAL
1590 parameter(nexternal=6)
1591 INTEGER NWAVEFUNCS, NCOLOR
1592 parameter(nwavefuncs=10, ncolor=1)
1593 real*8 zero
1594 parameter(zero=0d0)
1595 COMPLEX*16 IMAG1
1596 parameter(imag1=(0d0,1d0))
1597C
1598C ARGUMENTS
1599C
1600 real*8 p(0:3,nexternal)
1601 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1602C
1603C LOCAL VARIABLES
1604C
1605 INTEGER I,J
1606 COMPLEX*16 ZTEMP
1607 real*8 denom(ncolor), cf(ncolor,ncolor)
1608 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1609 COMPLEX*16 W(18,NWAVEFUNCS)
1610 COMPLEX*16 DUM0,DUM1
1611 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1612C
1613C GLOBAL VARIABLES
1614C
1615 include 'coupl.inc'
1616C
1617C COLOR DATA
1618C
1619 DATA denom(1)/1/
1620 DATA (cf(i, 1),i= 1, 1) / 9/
1621C 1 T(2,1) T(3,4)
1622
1623 INTEGER H1,H2
1624 real*8 matrix
1625 matrix_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
1629C ----------
1630C BEGIN CODE
1631C ----------
1632 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1633 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1634 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1635 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1636 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1637 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1638 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1639 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1640 CALL ffv2_1(w(1,6),w(1,7),gc_108,zero,zero,w(1,9))
1641C 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))
1644C 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))
1647C 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))
1650C 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))
1653C Amplitude(s) for diagram number 5
1654 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1655C 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))
1657C 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))
1660C 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))
1663C Amplitude(s) for diagram number 9
1664 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1665C Amplitude(s) for diagram number 10
1666 CALL ffv2_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))
1668C Amplitude(s) for diagram number 11
1669 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_100,amp(11))
1670 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1671 $ -amp(9)-amp(10)-amp(11)
1672
1673 matrix = 0.d0
1674 DO i = 1, ncolor
1675 ztemp = (0.d0,0.d0)
1676 DO j = 1, ncolor
1677 ztemp = ztemp + cf(j,i)*jamp(j)
1678 ENDDO
1679 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1680 ENDDO
1681 matrix_scx_dcx_noh=matrix
1682 ENDIF
1683 ENDIF
1684 END
1685
1686C ----- begin subprocesses SCX->SUX NO Higgs
1687
1688 SUBROUTINE scx_sux_noh(P,H1,H2,ANS)
1689C
1690C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1691C By the MadGraph Development Team
1692C Please visit us at https://launchpad.net/madgraph5
1693C
1694C MadGraph StandAlone Version
1695C
1696C Returns amplitude squared summed/avg over colors
1697C and helicities
1698C for the point in phase space P(0:3,NEXTERNAL)
1699C
1700C Process: s c~ > s u~ ta+ ta- / h QED=4
1701C
1702 IMPLICIT NONE
1703C
1704C CONSTANTS
1705C
1706 INTEGER NEXTERNAL
1707 parameter(nexternal=6)
1708 INTEGER NCOMB
1709 parameter( ncomb=64)
1710C
1711C ARGUMENTS
1712C
1713 real*8 p(0:3,nexternal),ans
1714 INTEGER H1,H2
1715C
1716C LOCAL VARIABLES
1717C
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/
1791C ----------
1792C BEGIN CODE
1793C ----------
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)
1807C
1808C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
1809C By the MadGraph Development Team
1810C Please visit us at https://launchpad.net/madgraph5
1811C
1812C Returns amplitude squared summed/avg over colors
1813C for the point with external lines W(0:6,NEXTERNAL)
1814C
1815C Process: s c~ > s u~ ta+ ta- / h QED=4
1816C
1817 IMPLICIT NONE
1818C
1819C CONSTANTS
1820C
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))
1831C
1832C ARGUMENTS
1833C
1834 real*8 p(0:3,nexternal)
1835 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1836C
1837C LOCAL VARIABLES
1838C
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)/
1846C
1847C GLOBAL VARIABLES
1848C
1849 include 'coupl.inc'
1850C
1851C COLOR DATA
1852C
1853 DATA denom(1)/1/
1854 DATA (cf(i, 1),i= 1, 1) / 9/
1855C 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
1863C ----------
1864C BEGIN CODE
1865C ----------
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))
1875C 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))
1878C 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))
1881C 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))
1884C 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))
1887C Amplitude(s) for diagram number 5
1888 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(5))
1889C 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))
1891C 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))
1894C 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))
1897C Amplitude(s) for diagram number 9
1898 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1899C 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))
1902C 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
1922C DONE 26.10.2015 STILL S=2 MISSING