The Hawking temperature of a Schwarzschild black hole is given in SI units as $$T_{H}=\frac{\hbar c^3}{8 \pi G k_{B} M},$$ where $k_{B}$ is the Boltzmann constant. I would like to know how $\hbar$ and $k_{B}$ show up in the temperature. I mean where in the original derivation by Hawking do these constants show up?
I have looked into the original paper by Hawking, "Particle Creation by Black Holes". There he begins the calculation by writing down the massless scalar wave equation in curved background $$ \frac{1}{\sqrt{-g}}\partial_{\mu} \left(\sqrt{-g}\,g^{\mu\nu} \partial_{\nu} \phi \right)=0.$$
Now as far as I can understand Hawking temperature shows up in the exponent when the modes $\sim\phi$ are traced from the surface of collapsing body to the future infinity. Alternative calculations without invoking the collapse geometry suggest modes tunnel through the horizon. So at first I thought $\hbar$ naturally shows because from quantum mechanics we have $\phi \sim e^{\frac{i}{\hbar}S}$. But what is bothering me here is that above wave equation does not have any $\hbar$ in it. In fact such wave equation comes from a Lagrangian of the form $$ I\left[\phi\right] = \int \left[\frac{1}{2}g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi \right]\sqrt{-g} d^4x,$$ where $\hbar$ does not show up. And I am not sure whether such Lagrangian should come with a $\hbar$ based on dimensional analysis.
My second source of confusion is related to the Boltzmann constant. Again, I have no idea how and where $k_{B}$ emerges in the derivation. Without the notion of temperature, $k_{B}$ seems unrelated to such calculations which involve wave equations, Bogoliubov transformations etc...