C++ Interface to Tauola
SD.f
1C JK. done 27.10.2015
2c------------------begin subprocess initiated by SD----------
3 SUBROUTINE sd(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 sd_ds_h(p,h1,h2,ans)
15 IF(i3.EQ.0 .AND. i4.EQ.0) CALL sd_ds_h(p,h1,h2,ans)
16 ELSE IF(key.EQ.0) THEN
17 IF(i3.EQ.1 .AND. i4.EQ.3) CALL sd_ds_noh(p,h1,h2,ans)
18 IF(i3.EQ.0 .AND. i4.EQ.0) CALL sd_ds_noh(p,h1,h2,ans)
19 ELSE
20 WRITE(*,*) 'spin=2 NOT FINISHED'
21 stop
22 ENDIF
23 END ! SUBROUTINE SD
24
25C ---------begin subprocess SD->jjH with H-> tautau
26c ---------------------jj=DS only -------------
27
28 SUBROUTINE sd_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: s d > 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_sd_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_sd_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_sd_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: s d > 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,2) T(4,1)
198
199 INTEGER H1,H2
200 real*8 matrix
201 matrix_sd_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,4),gc_50,gc_58,mz,wz,w(1,6))
216 CALL ffv2_3_3(w(1,2),w(1,3),gc_50,gc_58,mz,wz,w(1,4))
217C Amplitude(s) for diagram number 1
218 CALL vvs1_0(w(1,6),w(1,4),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_sd_ds_h=matrix
230 ENDIF ! CLOSES H1 IF
231 ENDIF ! CLOSES H2 IF
232 END
233C-----------------------end subprocess SD->ddH-------------
234
235
236
237
238c-----------------begin subprocess SD->jj_noH--------------
239c jj=DS only
240
241 SUBROUTINE sd_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: s d > 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_sd_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_sd_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_sd_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: s d > 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_sd_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
419C ----------
420C BEGIN CODE
421C ----------
422 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
423 CALL ixxxxx(p(0,2),zero,nhel(2),+1*ic(2),w(1,2))
424 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
425 CALL oxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
426 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
427 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
428 CALL ffv1p0_3(w(1,1),w(1,4),gc_1,zero,zero,w(1,7))
429 CALL ffv1p0_3(w(1,2),w(1,3),gc_1,zero,zero,w(1,8))
430 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
431C Amplitude(s) for diagram number 1
432 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
433 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
434C Amplitude(s) for diagram number 2
435 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
436 CALL ffv2_3_3(w(1,2),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
437C Amplitude(s) for diagram number 3
438 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
439C Amplitude(s) for diagram number 4
440 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
441 CALL ffv2_3_3(w(1,1),w(1,4),gc_50,gc_58,mz,wz,w(1,10))
442 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
443C Amplitude(s) for diagram number 5
444 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
445 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
446C Amplitude(s) for diagram number 6
447 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
448C Amplitude(s) for diagram number 7
449 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
450C Amplitude(s) for diagram number 8
451 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
452 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
453 CALL ffv1_2(w(1,2),w(1,7),gc_1,zero,zero,w(1,9))
454C Amplitude(s) for diagram number 9
455 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(9))
456 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,13))
457C Amplitude(s) for diagram number 10
458 CALL ffv1_0(w(1,2),w(1,13),w(1,12),gc_1,amp(10))
459 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
460C Amplitude(s) for diagram number 11
461 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(11))
462C Amplitude(s) for diagram number 12
463 CALL ffv2_3_0(w(1,2),w(1,13),w(1,7),gc_50,gc_58,amp(12))
464 CALL ffv1p0_3(w(1,1),w(1,4),gc_11,zero,zero,w(1,13))
465 CALL ffv1_2(w(1,2),w(1,13),gc_11,zero,zero,w(1,9))
466C Amplitude(s) for diagram number 13
467 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(13))
468 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,6))
469C Amplitude(s) for diagram number 14
470 CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(14))
471C Amplitude(s) for diagram number 15
472 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(15))
473C Amplitude(s) for diagram number 16
474 CALL ffv2_3_0(w(1,2),w(1,6),w(1,7),gc_50,gc_58,amp(16))
475 CALL ffv2_3_2(w(1,2),w(1,10),gc_50,gc_58,zero,zero,w(1,6))
476C Amplitude(s) for diagram number 17
477 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_1,amp(17))
478 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
479C Amplitude(s) for diagram number 18
480 CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
481C Amplitude(s) for diagram number 19
482 CALL ffv2_3_0(w(1,6),w(1,3),w(1,7),gc_50,gc_58,amp(19))
483C Amplitude(s) for diagram number 20
484 CALL ffv2_3_0(w(1,2),w(1,9),w(1,7),gc_50,gc_58,amp(20))
485 CALL ffv1_2(w(1,1),w(1,8),gc_1,zero,zero,w(1,9))
486C Amplitude(s) for diagram number 21
487 CALL ffv1_0(w(1,9),w(1,4),w(1,12),gc_1,amp(21))
488 CALL ffv1_2(w(1,1),w(1,12),gc_1,zero,zero,w(1,6))
489C Amplitude(s) for diagram number 22
490 CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_1,amp(22))
491C Amplitude(s) for diagram number 23
492 CALL ffv2_3_0(w(1,9),w(1,4),w(1,7),gc_50,gc_58,amp(23))
493 CALL ffv2_3_2(w(1,1),w(1,7),gc_50,gc_58,zero,zero,w(1,9))
494C Amplitude(s) for diagram number 24
495 CALL ffv1_0(w(1,9),w(1,4),w(1,8),gc_1,amp(24))
496 CALL ffv1p0_3(w(1,2),w(1,3),gc_11,zero,zero,w(1,8))
497 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,3))
498C Amplitude(s) for diagram number 25
499 CALL ffv1_0(w(1,3),w(1,4),w(1,12),gc_1,amp(25))
500C Amplitude(s) for diagram number 26
501 CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_11,amp(26))
502C Amplitude(s) for diagram number 27
503 CALL ffv2_3_0(w(1,3),w(1,4),w(1,7),gc_50,gc_58,amp(27))
504C Amplitude(s) for diagram number 28
505 CALL ffv1_0(w(1,9),w(1,4),w(1,8),gc_11,amp(28))
506 CALL ffv2_3_2(w(1,1),w(1,11),gc_50,gc_58,zero,zero,w(1,8))
507C Amplitude(s) for diagram number 29
508 CALL ffv1_0(w(1,8),w(1,4),w(1,12),gc_1,amp(29))
509C Amplitude(s) for diagram number 30
510 CALL ffv2_3_0(w(1,6),w(1,4),w(1,11),gc_50,gc_58,amp(30))
511C Amplitude(s) for diagram number 31
512 CALL ffv2_3_0(w(1,8),w(1,4),w(1,7),gc_50,gc_58,amp(31))
513C Amplitude(s) for diagram number 32
514 CALL ffv2_3_0(w(1,9),w(1,4),w(1,11),gc_50,gc_58,amp(32))
515 jamp(1)=+1./2.*(+amp(13)+amp(14)+amp(15)+amp(16)+amp(25)+amp(26)
516 $ +amp(27)+amp(28))
517 jamp(2)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
518 $ +amp(9)+amp(10)+amp(11)+amp(12)-1./6.*amp(13)-1./6.*amp(14)
519 $ -1./6.*amp(15)-1./6.*amp(16)+amp(17)+amp(18)+amp(19)+amp(20)
520 $ +amp(21)+amp(22)+amp(23)+amp(24)-1./6.*amp(25)-1./6.*amp(26)
521 $ -1./6.*amp(27)-1./6.*amp(28)+amp(29)+amp(30)+amp(31)+amp(32)
522
523 matrix = 0.d0
524 DO i = 1, ncolor
525 ztemp = (0.d0,0.d0)
526 DO j = 1, ncolor
527 ztemp = ztemp + cf(j,i)*jamp(j)
528 ENDDO
529 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
530 ENDDO
531 matrix_sd_ds_noh=matrix
532 ENDIF
533 ENDIF
534 END
535C---------------------end subprocess initiated by SD-------
536C DONE 27.10.2015