5 C*******************************************************************
6 C make boost and rotation to entries from IMIN to IMAX
7 C*******************************************************************
8 SUBROUTINE ATROBO(THE,PHI,BEX,BEY,BEZ,IMIN,IMAX,IERROR)
10 #include "lujets_hijing.inc"
11 DIMENSION ROT(3,3),PV(3)
12 DOUBLE PRECISION DP(4),DBEX,DBEY,DBEZ,DGA,DGA2,DBEP,DGABEP
16 IF(IMIN.LE.0 .OR. IMAX.GT.N .OR. IMIN.GT.IMAX) RETURN
18 IF(THE**2+PHI**2.GT.1E-20) THEN
19 C...ROTATE (TYPICALLY FROM Z AXIS TO DIRECTION THETA,PHI)
20 ROT(1,1)=COS(THE)*COS(PHI)
22 ROT(1,3)=SIN(THE)*COS(PHI)
23 ROT(2,1)=COS(THE)*SIN(PHI)
25 ROT(2,3)=SIN(THE)*SIN(PHI)
30 C************** IF(MOD(K(I,1)/10000,10).GE.6) GOTO 120
34 110 P(I,J)=ROT(J,1)*PV(1)+ROT(J,2)*PV(2)
39 IF(BEX**2+BEY**2+BEZ**2.GT.1E-20) THEN
40 C...LORENTZ BOOST (TYPICALLY FROM REST TO MOMENTUM/ENERGY=BETA)
44 DGA2=1D0-DBEX**2-DBEY**2-DBEZ**2
51 C************* IF(MOD(K(I,1)/10000,10).GE.6) GOTO 140
54 DBEP=DBEX*DP(1)+DBEY*DP(2)+DBEZ*DP(3)
55 DGABEP=DGA*(DGA*DBEP/(1D0+DGA)+DP(4))
56 P(I,1)=DP(1)+DGABEP*DBEX
57 P(I,2)=DP(2)+DGABEP*DBEY
58 P(I,3)=DP(3)+DGABEP*DBEZ
59 P(I,4)=DGA*(DP(4)+DBEP)