HIJING/hijing1_36/atrobo.F

   1 * $Id$
   2
   3
   4 C
   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)
   9 #include "lujets_hijing.inc"
  10         DIMENSION ROT(3,3),PV(3)
  11         DOUBLE PRECISION DP(4),DBEX,DBEY,DBEZ,DGA,DGA2,DBEP,DGABEP
  12         SAVE
  13         IERROR=0
  14
  15         IF(IMIN.LE.0 .OR. IMAX.GT.N .OR. IMIN.GT.IMAX) RETURN
  16
  17         IF(THE**2+PHI**2.GT.1E-20) THEN
  18 C...ROTATE (TYPICALLY FROM Z AXIS TO DIRECTION THETA,PHI)
  19            ROT(1,1)=COS(THE)*COS(PHI)
  20            ROT(1,2)=-SIN(PHI)
  21            ROT(1,3)=SIN(THE)*COS(PHI)
  22            ROT(2,1)=COS(THE)*SIN(PHI)
  23            ROT(2,2)=COS(PHI)
  24            ROT(2,3)=SIN(THE)*SIN(PHI)
  25            ROT(3,1)=-SIN(THE)
  26            ROT(3,2)=0.
  27            ROT(3,3)=COS(THE)
  28            DO 120 I=IMIN,IMAX
  29 C**************    IF(MOD(K(I,1)/10000,10).GE.6) GOTO 120
  30               DO 100 J=1,3
  31  100             PV(J)=P(I,J)
  32                  DO 110 J=1,3
  33  110                P(I,J)=ROT(J,1)*PV(1)+ROT(J,2)*PV(2)
  34      &                     +ROT(J,3)*PV(3)
  35  120             CONTINUE
  36         ENDIF
  37
  38         IF(BEX**2+BEY**2+BEZ**2.GT.1E-20) THEN
  39 C...LORENTZ BOOST (TYPICALLY FROM REST TO MOMENTUM/ENERGY=BETA)
  40                 DBEX=BEX
  41                 DBEY=BEY
  42                 DBEZ=BEZ
  43                 DGA2=1D0-DBEX**2-DBEY**2-DBEZ**2
  44                 IF(DGA2.LE.0D0) THEN
  45                         IERROR=1
  46                         RETURN
  47                 ENDIF
  48                 DGA=1D0/DSQRT(DGA2)
  49                 DO 140 I=IMIN,IMAX
  50 C*************     IF(MOD(K(I,1)/10000,10).GE.6) GOTO 140
  51                    DO 130 J=1,4
  52 130                DP(J)=P(I,J)
  53                    DBEP=DBEX*DP(1)+DBEY*DP(2)+DBEZ*DP(3)
  54                    DGABEP=DGA*(DGA*DBEP/(1D0+DGA)+DP(4))
  55                    P(I,1)=DP(1)+DGABEP*DBEX
  56                    P(I,2)=DP(2)+DGABEP*DBEY
  57                    P(I,3)=DP(3)+DGABEP*DBEZ
  58                    P(I,4)=DGA*(DP(4)+DBEP)
  59 140             CONTINUE
  60         ENDIF
  61
  62         RETURN
  63         END