#include "isajet/pilot.h" FUNCTION ZZALL(T,U,T1,U1,T3,U3,P1,P2) C DECAY DISTRIBUTION FOR Z0 Z0 PAIRS FROM SCHOONSCHIP(1980). C INCLUDES TT, TU, AND UU TERMS. #include "isajet/itapes.inc" #include "isajet/wwpar.inc" DIMENSION P1(4),P2(4) #if defined(CERNLIB_DOUBLE) DOUBLE PRECISION ZZALL DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2 DOUBLE PRECISION CVAVA,DVAVA,CVVVV,CVVAA,TT,TU,UU,EPF #endif TTUU=(T*U)**2 CVAVA=CV*CA*(CV1*CA3+CA1*CV3)/TTUU DVAVA=CV*CA*(CV1*CA3-CA1*CV3)/TTUU CVVVV=(CV**2+CA**2)*CV1*CV3/TTUU CVVAA=(CV**2+CA**2)*CA1*CA3/TTUU TT=T**2 TU=T*U UU=U**2 ZZALL= 1 +CVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3+16.*ZM2*T1*T3*U3- 1 16.*ZM2*U1*T3*U3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1**2+16.*ZM2**2*U1 1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2-16.*ZM2** 1 2*U3*S13-8.*ZM2**2*U3**2) 1 +CVAVA*T*U**2*(16.*ZM2*T1*S13-16.*ZM2*T3*S13) 1 +CVAVA*T**2*(16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-8.*ZM2**2*U1**2+8 1 .*ZM2**2*U3**2) 1 +CVAVA*T**2*U*(-16.*ZM2*U1*S13+16.*ZM2*U3*S13) 1 +CVAVA*U**2*(-16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+8.*ZM2**2*T1**2- 1 8.*ZM2**2*T3**2) ZZALL=ZZALL 1 +DVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3+ 1 16.*ZM2*U1*T3*U3-32.*ZM2**2*T1*T3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1 1 **2+32.*ZM2**2*U1*U3+16.*ZM2**2*U1*S13+8.*ZM2**2*U1**2-16.*ZM2** 1 2*T3*S13-8.*ZM2**2*T3**2+16.*ZM2**2*U3*S13+8.*ZM2**2*U3**2-32.*Z 1 M2**3*T1+32.*ZM2**3*U1-32.*ZM2**3*T3+32.*ZM2**3*U3) 1 +DVAVA*T*U**2*(16.*ZM2*T1*S13+16.*ZM2*T3*S13-16.*ZM2**2*U1-16.*Z 1 M2**2*U3-32.*ZM2**3) 1 +DVAVA*T*U**3*(8.*ZM2**2) ZZALL=ZZALL 1 +DVAVA*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-32.*ZM2**2*U1*U3 1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-16.*ZM2**3*U1-16.*ZM2**3*U3-8.* 1 ZM2**4) 1 +DVAVA*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13+16.*ZM2**2*T1+16.* 1 ZM2**2*T3+32.*ZM2**3) 1 +DVAVA*T**3*U*(-8.*ZM2**2) 1 +DVAVA*U**2*(16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+32.*ZM2**2*T1*T3+ 1 8.*ZM2**2*T1**2+8.*ZM2**2*T3**2+16.*ZM2**3*T1+16.*ZM2**3*T3+8.*Z 1 M2**4) ZZALL=ZZALL 1 +CVVVV*T*U*(-16.*ZM2*T1*U1*T3-16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3- 1 16.*ZM2*U1*T3*U3+16.*ZM2**2*T1*S13+8.*ZM2**2*T1**2+16.*ZM2**2*U1 1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2+16.*ZM2** 1 2*U3*S13+8.*ZM2**2*U3**2+16.*ZM2**2*S13**2+16.*ZM2**3*T1+16.*ZM2 1 **3*U1+16.*ZM2**3*T3+16.*ZM2**3*U3+32.*ZM2**3*S13+32.*ZM2**4-32. 1 *T1*U1*T3*U3) 1 +CVVVV*T*U**2*(-16.*ZM2*T1*S13-16.*ZM2*T3*S13-8.*ZM2**2*U1-8.*ZM 1 2**2*U3-16.*ZM2**2*S13-16.*ZM2**3-32.*T1*T3*S13) ZZALL=ZZALL 1 +CVVVV*T*U**3*(4.*ZM2**2) 1 +CVVVV*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-16.*ZM2**2*U1*U3 1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-8.*ZM2**3*U1-8.*ZM2**3*U3-4.*ZM 1 2**4-16.*U1**2*U3**2) 1 +CVVVV*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13-8.*ZM2**2*T1-8.*ZM 1 2**2*T3-16.*ZM2**2*S13-16.*ZM2**3-32.*U1*U3*S13) 1 +CVVVV*T**2*U**2*(-16.*S13**2) 1 +CVVVV*T**3*U*(4.*ZM2**2) ZZALL=ZZALL 1 +CVVVV*U**2*(-16.*ZM2*T1*T3**2-16.*ZM2*T1**2*T3-16.*ZM2**2*T1*T3 1 -8.*ZM2**2*T1**2-8.*ZM2**2*T3**2-8.*ZM2**3*T1-8.*ZM2**3*T3-4.*ZM 1 2**4-16.*T1**2*T3**2) 1 +CVVAA*T*U*(-16.*ZM2**3*T1-16.*ZM2**3*U1-16.*ZM2**3*T3-16.*ZM2** 1 3*U3-32.*ZM2**3*S13-32.*ZM2**4) 1 +CVVAA*T*U**2*(8.*ZM2**2*U1+8.*ZM2**2*U3+16.*ZM2**2*S13+16.*ZM2* 1 *3) 1 +CVVAA*T*U**3*(-4.*ZM2**2) ZZALL=ZZALL 1 +CVVAA*T**2*(16.*ZM2**2*U1*U3+8.*ZM2**3*U1+8.*ZM2**3*U3+4.*ZM2** 1 4) 1 +CVVAA*T**2*U*(8.*ZM2**2*T1+8.*ZM2**2*T3+16.*ZM2**2*S13+16.*ZM2* 1 *3) 1 +CVVAA*T**3*U*(-4.*ZM2**2) 1 +CVVAA*U**2*(16.*ZM2**2*T1*T3+8.*ZM2**3*T1+8.*ZM2**3*T3+4.*ZM2** 1 4) C THE EPF TERMS FROM SCHOONSCHIP DID NOT CONTAIN CVAVA, ETC., C BUT DID CONTAIN AN EXPLICIT T*U. THIS WAS REPLACED BY HAND C BY 1./TU. ZZALL=ZZALL 1 +EPF(P1,P2,P3,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*S13) 1 +EPF(P1,P2,Q1,Q3)/TU*(-32.*ZM2*CV*CA*CV1*CV3*T1-32.*CV*CA*CV1*C 1 V3*T1*T3-32.*CV*CA*CV1*CV3*T1*U3) 1 +EPF(P1,P3,Q1,Q3)/TU*(32.*CV*CA*CV1*CV3*T1*T3) 1 +EPF(P2,P3,Q1,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*U3)+0. ZZALL=2.*ZZALL RETURN END