← 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_gggg/matrix1.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_gggg/matrix1.f
 
1
      SUBROUTINE SMATRIX1(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~ > g g g g WEIGHTED=4
 
14
C     Process: s s~ > g g g g 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=159)
 
28
      INTEGER    NDIAGS
 
29
      PARAMETER (NDIAGS=123)
 
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,MATRIX1
 
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/864/
 
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)=24
 
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=MATRIX1(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=MATRIX1(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(1))/XTOT
 
241
        ELSE
 
242
          ANS=0D0
 
243
        ENDIF
 
244
      ENDIF
 
245
      ANS=ANS/DBLE(IDEN)
 
246
      END
 
247
 
 
248
 
 
249
      REAL*8 FUNCTION MATRIX1(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~ > g g g g WEIGHTED=4
 
259
C     Process: s s~ > g g g g WEIGHTED=4
 
260
C     
 
261
      IMPLICIT NONE
 
262
C     
 
263
C     CONSTANTS
 
264
C     
 
265
      INTEGER    NGRAPHS
 
266
      PARAMETER (NGRAPHS=159)
 
267
      INCLUDE 'genps.inc'
 
268
      INCLUDE 'nexternal.inc'
 
269
      INCLUDE 'maxamps.inc'
 
270
      INTEGER    NWAVEFUNCS,     NCOLOR
 
271
      PARAMETER (NWAVEFUNCS=24, NCOLOR=24)
 
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)/54/
 
302
      DATA (CF(I,  1),I=  1,  6) /  512,  -64,  -64,    8,    8,   80/
 
303
      DATA (CF(I,  1),I=  7, 12) /  -64,    8,    8,   -1,   -1,  -10/
 
304
      DATA (CF(I,  1),I= 13, 18) /    8,   -1,   80,  -10,   71,   62/
 
305
      DATA (CF(I,  1),I= 19, 24) /   -1,  -10,  -10,   62,   62,  -28/
 
306
C     1 T(3,4,5,6,2,1)
 
307
      DATA DENOM(2)/54/
 
308
      DATA (CF(I,  2),I=  1,  6) /  -64,  512,    8,   80,  -64,    8/
 
309
      DATA (CF(I,  2),I=  7, 12) /    8,  -64,   -1,  -10,    8,   -1/
 
310
      DATA (CF(I,  2),I= 13, 18) /   -1,  -10,  -10,   62,   62,  -28/
 
311
      DATA (CF(I,  2),I= 19, 24) /    8,   -1,   80,  -10,   71,   62/
 
312
C     1 T(3,4,6,5,2,1)
 
313
      DATA DENOM(3)/54/
 
314
      DATA (CF(I,  3),I=  1,  6) /  -64,    8,  512,  -64,   80,    8/
 
315
      DATA (CF(I,  3),I=  7, 12) /    8,   -1,   80,  -10,   71,   62/
 
316
      DATA (CF(I,  3),I= 13, 18) /  -64,    8,    8,   -1,   -1,  -10/
 
317
      DATA (CF(I,  3),I= 19, 24) /  -10,   -1,   62,  -28,  -10,   62/
 
318
C     1 T(3,5,4,6,2,1)
 
319
      DATA DENOM(4)/54/
 
320
      DATA (CF(I,  4),I=  1,  6) /    8,   80,  -64,  512,    8,  -64/
 
321
      DATA (CF(I,  4),I=  7, 12) /   -1,  -10,  -10,   62,   62,  -28/
 
322
      DATA (CF(I,  4),I= 13, 18) /    8,  -64,   -1,  -10,    8,   -1/
 
323
      DATA (CF(I,  4),I= 19, 24) /   -1,    8,   71,   62,   80,  -10/
 
324
C     1 T(3,5,6,4,2,1)
 
325
      DATA DENOM(5)/54/
 
326
      DATA (CF(I,  5),I=  1,  6) /    8,  -64,   80,    8,  512,  -64/
 
327
      DATA (CF(I,  5),I=  7, 12) /   -1,    8,   71,   62,   80,  -10/
 
328
      DATA (CF(I,  5),I= 13, 18) /  -10,   -1,   62,  -28,  -10,   62/
 
329
      DATA (CF(I,  5),I= 19, 24) /  -64,    8,    8,   -1,   -1,  -10/
 
330
C     1 T(3,6,4,5,2,1)
 
331
      DATA DENOM(6)/54/
 
332
      DATA (CF(I,  6),I=  1,  6) /   80,    8,    8,  -64,  -64,  512/
 
333
      DATA (CF(I,  6),I=  7, 12) /  -10,   -1,   62,  -28,  -10,   62/
 
334
      DATA (CF(I,  6),I= 13, 18) /   -1,    8,   71,   62,   80,  -10/
 
335
      DATA (CF(I,  6),I= 19, 24) /    8,  -64,   -1,  -10,    8,   -1/
 
336
C     1 T(3,6,5,4,2,1)
 
337
      DATA DENOM(7)/54/
 
338
      DATA (CF(I,  7),I=  1,  6) /  -64,    8,    8,   -1,   -1,  -10/
 
339
      DATA (CF(I,  7),I=  7, 12) /  512,  -64,  -64,    8,    8,   80/
 
340
      DATA (CF(I,  7),I= 13, 18) /   80,  -10,    8,   -1,   62,   71/
 
341
      DATA (CF(I,  7),I= 19, 24) /  -10,   62,   -1,  -10,  -28,   62/
 
342
C     1 T(4,3,5,6,2,1)
 
343
      DATA DENOM(8)/54/
 
344
      DATA (CF(I,  8),I=  1,  6) /    8,  -64,   -1,  -10,    8,   -1/
 
345
      DATA (CF(I,  8),I=  7, 12) /  -64,  512,    8,   80,  -64,    8/
 
346
      DATA (CF(I,  8),I= 13, 18) /  -10,   62,   -1,  -10,  -28,   62/
 
347
      DATA (CF(I,  8),I= 19, 24) /   80,  -10,    8,   -1,   62,   71/
 
348
C     1 T(4,3,6,5,2,1)
 
349
      DATA DENOM(9)/54/
 
350
      DATA (CF(I,  9),I=  1,  6) /    8,   -1,   80,  -10,   71,   62/
 
351
      DATA (CF(I,  9),I=  7, 12) /  -64,    8,  512,  -64,   80,    8/
 
352
      DATA (CF(I,  9),I= 13, 18) /    8,   -1,  -64,    8,  -10,   -1/
 
353
      DATA (CF(I,  9),I= 19, 24) /   62,  -28,  -10,   -1,   62,  -10/
 
354
C     1 T(4,5,3,6,2,1)
 
355
      DATA DENOM(10)/54/
 
356
      DATA (CF(I, 10),I=  1,  6) /   -1,  -10,  -10,   62,   62,  -28/
 
357
      DATA (CF(I, 10),I=  7, 12) /    8,   80,  -64,  512,    8,  -64/
 
358
      DATA (CF(I, 10),I= 13, 18) /   -1,  -10,    8,  -64,   -1,    8/
 
359
      DATA (CF(I, 10),I= 19, 24) /   71,   62,   -1,    8,  -10,   80/
 
360
C     1 T(4,5,6,3,2,1)
 
361
      DATA DENOM(11)/54/
 
362
      DATA (CF(I, 11),I=  1,  6) /   -1,    8,   71,   62,   80,  -10/
 
363
      DATA (CF(I, 11),I=  7, 12) /    8,  -64,   80,    8,  512,  -64/
 
364
      DATA (CF(I, 11),I= 13, 18) /   62,  -28,  -10,   -1,   62,  -10/
 
365
      DATA (CF(I, 11),I= 19, 24) /    8,   -1,  -64,    8,  -10,   -1/
 
366
C     1 T(4,6,3,5,2,1)
 
367
      DATA DENOM(12)/54/
 
368
      DATA (CF(I, 12),I=  1,  6) /  -10,   -1,   62,  -28,  -10,   62/
 
369
      DATA (CF(I, 12),I=  7, 12) /   80,    8,    8,  -64,  -64,  512/
 
370
      DATA (CF(I, 12),I= 13, 18) /   71,   62,   -1,    8,  -10,   80/
 
371
      DATA (CF(I, 12),I= 19, 24) /   -1,  -10,    8,  -64,   -1,    8/
 
372
C     1 T(4,6,5,3,2,1)
 
373
      DATA DENOM(13)/54/
 
374
      DATA (CF(I, 13),I=  1,  6) /    8,   -1,  -64,    8,  -10,   -1/
 
375
      DATA (CF(I, 13),I=  7, 12) /   80,  -10,    8,   -1,   62,   71/
 
376
      DATA (CF(I, 13),I= 13, 18) /  512,  -64,  -64,    8,    8,   80/
 
377
      DATA (CF(I, 13),I= 19, 24) /   62,  -10,  -28,   62,   -1,  -10/
 
378
C     1 T(5,3,4,6,2,1)
 
379
      DATA DENOM(14)/54/
 
380
      DATA (CF(I, 14),I=  1,  6) /   -1,  -10,    8,  -64,   -1,    8/
 
381
      DATA (CF(I, 14),I=  7, 12) /  -10,   62,   -1,  -10,  -28,   62/
 
382
      DATA (CF(I, 14),I= 13, 18) /  -64,  512,    8,   80,  -64,    8/
 
383
      DATA (CF(I, 14),I= 19, 24) /  -10,   80,   62,   71,    8,   -1/
 
384
C     1 T(5,3,6,4,2,1)
 
385
      DATA DENOM(15)/54/
 
386
      DATA (CF(I, 15),I=  1,  6) /   80,  -10,    8,   -1,   62,   71/
 
387
      DATA (CF(I, 15),I=  7, 12) /    8,   -1,  -64,    8,  -10,   -1/
 
388
      DATA (CF(I, 15),I= 13, 18) /  -64,    8,  512,  -64,   80,    8/
 
389
      DATA (CF(I, 15),I= 19, 24) /  -28,   62,   62,  -10,  -10,   -1/
 
390
C     1 T(5,4,3,6,2,1)
 
391
      DATA DENOM(16)/54/
 
392
      DATA (CF(I, 16),I=  1,  6) /  -10,   62,   -1,  -10,  -28,   62/
 
393
      DATA (CF(I, 16),I=  7, 12) /   -1,  -10,    8,  -64,   -1,    8/
 
394
      DATA (CF(I, 16),I= 13, 18) /    8,   80,  -64,  512,    8,  -64/
 
395
      DATA (CF(I, 16),I= 19, 24) /   62,   71,  -10,   80,   -1,    8/
 
396
C     1 T(5,4,6,3,2,1)
 
397
      DATA DENOM(17)/54/
 
398
      DATA (CF(I, 17),I=  1,  6) /   71,   62,   -1,    8,  -10,   80/
 
399
      DATA (CF(I, 17),I=  7, 12) /   62,  -28,  -10,   -1,   62,  -10/
 
400
      DATA (CF(I, 17),I= 13, 18) /    8,  -64,   80,    8,  512,  -64/
 
401
      DATA (CF(I, 17),I= 19, 24) /   -1,    8,  -10,   -1,  -64,    8/
 
402
C     1 T(5,6,3,4,2,1)
 
403
      DATA DENOM(18)/54/
 
404
      DATA (CF(I, 18),I=  1,  6) /   62,  -28,  -10,   -1,   62,  -10/
 
405
      DATA (CF(I, 18),I=  7, 12) /   71,   62,   -1,    8,  -10,   80/
 
406
      DATA (CF(I, 18),I= 13, 18) /   80,    8,    8,  -64,  -64,  512/
 
407
      DATA (CF(I, 18),I= 19, 24) /  -10,   -1,   -1,    8,    8,  -64/
 
408
C     1 T(5,6,4,3,2,1)
 
409
      DATA DENOM(19)/54/
 
410
      DATA (CF(I, 19),I=  1,  6) /   -1,    8,  -10,   -1,  -64,    8/
 
411
      DATA (CF(I, 19),I=  7, 12) /  -10,   80,   62,   71,    8,   -1/
 
412
      DATA (CF(I, 19),I= 13, 18) /   62,  -10,  -28,   62,   -1,  -10/
 
413
      DATA (CF(I, 19),I= 19, 24) /  512,  -64,  -64,    8,    8,   80/
 
414
C     1 T(6,3,4,5,2,1)
 
415
      DATA DENOM(20)/54/
 
416
      DATA (CF(I, 20),I=  1,  6) /  -10,   -1,   -1,    8,    8,  -64/
 
417
      DATA (CF(I, 20),I=  7, 12) /   62,  -10,  -28,   62,   -1,  -10/
 
418
      DATA (CF(I, 20),I= 13, 18) /  -10,   80,   62,   71,    8,   -1/
 
419
      DATA (CF(I, 20),I= 19, 24) /  -64,  512,    8,   80,  -64,    8/
 
420
C     1 T(6,3,5,4,2,1)
 
421
      DATA DENOM(21)/54/
 
422
      DATA (CF(I, 21),I=  1,  6) /  -10,   80,   62,   71,    8,   -1/
 
423
      DATA (CF(I, 21),I=  7, 12) /   -1,    8,  -10,   -1,  -64,    8/
 
424
      DATA (CF(I, 21),I= 13, 18) /  -28,   62,   62,  -10,  -10,   -1/
 
425
      DATA (CF(I, 21),I= 19, 24) /  -64,    8,  512,  -64,   80,    8/
 
426
C     1 T(6,4,3,5,2,1)
 
427
      DATA DENOM(22)/54/
 
428
      DATA (CF(I, 22),I=  1,  6) /   62,  -10,  -28,   62,   -1,  -10/
 
429
      DATA (CF(I, 22),I=  7, 12) /  -10,   -1,   -1,    8,    8,  -64/
 
430
      DATA (CF(I, 22),I= 13, 18) /   62,   71,  -10,   80,   -1,    8/
 
431
      DATA (CF(I, 22),I= 19, 24) /    8,   80,  -64,  512,    8,  -64/
 
432
C     1 T(6,4,5,3,2,1)
 
433
      DATA DENOM(23)/54/
 
434
      DATA (CF(I, 23),I=  1,  6) /   62,   71,  -10,   80,   -1,    8/
 
435
      DATA (CF(I, 23),I=  7, 12) /  -28,   62,   62,  -10,  -10,   -1/
 
436
      DATA (CF(I, 23),I= 13, 18) /   -1,    8,  -10,   -1,  -64,    8/
 
437
      DATA (CF(I, 23),I= 19, 24) /    8,  -64,   80,    8,  512,  -64/
 
438
C     1 T(6,5,3,4,2,1)
 
439
      DATA DENOM(24)/54/
 
440
      DATA (CF(I, 24),I=  1,  6) /  -28,   62,   62,  -10,  -10,   -1/
 
441
      DATA (CF(I, 24),I=  7, 12) /   62,   71,  -10,   80,   -1,    8/
 
442
      DATA (CF(I, 24),I= 13, 18) /  -10,   -1,   -1,    8,    8,  -64/
 
443
      DATA (CF(I, 24),I= 19, 24) /   80,    8,    8,  -64,  -64,  512/
 
444
C     1 T(6,5,4,3,2,1)
 
445
C     ----------
 
446
C     BEGIN CODE
 
447
C     ----------
 
448
      CALL IXXXXX(P(0,1),ZERO,NHEL(1),+1*IC(1),W(1,1))
 
449
      CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
 
450
      CALL VXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
 
451
      CALL VXXXXX(P(0,4),ZERO,NHEL(4),+1*IC(4),W(1,4))
 
452
      CALL VXXXXX(P(0,5),ZERO,NHEL(5),+1*IC(5),W(1,5))
 
453
      CALL VXXXXX(P(0,6),ZERO,NHEL(6),+1*IC(6),W(1,6))
 
454
      CALL VVV1_1(W(1,3),W(1,4),GC_10,ZERO,ZERO,W(1,7))
 
455
      CALL VVV1_1(W(1,5),W(1,6),GC_10,ZERO,ZERO,W(1,8))
 
456
      CALL FFV1_1(W(1,2),W(1,7),GC_11,ZERO,ZERO,W(1,9))
 
457
C     Amplitude(s) for diagram number 1
 
458
      CALL FFV1_0(W(1,1),W(1,9),W(1,8),GC_11,AMP(1))
 
459
      CALL FFV1_2(W(1,1),W(1,7),GC_11,ZERO,ZERO,W(1,10))
 
460
C     Amplitude(s) for diagram number 2
 
461
      CALL FFV1_0(W(1,10),W(1,2),W(1,8),GC_11,AMP(2))
 
462
      CALL FFV1_3(W(1,1),W(1,2),GC_11,ZERO,ZERO,W(1,11))
 
463
C     Amplitude(s) for diagram number 3
 
464
      CALL VVV1_0(W(1,7),W(1,8),W(1,11),GC_10,AMP(3))
 
465
      CALL FFV1_1(W(1,2),W(1,5),GC_11,ZERO,ZERO,W(1,12))
 
466
      CALL VVV1_1(W(1,7),W(1,6),GC_10,ZERO,ZERO,W(1,13))
 
467
C     Amplitude(s) for diagram number 4
 
468
      CALL FFV1_0(W(1,1),W(1,12),W(1,13),GC_11,AMP(4))
 
469
C     Amplitude(s) for diagram number 5
 
470
      CALL FFV1_0(W(1,10),W(1,12),W(1,6),GC_11,AMP(5))
 
471
      CALL FFV1_2(W(1,1),W(1,6),GC_11,ZERO,ZERO,W(1,14))
 
472
C     Amplitude(s) for diagram number 6
 
473
      CALL FFV1_0(W(1,14),W(1,12),W(1,7),GC_11,AMP(6))
 
474
      CALL FFV1_2(W(1,1),W(1,5),GC_11,ZERO,ZERO,W(1,15))
 
475
C     Amplitude(s) for diagram number 7
 
476
      CALL FFV1_0(W(1,15),W(1,2),W(1,13),GC_11,AMP(7))
 
477
C     Amplitude(s) for diagram number 8
 
478
      CALL FFV1_0(W(1,15),W(1,9),W(1,6),GC_11,AMP(8))
 
479
      CALL FFV1_1(W(1,2),W(1,6),GC_11,ZERO,ZERO,W(1,16))
 
480
C     Amplitude(s) for diagram number 9
 
481
      CALL FFV1_0(W(1,15),W(1,16),W(1,7),GC_11,AMP(9))
 
482
      CALL VVV1_1(W(1,7),W(1,5),GC_10,ZERO,ZERO,W(1,17))
 
483
C     Amplitude(s) for diagram number 10
 
484
      CALL FFV1_0(W(1,1),W(1,16),W(1,17),GC_11,AMP(10))
 
485
C     Amplitude(s) for diagram number 11
 
486
      CALL FFV1_0(W(1,10),W(1,16),W(1,5),GC_11,AMP(11))
 
487
C     Amplitude(s) for diagram number 12
 
488
      CALL FFV1_0(W(1,14),W(1,2),W(1,17),GC_11,AMP(12))
 
489
C     Amplitude(s) for diagram number 13
 
490
      CALL FFV1_0(W(1,14),W(1,9),W(1,5),GC_11,AMP(13))
 
491
C     Amplitude(s) for diagram number 14
 
492
      CALL VVVV1_0(W(1,7),W(1,5),W(1,6),W(1,11),GC_12,AMP(14))
 
493
      CALL VVVV3_0(W(1,7),W(1,5),W(1,6),W(1,11),GC_12,AMP(15))
 
494
      CALL VVVV4_0(W(1,7),W(1,5),W(1,6),W(1,11),GC_12,AMP(16))
 
495
C     Amplitude(s) for diagram number 15
 
496
      CALL VVV1_0(W(1,6),W(1,11),W(1,17),GC_10,AMP(17))
 
497
C     Amplitude(s) for diagram number 16
 
498
      CALL VVV1_0(W(1,5),W(1,11),W(1,13),GC_10,AMP(18))
 
499
      CALL VVV1_1(W(1,3),W(1,5),GC_10,ZERO,ZERO,W(1,13))
 
500
      CALL VVV1_1(W(1,4),W(1,6),GC_10,ZERO,ZERO,W(1,17))
 
501
      CALL FFV1_1(W(1,2),W(1,13),GC_11,ZERO,ZERO,W(1,7))
 
502
C     Amplitude(s) for diagram number 17
 
503
      CALL FFV1_0(W(1,1),W(1,7),W(1,17),GC_11,AMP(19))
 
504
      CALL FFV1_2(W(1,1),W(1,13),GC_11,ZERO,ZERO,W(1,9))
 
505
C     Amplitude(s) for diagram number 18
 
506
      CALL FFV1_0(W(1,9),W(1,2),W(1,17),GC_11,AMP(20))
 
507
C     Amplitude(s) for diagram number 19
 
508
      CALL VVV1_0(W(1,13),W(1,17),W(1,11),GC_10,AMP(21))
 
509
      CALL FFV1_1(W(1,2),W(1,4),GC_11,ZERO,ZERO,W(1,10))
 
510
      CALL VVV1_1(W(1,13),W(1,6),GC_10,ZERO,ZERO,W(1,18))
 
511
C     Amplitude(s) for diagram number 20
 
512
      CALL FFV1_0(W(1,1),W(1,10),W(1,18),GC_11,AMP(22))
 
513
C     Amplitude(s) for diagram number 21
 
514
      CALL FFV1_0(W(1,9),W(1,10),W(1,6),GC_11,AMP(23))
 
515
C     Amplitude(s) for diagram number 22
 
516
      CALL FFV1_0(W(1,14),W(1,10),W(1,13),GC_11,AMP(24))
 
517
      CALL FFV1_2(W(1,1),W(1,4),GC_11,ZERO,ZERO,W(1,19))
 
518
C     Amplitude(s) for diagram number 23
 
519
      CALL FFV1_0(W(1,19),W(1,2),W(1,18),GC_11,AMP(25))
 
520
C     Amplitude(s) for diagram number 24
 
521
      CALL FFV1_0(W(1,19),W(1,7),W(1,6),GC_11,AMP(26))
 
522
C     Amplitude(s) for diagram number 25
 
523
      CALL FFV1_0(W(1,19),W(1,16),W(1,13),GC_11,AMP(27))
 
524
      CALL VVV1_1(W(1,13),W(1,4),GC_10,ZERO,ZERO,W(1,20))
 
525
C     Amplitude(s) for diagram number 26
 
526
      CALL FFV1_0(W(1,1),W(1,16),W(1,20),GC_11,AMP(28))
 
527
C     Amplitude(s) for diagram number 27
 
528
      CALL FFV1_0(W(1,9),W(1,16),W(1,4),GC_11,AMP(29))
 
529
C     Amplitude(s) for diagram number 28
 
530
      CALL FFV1_0(W(1,14),W(1,2),W(1,20),GC_11,AMP(30))
 
531
C     Amplitude(s) for diagram number 29
 
532
      CALL FFV1_0(W(1,14),W(1,7),W(1,4),GC_11,AMP(31))
 
533
C     Amplitude(s) for diagram number 30
 
534
      CALL VVVV1_0(W(1,13),W(1,4),W(1,6),W(1,11),GC_12,AMP(32))
 
535
      CALL VVVV3_0(W(1,13),W(1,4),W(1,6),W(1,11),GC_12,AMP(33))
 
536
      CALL VVVV4_0(W(1,13),W(1,4),W(1,6),W(1,11),GC_12,AMP(34))
 
537
C     Amplitude(s) for diagram number 31
 
538
      CALL VVV1_0(W(1,6),W(1,11),W(1,20),GC_10,AMP(35))
 
539
C     Amplitude(s) for diagram number 32
 
540
      CALL VVV1_0(W(1,4),W(1,11),W(1,18),GC_10,AMP(36))
 
541
      CALL VVV1_1(W(1,3),W(1,6),GC_10,ZERO,ZERO,W(1,18))
 
542
      CALL VVV1_1(W(1,4),W(1,5),GC_10,ZERO,ZERO,W(1,20))
 
543
      CALL FFV1_1(W(1,2),W(1,18),GC_11,ZERO,ZERO,W(1,13))
 
544
C     Amplitude(s) for diagram number 33
 
545
      CALL FFV1_0(W(1,1),W(1,13),W(1,20),GC_11,AMP(37))
 
546
      CALL FFV1_2(W(1,1),W(1,18),GC_11,ZERO,ZERO,W(1,7))
 
547
C     Amplitude(s) for diagram number 34
 
548
      CALL FFV1_0(W(1,7),W(1,2),W(1,20),GC_11,AMP(38))
 
549
C     Amplitude(s) for diagram number 35
 
550
      CALL VVV1_0(W(1,18),W(1,20),W(1,11),GC_10,AMP(39))
 
551
      CALL VVV1_1(W(1,18),W(1,5),GC_10,ZERO,ZERO,W(1,9))
 
552
C     Amplitude(s) for diagram number 36
 
553
      CALL FFV1_0(W(1,1),W(1,10),W(1,9),GC_11,AMP(40))
 
554
C     Amplitude(s) for diagram number 37
 
555
      CALL FFV1_0(W(1,7),W(1,10),W(1,5),GC_11,AMP(41))
 
556
C     Amplitude(s) for diagram number 38
 
557
      CALL FFV1_0(W(1,15),W(1,10),W(1,18),GC_11,AMP(42))
 
558
C     Amplitude(s) for diagram number 39
 
559
      CALL FFV1_0(W(1,19),W(1,2),W(1,9),GC_11,AMP(43))
 
560
C     Amplitude(s) for diagram number 40
 
561
      CALL FFV1_0(W(1,19),W(1,13),W(1,5),GC_11,AMP(44))
 
562
C     Amplitude(s) for diagram number 41
 
563
      CALL FFV1_0(W(1,19),W(1,12),W(1,18),GC_11,AMP(45))
 
564
      CALL VVV1_1(W(1,18),W(1,4),GC_10,ZERO,ZERO,W(1,21))
 
565
C     Amplitude(s) for diagram number 42
 
566
      CALL FFV1_0(W(1,1),W(1,12),W(1,21),GC_11,AMP(46))
 
567
C     Amplitude(s) for diagram number 43
 
568
      CALL FFV1_0(W(1,7),W(1,12),W(1,4),GC_11,AMP(47))
 
569
C     Amplitude(s) for diagram number 44
 
570
      CALL FFV1_0(W(1,15),W(1,2),W(1,21),GC_11,AMP(48))
 
571
C     Amplitude(s) for diagram number 45
 
572
      CALL FFV1_0(W(1,15),W(1,13),W(1,4),GC_11,AMP(49))
 
573
C     Amplitude(s) for diagram number 46
 
574
      CALL VVVV1_0(W(1,18),W(1,4),W(1,5),W(1,11),GC_12,AMP(50))
 
575
      CALL VVVV3_0(W(1,18),W(1,4),W(1,5),W(1,11),GC_12,AMP(51))
 
576
      CALL VVVV4_0(W(1,18),W(1,4),W(1,5),W(1,11),GC_12,AMP(52))
 
577
C     Amplitude(s) for diagram number 47
 
578
      CALL VVV1_0(W(1,5),W(1,11),W(1,21),GC_10,AMP(53))
 
579
C     Amplitude(s) for diagram number 48
 
580
      CALL VVV1_0(W(1,4),W(1,11),W(1,9),GC_10,AMP(54))
 
581
      CALL FFV1_1(W(1,2),W(1,3),GC_11,ZERO,ZERO,W(1,9))
 
582
      CALL FFV1_1(W(1,9),W(1,6),GC_11,ZERO,ZERO,W(1,21))
 
583
C     Amplitude(s) for diagram number 49
 
584
      CALL FFV1_0(W(1,1),W(1,21),W(1,20),GC_11,AMP(55))
 
585
      CALL FFV1_3(W(1,1),W(1,9),GC_11,ZERO,ZERO,W(1,18))
 
586
C     Amplitude(s) for diagram number 50
 
587
      CALL VVV1_0(W(1,20),W(1,6),W(1,18),GC_10,AMP(56))
 
588
C     Amplitude(s) for diagram number 51
 
589
      CALL FFV1_0(W(1,14),W(1,9),W(1,20),GC_11,AMP(57))
 
590
      CALL FFV1_1(W(1,9),W(1,5),GC_11,ZERO,ZERO,W(1,13))
 
591
C     Amplitude(s) for diagram number 52
 
592
      CALL FFV1_0(W(1,1),W(1,13),W(1,17),GC_11,AMP(58))
 
593
C     Amplitude(s) for diagram number 53
 
594
      CALL VVV1_0(W(1,17),W(1,5),W(1,18),GC_10,AMP(59))
 
595
C     Amplitude(s) for diagram number 54
 
596
      CALL FFV1_0(W(1,15),W(1,9),W(1,17),GC_11,AMP(60))
 
597
C     Amplitude(s) for diagram number 55
 
598
      CALL FFV1_0(W(1,19),W(1,13),W(1,6),GC_11,AMP(61))
 
599
C     Amplitude(s) for diagram number 56
 
600
      CALL FFV1_0(W(1,19),W(1,21),W(1,5),GC_11,AMP(62))
 
601
C     Amplitude(s) for diagram number 57
 
602
      CALL FFV1_0(W(1,19),W(1,9),W(1,8),GC_11,AMP(63))
 
603
      CALL FFV1_1(W(1,9),W(1,4),GC_11,ZERO,ZERO,W(1,7))
 
604
C     Amplitude(s) for diagram number 58
 
605
      CALL FFV1_0(W(1,1),W(1,7),W(1,8),GC_11,AMP(64))
 
606
C     Amplitude(s) for diagram number 59
 
607
      CALL VVV1_0(W(1,4),W(1,8),W(1,18),GC_10,AMP(65))
 
608
C     Amplitude(s) for diagram number 60
 
609
      CALL FFV1_0(W(1,15),W(1,7),W(1,6),GC_11,AMP(66))
 
610
C     Amplitude(s) for diagram number 61
 
611
      CALL FFV1_0(W(1,15),W(1,21),W(1,4),GC_11,AMP(67))
 
612
C     Amplitude(s) for diagram number 62
 
613
      CALL FFV1_0(W(1,14),W(1,7),W(1,5),GC_11,AMP(68))
 
614
C     Amplitude(s) for diagram number 63
 
615
      CALL FFV1_0(W(1,14),W(1,13),W(1,4),GC_11,AMP(69))
 
616
      CALL VVVV1_1(W(1,4),W(1,5),W(1,6),GC_12,ZERO,ZERO,W(1,13))
 
617
      CALL VVVV3_1(W(1,4),W(1,5),W(1,6),GC_12,ZERO,ZERO,W(1,7))
 
618
      CALL VVVV4_1(W(1,4),W(1,5),W(1,6),GC_12,ZERO,ZERO,W(1,21))
 
619
C     Amplitude(s) for diagram number 64
 
620
      CALL FFV1_0(W(1,1),W(1,9),W(1,13),GC_11,AMP(70))
 
621
      CALL FFV1_0(W(1,1),W(1,9),W(1,7),GC_11,AMP(71))
 
622
      CALL FFV1_0(W(1,1),W(1,9),W(1,21),GC_11,AMP(72))
 
623
      CALL FFV1_2(W(1,1),W(1,3),GC_11,ZERO,ZERO,W(1,9))
 
624
      CALL FFV1_2(W(1,9),W(1,6),GC_11,ZERO,ZERO,W(1,18))
 
625
C     Amplitude(s) for diagram number 65
 
626
      CALL FFV1_0(W(1,18),W(1,2),W(1,20),GC_11,AMP(73))
 
627
      CALL FFV1_3(W(1,9),W(1,2),GC_11,ZERO,ZERO,W(1,22))
 
628
C     Amplitude(s) for diagram number 66
 
629
      CALL VVV1_0(W(1,20),W(1,6),W(1,22),GC_10,AMP(74))
 
630
C     Amplitude(s) for diagram number 67
 
631
      CALL FFV1_0(W(1,9),W(1,16),W(1,20),GC_11,AMP(75))
 
632
      CALL FFV1_2(W(1,9),W(1,5),GC_11,ZERO,ZERO,W(1,23))
 
633
C     Amplitude(s) for diagram number 68
 
634
      CALL FFV1_0(W(1,23),W(1,2),W(1,17),GC_11,AMP(76))
 
635
C     Amplitude(s) for diagram number 69
 
636
      CALL VVV1_0(W(1,17),W(1,5),W(1,22),GC_10,AMP(77))
 
637
C     Amplitude(s) for diagram number 70
 
638
      CALL FFV1_0(W(1,9),W(1,12),W(1,17),GC_11,AMP(78))
 
639
C     Amplitude(s) for diagram number 71
 
640
      CALL FFV1_0(W(1,23),W(1,10),W(1,6),GC_11,AMP(79))
 
641
C     Amplitude(s) for diagram number 72
 
642
      CALL FFV1_0(W(1,18),W(1,10),W(1,5),GC_11,AMP(80))
 
643
C     Amplitude(s) for diagram number 73
 
644
      CALL FFV1_0(W(1,9),W(1,10),W(1,8),GC_11,AMP(81))
 
645
      CALL FFV1_2(W(1,9),W(1,4),GC_11,ZERO,ZERO,W(1,24))
 
646
C     Amplitude(s) for diagram number 74
 
647
      CALL FFV1_0(W(1,24),W(1,2),W(1,8),GC_11,AMP(82))
 
648
C     Amplitude(s) for diagram number 75
 
649
      CALL VVV1_0(W(1,4),W(1,8),W(1,22),GC_10,AMP(83))
 
650
C     Amplitude(s) for diagram number 76
 
651
      CALL FFV1_0(W(1,24),W(1,12),W(1,6),GC_11,AMP(84))
 
652
C     Amplitude(s) for diagram number 77
 
653
      CALL FFV1_0(W(1,18),W(1,12),W(1,4),GC_11,AMP(85))
 
654
C     Amplitude(s) for diagram number 78
 
655
      CALL FFV1_0(W(1,24),W(1,16),W(1,5),GC_11,AMP(86))
 
656
C     Amplitude(s) for diagram number 79
 
657
      CALL FFV1_0(W(1,23),W(1,16),W(1,4),GC_11,AMP(87))
 
658
C     Amplitude(s) for diagram number 80
 
659
      CALL FFV1_0(W(1,9),W(1,2),W(1,13),GC_11,AMP(88))
 
660
      CALL FFV1_0(W(1,9),W(1,2),W(1,7),GC_11,AMP(89))
 
661
      CALL FFV1_0(W(1,9),W(1,2),W(1,21),GC_11,AMP(90))
 
662
      CALL VVV1_1(W(1,3),W(1,20),GC_10,ZERO,ZERO,W(1,9))
 
663
C     Amplitude(s) for diagram number 81
 
664
      CALL FFV1_0(W(1,1),W(1,16),W(1,9),GC_11,AMP(91))
 
665
      CALL FFV1_1(W(1,16),W(1,3),GC_11,ZERO,ZERO,W(1,23))
 
666
C     Amplitude(s) for diagram number 82
 
667
      CALL FFV1_0(W(1,1),W(1,23),W(1,20),GC_11,AMP(92))
 
668
C     Amplitude(s) for diagram number 83
 
669
      CALL FFV1_0(W(1,14),W(1,2),W(1,9),GC_11,AMP(93))
 
670
      CALL FFV1_2(W(1,14),W(1,3),GC_11,ZERO,ZERO,W(1,24))
 
671
C     Amplitude(s) for diagram number 84
 
672
      CALL FFV1_0(W(1,24),W(1,2),W(1,20),GC_11,AMP(94))
 
673
C     Amplitude(s) for diagram number 85
 
674
      CALL VVVV1_0(W(1,3),W(1,20),W(1,6),W(1,11),GC_12,AMP(95))
 
675
      CALL VVVV3_0(W(1,3),W(1,20),W(1,6),W(1,11),GC_12,AMP(96))
 
676
      CALL VVVV4_0(W(1,3),W(1,20),W(1,6),W(1,11),GC_12,AMP(97))
 
677
C     Amplitude(s) for diagram number 86
 
678
      CALL VVV1_0(W(1,6),W(1,11),W(1,9),GC_10,AMP(98))
 
679
      CALL VVV1_1(W(1,3),W(1,11),GC_10,ZERO,ZERO,W(1,9))
 
680
C     Amplitude(s) for diagram number 87
 
681
      CALL VVV1_0(W(1,20),W(1,6),W(1,9),GC_10,AMP(99))
 
682
      CALL VVV1_1(W(1,3),W(1,17),GC_10,ZERO,ZERO,W(1,20))
 
683
C     Amplitude(s) for diagram number 88
 
684
      CALL FFV1_0(W(1,1),W(1,12),W(1,20),GC_11,AMP(100))
 
685
      CALL FFV1_1(W(1,12),W(1,3),GC_11,ZERO,ZERO,W(1,18))
 
686
C     Amplitude(s) for diagram number 89
 
687
      CALL FFV1_0(W(1,1),W(1,18),W(1,17),GC_11,AMP(101))
 
688
C     Amplitude(s) for diagram number 90
 
689
      CALL FFV1_0(W(1,15),W(1,2),W(1,20),GC_11,AMP(102))
 
690
      CALL FFV1_2(W(1,15),W(1,3),GC_11,ZERO,ZERO,W(1,22))
 
691
C     Amplitude(s) for diagram number 91
 
692
      CALL FFV1_0(W(1,22),W(1,2),W(1,17),GC_11,AMP(103))
 
693
C     Amplitude(s) for diagram number 92
 
694
      CALL VVVV1_0(W(1,3),W(1,17),W(1,5),W(1,11),GC_12,AMP(104))
 
695
      CALL VVVV3_0(W(1,3),W(1,17),W(1,5),W(1,11),GC_12,AMP(105))
 
696
      CALL VVVV4_0(W(1,3),W(1,17),W(1,5),W(1,11),GC_12,AMP(106))
 
697
C     Amplitude(s) for diagram number 93
 
698
      CALL VVV1_0(W(1,5),W(1,11),W(1,20),GC_10,AMP(107))
 
699
C     Amplitude(s) for diagram number 94
 
700
      CALL VVV1_0(W(1,17),W(1,5),W(1,9),GC_10,AMP(108))
 
701
      CALL FFV1_1(W(1,10),W(1,3),GC_11,ZERO,ZERO,W(1,17))
 
702
C     Amplitude(s) for diagram number 95
 
703
      CALL FFV1_0(W(1,1),W(1,17),W(1,8),GC_11,AMP(109))
 
704
      CALL VVV1_1(W(1,3),W(1,8),GC_10,ZERO,ZERO,W(1,20))
 
705
C     Amplitude(s) for diagram number 96
 
706
      CALL FFV1_0(W(1,1),W(1,10),W(1,20),GC_11,AMP(110))
 
707
C     Amplitude(s) for diagram number 97
 
708
      CALL FFV1_0(W(1,15),W(1,17),W(1,6),GC_11,AMP(111))
 
709
C     Amplitude(s) for diagram number 98
 
710
      CALL FFV1_0(W(1,22),W(1,10),W(1,6),GC_11,AMP(112))
 
711
C     Amplitude(s) for diagram number 99
 
712
      CALL FFV1_0(W(1,14),W(1,17),W(1,5),GC_11,AMP(113))
 
713
C     Amplitude(s) for diagram number 100
 
714
      CALL FFV1_0(W(1,24),W(1,10),W(1,5),GC_11,AMP(114))
 
715
      CALL FFV1_2(W(1,19),W(1,3),GC_11,ZERO,ZERO,W(1,17))
 
716
C     Amplitude(s) for diagram number 101
 
717
      CALL FFV1_0(W(1,17),W(1,2),W(1,8),GC_11,AMP(115))
 
718
C     Amplitude(s) for diagram number 102
 
719
      CALL FFV1_0(W(1,19),W(1,2),W(1,20),GC_11,AMP(116))
 
720
C     Amplitude(s) for diagram number 103
 
721
      CALL FFV1_0(W(1,17),W(1,12),W(1,6),GC_11,AMP(117))
 
722
C     Amplitude(s) for diagram number 104
 
723
      CALL FFV1_0(W(1,19),W(1,18),W(1,6),GC_11,AMP(118))
 
724
C     Amplitude(s) for diagram number 105
 
725
      CALL FFV1_0(W(1,17),W(1,16),W(1,5),GC_11,AMP(119))
 
726
C     Amplitude(s) for diagram number 106
 
727
      CALL FFV1_0(W(1,19),W(1,23),W(1,5),GC_11,AMP(120))
 
728
C     Amplitude(s) for diagram number 107
 
729
      CALL VVVV1_0(W(1,3),W(1,4),W(1,8),W(1,11),GC_12,AMP(121))
 
730
      CALL VVVV3_0(W(1,3),W(1,4),W(1,8),W(1,11),GC_12,AMP(122))
 
731
      CALL VVVV4_0(W(1,3),W(1,4),W(1,8),W(1,11),GC_12,AMP(123))
 
732
C     Amplitude(s) for diagram number 108
 
733
      CALL VVV1_0(W(1,4),W(1,11),W(1,20),GC_10,AMP(124))
 
734
C     Amplitude(s) for diagram number 109
 
735
      CALL VVV1_0(W(1,4),W(1,8),W(1,9),GC_10,AMP(125))
 
736
C     Amplitude(s) for diagram number 110
 
737
      CALL FFV1_0(W(1,14),W(1,18),W(1,4),GC_11,AMP(126))
 
738
C     Amplitude(s) for diagram number 111
 
739
      CALL FFV1_0(W(1,24),W(1,12),W(1,4),GC_11,AMP(127))
 
740
C     Amplitude(s) for diagram number 112
 
741
      CALL FFV1_0(W(1,22),W(1,16),W(1,4),GC_11,AMP(128))
 
742
C     Amplitude(s) for diagram number 113
 
743
      CALL FFV1_0(W(1,15),W(1,23),W(1,4),GC_11,AMP(129))
 
744
      CALL VVVV1_1(W(1,3),W(1,4),W(1,5),GC_12,ZERO,ZERO,W(1,23))
 
745
      CALL VVVV3_1(W(1,3),W(1,4),W(1,5),GC_12,ZERO,ZERO,W(1,22))
 
746
      CALL VVVV4_1(W(1,3),W(1,4),W(1,5),GC_12,ZERO,ZERO,W(1,24))
 
747
C     Amplitude(s) for diagram number 114
 
748
      CALL FFV1_0(W(1,1),W(1,16),W(1,23),GC_11,AMP(130))
 
749
      CALL FFV1_0(W(1,1),W(1,16),W(1,22),GC_11,AMP(131))
 
750
      CALL FFV1_0(W(1,1),W(1,16),W(1,24),GC_11,AMP(132))
 
751
C     Amplitude(s) for diagram number 115
 
752
      CALL FFV1_0(W(1,14),W(1,2),W(1,23),GC_11,AMP(133))
 
753
      CALL FFV1_0(W(1,14),W(1,2),W(1,22),GC_11,AMP(134))
 
754
      CALL FFV1_0(W(1,14),W(1,2),W(1,24),GC_11,AMP(135))
 
755
C     Amplitude(s) for diagram number 116
 
756
      CALL VVV1_0(W(1,23),W(1,6),W(1,11),GC_10,AMP(136))
 
757
      CALL VVV1_0(W(1,22),W(1,6),W(1,11),GC_10,AMP(137))
 
758
      CALL VVV1_0(W(1,24),W(1,6),W(1,11),GC_10,AMP(138))
 
759
      CALL VVVV1_1(W(1,3),W(1,4),W(1,6),GC_12,ZERO,ZERO,W(1,24))
 
760
      CALL VVVV3_1(W(1,3),W(1,4),W(1,6),GC_12,ZERO,ZERO,W(1,22))
 
761
      CALL VVVV4_1(W(1,3),W(1,4),W(1,6),GC_12,ZERO,ZERO,W(1,23))
 
762
C     Amplitude(s) for diagram number 117
 
763
      CALL FFV1_0(W(1,1),W(1,12),W(1,24),GC_11,AMP(139))
 
764
      CALL FFV1_0(W(1,1),W(1,12),W(1,22),GC_11,AMP(140))
 
765
      CALL FFV1_0(W(1,1),W(1,12),W(1,23),GC_11,AMP(141))
 
766
C     Amplitude(s) for diagram number 118
 
767
      CALL FFV1_0(W(1,15),W(1,2),W(1,24),GC_11,AMP(142))
 
768
      CALL FFV1_0(W(1,15),W(1,2),W(1,22),GC_11,AMP(143))
 
769
      CALL FFV1_0(W(1,15),W(1,2),W(1,23),GC_11,AMP(144))
 
770
C     Amplitude(s) for diagram number 119
 
771
      CALL VVV1_0(W(1,24),W(1,5),W(1,11),GC_10,AMP(145))
 
772
      CALL VVV1_0(W(1,22),W(1,5),W(1,11),GC_10,AMP(146))
 
773
      CALL VVV1_0(W(1,23),W(1,5),W(1,11),GC_10,AMP(147))
 
774
      CALL VVVV1_1(W(1,3),W(1,5),W(1,6),GC_12,ZERO,ZERO,W(1,23))
 
775
      CALL VVVV3_1(W(1,3),W(1,5),W(1,6),GC_12,ZERO,ZERO,W(1,22))
 
776
      CALL VVVV4_1(W(1,3),W(1,5),W(1,6),GC_12,ZERO,ZERO,W(1,24))
 
777
C     Amplitude(s) for diagram number 120
 
778
      CALL FFV1_0(W(1,1),W(1,10),W(1,23),GC_11,AMP(148))
 
779
      CALL FFV1_0(W(1,1),W(1,10),W(1,22),GC_11,AMP(149))
 
780
      CALL FFV1_0(W(1,1),W(1,10),W(1,24),GC_11,AMP(150))
 
781
C     Amplitude(s) for diagram number 121
 
782
      CALL FFV1_0(W(1,19),W(1,2),W(1,23),GC_11,AMP(151))
 
783
      CALL FFV1_0(W(1,19),W(1,2),W(1,22),GC_11,AMP(152))
 
784
      CALL FFV1_0(W(1,19),W(1,2),W(1,24),GC_11,AMP(153))
 
785
C     Amplitude(s) for diagram number 122
 
786
      CALL VVV1_0(W(1,23),W(1,4),W(1,11),GC_10,AMP(154))
 
787
      CALL VVV1_0(W(1,22),W(1,4),W(1,11),GC_10,AMP(155))
 
788
      CALL VVV1_0(W(1,24),W(1,4),W(1,11),GC_10,AMP(156))
 
789
C     Amplitude(s) for diagram number 123
 
790
      CALL VVV1_0(W(1,3),W(1,13),W(1,11),GC_10,AMP(157))
 
791
      CALL VVV1_0(W(1,3),W(1,7),W(1,11),GC_10,AMP(158))
 
792
      CALL VVV1_0(W(1,3),W(1,21),W(1,11),GC_10,AMP(159))
 
793
      JAMP(1)=-AMP(1)+IMAG1*AMP(3)-AMP(12)-IMAG1*AMP(13)+IMAG1*AMP(14)
 
794
     $ -IMAG1*AMP(16)+IMAG1*AMP(17)-AMP(56)-IMAG1*AMP(57)-IMAG1
 
795
     $ *AMP(64)-AMP(65)+AMP(68)+AMP(72)-AMP(70)-AMP(93)+IMAG1*AMP(95)
 
796
     $ -IMAG1*AMP(97)+IMAG1*AMP(98)-IMAG1*AMP(99)+IMAG1*AMP(121)
 
797
     $ -IMAG1*AMP(123)-IMAG1*AMP(125)+AMP(135)-AMP(133)+IMAG1*AMP(136)
 
798
     $ -IMAG1*AMP(138)+IMAG1*AMP(157)-IMAG1*AMP(159)
 
799
      JAMP(2)=+AMP(1)-IMAG1*AMP(3)-AMP(7)-IMAG1*AMP(8)+IMAG1*AMP(15)
 
800
     $ +IMAG1*AMP(16)+IMAG1*AMP(18)-AMP(59)-IMAG1*AMP(60)+IMAG1
 
801
     $ *AMP(64)+AMP(65)+AMP(66)+AMP(71)+AMP(70)-AMP(102)+IMAG1
 
802
     $ *AMP(104)-IMAG1*AMP(106)+IMAG1*AMP(107)-IMAG1*AMP(108)
 
803
     $ -IMAG1*AMP(121)+IMAG1*AMP(123)+IMAG1*AMP(125)+AMP(144)-AMP(142)
 
804
     $ +IMAG1*AMP(145)-IMAG1*AMP(147)-IMAG1*AMP(158)-IMAG1*AMP(157)
 
805
      JAMP(3)=-AMP(19)+IMAG1*AMP(21)-AMP(30)-IMAG1*AMP(31)+IMAG1
 
806
     $ *AMP(32)-IMAG1*AMP(34)+IMAG1*AMP(35)+AMP(56)+IMAG1*AMP(57)
 
807
     $ -IMAG1*AMP(58)+AMP(59)+AMP(69)-AMP(72)-AMP(71)+AMP(93)
 
808
     $ -IMAG1*AMP(95)+IMAG1*AMP(97)-IMAG1*AMP(98)+IMAG1*AMP(99)
 
809
     $ +IMAG1*AMP(105)+IMAG1*AMP(106)+IMAG1*AMP(108)+AMP(134)
 
810
     $ +AMP(133)-IMAG1*AMP(137)-IMAG1*AMP(136)+IMAG1*AMP(158)
 
811
     $ +IMAG1*AMP(159)
 
812
      JAMP(4)=+AMP(19)-IMAG1*AMP(21)-AMP(25)-IMAG1*AMP(26)+IMAG1
 
813
     $ *AMP(33)+IMAG1*AMP(34)+IMAG1*AMP(36)+IMAG1*AMP(58)-AMP(59)
 
814
     $ +AMP(61)-IMAG1*AMP(63)+AMP(65)+AMP(71)+AMP(70)-IMAG1*AMP(105)
 
815
     $ -IMAG1*AMP(106)-IMAG1*AMP(108)-AMP(116)+IMAG1*AMP(122)
 
816
     $ +IMAG1*AMP(123)+IMAG1*AMP(124)+IMAG1*AMP(125)+AMP(153)-AMP(151)
 
817
     $ +IMAG1*AMP(154)-IMAG1*AMP(156)-IMAG1*AMP(158)-IMAG1*AMP(157)
 
818
      JAMP(5)=-AMP(37)+IMAG1*AMP(39)-AMP(48)-IMAG1*AMP(49)+IMAG1
 
819
     $ *AMP(50)-IMAG1*AMP(52)+IMAG1*AMP(53)-IMAG1*AMP(55)+AMP(56)
 
820
     $ +AMP(59)+IMAG1*AMP(60)+AMP(67)-AMP(72)-AMP(71)+IMAG1*AMP(96)
 
821
     $ +IMAG1*AMP(97)+IMAG1*AMP(99)+AMP(102)-IMAG1*AMP(104)+IMAG1
 
822
     $ *AMP(106)-IMAG1*AMP(107)+IMAG1*AMP(108)+AMP(143)+AMP(142)
 
823
     $ -IMAG1*AMP(146)-IMAG1*AMP(145)+IMAG1*AMP(158)+IMAG1*AMP(159)
 
824
      JAMP(6)=+AMP(37)-IMAG1*AMP(39)-AMP(43)-IMAG1*AMP(44)+IMAG1
 
825
     $ *AMP(51)+IMAG1*AMP(52)+IMAG1*AMP(54)+IMAG1*AMP(55)-AMP(56)
 
826
     $ +AMP(62)+IMAG1*AMP(63)-AMP(65)+AMP(72)-AMP(70)-IMAG1*AMP(96)
 
827
     $ -IMAG1*AMP(97)-IMAG1*AMP(99)+AMP(116)-IMAG1*AMP(122)-IMAG1
 
828
     $ *AMP(123)-IMAG1*AMP(124)-IMAG1*AMP(125)+AMP(152)+AMP(151)
 
829
     $ -IMAG1*AMP(155)-IMAG1*AMP(154)+IMAG1*AMP(157)-IMAG1*AMP(159)
 
830
      JAMP(7)=+AMP(1)-IMAG1*AMP(3)+AMP(12)+IMAG1*AMP(13)-IMAG1*AMP(14)
 
831
     $ +IMAG1*AMP(16)-IMAG1*AMP(17)-AMP(22)-IMAG1*AMP(24)+AMP(30)
 
832
     $ -IMAG1*AMP(32)-IMAG1*AMP(33)-IMAG1*AMP(35)-IMAG1*AMP(36)
 
833
     $ -IMAG1*AMP(109)-AMP(110)+AMP(113)-IMAG1*AMP(121)-IMAG1*AMP(122)
 
834
     $ -IMAG1*AMP(124)-AMP(135)-AMP(134)+IMAG1*AMP(137)+IMAG1*AMP(138)
 
835
     $ +AMP(150)-AMP(148)-IMAG1*AMP(154)+IMAG1*AMP(156)
 
836
      JAMP(8)=-AMP(1)+IMAG1*AMP(3)+AMP(7)+IMAG1*AMP(8)-IMAG1*AMP(15)
 
837
     $ -IMAG1*AMP(16)-IMAG1*AMP(18)-AMP(40)-IMAG1*AMP(42)+AMP(48)
 
838
     $ -IMAG1*AMP(50)-IMAG1*AMP(51)-IMAG1*AMP(53)-IMAG1*AMP(54)
 
839
     $ +IMAG1*AMP(109)+AMP(110)+AMP(111)+IMAG1*AMP(121)+IMAG1*AMP(122)
 
840
     $ +IMAG1*AMP(124)-AMP(144)-AMP(143)+IMAG1*AMP(146)+IMAG1*AMP(147)
 
841
     $ +AMP(149)+AMP(148)+IMAG1*AMP(155)+IMAG1*AMP(154)
 
842
      JAMP(9)=+AMP(22)+IMAG1*AMP(24)-AMP(30)+IMAG1*AMP(32)+IMAG1
 
843
     $ *AMP(33)+IMAG1*AMP(35)+IMAG1*AMP(36)-AMP(38)-IMAG1*AMP(39)
 
844
     $ +AMP(40)-IMAG1*AMP(41)+IMAG1*AMP(51)+IMAG1*AMP(52)+IMAG1
 
845
     $ *AMP(54)+AMP(93)-IMAG1*AMP(94)-IMAG1*AMP(95)-IMAG1*AMP(96)
 
846
     $ -IMAG1*AMP(98)+AMP(114)+AMP(134)+AMP(133)-IMAG1*AMP(137)
 
847
     $ -IMAG1*AMP(136)-AMP(150)-AMP(149)-IMAG1*AMP(155)-IMAG1*AMP(156)
 
848
      JAMP(10)=+AMP(38)+IMAG1*AMP(39)-AMP(40)+IMAG1*AMP(41)-IMAG1
 
849
     $ *AMP(51)-IMAG1*AMP(52)-IMAG1*AMP(54)-IMAG1*AMP(73)-AMP(74)
 
850
     $ +AMP(80)-IMAG1*AMP(81)-AMP(83)+AMP(90)-AMP(88)+IMAG1*AMP(96)
 
851
     $ +IMAG1*AMP(97)+IMAG1*AMP(99)+AMP(110)+IMAG1*AMP(122)+IMAG1
 
852
     $ *AMP(123)+IMAG1*AMP(124)+IMAG1*AMP(125)+AMP(149)+AMP(148)
 
853
     $ +IMAG1*AMP(155)+IMAG1*AMP(154)-IMAG1*AMP(157)+IMAG1*AMP(159)
 
854
      JAMP(11)=-AMP(20)-IMAG1*AMP(21)+AMP(22)-IMAG1*AMP(23)+IMAG1
 
855
     $ *AMP(33)+IMAG1*AMP(34)+IMAG1*AMP(36)+AMP(40)+IMAG1*AMP(42)
 
856
     $ -AMP(48)+IMAG1*AMP(50)+IMAG1*AMP(51)+IMAG1*AMP(53)+IMAG1
 
857
     $ *AMP(54)+AMP(102)-IMAG1*AMP(103)-IMAG1*AMP(104)-IMAG1*AMP(105)
 
858
     $ -IMAG1*AMP(107)+AMP(112)+AMP(143)+AMP(142)-IMAG1*AMP(146)
 
859
     $ -IMAG1*AMP(145)-AMP(150)-AMP(149)-IMAG1*AMP(155)-IMAG1*AMP(156)
 
860
      JAMP(12)=+AMP(20)+IMAG1*AMP(21)-AMP(22)+IMAG1*AMP(23)-IMAG1
 
861
     $ *AMP(33)-IMAG1*AMP(34)-IMAG1*AMP(36)-IMAG1*AMP(76)-AMP(77)
 
862
     $ +AMP(79)+IMAG1*AMP(81)+AMP(83)+AMP(89)+AMP(88)+IMAG1*AMP(105)
 
863
     $ +IMAG1*AMP(106)+IMAG1*AMP(108)-AMP(110)-IMAG1*AMP(122)
 
864
     $ -IMAG1*AMP(123)-IMAG1*AMP(124)-IMAG1*AMP(125)+AMP(150)-AMP(148)
 
865
     $ -IMAG1*AMP(154)+IMAG1*AMP(156)+IMAG1*AMP(158)+IMAG1*AMP(157)
 
866
      JAMP(13)=-AMP(4)-IMAG1*AMP(6)+AMP(12)-IMAG1*AMP(14)-IMAG1
 
867
     $ *AMP(15)-IMAG1*AMP(17)-IMAG1*AMP(18)+AMP(19)-IMAG1*AMP(21)
 
868
     $ +AMP(30)+IMAG1*AMP(31)-IMAG1*AMP(32)+IMAG1*AMP(34)-IMAG1
 
869
     $ *AMP(35)-AMP(100)-IMAG1*AMP(101)-IMAG1*AMP(104)-IMAG1*AMP(105)
 
870
     $ -IMAG1*AMP(107)+AMP(126)-AMP(135)-AMP(134)+IMAG1*AMP(137)
 
871
     $ +IMAG1*AMP(138)+AMP(141)-AMP(139)-IMAG1*AMP(145)+IMAG1*AMP(147)
 
872
      JAMP(14)=-AMP(19)+IMAG1*AMP(21)+AMP(25)+IMAG1*AMP(26)-IMAG1
 
873
     $ *AMP(33)-IMAG1*AMP(34)-IMAG1*AMP(36)+AMP(43)-IMAG1*AMP(45)
 
874
     $ -AMP(46)-IMAG1*AMP(50)-IMAG1*AMP(51)-IMAG1*AMP(53)-IMAG1
 
875
     $ *AMP(54)+AMP(100)+IMAG1*AMP(101)+IMAG1*AMP(104)+IMAG1*AMP(105)
 
876
     $ +IMAG1*AMP(107)+AMP(118)+AMP(140)+AMP(139)+IMAG1*AMP(146)
 
877
     $ +IMAG1*AMP(145)-AMP(153)-AMP(152)+IMAG1*AMP(155)+IMAG1*AMP(156)
 
878
      JAMP(15)=+AMP(4)+IMAG1*AMP(6)-AMP(12)+IMAG1*AMP(14)+IMAG1
 
879
     $ *AMP(15)+IMAG1*AMP(17)+IMAG1*AMP(18)+AMP(38)+IMAG1*AMP(39)
 
880
     $ +AMP(46)-IMAG1*AMP(47)+IMAG1*AMP(50)-IMAG1*AMP(52)+IMAG1
 
881
     $ *AMP(53)-AMP(93)+IMAG1*AMP(94)+IMAG1*AMP(95)+IMAG1*AMP(96)
 
882
     $ +IMAG1*AMP(98)+AMP(127)+AMP(135)-AMP(133)+IMAG1*AMP(136)
 
883
     $ -IMAG1*AMP(138)-AMP(141)-AMP(140)-IMAG1*AMP(146)-IMAG1*AMP(147)
 
884
      JAMP(16)=-AMP(38)-IMAG1*AMP(39)-AMP(46)+IMAG1*AMP(47)-IMAG1
 
885
     $ *AMP(50)+IMAG1*AMP(52)-IMAG1*AMP(53)+IMAG1*AMP(73)+AMP(74)
 
886
     $ +AMP(77)-IMAG1*AMP(78)+AMP(85)-AMP(90)-AMP(89)-IMAG1*AMP(96)
 
887
     $ -IMAG1*AMP(97)-IMAG1*AMP(99)+AMP(100)+IMAG1*AMP(104)-IMAG1
 
888
     $ *AMP(106)+IMAG1*AMP(107)-IMAG1*AMP(108)+AMP(140)+AMP(139)
 
889
     $ +IMAG1*AMP(146)+IMAG1*AMP(145)-IMAG1*AMP(158)-IMAG1*AMP(159)
 
890
      JAMP(17)=-AMP(2)-IMAG1*AMP(3)+AMP(4)-IMAG1*AMP(5)+IMAG1*AMP(15)
 
891
     $ +IMAG1*AMP(16)+IMAG1*AMP(18)-AMP(43)+IMAG1*AMP(45)+AMP(46)
 
892
     $ +IMAG1*AMP(50)+IMAG1*AMP(51)+IMAG1*AMP(53)+IMAG1*AMP(54)
 
893
     $ -IMAG1*AMP(115)+AMP(116)+AMP(117)-IMAG1*AMP(121)-IMAG1*AMP(122)
 
894
     $ -IMAG1*AMP(124)-AMP(141)-AMP(140)-IMAG1*AMP(146)-IMAG1*AMP(147)
 
895
     $ +AMP(152)+AMP(151)-IMAG1*AMP(155)-IMAG1*AMP(154)
 
896
      JAMP(18)=+AMP(2)+IMAG1*AMP(3)-AMP(4)+IMAG1*AMP(5)-IMAG1*AMP(15)
 
897
     $ -IMAG1*AMP(16)-IMAG1*AMP(18)-AMP(77)+IMAG1*AMP(78)-IMAG1
 
898
     $ *AMP(82)+AMP(83)+AMP(84)+AMP(89)+AMP(88)-AMP(100)-IMAG1
 
899
     $ *AMP(104)+IMAG1*AMP(106)-IMAG1*AMP(107)+IMAG1*AMP(108)
 
900
     $ +IMAG1*AMP(121)-IMAG1*AMP(123)-IMAG1*AMP(125)+AMP(141)-AMP(139)
 
901
     $ -IMAG1*AMP(145)+IMAG1*AMP(147)+IMAG1*AMP(158)+IMAG1*AMP(157)
 
902
      JAMP(19)=+AMP(7)-IMAG1*AMP(9)-AMP(10)-IMAG1*AMP(14)-IMAG1
 
903
     $ *AMP(15)-IMAG1*AMP(17)-IMAG1*AMP(18)+AMP(37)-IMAG1*AMP(39)
 
904
     $ +AMP(48)+IMAG1*AMP(49)-IMAG1*AMP(50)+IMAG1*AMP(52)-IMAG1
 
905
     $ *AMP(53)-AMP(91)-IMAG1*AMP(92)-IMAG1*AMP(95)-IMAG1*AMP(96)
 
906
     $ -IMAG1*AMP(98)+AMP(129)+AMP(132)-AMP(130)-IMAG1*AMP(136)
 
907
     $ +IMAG1*AMP(138)-AMP(144)-AMP(143)+IMAG1*AMP(146)+IMAG1*AMP(147)
 
908
      JAMP(20)=+AMP(25)-IMAG1*AMP(27)-AMP(28)-IMAG1*AMP(32)-IMAG1
 
909
     $ *AMP(33)-IMAG1*AMP(35)-IMAG1*AMP(36)-AMP(37)+IMAG1*AMP(39)
 
910
     $ +AMP(43)+IMAG1*AMP(44)-IMAG1*AMP(51)-IMAG1*AMP(52)-IMAG1
 
911
     $ *AMP(54)+AMP(91)+IMAG1*AMP(92)+IMAG1*AMP(95)+IMAG1*AMP(96)
 
912
     $ +IMAG1*AMP(98)+AMP(120)+AMP(131)+AMP(130)+IMAG1*AMP(137)
 
913
     $ +IMAG1*AMP(136)-AMP(153)-AMP(152)+IMAG1*AMP(155)+IMAG1*AMP(156)
 
914
      JAMP(21)=-AMP(7)+IMAG1*AMP(9)+AMP(10)+IMAG1*AMP(14)+IMAG1
 
915
     $ *AMP(15)+IMAG1*AMP(17)+IMAG1*AMP(18)+AMP(20)+IMAG1*AMP(21)
 
916
     $ +AMP(28)-IMAG1*AMP(29)+IMAG1*AMP(32)-IMAG1*AMP(34)+IMAG1
 
917
     $ *AMP(35)-AMP(102)+IMAG1*AMP(103)+IMAG1*AMP(104)+IMAG1*AMP(105)
 
918
     $ +IMAG1*AMP(107)+AMP(128)-AMP(132)-AMP(131)-IMAG1*AMP(137)
 
919
     $ -IMAG1*AMP(138)+AMP(144)-AMP(142)+IMAG1*AMP(145)-IMAG1*AMP(147)
 
920
      JAMP(22)=-AMP(20)-IMAG1*AMP(21)-AMP(28)+IMAG1*AMP(29)-IMAG1
 
921
     $ *AMP(32)+IMAG1*AMP(34)-IMAG1*AMP(35)+AMP(74)-IMAG1*AMP(75)
 
922
     $ +IMAG1*AMP(76)+AMP(77)+AMP(87)-AMP(90)-AMP(89)+AMP(91)
 
923
     $ +IMAG1*AMP(95)-IMAG1*AMP(97)+IMAG1*AMP(98)-IMAG1*AMP(99)
 
924
     $ -IMAG1*AMP(105)-IMAG1*AMP(106)-IMAG1*AMP(108)+AMP(131)+AMP(130)
 
925
     $ +IMAG1*AMP(137)+IMAG1*AMP(136)-IMAG1*AMP(158)-IMAG1*AMP(159)
 
926
      JAMP(23)=+AMP(2)+IMAG1*AMP(3)+AMP(10)-IMAG1*AMP(11)+IMAG1
 
927
     $ *AMP(14)-IMAG1*AMP(16)+IMAG1*AMP(17)-AMP(25)+IMAG1*AMP(27)
 
928
     $ +AMP(28)+IMAG1*AMP(32)+IMAG1*AMP(33)+IMAG1*AMP(35)+IMAG1
 
929
     $ *AMP(36)+IMAG1*AMP(115)-AMP(116)+AMP(119)+IMAG1*AMP(121)
 
930
     $ +IMAG1*AMP(122)+IMAG1*AMP(124)-AMP(132)-AMP(131)-IMAG1*AMP(137)
 
931
     $ -IMAG1*AMP(138)+AMP(153)-AMP(151)+IMAG1*AMP(154)-IMAG1*AMP(156)
 
932
      JAMP(24)=-AMP(2)-IMAG1*AMP(3)-AMP(10)+IMAG1*AMP(11)-IMAG1
 
933
     $ *AMP(14)+IMAG1*AMP(16)-IMAG1*AMP(17)-AMP(74)+IMAG1*AMP(75)
 
934
     $ +IMAG1*AMP(82)-AMP(83)+AMP(86)+AMP(90)-AMP(88)-AMP(91)
 
935
     $ -IMAG1*AMP(95)+IMAG1*AMP(97)-IMAG1*AMP(98)+IMAG1*AMP(99)
 
936
     $ -IMAG1*AMP(121)+IMAG1*AMP(123)+IMAG1*AMP(125)+AMP(132)-AMP(130)
 
937
     $ -IMAG1*AMP(136)+IMAG1*AMP(138)-IMAG1*AMP(157)+IMAG1*AMP(159)
 
938
      MATRIX1 = 0.D0
 
939
      DO I = 1, NCOLOR
 
940
        ZTEMP = (0.D0,0.D0)
 
941
        DO J = 1, NCOLOR
 
942
          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
 
943
        ENDDO
 
944
        MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
 
945
      ENDDO
 
946
      AMP2(1)=AMP2(1)+AMP(1)*DCONJG(AMP(1))
 
947
      AMP2(2)=AMP2(2)+AMP(2)*DCONJG(AMP(2))
 
948
      AMP2(3)=AMP2(3)+AMP(3)*DCONJG(AMP(3))
 
949
      AMP2(4)=AMP2(4)+AMP(4)*DCONJG(AMP(4))
 
950
      AMP2(5)=AMP2(5)+AMP(5)*DCONJG(AMP(5))
 
951
      AMP2(6)=AMP2(6)+AMP(6)*DCONJG(AMP(6))
 
952
      AMP2(7)=AMP2(7)+AMP(7)*DCONJG(AMP(7))
 
953
      AMP2(8)=AMP2(8)+AMP(8)*DCONJG(AMP(8))
 
954
      AMP2(9)=AMP2(9)+AMP(9)*DCONJG(AMP(9))
 
955
      AMP2(10)=AMP2(10)+AMP(10)*DCONJG(AMP(10))
 
956
      AMP2(11)=AMP2(11)+AMP(11)*DCONJG(AMP(11))
 
957
      AMP2(12)=AMP2(12)+AMP(12)*DCONJG(AMP(12))
 
958
      AMP2(13)=AMP2(13)+AMP(13)*DCONJG(AMP(13))
 
959
      AMP2(15)=AMP2(15)+AMP(17)*DCONJG(AMP(17))
 
960
      AMP2(16)=AMP2(16)+AMP(18)*DCONJG(AMP(18))
 
961
      AMP2(17)=AMP2(17)+AMP(19)*DCONJG(AMP(19))
 
962
      AMP2(18)=AMP2(18)+AMP(20)*DCONJG(AMP(20))
 
963
      AMP2(19)=AMP2(19)+AMP(21)*DCONJG(AMP(21))
 
964
      AMP2(20)=AMP2(20)+AMP(22)*DCONJG(AMP(22))
 
965
      AMP2(21)=AMP2(21)+AMP(23)*DCONJG(AMP(23))
 
966
      AMP2(22)=AMP2(22)+AMP(24)*DCONJG(AMP(24))
 
967
      AMP2(23)=AMP2(23)+AMP(25)*DCONJG(AMP(25))
 
968
      AMP2(24)=AMP2(24)+AMP(26)*DCONJG(AMP(26))
 
969
      AMP2(25)=AMP2(25)+AMP(27)*DCONJG(AMP(27))
 
970
      AMP2(26)=AMP2(26)+AMP(28)*DCONJG(AMP(28))
 
971
      AMP2(27)=AMP2(27)+AMP(29)*DCONJG(AMP(29))
 
972
      AMP2(28)=AMP2(28)+AMP(30)*DCONJG(AMP(30))
 
973
      AMP2(29)=AMP2(29)+AMP(31)*DCONJG(AMP(31))
 
974
      AMP2(31)=AMP2(31)+AMP(35)*DCONJG(AMP(35))
 
975
      AMP2(32)=AMP2(32)+AMP(36)*DCONJG(AMP(36))
 
976
      AMP2(33)=AMP2(33)+AMP(37)*DCONJG(AMP(37))
 
977
      AMP2(34)=AMP2(34)+AMP(38)*DCONJG(AMP(38))
 
978
      AMP2(35)=AMP2(35)+AMP(39)*DCONJG(AMP(39))
 
979
      AMP2(36)=AMP2(36)+AMP(40)*DCONJG(AMP(40))
 
980
      AMP2(37)=AMP2(37)+AMP(41)*DCONJG(AMP(41))
 
981
      AMP2(38)=AMP2(38)+AMP(42)*DCONJG(AMP(42))
 
982
      AMP2(39)=AMP2(39)+AMP(43)*DCONJG(AMP(43))
 
983
      AMP2(40)=AMP2(40)+AMP(44)*DCONJG(AMP(44))
 
984
      AMP2(41)=AMP2(41)+AMP(45)*DCONJG(AMP(45))
 
985
      AMP2(42)=AMP2(42)+AMP(46)*DCONJG(AMP(46))
 
986
      AMP2(43)=AMP2(43)+AMP(47)*DCONJG(AMP(47))
 
987
      AMP2(44)=AMP2(44)+AMP(48)*DCONJG(AMP(48))
 
988
      AMP2(45)=AMP2(45)+AMP(49)*DCONJG(AMP(49))
 
989
      AMP2(47)=AMP2(47)+AMP(53)*DCONJG(AMP(53))
 
990
      AMP2(48)=AMP2(48)+AMP(54)*DCONJG(AMP(54))
 
991
      AMP2(49)=AMP2(49)+AMP(55)*DCONJG(AMP(55))
 
992
      AMP2(50)=AMP2(50)+AMP(56)*DCONJG(AMP(56))
 
993
      AMP2(51)=AMP2(51)+AMP(57)*DCONJG(AMP(57))
 
994
      AMP2(52)=AMP2(52)+AMP(58)*DCONJG(AMP(58))
 
995
      AMP2(53)=AMP2(53)+AMP(59)*DCONJG(AMP(59))
 
996
      AMP2(54)=AMP2(54)+AMP(60)*DCONJG(AMP(60))
 
997
      AMP2(55)=AMP2(55)+AMP(61)*DCONJG(AMP(61))
 
998
      AMP2(56)=AMP2(56)+AMP(62)*DCONJG(AMP(62))
 
999
      AMP2(57)=AMP2(57)+AMP(63)*DCONJG(AMP(63))
 
1000
      AMP2(58)=AMP2(58)+AMP(64)*DCONJG(AMP(64))
 
1001
      AMP2(59)=AMP2(59)+AMP(65)*DCONJG(AMP(65))
 
1002
      AMP2(60)=AMP2(60)+AMP(66)*DCONJG(AMP(66))
 
1003
      AMP2(61)=AMP2(61)+AMP(67)*DCONJG(AMP(67))
 
1004
      AMP2(62)=AMP2(62)+AMP(68)*DCONJG(AMP(68))
 
1005
      AMP2(63)=AMP2(63)+AMP(69)*DCONJG(AMP(69))
 
1006
      AMP2(65)=AMP2(65)+AMP(73)*DCONJG(AMP(73))
 
1007
      AMP2(66)=AMP2(66)+AMP(74)*DCONJG(AMP(74))
 
1008
      AMP2(67)=AMP2(67)+AMP(75)*DCONJG(AMP(75))
 
1009
      AMP2(68)=AMP2(68)+AMP(76)*DCONJG(AMP(76))
 
1010
      AMP2(69)=AMP2(69)+AMP(77)*DCONJG(AMP(77))
 
1011
      AMP2(70)=AMP2(70)+AMP(78)*DCONJG(AMP(78))
 
1012
      AMP2(71)=AMP2(71)+AMP(79)*DCONJG(AMP(79))
 
1013
      AMP2(72)=AMP2(72)+AMP(80)*DCONJG(AMP(80))
 
1014
      AMP2(73)=AMP2(73)+AMP(81)*DCONJG(AMP(81))
 
1015
      AMP2(74)=AMP2(74)+AMP(82)*DCONJG(AMP(82))
 
1016
      AMP2(75)=AMP2(75)+AMP(83)*DCONJG(AMP(83))
 
1017
      AMP2(76)=AMP2(76)+AMP(84)*DCONJG(AMP(84))
 
1018
      AMP2(77)=AMP2(77)+AMP(85)*DCONJG(AMP(85))
 
1019
      AMP2(78)=AMP2(78)+AMP(86)*DCONJG(AMP(86))
 
1020
      AMP2(79)=AMP2(79)+AMP(87)*DCONJG(AMP(87))
 
1021
      AMP2(81)=AMP2(81)+AMP(91)*DCONJG(AMP(91))
 
1022
      AMP2(82)=AMP2(82)+AMP(92)*DCONJG(AMP(92))
 
1023
      AMP2(83)=AMP2(83)+AMP(93)*DCONJG(AMP(93))
 
1024
      AMP2(84)=AMP2(84)+AMP(94)*DCONJG(AMP(94))
 
1025
      AMP2(86)=AMP2(86)+AMP(98)*DCONJG(AMP(98))
 
1026
      AMP2(87)=AMP2(87)+AMP(99)*DCONJG(AMP(99))
 
1027
      AMP2(88)=AMP2(88)+AMP(100)*DCONJG(AMP(100))
 
1028
      AMP2(89)=AMP2(89)+AMP(101)*DCONJG(AMP(101))
 
1029
      AMP2(90)=AMP2(90)+AMP(102)*DCONJG(AMP(102))
 
1030
      AMP2(91)=AMP2(91)+AMP(103)*DCONJG(AMP(103))
 
1031
      AMP2(93)=AMP2(93)+AMP(107)*DCONJG(AMP(107))
 
1032
      AMP2(94)=AMP2(94)+AMP(108)*DCONJG(AMP(108))
 
1033
      AMP2(95)=AMP2(95)+AMP(109)*DCONJG(AMP(109))
 
1034
      AMP2(96)=AMP2(96)+AMP(110)*DCONJG(AMP(110))
 
1035
      AMP2(97)=AMP2(97)+AMP(111)*DCONJG(AMP(111))
 
1036
      AMP2(98)=AMP2(98)+AMP(112)*DCONJG(AMP(112))
 
1037
      AMP2(99)=AMP2(99)+AMP(113)*DCONJG(AMP(113))
 
1038
      AMP2(100)=AMP2(100)+AMP(114)*DCONJG(AMP(114))
 
1039
      AMP2(101)=AMP2(101)+AMP(115)*DCONJG(AMP(115))
 
1040
      AMP2(102)=AMP2(102)+AMP(116)*DCONJG(AMP(116))
 
1041
      AMP2(103)=AMP2(103)+AMP(117)*DCONJG(AMP(117))
 
1042
      AMP2(104)=AMP2(104)+AMP(118)*DCONJG(AMP(118))
 
1043
      AMP2(105)=AMP2(105)+AMP(119)*DCONJG(AMP(119))
 
1044
      AMP2(106)=AMP2(106)+AMP(120)*DCONJG(AMP(120))
 
1045
      AMP2(108)=AMP2(108)+AMP(124)*DCONJG(AMP(124))
 
1046
      AMP2(109)=AMP2(109)+AMP(125)*DCONJG(AMP(125))
 
1047
      AMP2(110)=AMP2(110)+AMP(126)*DCONJG(AMP(126))
 
1048
      AMP2(111)=AMP2(111)+AMP(127)*DCONJG(AMP(127))
 
1049
      AMP2(112)=AMP2(112)+AMP(128)*DCONJG(AMP(128))
 
1050
      AMP2(113)=AMP2(113)+AMP(129)*DCONJG(AMP(129))
 
1051
      DO I = 1, NCOLOR
 
1052
        JAMP2(I)=JAMP2(I)+JAMP(I)*DCONJG(JAMP(I))
 
1053
      ENDDO
 
1054
 
 
1055
      END
 
1056