Merging the VirtualMC branch to the main development branch (HEAD)
[u/mrichter/AliRoot.git] / HIJING / hipyset1_35 / luptdi_hijing.F
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