C++ Interface to Tauola
CC.f
1C JK. corrected 26.10.2015
2c------------------begin subprocess initiated by CC----------
3
4 SUBROUTINE cc(P,I3,I4,H1,H2,KEY,ANS)
5 IMPLICIT NONE
6
7 INTEGER I3,I4, H1,H2, KEY
8 real*8 p(0:3,6), ans
9
10C I3=0,4 I4=0,4 ! only cc->cc
11
12 ans=0.d0
13
14 IF(key.EQ.1) THEN
15 IF(i3.EQ.4 .AND. i4.EQ.4) CALL cc_cc_h(p,h1,h2,ans)
16 IF(i3.EQ.0 .AND. i4.EQ.0) CALL cc_cc_h(p,h1,h2,ans)
17 ELSE IF(key.EQ.0) THEN
18 IF(i3.EQ.4 .AND. i4.EQ.4) CALL cc_cc_noh(p,h1,h2,ans)
19 IF(i3.EQ.0 .AND. i4.EQ.0) CALL cc_cc_noh(p,h1,h2,ans)
20 ELSE
21 WRITE(*,*) 'spin=2 NOT FINISHED'
22 stop
23 ENDIF
24 END ! SUBROUTINE CC
25
26C ---------begin subprocess CC->jjH with H-> tautau
27c ---------------------jj=uu only -------------
28
29 SUBROUTINE cc_cc_h(P,H1,H2,ANS)
30C
31C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
32C By the MadGraph Development Team
33C Please visit us at https://launchpad.net/madgraph5
34C
35C MadGraph StandAlone Version
36C
37C Returns amplitude squared summed/avg over colors
38C and helicities
39C for the point in phase space P(0:3,NEXTERNAL)
40C
41C Process: c c > c c h WEIGHTED=6
42C * Decay: h > ta+ ta- WEIGHTED=2
43C
44 IMPLICIT NONE
45C
46C CONSTANTS
47C
48 INTEGER NEXTERNAL
49 parameter(nexternal=6)
50 INTEGER NCOMB
51 parameter( ncomb=64)
52C
53C ARGUMENTS
54C
55 real*8 p(0:3,nexternal),ans
56 INTEGER H1,H2
57C
58C LOCAL VARIABLES
59C
60 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
61 real*8 t
62 real*8 matrix_cc_cc_h
63 INTEGER IHEL,IDEN, I
64 INTEGER JC(NEXTERNAL)
65 LOGICAL GOODHEL(NCOMB)
66 DATA ntry/0/
67 DATA goodhel/ncomb*.false./
68 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
69 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
70 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
71 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
72 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
73 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
74 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
75 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
76 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
77 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
78 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
79 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
80 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
81 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
82 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
83 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
84 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
85 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
86 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
87 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
88 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
89 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
90 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
91 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
92 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
93 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
94 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
95 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
96 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
97 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
98 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
99 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
100 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
101 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
102 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
103 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
104 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
105 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
106 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
107 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
108 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
109 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
110 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
111 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
112 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
113 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
114 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
115 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
116 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
117 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
118 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
119 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
120 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
121 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
122 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
123 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
124 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
125 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
126 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
127 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
128 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
129 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
130 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
131 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
132 DATA iden/72/
133C ----------
134C BEGIN CODE
135C ----------
136 DO ihel=1,nexternal
137 jc(ihel) = +1
138 ENDDO
139 ans = 0d0
140 DO ihel=1,ncomb
141 t=matrix_cc_cc_h(p ,h1,h2,nhel(1,ihel),jc(1))
142 ans=ans+t
143 ENDDO
144 ans=ans/dble(iden)
145 END
146
147
148 REAL*8 FUNCTION matrix_cc_cc_h(P,H1,H2,NHEL,IC)
149C
150C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
151C By the MadGraph Development Team
152C Please visit us at https://launchpad.net/madgraph5
153C
154C Returns amplitude squared summed/avg over colors
155C for the point with external lines W(0:6,NEXTERNAL)
156C
157C Process: c c > c c h WEIGHTED=6
158C * Decay: h > ta+ ta- WEIGHTED=2
159C
160 IMPLICIT NONE
161C
162C CONSTANTS
163C
164 INTEGER NGRAPHS
165 parameter(ngraphs=2)
166 INTEGER NEXTERNAL
167 parameter(nexternal=6)
168 INTEGER NWAVEFUNCS, NCOLOR
169 parameter(nwavefuncs=7, ncolor=2)
170 real*8 zero
171 parameter(zero=0d0)
172 COMPLEX*16 IMAG1
173 parameter(imag1=(0d0,1d0))
174C
175C ARGUMENTS
176C
177 real*8 p(0:3,nexternal)
178 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
179C
180C LOCAL VARIABLES
181C
182 INTEGER I,J
183 COMPLEX*16 ZTEMP
184 real*8 denom(ncolor), cf(ncolor,ncolor)
185 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
186 COMPLEX*16 W(18,NWAVEFUNCS)
187 COMPLEX*16 DUM0,DUM1
188 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
189C
190C GLOBAL VARIABLES
191C
192 include 'coupl.inc'
193C
194C COLOR DATA
195C
196 DATA denom(1)/1/
197 DATA (cf(i, 1),i= 1, 2) / 9, 3/
198C 1 T(3,1) T(4,2)
199 DATA denom(2)/1/
200 DATA (cf(i, 2),i= 1, 2) / 3, 9/
201C 1 T(3,2) T(4,1)
202
203 INTEGER H1,H2
204 real*8 matrix
205 matrix_cc_cc_h=0.d0
206 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
207 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
208
209
210C ----------
211C BEGIN CODE
212C ----------
213 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
214 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
215 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
216 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
217 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
218 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
219 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
220 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,6))
221 CALL ffv2_5_3(w(1,2),w(1,4),gc_51,gc_58,mz,wz,w(1,5))
222C Amplitude(s) for diagram number 1
223 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
224 CALL ffv2_5_3(w(1,1),w(1,4),gc_51,gc_58,mz,wz,w(1,5))
225 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,4))
226C Amplitude(s) for diagram number 2
227 CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_81,amp(2))
228 jamp(1)=-amp(1)
229 jamp(2)=+amp(2)
230
231 matrix = 0.d0
232 DO i = 1, ncolor
233 ztemp = (0.d0,0.d0)
234 DO j = 1, ncolor
235 ztemp = ztemp + cf(j,i)*jamp(j)
236 ENDDO
237 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
238 ENDDO
239 matrix_cc_cc_h=matrix
240 ENDIF
241 ENDIF
242 END
243C----------------end subprocess CC->jjH---------
244
245c----------------begin subprocess CC->jj_noH
246
247
248
249 SUBROUTINE cc_cc_noh(P,H1,H2,ANS)
250C
251C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
252C By the MadGraph Development Team
253C Please visit us at https://launchpad.net/madgraph5
254C
255C MadGraph StandAlone Version
256C
257C Returns amplitude squared summed/avg over colors
258C and helicities
259C for the point in phase space P(0:3,NEXTERNAL)
260C
261C Process: c c > c c ta+ ta- / h QED=4
262C
263 IMPLICIT NONE
264C
265C CONSTANTS
266C
267 INTEGER NEXTERNAL
268 parameter(nexternal=6)
269 INTEGER NCOMB
270 parameter( ncomb=64)
271C
272C ARGUMENTS
273C
274 real*8 p(0:3,nexternal),ans
275 INTEGER H1,H2
276C
277C LOCAL VARIABLES
278C
279 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
280 real*8 t
281 real*8 matrix_cc_cc_noh
282 INTEGER IHEL,IDEN, I
283 INTEGER JC(NEXTERNAL)
284 LOGICAL GOODHEL(NCOMB)
285 DATA ntry/0/
286 DATA goodhel/ncomb*.false./
287 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
288 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
289 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
290 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
291 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
292 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
293 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
294 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
295 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
296 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
297 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
298 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
299 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
300 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
301 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
302 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
303 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
304 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
305 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
306 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
307 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
308 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
309 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
310 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
311 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
312 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
313 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
314 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
315 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
316 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
317 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
318 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
319 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
320 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
321 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
322 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
323 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
324 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
325 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
326 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
327 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
328 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
329 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
330 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
331 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
332 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
333 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
334 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
335 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
336 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
337 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
338 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
339 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
340 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
341 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
342 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
343 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
344 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
345 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
346 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
347 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
348 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
349 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
350 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
351 DATA iden/72/
352C ----------
353C BEGIN CODE
354C ----------
355 DO ihel=1,nexternal
356 jc(ihel) = +1
357 ENDDO
358 ans = 0d0
359 DO ihel=1,ncomb
360 t=matrix_cc_cc_noh(p ,h1,h2,nhel(1,ihel),jc(1))
361 ans=ans+t
362 ENDDO
363 ans=ans/dble(iden)
364 END
365
366
367 REAL*8 FUNCTION matrix_cc_cc_noh(P,H1,H2,NHEL,IC)
368C
369C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
370C By the MadGraph Development Team
371C Please visit us at https://launchpad.net/madgraph5
372C
373C Returns amplitude squared summed/avg over colors
374C for the point with external lines W(0:6,NEXTERNAL)
375C
376C Process: c c > c c ta+ ta- / h QED=4
377C
378 IMPLICIT NONE
379C
380C CONSTANTS
381C
382 INTEGER NGRAPHS
383 parameter(ngraphs=64)
384 INTEGER NEXTERNAL
385 parameter(nexternal=6)
386 INTEGER NWAVEFUNCS, NCOLOR
387 parameter(nwavefuncs=16, ncolor=2)
388 real*8 zero
389 parameter(zero=0d0)
390 COMPLEX*16 IMAG1
391 parameter(imag1=(0d0,1d0))
392C
393C ARGUMENTS
394C
395 real*8 p(0:3,nexternal)
396 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
397C
398C LOCAL VARIABLES
399C
400 INTEGER I,J
401 COMPLEX*16 ZTEMP
402 real*8 denom(ncolor), cf(ncolor,ncolor)
403 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
404 COMPLEX*16 W(18,NWAVEFUNCS)
405 COMPLEX*16 DUM0,DUM1
406 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
407C
408C GLOBAL VARIABLES
409C
410 include 'coupl.inc'
411C
412C COLOR DATA
413C
414 DATA denom(1)/1/
415 DATA (cf(i, 1),i= 1, 2) / 9, 3/
416C 1 T(3,1) T(4,2)
417 DATA denom(2)/1/
418 DATA (cf(i, 2),i= 1, 2) / 3, 9/
419C 1 T(3,2) T(4,1)
420
421 INTEGER H1,H2
422 real*8 matrix
423 matrix_cc_cc_noh=0.d0
424 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
425 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
426
427
428C ----------
429C BEGIN CODE
430C ----------
431 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
432 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
433 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
434 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
435 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
436 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
437 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,7))
438 CALL ffv1p0_3(w(1,2),w(1,4),gc_2,zero,zero,w(1,8))
439 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
440C Amplitude(s) for diagram number 1
441 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
442 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
443C Amplitude(s) for diagram number 2
444 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
445 CALL ffv2_5_3(w(1,2),w(1,4),gc_51,gc_58,mz,wz,w(1,11))
446C Amplitude(s) for diagram number 3
447 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
448C Amplitude(s) for diagram number 4
449 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
450 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,10))
451 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
452C Amplitude(s) for diagram number 5
453 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
454 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
455C Amplitude(s) for diagram number 6
456 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
457C Amplitude(s) for diagram number 7
458 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
459C Amplitude(s) for diagram number 8
460 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
461 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
462 CALL ffv1_2(w(1,2),w(1,7),gc_2,zero,zero,w(1,9))
463C Amplitude(s) for diagram number 9
464 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_2,amp(9))
465 CALL ffv1_1(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
466C Amplitude(s) for diagram number 10
467 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_2,amp(10))
468 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
469C Amplitude(s) for diagram number 11
470 CALL ffv2_5_0(w(1,9),w(1,4),w(1,7),gc_51,gc_58,amp(11))
471C Amplitude(s) for diagram number 12
472 CALL ffv2_5_0(w(1,2),w(1,13),w(1,7),gc_51,gc_58,amp(12))
473 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
474 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
475C Amplitude(s) for diagram number 13
476 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_2,amp(13))
477 CALL ffv1_1(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
478C Amplitude(s) for diagram number 14
479 CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_2,amp(14))
480C Amplitude(s) for diagram number 15
481 CALL ffv2_5_0(w(1,9),w(1,4),w(1,7),gc_51,gc_58,amp(15))
482C Amplitude(s) for diagram number 16
483 CALL ffv2_5_0(w(1,2),w(1,14),w(1,7),gc_51,gc_58,amp(16))
484 CALL ffv2_5_2(w(1,2),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
485C Amplitude(s) for diagram number 17
486 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_2,amp(17))
487 CALL ffv2_5_1(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
488C Amplitude(s) for diagram number 18
489 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(18))
490C Amplitude(s) for diagram number 19
491 CALL ffv2_5_0(w(1,14),w(1,4),w(1,7),gc_51,gc_58,amp(19))
492C Amplitude(s) for diagram number 20
493 CALL ffv2_5_0(w(1,2),w(1,9),w(1,7),gc_51,gc_58,amp(20))
494 CALL ffv1p0_3(w(1,1),w(1,4),gc_2,zero,zero,w(1,9))
495 CALL ffv1p0_3(w(1,2),w(1,3),gc_2,zero,zero,w(1,14))
496 CALL ffv1_2(w(1,5),w(1,9),gc_3,mta,zero,w(1,10))
497C Amplitude(s) for diagram number 21
498 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(21))
499 CALL ffv1_1(w(1,6),w(1,9),gc_3,mta,zero,w(1,13))
500C Amplitude(s) for diagram number 22
501 CALL ffv1_0(w(1,5),w(1,13),w(1,14),gc_3,amp(22))
502 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,15))
503C Amplitude(s) for diagram number 23
504 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(23))
505C Amplitude(s) for diagram number 24
506 CALL ffv2_4_0(w(1,5),w(1,13),w(1,15),gc_50,gc_59,amp(24))
507 CALL ffv2_5_3(w(1,1),w(1,4),gc_51,gc_58,mz,wz,w(1,13))
508 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,10))
509C Amplitude(s) for diagram number 25
510 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(25))
511 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,16))
512C Amplitude(s) for diagram number 26
513 CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
514C Amplitude(s) for diagram number 27
515 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(27))
516C Amplitude(s) for diagram number 28
517 CALL ffv2_4_0(w(1,5),w(1,16),w(1,15),gc_50,gc_59,amp(28))
518 CALL ffv1_2(w(1,2),w(1,9),gc_2,zero,zero,w(1,16))
519C Amplitude(s) for diagram number 29
520 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_2,amp(29))
521 CALL ffv1_1(w(1,3),w(1,9),gc_2,zero,zero,w(1,5))
522C Amplitude(s) for diagram number 30
523 CALL ffv1_0(w(1,2),w(1,5),w(1,12),gc_2,amp(30))
524C Amplitude(s) for diagram number 31
525 CALL ffv2_5_0(w(1,16),w(1,3),w(1,7),gc_51,gc_58,amp(31))
526C Amplitude(s) for diagram number 32
527 CALL ffv2_5_0(w(1,2),w(1,5),w(1,7),gc_51,gc_58,amp(32))
528 CALL ffv1p0_3(w(1,1),w(1,4),gc_11,zero,zero,w(1,5))
529 CALL ffv1_2(w(1,2),w(1,5),gc_11,zero,zero,w(1,16))
530C Amplitude(s) for diagram number 33
531 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_2,amp(33))
532 CALL ffv1_1(w(1,3),w(1,5),gc_11,zero,zero,w(1,9))
533C Amplitude(s) for diagram number 34
534 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(34))
535C Amplitude(s) for diagram number 35
536 CALL ffv2_5_0(w(1,16),w(1,3),w(1,7),gc_51,gc_58,amp(35))
537C Amplitude(s) for diagram number 36
538 CALL ffv2_5_0(w(1,2),w(1,9),w(1,7),gc_51,gc_58,amp(36))
539 CALL ffv2_5_2(w(1,2),w(1,13),gc_51,gc_58,zero,zero,w(1,9))
540C Amplitude(s) for diagram number 37
541 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(37))
542 CALL ffv2_5_1(w(1,3),w(1,13),gc_51,gc_58,zero,zero,w(1,16))
543C Amplitude(s) for diagram number 38
544 CALL ffv1_0(w(1,2),w(1,16),w(1,12),gc_2,amp(38))
545C Amplitude(s) for diagram number 39
546 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(39))
547C Amplitude(s) for diagram number 40
548 CALL ffv2_5_0(w(1,2),w(1,16),w(1,7),gc_51,gc_58,amp(40))
549 CALL ffv1_2(w(1,1),w(1,14),gc_2,zero,zero,w(1,16))
550C Amplitude(s) for diagram number 41
551 CALL ffv1_0(w(1,16),w(1,4),w(1,12),gc_2,amp(41))
552 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,9))
553C Amplitude(s) for diagram number 42
554 CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_2,amp(42))
555C Amplitude(s) for diagram number 43
556 CALL ffv2_5_0(w(1,16),w(1,4),w(1,7),gc_51,gc_58,amp(43))
557 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,16))
558C Amplitude(s) for diagram number 44
559 CALL ffv1_0(w(1,16),w(1,4),w(1,14),gc_2,amp(44))
560 CALL ffv1p0_3(w(1,2),w(1,3),gc_11,zero,zero,w(1,14))
561 CALL ffv1_2(w(1,1),w(1,14),gc_11,zero,zero,w(1,13))
562C Amplitude(s) for diagram number 45
563 CALL ffv1_0(w(1,13),w(1,4),w(1,12),gc_2,amp(45))
564C Amplitude(s) for diagram number 46
565 CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_11,amp(46))
566C Amplitude(s) for diagram number 47
567 CALL ffv2_5_0(w(1,13),w(1,4),w(1,7),gc_51,gc_58,amp(47))
568C Amplitude(s) for diagram number 48
569 CALL ffv1_0(w(1,16),w(1,4),w(1,14),gc_11,amp(48))
570 CALL ffv2_5_2(w(1,1),w(1,15),gc_51,gc_58,zero,zero,w(1,14))
571C Amplitude(s) for diagram number 49
572 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_2,amp(49))
573C Amplitude(s) for diagram number 50
574 CALL ffv2_5_0(w(1,9),w(1,4),w(1,15),gc_51,gc_58,amp(50))
575C Amplitude(s) for diagram number 51
576 CALL ffv2_5_0(w(1,14),w(1,4),w(1,7),gc_51,gc_58,amp(51))
577C Amplitude(s) for diagram number 52
578 CALL ffv2_5_0(w(1,16),w(1,4),w(1,15),gc_51,gc_58,amp(52))
579 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,15))
580C Amplitude(s) for diagram number 53
581 CALL ffv1_0(w(1,15),w(1,3),w(1,12),gc_2,amp(53))
582C Amplitude(s) for diagram number 54
583 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_2,amp(54))
584C Amplitude(s) for diagram number 55
585 CALL ffv2_5_0(w(1,15),w(1,3),w(1,7),gc_51,gc_58,amp(55))
586C Amplitude(s) for diagram number 56
587 CALL ffv1_0(w(1,16),w(1,3),w(1,8),gc_2,amp(56))
588 CALL ffv1p0_3(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
589 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
590C Amplitude(s) for diagram number 57
591 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(57))
592C Amplitude(s) for diagram number 58
593 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(58))
594C Amplitude(s) for diagram number 59
595 CALL ffv2_5_0(w(1,4),w(1,3),w(1,7),gc_51,gc_58,amp(59))
596C Amplitude(s) for diagram number 60
597 CALL ffv1_0(w(1,16),w(1,3),w(1,8),gc_11,amp(60))
598 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
599C Amplitude(s) for diagram number 61
600 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(61))
601C Amplitude(s) for diagram number 62
602 CALL ffv2_5_0(w(1,9),w(1,3),w(1,11),gc_51,gc_58,amp(62))
603C Amplitude(s) for diagram number 63
604 CALL ffv2_5_0(w(1,8),w(1,3),w(1,7),gc_51,gc_58,amp(63))
605C Amplitude(s) for diagram number 64
606 CALL ffv2_5_0(w(1,16),w(1,3),w(1,11),gc_51,gc_58,amp(64))
607 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
608 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
609 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
610 $ +1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)+1./2.*amp(36)
611 $ +1./2.*amp(45)+1./2.*amp(46)+1./2.*amp(47)+1./2.*amp(48)
612 $ -amp(53)-amp(54)-amp(55)-amp(56)+1./6.*amp(57)+1./6.*amp(58)
613 $ +1./6.*amp(59)+1./6.*amp(60)-amp(61)-amp(62)-amp(63)-amp(64)
614 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
615 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
616 $ +amp(28)+amp(29)+amp(30)+amp(31)+amp(32)-1./6.*amp(33)
617 $ -1./6.*amp(34)-1./6.*amp(35)-1./6.*amp(36)+amp(37)+amp(38)
618 $ +amp(39)+amp(40)+amp(41)+amp(42)+amp(43)+amp(44)-1./6.*amp(45)
619 $ -1./6.*amp(46)-1./6.*amp(47)-1./6.*amp(48)+amp(49)+amp(50)
620 $ +amp(51)+amp(52)-1./2.*amp(57)-1./2.*amp(58)-1./2.*amp(59)
621 $ -1./2.*amp(60)
622
623 matrix = 0.d0
624 DO i = 1, ncolor
625 ztemp = (0.d0,0.d0)
626 DO j = 1, ncolor
627 ztemp = ztemp + cf(j,i)*jamp(j)
628 ENDDO
629 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
630 ENDDO
631 matrix_cc_cc_noh=matrix
632 ENDIF
633 ENDIF
634 END
635
636C -------------end subprocess initiated by UU---------
637C DONE 26.10.2015
638