Abstract: Precise modeling of a material is a key to understand its underlying interactions and physics but also revealing the competing phases in the nearby interaction space. Highly frustrated systems are important due to the richness in physics and diversity of phases including spin liquids with exotic topological states they display. Here, we present a machine learning workflow to fit multi-experimental datasets to find an optimal Hamiltonian while undertaking phase classification and extracting information about the topography around the region of interest. Experimental data from a spin-ice material, Dy₂Ti₂O₇ including diffuse neutron scattering, heat capacity and susceptibility are utilized. This approach is shown to provide the best model in an efficient and effective way but also is powerful at planning the best experimental strategies.