We want to be able to use multiple reference simulations, each potentially run with different force field parameters or sequences, as input to the reweighting framework. In addition, we want to allow for each reference simulation to have a different set of temperatures, since optimizing the temperature spacing for constant entropy increase leads to different temperature sets for each system. This means we need to evaluate the energies at all possible temperatures, and also their finite difference intervals.

This will be very useful for increasing configurational sampling across ranges of parameters / sequences. For example, if we run a homopolymer A, homopolymer B, and a few AB copolymers with different LJ epsilon parameters in A and B sidechains, a lot more configurational space is sampled than just using homopolymer A as a reference.