1
SUBROUTINE SMATRIX5(P,ANS)
3
C Generated by MadGraph 5 v. 2.0.0.beta3, 2013-02-14
4
C By the MadGraph Development Team
5
C Please visit us at https://launchpad.net/madgraph5
7
C MadGraph for Madevent Version
9
C Returns amplitude squared summed/avg over colors
11
C for the point in phase space P(0:3,NEXTERNAL)
13
C Process: u d~ > u d~ WEIGHTED=2
14
C Process: u s~ > u s~ WEIGHTED=2
15
C Process: s d~ > s d~ WEIGHTED=2
16
C Process: s u~ > s u~ WEIGHTED=2
23
INCLUDE 'maxconfigs.inc'
24
INCLUDE 'nexternal.inc'
33
PARAMETER (THEL=2*NCOMB)
37
REAL*8 P(0:3,NEXTERNAL),ANS
41
INTEGER NHEL(NEXTERNAL,NCOMB),NTRY(2)
44
REAL*8 R,SUMHEL,TS(NCOMB)
46
INTEGER JC(NEXTERNAL),II
47
LOGICAL GOODHEL(NCOMB,2)
48
REAL*8 HWGT, XTOT, XTRY, XREJ, XR, YFRAC(0:NCOMB)
49
INTEGER NGOOD(2), IGOOD(NCOMB,2)
50
INTEGER JHEL(2), J, JJ
54
DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
55
COMMON/TO_AMPS/ AMP2, JAMP2
57
CHARACTER*101 HEL_BUFF
58
COMMON/TO_HELICITY/ HEL_BUFF
61
COMMON/TO_MIRROR/ IMIRROR
64
COMMON/TO_POLARIZATION/ POL
68
COMMON/TO_MATRIX/ISUM_HEL, MULTI_CHANNEL
69
INTEGER MAPCONFIG(0:LMAXCONFIGS), ICONFIG
70
COMMON/TO_MCONFIGS/MAPCONFIG, ICONFIG
71
INTEGER SUBDIAG(MAXSPROC),IB(2)
72
COMMON/TO_SUB_DIAG/SUBDIAG,IB
77
SAVE YFRAC, IGOOD, JHEL
78
DATA GOODHEL/THEL*.FALSE./
79
DATA (NHEL(I, 1),I=1,4) /-1,-1,-1,-1/
80
DATA (NHEL(I, 2),I=1,4) /-1,-1,-1, 1/
81
DATA (NHEL(I, 3),I=1,4) /-1,-1, 1,-1/
82
DATA (NHEL(I, 4),I=1,4) /-1,-1, 1, 1/
83
DATA (NHEL(I, 5),I=1,4) /-1, 1,-1,-1/
84
DATA (NHEL(I, 6),I=1,4) /-1, 1,-1, 1/
85
DATA (NHEL(I, 7),I=1,4) /-1, 1, 1,-1/
86
DATA (NHEL(I, 8),I=1,4) /-1, 1, 1, 1/
87
DATA (NHEL(I, 9),I=1,4) / 1,-1,-1,-1/
88
DATA (NHEL(I, 10),I=1,4) / 1,-1,-1, 1/
89
DATA (NHEL(I, 11),I=1,4) / 1,-1, 1,-1/
90
DATA (NHEL(I, 12),I=1,4) / 1,-1, 1, 1/
91
DATA (NHEL(I, 13),I=1,4) / 1, 1,-1,-1/
92
DATA (NHEL(I, 14),I=1,4) / 1, 1,-1, 1/
93
DATA (NHEL(I, 15),I=1,4) / 1, 1, 1,-1/
94
DATA (NHEL(I, 16),I=1,4) / 1, 1, 1, 1/
99
NTRY(IMIRROR)=NTRY(IMIRROR)+1
104
IF (MULTI_CHANNEL) THEN
114
WRITE(HEL_BUFF,'(20I5)') (0,I=1,NEXTERNAL)
118
IF (ISHEL(IMIRROR) .EQ. 0 .OR. NTRY(IMIRROR) .LE. MAXTRIES) THEN
120
IF (GOODHEL(I,IMIRROR) .OR. NTRY(IMIRROR).LE.MAXTRIES) THEN
121
T=MATRIX5(P ,NHEL(1,I),JC(1))
123
IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0
126
ELSE IF(POL(JJ).NE.1D0)THEN
127
T=T*(2D0-ABS(POL(JJ)))
135
IF(NTRY(IMIRROR).LE.MAXTRIES)THEN
137
IF (.NOT.GOODHEL(I,IMIRROR) .AND. (DABS(TS(I)).GT.ANS
138
$ *LIMHEL/NCOMB)) THEN
139
GOODHEL(I,IMIRROR)=.TRUE.
140
NGOOD(IMIRROR) = NGOOD(IMIRROR) +1
141
IGOOD(NGOOD(IMIRROR),IMIRROR) = I
142
PRINT *,'Added good helicity ',I,TS(I)*NCOMB/ANS
143
$ ,' in event ',NTRY(IMIRROR)
147
IF(NTRY(IMIRROR).EQ.MAXTRIES)THEN
148
ISHEL(IMIRROR)=MIN(ISUM_HEL,NGOOD(IMIRROR))
150
ELSE !LOOP OVER GOOD HELICITIES
151
DO J=1,ISHEL(IMIRROR)
152
JHEL(IMIRROR)=JHEL(IMIRROR)+1
153
IF (JHEL(IMIRROR) .GT. NGOOD(IMIRROR)) JHEL(IMIRROR)=1
154
HWGT = REAL(NGOOD(IMIRROR))/REAL(ISHEL(IMIRROR))
155
I = IGOOD(JHEL(IMIRROR),IMIRROR)
156
T=MATRIX5(P ,NHEL(1,I),JC(1))
158
IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0,POL(JJ)))
161
ELSE IF(POL(JJ).NE.1D0)THEN
162
T=T*(2D0-ABS(POL(JJ)))
168
IF (ISHEL(IMIRROR) .EQ. 1) THEN
169
WRITE(HEL_BUFF,'(20i5)')(NHEL(II,I),II=1,NEXTERNAL)
170
C Set right sign for ANS, based on sign of chosen helicity
174
IF (ISHEL(IMIRROR) .NE. 1) THEN
178
SUMHEL=SUMHEL+DABS(TS(I))/ANS
180
WRITE(HEL_BUFF,'(20i5)')(NHEL(II,I),II=1,NEXTERNAL)
181
C Set right sign for ANS, based on sign of chosen helicity
188
IF (MULTI_CHANNEL) THEN
193
IF (XTOT.NE.0D0) THEN
194
ANS=ANS*AMP2(SUBDIAG(5))/XTOT
203
REAL*8 FUNCTION MATRIX5(P,NHEL,IC)
205
C Generated by MadGraph 5 v. 2.0.0.beta3, 2013-02-14
206
C By the MadGraph Development Team
207
C Please visit us at https://launchpad.net/madgraph5
209
C Returns amplitude squared summed/avg over colors
210
C for the point with external lines W(0:6,NEXTERNAL)
212
C Process: u d~ > u d~ WEIGHTED=2
213
C Process: u s~ > u s~ WEIGHTED=2
214
C Process: s d~ > s d~ WEIGHTED=2
215
C Process: s u~ > s u~ WEIGHTED=2
222
PARAMETER (NGRAPHS=1)
224
INCLUDE 'nexternal.inc'
225
INCLUDE 'maxamps.inc'
226
INTEGER NWAVEFUNCS, NCOLOR
227
PARAMETER (NWAVEFUNCS=5, NCOLOR=2)
231
PARAMETER (IMAG1=(0D0,1D0))
235
REAL*8 P(0:3,NEXTERNAL)
236
INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
242
REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
243
COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
244
COMPLEX*16 W(6,NWAVEFUNCS)
245
C Needed for v4 models
247
DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
251
DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
252
COMMON/TO_AMPS/ AMP2, JAMP2
258
DATA (CF(I, 1),I= 1, 2) / 9, 3/
261
DATA (CF(I, 2),I= 1, 2) / 3, 9/
266
CALL IXXXXX(P(0,1),ZERO,NHEL(1),+1*IC(1),W(1,1))
267
CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
268
CALL OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
269
CALL IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
270
CALL FFV1_3(W(1,1),W(1,3),GC_11,ZERO,ZERO,W(1,5))
271
C Amplitude(s) for diagram number 1
272
CALL FFV1_0(W(1,4),W(1,2),W(1,5),GC_11,AMP(1))
273
JAMP(1)=+1D0/2D0*(+AMP(1))
274
JAMP(2)=+1D0/2D0*(-1D0/3D0*AMP(1))
279
ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
281
MATRIX5=MATRIX5+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
283
AMP2(1)=AMP2(1)+AMP(1)*DCONJG(AMP(1))
285
JAMP2(I)=JAMP2(I)+JAMP(I)*DCONJG(JAMP(I))