Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:15:43

0001     
0002 C*********************************************************************  
0003     
0004       SUBROUTINE LUPTDI(KFL,PX,PY)  
0005     
0006 C...Purpose: to generate transverse momentum according to a Gaussian.   
0007       COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200) 
0008       SAVE /LUDAT1/ 
0009     
0010 C...Generate p_T and azimuthal angle, gives p_x and p_y.    
0011       KFLA=IABS(KFL)    
0012       PT=PARJ(21)*SQRT(-LOG(MAX(1E-10,RLU(0)))) 
0013       IF(MSTJ(91).EQ.1) PT=PARJ(22)*PT  
0014       IF(KFLA.EQ.0.AND.MSTJ(13).LE.0) PT=0. 
0015       PHI=PARU(2)*RLU(0)    
0016       PX=PT*COS(PHI)    
0017       PY=PT*SIN(PHI)    
0018     
0019       RETURN    
0020       END