ISAJET/code/zzall.F

   1 #include "isajet/pilot.h"
   2       FUNCTION ZZALL(T,U,T1,U1,T3,U3,P1,P2)
   3 C          DECAY DISTRIBUTION FOR Z0 Z0 PAIRS FROM SCHOONSCHIP(1980).
   4 C          INCLUDES TT, TU, AND UU 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 ZZALL
  10       DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
  11       DOUBLE PRECISION CVAVA,DVAVA,CVVVV,CVVAA,TT,TU,UU,EPF
  12 #endif
  13       TTUU=(T*U)**2
  14       CVAVA=CV*CA*(CV1*CA3+CA1*CV3)/TTUU
  15       DVAVA=CV*CA*(CV1*CA3-CA1*CV3)/TTUU
  16       CVVVV=(CV**2+CA**2)*CV1*CV3/TTUU
  17       CVVAA=(CV**2+CA**2)*CA1*CA3/TTUU
  18       TT=T**2
  19       TU=T*U
  20       UU=U**2
  21       ZZALL=
  22      1 +CVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3+16.*ZM2*T1*T3*U3-
  23      1 16.*ZM2*U1*T3*U3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1**2+16.*ZM2**2*U1
  24      1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2-16.*ZM2**
  25      1 2*U3*S13-8.*ZM2**2*U3**2)
  26      1 +CVAVA*T*U**2*(16.*ZM2*T1*S13-16.*ZM2*T3*S13)
  27      1 +CVAVA*T**2*(16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-8.*ZM2**2*U1**2+8
  28      1 .*ZM2**2*U3**2)
  29      1 +CVAVA*T**2*U*(-16.*ZM2*U1*S13+16.*ZM2*U3*S13)
  30      1 +CVAVA*U**2*(-16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+8.*ZM2**2*T1**2-
  31      1 8.*ZM2**2*T3**2)
  32       ZZALL=ZZALL
  33      1 +DVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3+
  34      1 16.*ZM2*U1*T3*U3-32.*ZM2**2*T1*T3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1
  35      1 **2+32.*ZM2**2*U1*U3+16.*ZM2**2*U1*S13+8.*ZM2**2*U1**2-16.*ZM2**
  36      1 2*T3*S13-8.*ZM2**2*T3**2+16.*ZM2**2*U3*S13+8.*ZM2**2*U3**2-32.*Z
  37      1 M2**3*T1+32.*ZM2**3*U1-32.*ZM2**3*T3+32.*ZM2**3*U3)
  38      1 +DVAVA*T*U**2*(16.*ZM2*T1*S13+16.*ZM2*T3*S13-16.*ZM2**2*U1-16.*Z
  39      1 M2**2*U3-32.*ZM2**3)
  40      1 +DVAVA*T*U**3*(8.*ZM2**2)
  41       ZZALL=ZZALL
  42      1 +DVAVA*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-32.*ZM2**2*U1*U3
  43      1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-16.*ZM2**3*U1-16.*ZM2**3*U3-8.*
  44      1 ZM2**4)
  45      1 +DVAVA*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13+16.*ZM2**2*T1+16.*
  46      1 ZM2**2*T3+32.*ZM2**3)
  47      1 +DVAVA*T**3*U*(-8.*ZM2**2)
  48      1 +DVAVA*U**2*(16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+32.*ZM2**2*T1*T3+
  49      1 8.*ZM2**2*T1**2+8.*ZM2**2*T3**2+16.*ZM2**3*T1+16.*ZM2**3*T3+8.*Z
  50      1 M2**4)
  51       ZZALL=ZZALL
  52      1 +CVVVV*T*U*(-16.*ZM2*T1*U1*T3-16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3-
  53      1 16.*ZM2*U1*T3*U3+16.*ZM2**2*T1*S13+8.*ZM2**2*T1**2+16.*ZM2**2*U1
  54      1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2+16.*ZM2**
  55      1 2*U3*S13+8.*ZM2**2*U3**2+16.*ZM2**2*S13**2+16.*ZM2**3*T1+16.*ZM2
  56      1 **3*U1+16.*ZM2**3*T3+16.*ZM2**3*U3+32.*ZM2**3*S13+32.*ZM2**4-32.
  57      1 *T1*U1*T3*U3)
  58      1 +CVVVV*T*U**2*(-16.*ZM2*T1*S13-16.*ZM2*T3*S13-8.*ZM2**2*U1-8.*ZM
  59      1 2**2*U3-16.*ZM2**2*S13-16.*ZM2**3-32.*T1*T3*S13)
  60       ZZALL=ZZALL
  61      1 +CVVVV*T*U**3*(4.*ZM2**2)
  62      1 +CVVVV*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-16.*ZM2**2*U1*U3
  63      1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-8.*ZM2**3*U1-8.*ZM2**3*U3-4.*ZM
  64      1 2**4-16.*U1**2*U3**2)
  65      1 +CVVVV*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13-8.*ZM2**2*T1-8.*ZM
  66      1 2**2*T3-16.*ZM2**2*S13-16.*ZM2**3-32.*U1*U3*S13)
  67      1 +CVVVV*T**2*U**2*(-16.*S13**2)
  68      1 +CVVVV*T**3*U*(4.*ZM2**2)
  69       ZZALL=ZZALL
  70      1 +CVVVV*U**2*(-16.*ZM2*T1*T3**2-16.*ZM2*T1**2*T3-16.*ZM2**2*T1*T3
  71      1 -8.*ZM2**2*T1**2-8.*ZM2**2*T3**2-8.*ZM2**3*T1-8.*ZM2**3*T3-4.*ZM
  72      1 2**4-16.*T1**2*T3**2)
  73      1 +CVVAA*T*U*(-16.*ZM2**3*T1-16.*ZM2**3*U1-16.*ZM2**3*T3-16.*ZM2**
  74      1 3*U3-32.*ZM2**3*S13-32.*ZM2**4)
  75      1 +CVVAA*T*U**2*(8.*ZM2**2*U1+8.*ZM2**2*U3+16.*ZM2**2*S13+16.*ZM2*
  76      1 *3)
  77      1 +CVVAA*T*U**3*(-4.*ZM2**2)
  78       ZZALL=ZZALL
  79      1 +CVVAA*T**2*(16.*ZM2**2*U1*U3+8.*ZM2**3*U1+8.*ZM2**3*U3+4.*ZM2**
  80      1 4)
  81      1 +CVVAA*T**2*U*(8.*ZM2**2*T1+8.*ZM2**2*T3+16.*ZM2**2*S13+16.*ZM2*
  82      1 *3)
  83      1 +CVVAA*T**3*U*(-4.*ZM2**2)
  84      1 +CVVAA*U**2*(16.*ZM2**2*T1*T3+8.*ZM2**3*T1+8.*ZM2**3*T3+4.*ZM2**
  85      1 4)
  86 C          THE EPF TERMS FROM SCHOONSCHIP DID NOT CONTAIN CVAVA, ETC.,
  87 C          BUT DID CONTAIN AN EXPLICIT T*U. THIS WAS REPLACED BY HAND
  88 C          BY 1./TU.
  89       ZZALL=ZZALL
  90      1 +EPF(P1,P2,P3,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*S13)
  91      1 +EPF(P1,P2,Q1,Q3)/TU*(-32.*ZM2*CV*CA*CV1*CV3*T1-32.*CV*CA*CV1*C
  92      1 V3*T1*T3-32.*CV*CA*CV1*CV3*T1*U3)
  93      1 +EPF(P1,P3,Q1,Q3)/TU*(32.*CV*CA*CV1*CV3*T1*T3)
  94      1 +EPF(P2,P3,Q1,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*U3)+0.
  95       ZZALL=2.*ZZALL
  96       RETURN
  97       END