HIJING/hijing1_36/hiboost.F

   1 * $Id$
   2 C
   3 C
   4 C
   5 C
   6         SUBROUTINE HIBOOST
   7         IMPLICIT DOUBLE PRECISION(D)
   8 #include "lujets_hijing.inc"
   9 #include "ludat1_hijing.inc"
  10 #include "hiparnt.inc"
  11         SAVE
  12         DO 100 I=1,N
  13            DBETA=P(I,3)/P(I,4)
  14            IF(ABS(DBETA).GE.1.D0) THEN
  15               DB=HINT1(2)
  16               IF(DB.GT.0.99999999D0) THEN
  17 C               ********Rescale boost vector if too close to unity.
  18                  WRITE(6,*) '(HIBOOT:) boost vector too large'
  19                  DB=0.99999999D0
  20               ENDIF
  21               DGA=1D0/SQRT(1D0-DB**2)
  22               DP3=P(I,3)
  23               DP4=P(I,4)
  24               P(I,3)=(DP3+DB*DP4)*DGA
  25               P(I,4)=(DP4+DB*DP3)*DGA
  26               GO TO 100
  27            ENDIF
  28            Y=0.5*DLOG((1.D0+DBETA)/(1.D0-DBETA))
  29            AMT=SQRT(P(I,1)**2+P(I,2)**2+P(I,5)**2)
  30            P(I,3)=AMT*SINH(Y+HINT1(3))
  31            P(I,4)=AMT*COSH(Y+HINT1(3))
  32 100     CONTINUE
  33         RETURN
  34         END