We study an ensemble of quantum pure states, the thermalization resilient ensemble (TRE), providing the statistical characterization of the thermal equilibrium of 5 isolated quantum systems. Following a previous work where the ensemble was defined based on the invariance of the average populations upon thermal contact of identical systems, here we introduce a general methodology to generate quantum states according to the TRE statistic. The sampling is employed to characterize the ensemble distribution of thermodynamic functions like the entropy, internal energy, and temperature. The possibility of defining the temperature also for isolated quantum systems with a limited number of degrees of freedom is a distinctive feature of the TRE statistic which has no counterpart in other quantum statistical ensembles. The results are illustrated by explicit calculations for spin model systems.