 Rotating Mestel disk distribution function

In the midplane Mestel disk potential $\Phi =v_c^2\log (R/R_0)$, take specific energy $E=(v_R^2+v_{\varphi }^2)/2+\Phi$ and $L_z=R{v}_{\varphi }$. A prograde galactic distribution function $f=A{L}_z^q{e}^{-E/s}$ for $L_z>0$, with $s>0$, has surface density proportional to $R^{q-v_c^2/s}$. Matching $\Sigma ={\Sigma }_0{R}_0/R$ gives $q=v_c^2/s-1>-1$ andThe denominator is the product of the full radial Gaussian integral and the positive-azimuthal half-line moment. Here $s$ is a squared-velocity scale, not a velocity dispersion itself. Changing the potential's additive constant changes $A$. With positive exponent $e^{\epsilon /s}$, $\epsilon$ must instead be the negative of the specific total energy.