C++ Interface to Tauola
UU.f
1C JK. corrected 20.11.2014
2c------------------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
9C 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
25C ---------begin subprocess UU->jjH with H-> tautau
26c ---------------------jj=uu only -------------
27
28 SUBROUTINE uu_uu_h(P,H1,H2,ANS)
29C
30C
31C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
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: u u > u u 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_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/
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_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)
149C
150C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
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: u u > u u 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_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
209C ----------
210C BEGIN CODE
211C ----------
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))
221C 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))
225C 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
242C----------------end subprocess uu->jjH---------
243
244c----------------begin subprocess uu->jj_noH
245
246
247
248 SUBROUTINE uu_uu_noh(P,H1,H2,ANS)
249C
250C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
251C By the MadGraph Development Team
252C Please visit us at https://launchpad.net/madgraph5
253C
254C MadGraph StandAlone Version
255C
256C Returns amplitude squared summed/avg over colors
257C and helicities
258C for the point in phase space P(0:3,NEXTERNAL)
259C
260C Process: u u > u u ta+ ta- / h QED=4
261C
262 IMPLICIT NONE
263C
264C CONSTANTS
265C
266 INTEGER NEXTERNAL
267 parameter(nexternal=6)
268 INTEGER NCOMB
269 parameter( ncomb=64)
270C
271C ARGUMENTS
272C
273 real*8 p(0:3,nexternal),ans
274 INTEGER H1,H2
275C
276C LOCAL VARIABLES
277C
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/
351C ----------
352C BEGIN CODE
353C ----------
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)
367C
368C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
369C By the MadGraph Development Team
370C Please visit us at https://launchpad.net/madgraph5
371C
372C Returns amplitude squared summed/avg over colors
373C for the point with external lines W(0:6,NEXTERNAL)
374C
375C Process: u u > u u ta+ ta- / h QED=4
376C
377 IMPLICIT NONE
378C
379C CONSTANTS
380C
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))
391C
392C ARGUMENTS
393C
394 real*8 p(0:3,nexternal)
395 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
396C
397C LOCAL VARIABLES
398C
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)/
406C
407C GLOBAL VARIABLES
408C
409 include 'coupl.inc'
410C
411C COLOR DATA
412C
413 DATA denom(1)/1/
414 DATA (cf(i, 1),i= 1, 2) / 9, 3/
415C 1 T(3,1) T(4,2)
416 DATA denom(2)/1/
417 DATA (cf(i, 2),i= 1, 2) / 3, 9/
418C 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
426C ----------
427C BEGIN CODE
428C ----------
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))
438C 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))
441C 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))
444C 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))
446C 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))
450C 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))
453C Amplitude(s) for diagram number 6
454 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
455C 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))
457C 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))
461C 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))
464C 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))
467C 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))
469C 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))
473C 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))
476C Amplitude(s) for diagram number 14
477 CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_2,amp(14))
478C 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))
480C 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))
483C 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))
486C Amplitude(s) for diagram number 18
487 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(18))
488C 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))
490C 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))
495C 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))
498C 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))
501C 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))
503C 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))
507C 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))
510C Amplitude(s) for diagram number 26
511 CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
512C 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))
514C 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))
517C 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))
520C Amplitude(s) for diagram number 30
521 CALL ffv1_0(w(1,2),w(1,5),w(1,12),gc_2,amp(30))
522C 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))
524C 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))
528C 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))
531C Amplitude(s) for diagram number 34
532 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_2,amp(34))
533C 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))
535C 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))
538C 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))
541C Amplitude(s) for diagram number 38
542 CALL ffv1_0(w(1,2),w(1,16),w(1,12),gc_2,amp(38))
543C 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))
545C 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))
548C 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))
551C Amplitude(s) for diagram number 42
552 CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_2,amp(42))
553C 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))
556C 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))
560C Amplitude(s) for diagram number 45
561 CALL ffv1_0(w(1,13),w(1,4),w(1,12),gc_2,amp(45))
562C Amplitude(s) for diagram number 46
563 CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_11,amp(46))
564C 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))
566C 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))
569C Amplitude(s) for diagram number 49
570 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_2,amp(49))
571C 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))
573C 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))
575C 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))
578C Amplitude(s) for diagram number 53
579 CALL ffv1_0(w(1,15),w(1,3),w(1,12),gc_2,amp(53))
580C Amplitude(s) for diagram number 54
581 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_2,amp(54))
582C 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))
584C 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))
588C Amplitude(s) for diagram number 57
589 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_2,amp(57))
590C Amplitude(s) for diagram number 58
591 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(58))
592C 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))
594C 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))
597C Amplitude(s) for diagram number 61
598 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_2,amp(61))
599C 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))
601C 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))
603C 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
634C -------------end subprocess initiated by UU---------