[u/mrichter/AliRoot.git] / PHOS / shaker / shtolu.f

*CMZ :          17/07/98  15.49.05  by  Federico Carminati
*-- Author :
      SUBROUTINE SHTOLU(KF,Y,PT,PHI,W)
c	================================

c	Load Particle in /LUJETS/

*KEEP,LUJETS.
      COMMON /LUJETS/ N,K(200000,5),P(200000,5),V(200000,5)
      SAVE /LUJETS/
*KEEP,SHWATE.
      COMMON /SHWATE/ WEI(200000)

*KEND.
      N=N+1
      K(N,1) = 1
      K(N,2) = KF
      K(N,3) = 0
      K(N,4) = 0
      K(N,5) = 0

      AM = ULMASS(KF)
      AMT = SQRT(AM**2+PT**2)
      P(N,1) = PT*COS(PHI)
      P(N,2) = PT*SIN(PHI)
      P(N,3) = AMT*SINH(Y)
      P(N,4) = AMT*COSH(Y)
      P(N,5) = AM
      DO 50 JJ=1,5
        V(N,JJ) = 0.
50	CONTINUE
      WEI(N) = W

      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*CMZ : 17/07/98 15.49.05 by Federico Carminati
	2	*-- Author :
	3	SUBROUTINE SHTOLU(KF,Y,PT,PHI,W)
	4	c ================================
	5
	6	c Load Particle in /LUJETS/
	7
	8	*KEEP,LUJETS.
	9	COMMON /LUJETS/ N,K(200000,5),P(200000,5),V(200000,5)
	10	SAVE /LUJETS/
	11	*KEEP,SHWATE.
	12	COMMON /SHWATE/ WEI(200000)
	13
	14	*KEND.
	15	N=N+1
	16	K(N,1) = 1
	17	K(N,2) = KF
	18	K(N,3) = 0
	19	K(N,4) = 0
	20	K(N,5) = 0
	21
	22	AM = ULMASS(KF)
	23	AMT = SQRT(AM2+PT2)
	24	P(N,1) = PT*COS(PHI)
	25	P(N,2) = PT*SIN(PHI)
	26	P(N,3) = AMT*SINH(Y)
	27	P(N,4) = AMT*COSH(Y)
	28	P(N,5) = AM
	29	DO 50 JJ=1,5
	30	V(N,JJ) = 0.
	31	50 CONTINUE
	32	WEI(N) = W
	33
	34	RETURN
	35	END