C++ Interface to Tauola
DS.f
1 C JK. done 27.10.2015
2 c------------------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 
9 C 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 
25 C ---------begin subprocess DS->jjH with H-> tautau
26 c ---------------------jj=DS only -------------
27 
28  SUBROUTINE ds_ds_h(P,H1,H2,ANS)
29 C
30 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
31 C By the MadGraph Development Team
32 C Please visit us at https://launchpad.net/madgraph5
33 C
34 C MadGraph StandAlone Version
35 C
36 C Returns amplitude squared summed/avg over colors
37 C and helicities
38 C for the point in phase space P(0:3,NEXTERNAL)
39 C
40 C Process: d s > d s h WEIGHTED=6
41 C * Decay: h > ta+ ta- WEIGHTED=2
42 C
43  IMPLICIT NONE
44 C
45 C CONSTANTS
46 C
47  INTEGER NEXTERNAL
48  parameter(nexternal=6)
49  INTEGER NCOMB
50  parameter( ncomb=64)
51 C
52 C ARGUMENTS
53 C
54  real*8 p(0:3,nexternal),ans
55  INTEGER H1,H2
56 C
57 C LOCAL VARIABLES
58 C
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/
132 C ----------
133 C BEGIN CODE
134 C ----------
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)
148 C
149 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
150 C By the MadGraph Development Team
151 C Please visit us at https://launchpad.net/madgraph5
152 C
153 C Returns amplitude squared summed/avg over colors
154 C for the point with external lines W(0:6,NEXTERNAL)
155 C
156 C Process: d s > d s h WEIGHTED=6
157 C * Decay: h > ta+ ta- WEIGHTED=2
158 C
159  IMPLICIT NONE
160 C
161 C CONSTANTS
162 C
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))
173 C
174 C ARGUMENTS
175 C
176  real*8 p(0:3,nexternal)
177  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
178 C
179 C LOCAL VARIABLES
180 C
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)/
188 C
189 C GLOBAL VARIABLES
190 C
191  include 'coupl.inc'
192 C
193 C COLOR DATA
194 C
195  DATA denom(1)/1/
196  DATA (cf(i, 1),i= 1, 1) / 9/
197 C 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 
205 C ----------
206 C BEGIN CODE
207 C ----------
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))
217 C 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
233 C-----------------------end subprocess DS->ddH-------------
234 
235 
236 
237 
238 c-----------------begin subprocess DS->jj_noH--------------
239 c jj=dd only
240 
241  SUBROUTINE ds_ds_noh(P,H1,H2,ANS)
242 C
243 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
244 C By the MadGraph Development Team
245 C Please visit us at https://launchpad.net/madgraph5
246 C
247 C MadGraph StandAlone Version
248 C
249 C Returns amplitude squared summed/avg over colors
250 C and helicities
251 C for the point in phase space P(0:3,NEXTERNAL)
252 C
253 C Process: d s > d s ta+ ta- / h QED=4
254 C
255  IMPLICIT NONE
256 C
257 C CONSTANTS
258 C
259  INTEGER NEXTERNAL
260  parameter(nexternal=6)
261  INTEGER NCOMB
262  parameter( ncomb=64)
263 C
264 C ARGUMENTS
265 C
266  real*8 p(0:3,nexternal),ans
267  INTEGER H1,H2
268 C
269 C LOCAL VARIABLES
270 C
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/
344 C ----------
345 C BEGIN CODE
346 C ----------
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)
360 C
361 C Generated by MadGraph 5 v. 1.5.15, 2013-12-11
362 C By the MadGraph Development Team
363 C Please visit us at https://launchpad.net/madgraph5
364 C
365 C Returns amplitude squared summed/avg over colors
366 C for the point with external lines W(0:6,NEXTERNAL)
367 C
368 C Process: d s > d s ta+ ta- / h QED=4
369 C
370  IMPLICIT NONE
371 C
372 C CONSTANTS
373 C
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))
384 C
385 C ARGUMENTS
386 C
387  real*8 p(0:3,nexternal)
388  INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
389 C
390 C LOCAL VARIABLES
391 C
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)/
399 C
400 C GLOBAL VARIABLES
401 C
402  include 'coupl.inc'
403 C
404 C COLOR DATA
405 C
406  DATA denom(1)/1/
407  DATA (cf(i, 1),i= 1, 2) / 9, 3/
408 C 1 T(3,1) T(4,2)
409  DATA denom(2)/1/
410  DATA (cf(i, 2),i= 1, 2) / 3, 9/
411 C 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 
420 C ----------
421 C BEGIN CODE
422 C ----------
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))
432 C 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))
435 C 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))
438 C 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))
440 C 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))
444 C 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))
447 C Amplitude(s) for diagram number 6
448  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
449 C 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))
451 C 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))
455 C 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))
458 C 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))
461 C 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))
463 C 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))
467 C 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))
470 C Amplitude(s) for diagram number 14
471  CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(14))
472 C 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))
474 C 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))
477 C 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))
480 C Amplitude(s) for diagram number 18
481  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
482 C 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))
484 C 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))
487 C 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))
490 C Amplitude(s) for diagram number 22
491  CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_1,amp(22))
492 C 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))
495 C 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))
499 C Amplitude(s) for diagram number 25
500  CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_1,amp(25))
501 C Amplitude(s) for diagram number 26
502  CALL ffv1_0(w(1,6),w(1,3),w(1,8),gc_11,amp(26))
503 C 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))
505 C 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))
508 C Amplitude(s) for diagram number 29
509  CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_1,amp(29))
510 C 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))
512 C 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))
514 C 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
536 C---------------------end subprocess initiated by DS-------
537 C DONE 27.10.2015
Generated by doxygen 1.9.1