ISAJET/isasusy/sshww1.F

   1 #include "isajet/pilot.h"
   2       DOUBLE PRECISION FUNCTION SSHWW1(XX)
   3 C-----------------------------------------------------------------------
   4 C     SSHWW: hi -> w + w*
   5 C     Bisset's FUNWW
   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,MW,DELTAW,MH
  15 C
  16       MW=AMW
  17       DELTAW=GAMW
  18       MH=TMP(1)
  19 C
  20       EPLN=MW/MH
  21       PROP=(1.D0-XX)**2
  22       PROP=PROP+(EPLN**2)*DELTAW**2/MH**2
  23       TEMP=XX**2-12*XX*EPLN**2
  24       TEMP=TEMP+8*EPLN**2+12*EPLN**4
  25       AAA=XX**2-4*EPLN**2
  26       IF(AAA.LT.0) THEN
  27          AAA=0
  28       ENDIF
  29       FN=TEMP*SQRT(AAA)/PROP
  30       SSHWW1=FN
  31       RETURN
  32       END