4 *CMZ : 20/08/91 22.09.33 by Federico Carminati
6 *-- Author : Zoltan Kunszt, modified by Mike Seymour
8 C-----------------------------------------------------------------------
10 SUBROUTINE HWHEW5(N1,N2,N3,N4,N5,N6,HELSUM,HELCTY,ID1,ID2)
12 C-----------------------------------------------------------------------
14 C RECALL THAT N1,N3,N5 MUST BE OUTGOING FERMIONS, AND N2,N4,N6 MUST BE
16 C OUTGOING ANTI-FERMIONS; 3,4 FOR Z0, 5,6 FOR Z0
20 C EQ1 AND T31 ARE FOR OUTOING INITIAL QUARK
22 C CHOOSE APPROPRIATE CASE ACCORDING TO NUPDN
24 C NUPDN=1 FOR UUBAR COLLISIONS, NUPDN=2 FOR DDBAR COLLISIONS
26 C NFINAL CHOOSES THE FINAL DECAYS, 1 FOR DOUBLE LEPTON, 2 FOR 1 FLAVOR
28 C LEPTON+2FAMILIES OF QUARKS, 3 THE SAME, 4 FOR DOUBLE 2FAM3COLOR QUARKS
32 C NOTE: EXTERNAL FACTOR OF COLOR AVERAGE AND SPIN AVERAGE AND
34 C COUPLING (E**8/4/9) MUST BE INCLUDED AS WELL AS COMPENSATION
36 C FOR ON POLE APPROXIMATION AS DESIRED.
40 C---SLIGHTLY MODIFIED BY MHS, SO THAT HELCTY REFERS TO THE FINAL STATE
42 C INDICATED BY ID1,ID2
44 C-----------------------------------------------------------------------
48 COMPLEX HWHEW4,ZH,ZCH,ZD,ZAMM(8),ZS134,ZS156,ZS234,ZS256,ZTWO
50 DOUBLE PRECISION CPFAC,CPALL,HELSUM,HELCTY,AMM
52 INTEGER N1,N2,N3,N4,N5,N6,ID1,ID2,I
56 COMMON/HWHEWQ/ZH(7,7),ZCH(7,7),ZD(7,7)
58 COMMON/HWHEWR/CPFAC(12,12,8),CPALL(8)
62 C THE MATRIX ELEMENT DEPENDS ON
64 ZS134=(ZD(N1,N3)+ZD(N1,N4)+ZD(N3,N4))*ZTWO
66 ZS156=(ZD(N1,N5)+ZD(N1,N6)+ZD(N5,N6))*ZTWO
68 ZS234=(ZD(N2,N3)+ZD(N2,N4)+ZD(N3,N4))*ZTWO
70 ZS256=(ZD(N2,N5)+ZD(N2,N6)+ZD(N5,N6))*ZTWO
72 ZAMM(1)=HWHEW4(N1,N2,N3,N4,N5,N6)/ZS134+
74 > HWHEW4(N1,N2,N5,N6,N3,N4)/ZS156
76 ZAMM(2)=HWHEW4(N1,N2,N4,N3,N5,N6)/ZS134+
78 > HWHEW4(N1,N2,N5,N6,N4,N3)/ZS156
80 ZAMM(3)=HWHEW4(N1,N2,N3,N4,N6,N5)/ZS134+
82 > HWHEW4(N1,N2,N6,N5,N3,N4)/ZS156
84 ZAMM(4)=HWHEW4(N1,N2,N4,N3,N6,N5)/ZS134+
86 > HWHEW4(N1,N2,N6,N5,N4,N3)/ZS156
88 ZAMM(5)=HWHEW4(N2,N1,N3,N4,N5,N6)/ZS234+
90 > HWHEW4(N2,N1,N5,N6,N3,N4)/ZS256
92 ZAMM(6)=HWHEW4(N2,N1,N4,N3,N5,N6)/ZS234+
94 > HWHEW4(N2,N1,N5,N6,N4,N3)/ZS256
96 ZAMM(7)=HWHEW4(N2,N1,N3,N4,N6,N5)/ZS234+
98 > HWHEW4(N2,N1,N6,N5,N3,N4)/ZS256
100 ZAMM(8)=HWHEW4(N2,N1,N4,N3,N6,N5)/ZS234+
102 > HWHEW4(N2,N1,N6,N5,N4,N3)/ZS256
110 AMM=REAL(ZAMM(I)*CONJG(ZAMM(I)))
112 HELSUM=HELSUM+CPALL(I)*AMM
114 HELCTY=HELCTY+CPFAC(ID1,ID2,I)*AMM