HERWIG/src/hwupcm.f

   1
   2 CDECK  ID>, HWUPCM.
   3
   4 *CMZ :-        -26/04/91  11.11.56  by  Bryan Webber
   5
   6 *-- Author :    Bryan Webber
   7
   8 C-----------------------------------------------------------------------
   9
  10       FUNCTION HWUPCM(EM0,EM1,EM2)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     C.M. MOMENTUM FOR DECAY MASSES EM0 -> EM1 + EM2
  15
  16 C     SET TO -1 BELOW THRESHOLD
  17
  18 C-----------------------------------------------------------------------
  19
  20       DOUBLE PRECISION HWUPCM,EM0,EM1,EM2,EMS,EMD
  21
  22       EMS=ABS(EM1+EM2)
  23
  24       EMD=ABS(EM1-EM2)
  25
  26       IF (EM0.LT.EMS.OR.EM0.LT.EMD) THEN
  27
  28         HWUPCM=-1.
  29
  30       ELSEIF (EM0.EQ.EMS.OR.EM0.EQ.EMD) THEN
  31
  32         HWUPCM=0.
  33
  34       ELSE
  35
  36         HWUPCM=SQRT((EM0+EMD)*(EM0-EMD)*
  37
  38      &              (EM0+EMS)*(EM0-EMS))*.5/EM0
  39
  40       ENDIF
  41
  42       END