Inferring flavor mixtures in multijet events

Abstract

Multijet events with heavy-flavors are of central importance at the LHC since many relevant processes -- such as $t\bar{t}$, $hh$, $t\bar{t}h$ and others -- have a preferred branching ratio for this final state. Current techniques for tackling these processes use hard-assignment selections through $b$-tagging working points, and suffer from systematic uncertainties because of the difficulties in Monte Carlo simulations. We develop a flexible Bayesian mixture model approach to simultaneously infer $b$-tagging score distributions and the flavor mixture composition in the dataset. We model multidimensional jet events, and to enhance estimation efficiency, we design structured priors that leverages the continuity and unimodality of the $b$-tagging score distributions. Remarkably, our method eliminates the need for a parametric assumption and is robust against model misspecification -- It works for arbitrarily flexible continuous curves and is better if they are unimodal. We have run a toy inferential process with signal $bbbb$ and backgrounds $bbcc$ and $cccc$, and we find that with a few hundred events we can recover the true mixture fractions of the signal and backgrounds, as well as the true $b$-tagging score distribution curves, despite their arbitrariness and nonparametric shapes. We discuss prospects for taking these findings into a realistic scenario in a physics analysis. The presented results could be a starting point for a different and novel kind of analysis in multijet events, with a scope competitive with current state-of-the-art analyses. We also discuss the possibility of using these results in general cases of signals and backgrounds with approximately known continuous distributions and/or expected unimodality.

23 pages, 13 figures. Response to Referees in journal https://scipost.org/submission/2404.01387v2/. No change in reasoning nor conclusion, but a new section comparing new method to usual methods, including a ROC curve comparison, has been added. Also a more realistic case with 1% signal is shown

References in corpus (12)

• Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC
• Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC
• Rank-normalization, folding, and localization: An improved $\hat{R}$ for assessing convergence of MCMC
• Expected Performance of the ATLAS Experiment - Detector, Trigger and Physics
• ATLAS $b$-jet identification performance and efficiency measurement with $t\bar{t}$ events in $pp$ collisions at $\sqrt{s}=13$ TeV
• Bayesian Workflow
• Challenges in Monte Carlo event generator software for High-Luminosity LHC
• Search for nonresonant pair production of Higgs bosons in the $b\bar{b}b\bar{b}$ final state in $pp$ collisions at $\sqrt{s}=13$ TeV with the ATLAS detector
• Probing the $CP$ nature of the top-Higgs Yukawa coupling in $t\bar{t}H$ and $tH$ events with $H\to b\bar{b}$ decays using the ATLAS detector at the LHC
• Bayesian inference to study a signal with two or more decaying particles in a non-resonant background
• Improvement and generalization of ABCD method with Bayesian inference
• On the Modelling of Energetic Multi-jet QCD Events

Cited by in corpus (1)

• Inferring correlated distributions: boosted top jets