CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-EXO-24-039
Search for resonant and nonresonant production of pairs of dijet resonances with b jets in the final state in proton-proton collisions at $ \sqrt{s} = $ 13 TeV
Abstract: A search is presented for pairs of dijet resonances with identical mass, each containing one jet originating from a b quark and one jet from a light flavor quark. The search uses a data sample corresponding to an integrated luminosity of 138 fb$ ^{-1} $ collected by the CMS detector in proton-proton collisions at $ \sqrt{s}= $ 13 TeV. A machine learning based tagger is used for the identification of at least one b jet within each dijet resonance. Results are presented separately for the case where the four-jet production proceeds via an intermediate resonant state, and for nonresonant production. Upper limits at 95$ % $ confidence level are reported on the production of four-jet and dijet resonances with b jets in the final state. These are the first LHC limits on resonant dijet pair production involving b jets that are outside the standard model. The results are also used to set limits at the 95$ % $ confidence level on the diquark-to-diquark model with masses between 2 and 7 TeV. The nonresonant search excludes pair production of top squarks, with masses between 0.5 and 0.8 TeV, extending previous searches for such $ R $-parity-violating decays to s and b quarks.
Figures Summary References CMS Publications
Figures

png pdf
Figure 1:
(Left) Nonresonant production of pairs of dijet resonances, $ \mathrm{X} $. (Right) Resonant production of pairs of dijet resonances, $ \mathrm{X} $, via a massive resonance, $\mathrm{Y}$. In both scenarios, each dijet resonance, $ \mathrm{X} $, decays to a b jet and a light flavor jet.

png pdf
Figure 1-a:
(Left) Nonresonant production of pairs of dijet resonances, $ \mathrm{X} $. (Right) Resonant production of pairs of dijet resonances, $ \mathrm{X} $, via a massive resonance, $\mathrm{Y}$. In both scenarios, each dijet resonance, $ \mathrm{X} $, decays to a b jet and a light flavor jet.

png pdf
Figure 1-b:
(Left) Nonresonant production of pairs of dijet resonances, $ \mathrm{X} $. (Right) Resonant production of pairs of dijet resonances, $ \mathrm{X} $, via a massive resonance, $\mathrm{Y}$. In both scenarios, each dijet resonance, $ \mathrm{X} $, decays to a b jet and a light flavor jet.

png pdf
Figure 2:
Number of events observed (color scale) within bins of the invariant mass of the four jets, $ m_{\mathrm{4j}} $, and the average dijet mass, $ \overline{m}_{\mathrm{2j}} $.

png pdf
Figure 3:
(Left) Number of events observed (color scale) within bins of $ \overline{m}_{\mathrm{2j}} $ and the ratio $ \alpha $ of that mass to the $ m_{\mathrm{4j}} $. (Right) Number of events predicted in the same bins by a simulation of the production and R-parity violating (RPV) decay of a pair of top squarks with a mass of 1 TeV. The distribution of the simulated signal events arises from the performance of the jet pairing algorithm, and is explained in further detail in Section 6. The dashed lines indicate the lower edges of the three $ \alpha $ bins in which the analysis is conducted.

png pdf
Figure 3-a:
(Left) Number of events observed (color scale) within bins of $ \overline{m}_{\mathrm{2j}} $ and the ratio $ \alpha $ of that mass to the $ m_{\mathrm{4j}} $. (Right) Number of events predicted in the same bins by a simulation of the production and R-parity violating (RPV) decay of a pair of top squarks with a mass of 1 TeV. The distribution of the simulated signal events arises from the performance of the jet pairing algorithm, and is explained in further detail in Section 6. The dashed lines indicate the lower edges of the three $ \alpha $ bins in which the analysis is conducted.

png pdf
Figure 3-b:
(Left) Number of events observed (color scale) within bins of $ \overline{m}_{\mathrm{2j}} $ and the ratio $ \alpha $ of that mass to the $ m_{\mathrm{4j}} $. (Right) Number of events predicted in the same bins by a simulation of the production and R-parity violating (RPV) decay of a pair of top squarks with a mass of 1 TeV. The distribution of the simulated signal events arises from the performance of the jet pairing algorithm, and is explained in further detail in Section 6. The dashed lines indicate the lower edges of the three $ \alpha $ bins in which the analysis is conducted.

png pdf
Figure 4:
(Left) Number of events observed (color scale) within bins of $ m_{\mathrm{4j}} $ and the ratio $ \alpha $. (Right) Number of events predicted in the same bins by a simulation of a diquark with a mass of 6 TeV, decaying to a pair of color-triplet scalar diquarks, each with a mass of 1.5 TeV. The dashed lines indicate the lower edges of the twelve $ \alpha $ bins in which the analysis is conducted.

png pdf
Figure 4-a:
(Left) Number of events observed (color scale) within bins of $ m_{\mathrm{4j}} $ and the ratio $ \alpha $. (Right) Number of events predicted in the same bins by a simulation of a diquark with a mass of 6 TeV, decaying to a pair of color-triplet scalar diquarks, each with a mass of 1.5 TeV. The dashed lines indicate the lower edges of the twelve $ \alpha $ bins in which the analysis is conducted.

png pdf
Figure 4-b:
(Left) Number of events observed (color scale) within bins of $ m_{\mathrm{4j}} $ and the ratio $ \alpha $. (Right) Number of events predicted in the same bins by a simulation of a diquark with a mass of 6 TeV, decaying to a pair of color-triplet scalar diquarks, each with a mass of 1.5 TeV. The dashed lines indicate the lower edges of the twelve $ \alpha $ bins in which the analysis is conducted.

png pdf
Figure 5:
(Left) The $ \overline{m}_{\mathrm{2j}} $ distributions for different $ R $-parity-violating top squark signal masses, with the kinematic and b tagging selections applied. (Right) The $ m_{\mathrm{4j}} $ distributions for different $ \mathrm{S_{uu}} $ diquark masses with $ \alpha_{\mathrm{true}} = $ 0.25, with the kinematic and b tagging selections applied. Shapes are shown for the $ \alpha $ bins with the largest yield for each search, and have been normalized to have an integral of unity.

png pdf
Figure 5-a:
(Left) The $ \overline{m}_{\mathrm{2j}} $ distributions for different $ R $-parity-violating top squark signal masses, with the kinematic and b tagging selections applied. (Right) The $ m_{\mathrm{4j}} $ distributions for different $ \mathrm{S_{uu}} $ diquark masses with $ \alpha_{\mathrm{true}} = $ 0.25, with the kinematic and b tagging selections applied. Shapes are shown for the $ \alpha $ bins with the largest yield for each search, and have been normalized to have an integral of unity.

png pdf
Figure 5-b:
(Left) The $ \overline{m}_{\mathrm{2j}} $ distributions for different $ R $-parity-violating top squark signal masses, with the kinematic and b tagging selections applied. (Right) The $ m_{\mathrm{4j}} $ distributions for different $ \mathrm{S_{uu}} $ diquark masses with $ \alpha_{\mathrm{true}} = $ 0.25, with the kinematic and b tagging selections applied. Shapes are shown for the $ \alpha $ bins with the largest yield for each search, and have been normalized to have an integral of unity.

png pdf
Figure 6:
The products of acceptance and efficiency for a nonresonant signal vs. the top squark mass (left), and a resonant signal vs. the diquark mass (right) inclusively, i.e., for all $ \alpha $ values, and for the three $ \alpha $ bins that contain the majority ($ \geq $ 85%) of the signal. The case where the efficiency of the mass selection is unity is shown as a solid line, and the case where both the mass selection and b tagging efficiencies are unity is shown as a dashed black line.

png pdf
Figure 6-a:
The products of acceptance and efficiency for a nonresonant signal vs. the top squark mass (left), and a resonant signal vs. the diquark mass (right) inclusively, i.e., for all $ \alpha $ values, and for the three $ \alpha $ bins that contain the majority ($ \geq $ 85%) of the signal. The case where the efficiency of the mass selection is unity is shown as a solid line, and the case where both the mass selection and b tagging efficiencies are unity is shown as a dashed black line.

png pdf
Figure 6-b:
The products of acceptance and efficiency for a nonresonant signal vs. the top squark mass (left), and a resonant signal vs. the diquark mass (right) inclusively, i.e., for all $ \alpha $ values, and for the three $ \alpha $ bins that contain the majority ($ \geq $ 85%) of the signal. The case where the efficiency of the mass selection is unity is shown as a solid line, and the case where both the mass selection and b tagging efficiencies are unity is shown as a dashed black line.

png pdf
Figure 7:
(Upper panel) The measured $ \overline{m}_{\mathrm{2j}} $ distributions (points) and background fits (red curves) in each slice of $ \alpha $, in the nonresonant search. (Lower panel) Pulls for the Dijet-3p function, calculated using the statistical uncertainty of the data. In both panels, examples of predicted top squark pair production signals are shown, with cross sections equal to the observed upper limits at 95% confidence level, for top squark masses of 0.6 (blue solid) and 1.0 TeV (blue dashed).

png pdf
Figure 8:
(Upper panel) The measured $ m_{\mathrm{4j}} $ distributions (points) and background fits (red curves) in three out of the twelve slices of $ \alpha $ in the resonant search. (Lower panel) Pulls for the Dijet-3p function, calculated using the statistical uncertainty of the data. In both panels, examples of predicted D2D signals are shown, with cross sections equal to the observed upper limits at 95% confidence level, for $ \mathrm{S_{uu}} $ diquark masses of 2.0 (blue solid), 4.0 (blue dashed), and 6.0 TeV (blue dotted) and for M( $ \mathrm{S_{\frac{2}{3}}} $)/M( $ \mathrm{S_{uu}} $) = 0.25.

png pdf
Figure 9:
Pulls from the fit of the Dijet-3p function to the $ m_{\mathrm{4j}} $ distributions, calculated using the statistical uncertainty of the data, in all twelve $ \alpha $ slices used in the resonant search. The reduced chi-squared of the fit ($ \chi^{2} $ /NDF) for each $ \alpha $ slice is also indicated in each panel.

png pdf
Figure 10:
The observed 95% CL upper limits (black lines with points) on the product of the cross section, branching fraction, and acceptance in the nonresonant search for $ \mathrm{X} \mathrm{X} \rightarrow (\mathrm{b} j)(\mathrm{b} j) $. The expected limits (dashed lines) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to the predicted cross section of the RPV SUSY model [52,53] (dot-dashed).

png pdf
Figure 11:
The observed 95% CL upper limits (black lines with points) on the product of the cross section, branching fraction, and acceptance in the resonant search for $ \mathrm{Y} \rightarrow \mathrm{X} \mathrm{X} \rightarrow (\mathrm{b} j)(\mathrm{b} j) $ with $ \alpha_{\mathrm{true}} = M(\mathrm{X})/M( \mathrm{Y} ) = $ 0.25. The expected limits (dashed lines) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to the predicted cross section of the D2D diquark model [12] (dot-dashed) for M( $ \mathrm{S_{\frac{2}{3}}} $)/M( $ \mathrm{S_{uu}} $) = 0.25.

png pdf
Figure 12:
The observed 95% CL upper limits (black lines with points) on the product of the cross section, branching fraction, and acceptance in the resonant search for $ \mathrm{Y} \rightarrow \mathrm{X} \mathrm{X} \rightarrow (\mathrm{b} j)(\mathrm{b} j) $ for the other twelve values of $ \alpha_{\mathrm{true}} = M(\mathrm{X})/M( \mathrm{Y} ) $. The expected limits (dashed lines) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to the predicted cross section of the D2D diquark model [12] (dot-dashed).

png pdf
Figure 13:
The observed $p$-values corresponding to the nonresonant production of pairs of dijet resonances, each decaying to a b jet and a light flavor jet. The vertical axis represents the local $p$-value for a signal across all $ \alpha $ bins. Dashed lines indicate the corresponding levels of local significance, expressed in units of standard deviation (s.d.).

png pdf
Figure 14:
The observed $p$-values are shown for $ \mathrm{Y} \rightarrow \mathrm{X} \mathrm{X} \rightarrow (\mathrm{b} j)(\mathrm{b} j) $ for thirteen values of $ \alpha_{\mathrm{true}} = M(\mathrm{X})/M( \mathrm{Y} ) $. The vertical axis represents the local $p$-value for a signal across all $ \alpha $ slices. Dashed lines denote the corresponding levels of local significance, expressed in units of standard deviation (s.d.).

png pdf
Figure 15:
Three-dimensional display of the event with the highest $ m_{\mathrm{4j}} = $ 5.5 TeV and highest $ \overline{m}_{\mathrm{2j}} = $ 1.5 TeV. The display shows the energy deposited in the electromagnetic (red) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (green). The two b jets in the event have passing b tagging scores, shown in green, and the two other jets in the event have failing b tagging scores, shown in red. The grouping of the four observed jets into two dijet pairs (purple box) is discussed in the text.
Summary
A search for nonresonant and resonant production of pairs of dijet resonances, with a b jet within each dijet, has been performed. Data from proton-proton collisions at $ \sqrt{s}= $ 13 TeV were used in this search, collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$ ^{-1} $. The steeply falling dijet and four-jet mass distributions are fit with empirical background functions and the simulated shapes of resonance signals. We observe no significant evidence for new resonances. Upper limits at the 95% confidence level (CL) are provided for the production cross section multiplied by the branching fraction and acceptance, for models predicting pair produced dijet resonances, decaying to quarks. These are the first LHC limits on resonant production of pairs of dijet resonances involving b jets, and the first to extend the previous search for nonresonant dijet production to 138 fb$ ^{-1} $ data. For nonresonant production, the limits are presented as a function of the mass of the spin-0 dijet resonance, in the range 0.5 to 2.0 TeV. For resonant production, the limits are given for values of four-jet resonance mass between 2 and 8 TeV. Limits from the nonresonant search are compared to a model [11,10] of $ R $-parity-violating supersymmetry, with pair produced top squarks decaying to bottom and strange quarks. Top squarks with masses between 0.5 and 0.8 TeV are excluded at 95% CL. This significantly extends the prior best mass limit on these top squark decays, which was roughly 0.5 TeV [9]. The most significant signal hypothesis seen in the nonresonant search occurs at a dijet resonance mass of 1.0 TeV and has a local significance of 2.1 standard deviations. Limits from the resonant search are compared to a recent model [12] of massive intermediate scalar diquarks ( $ \mathrm{S_{uu}} $) decaying to pairs of final state diquarks ( $ \mathrm{S_{\frac{2}{3}}} $), and exclude the $ \mathrm{S_{uu}} $ for masses between 2 and 7 TeV, for nearly all values of the $ \mathrm{S_{\frac{2}{3}}} $ mass considered. We note that these are the first limits on this model of $ \mathrm{S_{uu}} $ decays, and that the $ \mathrm{S_{\frac{2}{3}}} $ is indistinguishable from an $ R $-parity-violating top squark. The most significant signal hypothesis seen in the resonant search occurs at an $ \mathrm{S_{uu}} $ mass of 3.1 TeV, and an $ \mathrm{S_{\frac{2}{3}}} $ mass of 1.3 TeV, and has a local significance of 2.6 standard deviations.
References
1 C. Kilic, T. Okui, and R. Sundrum Colored resonances at the Tevatron: Phenomenology and discovery potential in multijets JHEP 07 (2008) 038 0802.2568
2 C. T. Hill Topcolor: Top quark condensation in a gauge extension of the standard model PLB 266 (1991) 419
3 G. D. Kribs, E. Poppitz, and N. Weiner Flavor in supersymmetry with an extended R-symmetry PRD 78 (2008) 055010 0712.2039
4 L. Randall and R. Sundrum A large mass hierarchy from a small extra dimension PRL 83 (1999) 3370 hep-ph/9905221
5 B. A. Dobrescu, R. M. Harris, and J. Isaacson Ultraheavy resonances at the LHC: beyond the QCD background 1810.09429
6 I. Duminica et al. Ultraheavy diquark decaying into vectorlike quarks at the LHC PRD 111 (2025) 115025 2503.17031
7 CMS Collaboration Search for resonant and nonresonant production of pairs of dijet resonances in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 07 (2023) 161 CMS-EXO-21-010
2206.09997
8 ATLAS Collaboration Pursuit of paired dijet resonances in the Run 2 dataset with ATLAS PRD 108 (2023) 112005 2307.14944
9 CMS Collaboration Search for pair-produced resonances decaying to quark pairs in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PRD 98 (2018) 112014 CMS-EXO-17-021
1808.03124
10 R. Barbier et al. R-parity violating supersymmetry Phys. Rept. 420 (2005) 1 hep-ph/0406039
11 J. A. Evans and Y. Kats LHC coverage of RPV MSSM with light stops JHEP 04 (2013) 028 1209.0764
12 B. A. Dobrescu LHC probes of the 10 TeV scale 1912.13155
13 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
14 CMS Collaboration Development of the CMS detector for the CERN LHC Run 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
15 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
16 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
17 CMS Collaboration Performance of the CMS high-level trigger during LHC Run 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
18 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC JINST 16 (2021) P05014 CMS-EGM-17-001
2012.06888
19 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
20 CMS Collaboration Description and performance of track and primary-vertex reconstruction with the CMS tracker JINST 9 (2014) P10009 CMS-TRK-11-001
1405.6569
21 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
22 CMS Collaboration Performance of reconstruction and identification of $ \tau $ leptons decaying to hadrons and $ \nu_\tau $ in pp collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P10005 CMS-TAU-16-003
1809.02816
23 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
24 CMS Collaboration Performance of missing transverse momentum reconstruction in proton-proton collisions at $ \sqrt{s} = $ 13 TeV using the CMS detector JINST 14 (2019) P07004 CMS-JME-17-001
1903.06078
25 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
26 P. Skands, S. Carrazza, and J. Rojo Tuning PYTHIA 8.1: the Monash 2013 tune EPJC 74 (2014) 3024 1404.5630
27 NNPDF Collaboration Parton distributions with LHC data NPB 867 (2013) 244 1207.1303
28 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
29 J. Alwall et al. Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions EPJC 53 (2008) 473 0706.2569
30 T. Sjöstrand et al. An introduction to PYTHIA 8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
31 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
32 NNPDF Collaboration Parton distributions for the LHC Run II JHEP 04 (2015) 040 1410.8849
33 GEANT4 Collaboration Geant4 --- a simulation toolkit Nucl. Instr. Meth. A 506 (2003) 250
34 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
35 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
36 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
37 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
38 CMS Collaboration Jet algorithms performance in 13 TeV data CMS Physics Analysis Summary, CERN, Geneva, 2017
CMS-PAS-JME-16-003
CMS-PAS-JME-16-003
39 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
40 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020) P12012 2008.10519
41 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 Summary CMS-DP-2018-058, 2018
CDS
42 CMS Collaboration Search for narrow resonances decaying to dijets in proton-proton collisions at $ \sqrt{s} = $ 13 TeV PRL 116 (2016) 071801 CMS-EXO-15-001
1512.01224
43 R. A. Fisher On the interpretation of $ \chi^{2} $ from contingency tables, and the calculation of P J. R. Stat. Soc. 85 (1922) 87
44 P. D. Dauncey, M. Kenzie, N. Wardle, and G. J. Davies Handling uncertainties in background shapes: the discrete profiling method JINST 10 (2015) P04015 1408.6865
45 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
46 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s} = $ 13 TeV CMS Physics Analysis Summary, 2018
link
CMS-PAS-LUM-17-004
47 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} = $ 13 TeV CMS Physics Analysis Summary, 2019
link
CMS-PAS-LUM-18-002
48 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
49 ATLAS and CMS Collaborations, and LHC Higgs Combination Group Procedure for the LHC Higgs boson search combination in Summer 2011 Technical Report CMS-NOTE-2011-005, ATL-PHYS-PUB-2011-11, 2011
50 T. Junk Confidence level computation for combining searches with small statistics Nucl. Instr. Meth. A 434 (1999) 435 hep-ex/9902006
51 A. L. Read Presentation of search results: the CLs technique JPG 28 (2002) 2693
52 C. Borschensky et al. Squark and gluino production cross sections in pp collisions at 13, 14, 33 and 100 TeV EPJC 74 (2014) 3174 1407.5066
53 W. Beenakker et al. NLO+NLL squark and gluino production cross-sections with threshold-improved parton distributions EPJC 76 (2016) 53 1510.00375
Compact Muon Solenoid
LHC, CERN