PeterDonis said:
Yes.
No, you can't, because there is no way to distinguish "the already existing radiation" from "additional radiation" that you added. The radiation isn't going to have a well-defined photon number anyway since it's not going to be in a Fock state; it will basically be black-body radiation. But the radiation will have a well-defined temperature, and that temperature will increase when you apply the blowtorch. So it certainly seems appropriate to say that you heated the radiation inside.
Just to add a bit more: If you have a container and heat the walls, the electrons contained in the walls are rattling more and more, and this random motion leads to em. radiation. Keeping the walls at a certain constant temperature the corresponding radiation inside will also come to thermal equilibrium. This process is easily understood in a kinetic-theory way: There are photons created and absorbed all the time until everything comes to thermal equilibrium, where the absorption and creation rate becomes equal. Fortunately we don't need to solve this kinetic problem to know the equilibrium state of the radiation, but that's given by the maximum-entropy principle (equilibrium is the state of maximal entropy at the given constraints; in our discussed case it's the given temperature of the walls which equivalently means that the mean energy density of the radiation is fixed).
For simplicity let's make the quantization volume a cube of length ##L## and take it as a subvolume within the very big total volume of the cavity. Then we can impose periodic boundary conditions, which makes the calculation a bit simpler. This means a complete basis of the electromagnetic field are bosonic Fock states with the single-particle states chosen as momentum-helicity eigenstates. The momenta ##\vec{p} \in \frac{2 \pi}{L} \mathbb{Z}## and the helicities take the values ##\lambda=\pm 1## (I use natural units with ##\hbar=c=k_{\text{B}}=1##). The Hamiltonian is given by
$$\hat{H}=\sum_{\vec{p},\lambda} E_{\vec{p}} \hat{N}(\vec{p},\lambda),$$
Where ##\hat{N}(\vec{p},\lambda)## are the number operators for photons in the wave mode given by ##\vec{p}## and ##\lambda##. The energies are ##E(\vec{p})=|\vec{p}|## since photons are massless.
The partition sum thus is given by (with ##\beta=1/T## the inverse temperature)
$$Z=\mathrm{Tr} \exp(-\beta \hat{H}).$$
For the following it's more convenient to first calculate a somewhat more general quantity in order to be able to calculate various interesting expectation values, namely
\begin{equation*}
\begin{split}
Z[\alpha]&=\mathrm{Tr} \exp(-\beta \hat{H} + \sum_{\vec{p},\lambda} \alpha(\vec{p},\lambda) \hat{N}(\vec{p},\lambda) = \mathrm{Tr} \exp(\sum_{\vec{p},\lambda}[-\beta E(\vec{p})+\alpha(\vec{p},\lambda)] \hat{N}(\vec{p},\lambda)] \\
&=\prod_{\vec{p},\lambda} \sum_{N(\vec{p},\lambda)=0}^{\infty} \exp[(-\beta E(\vec{p})+\alpha(\vec{p},\lambda))N(\vec{p},\lambda)] \\
&=\prod_{\vec{p},\lambda} \frac{1}{1-\exp(-\beta E(\vec{p})+\alpha(\vec{p},\lambda))}.
\end{split}
\end{equation*}
Now it's easy to calculate the mean number of photons in each wave mode. Given that the Statistical operator is
$$\hat{\rho}=\frac{1}{Z} \exp(-\beta \hat{H})=\frac{1}{Z} \left [\frac{\partial Z[\alpha]}{\partial \alpha(\vec{p},\lambda)} \right]_{\alpha(\vec{p},\lambda)=0} = \left [ \frac{\partial}{\partial \alpha(\vec{p},\lambda)} \ln Z[\alpha] \right]_{\alpha(\vec{p},\lambda)=0}.$$
Which leads after some algebra to
$$n(\vec{p},\lambda)=\langle N(\vec{p},\lambda) \rangle=\frac{1}{\exp(\beta E(\vec{p}))-1},$$
i.e., the Planck black-body distribution as expected.
This doesn't depend on ##\lambda##, i.e., the radiation is completely unpolarized (which is also expected from the equlibrium state being the state of maximal disorder).
BTW: A perfectly reflecting wall won't lead to thermal equilibrium at all since no em. radiation (or photons) can be absorbed. It's perfectly isolating and thus you can't "heat up the radiation", and you'd just have the one "green" cavity mode excited all the time. On the other hand, if the box is not a perfect mirror inside (and that's the case of course for any real-world cavity, although you can indeed make quite perfectly reflecting cavities nowadays) the initially excited radiation will get partially absorbed and the cavity walls heat up and at the end you have Planck radiation (the better the walls are reflecting the longer this will take) with the temperature determined by the total energy within the cavity and the walls.