ISAJET/code/zjj4.F

   1 #include "isajet/pilot.h"
   2       REAL*8 FUNCTION ZJJ4(P1, P2, P3, P4, P5,NHEL,IM)
   3 C
   4 C Function generated by Madgraph + hand coding
   5 C Returns amplitude squared summed/ave over colors
   6 C for the point in phase space P1,P2,P3,P4,P5
   7 C and helicity NHEL(1)...
   8 C for the process: q(im) qb(im) -> z q(im) qb(im)
   9 C with Madgraph code IM
  10 C
  11 #if defined(CERNLIB_IMPNONE)
  12       IMPLICIT NONE
  13 #endif
  14 C
  15 C CONSTANTS
  16 C
  17       INTEGER    NGRAPHS,    NEIGEN,    NEXTERNAL
  18       PARAMETER (NGRAPHS=  8,NEIGEN=  2,NEXTERNAL=5)
  19       REAL*8     ZERO
  20       PARAMETER (ZERO=0D0)
  21 C
  22 C ARGUMENTS
  23 C
  24       REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
  25       INTEGER NHEL(NEXTERNAL),IM
  26 C
  27 C LOCAL VARIABLES
  28 C
  29       INTEGER I,J
  30       REAL*8 EIGEN_VAL(NEIGEN), EIGEN_VEC(NGRAPHS,NEIGEN)
  31       COMPLEX*16 ZTEMP
  32       COMPLEX*16 AMP(NGRAPHS)
  33       COMPLEX*16 W1(6)  , W2(6)  , W3(6)  , W4(6)  , W5(6)
  34       COMPLEX*16 W6(6)  , W7(6)  , W8(6)  , W9(6)  , W10(6)
  35       COMPLEX*16 W11(6) , W12(6) , W13(6) , W14(6) , W15(6)
  36       COMPLEX*16 W16(6)
  37       REAL*8 GZUD(2)
  38 C
  39 C GLOBAL VARIABLES
  40 C
  41 #include "isajet/mgcoms.inc"
  42 C
  43 C COLOR DATA
  44 C
  45       DATA EIGEN_VAL(1  )/       5.9259259259259234D-01 /
  46       DATA EIGEN_VEC(1  ,1  )/  -3.5355339059327379D-01 /
  47       DATA EIGEN_VEC(2  ,1  )/   3.5355339059327379D-01 /
  48       DATA EIGEN_VEC(3  ,1  )/  -3.5355339059327379D-01 /
  49       DATA EIGEN_VEC(4  ,1  )/  -3.5355339059327379D-01 /
  50       DATA EIGEN_VEC(5  ,1  )/   3.5355339059327379D-01 /
  51       DATA EIGEN_VEC(6  ,1  )/  -3.5355339059327379D-01 /
  52       DATA EIGEN_VEC(7  ,1  )/   3.5355339059327379D-01 /
  53       DATA EIGEN_VEC(8  ,1  )/   3.5355339059327379D-01 /
  54       DATA EIGEN_VAL(2  )/       1.1851851851851847D+00 /
  55       DATA EIGEN_VEC(1  ,2  )/  -3.5355339059327379D-01 /
  56       DATA EIGEN_VEC(2  ,2  )/  -3.5355339059327379D-01 /
  57       DATA EIGEN_VEC(3  ,2  )/  -3.5355339059327379D-01 /
  58       DATA EIGEN_VEC(4  ,2  )/  -3.5355339059327379D-01 /
  59       DATA EIGEN_VEC(5  ,2  )/  -3.5355339059327379D-01 /
  60       DATA EIGEN_VEC(6  ,2  )/  -3.5355339059327379D-01 /
  61       DATA EIGEN_VEC(7  ,2  )/  -3.5355339059327379D-01 /
  62       DATA EIGEN_VEC(8  ,2  )/  -3.5355339059327379D-01 /
  63 C ----------
  64 C BEGIN CODE
  65 C ----------
  66       IF(IM.EQ.4.OR.IM.EQ.8.OR.IM.EQ.12) THEN
  67         GZUD(1)=GZD(1)
  68         GZUD(2)=GZD(2)
  69       ELSE
  70         GZUD(1)=GZU(1)
  71         GZUD(2)=GZU(2)
  72       ENDIF
  73 C
  74       CALL IXXXXX(P1  ,FMASS(IM),NHEL(1  ), 1,W1  )
  75       CALL OXXXXX(P2  ,FMASS(IM),NHEL(2  ),-1,W2  )
  76       CALL VXXXXX(P3  ,ZMASS,NHEL(3  ), 1,W3  )
  77       CALL OXXXXX(P4  ,FMASS(IM),NHEL(4  ), 1,W4  )
  78       CALL IXXXXX(P5  ,FMASS(IM),NHEL(5  ),-1,W5  )
  79       CALL JIOXXX(W1  ,W4  ,GG,ZERO,ZERO,W6  )
  80       CALL FVIXXX(W5  ,W3  ,GZUD,FMASS(IM),FWIDTH(IM),W7  )
  81       CALL IOVXXX(W7  ,W2  ,W6  ,GG,AMP(1  ))
  82       CALL FVOXXX(W2  ,W3  ,GZUD,FMASS(IM),FWIDTH(IM),W8  )
  83       CALL JIOXXX(W1  ,W8  ,GG,ZERO,ZERO,W9  )
  84       CALL IOVXXX(W5  ,W4  ,W9  ,GG,AMP(2  ))
  85       CALL IOVXXX(W5  ,W8  ,W6  ,GG,AMP(3  ))
  86       CALL FVIXXX(W1  ,W3  ,GZUD,FMASS(IM),FWIDTH(IM),W10 )
  87       CALL JIOXXX(W10 ,W4  ,GG,ZERO,ZERO,W11 )
  88       CALL IOVXXX(W5  ,W2  ,W11 ,GG,AMP(4  ))
  89       CALL JIOXXX(W10 ,W2  ,GG,ZERO,ZERO,W12 )
  90       CALL IOVXXX(W5  ,W4  ,W12 ,GG,AMP(5  ))
  91       CALL FVOXXX(W4  ,W3  ,GZUD,FMASS(IM),FWIDTH(IM),W13 )
  92       CALL JIOXXX(W5  ,W2  ,GG,ZERO,ZERO,W14 )
  93       CALL IOVXXX(W1  ,W13 ,W14 ,GG,AMP(6  ))
  94       CALL JIOXXX(W1  ,W2  ,GG,ZERO,ZERO,W15 )
  95       CALL FVOXXX(W4  ,W15 ,GG,FMASS(IM),FWIDTH(IM),W16 )
  96       CALL IOVXXX(W5  ,W16 ,W3  ,GZUD,AMP(7  ))
  97       CALL IOVXXX(W5  ,W13 ,W15 ,GG,AMP(8  ))
  98       ZJJ4 = 0.D0
  99       DO I = 1, NEIGEN
 100           ZTEMP = (0.D0,0.D0)
 101           DO J = 1, NGRAPHS
 102               ZTEMP = ZTEMP + EIGEN_VEC(J,I)*AMP(J)
 103           ENDDO
 104           ZJJ4 =ZJJ4+ZTEMP*EIGEN_VAL(I)*CONJG(ZTEMP)
 105       ENDDO
 106 C      CALL GAUGECHECK(AMP,ZTEMP,EIGEN_VEC,EIGEN_VAL,NGRAPHS,NEIGEN)
 107       END