N-dimensional maximum-entropy tomography via particle sampling

Austin Hoover
2024-09-26
Abstract:This paper proposes an alternative implementation of the MENT algorithm, an exact maximum-entropy algorithm used to infer a phase space distribution from its projections. A key step in the MENT algorithm is to compute the distribution's projections via numerical integration. In this approach, the run time scales quickly with the phase space dimension and measurement resolution. The proposed MENT implementation computes the projections via particle sampling, rather than numerical integration, eliminating the dependence on the measurement resolution. Furthermore, with the appropriate sampling algorithm, the particle-based approach scales to $N$-dimensions without computer memory limitations. Using synthetic data, we demonstrate MENT convergence in four dimensions using a grid-based sampling method and in six dimensions using Markov Chain Monte Carlo (MCMC) sampling.
Accelerator Physics,Data Analysis, Statistics and Probability
What problem does this paper attempt to address?