C++ Interface to Tauola
CU.f
1 C JK. poprawka 10.11.2015
2 c------------------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 
10 C 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 
26 C ---------begin subprocess CU->jjH with H-> tautau
27 c ---------------------jj=CU only -------------
28 
29  SUBROUTINE cu_uc_h(P,H1,H2,ANS)
30 C
31 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
32 C By the MadGraph Development Team
33 C Please visit us at https://launchpad.net/madgraph5
34 C
35 C MadGraph StandAlone Version
36 C
37 C Returns amplitude squared summed/avg over colors
38 C and helicities
39 C for the point in phase space P(0:3,NEXTERNAL)
40 C
41 C Process: c u > u c h WEIGHTED=6
42 C * Decay: h > ta+ ta- WEIGHTED=2
43 C
44  IMPLICIT NONE
45 C
46 C CONSTANTS
47 C
48  INTEGER NEXTERNAL
49  parameter(nexternal=6)
50  INTEGER NCOMB
51  parameter( ncomb=64)
52 C
53 C ARGUMENTS
54 C
55  real*8 p(0:3,nexternal),ans
56  INTEGER H1,H2
57 C
58 C LOCAL VARIABLES
59 C
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/
133 C ----------
134 C BEGIN CODE
135 C ----------
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)
149 C
150 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
151 C By the MadGraph Development Team
152 C Please visit us at https://launchpad.net/madgraph5
153 C
154 C Returns amplitude squared summed/avg over colors
155 C for the point with external lines W(0:6,NEXTERNAL)
156 C
157 C Process: c u > u c h WEIGHTED=6
158 C * Decay: h > ta+ ta- WEIGHTED=2
159 C
160  IMPLICIT NONE
161 C
162 C CONSTANTS
163 C
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))
174 C
175 C ARGUMENTS
176 C
177  real*8 p(0:3,nexternal)
178  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
179 C
180 C LOCAL VARIABLES
181 C
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)/
189 C
190 C GLOBAL VARIABLES
191 C
192  include 'coupl.inc'
193 C
194 C COLOR DATA
195 C
196  DATA denom(1)/1/
197  DATA (cf(i, 1),i= 1, 1) / 9/
198 C 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 
206 C ----------
207 C BEGIN CODE
208 C ----------
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))
218 C 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
234 C----------------end subprocess CU->jjH---------
235 
236 c----------------begin subprocess CU->jj_noH
237 
238 
239 
240  SUBROUTINE cu_uc_noh(P,H1,H2,ANS)
241 C
242 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
243 C By the MadGraph Development Team
244 C Please visit us at https://launchpad.net/madgraph5
245 C
246 C MadGraph StandAlone Version
247 C
248 C Returns amplitude squared summed/avg over colors
249 C and helicities
250 C for the point in phase space P(0:3,NEXTERNAL)
251 C
252 C Process: c u > u c ta+ ta- / h QED=4
253 C
254  IMPLICIT NONE
255 C
256 C CONSTANTS
257 C
258  INTEGER NEXTERNAL
259  parameter(nexternal=6)
260  INTEGER NCOMB
261  parameter( ncomb=64)
262 C
263 C ARGUMENTS
264 C
265  real*8 p(0:3,nexternal),ans
266  INTEGER H1,H2
267 C
268 C LOCAL VARIABLES
269 C
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/
343 C ----------
344 C BEGIN CODE
345 C ----------
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)
360 C
361 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
362 C By the MadGraph Development Team
363 C Please visit us at https://launchpad.net/madgraph5
364 C
365 C Returns amplitude squared summed/avg over colors
366 C for the point with external lines W(0:6,NEXTERNAL)
367 C
368 C Process: c u > u c ta+ ta- / h QED=4
369 C
370  IMPLICIT NONE
371 C
372 C CONSTANTS
373 C
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))
384 C
385 C ARGUMENTS
386 C
387  real*8 p(0:3,nexternal)
388  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
389 C
390 C LOCAL VARIABLES
391 C
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)/
399 C
400 C GLOBAL VARIABLES
401 C
402  include 'coupl.inc'
403 C
404 C COLOR DATA
405 C
406  DATA denom(1)/1/
407  DATA (cf(i, 1),i= 1, 2) / 9, 3/
408 C 1 T(3,1) T(4,2)
409  DATA denom(2)/1/
410  DATA (cf(i, 2),i= 1, 2) / 3, 9/
411 C 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 
419 C ----------
420 C BEGIN CODE
421 C ----------
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))
431 C 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))
434 C 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))
437 C 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))
439 C 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))
443 C 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))
446 C Amplitude(s) for diagram number 6
447  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
448 C 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))
450 C 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))
454 C 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))
457 C 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))
460 C 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))
462 C 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))
466 C 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))
469 C Amplitude(s) for diagram number 14
470  CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_2,amp(14))
471 C 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))
473 C 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))
476 C 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))
479 C Amplitude(s) for diagram number 18
480  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(18))
481 C 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))
483 C 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))
486 C 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))
489 C Amplitude(s) for diagram number 22
490  CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_2,amp(22))
491 C 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))
494 C 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))
498 C Amplitude(s) for diagram number 25
499  CALL ffv1_0(w(1,3),w(1,4),w(1,12),gc_2,amp(25))
500 C Amplitude(s) for diagram number 26
501  CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_11,amp(26))
502 C 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))
504 C 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))
507 C Amplitude(s) for diagram number 29
508  CALL ffv1_0(w(1,8),w(1,4),w(1,12),gc_2,amp(29))
509 C 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))
511 C 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))
513 C 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 
536 C -------------end subprocess initiated by CU---------
537 C DONE 26.10.2015, poprawka 10.11.2015
538 
Generated by doxygen 1.9.1