[u/mrichter/AliRoot.git] / HIJING / hipyset1_35 / luptdi_hijing.F

* $Id$
    
C*********************************************************************  
    
      SUBROUTINE LUPTDI_HIJING(KFL,PX,PY)  
    
C...Purpose: to generate transverse momentum according to a Gaussian.   
#include "ludat1_hijing.inc"
    
C...Generate p_T and azimuthal angle, gives p_x and p_y.    
      KFLA=IABS(KFL)    
      PT=PARJ(21)*SQRT(-LOG(MAX(1E-10,RLU_HIJING(0)))) 
      IF(MSTJ(91).EQ.1) PT=PARJ(22)*PT  
      IF(KFLA.EQ.0.AND.MSTJ(13).LE.0) PT=0. 
      PHI=PARU(2)*RLU_HIJING(0)    
      PX=PT*COS(PHI)    
      PY=PT*SIN(PHI)    
    
      RETURN    
      END
Commit	Line	Data
e74335a4	1	* $Id$
	2
	3	C*********************************************************************
	4
	5	SUBROUTINE LUPTDI_HIJING(KFL,PX,PY)
	6
	7	C...Purpose: to generate transverse momentum according to a Gaussian.
	8	#include "ludat1_hijing.inc"
	9
	10	C...Generate p_T and azimuthal angle, gives p_x and p_y.
	11	KFLA=IABS(KFL)
	12	PT=PARJ(21)*SQRT(-LOG(MAX(1E-10,RLU_HIJING(0))))
	13	IF(MSTJ(91).EQ.1) PT=PARJ(22)*PT
	14	IF(KFLA.EQ.0.AND.MSTJ(13).LE.0) PT=0.
	15	PHI=PARU(2)*RLU_HIJING(0)
	16	PX=PT*COS(PHI)
	17	PY=PT*SIN(PHI)
	18
	19	RETURN
	20	END