We have a molecule that is emitting/absorbing photons. We know the Hamiltonian and that there are several levels. We count the emitted photons at different angles and frequencies. We can also do scattering with a beam of photons. From the absorption data, can we determine if the system occupies each energy level with a thermal probability distribution? I.e. I mean that an energy eigenstate $i$ with energy $E_i$ is occupied with probability proportional to $e^{-\beta E_i}$. How do we do this experimentally?
As a simple example, how can we tell if a hydrogen atom is in a thermal state? Can we generalize this to bigger molecules?