With the effective variable reduction per spectrum, it is now interesting to see how AIC, BIC, and MSE behave when switching from individual fit solutions to one global fit solution. Table shows that the MSE increases up to a factor of 2.2 for complex 1; however, the common AIC and BIC also improve up to a factor of 1.98. Now, with the comparison of the individual versus global fitting results, the user faces the challenge of determining which model is more chemically reasonable. With the global fitting model, the user has to consider that in the presence of strong covalency, the absolute energies may not be the best indicator of the electronic structure (Figure ). This topic has been discussed in the literature. For complex 2 and its MMCT at 7115.4 eV, the peak with the same energy but different amplitude (1: 0.132 a.u. vs. 2: 0.156 a.u. -> ~15% increase in amplitude) is also an important piece to describe the spectra of complex 1 via a fitting model. Finally, we want to emphasize how user-friendly (Formally: <<y_column="2Fe_2C_norm" or "2Fe_2C_DMAP_norm">>, <<now: y_column=["2Fe_2C_norm", "2Fe_2C_DMAP_norm"]>>.) and safe is it is to change from regular fit to global fit by just providing a list of spectra. Figure , along with SI-Advanced Usage III, highlights this aspect. The Open-Source Community aspect welcomes feedback for further improvements and extensions.