← Back to branch summary

~maddevelopers/mg5amcnlo/WWW5_caching

« back to all changes in this revision

Viewing changes to users/mardelcourt/PROC_970226/PROC_970226/SubProcesses/P0_qq_qqqq/matrix3.f

  • Committer: John Doe
  • Date: 2013-03-25 20:27:02 UTC
  • Revision ID: john.doe@gmail.com-20130325202702-5sk3t1r8h33ca4p4
first clean version

Show diffs side-by-side

added added

removed removed

Lines of Context:
users/mardelcourt/PROC_970226/PROC_970226/SubProcesses/P0_qq_qqqq/matrix3.f
 
1
      SUBROUTINE SMATRIX3(P,ANS)
 
2
C     
 
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
 
6
C     
 
7
C     MadGraph for Madevent Version
 
8
C     
 
9
C     Returns amplitude squared summed/avg over colors
 
10
C     and helicities
 
11
C     for the point in phase space P(0:3,NEXTERNAL)
 
12
C     
 
13
C     Process: u~ u~ > u u~ u~ u~ WEIGHTED=4
 
14
C     Process: s~ s~ > s s~ s~ s~ WEIGHTED=4
 
15
C     
 
16
      IMPLICIT NONE
 
17
C     
 
18
C     CONSTANTS
 
19
C     
 
20
      INCLUDE 'genps.inc'
 
21
      INCLUDE 'maxconfigs.inc'
 
22
      INCLUDE 'nexternal.inc'
 
23
      INCLUDE 'maxamps.inc'
 
24
      INTEGER                 NCOMB
 
25
      PARAMETER (             NCOMB=64)
 
26
      INTEGER    NGRAPHS
 
27
      PARAMETER (NGRAPHS=42)
 
28
      INTEGER    NDIAGS
 
29
      PARAMETER (NDIAGS=42)
 
30
      INTEGER    THEL
 
31
      PARAMETER (THEL=2*NCOMB)
 
32
C     
 
33
C     ARGUMENTS 
 
34
C     
 
35
      REAL*8 P(0:3,NEXTERNAL),ANS
 
36
C     
 
37
C     LOCAL VARIABLES 
 
38
C     
 
39
      INTEGER NHEL(NEXTERNAL,NCOMB),NTRY(2)
 
40
      INTEGER ISHEL(2)
 
41
      REAL*8 T,MATRIX3
 
42
      REAL*8 R,SUMHEL,TS(NCOMB)
 
43
      INTEGER I,IDEN
 
44
      INTEGER JC(NEXTERNAL),II
 
45
      LOGICAL GOODHEL(NCOMB,2)
 
46
      REAL*8 HWGT, XTOT, XTRY, XREJ, XR, YFRAC(0:NCOMB)
 
47
      INTEGER NGOOD(2), IGOOD(NCOMB,2)
 
48
      INTEGER JHEL(2), J, JJ
 
49
C     
 
50
C     GLOBAL VARIABLES
 
51
C     
 
52
      DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
 
53
      COMMON/TO_AMPS/  AMP2,       JAMP2
 
54
 
 
55
      CHARACTER*101         HEL_BUFF
 
56
      COMMON/TO_HELICITY/  HEL_BUFF
 
57
 
 
58
      INTEGER IMIRROR
 
59
      COMMON/TO_MIRROR/ IMIRROR
 
60
 
 
61
      REAL*8 POL(2)
 
62
      COMMON/TO_POLARIZATION/ POL
 
63
 
 
64
      INTEGER          ISUM_HEL
 
65
      LOGICAL                    MULTI_CHANNEL
 
66
      COMMON/TO_MATRIX/ISUM_HEL, MULTI_CHANNEL
 
67
      INTEGER MAPCONFIG(0:LMAXCONFIGS), ICONFIG
 
68
      COMMON/TO_MCONFIGS/MAPCONFIG, ICONFIG
 
69
      INTEGER SUBDIAG(MAXSPROC),IB(2)
 
70
      COMMON/TO_SUB_DIAG/SUBDIAG,IB
 
71
      DATA XTRY, XREJ /0,0/
 
72
      DATA NTRY /0,0/
 
73
      DATA NGOOD /0,0/
 
74
      DATA ISHEL/0,0/
 
75
      SAVE YFRAC, IGOOD, JHEL
 
76
      DATA GOODHEL/THEL*.FALSE./
 
77
      DATA (NHEL(I,   1),I=1,6) /-1,-1,-1,-1,-1,-1/
 
78
      DATA (NHEL(I,   2),I=1,6) /-1,-1,-1,-1,-1, 1/
 
79
      DATA (NHEL(I,   3),I=1,6) /-1,-1,-1,-1, 1,-1/
 
80
      DATA (NHEL(I,   4),I=1,6) /-1,-1,-1,-1, 1, 1/
 
81
      DATA (NHEL(I,   5),I=1,6) /-1,-1,-1, 1,-1,-1/
 
82
      DATA (NHEL(I,   6),I=1,6) /-1,-1,-1, 1,-1, 1/
 
83
      DATA (NHEL(I,   7),I=1,6) /-1,-1,-1, 1, 1,-1/
 
84
      DATA (NHEL(I,   8),I=1,6) /-1,-1,-1, 1, 1, 1/
 
85
      DATA (NHEL(I,   9),I=1,6) /-1,-1, 1,-1,-1,-1/
 
86
      DATA (NHEL(I,  10),I=1,6) /-1,-1, 1,-1,-1, 1/
 
87
      DATA (NHEL(I,  11),I=1,6) /-1,-1, 1,-1, 1,-1/
 
88
      DATA (NHEL(I,  12),I=1,6) /-1,-1, 1,-1, 1, 1/
 
89
      DATA (NHEL(I,  13),I=1,6) /-1,-1, 1, 1,-1,-1/
 
90
      DATA (NHEL(I,  14),I=1,6) /-1,-1, 1, 1,-1, 1/
 
91
      DATA (NHEL(I,  15),I=1,6) /-1,-1, 1, 1, 1,-1/
 
92
      DATA (NHEL(I,  16),I=1,6) /-1,-1, 1, 1, 1, 1/
 
93
      DATA (NHEL(I,  17),I=1,6) /-1, 1,-1,-1,-1,-1/
 
94
      DATA (NHEL(I,  18),I=1,6) /-1, 1,-1,-1,-1, 1/
 
95
      DATA (NHEL(I,  19),I=1,6) /-1, 1,-1,-1, 1,-1/
 
96
      DATA (NHEL(I,  20),I=1,6) /-1, 1,-1,-1, 1, 1/
 
97
      DATA (NHEL(I,  21),I=1,6) /-1, 1,-1, 1,-1,-1/
 
98
      DATA (NHEL(I,  22),I=1,6) /-1, 1,-1, 1,-1, 1/
 
99
      DATA (NHEL(I,  23),I=1,6) /-1, 1,-1, 1, 1,-1/
 
100
      DATA (NHEL(I,  24),I=1,6) /-1, 1,-1, 1, 1, 1/
 
101
      DATA (NHEL(I,  25),I=1,6) /-1, 1, 1,-1,-1,-1/
 
102
      DATA (NHEL(I,  26),I=1,6) /-1, 1, 1,-1,-1, 1/
 
103
      DATA (NHEL(I,  27),I=1,6) /-1, 1, 1,-1, 1,-1/
 
104
      DATA (NHEL(I,  28),I=1,6) /-1, 1, 1,-1, 1, 1/
 
105
      DATA (NHEL(I,  29),I=1,6) /-1, 1, 1, 1,-1,-1/
 
106
      DATA (NHEL(I,  30),I=1,6) /-1, 1, 1, 1,-1, 1/
 
107
      DATA (NHEL(I,  31),I=1,6) /-1, 1, 1, 1, 1,-1/
 
108
      DATA (NHEL(I,  32),I=1,6) /-1, 1, 1, 1, 1, 1/
 
109
      DATA (NHEL(I,  33),I=1,6) / 1,-1,-1,-1,-1,-1/
 
110
      DATA (NHEL(I,  34),I=1,6) / 1,-1,-1,-1,-1, 1/
 
111
      DATA (NHEL(I,  35),I=1,6) / 1,-1,-1,-1, 1,-1/
 
112
      DATA (NHEL(I,  36),I=1,6) / 1,-1,-1,-1, 1, 1/
 
113
      DATA (NHEL(I,  37),I=1,6) / 1,-1,-1, 1,-1,-1/
 
114
      DATA (NHEL(I,  38),I=1,6) / 1,-1,-1, 1,-1, 1/
 
115
      DATA (NHEL(I,  39),I=1,6) / 1,-1,-1, 1, 1,-1/
 
116
      DATA (NHEL(I,  40),I=1,6) / 1,-1,-1, 1, 1, 1/
 
117
      DATA (NHEL(I,  41),I=1,6) / 1,-1, 1,-1,-1,-1/
 
118
      DATA (NHEL(I,  42),I=1,6) / 1,-1, 1,-1,-1, 1/
 
119
      DATA (NHEL(I,  43),I=1,6) / 1,-1, 1,-1, 1,-1/
 
120
      DATA (NHEL(I,  44),I=1,6) / 1,-1, 1,-1, 1, 1/
 
121
      DATA (NHEL(I,  45),I=1,6) / 1,-1, 1, 1,-1,-1/
 
122
      DATA (NHEL(I,  46),I=1,6) / 1,-1, 1, 1,-1, 1/
 
123
      DATA (NHEL(I,  47),I=1,6) / 1,-1, 1, 1, 1,-1/
 
124
      DATA (NHEL(I,  48),I=1,6) / 1,-1, 1, 1, 1, 1/
 
125
      DATA (NHEL(I,  49),I=1,6) / 1, 1,-1,-1,-1,-1/
 
126
      DATA (NHEL(I,  50),I=1,6) / 1, 1,-1,-1,-1, 1/
 
127
      DATA (NHEL(I,  51),I=1,6) / 1, 1,-1,-1, 1,-1/
 
128
      DATA (NHEL(I,  52),I=1,6) / 1, 1,-1,-1, 1, 1/
 
129
      DATA (NHEL(I,  53),I=1,6) / 1, 1,-1, 1,-1,-1/
 
130
      DATA (NHEL(I,  54),I=1,6) / 1, 1,-1, 1,-1, 1/
 
131
      DATA (NHEL(I,  55),I=1,6) / 1, 1,-1, 1, 1,-1/
 
132
      DATA (NHEL(I,  56),I=1,6) / 1, 1,-1, 1, 1, 1/
 
133
      DATA (NHEL(I,  57),I=1,6) / 1, 1, 1,-1,-1,-1/
 
134
      DATA (NHEL(I,  58),I=1,6) / 1, 1, 1,-1,-1, 1/
 
135
      DATA (NHEL(I,  59),I=1,6) / 1, 1, 1,-1, 1,-1/
 
136
      DATA (NHEL(I,  60),I=1,6) / 1, 1, 1,-1, 1, 1/
 
137
      DATA (NHEL(I,  61),I=1,6) / 1, 1, 1, 1,-1,-1/
 
138
      DATA (NHEL(I,  62),I=1,6) / 1, 1, 1, 1,-1, 1/
 
139
      DATA (NHEL(I,  63),I=1,6) / 1, 1, 1, 1, 1,-1/
 
140
      DATA (NHEL(I,  64),I=1,6) / 1, 1, 1, 1, 1, 1/
 
141
      DATA IDEN/216/
 
142
C     ----------
 
143
C     BEGIN CODE
 
144
C     ----------
 
145
      NTRY(IMIRROR)=NTRY(IMIRROR)+1
 
146
      DO I=1,NEXTERNAL
 
147
        JC(I) = +1
 
148
      ENDDO
 
149
 
 
150
      IF (MULTI_CHANNEL) THEN
 
151
        DO I=1,NDIAGS
 
152
          AMP2(I)=0D0
 
153
        ENDDO
 
154
        JAMP2(0)=6
 
155
        DO I=1,INT(JAMP2(0))
 
156
          JAMP2(I)=0D0
 
157
        ENDDO
 
158
      ENDIF
 
159
      ANS = 0D0
 
160
      WRITE(HEL_BUFF,'(20I5)') (0,I=1,NEXTERNAL)
 
161
      DO I=1,NCOMB
 
162
        TS(I)=0D0
 
163
      ENDDO
 
164
      IF (ISHEL(IMIRROR) .EQ. 0 .OR. NTRY(IMIRROR) .LE. MAXTRIES) THEN
 
165
        DO I=1,NCOMB
 
166
          IF (GOODHEL(I,IMIRROR) .OR. NTRY(IMIRROR).LE.MAXTRIES) THEN
 
167
            T=MATRIX3(P ,NHEL(1,I),JC(1))
 
168
            DO JJ=1,NINCOMING
 
169
              IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0
 
170
     $         ,POL(JJ)))) THEN
 
171
                T=T*ABS(POL(JJ))
 
172
              ELSE IF(POL(JJ).NE.1D0)THEN
 
173
                T=T*(2D0-ABS(POL(JJ)))
 
174
              ENDIF
 
175
            ENDDO
 
176
            ANS=ANS+DABS(T)
 
177
            TS(I)=T
 
178
          ENDIF
 
179
        ENDDO
 
180
        JHEL(IMIRROR) = 1
 
181
        IF(NTRY(IMIRROR).LE.MAXTRIES)THEN
 
182
          DO I=1,NCOMB
 
183
            IF (.NOT.GOODHEL(I,IMIRROR) .AND. (DABS(TS(I)).GT.ANS
 
184
     $       *LIMHEL/NCOMB)) THEN
 
185
              GOODHEL(I,IMIRROR)=.TRUE.
 
186
              NGOOD(IMIRROR) = NGOOD(IMIRROR) +1
 
187
              IGOOD(NGOOD(IMIRROR),IMIRROR) = I
 
188
              PRINT *,'Added good helicity ',I,TS(I)*NCOMB/ANS
 
189
     $         ,' in event ',NTRY(IMIRROR)
 
190
            ENDIF
 
191
          ENDDO
 
192
        ENDIF
 
193
        IF(NTRY(IMIRROR).EQ.MAXTRIES)THEN
 
194
          ISHEL(IMIRROR)=MIN(ISUM_HEL,NGOOD(IMIRROR))
 
195
        ENDIF
 
196
      ELSE  !LOOP OVER GOOD HELICITIES
 
197
        DO J=1,ISHEL(IMIRROR)
 
198
          JHEL(IMIRROR)=JHEL(IMIRROR)+1
 
199
          IF (JHEL(IMIRROR) .GT. NGOOD(IMIRROR)) JHEL(IMIRROR)=1
 
200
          HWGT = REAL(NGOOD(IMIRROR))/REAL(ISHEL(IMIRROR))
 
201
          I = IGOOD(JHEL(IMIRROR),IMIRROR)
 
202
          T=MATRIX3(P ,NHEL(1,I),JC(1))
 
203
          DO JJ=1,NINCOMING
 
204
            IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0,POL(JJ)))
 
205
     $       ) THEN
 
206
              T=T*ABS(POL(JJ))
 
207
            ELSE IF(POL(JJ).NE.1D0)THEN
 
208
              T=T*(2D0-ABS(POL(JJ)))
 
209
            ENDIF
 
210
          ENDDO
 
211
          ANS=ANS+DABS(T)*HWGT
 
212
          TS(I)=T*HWGT
 
213
        ENDDO
 
214
        IF (ISHEL(IMIRROR) .EQ. 1) THEN
 
215
          WRITE(HEL_BUFF,'(20i5)')(NHEL(II,I),II=1,NEXTERNAL)
 
216
C         Set right sign for ANS, based on sign of chosen helicity
 
217
          ANS=DSIGN(ANS,TS(I))
 
218
        ENDIF
 
219
      ENDIF
 
220
      IF (ISHEL(IMIRROR) .NE. 1) THEN
 
221
        CALL RANMAR(R)
 
222
        SUMHEL=0D0
 
223
        DO I=1,NCOMB
 
224
          SUMHEL=SUMHEL+DABS(TS(I))/ANS
 
225
          IF(R.LT.SUMHEL)THEN
 
226
            WRITE(HEL_BUFF,'(20i5)')(NHEL(II,I),II=1,NEXTERNAL)
 
227
C           Set right sign for ANS, based on sign of chosen helicity
 
228
            ANS=DSIGN(ANS,TS(I))
 
229
            GOTO 10
 
230
          ENDIF
 
231
        ENDDO
 
232
 10     CONTINUE
 
233
      ENDIF
 
234
      IF (MULTI_CHANNEL) THEN
 
235
        XTOT=0D0
 
236
        DO I=1,NDIAGS
 
237
          XTOT=XTOT+AMP2(I)
 
238
        ENDDO
 
239
        IF (XTOT.NE.0D0) THEN
 
240
          ANS=ANS*AMP2(SUBDIAG(3))/XTOT
 
241
        ELSE
 
242
          ANS=0D0
 
243
        ENDIF
 
244
      ENDIF
 
245
      ANS=ANS/DBLE(IDEN)
 
246
      END
 
247
 
 
248
 
 
249
      REAL*8 FUNCTION MATRIX3(P,NHEL,IC)
 
250
C     
 
251
C     Generated by MadGraph 5 v. 2.0.0.beta3, 2013-02-14
 
252
C     By the MadGraph Development Team
 
253
C     Please visit us at https://launchpad.net/madgraph5
 
254
C     
 
255
C     Returns amplitude squared summed/avg over colors
 
256
C     for the point with external lines W(0:6,NEXTERNAL)
 
257
C     
 
258
C     Process: u~ u~ > u u~ u~ u~ WEIGHTED=4
 
259
C     Process: s~ s~ > s s~ s~ s~ WEIGHTED=4
 
260
C     
 
261
      IMPLICIT NONE
 
262
C     
 
263
C     CONSTANTS
 
264
C     
 
265
      INTEGER    NGRAPHS
 
266
      PARAMETER (NGRAPHS=42)
 
267
      INCLUDE 'genps.inc'
 
268
      INCLUDE 'nexternal.inc'
 
269
      INCLUDE 'maxamps.inc'
 
270
      INTEGER    NWAVEFUNCS,     NCOLOR
 
271
      PARAMETER (NWAVEFUNCS=16, NCOLOR=6)
 
272
      REAL*8     ZERO
 
273
      PARAMETER (ZERO=0D0)
 
274
      COMPLEX*16 IMAG1
 
275
      PARAMETER (IMAG1=(0D0,1D0))
 
276
C     
 
277
C     ARGUMENTS 
 
278
C     
 
279
      REAL*8 P(0:3,NEXTERNAL)
 
280
      INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
 
281
C     
 
282
C     LOCAL VARIABLES 
 
283
C     
 
284
      INTEGER I,J
 
285
      COMPLEX*16 ZTEMP
 
286
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
 
287
      COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
 
288
      COMPLEX*16 W(6,NWAVEFUNCS)
 
289
C     Needed for v4 models
 
290
      COMPLEX*16 DUM0,DUM1
 
291
      DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
 
292
C     
 
293
C     GLOBAL VARIABLES
 
294
C     
 
295
      DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXFLOW)
 
296
      COMMON/TO_AMPS/  AMP2,       JAMP2
 
297
      INCLUDE 'coupl.inc'
 
298
C     
 
299
C     COLOR DATA
 
300
C     
 
301
      DATA DENOM(1)/1/
 
302
      DATA (CF(I,  1),I=  1,  6) /   27,    9,    9,    3,    3,    9/
 
303
C     1 T(1,4) T(2,5) T(3,6)
 
304
      DATA DENOM(2)/1/
 
305
      DATA (CF(I,  2),I=  1,  6) /    9,   27,    3,    9,    9,    3/
 
306
C     1 T(1,4) T(2,6) T(3,5)
 
307
      DATA DENOM(3)/1/
 
308
      DATA (CF(I,  3),I=  1,  6) /    9,    3,   27,    9,    9,    3/
 
309
C     1 T(1,5) T(2,4) T(3,6)
 
310
      DATA DENOM(4)/1/
 
311
      DATA (CF(I,  4),I=  1,  6) /    3,    9,    9,   27,    3,    9/
 
312
C     1 T(1,5) T(2,6) T(3,4)
 
313
      DATA DENOM(5)/1/
 
314
      DATA (CF(I,  5),I=  1,  6) /    3,    9,    9,    3,   27,    9/
 
315
C     1 T(1,6) T(2,4) T(3,5)
 
316
      DATA DENOM(6)/1/
 
317
      DATA (CF(I,  6),I=  1,  6) /    9,    3,    3,    9,    9,   27/
 
318
C     1 T(1,6) T(2,5) T(3,4)
 
319
C     ----------
 
320
C     BEGIN CODE
 
321
C     ----------
 
322
      CALL OXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
 
323
      CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
 
324
      CALL OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
 
325
      CALL IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
 
326
      CALL IXXXXX(P(0,5),ZERO,NHEL(5),-1*IC(5),W(1,5))
 
327
      CALL IXXXXX(P(0,6),ZERO,NHEL(6),-1*IC(6),W(1,6))
 
328
      CALL FFV1_3(W(1,4),W(1,1),GC_11,ZERO,ZERO,W(1,7))
 
329
      CALL FFV1_3(W(1,5),W(1,2),GC_11,ZERO,ZERO,W(1,8))
 
330
      CALL FFV1_1(W(1,3),W(1,7),GC_11,ZERO,ZERO,W(1,9))
 
331
C     Amplitude(s) for diagram number 1
 
332
      CALL FFV1_0(W(1,6),W(1,9),W(1,8),GC_11,AMP(1))
 
333
      CALL FFV1_2(W(1,6),W(1,7),GC_11,ZERO,ZERO,W(1,10))
 
334
C     Amplitude(s) for diagram number 2
 
335
      CALL FFV1_0(W(1,10),W(1,3),W(1,8),GC_11,AMP(2))
 
336
      CALL FFV1_3(W(1,6),W(1,3),GC_11,ZERO,ZERO,W(1,11))
 
337
C     Amplitude(s) for diagram number 3
 
338
      CALL VVV1_0(W(1,7),W(1,8),W(1,11),GC_10,AMP(3))
 
339
      CALL FFV1_3(W(1,5),W(1,3),GC_11,ZERO,ZERO,W(1,12))
 
340
      CALL FFV1_1(W(1,2),W(1,7),GC_11,ZERO,ZERO,W(1,13))
 
341
C     Amplitude(s) for diagram number 4
 
342
      CALL FFV1_0(W(1,6),W(1,13),W(1,12),GC_11,AMP(4))
 
343
C     Amplitude(s) for diagram number 5
 
344
      CALL FFV1_0(W(1,10),W(1,2),W(1,12),GC_11,AMP(5))
 
345
      CALL FFV1_3(W(1,6),W(1,2),GC_11,ZERO,ZERO,W(1,10))
 
346
C     Amplitude(s) for diagram number 6
 
347
      CALL VVV1_0(W(1,7),W(1,12),W(1,10),GC_10,AMP(6))
 
348
      CALL FFV1_2(W(1,5),W(1,7),GC_11,ZERO,ZERO,W(1,14))
 
349
C     Amplitude(s) for diagram number 7
 
350
      CALL FFV1_0(W(1,14),W(1,3),W(1,10),GC_11,AMP(7))
 
351
C     Amplitude(s) for diagram number 8
 
352
      CALL FFV1_0(W(1,5),W(1,9),W(1,10),GC_11,AMP(8))
 
353
C     Amplitude(s) for diagram number 9
 
354
      CALL FFV1_0(W(1,14),W(1,2),W(1,11),GC_11,AMP(9))
 
355
C     Amplitude(s) for diagram number 10
 
356
      CALL FFV1_0(W(1,5),W(1,13),W(1,11),GC_11,AMP(10))
 
357
      CALL FFV1_3(W(1,4),W(1,2),GC_11,ZERO,ZERO,W(1,13))
 
358
      CALL FFV1_3(W(1,5),W(1,1),GC_11,ZERO,ZERO,W(1,14))
 
359
      CALL FFV1_1(W(1,3),W(1,13),GC_11,ZERO,ZERO,W(1,9))
 
360
C     Amplitude(s) for diagram number 11
 
361
      CALL FFV1_0(W(1,6),W(1,9),W(1,14),GC_11,AMP(11))
 
362
      CALL FFV1_2(W(1,6),W(1,13),GC_11,ZERO,ZERO,W(1,7))
 
363
C     Amplitude(s) for diagram number 12
 
364
      CALL FFV1_0(W(1,7),W(1,3),W(1,14),GC_11,AMP(12))
 
365
C     Amplitude(s) for diagram number 13
 
366
      CALL VVV1_0(W(1,13),W(1,14),W(1,11),GC_10,AMP(13))
 
367
      CALL FFV1_1(W(1,1),W(1,13),GC_11,ZERO,ZERO,W(1,15))
 
368
C     Amplitude(s) for diagram number 14
 
369
      CALL FFV1_0(W(1,6),W(1,15),W(1,12),GC_11,AMP(14))
 
370
C     Amplitude(s) for diagram number 15
 
371
      CALL FFV1_0(W(1,7),W(1,1),W(1,12),GC_11,AMP(15))
 
372
      CALL FFV1_3(W(1,6),W(1,1),GC_11,ZERO,ZERO,W(1,7))
 
373
C     Amplitude(s) for diagram number 16
 
374
      CALL VVV1_0(W(1,13),W(1,12),W(1,7),GC_10,AMP(16))
 
375
      CALL FFV1_2(W(1,5),W(1,13),GC_11,ZERO,ZERO,W(1,16))
 
376
C     Amplitude(s) for diagram number 17
 
377
      CALL FFV1_0(W(1,16),W(1,3),W(1,7),GC_11,AMP(17))
 
378
C     Amplitude(s) for diagram number 18
 
379
      CALL FFV1_0(W(1,5),W(1,9),W(1,7),GC_11,AMP(18))
 
380
C     Amplitude(s) for diagram number 19
 
381
      CALL FFV1_0(W(1,16),W(1,1),W(1,11),GC_11,AMP(19))
 
382
C     Amplitude(s) for diagram number 20
 
383
      CALL FFV1_0(W(1,5),W(1,15),W(1,11),GC_11,AMP(20))
 
384
      CALL FFV1_3(W(1,4),W(1,3),GC_11,ZERO,ZERO,W(1,15))
 
385
      CALL FFV1_1(W(1,2),W(1,15),GC_11,ZERO,ZERO,W(1,16))
 
386
C     Amplitude(s) for diagram number 21
 
387
      CALL FFV1_0(W(1,6),W(1,16),W(1,14),GC_11,AMP(21))
 
388
      CALL FFV1_2(W(1,6),W(1,15),GC_11,ZERO,ZERO,W(1,9))
 
389
C     Amplitude(s) for diagram number 22
 
390
      CALL FFV1_0(W(1,9),W(1,2),W(1,14),GC_11,AMP(22))
 
391
C     Amplitude(s) for diagram number 23
 
392
      CALL VVV1_0(W(1,15),W(1,14),W(1,10),GC_10,AMP(23))
 
393
      CALL FFV1_1(W(1,1),W(1,15),GC_11,ZERO,ZERO,W(1,13))
 
394
C     Amplitude(s) for diagram number 24
 
395
      CALL FFV1_0(W(1,6),W(1,13),W(1,8),GC_11,AMP(24))
 
396
C     Amplitude(s) for diagram number 25
 
397
      CALL FFV1_0(W(1,9),W(1,1),W(1,8),GC_11,AMP(25))
 
398
C     Amplitude(s) for diagram number 26
 
399
      CALL VVV1_0(W(1,15),W(1,8),W(1,7),GC_10,AMP(26))
 
400
      CALL FFV1_2(W(1,5),W(1,15),GC_11,ZERO,ZERO,W(1,9))
 
401
C     Amplitude(s) for diagram number 27
 
402
      CALL FFV1_0(W(1,9),W(1,2),W(1,7),GC_11,AMP(27))
 
403
C     Amplitude(s) for diagram number 28
 
404
      CALL FFV1_0(W(1,5),W(1,16),W(1,7),GC_11,AMP(28))
 
405
C     Amplitude(s) for diagram number 29
 
406
      CALL FFV1_0(W(1,9),W(1,1),W(1,10),GC_11,AMP(29))
 
407
C     Amplitude(s) for diagram number 30
 
408
      CALL FFV1_0(W(1,5),W(1,13),W(1,10),GC_11,AMP(30))
 
409
      CALL FFV1_2(W(1,4),W(1,14),GC_11,ZERO,ZERO,W(1,13))
 
410
C     Amplitude(s) for diagram number 31
 
411
      CALL FFV1_0(W(1,13),W(1,3),W(1,10),GC_11,AMP(31))
 
412
      CALL FFV1_2(W(1,4),W(1,10),GC_11,ZERO,ZERO,W(1,5))
 
413
C     Amplitude(s) for diagram number 32
 
414
      CALL FFV1_0(W(1,5),W(1,3),W(1,14),GC_11,AMP(32))
 
415
C     Amplitude(s) for diagram number 33
 
416
      CALL FFV1_0(W(1,13),W(1,2),W(1,11),GC_11,AMP(33))
 
417
      CALL FFV1_2(W(1,4),W(1,11),GC_11,ZERO,ZERO,W(1,13))
 
418
C     Amplitude(s) for diagram number 34
 
419
      CALL FFV1_0(W(1,13),W(1,2),W(1,14),GC_11,AMP(34))
 
420
      CALL FFV1_2(W(1,4),W(1,8),GC_11,ZERO,ZERO,W(1,14))
 
421
C     Amplitude(s) for diagram number 35
 
422
      CALL FFV1_0(W(1,14),W(1,3),W(1,7),GC_11,AMP(35))
 
423
      CALL FFV1_2(W(1,4),W(1,7),GC_11,ZERO,ZERO,W(1,9))
 
424
C     Amplitude(s) for diagram number 36
 
425
      CALL FFV1_0(W(1,9),W(1,3),W(1,8),GC_11,AMP(36))
 
426
C     Amplitude(s) for diagram number 37
 
427
      CALL FFV1_0(W(1,14),W(1,1),W(1,11),GC_11,AMP(37))
 
428
C     Amplitude(s) for diagram number 38
 
429
      CALL FFV1_0(W(1,13),W(1,1),W(1,8),GC_11,AMP(38))
 
430
      CALL FFV1_2(W(1,4),W(1,12),GC_11,ZERO,ZERO,W(1,13))
 
431
C     Amplitude(s) for diagram number 39
 
432
      CALL FFV1_0(W(1,13),W(1,2),W(1,7),GC_11,AMP(39))
 
433
C     Amplitude(s) for diagram number 40
 
434
      CALL FFV1_0(W(1,9),W(1,2),W(1,12),GC_11,AMP(40))
 
435
C     Amplitude(s) for diagram number 41
 
436
      CALL FFV1_0(W(1,13),W(1,1),W(1,10),GC_11,AMP(41))
 
437
C     Amplitude(s) for diagram number 42
 
438
      CALL FFV1_0(W(1,5),W(1,1),W(1,12),GC_11,AMP(42))
 
439
      JAMP(1)=+1D0/4D0*(+1D0/9D0*AMP(1)+1D0/9D0*AMP(2)+1D0/3D0*AMP(4)
 
440
     $ +1D0/3D0*AMP(5)+1D0/3D0*AMP(7)+1D0/3D0*AMP(8)+1D0/9D0*AMP(9)
 
441
     $ +1D0/9D0*AMP(10)+AMP(14)-IMAG1*AMP(16)+AMP(17)+1D0/3D0*AMP(19)
 
442
     $ +1D0/3D0*AMP(20)+AMP(22)+IMAG1*AMP(23)+1D0/3D0*AMP(24)
 
443
     $ +1D0/3D0*AMP(25)+AMP(30)+AMP(31)+1D0/3D0*AMP(33)+1D0/3D0
 
444
     $ *AMP(34)+1D0/3D0*AMP(35)+1D0/3D0*AMP(36)+1D0/9D0*AMP(37)
 
445
     $ +1D0/9D0*AMP(38)+AMP(40))
 
446
      JAMP(2)=+1D0/4D0*(-1D0/3D0*AMP(1)-1D0/3D0*AMP(2)-1D0/9D0*AMP(4)
 
447
     $ -1D0/9D0*AMP(5)-1D0/9D0*AMP(7)-1D0/9D0*AMP(8)-1D0/3D0*AMP(9)
 
448
     $ -1D0/3D0*AMP(10)-AMP(12)-IMAG1*AMP(13)-1D0/3D0*AMP(14)
 
449
     $ -1D0/3D0*AMP(15)-AMP(20)-AMP(24)+IMAG1*AMP(26)-AMP(27)
 
450
     $ -1D0/3D0*AMP(29)-1D0/3D0*AMP(30)-1D0/3D0*AMP(31)-1D0/3D0
 
451
     $ *AMP(32)-AMP(33)-AMP(36)-1D0/3D0*AMP(39)-1D0/3D0*AMP(40)
 
452
     $ -1D0/9D0*AMP(41)-1D0/9D0*AMP(42))
 
453
      JAMP(3)=+1D0/4D0*(-AMP(4)+IMAG1*AMP(6)-AMP(7)-1D0/3D0*AMP(9)
 
454
     $ -1D0/3D0*AMP(10)-1D0/9D0*AMP(11)-1D0/9D0*AMP(12)-1D0/3D0
 
455
     $ *AMP(14)-1D0/3D0*AMP(15)-1D0/3D0*AMP(17)-1D0/3D0*AMP(18)
 
456
     $ -1D0/9D0*AMP(19)-1D0/9D0*AMP(20)-1D0/3D0*AMP(21)-1D0/3D0
 
457
     $ *AMP(22)-AMP(25)-IMAG1*AMP(26)-AMP(28)-1D0/3D0*AMP(31)
 
458
     $ -1D0/3D0*AMP(32)-1D0/9D0*AMP(33)-1D0/9D0*AMP(34)-AMP(35)
 
459
     $ -1D0/3D0*AMP(37)-1D0/3D0*AMP(38)-AMP(42))
 
460
      JAMP(4)=+1D0/4D0*(+AMP(1)-IMAG1*AMP(3)+1D0/3D0*AMP(7)+1D0/3D0
 
461
     $ *AMP(8)+AMP(9)+1D0/3D0*AMP(11)+1D0/3D0*AMP(12)+AMP(15)
 
462
     $ +IMAG1*AMP(16)+AMP(18)+1D0/9D0*AMP(21)+1D0/9D0*AMP(22)
 
463
     $ +1D0/3D0*AMP(24)+1D0/3D0*AMP(25)+1D0/3D0*AMP(27)+1D0/3D0
 
464
     $ *AMP(28)+1D0/9D0*AMP(29)+1D0/9D0*AMP(30)+1D0/9D0*AMP(31)
 
465
     $ +1D0/9D0*AMP(32)+1D0/3D0*AMP(33)+1D0/3D0*AMP(34)+AMP(38)
 
466
     $ +AMP(39)+1D0/3D0*AMP(41)+1D0/3D0*AMP(42))
 
467
      JAMP(5)=+1D0/4D0*(+AMP(2)+IMAG1*AMP(3)+1D0/3D0*AMP(4)+1D0/3D0
 
468
     $ *AMP(5)+AMP(10)+1D0/3D0*AMP(11)+1D0/3D0*AMP(12)+1D0/9D0*AMP(14)
 
469
     $ +1D0/9D0*AMP(15)+1D0/9D0*AMP(17)+1D0/9D0*AMP(18)+1D0/3D0
 
470
     $ *AMP(19)+1D0/3D0*AMP(20)+AMP(21)-IMAG1*AMP(23)+1D0/3D0*AMP(27)
 
471
     $ +1D0/3D0*AMP(28)+AMP(29)+AMP(32)+1D0/3D0*AMP(35)+1D0/3D0
 
472
     $ *AMP(36)+AMP(37)+1D0/9D0*AMP(39)+1D0/9D0*AMP(40)+1D0/3D0
 
473
     $ *AMP(41)+1D0/3D0*AMP(42))
 
474
      JAMP(6)=+1D0/4D0*(-1D0/3D0*AMP(1)-1D0/3D0*AMP(2)-AMP(5)
 
475
     $ -IMAG1*AMP(6)-AMP(8)-AMP(11)+IMAG1*AMP(13)-1D0/3D0*AMP(17)
 
476
     $ -1D0/3D0*AMP(18)-AMP(19)-1D0/3D0*AMP(21)-1D0/3D0*AMP(22)
 
477
     $ -1D0/9D0*AMP(24)-1D0/9D0*AMP(25)-1D0/9D0*AMP(27)-1D0/9D0
 
478
     $ *AMP(28)-1D0/3D0*AMP(29)-1D0/3D0*AMP(30)-AMP(34)-1D0/9D0
 
479
     $ *AMP(35)-1D0/9D0*AMP(36)-1D0/3D0*AMP(37)-1D0/3D0*AMP(38)
 
480
     $ -1D0/3D0*AMP(39)-1D0/3D0*AMP(40)-AMP(41))
 
481
      MATRIX3 = 0.D0
 
482
      DO I = 1, NCOLOR
 
483
        ZTEMP = (0.D0,0.D0)
 
484
        DO J = 1, NCOLOR
 
485
          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
 
486
        ENDDO
 
487
        MATRIX3=MATRIX3+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
 
488
      ENDDO
 
489
      AMP2(1)=AMP2(1)+AMP(1)*DCONJG(AMP(1))
 
490
      AMP2(2)=AMP2(2)+AMP(2)*DCONJG(AMP(2))
 
491
      AMP2(3)=AMP2(3)+AMP(3)*DCONJG(AMP(3))
 
492
      AMP2(10)=AMP2(10)+AMP(10)*DCONJG(AMP(10))
 
493
      AMP2(9)=AMP2(9)+AMP(9)*DCONJG(AMP(9))
 
494
      AMP2(12)=AMP2(12)+AMP(12)*DCONJG(AMP(12))
 
495
      AMP2(11)=AMP2(11)+AMP(11)*DCONJG(AMP(11))
 
496
      AMP2(13)=AMP2(13)+AMP(13)*DCONJG(AMP(13))
 
497
      AMP2(34)=AMP2(34)+AMP(34)*DCONJG(AMP(34))
 
498
      AMP2(33)=AMP2(33)+AMP(33)*DCONJG(AMP(33))
 
499
      AMP2(20)=AMP2(20)+AMP(20)*DCONJG(AMP(20))
 
500
      AMP2(19)=AMP2(19)+AMP(19)*DCONJG(AMP(19))
 
501
      AMP2(38)=AMP2(38)+AMP(38)*DCONJG(AMP(38))
 
502
      AMP2(37)=AMP2(37)+AMP(37)*DCONJG(AMP(37))
 
503
      AMP2(4)=AMP2(4)+AMP(4)*DCONJG(AMP(4))
 
504
      AMP2(5)=AMP2(5)+AMP(5)*DCONJG(AMP(5))
 
505
      AMP2(6)=AMP2(6)+AMP(6)*DCONJG(AMP(6))
 
506
      AMP2(7)=AMP2(7)+AMP(7)*DCONJG(AMP(7))
 
507
      AMP2(8)=AMP2(8)+AMP(8)*DCONJG(AMP(8))
 
508
      AMP2(42)=AMP2(42)+AMP(42)*DCONJG(AMP(42))
 
509
      AMP2(41)=AMP2(41)+AMP(41)*DCONJG(AMP(41))
 
510
      AMP2(14)=AMP2(14)+AMP(14)*DCONJG(AMP(14))
 
511
      AMP2(15)=AMP2(15)+AMP(15)*DCONJG(AMP(15))
 
512
      AMP2(16)=AMP2(16)+AMP(16)*DCONJG(AMP(16))
 
513
      AMP2(17)=AMP2(17)+AMP(17)*DCONJG(AMP(17))
 
514
      AMP2(18)=AMP2(18)+AMP(18)*DCONJG(AMP(18))
 
515
      AMP2(21)=AMP2(21)+AMP(21)*DCONJG(AMP(21))
 
516
      AMP2(22)=AMP2(22)+AMP(22)*DCONJG(AMP(22))
 
517
      AMP2(23)=AMP2(23)+AMP(23)*DCONJG(AMP(23))
 
518
      AMP2(24)=AMP2(24)+AMP(24)*DCONJG(AMP(24))
 
519
      AMP2(25)=AMP2(25)+AMP(25)*DCONJG(AMP(25))
 
520
      AMP2(26)=AMP2(26)+AMP(26)*DCONJG(AMP(26))
 
521
      AMP2(27)=AMP2(27)+AMP(27)*DCONJG(AMP(27))
 
522
      AMP2(28)=AMP2(28)+AMP(28)*DCONJG(AMP(28))
 
523
      AMP2(29)=AMP2(29)+AMP(29)*DCONJG(AMP(29))
 
524
      AMP2(30)=AMP2(30)+AMP(30)*DCONJG(AMP(30))
 
525
      AMP2(31)=AMP2(31)+AMP(31)*DCONJG(AMP(31))
 
526
      AMP2(32)=AMP2(32)+AMP(32)*DCONJG(AMP(32))
 
527
      AMP2(35)=AMP2(35)+AMP(35)*DCONJG(AMP(35))
 
528
      AMP2(36)=AMP2(36)+AMP(36)*DCONJG(AMP(36))
 
529
      AMP2(39)=AMP2(39)+AMP(39)*DCONJG(AMP(39))
 
530
      AMP2(40)=AMP2(40)+AMP(40)*DCONJG(AMP(40))
 
531
      DO I = 1, NCOLOR
 
532
        JAMP2(I)=JAMP2(I)+JAMP(I)*DCONJG(JAMP(I))
 
533
      ENDDO
 
534
 
 
535
      END
 
536