HERWIG/src/hwumbw.f

   1
   2 CDECK  ID>, HWUMBW.
   3
   4 *CMZ :-        -21/02/98  11.11.56  by  Bryan Webber
   5
   6 *-- Author :    Bryan Webber
   7
   8 C-----------------------------------------------------------------------
   9
  10       FUNCTION HWUMBW(ID)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     CHOOSES MASS ACCORDING TO BREIT-WIGNER DISTRIBUTION
  15
  16 C-----------------------------------------------------------------------
  17
  18       INCLUDE 'HERWIG61.INC'
  19
  20       DOUBLE PRECISION HWUMBW,HWR,WMX,TAU,T
  21
  22       INTEGER ID
  23
  24 C--WMX IS MAX NUMBER OF WIDTHS FROM NOMINAL MASS
  25
  26       DATA WMX/10D0/
  27
  28       SAVE WMX
  29
  30       HWUMBW=RMASS(ID)
  31
  32       TAU=RLTIM(ID)
  33
  34       IF (TAU.EQ.ZERO.OR.TAU.GT.1D-18) RETURN
  35
  36  1    T=HALF*TAN(PIFAC*(HWR()-HALF))
  37
  38       IF (ABS(T).GT.WMX) GO TO 1
  39
  40       HWUMBW=HWUMBW+HBAR*T/TAU
  41
  42       IF(HWUMBW.LT.ZERO) GOTO 1
  43
  44       END