ISAJET/isasusy/sshww2.F

   1 #include "isajet/pilot.h"
   2       DOUBLE PRECISION FUNCTION SSHWW2(XX)
   3 C-----------------------------------------------------------------------
   4 C     SSHWW: hi -> z + z*
   5 C     Bisset's FUNZZ
   6 C-----------------------------------------------------------------------
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 #include "isajet/sspar.inc"
  11 #include "isajet/sssm.inc"
  12 #include "isajet/sstmp.inc"
  13 C
  14       DOUBLE PRECISION XX,EPLN,PROP,TEMP,FN,AAA,MZ,DELTAZ,MH
  15 C
  16 C
  17       MZ=AMZ
  18       DELTAZ=GAMZ
  19       MH=TMP(1)
  20 C
  21       EPLN=MZ/MH
  22       PROP=(1.D0-XX)**2
  23       PROP=PROP+(EPLN**2)*DELTAZ**2/MH**2
  24       TEMP=XX**2-12*XX*EPLN**2
  25       TEMP=TEMP+8*EPLN**2+12*EPLN**4
  26       AAA=XX**2-4*EPLN**2
  27       IF(AAA.LT.0) THEN
  28         AAA=0
  29       ENDIF
  30       FN=TEMP*SQRT(AAA)/PROP
  31       SSHWW2=FN
  32       RETURN
  33       END