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