ISAJET/code/wzsu.F

   1 #include "isajet/pilot.h"
   2       FUNCTION WZSU(T,U,T1,U1,T3,U3,P1,P2)
   3 C          DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
   4 C          INTERFERENCE OF S AND U GRAPHS.
   5 #include "isajet/itapes.inc"
   6 #include "isajet/wwpar.inc"
   7       DIMENSION P1(4),P2(4)
   8 #if defined(CERNLIB_DOUBLE)
   9       DOUBLE PRECISION WZSU
  10       DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
  11       DOUBLE PRECISION WM4,ZM4,WZM2,CSXCU,EPF
  12 #endif
  13       WM4=WM2**2
  14       ZM4=ZM2**2
  15       WZM2=WM2*ZM2
  16       CSXCU=CS*CU/U
  17       WZSU=
  18      1 +CSXCU*CV3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-64.*WM2*T1*ZM4+128.
  19      1 *WM2*U1*T3*U3-64.*WM2*T3*ZM4-64.*WM2*S13*ZM4+64.*ZM2*T1*U1*T3+64
  20      1 .*ZM2*T1*U1*U3-64.*ZM2*T1*WM4-128.*ZM2*U1**2*U3-64.*ZM2*T3*WM4-6
  21      1 4.*ZM2*S13*WM4+128.*T1*U1*T3*U3-64.*T1*T3*WZM2-64.*T1*S13*WZM2-3
  22      1 2.*T1**2*WZM2+64.*U1*U3*WZM2-64.*U1*S13*WZM2-128.*U1**2*U3**2-32
  23      1 .*U1**2*WZM2-64.*U1**2*ZM4-128.*T3*S13*WZM2-64.*T3**2*WZM2-64.*S
  24      1 13**2*WZM2-64.*WM4*ZM4)
  25       WZSU=WZSU
  26      1 +CSXCU*CV3*T*(-32.*WM2*U1*U3+64.*WM2*ZM4-32.*ZM2*U1*T3-32.*ZM2*U
  27      1 1*U3+32.*ZM2*U1*S13+32.*ZM2*U1**2+64.*ZM2*WM4-32.*T1*U1*U3+64.*T
  28      1 1*WZM2-64.*U1*T3*U3+64.*U1*U3*S13+32.*U1**2*U3+64.*T3*WZM2+64.*S
  29      1 13*WZM2)
  30      1 +CSXCU*CV3*T*U*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-3
  31      1 2.*T1*T3+32.*T1*S13+64.*U1*T3+64.*U1*U3+32.*U1*S13+64.*T3*S13+64
  32      1 .*S13**2-32.*WZM2)
  33       WZSU=WZSU
  34      1 +CSXCU*CV3*T**2*(32.*U1*U3-32.*WZM2)
  35      1 +CSXCU*CV3*T**2*U*(-32.*S13)
  36       WZSU=WZSU
  37      1 +CSXCU*CV3*U*(64.*WM2*T1*T3-32.*WM2*U1*T3-64.*WM2*U1*U3+64.*WM2*
  38      1 T3*S13+64.*WM2*T3**2+32.*WM2*ZM4+64.*ZM2*T1*T3+32.*ZM2*T1*S13+32
  39      1 .*ZM2*T1**2-64.*ZM2*U1*T3-64.*ZM2*U1*U3-64.*ZM2*U1*S13+64.*ZM2*U
  40      1 1**2+32.*ZM2*WM4-32.*T1*U1*T3-64.*T1*U1*U3+64.*T1*T3*S13+64.*T1*
  41      1 T3**2+64.*T1*WZM2+32.*T1*ZM4+32.*T1**2*T3-128.*U1*T3*U3-128.*U1*
  42      1 U3*S13+64.*U1**2*U3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM2)
  43       WZSU=WZSU
  44      1 +CSXCU*CV3*U**2*(64.*U1*T3+64.*U1*S13)
  45      1 +CSXCU*CA3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+64.*WM2*T1*ZM4+64.*W
  46      1 M2*T3*ZM4+64.*WM2*S13*ZM4+64.*ZM2*T1*U1*T3-64.*ZM2*T1*U1*U3+64.*
  47      1 ZM2*T1*WM4-128.*ZM2*U1**2*U3+64.*ZM2*T3*WM4+64.*ZM2*S13*WM4+64.*
  48      1 T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2-64.*U1*U3*WZM2-64.*U1*
  49      1 S13*WZM2-32.*U1**2*WZM2-64.*U1**2*ZM4+64.*WM4*ZM4)
  50       WZSU=WZSU
  51      1 +CSXCU*CA3*T*(32.*WM2*U1*U3-64.*WM2*ZM4-32.*ZM2*U1*T3+32.*ZM2*U1
  52      1 *U3+32.*ZM2*U1*S13+32.*ZM2*U1**2-64.*ZM2*WM4+32.*T1*U1*U3-64.*T1
  53      1 *WZM2+32.*U1**2*U3-64.*T3*WZM2-64.*S13*WZM2)
  54      1 +CSXCU*CA3*T*U*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+32
  55      1 .*T1*T3-32.*T1*S13+64.*U1*T3-64.*U1*U3+32.*U1*S13+32.*WZM2)
  56      1 +CSXCU*CA3*T**2*(-32.*U1*U3+32.*WZM2)
  57       WZSU=WZSU
  58      1 +CSXCU*CA3*T**2*U*(32.*S13)
  59      1 +CSXCU*CA3*U*(-64.*WM2*T1*T3-32.*WM2*U1*T3+64.*WM2*U1*U3-32.*WM2
  60      1 *ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-64.*ZM2*U1*T3+64
  61      1 .*ZM2*U1*U3-64.*ZM2*U1*S13+64.*ZM2*U1**2-32.*ZM2*WM4-32.*T1*U1*T
  62      1 3+64.*T1*U1*U3-64.*T1*WZM2-32.*T1*ZM4-32.*T1**2*T3+64.*U1**2*U3-
  63      1 64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
  64       WZSU=WZSU
  65      1 +CSXCU*CA3*U**2*(64.*U1*T3+64.*U1*S13)
  66      1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*(32.*ZM2*T1+32.*ZM2*U1+32.*WZM2)
  67      1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*T*(-16.*ZM2)
  68      1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*U*(-16.*ZM2)
  69      1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*(64.*WM2*T1+64.*WM2*U1+64.*WM2*S13+6
  70      1 4.*T1*U1+64.*T1*S13+32.*T1**2+64.*U1*S13+32.*U1**2+32.*WM4)
  71       WZSU=WZSU
  72      1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*T*(-16.*WM2-16.*T1-16.*U1-32.*S13)
  73      1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*U*(-16.*WM2-16.*T1-16.*U1-32.*S13)
  74      1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*(64.*WM2*T3+64.*WM2*U3+64.*ZM2*T1+64
  75      1 .*ZM2*U1+64.*T1*T3+64.*T1*U3+64.*U1*T3+64.*U1*U3+64.*WZM2)
  76      1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T*(-32.*WM2-32.*ZM2-32.*T1-32.*U1-32
  77      1 .*T3-32.*U3)
  78       WZSU=WZSU
  79      1 +32.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T*U
  80      1 +16.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T**2
  81      1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*U*(-32.*WM2-32.*ZM2-32.*T1-32.*U1-32
  82      1 .*T3-32.*U3)
  83      1 +16.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*U**2
  84      1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*(-64.*WM2*T3-32.*ZM2*T1-32.*ZM2*U1-6
  85      1 4.*T1*T3-64.*U1*T3-32.*WZM2)
  86       WZSU=WZSU
  87      1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T*(32.*WM2+16.*ZM2+32.*T1+32.*U1+32.
  88      1 *T3)
  89      1 -16.*EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T*U
  90      1 -16.*EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T**2
  91      1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*U*(16.*ZM2+32.*T3)
  92      1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*(64.*WM2*U3+32.*ZM2*T1+32.*ZM2*U1+64
  93      1 .*T1*U3+64.*U1*U3+32.*WZM2)
  94       WZSU=WZSU
  95      1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*T*(-16.*ZM2-32.*U3)
  96      1 +16.*EPF(P2,P3,Q1,Q3)*CSXCU*CV3*T*U
  97      1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*U*(-32.*WM2-16.*ZM2-32.*T1-32.*U1-32
  98      1 .*U3)
  99      1 +16.*EPF(P2,P3,Q1,Q3)*CSXCU*CV3*U**2
 100       RETURN
 101       END