Preliminary tests demonstrated that this provided the best setup to compare g(r) curves as a whole in the most robust way. This modification also allows mitigating potential problems emerging from the fact that classical three-site water models usually have difficulty reproducing the first peak of the RDFs (due to the fact that quantum effects are not included in the description of the system). In particular, the scoring function used in the optimization presented the first subsection of the results is: