ISAJET/code/zjj7.F

   1 #include "isajet/pilot.h"
   2       REAL*8 FUNCTION ZJJ7(P1, P2, P3, P4, P5,NHEL,IM1)
   3 C
   4 C FUNCTION GENERATED BY MADGRAPH
   5 C RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
   6 C FOR THE POINT IN PHASE SPACE P1,P2,P3,P4,...
   7 C AND HELICITY NHEL(1),NHEL(2),....
   8 C
   9 C FOR PROCESS : g u  -> z g u
  10 C
  11       IMPLICIT NONE
  12 C
  13 C CONSTANTS
  14 C
  15       INTEGER    NGRAPHS,    NEIGEN,    NEXTERNAL
  16       PARAMETER (NGRAPHS=  8,NEIGEN=  2,NEXTERNAL=5)
  17       REAL*8     ZERO
  18       PARAMETER (ZERO=0D0)
  19 C
  20 C ARGUMENTS
  21 C
  22       REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
  23       INTEGER NHEL(NEXTERNAL)
  24       INTEGER IM1
  25 C
  26 C LOCAL VARIABLES
  27 C
  28       INTEGER I,J
  29       REAL*8 EIGEN_VAL(NEIGEN), EIGEN_VEC(NGRAPHS,NEIGEN)
  30       COMPLEX*16 ZTEMP
  31       COMPLEX*16 AMP(NGRAPHS)
  32       COMPLEX*16 W1(6)  , W2(6)  , W3(6)  , W4(6)  , W5(6)
  33       COMPLEX*16 W6(6)  , W7(6)  , W8(6)  , W9(6)  , W10(6)
  34       COMPLEX*16 W11(6) , W12(6) , W13(6) , W14(6) , W15(6)
  35       REAL*8 GZUD(2)
  36 C
  37 C GLOBAL VARIABLES
  38 C
  39 #include "isajet/mgcoms.inc"
  40 C
  41 C COLOR DATA
  42 C
  43       DATA EIGEN_VAL(1  )/       5.8333333333333270D-01 /
  44       DATA EIGEN_VEC(1  ,1  )/   4.0824829046386285D-01 /
  45       DATA EIGEN_VEC(2  ,1  )/   0.0000000000000000D+00 /
  46       DATA EIGEN_VEC(3  ,1  )/   4.0824829046386318D-01 /
  47       DATA EIGEN_VEC(4  ,1  )/   0.0000000000000000D+00 /
  48       DATA EIGEN_VEC(5  ,1  )/   4.0824829046386285D-01 /
  49       DATA EIGEN_VEC(6  ,1  )/   4.0824829046386285D-01 /
  50       DATA EIGEN_VEC(7  ,1  )/   4.0824829046386318D-01 /
  51       DATA EIGEN_VEC(8  ,1  )/   4.0824829046386318D-01 /
  52       DATA EIGEN_VAL(2  )/       1.7499999999999991D+00 /
  53       DATA EIGEN_VEC(1  ,2  )/   2.6726124191242445D-01 /
  54       DATA EIGEN_VEC(2  ,2  )/  -5.3452248382484879D-01 /
  55       DATA EIGEN_VEC(3  ,2  )/  -2.6726124191242434D-01 /
  56       DATA EIGEN_VEC(4  ,2  )/  -5.3452248382484879D-01 /
  57       DATA EIGEN_VEC(5  ,2  )/   2.6726124191242445D-01 /
  58       DATA EIGEN_VEC(6  ,2  )/   2.6726124191242445D-01 /
  59       DATA EIGEN_VEC(7  ,2  )/  -2.6726124191242434D-01 /
  60       DATA EIGEN_VEC(8  ,2  )/  -2.6726124191242434D-01 /
  61 C ----------
  62 C BEGIN CODE
  63 C ----------
  64       IF(IM1.EQ.3.OR.IM1.EQ.7.OR.IM1.EQ.11) THEN
  65         GZUD(1)=GZD(1)
  66         GZUD(2)=GZD(2)
  67       ELSE
  68         GZUD(1)=GZU(1)
  69         GZUD(2)=GZU(2)
  70       ENDIF
  71 C
  72       CALL VXXXXX(P1  , ZERO,NHEL(1  ),-1,W1  )
  73       CALL IXXXXX(P2  ,FMASS(IM1),NHEL(2  ), 1,W2  )
  74       CALL VXXXXX(P3  ,ZMASS,NHEL(3  ), 1,W3  )
  75       CALL VXXXXX(P4  , ZERO,NHEL(4  ), 1,W4  )
  76       CALL OXXXXX(P5  ,FMASS(IM1),NHEL(5  ), 1,W5  )
  77       CALL FVIXXX(W2  ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W6  )
  78       CALL FVIXXX(W6  ,W4  ,GG,FMASS(IM1),FWIDTH(IM1),W7  )
  79       CALL IOVXXX(W7  ,W5  ,W1  ,GG,AMP(1  ))
  80       CALL JGGXXX(W4  ,W1  ,G,W8  )
  81       CALL FVOXXX(W5  ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W9  )
  82       CALL IOVXXX(W2  ,W9  ,W8  ,GG,AMP(2  ))
  83       CALL FVIXXX(W6  ,W1  ,GG,FMASS(IM1),FWIDTH(IM1),W10 )
  84       CALL IOVXXX(W10 ,W5  ,W4  ,GG,AMP(3  ))
  85       CALL IOVXXX(W6  ,W5  ,W8  ,GG,AMP(4  ))
  86       CALL FVIXXX(W2  ,W4  ,GG,FMASS(IM1),FWIDTH(IM1),W11 )
  87       CALL FVOXXX(W5  ,W1  ,GG,FMASS(IM1),FWIDTH(IM1),W12 )
  88       CALL IOVXXX(W11 ,W12 ,W3  ,GZUD,AMP(5  ))
  89       CALL IOVXXX(W11 ,W9  ,W1  ,GG,AMP(6  ))
  90       CALL FVIXXX(W2  ,W1  ,GG,FMASS(IM1),FWIDTH(IM1),W13 )
  91       CALL FVIXXX(W13 ,W4  ,GG,FMASS(IM1),FWIDTH(IM1),W14 )
  92       CALL IOVXXX(W14 ,W5  ,W3  ,GZUD,AMP(7  ))
  93       CALL FVIXXX(W13 ,W3  ,GZUD,FMASS(IM1),FWIDTH(IM1),W15 )
  94       CALL IOVXXX(W15 ,W5  ,W4  ,GG,AMP(8  ))
  95       ZJJ7 = 0.D0
  96       DO I = 1, NEIGEN
  97           ZTEMP = (0.D0,0.D0)
  98           DO J = 1, NGRAPHS
  99               ZTEMP = ZTEMP + EIGEN_VEC(J,I)*AMP(J)
 100           ENDDO
 101           ZJJ7 =ZJJ7+ZTEMP*EIGEN_VAL(I)*CONJG(ZTEMP)
 102       ENDDO
 103 C      CALL GAUGECHECK(AMP,ZTEMP,EIGEN_VEC,EIGEN_VAL,NGRAPHS,NEIGEN)
 104       END