Metal hydrides are promising candidates for hydrogen storage applications. From a materials discovery perspective, an accurate, efficient computational workflow is urgently required that can rapidly analyze/predict thermodynamic properties of these materials. In this paper, we develop a thermodynamic property estimation/phase diagram construction framework based on the lattice reverse Monte Carlo (RMC) method. We apply the technique to the nickel hydride (NiHx) system by calculating the pressure-composition-temperature (PCT) isotherm and constructing its phase diagram. Formation of nickel hydride entails significant volume expansion, strong interaction between hydrogen and the host atoms, lattice strain, and a phase transition. An attractive feature of our approach is that the entire phase diagram can be accurately constructed in few minutes by considering <10 configurations. In contrast, a popular technique based on grand canonical Monte Carlo would require sampling of several million configurations. The promising results obtained with our computational workflow suggest that the approach can be used in future for wider materials search and discovery.