Automatic differentiation for error analysis of Monte Carlo data
Abstract
Automatic Differentiation (AD) allows to determine exactly the Taylor series of any function truncated at any order. Here we propose to use AD techniques for Monte Carlo data analysis. We discuss how to estimate errors of a general function of measured observables in different Monte Carlo simulations. Our proposal combines the Γ-method with Automatic differentiation, allowing exact error propagation in arbitrary observables, even those defined via iterative algorithms. The case of special interest where we estimate the error in fit parameters is discussed in detail. We also present a freely available fortran reference implementation of the ideas discussed in this work.
- Publication:
-
Computer Physics Communications
- Pub Date:
- May 2019
- DOI:
- 10.1016/j.cpc.2018.12.020
- arXiv:
- arXiv:1809.01289
- Bibcode:
- 2019CoPhC.238...19R
- Keywords:
-
- Lattice QCD;
- Monte Carlo;
- Error analysis;
- High Energy Physics - Lattice;
- Physics - Computational Physics;
- Physics - Data Analysis;
- Statistics and Probability
- E-Print:
- LaTeX, 25 pags. More explicit notation. Software available in https://gitlab.ift.uam-csic.es/alberto/aderrors