1
SUBROUTINE SMATRIX_3(P,ANS)
3
C Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
4
C By the MadGraph5_aMC@NLO Development Team
5
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
7
C Returns amplitude squared summed/avg over colors
9
C for the point in phase space P(0:3,NEXTERNAL)
11
C Process: d~ g > w+ u~ WEIGHTED<=3 [ all = QCD ]
12
C Process: s~ g > w+ c~ WEIGHTED<=3 [ all = QCD ]
18
INCLUDE 'nexternal.inc'
24
REAL*8 P(0:3,NEXTERNAL),ANS
25
DOUBLE PRECISION WGT_ME_BORN,WGT_ME_REAL
26
COMMON /C_WGT_ME_TREE/ WGT_ME_BORN,WGT_ME_REAL
30
INTEGER IHEL,IDEN,I,T_IDENT(NCOMB)
32
REAL*8 T,T_SAVE(NCOMB)
34
INTEGER NHEL(NEXTERNAL,NCOMB)
35
DATA (NHEL(I, 1),I=1,4) /-1,-1,-1, 1/
36
DATA (NHEL(I, 2),I=1,4) /-1,-1,-1,-1/
37
DATA (NHEL(I, 3),I=1,4) /-1,-1, 0, 1/
38
DATA (NHEL(I, 4),I=1,4) /-1,-1, 0,-1/
39
DATA (NHEL(I, 5),I=1,4) /-1,-1, 1, 1/
40
DATA (NHEL(I, 6),I=1,4) /-1,-1, 1,-1/
41
DATA (NHEL(I, 7),I=1,4) /-1, 1,-1, 1/
42
DATA (NHEL(I, 8),I=1,4) /-1, 1,-1,-1/
43
DATA (NHEL(I, 9),I=1,4) /-1, 1, 0, 1/
44
DATA (NHEL(I, 10),I=1,4) /-1, 1, 0,-1/
45
DATA (NHEL(I, 11),I=1,4) /-1, 1, 1, 1/
46
DATA (NHEL(I, 12),I=1,4) /-1, 1, 1,-1/
47
DATA (NHEL(I, 13),I=1,4) / 1,-1,-1, 1/
48
DATA (NHEL(I, 14),I=1,4) / 1,-1,-1,-1/
49
DATA (NHEL(I, 15),I=1,4) / 1,-1, 0, 1/
50
DATA (NHEL(I, 16),I=1,4) / 1,-1, 0,-1/
51
DATA (NHEL(I, 17),I=1,4) / 1,-1, 1, 1/
52
DATA (NHEL(I, 18),I=1,4) / 1,-1, 1,-1/
53
DATA (NHEL(I, 19),I=1,4) / 1, 1,-1, 1/
54
DATA (NHEL(I, 20),I=1,4) / 1, 1,-1,-1/
55
DATA (NHEL(I, 21),I=1,4) / 1, 1, 0, 1/
56
DATA (NHEL(I, 22),I=1,4) / 1, 1, 0,-1/
57
DATA (NHEL(I, 23),I=1,4) / 1, 1, 1, 1/
58
DATA (NHEL(I, 24),I=1,4) / 1, 1, 1,-1/
59
LOGICAL GOODHEL(NCOMB)
60
DATA GOODHEL/NCOMB*.FALSE./
70
IF (GOODHEL(IHEL) .OR. NTRY .LT. 2) THEN
72
C for the first ps-point, check for helicities that give
73
C identical matrix elements
74
T=MATRIX_3(P ,NHEL(1,IHEL))
79
IF (T_SAVE(I).EQ.0D0) CYCLE
80
IF (ABS(T/T_SAVE(I)-1D0) .LT. 1D-12) THEN
81
C WRITE (*,*) 'FOUND IDENTICAL',T,IHEL,T_SAVE(I),I
86
IF (T_IDENT(IHEL).GT.0) THEN
87
C if two helicity states are identical, dont recompute
88
T=T_SAVE(T_IDENT(IHEL))
91
T=MATRIX_3(P ,NHEL(1,IHEL))
95
C add to the sum of helicities
97
IF (T .NE. 0D0 .AND. .NOT. GOODHEL(IHEL)) THEN
107
REAL*8 FUNCTION MATRIX_3(P,NHEL)
109
C Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
110
C By the MadGraph5_aMC@NLO Development Team
111
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
113
C Returns amplitude squared summed/avg over colors
114
C for the point with external lines W(0:6,NEXTERNAL)
116
C Process: d~ g > w+ u~ WEIGHTED<=3 [ all = QCD ]
117
C Process: s~ g > w+ c~ WEIGHTED<=3 [ all = QCD ]
124
PARAMETER (NGRAPHS=2)
125
INTEGER NWAVEFUNCS, NCOLOR
126
PARAMETER (NWAVEFUNCS=5, NCOLOR=1)
130
PARAMETER (IMAG1=(0D0,1D0))
131
INCLUDE 'nexternal.inc'
136
REAL*8 P(0:3,NEXTERNAL)
137
INTEGER NHEL(NEXTERNAL)
142
INTEGER IC(NEXTERNAL)
143
DATA IC /NEXTERNAL*1/
144
REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
145
COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR), W(8,NWAVEFUNCS)
150
DATA (CF(I, 1),I= 1, 1) / 4/
155
CALL OXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
156
CALL VXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
157
CALL VXXXXX(P(0,3),MDL_MW,NHEL(3),+1*IC(3),W(1,3))
158
CALL IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
159
CALL FFV1_1(W(1,1),W(1,2),GC_5,ZERO,ZERO,W(1,5))
160
C Amplitude(s) for diagram number 1
161
CALL FFV2_0(W(1,4),W(1,5),W(1,3),GC_47,AMP(1))
162
CALL FFV2_1(W(1,1),W(1,3),GC_47,ZERO,ZERO,W(1,5))
163
C Amplitude(s) for diagram number 2
164
CALL FFV1_0(W(1,4),W(1,5),W(1,2),GC_5,AMP(2))
165
JAMP(1)=-AMP(1)-AMP(2)
170
ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
172
MATRIX_3 = MATRIX_3+ZTEMP*DCONJG(JAMP(I))/DENOM(I)