HERWIG/src/hwhew5.f

   1
   2 CDECK  ID>, HWHEW5.
   3
   4 *CMZ :          20/08/91  22.09.33  by  Federico Carminati
   5
   6 *-- Author :    Zoltan Kunszt, modified by Mike Seymour
   7
   8 C-----------------------------------------------------------------------
   9
  10       SUBROUTINE HWHEW5(N1,N2,N3,N4,N5,N6,HELSUM,HELCTY,ID1,ID2)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C RECALL THAT N1,N3,N5 MUST BE OUTGOING FERMIONS, AND N2,N4,N6 MUST BE
  15
  16 C OUTGOING ANTI-FERMIONS; 3,4 FOR Z0, 5,6 FOR Z0
  17
  18 C
  19
  20 C EQ1 AND T31 ARE FOR OUTOING INITIAL QUARK
  21
  22 C CHOOSE APPROPRIATE CASE ACCORDING TO NUPDN
  23
  24 C NUPDN=1 FOR UUBAR COLLISIONS, NUPDN=2 FOR DDBAR COLLISIONS
  25
  26 C NFINAL CHOOSES THE FINAL DECAYS, 1 FOR DOUBLE LEPTON, 2 FOR 1 FLAVOR
  27
  28 C LEPTON+2FAMILIES OF QUARKS, 3 THE SAME, 4 FOR DOUBLE 2FAM3COLOR QUARKS
  29
  30 C
  31
  32 C NOTE: EXTERNAL FACTOR OF COLOR AVERAGE AND SPIN AVERAGE AND
  33
  34 C COUPLING (E**8/4/9) MUST BE INCLUDED AS WELL AS COMPENSATION
  35
  36 C FOR ON POLE APPROXIMATION AS DESIRED.
  37
  38 C
  39
  40 C---SLIGHTLY MODIFIED BY MHS, SO THAT HELCTY REFERS TO THE FINAL STATE
  41
  42 C   INDICATED BY ID1,ID2
  43
  44 C-----------------------------------------------------------------------
  45
  46       IMPLICIT NONE
  47
  48       COMPLEX HWHEW4,ZH,ZCH,ZD,ZAMM(8),ZS134,ZS156,ZS234,ZS256,ZTWO
  49
  50       DOUBLE PRECISION CPFAC,CPALL,HELSUM,HELCTY,AMM
  51
  52       INTEGER N1,N2,N3,N4,N5,N6,ID1,ID2,I
  53
  54       EXTERNAL HWHEW4
  55
  56       COMMON/HWHEWQ/ZH(7,7),ZCH(7,7),ZD(7,7)
  57
  58       COMMON/HWHEWR/CPFAC(12,12,8),CPALL(8)
  59
  60       DATA ZTWO/(2.0,0.0)/
  61
  62 C THE MATRIX ELEMENT DEPENDS ON
  63
  64       ZS134=(ZD(N1,N3)+ZD(N1,N4)+ZD(N3,N4))*ZTWO
  65
  66       ZS156=(ZD(N1,N5)+ZD(N1,N6)+ZD(N5,N6))*ZTWO
  67
  68       ZS234=(ZD(N2,N3)+ZD(N2,N4)+ZD(N3,N4))*ZTWO
  69
  70       ZS256=(ZD(N2,N5)+ZD(N2,N6)+ZD(N5,N6))*ZTWO
  71
  72       ZAMM(1)=HWHEW4(N1,N2,N3,N4,N5,N6)/ZS134+
  73
  74      >        HWHEW4(N1,N2,N5,N6,N3,N4)/ZS156
  75
  76       ZAMM(2)=HWHEW4(N1,N2,N4,N3,N5,N6)/ZS134+
  77
  78      >        HWHEW4(N1,N2,N5,N6,N4,N3)/ZS156
  79
  80       ZAMM(3)=HWHEW4(N1,N2,N3,N4,N6,N5)/ZS134+
  81
  82      >        HWHEW4(N1,N2,N6,N5,N3,N4)/ZS156
  83
  84       ZAMM(4)=HWHEW4(N1,N2,N4,N3,N6,N5)/ZS134+
  85
  86      >        HWHEW4(N1,N2,N6,N5,N4,N3)/ZS156
  87
  88       ZAMM(5)=HWHEW4(N2,N1,N3,N4,N5,N6)/ZS234+
  89
  90      >        HWHEW4(N2,N1,N5,N6,N3,N4)/ZS256
  91
  92       ZAMM(6)=HWHEW4(N2,N1,N4,N3,N5,N6)/ZS234+
  93
  94      >        HWHEW4(N2,N1,N5,N6,N4,N3)/ZS256
  95
  96       ZAMM(7)=HWHEW4(N2,N1,N3,N4,N6,N5)/ZS234+
  97
  98      >        HWHEW4(N2,N1,N6,N5,N3,N4)/ZS256
  99
 100       ZAMM(8)=HWHEW4(N2,N1,N4,N3,N6,N5)/ZS234+
 101
 102      >        HWHEW4(N2,N1,N6,N5,N4,N3)/ZS256
 103
 104       HELSUM=0.0
 105
 106       HELCTY=0.0
 107
 108       DO 1 I=1,8
 109
 110         AMM=REAL(ZAMM(I)*CONJG(ZAMM(I)))
 111
 112         HELSUM=HELSUM+CPALL(I)*AMM
 113
 114         HELCTY=HELCTY+CPFAC(ID1,ID2,I)*AMM
 115
 116  1    CONTINUE
 117
 118       RETURN
 119
 120       END