New functions and constructors added and some other fixes.
[u/mrichter/AliRoot.git] / PYTHIA / jetset / luptdi.F
1  
2 C********************************************************************* 
3  
4       SUBROUTINE LUPTDI(KFL,PX,PY) 
5  
6 C...Purpose: to generate transverse momentum according to a Gaussian. 
7       COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200) 
8       SAVE /LUDAT1/ 
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(0)))) 
13       IF(PARJ(23).GT.RLU(0)) PT=PARJ(24)*PT 
14       IF(MSTJ(91).EQ.1) PT=PARJ(22)*PT 
15       IF(KFLA.EQ.0.AND.MSTJ(13).LE.0) PT=0. 
16       PHI=PARU(2)*RLU(0) 
17       PX=PT*COS(PHI) 
18       PY=PT*SIN(PHI) 
19  
20       RETURN 
21       END