This is actually a wonderful question and I would encourage you to think more about it.
As Emilio said in the comments. The two really describe two different scenarios. In the case of a wavefunction, we consider the system of interest as isolated and there is no notion of temperature (The notion of temperature most notably comes when our system is in contact with a heat reservoir). In such a case, the system has a fixed average energy i.e $\langle\psi| H |\psi\rangle$, that is the expectation value of the Hamiltonian with respect to the wavefunction $|\psi\rangle$.
This is different than the case when we attach our system to a heat reservoir of temperature T. Here, we consider an ensemble of wavefunctions $|\psi_i\rangle$ each of average energy $\epsilon_i$. Calculations of various properties of the system can be done by using the density operator (a.k.a density matrix) $ \sum_i p_i |\psi_i\rangle\langle\psi_i|$, here $p_i $ will be the probability as per the Boltzmann distribution of $|\psi_i\rangle$, the wavefunctions present in the system. (note: one thing you should think is that ask yourself is that wither this different from a superposition of different wavefunctions). This is in a sense how one goes from quantum mechanics to quantum statistical mechanics.