Stay Positive: Neural Refinement of Sample Weights
arXiv:2505.03724 · doi:10.1103/yx16-h7n9
Abstract
Monte Carlo simulations are an essential tool in particle physics data analysis. Events are typically generated alongside weights that redistribute the cross section of the simulated process across the phase space. These weights can be negative, and several post-hoc methods have been developed to eliminate or mitigate the negative values. All of these methods share the common strategy of approximating the average weight as a function of phase space. We introduce an alternative approach, which, instead of reweighting to the average, refines the initial weights with a scaling transformation, utilizing a phase space-dependent factor. Since this new refinement method does not need to model the full weight distribution, it can be more accurate. High-dimensional and unbinned phase space is processed using neural networks for the refinement method. In addition to the refinement method, we introduce a new resampling protocol, which can be used in conjunction with any weight transformation to not only preserve the average weight but also the statistical uncertainties of the initial distribution. Using both realistic and synthetic examples, we show that the new neural refinement method is able to match or exceed the accuracy of similar weight transformations and that the new resampling protocol is simpler in implementation than previous methods while exhibiting equivalent statistical properties.
10 pages, 7 figures; enhanced readability, added resampling protocol
References in corpus (22)
- PYTHIA 6.4 Physics and Manual
- PYTHIA 6.2 Physics and Manual
- The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations
- An Introduction to PYTHIA 8.2
- FastJet user manual
- DELPHES 3, A modular framework for fast simulation of a generic collider experiment
- A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX
- Parton distributions with LHC data
- Dispelling the N^3 myth for the Kt jet-finder
- Three-jet cross sections to next-to-leading order
- CutTools: a program implementing the OPP reduction method to compute one-loop amplitudes
- Automation of one-loop QCD corrections
- OneLOop: for the evaluation of one-loop scalar functions
- A general approach to jet cross sections in QCD
- Automated one-loop calculations: a proof of concept
- Ninja: Automated Integrand Reduction via Laurent Expansion for One-Loop Amplitudes
- Integrand reduction of one-loop scattering amplitudes through Laurent series expansion
- A Neural Resampler for Monte Carlo Reweighting with Preserved Uncertainties
- On the reduction of negative weights in MC@NLO-type matching procedures
- A Positive Resampler for Monte Carlo Events with Negative Weights
- Unbiased Elimination of Negative Weights in Monte Carlo Samples
- Efficient negative-weight elimination in large high-multiplicity Monte Carlo event samples