In F.Mandl's book Statisical Physics. He derived Boltsmann's Distribution for a canonical system of N particles to be:
$$ P_r = \frac{e^{-\beta E_r}}{Z} $$
Where Z is the partition function of the system, and $\beta = \frac{1}{kT}$.
This is my uncertainty: My tutor told me that this probability is the probability of one of the mircrostates of Energy = $E_r$ in other words $P_r = \frac{1}{\Omega(E_r)} = \frac{1}{all \hspace{2mm} arrangements \hspace{2mm} with \hspace{2mm} energy \hspace{2mm} E_r}$. Therefore if the canonical system was made up of N subsystems instead of N particles, the partition function would be simply powered by N. Thus the new boltzmann's distribution is given by: $P_r = \frac{e^{-\beta E_r}}{Z^N}$.
However to me, from the derivation, $P_r$ looks more like the probability of finding the state of the system of interest to have energy $E_r$. I guessed this was so, since we started the derivation of $P_r$ in Mandl's textbook to be given by the ratio of statistical weight of the heatbath to have energy corresponding to total energy minus energy of system of interest, $E_{0}-E_r$, to the sum of all statistical weights of all energies.
Please help me clarify my doubts, is my tutor correct, or are my thoughts correct? Your guidance is much appreciated, thank you!