C++InterfacetoTauola
DD.f
1 C JK. done 30.07.2014
2 c------------------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 
9 C 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 
25 C ---------begin subprocess DD->jjH with H-> tautau
26 c ---------------------jj=dd only -------------
27 
28  SUBROUTINE dd_dd_h(P,H1,H2,ANS)
29 C
30 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
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 d > d d 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_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/
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_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)
148 C
149 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
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 d > d d 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=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))
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, 2) / 9, 3/
197 C 1 T(3,1) T(4,2)
198  DATA denom(2)/1/
199  DATA (cf(i, 2),i= 1, 2) / 3, 9/
200 C 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 
208 C ----------
209 C BEGIN CODE
210 C ----------
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))
220 C 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))
224 C 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
241 C-----------------------end subprocess dd->ddH-------------
242 
243 
244 
245 
246 c-----------------begin subprocess dd->jj_noH--------------
247 c jj=dd only
248 
249  SUBROUTINE dd_dd_noh(P,H1,H2,ANS)
250 C
251 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
252 C By the MadGraph Development Team
253 C Please visit us at https://launchpad.net/madgraph5
254 C
255 C MadGraph StandAlone Version
256 C
257 C Returns amplitude squared summed/avg over colors
258 C and helicities
259 C for the point in phase space P(0:3,NEXTERNAL)
260 C
261 C Process: d d > d d ta+ ta- / h QED=4
262 C
263  IMPLICIT NONE
264 C
265 C CONSTANTS
266 C
267  INTEGER nexternal
268  parameter(nexternal=6)
269  INTEGER ncomb
270  parameter( ncomb=64)
271 C
272 C ARGUMENTS
273 C
274  REAL*8 p(0:3,nexternal),ans
275  INTEGER h1,h2
276 C
277 C LOCAL VARIABLES
278 C
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/
352 C ----------
353 C BEGIN CODE
354 C ----------
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)
368 C
369 C Generated by MadGraph 5 v. 1.5.12, 2013-08-21
370 C By the MadGraph Development Team
371 C Please visit us at https://launchpad.net/madgraph5
372 C
373 C Returns amplitude squared summed/avg over colors
374 C for the point with external lines W(0:6,NEXTERNAL)
375 C
376 C Process: d d > d d ta+ ta- / h QED=4
377 C
378  IMPLICIT NONE
379 C
380 C CONSTANTS
381 C
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))
392 C
393 C ARGUMENTS
394 C
395  REAL*8 p(0:3,nexternal)
396  INTEGER nhel(nexternal), ic(nexternal)
397 C
398 C LOCAL VARIABLES
399 C
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)/
407 C
408 C GLOBAL VARIABLES
409 C
410  include 'coupl.inc'
411 C
412 C COLOR DATA
413 C
414  DATA denom(1)/1/
415  DATA (cf(i, 1),i= 1, 2) / 9, 3/
416 C 1 T(3,1) T(4,2)
417  DATA denom(2)/1/
418  DATA (cf(i, 2),i= 1, 2) / 3, 9/
419 C 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 
427 C ----------
428 C BEGIN CODE
429 C ----------
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))
439 C 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))
442 C 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))
445 C 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))
447 C 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))
451 C 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))
454 C Amplitude(s) for diagram number 6
455  CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
456 C 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))
458 C 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))
462 C 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))
465 C 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))
468 C 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))
470 C 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))
474 C 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))
477 C Amplitude(s) for diagram number 14
478  CALL ffv1_0(w(1,2),w(1,14),w(1,12),gc_1,amp(14))
479 C 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))
481 C 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))
484 C 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))
487 C Amplitude(s) for diagram number 18
488  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(18))
489 C 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))
491 C 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))
496 C 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))
499 C 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))
502 C 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))
504 C 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))
508 C 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))
511 C Amplitude(s) for diagram number 26
512  CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
513 C 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))
515 C 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))
518 C 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))
521 C Amplitude(s) for diagram number 30
522  CALL ffv1_0(w(1,2),w(1,5),w(1,12),gc_1,amp(30))
523 C 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))
525 C 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))
529 C 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))
532 C Amplitude(s) for diagram number 34
533  CALL ffv1_0(w(1,2),w(1,9),w(1,12),gc_1,amp(34))
534 C 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))
536 C 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))
539 C 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))
542 C Amplitude(s) for diagram number 38
543  CALL ffv1_0(w(1,2),w(1,16),w(1,12),gc_1,amp(38))
544 C 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))
546 C 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))
549 C 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))
552 C Amplitude(s) for diagram number 42
553  CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_1,amp(42))
554 C 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))
557 C 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))
561 C Amplitude(s) for diagram number 45
562  CALL ffv1_0(w(1,13),w(1,4),w(1,12),gc_1,amp(45))
563 C Amplitude(s) for diagram number 46
564  CALL ffv1_0(w(1,9),w(1,4),w(1,14),gc_11,amp(46))
565 C 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))
567 C 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))
570 C Amplitude(s) for diagram number 49
571  CALL ffv1_0(w(1,14),w(1,4),w(1,12),gc_1,amp(49))
572 C 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))
574 C 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))
576 C 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))
579 C Amplitude(s) for diagram number 53
580  CALL ffv1_0(w(1,15),w(1,3),w(1,12),gc_1,amp(53))
581 C Amplitude(s) for diagram number 54
582  CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_1,amp(54))
583 C 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))
585 C 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))
589 C Amplitude(s) for diagram number 57
590  CALL ffv1_0(w(1,4),w(1,3),w(1,12),gc_1,amp(57))
591 C Amplitude(s) for diagram number 58
592  CALL ffv1_0(w(1,9),w(1,3),w(1,8),gc_11,amp(58))
593 C 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))
595 C 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))
598 C Amplitude(s) for diagram number 61
599  CALL ffv1_0(w(1,8),w(1,3),w(1,12),gc_1,amp(61))
600 C 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))
602 C 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))
604 C 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
634 C---------------------end subprocess initiated by DD-------
Generated on Tue Mar 12 2024 22:57:50 for C++InterfacetoTauola by   doxygen 1.8.7