ISAJET/isasusy/sswzf1.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSWZF1(EE)
   3 C-----------------------------------------------------------------------
   4 C          SSWZBF: wiss -> zjss f fbar
   5 C          Baer's WIWFUN
   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 EE
  15       DOUBLE PRECISION C1,C2,MWI,MZI,MW,PROP,T1,T2,E
  16 C
  17       E=EE
  18       C1=TMP(1)
  19       C2=TMP(2)
  20       MWI=TMP(3)
  21       MZI=TMP(4)
  22       MW=AMW
  23 C
  24       PROP=(MWI**2+MZI**2-2*MWI*E-MW**2)**2
  25       T1=C1*(3*(MWI**2+MZI**2)*MWI*E-4*MWI**2*E*E-2*MWI**2*MZI**2)
  26       T2=-3*C2*MWI*MZI*(MWI**2+MZI**2-2*MWI*E)
  27       SSWZF1=SQRT(MAX(0.D0,E*E-MZI**2))/PROP*(T1+T2)
  28       RETURN
  29       END