#include "isajet/pilot.h" FUNCTION WZSU(T,U,T1,U1,T3,U3,P1,P2) C DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980). C INTERFERENCE OF S AND U GRAPHS. #include "isajet/itapes.inc" #include "isajet/wwpar.inc" DIMENSION P1(4),P2(4) #if defined(CERNLIB_DOUBLE) DOUBLE PRECISION WZSU DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2 DOUBLE PRECISION WM4,ZM4,WZM2,CSXCU,EPF #endif WM4=WM2**2 ZM4=ZM2**2 WZM2=WM2*ZM2 CSXCU=CS*CU/U WZSU= 1 +CSXCU*CV3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-64.*WM2*T1*ZM4+128. 1 *WM2*U1*T3*U3-64.*WM2*T3*ZM4-64.*WM2*S13*ZM4+64.*ZM2*T1*U1*T3+64 1 .*ZM2*T1*U1*U3-64.*ZM2*T1*WM4-128.*ZM2*U1**2*U3-64.*ZM2*T3*WM4-6 1 4.*ZM2*S13*WM4+128.*T1*U1*T3*U3-64.*T1*T3*WZM2-64.*T1*S13*WZM2-3 1 2.*T1**2*WZM2+64.*U1*U3*WZM2-64.*U1*S13*WZM2-128.*U1**2*U3**2-32 1 .*U1**2*WZM2-64.*U1**2*ZM4-128.*T3*S13*WZM2-64.*T3**2*WZM2-64.*S 1 13**2*WZM2-64.*WM4*ZM4) WZSU=WZSU 1 +CSXCU*CV3*T*(-32.*WM2*U1*U3+64.*WM2*ZM4-32.*ZM2*U1*T3-32.*ZM2*U 1 1*U3+32.*ZM2*U1*S13+32.*ZM2*U1**2+64.*ZM2*WM4-32.*T1*U1*U3+64.*T 1 1*WZM2-64.*U1*T3*U3+64.*U1*U3*S13+32.*U1**2*U3+64.*T3*WZM2+64.*S 1 13*WZM2) 1 +CSXCU*CV3*T*U*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-3 1 2.*T1*T3+32.*T1*S13+64.*U1*T3+64.*U1*U3+32.*U1*S13+64.*T3*S13+64 1 .*S13**2-32.*WZM2) WZSU=WZSU 1 +CSXCU*CV3*T**2*(32.*U1*U3-32.*WZM2) 1 +CSXCU*CV3*T**2*U*(-32.*S13) WZSU=WZSU 1 +CSXCU*CV3*U*(64.*WM2*T1*T3-32.*WM2*U1*T3-64.*WM2*U1*U3+64.*WM2* 1 T3*S13+64.*WM2*T3**2+32.*WM2*ZM4+64.*ZM2*T1*T3+32.*ZM2*T1*S13+32 1 .*ZM2*T1**2-64.*ZM2*U1*T3-64.*ZM2*U1*U3-64.*ZM2*U1*S13+64.*ZM2*U 1 1**2+32.*ZM2*WM4-32.*T1*U1*T3-64.*T1*U1*U3+64.*T1*T3*S13+64.*T1* 1 T3**2+64.*T1*WZM2+32.*T1*ZM4+32.*T1**2*T3-128.*U1*T3*U3-128.*U1* 1 U3*S13+64.*U1**2*U3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM2) WZSU=WZSU 1 +CSXCU*CV3*U**2*(64.*U1*T3+64.*U1*S13) 1 +CSXCU*CA3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+64.*WM2*T1*ZM4+64.*W 1 M2*T3*ZM4+64.*WM2*S13*ZM4+64.*ZM2*T1*U1*T3-64.*ZM2*T1*U1*U3+64.* 1 ZM2*T1*WM4-128.*ZM2*U1**2*U3+64.*ZM2*T3*WM4+64.*ZM2*S13*WM4+64.* 1 T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2-64.*U1*U3*WZM2-64.*U1* 1 S13*WZM2-32.*U1**2*WZM2-64.*U1**2*ZM4+64.*WM4*ZM4) WZSU=WZSU 1 +CSXCU*CA3*T*(32.*WM2*U1*U3-64.*WM2*ZM4-32.*ZM2*U1*T3+32.*ZM2*U1 1 *U3+32.*ZM2*U1*S13+32.*ZM2*U1**2-64.*ZM2*WM4+32.*T1*U1*U3-64.*T1 1 *WZM2+32.*U1**2*U3-64.*T3*WZM2-64.*S13*WZM2) 1 +CSXCU*CA3*T*U*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+32 1 .*T1*T3-32.*T1*S13+64.*U1*T3-64.*U1*U3+32.*U1*S13+32.*WZM2) 1 +CSXCU*CA3*T**2*(-32.*U1*U3+32.*WZM2) WZSU=WZSU 1 +CSXCU*CA3*T**2*U*(32.*S13) 1 +CSXCU*CA3*U*(-64.*WM2*T1*T3-32.*WM2*U1*T3+64.*WM2*U1*U3-32.*WM2 1 *ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-64.*ZM2*U1*T3+64 1 .*ZM2*U1*U3-64.*ZM2*U1*S13+64.*ZM2*U1**2-32.*ZM2*WM4-32.*T1*U1*T 1 3+64.*T1*U1*U3-64.*T1*WZM2-32.*T1*ZM4-32.*T1**2*T3+64.*U1**2*U3- 1 64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2) WZSU=WZSU 1 +CSXCU*CA3*U**2*(64.*U1*T3+64.*U1*S13) 1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*(32.*ZM2*T1+32.*ZM2*U1+32.*WZM2) 1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*T*(-16.*ZM2) 1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*U*(-16.*ZM2) 1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*(64.*WM2*T1+64.*WM2*U1+64.*WM2*S13+6 1 4.*T1*U1+64.*T1*S13+32.*T1**2+64.*U1*S13+32.*U1**2+32.*WM4) WZSU=WZSU 1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*T*(-16.*WM2-16.*T1-16.*U1-32.*S13) 1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*U*(-16.*WM2-16.*T1-16.*U1-32.*S13) 1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*(64.*WM2*T3+64.*WM2*U3+64.*ZM2*T1+64 1 .*ZM2*U1+64.*T1*T3+64.*T1*U3+64.*U1*T3+64.*U1*U3+64.*WZM2) 1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T*(-32.*WM2-32.*ZM2-32.*T1-32.*U1-32 1 .*T3-32.*U3) WZSU=WZSU 1 +32.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T*U 1 +16.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T**2 1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*U*(-32.*WM2-32.*ZM2-32.*T1-32.*U1-32 1 .*T3-32.*U3) 1 +16.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*U**2 1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*(-64.*WM2*T3-32.*ZM2*T1-32.*ZM2*U1-6 1 4.*T1*T3-64.*U1*T3-32.*WZM2) WZSU=WZSU 1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T*(32.*WM2+16.*ZM2+32.*T1+32.*U1+32. 1 *T3) 1 -16.*EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T*U 1 -16.*EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T**2 1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*U*(16.*ZM2+32.*T3) 1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*(64.*WM2*U3+32.*ZM2*T1+32.*ZM2*U1+64 1 .*T1*U3+64.*U1*U3+32.*WZM2) WZSU=WZSU 1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*T*(-16.*ZM2-32.*U3) 1 +16.*EPF(P2,P3,Q1,Q3)*CSXCU*CV3*T*U 1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*U*(-32.*WM2-16.*ZM2-32.*T1-32.*U1-32 1 .*U3) 1 +16.*EPF(P2,P3,Q1,Q3)*CSXCU*CV3*U**2 RETURN END