Using the gravitational path integral we can define the partition function as:
$$ Z(\beta) = \int\mathcal{D}g\mathcal{D}\phi e^{-I_E[g,\phi]} $$ with boundary conditions: $$ t_E \sim t_E + \beta, \quad g_{tt} \to 1 \text{ at infinity.} $$ Here, $t_E$ is the Euclidean time with period $\beta$, $g$ denotes the gravitational field, $\phi$ denotes the matter fields and $I_E[g,\phi]$ is the Euclidean action. In the Schwarzschild solution we can approximate the partition function as: $$ Z(\beta) \approx e^{-I_G[g_0]} $$ where $I_G$ is the Euclidean gravitational action and $g_0$ is the classical Schwarzschild solution.
In order to calculate the partition function in this approximation one needs to perform a Wick rotation in the time coordinate and define the Euclidean Schwarzschild metric. Doing so, one finds that the time coordinate must have a period $\beta = 8\pi M$. In the context of gravitational path integrals $\beta$ is the inverse temperature, so we have the temperature as $T = \frac{1}{8\pi M}$. This is usually interpreted as the temperature of the black hole.
Having found the partition function for the Schwarzschild solution one can calculate an entropy, which turns out to be the Bekenstein-Hawking entropy $S = A/4$. This result is usually interpreted as the entropy of the black hole.
Why do we interpret these results as the temperature and energy of the black hole? Couldn't some of the entropy have come from the Hawking radiation?
In particular I find it curious that the Euclidean metric being used in the calculations have coordinate limits: $$ r \geq 2M, \quad t_E \sim t_E + \beta $$ The gravitational action is then given as an integral where $r$ is in the interval $[2M, \infty]$. In other words, the partition function is given as an integral over the entire manifold except from the black hole! Why do we then think of the entropy as being located in the black hole?