ISAJET/code/lboost.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE LBOOST(PREST,N,P1,P2)
   3 C
   4 C            BOOST 4-VECTORS P1 TO PREST REST FRAME
   5 C            PUT RESULTING 4-VECTORS IN P2
   6 C
   7 #include "isajet/itapes.inc"
   8       DIMENSION PREST(4),P1(4,N),P2(4,N)
   9       DO 1 I=1,N
  10       WCN=SQRT(PREST(4)**2-PREST(1)**2-PREST(2)**2-PREST(3)**2)
  11       II=(I-1)*4
  12       P2(4,I)=(P1(4,I)*PREST(4)-P1(1,I)*PREST(1)-P1(2,I)*PREST(2)
  13      1         -P1(3,I)*PREST(3))/WCN
  14       FACT=(P2(4,I)+P1(4,I))/(WCN+PREST(4))
  15       DO 2 K=1,3
  16     2 P2(K,I)=P1(K,I)-FACT*PREST(K)
  17     1 CONTINUE
  18       RETURN
  19       END