ISAJET/code/wztu.F

   1 #include "isajet/pilot.h"
   2       FUNCTION WZTU(T,U,T1,U1,T3,U3,P1,P2)
   3 C          DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
   4 C          ALL T AND U GRAPH TERMS.
   5 #include "isajet/itapes.inc"
   6 #include "isajet/wwpar.inc"
   7       DIMENSION P1(4),P2(4)
   8 #if defined(CERNLIB_DOUBLE)
   9       DOUBLE PRECISION WZTU
  10       DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
  11       DOUBLE PRECISION WM4,ZM4,WZM2,CTXCT,CTXCU,CUXCU,EPF
  12 #endif
  13       WM4=WM2**2
  14       ZM4=ZM2**2
  15       WZM2=WM2*ZM2
  16       CTXCT=CT**2/T**2
  17       CTXCU=CT*CU/(T*U)
  18       CUXCU=CU**2/U**2
  19       WZTU=
  20      1 +CTXCT*CV3*(-128.*WM2*T1*T3**2-64.*WM2*T1*ZM4-64.*ZM2*T1**2*T3-6
  21      1 4.*ZM2*T3*WM4-128.*T1*T3*WZM2-64.*T1**2*T3**2-32.*T1**2*ZM4-64.*
  22      1 T3**2*WM4-32.*WM4*ZM4)
  23       WZTU=WZTU
  24      1 +CTXCT*CV3*T*(64.*WM2*T1*T3-32.*WM2*U1*T3-64.*WM2*T3*S13+64.*WM2
  25      1 *T3**2+32.*WM2*ZM4+64.*ZM2*T1*T3-32.*ZM2*T1*S13+32.*ZM2*T1**2+32
  26      1 .*ZM2*WM4-32.*T1*U1*T3-64.*T1*T3*S13+64.*T1*T3**2+64.*T1*WZM2+32
  27      1 .*T1*ZM4+32.*T1**2*T3+64.*T3*WZM2+32.*T3*WM4-32.*S13*WZM2)
  28      1 +CTXCT*CV3*T**2*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-
  29      1 32.*T1*T3+32.*T1*S13+32.*U1*T3+64.*T3*S13-32.*WZM2)
  30       WZTU=WZTU
  31      1 +CTXCT*CV3*T**3*(-32.*S13)
  32      1 +CTXCT*CA3*(64.*WM2*T1*ZM4+64.*ZM2*T1**2*T3+64.*ZM2*T3*WM4+128.*
  33      1 T1*T3*WZM2+32.*T1**2*ZM4+32.*WM4*ZM4)
  34      1 +CTXCT*CA3*T*(-64.*WM2*T1*T3-32.*WM2*U1*T3-32.*WM2*ZM4-64.*ZM2*T
  35      1 1*T3+32.*ZM2*T1*S13-32.*ZM2*T1**2-32.*ZM2*WM4-32.*T1*U1*T3-64.*T
  36      1 1*WZM2-32.*T1*ZM4-32.*T1**2*T3-64.*T3*WZM2-32.*T3*WM4+32.*S13*WZ
  37      1 M2)
  38       WZTU=WZTU
  39      1 +CTXCT*CA3*T**2*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+3
  40      1 2.*T1*T3-32.*T1*S13+32.*U1*T3+32.*WZM2)
  41      1 +CTXCT*CA3*T**3*(32.*S13)
  42       WZTU=WZTU
  43      1 +CTXCU*CV3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+64.*WM2*T1*ZM4-128.*
  44      1 WM2*U1*T3*U3+64.*WM2*T3*ZM4+64.*WM2*S13*ZM4-64.*ZM2*T1*U1*T3-64.
  45      1 *ZM2*T1*U1*U3+64.*ZM2*T1*WM4+64.*ZM2*T3*WM4+64.*ZM2*S13*WM4-128.
  46      1 *T1*U1*T3*U3+64.*T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2-64.*U
  47      1 1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**2*WZM2+128.*T3*S13*WZM2+64.*T3
  48      1 **2*WZM2+64.*S13**2*WZM2+64.*WM4*ZM4)
  49       WZTU=WZTU
  50      1 +CTXCU*CV3*T*(32.*WM2*U1*U3-64.*WM2*ZM4+32.*ZM2*U1*T3+32.*ZM2*U1
  51      1 *U3-32.*ZM2*U1*S13-32.*ZM2*U1**2-64.*ZM2*WM4+32.*T1*U1*U3-64.*T1
  52      1 *WZM2+64.*U1*T3*U3-64.*U1*U3*S13-32.*U1**2*U3-64.*T3*WZM2-64.*S1
  53      1 3*WZM2)
  54      1 +CTXCU*CV3*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-32.*U1*S13-64.*T3*S13-64.*S13**2+32
  56      1 .*WZM2)
  57       WZTU=WZTU
  58      1 +CTXCU*CV3*T**2*(-32.*U1*U3+32.*WZM2)
  59      1 +CTXCU*CV3*T**2*U*(32.*S13)
  60      1 +CTXCU*CV3*U*(-64.*WM2*T1*T3+32.*WM2*U1*T3-64.*WM2*T3*S13-64.*WM
  61      1 2*T3**2-32.*WM2*ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-3
  62      1 2.*ZM2*WM4+32.*T1*U1*T3-64.*T1*T3*S13-64.*T1*T3**2-64.*T1*WZM2-3
  63      1 2.*T1*ZM4-32.*T1**2*T3-64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
  64       WZTU=WZTU
  65      1 +CTXCU*CA3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-64.*WM2*T1*ZM4-64.*
  66      1 WM2*T3*ZM4-64.*WM2*S13*ZM4-64.*ZM2*T1*U1*T3+64.*ZM2*T1*U1*U3-64.
  67      1 *ZM2*T1*WM4-64.*ZM2*T3*WM4-64.*ZM2*S13*WM4-64.*T1*T3*WZM2-64.*T1
  68      1 *S13*WZM2-32.*T1**2*WZM2+64.*U1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**
  69      1 2*WZM2-64.*WM4*ZM4)
  70       WZTU=WZTU
  71      1 +CTXCU*CA3*T*(-32.*WM2*U1*U3+64.*WM2*ZM4+32.*ZM2*U1*T3-32.*ZM2*U
  72      1 1*U3-32.*ZM2*U1*S13-32.*ZM2*U1**2+64.*ZM2*WM4-32.*T1*U1*U3+64.*T
  73      1 1*WZM2-32.*U1**2*U3+64.*T3*WZM2+64.*S13*WZM2)
  74      1 +CTXCU*CA3*T*U*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-3
  75      1 2.*T1*T3+32.*T1*S13-64.*U1*T3-32.*U1*S13-32.*WZM2)
  76      1 +CTXCU*CA3*T**2*(32.*U1*U3-32.*WZM2)
  77       WZTU=WZTU
  78      1 +CTXCU*CA3*T**2*U*(-32.*S13)
  79      1 +CTXCU*CA3*U*(64.*WM2*T1*T3+32.*WM2*U1*T3+32.*WM2*ZM4+64.*ZM2*T1
  80      1 *T3+32.*ZM2*T1*S13+32.*ZM2*T1**2+32.*ZM2*WM4+32.*T1*U1*T3+64.*T1
  81      1 *WZM2+32.*T1*ZM4+32.*T1**2*T3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM
  82      1 2)
  83      1 +CUXCU*CV3*(-64.*ZM2*U1**2*U3-64.*U1**2*U3**2-32.*U1**2*ZM4)
  84       WZTU=WZTU
  85      1 +CUXCU*CV3*T*U*(32.*U1*U3)
  86      1 +CUXCU*CV3*U*(-32.*WM2*U1*U3-32.*ZM2*U1*T3-32.*ZM2*U1*U3-32.*ZM2
  87      1 *U1*S13+32.*ZM2*U1**2-32.*T1*U1*U3-64.*U1*T3*U3-64.*U1*U3*S13+32
  88      1 .*U1**2*U3)
  89      1 +CUXCU*CV3*U**2*(32.*U1*T3+32.*U1*S13)
  90      1 +CUXCU*CA3*(-64.*ZM2*U1**2*U3-32.*U1**2*ZM4)
  91       WZTU=WZTU
  92      1 +CUXCU*CA3*T*U*(-32.*U1*U3)
  93      1 +CUXCU*CA3*U*(32.*WM2*U1*U3-32.*ZM2*U1*T3+32.*ZM2*U1*U3-32.*ZM2*
  94      1 U1*S13+32.*ZM2*U1**2+32.*T1*U1*U3+32.*U1**2*U3)
  95      1 +CUXCU*CA3*U**2*(32.*U1*T3+32.*U1*S13)
  96      1 +EPF(P1,P2,P3,Q1)*CTXCU*CV3*(-32.*ZM2*T1-16.*WZM2)
  97      1 +EPF(P1,P2,P3,Q1)*CTXCU*CV3*T*(16.*ZM2)
  98       WZTU=WZTU
  99      1 +EPF(P1,P2,P3,Q3)*CTXCU*CV3*(-48.*WM2*T1-16.*WM2*U1-32.*WM2*S13-
 100      1 32.*T1*U1-64.*T1*S13-32.*T1**2-16.*WM4)
 101      1 +EPF(P1,P2,P3,Q3)*CTXCU*CV3*T*(16.*WM2+16.*T1+16.*U1+32.*S13)
 102      1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*(-32.*WM2*T3-32.*WM2*U3-64.*ZM2*T1-6
 103      1 4.*T1*T3-64.*T1*U3-32.*WZM2)
 104      1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T*(16.*WM2+32.*ZM2+32.*T1+32.*T3+32.
 105      1 *U3)
 106       WZTU=WZTU
 107      1 -16.*EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T*U
 108      1 -16.*EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T**2
 109      1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*U*(16.*WM2+32.*T1)
 110      1 +EPF(P1,P3,Q1,Q3)*CTXCU*CV3*(32.*WM2*T3+32.*ZM2*T1+64.*T1*T3+16.
 111      1 *WZM2)
 112      1 +EPF(P1,P3,Q1,Q3)*CTXCU*CV3*T*(-16.*WM2-16.*ZM2-32.*T1-32.*T3)
 113       WZTU=WZTU
 114      1 +16.*EPF(P1,P3,Q1,Q3)*CTXCU*CV3*T**2
 115      1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*(-32.*WM2*U3-32.*ZM2*T1-64.*T1*U3-16
 116      1 .*WZM2)
 117      1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*T*(16.*ZM2+32.*U3)
 118      1 -16.*EPF(P2,P3,Q1,Q3)*CTXCU*CV3*T*U
 119      1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*U*(16.*WM2+32.*T1)+0.
 120       RETURN
 121       END