ISAJET/code/zjj5.F

   1 #include "isajet/pilot.h"
   2       REAL*8 FUNCTION ZJJ5(P1, P2, P3, P4, P5,NHEL,IM1,IM2)
   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 : q(im1) q(im2)  -> z q(im1) q(im2)
  10 C with IM1 != IM2
  11 C
  12 #if defined(CERNLIB_IMPNONE)
  13       IMPLICIT NONE
  14 #endif
  15 C
  16 C CONSTANTS
  17 C
  18       INTEGER    NGRAPHS,    NEIGEN,    NEXTERNAL
  19       PARAMETER (NGRAPHS=  4,NEIGEN=  1,NEXTERNAL=5)
  20       REAL*8     ZERO
  21       PARAMETER (ZERO=0D0)
  22 C
  23 C ARGUMENTS
  24 C
  25       REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
  26       INTEGER NHEL(NEXTERNAL)
  27       INTEGER IM1,IM2
  28 C
  29 C LOCAL VARIABLES
  30 C
  31       INTEGER I,J
  32       REAL*8 EIGEN_VAL(NEIGEN), EIGEN_VEC(NGRAPHS,NEIGEN)
  33       COMPLEX*16 ZTEMP
  34       COMPLEX*16 AMP(NGRAPHS)
  35       COMPLEX*16 W1(6)  , W2(6)  , W3(6)  , W4(6)  , W5(6)
  36       COMPLEX*16 W6(6)  , W7(6)  , W8(6)  , W9(6)  , W10(6)
  37       COMPLEX*16 W11(6) , W12(6)
  38       REAL*8 GZUD1(2),GZUD2(2)
  39 C
  40 C GLOBAL VARIABLES
  41 C
  42 #include "isajet/mgcoms.inc"
  43 C
  44 C COLOR DATA
  45 C
  46       DATA EIGEN_VAL(1  )/       8.8888888888888884D-01 /
  47       DATA EIGEN_VEC(1  ,1  )/  -5.0000000000000000D-01 /
  48       DATA EIGEN_VEC(2  ,1  )/  -5.0000000000000000D-01 /
  49       DATA EIGEN_VEC(3  ,1  )/  -5.0000000000000000D-01 /
  50       DATA EIGEN_VEC(4  ,1  )/  -5.0000000000000000D-01 /
  51 C ----------
  52 C BEGIN CODE
  53 C ----------
  54       IF(IM1.EQ.3.OR.IM1.EQ.7.OR.IM1.EQ.11) THEN
  55         GZUD1(1)=GZD(1)
  56         GZUD1(2)=GZD(2)
  57       ELSE
  58         GZUD1(1)=GZU(1)
  59         GZUD1(2)=GZU(2)
  60       ENDIF
  61       IF(IM2.EQ.3.OR.IM2.EQ.7.OR.IM2.EQ.11) THEN
  62         GZUD2(1)=GZD(1)
  63         GZUD2(2)=GZD(2)
  64       ELSE
  65         GZUD2(1)=GZU(1)
  66         GZUD2(2)=GZU(2)
  67       ENDIF
  68 C
  69       CALL IXXXXX(P1  ,FMASS(IM1),NHEL(1  ), 1,W1  )
  70       CALL IXXXXX(P2  ,FMASS(IM2),NHEL(2  ), 1,W2  )
  71       CALL VXXXXX(P3  ,ZMASS,NHEL(3  ), 1,W3  )
  72       CALL OXXXXX(P4  ,FMASS(IM1),NHEL(4  ), 1,W4  )
  73       CALL OXXXXX(P5  ,FMASS(IM2),NHEL(5  ), 1,W5  )
  74       CALL JIOXXX(W1  ,W4  ,GG,ZERO,ZERO,W6  )
  75       CALL FVOXXX(W5  ,W3  ,GZUD2,FMASS(IM2),FWIDTH(IM2),W7  )
  76       CALL IOVXXX(W2  ,W7  ,W6  ,GG,AMP(1  ))
  77       CALL FVIXXX(W2  ,W3  ,GZUD2,FMASS(IM2),FWIDTH(IM2),W8  )
  78       CALL IOVXXX(W8  ,W5  ,W6  ,GG,AMP(2  ))
  79       CALL FVIXXX(W1  ,W3  ,GZUD1,FMASS(IM1),FWIDTH(IM1),W9  )
  80       CALL JIOXXX(W9  ,W4  ,GG,ZERO,ZERO,W10 )
  81       CALL IOVXXX(W2  ,W5  ,W10 ,GG,AMP(3  ))
  82       CALL FVOXXX(W4  ,W3  ,GZUD1,FMASS(IM1),FWIDTH(IM1),W11 )
  83       CALL JIOXXX(W2  ,W5  ,GG,ZERO,ZERO,W12 )
  84       CALL IOVXXX(W1  ,W11 ,W12 ,GG,AMP(4  ))
  85       ZJJ5 = 0.D0
  86       DO I = 1, NEIGEN
  87           ZTEMP = (0.D0,0.D0)
  88           DO J = 1, NGRAPHS
  89               ZTEMP = ZTEMP + EIGEN_VEC(J,I)*AMP(J)
  90           ENDDO
  91           ZJJ5 =ZJJ5+ZTEMP*EIGEN_VAL(I)*CONJG(ZTEMP)
  92       ENDDO
  93 C      CALL GAUGECHECK(AMP,ZTEMP,EIGEN_VEC,EIGEN_VAL,NGRAPHS,NEIGEN)
  94       END