[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / hiboost.F

* $Id$
C
C
C
C
	SUBROUTINE HIBOOST
      	IMPLICIT DOUBLE PRECISION(D)  
#include "lujets_hijing.inc"
#include "ludat1_hijing.inc"
#include "hiparnt.inc"
	SAVE
	DO 100 I=1,N
	   DBETA=P(I,3)/P(I,4)
	   IF(ABS(DBETA).GE.1.D0) THEN
	      DB=HINT1(2)
	      IF(DB.GT.0.99999999D0) THEN 
C		********Rescale boost vector if too close to unity. 
		 WRITE(6,*) '(HIBOOT:) boost vector too large' 
		 DB=0.99999999D0
	      ENDIF 
	      DGA=1D0/SQRT(1D0-DB**2)
	      DP3=P(I,3)
	      DP4=P(I,4)
	      P(I,3)=(DP3+DB*DP4)*DGA  
	      P(I,4)=(DP4+DB*DP3)*DGA  
	      GO TO 100
	   ENDIF
	   Y=0.5*DLOG((1.D0+DBETA)/(1.D0-DBETA))
	   AMT=SQRT(P(I,1)**2+P(I,2)**2+P(I,5)**2)
	   P(I,3)=AMT*SINH(Y+HINT1(3))
	   P(I,4)=AMT*COSH(Y+HINT1(3))
100	CONTINUE
	RETURN
	END
Commit	Line	Data
e74335a4	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+DBDP4)DGA
	25	P(I,4)=(DP4+DBDP3)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