HERWIG/src/hwurap.f

   1
   2 CDECK  ID>, HWURAP.
   3
   4 *CMZ :-        -26/04/91  11.11.56  by  Bryan Webber
   5
   6 *-- Author :    Bryan Webber
   7
   8 C-----------------------------------------------------------------------
   9
  10       FUNCTION HWURAP(P)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     LONGITUDINAL RAPIDITY (SET TO +/-1000 IF TOO LARGE)
  15
  16 C-----------------------------------------------------------------------
  17
  18       DOUBLE PRECISION HWURAP,EMT2,P(5),ZERO
  19
  20       PARAMETER (ZERO=0.D0)
  21
  22       EMT2=P(1)**2+P(2)**2+P(5)**2
  23
  24       IF (P(3).GT.ZERO) THEN
  25
  26         IF (EMT2.EQ.ZERO) THEN
  27
  28           HWURAP=1000.
  29
  30         ELSE
  31
  32           HWURAP= 0.5*LOG((P(3)+P(4))**2/EMT2)
  33
  34         ENDIF
  35
  36       ELSEIF (P(3).LT.ZERO) THEN
  37
  38         IF (EMT2.EQ.ZERO) THEN
  39
  40           HWURAP=-1000.
  41
  42         ELSE
  43
  44           HWURAP=-0.5*LOG((P(3)-P(4))**2/EMT2)
  45
  46         ENDIF
  47
  48       ELSE
  49
  50           HWURAP=0.
  51
  52       ENDIF
  53
  54       END