I am trying to understand the idea that gravity breaks down at the Planck scale, but I am confused by the use of natural units ($c = \hbar = 1$). The Einstein-Hilbert action in natural units is:
\begin{equation} S_{EH} = \frac{1}{16 \pi G}\int d^4 x \sqrt{-g} R. \end{equation}
I have seen different arguments, but the simplest one just notes that a perturbative expansion of $S_{EH}$ leads to $E^2/G$ terms so the expansion blows up when $E > \sqrt{G}$. See for example p. 172 of Zee.
In natural units, the Planck length $l_P$ is equal to $\sqrt{G}$. Therefore, it seems that beyond the Planck scale perturbation theory breaks down.
Now I would like to put back all of the constants. From the action itself, I can see that $c$ will be there somewhere. But where does $\hbar$ come from? Please don't tell me it's just dimensional analysis, because all dimensional analysis tells me is that we need some constant with the right dimensions not necessarily equal to $\hbar$. Is there a specific physical reason it is $\hbar$?
EDIT: Of course we expect $\hbar$ because we are interested in quantum gravity. But I want to know where specifically in this sort of calculation $\hbar$ appears - which equations are appealed to?