Error propagation and analysis of Monte Carlo data with the \( \Gamma \)-method and automatic differentiation. There are Julia and fortran implementations.
Error propagation is performed exactly, even in quantities that are defined via an iterative algorithm (i.e. fit parameters). The paper (arXiv link) explains the benefits of such an approach.
Currently there exist two implementations: