luptdi.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file luptdi.f

    
C*********************************************************************  
    
      SUBROUTINE luptdi(KFL,PX,PY)  
    
C...Purpose: to generate transverse momentum according to a Gaussian.   
      common/ludat1/mstu(200),paru(200),mstj(200),parj(200) 
      SAVE /ludat1/ 
    
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(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(0)    
      px=pt*cos(phi)    
      py=pt*sin(phi)    
    
      RETURN    
      END