#include "isajet/pilot.h"
      REAL*8 FUNCTION ZJJ6(P1, P2, P3, P4, P5,NHEL,IM1)
C  
C Function generated by Madgraph + hand coding
C Returns amplitude squared summed/avg over colors
C for the point in phase space p1,p2,p3,p4,...
C and helicity NHEL(1),NHEL(2),....
C  
C FOR PROCESS : q(im1) q(im1)  -> z q(im1) q(im1)  
C  
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
C  
C CONSTANTS
C  
      INTEGER    NGRAPHS,    NEIGEN,    NEXTERNAL       
      PARAMETER (NGRAPHS=  8,NEIGEN=  2,NEXTERNAL=5)    
      REAL*8     ZERO
      PARAMETER (ZERO=0D0)
C  
C ARGUMENTS 
C  
      REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
      INTEGER NHEL(NEXTERNAL)
      INTEGER IM1
C  
C LOCAL VARIABLES 
C  
      INTEGER I,J
      REAL*8 EIGEN_VAL(NEIGEN), EIGEN_VEC(NGRAPHS,NEIGEN)
      COMPLEX*16 ZTEMP
      COMPLEX*16 AMP(NGRAPHS)
      COMPLEX*16 W1(6)  , W2(6)  , W3(6)  , W4(6)  , W5(6)        
      COMPLEX*16 W6(6)  , W7(6)  , W8(6)  , W9(6)  , W10(6)       
      COMPLEX*16 W11(6) , W12(6) , W13(6) , W14(6) , W15(6) 
      REAL*8 GZUD(2)
C  
C GLOBAL VARIABLES
C  
#include "isajet/mgcoms.inc"
C  
C COLOR DATA
C  
      DATA EIGEN_VAL(1  )/       2.9629629629629617D-01 /                  
      DATA EIGEN_VEC(1  ,1  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(2  ,1  )/   3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(3  ,1  )/   3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(4  ,1  )/   3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(5  ,1  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(6  ,1  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(7  ,1  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(8  ,1  )/   3.5355339059327379D-01 /                  
      DATA EIGEN_VAL(2  )/       5.9259259259259234D-01 /                  
      DATA EIGEN_VEC(1  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(2  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(3  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(4  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(5  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(6  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(7  ,2  )/  -3.5355339059327379D-01 /                  
      DATA EIGEN_VEC(8  ,2  )/  -3.5355339059327379D-01 /                  
C ----------
C BEGIN CODE
C ----------
      IF(IM1.EQ.3.OR.IM1.EQ.7.OR.IM1.EQ.11) THEN
        GZUD(1)=GZD(1)
        GZUD(2)=GZD(2)
      ELSE
        GZUD(1)=GZU(1)
        GZUD(2)=GZU(2)
      ENDIF
C
      CALL IXXXXX(P1  ,FMASS(IM1),NHEL(1  ), 1,W1  )                       
      CALL IXXXXX(P2  ,FMASS(IM1),NHEL(2  ), 1,W2  )                       
      CALL VXXXXX(P3  ,ZMASS,NHEL(3  ), 1,W3  )                            
      CALL OXXXXX(P4  ,FMASS(IM1),NHEL(4  ), 1,W4  )                       
      CALL OXXXXX(P5  ,FMASS(IM1),NHEL(5  ), 1,W5  )                       
      CALL FVIXXX(W2  ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W6  )               
      CALL JIOXXX(W6  ,W4  ,GG,ZERO,ZERO,W7  )                             
      CALL IOVXXX(W1  ,W5  ,W7  ,GG,AMP(1  ))                              
      CALL JIOXXX(W1  ,W4  ,GG,ZERO,ZERO,W8  )                             
      CALL FVOXXX(W5  ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W9  )               
      CALL IOVXXX(W2  ,W9  ,W8  ,GG,AMP(2  ))                              
      CALL IOVXXX(W6  ,W5  ,W8  ,GG,AMP(3  ))                              
      CALL FVIXXX(W1  ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W10 )               
      CALL JIOXXX(W10 ,W4  ,GG,ZERO,ZERO,W11 )                             
      CALL IOVXXX(W2  ,W5  ,W11 ,GG,AMP(4  ))                              
      CALL JIOXXX(W2  ,W4  ,GG,ZERO,ZERO,W12 )                             
      CALL IOVXXX(W1  ,W9  ,W12 ,GG,AMP(5  ))                              
      CALL IOVXXX(W10 ,W5  ,W12 ,GG,AMP(6  ))                              
      CALL FVOXXX(W4  ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W13 )               
      CALL JIOXXX(W1  ,W5  ,GG,ZERO,ZERO,W14 )                             
      CALL IOVXXX(W2  ,W13 ,W14 ,GG,AMP(7  ))                              
      CALL JIOXXX(W2  ,W5  ,GG,ZERO,ZERO,W15 )                             
      CALL IOVXXX(W1  ,W13 ,W15 ,GG,AMP(8  ))                              
      ZJJ6 = 0.D0 
      DO I = 1, NEIGEN
          ZTEMP = (0.D0,0.D0)
          DO J = 1, NGRAPHS
              ZTEMP = ZTEMP + EIGEN_VEC(J,I)*AMP(J)
          ENDDO
          ZJJ6 =ZJJ6+ZTEMP*EIGEN_VAL(I)*CONJG(ZTEMP) 
      ENDDO
C      CALL GAUGECHECK(AMP,ZTEMP,EIGEN_VEC,EIGEN_VAL,NGRAPHS,NEIGEN)
      END