ISAJET/isasusy/sszzf4.F

   1 #include "isajet/pilot.h"
   2       REAL FUNCTION SSZZF4(E)
   3 C-----------------------------------------------------------------------
   4 C          Z_I -> Z_J +B +BBAR VIA HIGGS
   5 C-----------------------------------------------------------------------
   6 #if defined(CERNLIB_IMPNONE)
   7       IMPLICIT NONE
   8 #endif
   9 #include "isajet/sssm.inc"
  10 #include "isajet/sspar.inc"
  11 #include "isajet/sstmp.inc"
  12 C
  13       REAL E
  14       DOUBLE PRECISION MZ1,MZ2,SN,EE,MH1,MH2,P,MF,BF,TM
  15 C
  16       MZ2=TMP(1)
  17       MZ1=TMP(2)
  18       MH1=TMP(3)
  19       SN=TMP(4)
  20       MH2=TMP(5)
  21       MF=TMP(6)
  22       EE=E
  23 C
  24       P=SQRT(MAX(0.D0,EE**2-MZ1**2))
  25       TM=1.D0-4*MF**2/(MZ2**2+MZ1**2-2*E*MZ2)
  26       BF=DSQRT(MAX(0.D0,TM))
  27       SSZZF4=P*BF*(MZ2**2+MZ1**2-2*MZ2*EE-2*MF**2)*
  28      $(MZ2*EE+SN*MZ2*MZ1)/(MZ2**2+MZ1**2-2*MZ2*EE-MH1**2)/
  29      $(MZ2**2+MZ1**2-2*MZ2*EE-MH2**2)
  30       RETURN
  31       END