ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
luptdi.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file 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(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(0)
16  px=pt*cos(phi)
17  py=pt*sin(phi)
18 
19  RETURN
20  END