C++ Interface to Tauola
CC.f
1 C JK. corrected 26.10.2015
2 c------------------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 
10 C 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 
26 C ---------begin subprocess CC->jjH with H-> tautau
27 c ---------------------jj=uu only -------------
28 
29  SUBROUTINE cc_cc_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 c > c 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_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/
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_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)
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 c > c 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=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))
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, 2) / 9, 3/
198 C 1 T(3,1) T(4,2)
199  DATA denom(2)/1/
200  DATA (cf(i, 2),i= 1, 2) / 3, 9/
201 C 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 
210 C ----------
211 C BEGIN CODE
212 C ----------
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))
222 C 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))
226 C 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
243 C----------------end subprocess CC->jjH---------
244 
245 c----------------begin subprocess CC->jj_noH
246 
247 
248 
249  SUBROUTINE cc_cc_noh(P,H1,H2,ANS)
250 C
251 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
252 C By the MadGraph Development Team
253 C Please visit us at https://launchpad.net/madgraph5
254 C
255 C MadGraph StandAlone Version
256 C
257 C Returns amplitude squared summed/avg over colors
258 C and helicities
259 C for the point in phase space P(0:3,NEXTERNAL)
260 C
261 C Process: c c > c c ta+ ta- / h QED=4
262 C
263  IMPLICIT NONE
264 C
265 C CONSTANTS
266 C
267  INTEGER NEXTERNAL
268  parameter(nexternal=6)
269  INTEGER NCOMB
270  parameter( ncomb=64)
271 C
272 C ARGUMENTS
273 C
274  real*8 p(0:3,nexternal),ans
275  INTEGER H1,H2
276 C
277 C LOCAL VARIABLES
278 C
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/
352 C ----------
353 C BEGIN CODE
354 C ----------
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)
368 C
369 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
370 C By the MadGraph Development Team
371 C Please visit us at https://launchpad.net/madgraph5
372 C
373 C Returns amplitude squared summed/avg over colors
374 C for the point with external lines W(0:6,NEXTERNAL)
375 C
376 C Process: c c > c c ta+ ta- / h QED=4
377 C
378  IMPLICIT NONE
379 C
380 C CONSTANTS
381 C
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))
392 C
393 C ARGUMENTS
394 C
395  real*8 p(0:3,nexternal)
396  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
397 C
398 C LOCAL VARIABLES
399 C
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)/
407 C
408 C GLOBAL VARIABLES
409 C
410  include 'coupl.inc'
411 C
412 C COLOR DATA
413 C
414  DATA denom(1)/1/
415  DATA (cf(i, 1),i= 1, 2) / 9, 3/
416 C 1 T(3,1) T(4,2)
417  DATA denom(2)/1/
418  DATA (cf(i, 2),i= 1, 2) / 3, 9/
419 C 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 
428 C ----------
429 C BEGIN CODE
430 C ----------
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))
440 C 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))
443 C 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))
446 C 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))
448 C 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))
452 C 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))
455 C Amplitude(s) for diagram number 6
456  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
457 C 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))
459 C 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))
463 C 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))
466 C 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))
469 C 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))
471 C 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))
475 C 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))
478 C Amplitude(s) for diagram number 14
479  CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_2,amp(14))
480 C 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))
482 C 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))
485 C 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))
488 C Amplitude(s) for diagram number 18
489  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(18))
490 C 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))
492 C 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))
497 C 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))
500 C 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))
503 C 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))
505 C 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))
509 C 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))
512 C Amplitude(s) for diagram number 26
513  CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
514 C 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))
516 C 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))
519 C 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))
522 C Amplitude(s) for diagram number 30
523  CALL ffv1_0(w(1,2),w(1,5),w(1,12),gc_2,amp(30))
524 C 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))
526 C 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))
530 C 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))
533 C Amplitude(s) for diagram number 34
534  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(34))
535 C 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))
537 C 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))
540 C 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))
543 C Amplitude(s) for diagram number 38
544  CALL ffv1_0(w(1,2),w(1,16),w(1,12),gc_2,amp(38))
545 C 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))
547 C 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))
550 C 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))
553 C Amplitude(s) for diagram number 42
554  CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_2,amp(42))
555 C 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))
558 C 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))
562 C Amplitude(s) for diagram number 45
563  CALL ffv1_0(w(1,13),w(1,4),w(1,12),gc_2,amp(45))
564 C Amplitude(s) for diagram number 46
565  CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_11,amp(46))
566 C 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))
568 C 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))
571 C Amplitude(s) for diagram number 49
572  CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_2,amp(49))
573 C 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))
575 C 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))
577 C 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))
580 C Amplitude(s) for diagram number 53
581  CALL ffv1_0(w(1,15),w(1,3),w(1,12),gc_2,amp(53))
582 C Amplitude(s) for diagram number 54
583  CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_2,amp(54))
584 C 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))
586 C 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))
590 C Amplitude(s) for diagram number 57
591  CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(57))
592 C Amplitude(s) for diagram number 58
593  CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(58))
594 C 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))
596 C 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))
599 C Amplitude(s) for diagram number 61
600  CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(61))
601 C 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))
603 C 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))
605 C 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 
636 C -------------end subprocess initiated by UU---------
637 C DONE 26.10.2015
638 
Generated by doxygen 1.9.1