ISAJET/code/wwst.F

   1 #include "isajet/pilot.h"
   2       FUNCTION WWST(T,U,T1,U1,T3,U3,P1,P2)
   3 C          DECAY DISTRIBUTION FOR W+ W- PAIRS FROM SCHOONSCHIP(1980).
   4 C          INTERFERENCE OF T AND S 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 WWST
  10       DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
  11       DOUBLE PRECISION CVACQ,EPF
  12 #endif
  13       CVACQ=(CV+CA)*CQ
  14       WWST=
  15      1 +CVACQ*(-256.*T1*U1*T3*U3+256.*T1**2*T3**2)
  16      1 +CVACQ*T*(256.*T1*U1*T3+256.*T1*T3*S13-256.*T1*T3**2+128.*U1*T3*
  17      1 U3-128.*U1*U3*S13-128.*U1**2*U3)
  18      1 +CVACQ*T*U*(-256.*U1*T3-128.*U1*S13-128.*T3*S13-128.*S13**2)
  19      1 +CVACQ*T*WM2*(384.*U1*T3-128.*U1*S13-128.*U1**2+256.*T3*S13-256.
  20      1 *T3**2)
  21      1 +CVACQ*T**2*(-256.*U1*T3-256.*T3*S13)
  22      1 +CVACQ*U*(128.*T1*U1*T3-128.*T1*T3*S13-128.*T1*T3**2)
  23       WWST=WWST
  24      1 +CVACQ*U*WM2*(128.*U1*T3-128.*T3*S13-128.*T3**2)
  25      1 +CVACQ*WM2*(-256.*T1*U1*T3+512.*T1*T3**2-256.*U1*T3*U3)
  26      1 +CVACQ*WM2**2*(256.*U1*S13+128.*U1**2+256.*T3*S13+384.*T3**2+128
  27      1 .*S13**2)
  28      1 +EPF(P1,P2,P3,Q1)*CVACQ*(128.*T3*U3+128.*T3*S13+64.*T3**2+128.*U
  29      1 3*S13+64.*U3**2)
  30      1 +EPF(P1,P2,P3,Q1)*CVACQ*T*(-32.*T3-32.*U3-64.*S13)
  31      1 -32.*EPF(P1,P2,P3,Q1)*CVACQ*T*WM2
  32       WWST=WWST
  33      1 +EPF(P1,P2,P3,Q1)*CVACQ*U*(-32.*T3-32.*U3-64.*S13)
  34      1 -32.*EPF(P1,P2,P3,Q1)*CVACQ*U*WM2
  35      1 +EPF(P1,P2,P3,Q1)*CVACQ*WM2*(128.*T3+128.*U3+128.*S13)
  36      1 +64.*EPF(P1,P2,P3,Q1)*CVACQ*WM2**2
  37      1 -32.*EPF(P1,P2,P3,Q3)*CVACQ*T*WM2
  38      1 -32.*EPF(P1,P2,P3,Q3)*CVACQ*U*WM2
  39      1 +EPF(P1,P2,P3,Q3)*CVACQ*WM2*(64.*T3+64.*U3)
  40      1 +64.*EPF(P1,P2,P3,Q3)*CVACQ*WM2**2
  41       WWST=WWST
  42      1 +EPF(P1,P3,Q1,Q3)*CVACQ*(128.*U1*T3+128.*U1*U3)
  43      1 +EPF(P1,P3,Q1,Q3)*CVACQ*T*(-64.*U1)
  44      1 +32.*EPF(P1,P3,Q1,Q3)*CVACQ*T*U
  45      1 -32.*EPF(P1,P3,Q1,Q3)*CVACQ*T*WM2
  46      1 +EPF(P1,P3,Q1,Q3)*CVACQ*U*(-64.*U1-64.*T3-64.*U3)
  47      1 -96.*EPF(P1,P3,Q1,Q3)*CVACQ*U*WM2
  48      1 +32.*EPF(P1,P3,Q1,Q3)*CVACQ*U**2
  49      1 +EPF(P1,P3,Q1,Q3)*CVACQ*WM2*(128.*U1+64.*T3+64.*U3)
  50       WWST=WWST
  51      1 +64.*EPF(P1,P3,Q1,Q3)*CVACQ*WM2**2
  52      1 +EPF(P2,P3,Q1,Q3)*CVACQ*(-128.*T1*T3-128.*T1*U3)
  53      1 +EPF(P2,P3,Q1,Q3)*CVACQ*T*(64.*T1+64.*T3+64.*U3)
  54      1 -32.*EPF(P2,P3,Q1,Q3)*CVACQ*T*U
  55      1 +96.*EPF(P2,P3,Q1,Q3)*CVACQ*T*WM2
  56      1 -32.*EPF(P2,P3,Q1,Q3)*CVACQ*T**2
  57      1 +EPF(P2,P3,Q1,Q3)*CVACQ*U*(64.*T1)
  58      1 +32.*EPF(P2,P3,Q1,Q3)*CVACQ*U*WM2
  59      1 +EPF(P2,P3,Q1,Q3)*CVACQ*WM2*(-128.*T1-64.*T3-64.*U3)
  60      1 -64.*EPF(P2,P3,Q1,Q3)*CVACQ*WM2**2
  61       WWST=WWST/T
  62       WWST=2.*WWST
  63       RETURN
  64       END