\[g(\epsilon) = g_s\left(\frac{L}{2\pi\hbar}\right)^D \int{d^Dp\delta(\epsilon - \epsilon(p))} = g_s\left(\frac{L}{2\pi\hbar}\right)^D \int d\Omega_D \frac{p(\epsilon)^{D-1}}{\left|\frac{\partial \epsilon}{\partial p}(p(\epsilon))\right|}\]