1
SUBROUTINE SMATRIX1(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: g g > g g WEIGHTED=2
20
INCLUDE 'maxconfigs.inc'
21
INCLUDE 'nexternal.inc'
30
PARAMETER (THEL=2*NCOMB)
34
REAL*8 P(0:3,NEXTERNAL),ANS
38
INTEGER NHEL(NEXTERNAL,NCOMB),NTRY(2)
41
REAL*8 R,SUMHEL,TS(NCOMB)
43
INTEGER JC(NEXTERNAL),II
44
LOGICAL GOODHEL(NCOMB,2)
45
REAL*8 HWGT, XTOT, XTRY, XREJ, XR, YFRAC(0:NCOMB)
46
INTEGER NGOOD(2), IGOOD(NCOMB,2)
47
INTEGER JHEL(2), J, JJ
51
DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
52
COMMON/TO_AMPS/ AMP2, JAMP2
54
CHARACTER*101 HEL_BUFF
55
COMMON/TO_HELICITY/ HEL_BUFF
58
COMMON/TO_MIRROR/ IMIRROR
61
COMMON/TO_POLARIZATION/ POL
65
COMMON/TO_MATRIX/ISUM_HEL, MULTI_CHANNEL
66
INTEGER MAPCONFIG(0:LMAXCONFIGS), ICONFIG
67
COMMON/TO_MCONFIGS/MAPCONFIG, ICONFIG
68
INTEGER SUBDIAG(MAXSPROC),IB(2)
69
COMMON/TO_SUB_DIAG/SUBDIAG,IB
74
SAVE YFRAC, IGOOD, JHEL
75
DATA GOODHEL/THEL*.FALSE./
76
DATA (NHEL(I, 1),I=1,4) /-1,-1,-1,-1/
77
DATA (NHEL(I, 2),I=1,4) /-1,-1,-1, 1/
78
DATA (NHEL(I, 3),I=1,4) /-1,-1, 1,-1/
79
DATA (NHEL(I, 4),I=1,4) /-1,-1, 1, 1/
80
DATA (NHEL(I, 5),I=1,4) /-1, 1,-1,-1/
81
DATA (NHEL(I, 6),I=1,4) /-1, 1,-1, 1/
82
DATA (NHEL(I, 7),I=1,4) /-1, 1, 1,-1/
83
DATA (NHEL(I, 8),I=1,4) /-1, 1, 1, 1/
84
DATA (NHEL(I, 9),I=1,4) / 1,-1,-1,-1/
85
DATA (NHEL(I, 10),I=1,4) / 1,-1,-1, 1/
86
DATA (NHEL(I, 11),I=1,4) / 1,-1, 1,-1/
87
DATA (NHEL(I, 12),I=1,4) / 1,-1, 1, 1/
88
DATA (NHEL(I, 13),I=1,4) / 1, 1,-1,-1/
89
DATA (NHEL(I, 14),I=1,4) / 1, 1,-1, 1/
90
DATA (NHEL(I, 15),I=1,4) / 1, 1, 1,-1/
91
DATA (NHEL(I, 16),I=1,4) / 1, 1, 1, 1/
96
NTRY(IMIRROR)=NTRY(IMIRROR)+1
101
IF (MULTI_CHANNEL) THEN
111
WRITE(HEL_BUFF,'(20I5)') (0,I=1,NEXTERNAL)
115
IF (ISHEL(IMIRROR) .EQ. 0 .OR. NTRY(IMIRROR) .LE. MAXTRIES) THEN
117
IF (GOODHEL(I,IMIRROR) .OR. NTRY(IMIRROR).LE.MAXTRIES) THEN
118
T=MATRIX1(P ,NHEL(1,I),JC(1))
120
IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0
123
ELSE IF(POL(JJ).NE.1D0)THEN
124
T=T*(2D0-ABS(POL(JJ)))
132
IF(NTRY(IMIRROR).LE.MAXTRIES)THEN
134
IF (.NOT.GOODHEL(I,IMIRROR) .AND. (DABS(TS(I)).GT.ANS
135
$ *LIMHEL/NCOMB)) THEN
136
GOODHEL(I,IMIRROR)=.TRUE.
137
NGOOD(IMIRROR) = NGOOD(IMIRROR) +1
138
IGOOD(NGOOD(IMIRROR),IMIRROR) = I
139
PRINT *,'Added good helicity ',I,TS(I)*NCOMB/ANS
140
$ ,' in event ',NTRY(IMIRROR)
144
IF(NTRY(IMIRROR).EQ.MAXTRIES)THEN
145
ISHEL(IMIRROR)=MIN(ISUM_HEL,NGOOD(IMIRROR))
147
ELSE !LOOP OVER GOOD HELICITIES
148
DO J=1,ISHEL(IMIRROR)
149
JHEL(IMIRROR)=JHEL(IMIRROR)+1
150
IF (JHEL(IMIRROR) .GT. NGOOD(IMIRROR)) JHEL(IMIRROR)=1
151
HWGT = REAL(NGOOD(IMIRROR))/REAL(ISHEL(IMIRROR))
152
I = IGOOD(JHEL(IMIRROR),IMIRROR)
153
T=MATRIX1(P ,NHEL(1,I),JC(1))
155
IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0,POL(JJ)))
158
ELSE IF(POL(JJ).NE.1D0)THEN
159
T=T*(2D0-ABS(POL(JJ)))
165
IF (ISHEL(IMIRROR) .EQ. 1) THEN
166
WRITE(HEL_BUFF,'(20i5)')(NHEL(II,I),II=1,NEXTERNAL)
167
C Set right sign for ANS, based on sign of chosen helicity
171
IF (ISHEL(IMIRROR) .NE. 1) THEN
175
SUMHEL=SUMHEL+DABS(TS(I))/ANS
177
WRITE(HEL_BUFF,'(20i5)')(NHEL(II,I),II=1,NEXTERNAL)
178
C Set right sign for ANS, based on sign of chosen helicity
185
IF (MULTI_CHANNEL) THEN
190
IF (XTOT.NE.0D0) THEN
191
ANS=ANS*AMP2(SUBDIAG(1))/XTOT
200
REAL*8 FUNCTION MATRIX1(P,NHEL,IC)
202
C Generated by MadGraph 5 v. 2.0.0.beta3, 2013-02-14
203
C By the MadGraph Development Team
204
C Please visit us at https://launchpad.net/madgraph5
206
C Returns amplitude squared summed/avg over colors
207
C for the point with external lines W(0:6,NEXTERNAL)
209
C Process: g g > g g WEIGHTED=2
216
PARAMETER (NGRAPHS=6)
218
INCLUDE 'nexternal.inc'
219
INCLUDE 'maxamps.inc'
220
INTEGER NWAVEFUNCS, NCOLOR
221
PARAMETER (NWAVEFUNCS=5, NCOLOR=6)
225
PARAMETER (IMAG1=(0D0,1D0))
229
REAL*8 P(0:3,NEXTERNAL)
230
INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
236
REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
237
COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
238
COMPLEX*16 W(6,NWAVEFUNCS)
239
C Needed for v4 models
241
DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
245
DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
246
COMMON/TO_AMPS/ AMP2, JAMP2
252
DATA (CF(I, 1),I= 1, 6) / 19, -2, -2, -2, -2, 4/
255
DATA (CF(I, 2),I= 1, 6) / -2, 19, -2, 4, -2, -2/
258
DATA (CF(I, 3),I= 1, 6) / -2, -2, 19, -2, 4, -2/
261
DATA (CF(I, 4),I= 1, 6) / -2, 4, -2, 19, -2, -2/
264
DATA (CF(I, 5),I= 1, 6) / -2, -2, 4, -2, 19, -2/
267
DATA (CF(I, 6),I= 1, 6) / 4, -2, -2, -2, -2, 19/
272
CALL VXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
273
CALL VXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
274
CALL VXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
275
CALL VXXXXX(P(0,4),ZERO,NHEL(4),+1*IC(4),W(1,4))
276
C Amplitude(s) for diagram number 1
277
CALL VVVV1_0(W(1,1),W(1,2),W(1,3),W(1,4),GC_12,AMP(1))
278
CALL VVVV3_0(W(1,1),W(1,2),W(1,3),W(1,4),GC_12,AMP(2))
279
CALL VVVV4_0(W(1,1),W(1,2),W(1,3),W(1,4),GC_12,AMP(3))
280
CALL VVV1_1(W(1,1),W(1,2),GC_10,ZERO,ZERO,W(1,5))
281
C Amplitude(s) for diagram number 2
282
CALL VVV1_0(W(1,3),W(1,4),W(1,5),GC_10,AMP(4))
283
CALL VVV1_1(W(1,1),W(1,3),GC_10,ZERO,ZERO,W(1,5))
284
C Amplitude(s) for diagram number 3
285
CALL VVV1_0(W(1,2),W(1,4),W(1,5),GC_10,AMP(5))
286
CALL VVV1_1(W(1,1),W(1,4),GC_10,ZERO,ZERO,W(1,5))
287
C Amplitude(s) for diagram number 4
288
CALL VVV1_0(W(1,2),W(1,3),W(1,5),GC_10,AMP(6))
289
JAMP(1)=+2D0*(+AMP(3)-AMP(1)-AMP(4)+AMP(6))
290
JAMP(2)=+2D0*(+AMP(1)+AMP(2)+AMP(4)+AMP(5))
291
JAMP(3)=+2D0*(-AMP(3)-AMP(2)-AMP(5)-AMP(6))
292
JAMP(4)=+2D0*(+AMP(1)+AMP(2)+AMP(4)+AMP(5))
293
JAMP(5)=+2D0*(-AMP(3)-AMP(2)-AMP(5)-AMP(6))
294
JAMP(6)=+2D0*(+AMP(3)-AMP(1)-AMP(4)+AMP(6))
299
ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
301
MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
303
AMP2(2)=AMP2(2)+AMP(4)*DCONJG(AMP(4))
304
AMP2(3)=AMP2(3)+AMP(5)*DCONJG(AMP(5))
305
AMP2(4)=AMP2(4)+AMP(6)*DCONJG(AMP(6))
307
JAMP2(I)=JAMP2(I)+JAMP(I)*DCONJG(JAMP(I))