ISAJET/isasusy/sszzf2.F

   1 #include "isajet/pilot.h"
   2       REAL FUNCTION SSZZF2(E)
   3 C-----------------------------------------------------------------------
   4 C          SSWZBF: ziss -> zjss f fbar
   5 C          Baer's Z2ZFUN
   6 C-----------------------------------------------------------------------
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 #include "isajet/sssm.inc"
  11 #include "isajet/sspar.inc"
  12 #include "isajet/sstmp.inc"
  13 C
  14       REAL E
  15       DOUBLE PRECISION MZ1,MZ2,SN,M,PI,F1,F2,ZWID,MF,BF
  16       DATA PI/3.14159265D0/
  17 C
  18       MZ2=TMP(1)
  19       MZ1=TMP(2)
  20       SN=TMP(4)
  21       MF=TMP(6)
  22       M=AMZ
  23       ZWID=GAMZ
  24       BF=DSQRT(MAX(0.D0,1.D0-4*MF**2/(MZ2**2+MZ1**2-2*E*MZ2)))
  25 C
  26       F1=SQRT(MAX(0.D0,E**2-MZ1**2))/
  27      $((MZ1**2+MZ2**2-M**2-2*MZ2*E)**2+ZWID**2*M**2)
  28       F2=E*(MZ1**2+MZ2**2+2*SN*MZ1*MZ2)-MZ2*(E**2+MZ1**2)-
  29      $BF*MZ2*(E**2-MZ1**2)/3.D0-SN*MZ1*(MZ1**2+MZ2**2-2*MF**2)
  30       SSZZF2=MZ2*BF*F1*F2/2.D0/PI**3
  31       RETURN
  32       END