C++ Interface to Tauola
DD.f
1C JK. done 30.07.2014
2c------------------begin subprocess initiated by DD----------
3 SUBROUTINE dd(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,1 I4=0,1 ! only dd->dd
10
11 ans=0.d0
12
13 IF(key.EQ.1) THEN
14 IF(i3.EQ.1 .AND. i4.EQ.1) CALL dd_dd_h(p,h1,h2,ans)
15 IF(i3.EQ.0 .AND. i4.EQ.0) CALL dd_dd_h(p,h1,h2,ans)
16 ELSE IF(key.EQ.0) THEN
17 IF(i3.EQ.1 .AND. i4.EQ.1) CALL dd_dd_noh(p,h1,h2,ans)
18 IF(i3.EQ.0 .AND. i4.EQ.0) CALL dd_dd_noh(p,h1,h2,ans)
19 ELSE
20 WRITE(*,*) 'spin=2 NOT FINISHED'
21 stop
22 ENDIF
23 END SUBROUTINE dd
24
25C ---------begin subprocess DD->jjH with H-> tautau
26c ---------------------jj=dd only -------------
27
28 SUBROUTINE dd_dd_h(P,H1,H2,ANS)
29C
30C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
31C By the MadGraph Development Team
32C Please visit us at https://launchpad.net/madgraph5
33C
34C MadGraph StandAlone Version
35C
36C Returns amplitude squared summed/avg over colors
37C and helicities
38C for the point in phase space P(0:3,NEXTERNAL)
39C
40C Process: d d > d d h WEIGHTED=6
41C * Decay: h > ta+ ta- WEIGHTED=2
42C
43 IMPLICIT NONE
44C
45C CONSTANTS
46C
47 INTEGER NEXTERNAL
48 parameter(nexternal=6)
49 INTEGER NCOMB
50 parameter( ncomb=64)
51C
52C ARGUMENTS
53C
54 real*8 p(0:3,nexternal),ans
55 INTEGER H1,H2
56C
57C LOCAL VARIABLES
58C
59 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
60 real*8 t
61 real*8 matrix_dd_dd_h
62 INTEGER IHEL,IDEN, I
63 INTEGER JC(NEXTERNAL)
64 LOGICAL GOODHEL(NCOMB)
65 DATA ntry/0/
66 DATA goodhel/ncomb*.false./
67 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
68 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
69 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
70 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
71 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
72 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
73 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
74 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
75 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
76 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
77 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
78 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
79 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
80 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
81 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
82 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
83 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
84 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
85 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
86 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
87 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
88 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
89 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
90 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
91 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
92 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
93 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
94 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
95 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
96 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
97 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
98 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
99 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
100 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
101 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
102 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
103 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
104 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
105 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
106 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
107 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
108 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
109 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
110 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
111 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
112 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
113 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
114 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
115 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
116 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
117 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
118 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
119 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
120 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
121 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
122 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
123 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
124 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
125 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
126 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
127 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
128 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
129 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
130 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
131 DATA iden/72/
132C ----------
133C BEGIN CODE
134C ----------
135 DO ihel=1,nexternal
136 jc(ihel) = +1
137 ENDDO
138 ans = 0d0
139 DO ihel=1,ncomb
140 t=matrix_dd_dd_h(p ,h1,h2,nhel(1,ihel),jc(1))
141 ans=ans+t
142 ENDDO
143 ans=ans/dble(iden)
144 END
145
146
147 REAL*8 FUNCTION matrix_dd_dd_h(P,H1,H2,NHEL,IC)
148C
149C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
150C By the MadGraph Development Team
151C Please visit us at https://launchpad.net/madgraph5
152C
153C Returns amplitude squared summed/avg over colors
154C for the point with external lines W(0:6,NEXTERNAL)
155C
156C Process: d d > d d h WEIGHTED=6
157C * Decay: h > ta+ ta- WEIGHTED=2
158C
159 IMPLICIT NONE
160C
161C CONSTANTS
162C
163 INTEGER NGRAPHS
164 parameter(ngraphs=2)
165 INTEGER NEXTERNAL
166 parameter(nexternal=6)
167 INTEGER NWAVEFUNCS, NCOLOR
168 parameter(nwavefuncs=7, ncolor=2)
169 real*8 zero
170 parameter(zero=0d0)
171 COMPLEX*16 IMAG1
172 parameter(imag1=(0d0,1d0))
173C
174C ARGUMENTS
175C
176 real*8 p(0:3,nexternal)
177 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
178C
179C LOCAL VARIABLES
180C
181 INTEGER I,J
182 COMPLEX*16 ZTEMP
183 real*8 denom(ncolor), cf(ncolor,ncolor)
184 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
185 COMPLEX*16 W(18,NWAVEFUNCS)
186 COMPLEX*16 DUM0,DUM1
187 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
188C
189C GLOBAL VARIABLES
190C
191 include 'coupl.inc'
192C
193C COLOR DATA
194C
195 DATA denom(1)/1/
196 DATA (cf(i, 1),i= 1, 2) / 9, 3/
197C 1 T(3,1) T(4,2)
198 DATA denom(2)/1/
199 DATA (cf(i, 2),i= 1, 2) / 3, 9/
200C 1 T(3,2) T(4,1)
201
202 INTEGER H1,H2
203 real*8 matrix
204 matrix_dd_dd_h=0.d0
205 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
206 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
207
208C ----------
209C BEGIN CODE
210C ----------
211 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
212 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
213 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
214 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
215 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
216 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
217 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
218 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,6))
219 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
220C Amplitude(s) for diagram number 1
221 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
222 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,5))
223 CALL ffv2_3_3(w(1,2),w(1,3),gc_50,gc_58,mz,wz,w(1,4))
224C Amplitude(s) for diagram number 2
225 CALL vvs1_0(w(1,5),w(1,4),w(1,7),gc_81,amp(2))
226 jamp(1)=-amp(1)
227 jamp(2)=+amp(2)
228
229 matrix = 0.d0
230 DO i = 1, ncolor
231 ztemp = (0.d0,0.d0)
232 DO j = 1, ncolor
233 ztemp = ztemp + cf(j,i)*jamp(j)
234 ENDDO
235 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
236 ENDDO
237 matrix_dd_dd_h=matrix
238 ENDIF ! CLOSES H1 IF
239 ENDIF ! CLOSES H2 IF
240 END
241C-----------------------end subprocess dd->ddH-------------
242
243
244
245
246c-----------------begin subprocess dd->jj_noH--------------
247c jj=dd only
248
249 SUBROUTINE dd_dd_noh(P,H1,H2,ANS)
250C
251C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
252C By the MadGraph Development Team
253C Please visit us at https://launchpad.net/madgraph5
254C
255C MadGraph StandAlone Version
256C
257C Returns amplitude squared summed/avg over colors
258C and helicities
259C for the point in phase space P(0:3,NEXTERNAL)
260C
261C Process: d d > d d ta+ ta- / h QED=4
262C
263 IMPLICIT NONE
264C
265C CONSTANTS
266C
267 INTEGER NEXTERNAL
268 parameter(nexternal=6)
269 INTEGER NCOMB
270 parameter( ncomb=64)
271C
272C ARGUMENTS
273C
274 real*8 p(0:3,nexternal),ans
275 INTEGER H1,H2
276C
277C LOCAL VARIABLES
278C
279 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
280 real*8 t
281 real*8 matrix_dd_dd_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/
352C ----------
353C BEGIN CODE
354C ----------
355 DO ihel=1,nexternal
356 jc(ihel) = +1
357 ENDDO
358 ans = 0d0
359 DO ihel=1,ncomb
360 t=matrix_dd_dd_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_dd_dd_noh(P,H1,H2,NHEL,IC)
368C
369C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
370C By the MadGraph Development Team
371C Please visit us at https://launchpad.net/madgraph5
372C
373C Returns amplitude squared summed/avg over colors
374C for the point with external lines W(0:6,NEXTERNAL)
375C
376C Process: d d > d d ta+ ta- / h QED=4
377C
378 IMPLICIT NONE
379C
380C CONSTANTS
381C
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))
392C
393C ARGUMENTS
394C
395 real*8 p(0:3,nexternal)
396 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
397C
398C LOCAL VARIABLES
399C
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)/
407C
408C GLOBAL VARIABLES
409C
410 include 'coupl.inc'
411C
412C COLOR DATA
413C
414 DATA denom(1)/1/
415 DATA (cf(i, 1),i= 1, 2) / 9, 3/
416C 1 T(3,1) T(4,2)
417 DATA denom(2)/1/
418 DATA (cf(i, 2),i= 1, 2) / 3, 9/
419C 1 T(3,2) T(4,1)
420
421 INTEGER H1,H2
422 real*8 matrix
423 matrix_dd_dd_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
427C ----------
428C BEGIN CODE
429C ----------
430 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
431 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
432 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
433 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
434 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
435 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
436 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,7))
437 CALL ffv1p0_3(w(1,2),w(1,4),gc_1,zero,zero,w(1,8))
438 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
439C Amplitude(s) for diagram number 1
440 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
441 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
442C Amplitude(s) for diagram number 2
443 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
444 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,11))
445C Amplitude(s) for diagram number 3
446 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
447C Amplitude(s) for diagram number 4
448 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
449 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,10))
450 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
451C Amplitude(s) for diagram number 5
452 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
453 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
454C Amplitude(s) for diagram number 6
455 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
456C Amplitude(s) for diagram number 7
457 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
458C Amplitude(s) for diagram number 8
459 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
460 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
461 CALL ffv1_2(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
462C Amplitude(s) for diagram number 9
463 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(9))
464 CALL ffv1_1(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
465C Amplitude(s) for diagram number 10
466 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_1,amp(10))
467 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
468C Amplitude(s) for diagram number 11
469 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(11))
470C Amplitude(s) for diagram number 12
471 CALL ffv2_3_0(w(1,2),w(1,13),w(1,7),gc_50,gc_58,amp(12))
472 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
473 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
474C Amplitude(s) for diagram number 13
475 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(13))
476 CALL ffv1_1(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
477C Amplitude(s) for diagram number 14
478 CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_1,amp(14))
479C Amplitude(s) for diagram number 15
480 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(15))
481C Amplitude(s) for diagram number 16
482 CALL ffv2_3_0(w(1,2),w(1,14),w(1,7),gc_50,gc_58,amp(16))
483 CALL ffv2_3_2(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
484C Amplitude(s) for diagram number 17
485 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_1,amp(17))
486 CALL ffv2_3_1(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
487C Amplitude(s) for diagram number 18
488 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
489C Amplitude(s) for diagram number 19
490 CALL ffv2_3_0(w(1,14),w(1,4),w(1,7),gc_50,gc_58,amp(19))
491C Amplitude(s) for diagram number 20
492 CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(20))
493 CALL ffv1p0_3(w(1,1),w(1,4),gc_1,zero,zero,w(1,9))
494 CALL ffv1p0_3(w(1,2),w(1,3),gc_1,zero,zero,w(1,14))
495 CALL ffv1_2(w(1,5),w(1,9),gc_3,mta,zero,w(1,10))
496C Amplitude(s) for diagram number 21
497 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(21))
498 CALL ffv1_1(w(1,6),w(1,9),gc_3,mta,zero,w(1,13))
499C Amplitude(s) for diagram number 22
500 CALL ffv1_0(w(1,5),w(1,13),w(1,14),gc_3,amp(22))
501 CALL ffv2_3_3(w(1,2),w(1,3),gc_50,gc_58,mz,wz,w(1,15))
502C Amplitude(s) for diagram number 23
503 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(23))
504C Amplitude(s) for diagram number 24
505 CALL ffv2_4_0(w(1,5),w(1,13),w(1,15),gc_50,gc_59,amp(24))
506 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,13))
507 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,10))
508C Amplitude(s) for diagram number 25
509 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(25))
510 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,16))
511C Amplitude(s) for diagram number 26
512 CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
513C Amplitude(s) for diagram number 27
514 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(27))
515C Amplitude(s) for diagram number 28
516 CALL ffv2_4_0(w(1,5),w(1,16),w(1,15),gc_50,gc_59,amp(28))
517 CALL ffv1_2(w(1,2),w(1,9),gc_1,zero,zero,w(1,16))
518C Amplitude(s) for diagram number 29
519 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_1,amp(29))
520 CALL ffv1_1(w(1,3),w(1,9),gc_1,zero,zero,w(1,5))
521C Amplitude(s) for diagram number 30
522 CALL ffv1_0(w(1,2),w(1,5),w(1,12),gc_1,amp(30))
523C Amplitude(s) for diagram number 31
524 CALL ffv2_3_0(w(1,16),w(1,3),w(1,7),gc_50,gc_58,amp(31))
525C Amplitude(s) for diagram number 32
526 CALL ffv2_3_0(w(1,2),w(1,5),w(1,7),gc_50,gc_58,amp(32))
527 CALL ffv1p0_3(w(1,1),w(1,4),gc_11,zero,zero,w(1,5))
528 CALL ffv1_2(w(1,2),w(1,5),gc_11,zero,zero,w(1,16))
529C Amplitude(s) for diagram number 33
530 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_1,amp(33))
531 CALL ffv1_1(w(1,3),w(1,5),gc_11,zero,zero,w(1,9))
532C Amplitude(s) for diagram number 34
533 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(34))
534C Amplitude(s) for diagram number 35
535 CALL ffv2_3_0(w(1,16),w(1,3),w(1,7),gc_50,gc_58,amp(35))
536C Amplitude(s) for diagram number 36
537 CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(36))
538 CALL ffv2_3_2(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,9))
539C Amplitude(s) for diagram number 37
540 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(37))
541 CALL ffv2_3_1(w(1,3),w(1,13),gc_50,gc_58,zero,zero,w(1,16))
542C Amplitude(s) for diagram number 38
543 CALL ffv1_0(w(1,2),w(1,16),w(1,12),gc_1,amp(38))
544C Amplitude(s) for diagram number 39
545 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(39))
546C Amplitude(s) for diagram number 40
547 CALL ffv2_3_0(w(1,2),w(1,16),w(1,7),gc_50,gc_58,amp(40))
548 CALL ffv1_2(w(1,1),w(1,14),gc_1,zero,zero,w(1,16))
549C Amplitude(s) for diagram number 41
550 CALL ffv1_0(w(1,16),w(1,4),w(1,12),gc_1,amp(41))
551 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,9))
552C Amplitude(s) for diagram number 42
553 CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_1,amp(42))
554C Amplitude(s) for diagram number 43
555 CALL ffv2_3_0(w(1,16),w(1,4),w(1,7),gc_50,gc_58,amp(43))
556 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,16))
557C Amplitude(s) for diagram number 44
558 CALL ffv1_0(w(1,16),w(1,4),w(1,14),gc_1,amp(44))
559 CALL ffv1p0_3(w(1,2),w(1,3),gc_11,zero,zero,w(1,14))
560 CALL ffv1_2(w(1,1),w(1,14),gc_11,zero,zero,w(1,13))
561C Amplitude(s) for diagram number 45
562 CALL ffv1_0(w(1,13),w(1,4),w(1,12),gc_1,amp(45))
563C Amplitude(s) for diagram number 46
564 CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_11,amp(46))
565C Amplitude(s) for diagram number 47
566 CALL ffv2_3_0(w(1,13),w(1,4),w(1,7),gc_50,gc_58,amp(47))
567C Amplitude(s) for diagram number 48
568 CALL ffv1_0(w(1,16),w(1,4),w(1,14),gc_11,amp(48))
569 CALL ffv2_3_2(w(1,1),w(1,15),gc_50,gc_58,zero,zero,w(1,14))
570C Amplitude(s) for diagram number 49
571 CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_1,amp(49))
572C Amplitude(s) for diagram number 50
573 CALL ffv2_3_0(w(1,9),w(1,4),w(1,15),gc_50,gc_58,amp(50))
574C Amplitude(s) for diagram number 51
575 CALL ffv2_3_0(w(1,14),w(1,4),w(1,7),gc_50,gc_58,amp(51))
576C Amplitude(s) for diagram number 52
577 CALL ffv2_3_0(w(1,16),w(1,4),w(1,15),gc_50,gc_58,amp(52))
578 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,15))
579C Amplitude(s) for diagram number 53
580 CALL ffv1_0(w(1,15),w(1,3),w(1,12),gc_1,amp(53))
581C Amplitude(s) for diagram number 54
582 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_1,amp(54))
583C Amplitude(s) for diagram number 55
584 CALL ffv2_3_0(w(1,15),w(1,3),w(1,7),gc_50,gc_58,amp(55))
585C Amplitude(s) for diagram number 56
586 CALL ffv1_0(w(1,16),w(1,3),w(1,8),gc_1,amp(56))
587 CALL ffv1p0_3(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
588 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
589C Amplitude(s) for diagram number 57
590 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_1,amp(57))
591C Amplitude(s) for diagram number 58
592 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(58))
593C Amplitude(s) for diagram number 59
594 CALL ffv2_3_0(w(1,4),w(1,3),w(1,7),gc_50,gc_58,amp(59))
595C Amplitude(s) for diagram number 60
596 CALL ffv1_0(w(1,16),w(1,3),w(1,8),gc_11,amp(60))
597 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
598C Amplitude(s) for diagram number 61
599 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_1,amp(61))
600C Amplitude(s) for diagram number 62
601 CALL ffv2_3_0(w(1,9),w(1,3),w(1,11),gc_50,gc_58,amp(62))
602C Amplitude(s) for diagram number 63
603 CALL ffv2_3_0(w(1,8),w(1,3),w(1,7),gc_50,gc_58,amp(63))
604C Amplitude(s) for diagram number 64
605 CALL ffv2_3_0(w(1,16),w(1,3),w(1,11),gc_50,gc_58,amp(64))
606 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
607 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
608 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
609 $ +1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)+1./2.*amp(36)
610 $ +1./2.*amp(45)+1./2.*amp(46)+1./2.*amp(47)+1./2.*amp(48)
611 $ -amp(53)-amp(54)-amp(55)-amp(56)+1./6.*amp(57)+1./6.*amp(58)
612 $ +1./6.*amp(59)+1./6.*amp(60)-amp(61)-amp(62)-amp(63)-amp(64)
613 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
614 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
615 $ +amp(28)+amp(29)+amp(30)+amp(31)+amp(32)-1./6.*amp(33)
616 $ -1./6.*amp(34)-1./6.*amp(35)-1./6.*amp(36)+amp(37)+amp(38)
617 $ +amp(39)+amp(40)+amp(41)+amp(42)+amp(43)+amp(44)-1./6.*amp(45)
618 $ -1./6.*amp(46)-1./6.*amp(47)-1./6.*amp(48)+amp(49)+amp(50)
619 $ +amp(51)+amp(52)-1./2.*amp(57)-1./2.*amp(58)-1./2.*amp(59)
620 $ -1./2.*amp(60)
621
622 matrix = 0.d0
623 DO i = 1, ncolor
624 ztemp = (0.d0,0.d0)
625 DO j = 1, ncolor
626 ztemp = ztemp + cf(j,i)*jamp(j)
627 ENDDO
628 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
629 ENDDO
630 matrix_dd_dd_noh=matrix
631 ENDIF
632 ENDIF
633 END
634C---------------------end subprocess initiated by DD-------