C++ Interface to Tauola
CU.f
1C JK. poprawka 10.11.2015
2c------------------begin subprocess initiated by CU----------
3
4 SUBROUTINE cu(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,2 I4=0,4 ! only CU->UC
11
12 ans=0.d0
13
14 IF(key.EQ.1) THEN
15 IF(i3.EQ.2 .AND. i4.EQ.4) CALL cu_uc_h(p,h1,h2,ans)
16 IF(i3.EQ.0 .AND. i4.EQ.0) CALL cu_uc_h(p,h1,h2,ans)
17 ELSE IF(key.EQ.0) THEN
18 IF(i3.EQ.2 .AND. i4.EQ.4) CALL cu_uc_noh(p,h1,h2,ans)
19 IF(i3.EQ.0 .AND. i4.EQ.0) CALL cu_uc_noh(p,h1,h2,ans)
20 ELSE
21 WRITE(*,*) 'spin=2 NOT FINISHED'
22 stop
23 ENDIF
24 END ! SUBROUTINE CU
25
26C ---------begin subprocess CU->jjH with H-> tautau
27c ---------------------jj=CU only -------------
28
29 SUBROUTINE cu_uc_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 u > u 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_cu_uc_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/36/
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_cu_uc_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_cu_uc_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 u > u 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=1)
166 INTEGER NEXTERNAL
167 parameter(nexternal=6)
168 INTEGER NWAVEFUNCS, NCOLOR
169 parameter(nwavefuncs=7, ncolor=1)
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, 1) / 9/
198C 1 T(3,2) T(4,1)
199
200 INTEGER H1,H2
201 real*8 matrix
202 matrix_cu_uc_h=0.d0
203 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
204 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
205
206C ----------
207C BEGIN CODE
208C ----------
209 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
210 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
211 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
212 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
213 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
214 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
215 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
216 CALL ffv2_5_3(w(1,1),w(1,4),gc_51,gc_58,mz,wz,w(1,6))
217 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,4))
218C Amplitude(s) for diagram number 1
219 CALL vvs1_0(w(1,6),w(1,4),w(1,7),gc_81,amp(1))
220 jamp(1)=+amp(1)
221
222 matrix = 0.d0
223 DO i = 1, ncolor
224 ztemp = (0.d0,0.d0)
225 DO j = 1, ncolor
226 ztemp = ztemp + cf(j,i)*jamp(j)
227 ENDDO
228 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
229 ENDDO
230 matrix_cu_uc_h=matrix
231 ENDIF
232 ENDIF
233 END
234C----------------end subprocess CU->jjH---------
235
236c----------------begin subprocess CU->jj_noH
237
238
239
240 SUBROUTINE cu_uc_noh(P,H1,H2,ANS)
241C
242C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
243C By the MadGraph Development Team
244C Please visit us at https://launchpad.net/madgraph5
245C
246C MadGraph StandAlone Version
247C
248C Returns amplitude squared summed/avg over colors
249C and helicities
250C for the point in phase space P(0:3,NEXTERNAL)
251C
252C Process: c u > u c ta+ ta- / h QED=4
253C
254 IMPLICIT NONE
255C
256C CONSTANTS
257C
258 INTEGER NEXTERNAL
259 parameter(nexternal=6)
260 INTEGER NCOMB
261 parameter( ncomb=64)
262C
263C ARGUMENTS
264C
265 real*8 p(0:3,nexternal),ans
266 INTEGER H1,H2
267C
268C LOCAL VARIABLES
269C
270 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
271 real*8 t
272 real*8 matrix_cu_uc_noh
273 INTEGER IHEL,IDEN, I
274 INTEGER JC(NEXTERNAL)
275 LOGICAL GOODHEL(NCOMB)
276 DATA ntry/0/
277 DATA goodhel/ncomb*.false./
278 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
279 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
280 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
281 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
282 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
283 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
284 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
285 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
286 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
287 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
288 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
289 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
290 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
291 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
292 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
293 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
294 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
295 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
296 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
297 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
298 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
299 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
300 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
301 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
302 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
303 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
304 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
305 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
306 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
307 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
308 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
309 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
310 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
311 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
312 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
313 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
314 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
315 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
316 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
317 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
318 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
319 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
320 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
321 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
322 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
323 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
324 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
325 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
326 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
327 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
328 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
329 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
330 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
331 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
332 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
333 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
334 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
335 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
336 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
337 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
338 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
339 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
340 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
341 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
342 DATA iden/36/
343C ----------
344C BEGIN CODE
345C ----------
346
347 DO ihel=1,nexternal
348 jc(ihel) = +1
349 ENDDO
350 ans = 0d0
351 DO ihel=1,ncomb
352 t=matrix_cu_uc_noh(p,h1,h2,nhel(1,ihel),jc(1))
353 ans=ans+t
354 ENDDO
355 ans=ans/dble(iden)
356 END
357
358
359 REAL*8 FUNCTION matrix_cu_uc_noh(P,H1,H2,NHEL,IC)
360C
361C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
362C By the MadGraph Development Team
363C Please visit us at https://launchpad.net/madgraph5
364C
365C Returns amplitude squared summed/avg over colors
366C for the point with external lines W(0:6,NEXTERNAL)
367C
368C Process: c u > u c ta+ ta- / h QED=4
369C
370 IMPLICIT NONE
371C
372C CONSTANTS
373C
374 INTEGER NGRAPHS
375 parameter(ngraphs=32)
376 INTEGER NEXTERNAL
377 parameter(nexternal=6)
378 INTEGER NWAVEFUNCS, NCOLOR
379 parameter(nwavefuncs=13, ncolor=2)
380 real*8 zero
381 parameter(zero=0d0)
382 COMPLEX*16 IMAG1
383 parameter(imag1=(0d0,1d0))
384C
385C ARGUMENTS
386C
387 real*8 p(0:3,nexternal)
388 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
389C
390C LOCAL VARIABLES
391C
392 INTEGER I,J
393 COMPLEX*16 ZTEMP
394 real*8 denom(ncolor), cf(ncolor,ncolor)
395 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
396 COMPLEX*16 W(18,NWAVEFUNCS)
397 COMPLEX*16 DUM0,DUM1
398 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
399C
400C GLOBAL VARIABLES
401C
402 include 'coupl.inc'
403C
404C COLOR DATA
405C
406 DATA denom(1)/1/
407 DATA (cf(i, 1),i= 1, 2) / 9, 3/
408C 1 T(3,1) T(4,2)
409 DATA denom(2)/1/
410 DATA (cf(i, 2),i= 1, 2) / 3, 9/
411C 1 T(3,2) T(4,1)
412
413 INTEGER H1,H2
414 real*8 matrix
415 matrix_cu_uc_noh=0.d0
416 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
417 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
418
419C ----------
420C BEGIN CODE
421C ----------
422 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
423 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
424 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
425 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
426 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
427 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
428 CALL ffv1p0_3(w(1,1),w(1,4),gc_2,zero,zero,w(1,7))
429 CALL ffv1p0_3(w(1,2),w(1,3),gc_2,zero,zero,w(1,8))
430 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
431C Amplitude(s) for diagram number 1
432 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
433 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
434C Amplitude(s) for diagram number 2
435 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
436 CALL ffv2_5_3(w(1,2),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
437C Amplitude(s) for diagram number 3
438 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
439C Amplitude(s) for diagram number 4
440 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
441 CALL ffv2_5_3(w(1,1),w(1,4),gc_51,gc_58,mz,wz,w(1,10))
442 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
443C Amplitude(s) for diagram number 5
444 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
445 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
446C Amplitude(s) for diagram number 6
447 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
448C Amplitude(s) for diagram number 7
449 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
450C Amplitude(s) for diagram number 8
451 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
452 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
453 CALL ffv1_2(w(1,2),w(1,7),gc_2,zero,zero,w(1,9))
454C Amplitude(s) for diagram number 9
455 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(9))
456 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,13))
457C Amplitude(s) for diagram number 10
458 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_2,amp(10))
459 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
460C Amplitude(s) for diagram number 11
461 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(11))
462C Amplitude(s) for diagram number 12
463 CALL ffv2_5_0(w(1,2),w(1,13),w(1,7),gc_51,gc_58,amp(12))
464 CALL ffv1p0_3(w(1,1),w(1,4),gc_11,zero,zero,w(1,13))
465 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
466C Amplitude(s) for diagram number 13
467 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(13))
468 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,6))
469C Amplitude(s) for diagram number 14
470 CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_2,amp(14))
471C Amplitude(s) for diagram number 15
472 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(15))
473C Amplitude(s) for diagram number 16
474 CALL ffv2_5_0(w(1,2),w(1,6),w(1,7),gc_51,gc_58,amp(16))
475 CALL ffv2_5_2(w(1,2),w(1,10),gc_51,gc_58,zero,zero,w(1,6))
476C Amplitude(s) for diagram number 17
477 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_2,amp(17))
478 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
479C Amplitude(s) for diagram number 18
480 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(18))
481C Amplitude(s) for diagram number 19
482 CALL ffv2_5_0(w(1,6),w(1,3),w(1,7),gc_51,gc_58,amp(19))
483C Amplitude(s) for diagram number 20
484 CALL ffv2_5_0(w(1,2),w(1,9),w(1,7),gc_51,gc_58,amp(20))
485 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,9))
486C Amplitude(s) for diagram number 21
487 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_2,amp(21))
488 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,6))
489C Amplitude(s) for diagram number 22
490 CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_2,amp(22))
491C Amplitude(s) for diagram number 23
492 CALL ffv2_5_0(w(1,9),w(1,4),w(1,7),gc_51,gc_58,amp(23))
493 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,9))
494C Amplitude(s) for diagram number 24
495 CALL ffv1_0(w(1,9),w(1,4),w(1,8),gc_2,amp(24))
496 CALL ffv1p0_3(w(1,2),w(1,3),gc_11,zero,zero,w(1,8))
497 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,3))
498C Amplitude(s) for diagram number 25
499 CALL ffv1_0(w(1,3),w(1,4),w(1,12),gc_2,amp(25))
500C Amplitude(s) for diagram number 26
501 CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_11,amp(26))
502C Amplitude(s) for diagram number 27
503 CALL ffv2_5_0(w(1,3),w(1,4),w(1,7),gc_51,gc_58,amp(27))
504C Amplitude(s) for diagram number 28
505 CALL ffv1_0(w(1,9),w(1,4),w(1,8),gc_11,amp(28))
506 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
507C Amplitude(s) for diagram number 29
508 CALL ffv1_0(w(1,8),w(1,4),w(1,12),gc_2,amp(29))
509C Amplitude(s) for diagram number 30
510 CALL ffv2_5_0(w(1,6),w(1,4),w(1,11),gc_51,gc_58,amp(30))
511C Amplitude(s) for diagram number 31
512 CALL ffv2_5_0(w(1,8),w(1,4),w(1,7),gc_51,gc_58,amp(31))
513C Amplitude(s) for diagram number 32
514 CALL ffv2_5_0(w(1,9),w(1,4),w(1,11),gc_51,gc_58,amp(32))
515 jamp(1)=+1./2.*(+amp(13)+amp(14)+amp(15)+amp(16)+amp(25)+amp(26)
516 $ +amp(27)+amp(28))
517 jamp(2)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
518 $ +amp(9)+amp(10)+amp(11)+amp(12)-1./6.*amp(13)-1./6.*amp(14)
519 $ -1./6.*amp(15)-1./6.*amp(16)+amp(17)+amp(18)+amp(19)+amp(20)
520 $ +amp(21)+amp(22)+amp(23)+amp(24)-1./6.*amp(25)-1./6.*amp(26)
521 $ -1./6.*amp(27)-1./6.*amp(28)+amp(29)+amp(30)+amp(31)+amp(32)
522
523 matrix = 0.d0
524 DO i = 1, ncolor
525 ztemp = (0.d0,0.d0)
526 DO j = 1, ncolor
527 ztemp = ztemp + cf(j,i)*jamp(j)
528 ENDDO
529 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
530 ENDDO
531 matrix_cu_uc_noh=matrix
532 ENDIF
533 ENDIF
534 END
535
536C -------------end subprocess initiated by CU---------
537C DONE 26.10.2015, poprawka 10.11.2015
538