HERWIG/src/hwulf4.f

   1
   2 CDECK  ID>, HWULF4.
   3
   4 *CMZ :-        -05/11/95  19.33.42  by  Mike Seymour
   5
   6 *-- Author :    Adapted by Bryan Webber
   7
   8 C-----------------------------------------------------------------------
   9
  10       SUBROUTINE HWULF4(PS,PI,PF)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     TRANSFORMS PI (GIVEN IN LAB) INTO PF (IN REST FRAME OF PS)
  15
  16 C     N.B. P(1,2,3,4) = (PX,PY,PZ,E); PS(5)=M
  17
  18 C-----------------------------------------------------------------------
  19
  20       DOUBLE PRECISION PF4,FN,PS(5),PI(4),PF(4)
  21
  22       IF (PS(4).EQ.PS(5)) THEN
  23
  24         PF(1)= PI(1)
  25
  26         PF(2)= PI(2)
  27
  28         PF(3)= PI(3)
  29
  30         PF(4)= PI(4)
  31
  32       ELSE
  33
  34         PF4  = (PI(4)*PS(4)-PI(3)*PS(3)
  35
  36      &         -PI(2)*PS(2)-PI(1)*PS(1))/PS(5)
  37
  38         FN   = (PF4+PI(4)) / (PS(4)+PS(5))
  39
  40         PF(1)= PI(1) - FN*PS(1)
  41
  42         PF(2)= PI(2) - FN*PS(2)
  43
  44         PF(3)= PI(3) - FN*PS(3)
  45
  46         PF(4)= PF4
  47
  48       END IF
  49
  50       END