CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-SUS-23-001
Search for top squarks in final states with many light flavor jets and 0, 1, or 2 leptons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV
Abstract: Several new physics models including versions of supersymmetry (SUSY) characterized by $ R $-parity violation (RPV) or with additional hidden sectors predict the production of events with top quarks, low missing transverse momentum, and many additional quarks or gluons. The results of a search for top squarks decaying to two top quarks and six additional light-flavor quarks or gluons are reported. The search employs a novel machine learning method for data-driven background estimation using decorrelated discriminators referred to as the ABCDisCoTEC technique. The search is performed using events with 0, 1, or 2 electrons or muons in conjunction with at least six jets. No requirement is placed on the presence of missing transverse momentum. The result is based on a sample of proton-proton collisions at $ \sqrt{s} = $ 13 TeV corresponding to 138 fb$ ^{-1} $ of integrated luminosity collected with the CMS detector at the LHC in 2016-2018. The data are used to determine upper limits on the top squark pair production cross section in the frameworks of RPV and stealth SUSY. Models with top squark masses less than 700 (920) GeV are excluded at the 95% confidence level for RPV (stealth) SUSY scenarios.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Diagrams of top squark pair production with decays to top quarks and additional light-flavor quarks for the RPV SUSY model (left) and with decays to top quarks and gluons for the stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ model (right).

png pdf
Figure 1-a:
Diagrams of top squark pair production with decays to top quarks and additional light-flavor quarks for the RPV SUSY model (left) and with decays to top quarks and gluons for the stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ model (right).

png pdf
Figure 1-b:
Diagrams of top squark pair production with decays to top quarks and additional light-flavor quarks for the RPV SUSY model (left) and with decays to top quarks and gluons for the stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ model (right).

png pdf
Figure 2:
For the NN training for the RPV signal model and the 1$ \ell $ channel, 1D distributions of $ S_{\mathrm{NN},1} $ (upper left) and $ S_{\mathrm{NN},2} $ (upper right) for data (black solid line), pre-fit simulated SM backgrounds (stacked filled histograms) and two RPV signal models with $ m_{\tilde{\mathrm{t}}}= $ 400 and 800 GeV (dashed lines), as well as the 2D distribution of $ S_{\mathrm{NN},2} $ versus $ S_{\mathrm{NN},1} $ (lower) visualized using a Gaussian kernel density estimation (KDE) for the $ m_{\tilde{\mathrm{t}}} $ = 800 GeV RPV signal model (red) and simulated $ {\mathrm{t}\overline{\mathrm{t}}} {+}$jets background events (gray). The corresponding ABCD bin boundaries for this signal model and channel are also shown in the dashed vertical and horizontal lines.

png pdf
Figure 2-a:
For the NN training for the RPV signal model and the 1$ \ell $ channel, 1D distributions of $ S_{\mathrm{NN},1} $ (upper left) and $ S_{\mathrm{NN},2} $ (upper right) for data (black solid line), pre-fit simulated SM backgrounds (stacked filled histograms) and two RPV signal models with $ m_{\tilde{\mathrm{t}}}= $ 400 and 800 GeV (dashed lines), as well as the 2D distribution of $ S_{\mathrm{NN},2} $ versus $ S_{\mathrm{NN},1} $ (lower) visualized using a Gaussian kernel density estimation (KDE) for the $ m_{\tilde{\mathrm{t}}} $ = 800 GeV RPV signal model (red) and simulated $ {\mathrm{t}\overline{\mathrm{t}}} {+}$jets background events (gray). The corresponding ABCD bin boundaries for this signal model and channel are also shown in the dashed vertical and horizontal lines.

png pdf
Figure 2-b:
For the NN training for the RPV signal model and the 1$ \ell $ channel, 1D distributions of $ S_{\mathrm{NN},1} $ (upper left) and $ S_{\mathrm{NN},2} $ (upper right) for data (black solid line), pre-fit simulated SM backgrounds (stacked filled histograms) and two RPV signal models with $ m_{\tilde{\mathrm{t}}}= $ 400 and 800 GeV (dashed lines), as well as the 2D distribution of $ S_{\mathrm{NN},2} $ versus $ S_{\mathrm{NN},1} $ (lower) visualized using a Gaussian kernel density estimation (KDE) for the $ m_{\tilde{\mathrm{t}}} $ = 800 GeV RPV signal model (red) and simulated $ {\mathrm{t}\overline{\mathrm{t}}} {+}$jets background events (gray). The corresponding ABCD bin boundaries for this signal model and channel are also shown in the dashed vertical and horizontal lines.

png pdf
Figure 2-c:
For the NN training for the RPV signal model and the 1$ \ell $ channel, 1D distributions of $ S_{\mathrm{NN},1} $ (upper left) and $ S_{\mathrm{NN},2} $ (upper right) for data (black solid line), pre-fit simulated SM backgrounds (stacked filled histograms) and two RPV signal models with $ m_{\tilde{\mathrm{t}}}= $ 400 and 800 GeV (dashed lines), as well as the 2D distribution of $ S_{\mathrm{NN},2} $ versus $ S_{\mathrm{NN},1} $ (lower) visualized using a Gaussian kernel density estimation (KDE) for the $ m_{\tilde{\mathrm{t}}} $ = 800 GeV RPV signal model (red) and simulated $ {\mathrm{t}\overline{\mathrm{t}}} {+}$jets background events (gray). The corresponding ABCD bin boundaries for this signal model and channel are also shown in the dashed vertical and horizontal lines.

png pdf
Figure 3:
A comparison of the simulation-based closure correction $ \kappa $ (blue circle), non-closure $ {C}\hspace{-0.70em}/\kern 0.0em \ $ in data (solid purple triangle), and FSR systematic uncertainty (red circle) is shown for each $ N_\text{jets} $ category for the RPV model for the 0$ \ell $ (upper), 1$ \ell $ (middle), and 2$ \ell $ (lower panel) channels. The value of $ {C}\hspace{-0.70em}/\kern 0.0em \ $ shown is the maximum value of the stepping procedure described in the text. The non-closure in the lowest $ N_\text{jets} $ bin is equivalent to the data-based non-closure systematic for each of the channels. Open purple triangles show the statistical uncertainty on $ {C}\hspace{-0.70em}/\kern 0.0em \ $ in data for categories for which all VR and ABCD boundaries have a signal fraction above 5%. All data-based non-closure values are calculated after applying the simulation-based closure correction, such that an observed non-closure of zero signifies identical modeling of the $ S_{\mathrm{NN},1} $-$ S_{\mathrm{NN},2} $ correlation in simulation and data.

png pdf
Figure 4:
$ N_\text{jets} $ distributions from the background-only fits to data. Fits are run with the signal strength fixed to zero and all background process event yields are obtained from their post-fit values. The signal distribution overlaid is for the RPV model with $ m_{\tilde{\mathrm{t}}} = $ 400 GeV. The three plots shown are for the 0$ \ell $, 1$ \ell $, and 2$ \ell $ channels (top to bottom) and the four panels in each plot show the $ N_\text{jets} $ distributions for the $ \mathbf{A} $, $ \mathbf{B} $, $ \mathbf{C} $, and $ \mathbf{D} $ regions (left to right). The grey error band shows the combined statistical and systematic uncertainty from the background-only post-fit distributions.

png pdf
Figure 4-a:
$ N_\text{jets} $ distributions from the background-only fits to data. Fits are run with the signal strength fixed to zero and all background process event yields are obtained from their post-fit values. The signal distribution overlaid is for the RPV model with $ m_{\tilde{\mathrm{t}}} = $ 400 GeV. The three plots shown are for the 0$ \ell $, 1$ \ell $, and 2$ \ell $ channels (top to bottom) and the four panels in each plot show the $ N_\text{jets} $ distributions for the $ \mathbf{A} $, $ \mathbf{B} $, $ \mathbf{C} $, and $ \mathbf{D} $ regions (left to right). The grey error band shows the combined statistical and systematic uncertainty from the background-only post-fit distributions.

png pdf
Figure 4-b:
$ N_\text{jets} $ distributions from the background-only fits to data. Fits are run with the signal strength fixed to zero and all background process event yields are obtained from their post-fit values. The signal distribution overlaid is for the RPV model with $ m_{\tilde{\mathrm{t}}} = $ 400 GeV. The three plots shown are for the 0$ \ell $, 1$ \ell $, and 2$ \ell $ channels (top to bottom) and the four panels in each plot show the $ N_\text{jets} $ distributions for the $ \mathbf{A} $, $ \mathbf{B} $, $ \mathbf{C} $, and $ \mathbf{D} $ regions (left to right). The grey error band shows the combined statistical and systematic uncertainty from the background-only post-fit distributions.

png pdf
Figure 4-c:
$ N_\text{jets} $ distributions from the background-only fits to data. Fits are run with the signal strength fixed to zero and all background process event yields are obtained from their post-fit values. The signal distribution overlaid is for the RPV model with $ m_{\tilde{\mathrm{t}}} = $ 400 GeV. The three plots shown are for the 0$ \ell $, 1$ \ell $, and 2$ \ell $ channels (top to bottom) and the four panels in each plot show the $ N_\text{jets} $ distributions for the $ \mathbf{A} $, $ \mathbf{B} $, $ \mathbf{C} $, and $ \mathbf{D} $ regions (left to right). The grey error band shows the combined statistical and systematic uncertainty from the background-only post-fit distributions.

png pdf
Figure 5:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 5-a:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 5-b:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 5-c:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 5-d:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 5-e:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 5-f:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (left) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (right) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the 0$ \ell $ channel (upper), 1$ \ell $ channel (middle), and 2$ \ell $ channel (lower). The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 6:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (upper) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (lower) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the combination of all three channels. The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 6-a:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (upper) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (lower) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the combination of all three channels. The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.

png pdf
Figure 6-b:
95% CL upper limit on $ \sigma_{\tilde{\mathrm{t}}\bar{\tilde{\mathrm{t}}}} $ for the RPV (upper) and stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ (lower) SUSY signal models as a function of $ m_{\tilde{\mathrm{t}}} $ for the combination of all three channels. The median expected limit is shown in the dashed blue line, with the 68% and 95% intervals shown in light-blue and yellow, respectively. The observed limit is shown in black. The vertical dashed line denotes the transition from the low-mass optimization to high-mass optimization ABCD boundaries. Additionally, the theoretical cross section is shown in red with its uncertainty in light brown.
Tables

png pdf
Table 1:
The list of the search region selections for the 0$ \ell $, 1$ \ell $, and 2$ \ell $ channels. $ N_{\text{leptons}}^{\text{iso}} $ and $ N_{\text{muon}}^{\text{non-iso}} $ denote the number of isolated leptons and nonisolated muons, respectively.

png pdf
Table 2:
Lower bounds on $ (S_{\mathrm{NN},1},S_{\mathrm{NN},2}) $ defining the $ \mathbf{A} $ region for the low- and high-mass optimization for each signal model and search channel.

png pdf
Table 3:
Magnitude of systematic uncertainties for the 0$ \ell $, 1$ \ell $, and 2$ \ell $ channels based on the RPV-trained NN and ABCD bin boundaries optimized for low $ m_{\tilde{\mathrm{t}}} $. Reported values are in units of % and ranges correspond to the 16th and 84th percentile for the value of a systematic uncertainty across all applicable analysis regions (ABCD regions and $ N_\text{jets} $ categories). The maximum value for a given systematic uncertainty across these regions is shown in parentheses. Other backgrounds include the QCD multijet and minor background processes. The RPV signal model used has $ m_{\tilde{\mathrm{t}}} = $ 550 GeV. The systematic values based on the RPV-trained model are representative of the values obtained for the $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $-trained model.
Summary
A search for the pair production of top squarks with decays to top quarks and six additional gluons or light-flavor quarks via stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ and RPV SUSY decays is presented. The search is performed using proton-proton collision events at $ \sqrt{s} = $ 13 TeV, corresponding to an integrated luminosity of 138 fb$ ^{-1} $, collected by the CMS detector between 2016 and 2018. Events are selected in three search channels, defined as having at least six jets and zero, one, or two leptons. No requirement is placed on the presence of $ p_{\mathrm{T}}^\text{miss} $. This analysis is an extension of a previous search for these signatures [7], which observed a deviation with local significance of 2.8 standard deviations for a top squark mass of 400 GeV for the RPV model. The main improvements of this search are the addition of the zero- and two-lepton channels as well as the inclusion of a novel, neural-network-based background estimation method referred to as ABCDisCoTEC. The key feature of this new method is the creation of two uncorrelated neural network variables that can be used with an ABCD-style background estimation method. The backgrounds are estimated from a simultaneous binned likelihood fit to all search channels in several categories of jet multiplicity, with the contribution from the main $ {\mathrm{t}\overline{\mathrm{t}}} {+}$jets background constrained via the ABCD relationship that encodes the independence between the two neural network variables. This new background estimation method reduced the dependence of the analysis on uncertainties related to the modeling of the jet multiplicity spectrum compared to Ref. [7]. Overall, good agreement between the data and the background prediction is observed, and the deviation observed previously is not confirmed. The results are interpreted using two top squark decay topologies which generate low-$ p_{\mathrm{T}}^\text{miss} $ signatures: the stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ and R-parity violating models. Top squark masses up to 700 GeV for the RPV model with top squark decays to a top quark and neutralino with subsequent decay of the neutralino to three light-flavored quarks are excluded at 95% confidence level. An upper mass exclusion of 930 GeV is placed on the stealth $ \mathrm{S} \mathrm{Y} \overline{\mathrm{Y}} $ model where top squarks decay via a top quark, gluons, and a soft gravitino via a hidden sector.
References
1 P. Fayet and S. Ferrara Supersymmetry Phys. Rept. 32 (1977) 249
2 S. P. Martin A supersymmetry primer Adv. Ser. Direct. High Energy Phys. 21 (1997) 1 hep-ph/9709356
3 R. Barbier et al. $ R $-parity violating supersymmetry Phys. Rept. 420 (2005) 1 hep-ph/0406039
4 J. Fan, M. Reece, and J. T. Ruderman Stealth supersymmetry JHEP 11 (2011) 012 1105.5135
5 J. Fan, M. Reece, and J. T. Ruderman A stealth supersymmetry sampler JHEP 07 (2012) 196 1201.4875
6 J. Fan et al. Stealth supersymmetry simplified JHEP 07 (2016) 016 1512.05781
7 CMS Collaboration Search for top squarks in final states with two top quarks and several light-flavor jets in proton-proton collisions at $ \sqrt {s} = $ 13 TeV PRD 104 (2021) 032006 CMS-SUS-19-004
2102.06976
8 G. Kasieczka, B. Nachman, M. D. Schwartz, and D. Shih Automating the ABCD method with machine learning PRD 103 (2021) 035021 2007.14400
9 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s}= $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
10 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2018
CMS-PAS-LUM-17-004
CMS-PAS-LUM-17-004
11 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2019
CMS-PAS-LUM-18-002
CMS-PAS-LUM-18-002
12 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
13 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
14 CMS Collaboration Technical proposal for the Phase-II upgrade of the Compact Muon Solenoid CMS Technical Proposal CERN-LHCC-2015-010, CMS-TDR-15-02, 2015
CDS
15 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
16 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
17 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
18 CMS Collaboration Jet performance in pp collisions at $ \sqrt{s}= $ 7 TeV CMS Physics Analysis Summary, 2010
link
19 CMS Collaboration Jet algorithms performance in 13 TeV data CMS Physics Analysis Summary, 2017
CMS-PAS-JME-16-003
CMS-PAS-JME-16-003
20 CMS Collaboration Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV JINST 12 (2017) P02014 CMS-JME-13-004
1607.03663
21 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
22 CMS Collaboration Identification of heavy-flavour jets with the CMS detector in pp collisions at 13 TeV JINST 13 (2018) P05011 CMS-BTV-16-002
1712.07158
23 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020) P12012 2008.10519
24 CMS Collaboration Performance of the DeepJet b tagging algorithm using 41.9/fb of data from proton-proton collisions at 13 TeV with Phase 1 CMS detector CMS Detector Performance Note CMS-DP-2018-058, 2018
CDS
25 CMS Collaboration Performance of electron reconstruction and selection with the CMS detector in proton-proton collisions at $ \sqrt{s} = $ 8 TeV JINST 10 (2015) P06005 CMS-EGM-13-001
1502.02701
26 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
27 K. Rehermann and B. Tweedie Efficient identification of boosted semileptonic top quarks at the LHC JHEP 03 (2011) 059 1007.2221
28 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
29 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
30 S. Alioli, P. Nason, C. Oleari, and E. Re A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX JHEP 06 (2010) 043 1002.2581
31 S. Frixione, P. Nason, and G. Ridolfi A positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction JHEP 09 (2007) 126 0707.3088
32 E. Re Single-top Wt-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
33 S. Alioli, P. Nason, C. Oleari, and E. Re NLO single-top production matched with shower in POWHEG: s- and t-channel contributions JHEP 09 (2009) 111 0907.4076
34 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
35 T. Sjöstrand et al. An introduction to PYTHIA 8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
36 M. Czakon and A. Mitov Top++: A program for the calculation of the top-pair cross-section at hadron colliders Comput. Phys. Commun. 185 (2014) 2930 1112.5675
37 P. Kant et al. HATHOR for single top-quark production: Updated predictions and uncertainty estimates for single top-quark production in hadronic collisions Comput. Phys. Commun. 191 (2015) 74 1406.4403
38 M. Aliev et al. HATHOR: HAdronic Top and Heavy quarks crOss section calculatoR Comput. Phys. Commun. 182 (2011) 1034 1007.1327
39 T. Gehrmann et al. W$ ^+ $W$ ^- $ production at hadron colliders in next to next to leading order QCD PRL 113 (2014) 212001 1408.5243
40 J. M. Campbell and R. K. Ellis An update on vector boson pair production at hadron colliders PRD 60 (1999) 113006 hep-ph/9905386
41 J. M. Campbell, R. K. Ellis, and C. Williams Vector boson pair production at the LHC JHEP 07 (2011) 018 1105.0020
42 Y. Li and F. Petriello Combining QCD and electroweak corrections to dilepton production in FEWZ PRD 86 (2012) 094034 1208.5967
43 C. Borschensky et al. Squark and gluino production cross sections in pp collisions at $ \sqrt{s}= $ 13, 14, 33 and 100 TeV EPJC 74 (2014) 3174 1407.5066
44 W. Beenakker et al. NNLL-fast: predictions for coloured supersymmetric particle production at the LHC with threshold and Coulomb resummation JHEP 12 (2016) 133 1607.07741
45 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
46 CMS Collaboration CMS pythia 8 colour reconnection tunes based on underlying-event data EPJC 83 (2023) 587 CMS-GEN-17-002
2205.02905
47 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
48 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
49 CMS Collaboration Identification of heavy, energetic, hadronically decaying particles using machine-learning techniques JINST 15 (2020) P06005 CMS-JME-18-002
2004.08262
50 CMS Collaboration Search for top squark production in fully-hadronic final states in proton-proton collisions at $ \sqrt{s} = $ 13 TeV PRD 104 (2021) 052001 CMS-SUS-19-010
2103.01290
51 O. Behnke, K. Kröninger, G. Schott, and T. Schörner-Sadenius Data analysis in high energy physics: a practical guide to statistical methods Wiley-VCH, 2013
link
52 G. J. Székely and M. L. Rizzo Partial distance correlation with methods for dissimilarities Ann. Statist. 42 (2014) 2382 1310.2926
53 S. Odaka, N. Watanabe, and Y. Kurihara ME-PS matching in the simulation of multi-jet production in hadron collisions using a subtraction method PTEP 2015 (2015) 053B04 1409.3334
54 G. C. Fox and S. Wolfram Observables for the analysis of event shapes in e+ e- annihilation and other processes PRL 41 (1978) 1581
55 J. D. Bjorken and S. J. Brodsky Statistical model for electron-positron annihilation into hadrons PRD 1 (1970) 1416
56 CMS Collaboration Search for supersymmetry in hadronic final states using MT2 in $ pp $ collisions at $ \sqrt{s} = $ 7 TeV JHEP 10 (2012) 018 CMS-SUS-12-002
1207.1798
57 CMS Collaboration CMS technical design report, volume II: Physics performance JPG 34 (2007) 995
58 CMS Collaboration The CMS statistical analysis and combination tool: Combine Submitted to Comput. Softw. Big Sci, 2024 CMS-CAT-23-001
2404.06614
Compact Muon Solenoid
LHC, CERN