C++ Interface to Tauola
DS.f
1C JK. done 27.10.2015
2c------------------begin subprocess initiated by DS----------
3 SUBROUTINE ds(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,3 ! only DS->DS
10
11 ans=0.d0
12
13 IF(key.EQ.1) THEN
14 IF(i3.EQ.1 .AND. i4.EQ.3) CALL ds_ds_h(p,h1,h2,ans)
15 IF(i3.EQ.0 .AND. i4.EQ.0) CALL ds_ds_h(p,h1,h2,ans)
16 ELSE IF(key.EQ.0) THEN
17 IF(i3.EQ.1 .AND. i4.EQ.3) CALL ds_ds_noh(p,h1,h2,ans)
18 IF(i3.EQ.0 .AND. i4.EQ.0) CALL ds_ds_noh(p,h1,h2,ans)
19 ELSE
20 WRITE(*,*) 'spin=2 NOT FINISHED'
21 stop
22 ENDIF
23 END ! SUBROUTINE DS
24
25C ---------begin subprocess DS->jjH with H-> tautau
26c ---------------------jj=DS only -------------
27
28 SUBROUTINE ds_ds_h(P,H1,H2,ANS)
29C
30C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
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 s > d s 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_ds_ds_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/36/
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_ds_ds_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_ds_ds_h(P,H1,H2,NHEL,IC)
148C
149C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
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 s > d s h WEIGHTED=6
157C * Decay: h > ta+ ta- WEIGHTED=2
158C
159 IMPLICIT NONE
160C
161C CONSTANTS
162C
163 INTEGER NGRAPHS
164 parameter(ngraphs=1)
165 INTEGER NEXTERNAL
166 parameter(nexternal=6)
167 INTEGER NWAVEFUNCS, NCOLOR
168 parameter(nwavefuncs=7, ncolor=1)
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, 1) / 9/
197C 1 T(3,1) T(4,2)
198
199 INTEGER H1,H2
200 real*8 matrix
201 matrix_ds_ds_h=0.d0
202 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
203 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
204
205C ----------
206C BEGIN CODE
207C ----------
208 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
209 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
210 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
211 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
212 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
213 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
214 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
215 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,6))
216 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,3))
217C Amplitude(s) for diagram number 1
218 CALL vvs1_0(w(1,6),w(1,3),w(1,7),gc_81,amp(1))
219 jamp(1)=-amp(1)
220
221 matrix = 0.d0
222 DO i = 1, ncolor
223 ztemp = (0.d0,0.d0)
224 DO j = 1, ncolor
225 ztemp = ztemp + cf(j,i)*jamp(j)
226 ENDDO
227 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
228 ENDDO
229 matrix_ds_ds_h=matrix
230 ENDIF ! CLOSES H1 IF
231 ENDIF ! CLOSES H2 IF
232 END
233C-----------------------end subprocess DS->ddH-------------
234
235
236
237
238c-----------------begin subprocess DS->jj_noH--------------
239c jj=dd only
240
241 SUBROUTINE ds_ds_noh(P,H1,H2,ANS)
242C
243C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
244C By the MadGraph Development Team
245C Please visit us at https://launchpad.net/madgraph5
246C
247C MadGraph StandAlone Version
248C
249C Returns amplitude squared summed/avg over colors
250C and helicities
251C for the point in phase space P(0:3,NEXTERNAL)
252C
253C Process: d s > d s ta+ ta- / h QED=4
254C
255 IMPLICIT NONE
256C
257C CONSTANTS
258C
259 INTEGER NEXTERNAL
260 parameter(nexternal=6)
261 INTEGER NCOMB
262 parameter( ncomb=64)
263C
264C ARGUMENTS
265C
266 real*8 p(0:3,nexternal),ans
267 INTEGER H1,H2
268C
269C LOCAL VARIABLES
270C
271 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
272 real*8 t
273 real*8 matrix_ds_ds_noh
274 INTEGER IHEL,IDEN, I
275 INTEGER JC(NEXTERNAL)
276 LOGICAL GOODHEL(NCOMB)
277 DATA ntry/0/
278 DATA goodhel/ncomb*.false./
279 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
280 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
281 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
282 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
283 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
284 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
285 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
286 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
287 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
288 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
289 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
290 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
291 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
292 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
293 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
294 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
295 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
296 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
297 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
298 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
299 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
300 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
301 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
302 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
303 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
304 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
305 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
306 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
307 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
308 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
309 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
310 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
311 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
312 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
313 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
314 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
315 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
316 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
317 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
318 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
319 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
320 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
321 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
322 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
323 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
324 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
325 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
326 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
327 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
328 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
329 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
330 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
331 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
332 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
333 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
334 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
335 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
336 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
337 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
338 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
339 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
340 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
341 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
342 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
343 DATA iden/36/
344C ----------
345C BEGIN CODE
346C ----------
347 DO ihel=1,nexternal
348 jc(ihel) = +1
349 ENDDO
350 ans = 0d0
351 DO ihel=1,ncomb
352 t=matrix_ds_ds_noh(p ,h1,h2,nhel(1,ihel),jc(1))
353 ans=ans+t
354 ENDDO
355 ans=ans/dble(iden)
356 END
357
358
359 REAL*8 FUNCTION matrix_ds_ds_noh(P,H1,H2,NHEL,IC)
360C
361C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
362C By the MadGraph Development Team
363C Please visit us at https://launchpad.net/madgraph5
364C
365C Returns amplitude squared summed/avg over colors
366C for the point with external lines W(0:6,NEXTERNAL)
367C
368C Process: d s > d s ta+ ta- / h QED=4
369C
370 IMPLICIT NONE
371C
372C CONSTANTS
373C
374 INTEGER NGRAPHS
375 parameter(ngraphs=32)
376 INTEGER NEXTERNAL
377 parameter(nexternal=6)
378 INTEGER NWAVEFUNCS, NCOLOR
379 parameter(nwavefuncs=13, ncolor=2)
380 real*8 zero
381 parameter(zero=0d0)
382 COMPLEX*16 IMAG1
383 parameter(imag1=(0d0,1d0))
384C
385C ARGUMENTS
386C
387 real*8 p(0:3,nexternal)
388 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
389C
390C LOCAL VARIABLES
391C
392 INTEGER I,J
393 COMPLEX*16 ZTEMP
394 real*8 denom(ncolor), cf(ncolor,ncolor)
395 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
396 COMPLEX*16 W(18,NWAVEFUNCS)
397 COMPLEX*16 DUM0,DUM1
398 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
399C
400C GLOBAL VARIABLES
401C
402 include 'coupl.inc'
403C
404C COLOR DATA
405C
406 DATA denom(1)/1/
407 DATA (cf(i, 1),i= 1, 2) / 9, 3/
408C 1 T(3,1) T(4,2)
409 DATA denom(2)/1/
410 DATA (cf(i, 2),i= 1, 2) / 3, 9/
411C 1 T(3,2) T(4,1)
412
413 INTEGER H1,H2
414 real*8 matrix
415 matrix_ds_ds_noh=0.d0
416 IF(h1.EQ.0. or .h1.EQ.nhel(5)) THEN
417 IF(h2.EQ.0. or .h2.EQ.nhel(6)) THEN
418
419
420C ----------
421C BEGIN CODE
422C ----------
423 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
424 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
425 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
426 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
427 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
428 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
429 CALL ffv1p0_3(w(1,1),w(1,3),gc_1,zero,zero,w(1,7))
430 CALL ffv1p0_3(w(1,2),w(1,4),gc_1,zero,zero,w(1,8))
431 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
432C Amplitude(s) for diagram number 1
433 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
434 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
435C Amplitude(s) for diagram number 2
436 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
437 CALL ffv2_3_3(w(1,2),w(1,4),gc_50,gc_58,mz,wz,w(1,11))
438C Amplitude(s) for diagram number 3
439 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
440C Amplitude(s) for diagram number 4
441 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
442 CALL ffv2_3_3(w(1,1),w(1,3),gc_50,gc_58,mz,wz,w(1,10))
443 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
444C Amplitude(s) for diagram number 5
445 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
446 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
447C Amplitude(s) for diagram number 6
448 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
449C Amplitude(s) for diagram number 7
450 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
451C Amplitude(s) for diagram number 8
452 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
453 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
454 CALL ffv1_2(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
455C Amplitude(s) for diagram number 9
456 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(9))
457 CALL ffv1_1(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
458C Amplitude(s) for diagram number 10
459 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_1,amp(10))
460 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
461C Amplitude(s) for diagram number 11
462 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(11))
463C Amplitude(s) for diagram number 12
464 CALL ffv2_3_0(w(1,2),w(1,13),w(1,7),gc_50,gc_58,amp(12))
465 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,13))
466 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
467C Amplitude(s) for diagram number 13
468 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(13))
469 CALL ffv1_1(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
470C Amplitude(s) for diagram number 14
471 CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(14))
472C Amplitude(s) for diagram number 15
473 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(15))
474C Amplitude(s) for diagram number 16
475 CALL ffv2_3_0(w(1,2),w(1,6),w(1,7),gc_50,gc_58,amp(16))
476 CALL ffv2_3_2(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
477C Amplitude(s) for diagram number 17
478 CALL ffv1_0(w(1,6),w(1,4),w(1,12),gc_1,amp(17))
479 CALL ffv2_3_1(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
480C Amplitude(s) for diagram number 18
481 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
482C Amplitude(s) for diagram number 19
483 CALL ffv2_3_0(w(1,6),w(1,4),w(1,7),gc_50,gc_58,amp(19))
484C Amplitude(s) for diagram number 20
485 CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(20))
486 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
487C Amplitude(s) for diagram number 21
488 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(21))
489 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
490C Amplitude(s) for diagram number 22
491 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_1,amp(22))
492C Amplitude(s) for diagram number 23
493 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(23))
494 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
495C Amplitude(s) for diagram number 24
496 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_1,amp(24))
497 CALL ffv1p0_3(w(1,2),w(1,4),gc_11,zero,zero,w(1,8))
498 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
499C Amplitude(s) for diagram number 25
500 CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_1,amp(25))
501C Amplitude(s) for diagram number 26
502 CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(26))
503C Amplitude(s) for diagram number 27
504 CALL ffv2_3_0(w(1,4),w(1,3),w(1,7),gc_50,gc_58,amp(27))
505C Amplitude(s) for diagram number 28
506 CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(28))
507 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
508C Amplitude(s) for diagram number 29
509 CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_1,amp(29))
510C Amplitude(s) for diagram number 30
511 CALL ffv2_3_0(w(1,6),w(1,3),w(1,11),gc_50,gc_58,amp(30))
512C Amplitude(s) for diagram number 31
513 CALL ffv2_3_0(w(1,8),w(1,3),w(1,7),gc_50,gc_58,amp(31))
514C Amplitude(s) for diagram number 32
515 CALL ffv2_3_0(w(1,9),w(1,3),w(1,11),gc_50,gc_58,amp(32))
516 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
517 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
518 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
519 $ -amp(21)-amp(22)-amp(23)-amp(24)+1./6.*amp(25)+1./6.*amp(26)
520 $ +1./6.*amp(27)+1./6.*amp(28)-amp(29)-amp(30)-amp(31)-amp(32)
521 jamp(2)=+1./2.*(-amp(13)-amp(14)-amp(15)-amp(16)-amp(25)-amp(26)
522 $ -amp(27)-amp(28))
523
524 matrix = 0.d0
525 DO i = 1, ncolor
526 ztemp = (0.d0,0.d0)
527 DO j = 1, ncolor
528 ztemp = ztemp + cf(j,i)*jamp(j)
529 ENDDO
530 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
531 ENDDO
532 matrix_ds_ds_noh=matrix
533 ENDIF
534 ENDIF
535 END
536C---------------------end subprocess initiated by DS-------
537C DONE 27.10.2015