1 #include "isajet/pilot.h"
2 SUBROUTINE LBOOST(PREST,N,P1,P2)
4 C BOOST 4-VECTORS P1 TO PREST REST FRAME
5 C PUT RESULTING 4-VECTORS IN P2
7 #include "isajet/itapes.inc"
8 DIMENSION PREST(4),P1(4,N),P2(4,N)
10 WCN=SQRT(PREST(4)**2-PREST(1)**2-PREST(2)**2-PREST(3)**2)
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))
16 2 P2(K,I)=P1(K,I)-FACT*PREST(K)