Publications

Order by: publication date number of citations

• Variance reduction via simultaneous importance sampling and control variates techniques using vegas

  Prasanth Shyamsundar, Jacob L. Scott, Stephen Mrenna, Konstantin T. Matchev, Kyoungchul Kong

  SciPost Phys. Codebases 28 (2024) · published 4 March 2024 | Toggle abstract · pdf

  Monte Carlo (MC) integration is an important calculational technique in the physical sciences. Practical considerations require that the calculations are performed as accurately as possible for a given set of computational resources. To improve the accuracy of MC integration, a number of useful variance reduction algorithms have been developed, including importance sampling and control variates. In this work, we demonstrate how these two methods can be applied simultaneously, thus combining their benefits. We provide a python wrapper, named CoVVVR, which implements our approach in the VEGAS program. The improvements are quantified with several benchmark examples from the literature.

  1 citation

• Codebase release r1.4 for CoVVVR

  Prasanth Shyamsundar, Jacob L. Scott, Stephen Mrenna, Konstantin T. Matchev, Kyoungchul Kong

  SciPost Phys. Codebases 28-r1.4 (2024) · published 4 March 2024 | Toggle abstract · src

  Monte Carlo (MC) integration is an important calculational technique in the physical sciences. Practical considerations require that the calculations are performed as accurately as possible for a given set of computational resources. To improve the accuracy of MC integration, a number of useful variance reduction algorithms have been developed, including importance sampling and control variates. In this work, we demonstrate how these two methods can be applied simultaneously, thus combining their benefits. We provide a python wrapper, named CoVVVR, which implements our approach in the VEGAS program. The improvements are quantified with several benchmark examples from the literature.

  1 citation

• High-dimensional and permutation invariant anomaly detection

  Vinicius Mikuni, Benjamin Nachman

  SciPost Phys. 16, 062 (2024) · published 29 February 2024 | Toggle abstract · pdf

  Methods for anomaly detection of new physics processes are often limited to low-dimensional spaces due to the difficulty of learning high-dimensional probability densities. Particularly at the constituent level, incorporating desirable properties such as permutation invariance and variable-length inputs becomes difficult within popular density estimation methods. In this work, we introduce a permutation-invariant density estimator for particle physics data based on diffusion models, specifically designed to handle variable-length inputs. We demonstrate the efficacy of our methodology by utilizing the learned density as a permutation-invariant anomaly detection score, effectively identifying jets with low likelihood under the background-only hypothesis. To validate our density estimation method, we investigate the ratio of learned densities and compare to those obtained by a supervised classification algorithm.

  12 citations

• Codebase release 0.5 for mapyde

  Giordon Stark, Camila Aristimuno Ots, Mike Hance

  SciPost Phys. Codebases 27-r0.5 (2024) · published 28 February 2024 | Toggle abstract · src

  Searches for new physics at the Large Hadron Collider have constrained many models of physics beyond the Standard Model. Many searches also provide resources that allow them to be reinterpreted in the context of other models. We describe a reinterpretation pipeline that examines previously untested models of new physics using supplementary information from ATLAS Supersymmetry (SUSY) searches in a way that provides accurate constraints even for models that differ meaningfully from the benchmark models of the original analysis. The public analysis information, such as public analysis routines and serialized probability models, is combined with common event generation and simulation toolkits MADGRAPH, PYTHIA8, and DELPHES into workflows steered by TOML configuration files, and bundled into the mapyde Python package. The use of mapyde is demonstrated by constraining previously untested SUSY models with compressed sleptons and electroweakinos using ATLAS results.

  1 citation

• Reduce, reuse, reinterpret: An end-to-end pipeline for recycling particle physics results

  Giordon Stark, Camila Aristimuno Ots, Mike Hance

  SciPost Phys. Codebases 27 (2024) · published 28 February 2024 | Toggle abstract · pdf

  Searches for new physics at the Large Hadron Collider have constrained many models of physics beyond the Standard Model. Many searches also provide resources that allow them to be reinterpreted in the context of other models. We describe a reinterpretation pipeline that examines previously untested models of new physics using supplementary information from ATLAS Supersymmetry (SUSY) searches in a way that provides accurate constraints even for models that differ meaningfully from the benchmark models of the original analysis. The public analysis information, such as public analysis routines and serialized probability models, is combined with common event generation and simulation toolkits MADGRAPH, PYTHIA8, and DELPHES into workflows steered by TOML configuration files, and bundled into the mapyde Python package. The use of mapyde is demonstrated by constraining previously untested SUSY models with compressed sleptons and electroweakinos using ATLAS results.

  2 citations

• An unfolding method based on conditional invertible neural networks (cINN) using iterative training

  Mathias Backes, Anja Butter, Monica Dunford, Bogdan Malaescu

  SciPost Phys. Core 7, 007 (2024) · published 21 February 2024 | Toggle abstract · pdf

  The unfolding of detector effects is crucial for the comparison of data to theory predictions. While traditional methods are limited to representing the data in a low number of dimensions, machine learning has enabled new unfolding techniques while retaining the full dimensionality. Generative networks like invertible neural networks~(INN) enable a probabilistic unfolding, which map individual data events to their corresponding unfolded probability distribution. The accuracy of such methods is however limited by how well simulated training samples model the actual data that is unfolded. We introduce the iterative conditional INN (IcINN) for unfolding that adjusts for deviations between simulated training samples and data. The IcINN unfolding is first validated on toy data and then applied to pseudo-data for the $pp \to Z \gamma \gamma$ process.

  15 citations

• LHC tau-pair production constraints on $a_\tau$ and $d_\tau$

  Ulrich Haisch, Luc Schnell, Joachim Weiss

  SciPost Phys. 16, 048 (2024) · published 14 February 2024 | Toggle abstract · pdf

  We point out that relevant constraints on the anomalous magnetic ($a_\tau$) and electric ($d_\tau$) moment of the tau lepton can be derived from tau-pair production measurements performed at the LHC. Our conclusion is based on the observation that the leading relative deviations from the Standard Model prediction for $pp \to \tau^+ \tau^-$ due to $a_\tau$ and $d_\tau$ are enhanced at high energies. Less precise measurements at hadron colliders can therefore offer the same or better sensitivity to new physics with respect to high-precision low-energy measurements performed at lepton machines. We derive bounds on $a_\tau$ and $d_\tau$ using the full LHC Run II data set on tau-pair production and compare our findings with the current best limits on the tau anomalous moments.

  6 citations

• Back to the formula - LHC edition

  Anja Butter, Tilman Plehn, Nathalie Soybelman, Johann Brehmer

  SciPost Phys. 16, 037 (2024) · published 29 January 2024 | Toggle abstract · pdf

  While neural networks offer an attractive way to numerically encode functions, actual formulas remain the language of theoretical particle physics. We use symbolic regression trained on matrix-element information to extract, for instance, optimal LHC observables. This way we invert the usual simulation paradigm and extract easily interpretable formulas from complex simulated data. We introduce the method using the effect of a dimension-6 coefficient on associated ZH production. We then validate it for the known case of CP-violation in weak-boson-fusion Higgs production, including detector effects.

  6 citations

• To profile or to marginalize - A SMEFT case study

  Ilaria Brivio, Sebastian Bruggisser, Nina Elmer, Emma Geoffray, Michel Luchmann, Tilman Plehn

  SciPost Phys. 16, 035 (2024) · published 29 January 2024 | Toggle abstract · pdf

  Global SMEFT analyses have become a key interpretation framework for LHC physics, quantifying how well a large set of kinematic measurements agrees with the Standard Model. This agreement is encoded in measured Wilson coefficients and their uncertainties. A technical challenge of global analyses are correlations. We compare, for the first time, results from a profile likelihood and a Bayesian marginalization for a given data set with a comprehensive uncertainty treatment. Using the validated Bayesian framework we analyse a series of new kinematic measurements. For the updated dataset we find and explain differences between the marginalization and profile likelihood treatments.

  6 citations

• Gravitational Regge bounds

  Kelian Häring, Alexander Zhiboedov

  SciPost Phys. 16, 034 (2024) · published 26 January 2024 | Toggle abstract · pdf

  We review the basic assumptions and spell out the detailed arguments that lead to the bound on the Regge growth of gravitational scattering amplitudes. The minimal extra ingredient compared to the gapped case - in addition to unitarity, analyticity, subexponentiality, and crossing - is the assumption that scattering at large impact parameters is controlled by known semi-classical physics. We bound the Regge growth of amplitudes both with the fixed transferred momentum and smeared over it. Our basic conclusion is that gravitational scattering amplitudes admit dispersion relations with two subtractions. For a sub-class of smeared amplitudes, black hole formation reduces the number of subtractions to one. Finally, using dispersion relations with two subtractions we derive bounds on the local growth of relativistic scattering amplitudes. Schematically, the local bound states that the amplitude cannot grow faster than $s^2$. The results obtained in the paper are valid for $d> 4$ for which the $2\to2$ scattering amplitude is well-defined.

  21 citations