C++InterfacetoTauola
SD.f
1 C JK. done 27.10.2015
2 c------------------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 
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 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 
25 C ---------begin subprocess SD->jjH with H-> tautau
26 c ---------------------jj=DS only -------------
27 
28  SUBROUTINE sd_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: s d > 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_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/
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_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)
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: s d > 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,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 
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,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))
217 C 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
233 C-----------------------end subprocess SD->ddH-------------
234 
235 
236 
237 
238 c-----------------begin subprocess SD->jj_noH--------------
239 c jj=DS only
240 
241  SUBROUTINE sd_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: s d > 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_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/
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_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)
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: s d > 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_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 
419 C ----------
420 C BEGIN CODE
421 C ----------
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))
431 C 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))
434 C 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))
437 C 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))
439 C 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))
443 C 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))
446 C Amplitude(s) for diagram number 6
447  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
448 C 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))
450 C 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))
454 C 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))
457 C 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))
460 C 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))
462 C 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))
466 C 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))
469 C Amplitude(s) for diagram number 14
470  CALL ffv1_0(w(1,2),w(1,6),w(1,12),gc_1,amp(14))
471 C 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))
473 C 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))
476 C 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))
479 C Amplitude(s) for diagram number 18
480  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
481 C 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))
483 C 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))
486 C 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))
489 C Amplitude(s) for diagram number 22
490  CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_1,amp(22))
491 C 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))
494 C 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))
498 C Amplitude(s) for diagram number 25
499  CALL ffv1_0(w(1,3),w(1,4),w(1,12),gc_1,amp(25))
500 C Amplitude(s) for diagram number 26
501  CALL ffv1_0(w(1,6),w(1,4),w(1,8),gc_11,amp(26))
502 C 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))
504 C 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))
507 C Amplitude(s) for diagram number 29
508  CALL ffv1_0(w(1,8),w(1,4),w(1,12),gc_1,amp(29))
509 C 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))
511 C 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))
513 C 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
535 C---------------------end subprocess initiated by SD-------
536 C DONE 27.10.2015
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7