CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-EXO-19-021 ; CERN-EP-2020-202
Search for long-lived particles using displaced jets in proton-proton collisions at $\sqrt{s} = $ 13 TeV
Phys. Rev. D 104 (2021) 012015
Abstract: An inclusive search is presented for long-lived particles using displaced jets. The search uses a data sample collected with the CMS detector at the CERN LHC in 2017 and 2018, from proton-proton collisions at a center-of-mass energy of 13 TeV. The results of this search are combined with those of a previous search using a data sample collected with the CMS detector in 2016, yielding a total integrated luminosity of 132 fb$^{-1}$. The analysis searches for the distinctive topology of displaced tracks and displaced vertices associated with a dijet system. For a simplified model, where pair-produced long-lived neutral particles decay into quark-antiquark pairs, pair production cross sections larger than 0.07 fb are excluded at 95% confidence level (CL) for long-lived particle masses larger than 500 GeV and mean proper decay lengths between 2 and 250 mm. For a model where the standard model-like Higgs boson decays to two long-lived scalar particles that each decays to a quark-antiquark pair, branching fractions larger than 1% are excluded at 95% CL for mean proper decay lengths between 1 mm and 340 mm. A group of supersymmetric models with pair-produced long-lived gluinos or top squarks decaying into various final-state topologies containing displaced jets is also tested. Gluino masses up to 2500 GeV and top squark masses up to 1600 GeV are excluded at 95% CL for mean proper decay lengths between 3 and 300 mm. The highest lower bounds on mass reach 2600 GeV for long-lived gluinos and 1800 GeV for long-lived top squarks. These are the most stringent limits to date on these models.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
The Feynman diagrams for the different long-lived models considered, including the jet-jet model (upper left), models with an exotic decay of the SM-like Higgs boson (upper right), general gauge mediation models with ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ decay (second row, left), mini-split SUSY with ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ decay (second row, right), RPV SUSY with ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ decay (third row, left), RPV SUSY with $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ decay (third row, right), RPV SUSY with $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ decay (lower left), and dRPV SUSY with $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ decay (lower right).

png pdf
Figure 1-a:
Feynman diagram for the jet-jet model.

png pdf
Figure 1-b:
Feynman diagram for models with an exotic decay of the SM-like Higgs boson.

png pdf
Figure 1-c:
Feynman diagram for general gauge mediation models with ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ decay.

png pdf
Figure 1-d:
Feynman diagram for mini-split SUSY with ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ decay.

png pdf
Figure 1-e:
Feynman diagram for RPV SUSY with ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ decay.

png pdf
Figure 1-f:
Feynman diagram for RPV SUSY with $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ decay.

png pdf
Figure 1-g:
Feynman diagram for RPV SUSY with $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ decay.

png pdf
Figure 1-h:
Feynman diagram for dRPV SUSY with $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ decay.

png pdf
Figure 2:
Left : the NI-veto map based on the NI vertex reconstruction in the 2017 and 2018 data collected by the CMS detector, the map corresponds to the geometry of the CMS pixel detector used in 2017-2018 data taking [79]. The structures of the different pixel layers can be clearly seen. Right : the efficiency for a given vertex candidate to pass the NI-veto as a function of radius $r$.

png pdf
Figure 2-a:
The NI-veto map based on the NI vertex reconstruction in the 2017 and 2018 data collected by the CMS detector, the map corresponds to the geometry of the CMS pixel detector used in 2017-2018 data taking [79]. The structures of the different pixel layers can be clearly seen.

png pdf
Figure 2-b:
The efficiency for a given vertex candidate to pass the NI-veto as a function of radius $r$.

png pdf
Figure 3:
The distributions of the vertex track multiplicity (upper left), vertex $L_{xy}$ significance (upper right), cluster RMS (lower left), and the magnitude of the signed ${\mathrm {Sig[IP_{2D}]}}$ sum $ {| \kappa |}$ (lower right), for data, simulated QCD multijet events, and simulated signal events. Data and simulated events are selected with the displaced-jet triggers and with the offline ${H_{\mathrm {T}}}$, jets ${p_{\mathrm {T}}}$, and $\eta $ selections applied. For a given event, if there is more than one SV candidate being reconstructed, the one with the largest vertex track multiplicity is chosen. If the track multiplicities are the same, the one with the smallest $\chi ^{2}/\mathrm {n_{dof}}$ is chosen. The lower panels show the ratios between the data and the simulated QCD multijet events. The blue shaded error bands and vertical bars represent the statistical uncertainties. Three benchmark signal distributions are shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and varying $c\tau _{0}$. For visualization purposes, each signal process is given a cross section that yields 10$^{6}$ events produced in the analyzed data sample.

png pdf
Figure 3-a:
Distribution of the vertex track multiplicity, for data, simulated QCD multijet events, and simulated signal events. Data and simulated events are selected with the displaced-jet triggers and with the offline ${H_{\mathrm {T}}}$, jets ${p_{\mathrm {T}}}$, and $\eta $ selections applied. For a given event, if there is more than one SV candidate being reconstructed, the one with the largest vertex track multiplicity is chosen. If the track multiplicities are the same, the one with the smallest $\chi ^{2}/\mathrm {n_{dof}}$ is chosen. The lower panel shows the ratios between the data and the simulated QCD multijet events. The blue shaded error bands and vertical bars represent the statistical uncertainties. Three benchmark signal distributions are shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and varying $c\tau _{0}$. For visualization purposes, each signal process is given a cross section that yields 10$^{6}$ events produced in the analyzed data sample.

png pdf
Figure 3-b:
Distribution of the vertex $L_{xy}$ significance, for data, simulated QCD multijet events, and simulated signal events. Data and simulated events are selected with the displaced-jet triggers and with the offline ${H_{\mathrm {T}}}$, jets ${p_{\mathrm {T}}}$, and $\eta $ selections applied. For a given event, if there is more than one SV candidate being reconstructed, the one with the largest vertex track multiplicity is chosen. If the track multiplicities are the same, the one with the smallest $\chi ^{2}/\mathrm {n_{dof}}$ is chosen. The lower panel shows the ratios between the data and the simulated QCD multijet events. The blue shaded error bands and vertical bars represent the statistical uncertainties. Three benchmark signal distributions are shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and varying $c\tau _{0}$. For visualization purposes, each signal process is given a cross section that yields 10$^{6}$ events produced in the analyzed data sample.

png pdf
Figure 3-c:
Distribution of the cluster RMS, for data, simulated QCD multijet events, and simulated signal events. Data and simulated events are selected with the displaced-jet triggers and with the offline ${H_{\mathrm {T}}}$, jets ${p_{\mathrm {T}}}$, and $\eta $ selections applied. For a given event, if there is more than one SV candidate being reconstructed, the one with the largest vertex track multiplicity is chosen. If the track multiplicities are the same, the one with the smallest $\chi ^{2}/\mathrm {n_{dof}}$ is chosen. The lower panel shows the ratios between the data and the simulated QCD multijet events. The blue shaded error bands and vertical bars represent the statistical uncertainties. Three benchmark signal distributions are shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and varying $c\tau _{0}$. For visualization purposes, each signal process is given a cross section that yields 10$^{6}$ events produced in the analyzed data sample.

png pdf
Figure 3-d:
Distribution of the magnitude of the signed ${\mathrm {Sig[IP_{2D}]}}$ sum $ {| \kappa |}$, for data, simulated QCD multijet events, and simulated signal events. Data and simulated events are selected with the displaced-jet triggers and with the offline ${H_{\mathrm {T}}}$, jets ${p_{\mathrm {T}}}$, and $\eta $ selections applied. For a given event, if there is more than one SV candidate being reconstructed, the one with the largest vertex track multiplicity is chosen. If the track multiplicities are the same, the one with the smallest $\chi ^{2}/\mathrm {n_{dof}}$ is chosen. The lower panel shows the ratios between the data and the simulated QCD multijet events. The blue shaded error bands and vertical bars represent the statistical uncertainties. Three benchmark signal distributions are shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and varying $c\tau _{0}$. For visualization purposes, each signal process is given a cross section that yields 10$^{6}$ events produced in the analyzed data sample.

png pdf
Figure 4:
The distributions of the GBDT output score for data, simulated QCD multijet events, and simulated signal events. Data and simulated events are selected with the displaced-jet triggers and with the offline ${H_{\mathrm {T}}}$, jets ${p_{\mathrm {T}}}$, and $\eta $ selections applied. For a given event, if there is more than one SV candidate being reconstructed, the one with the largest vertex track multiplicity is chosen. If the track multiplicities are the same, the one with the smallest $\chi ^{2}/\mathrm {n_{dof}}$ is chosen. The lower panel shows the ratio between the data and the simulated QCD multijet events. The blue shaded error bands and vertical bars represent the statistical uncertainties. Three benchmark signal distributions are shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and varying $c\tau _{0}$. For visualization purposes, each signal process is given a cross section corresponding to 10$^{6}$ events produced in the analyzed data sample. The signal events shown in this plot are not used in the GBDT training.

png pdf
Figure 5:
The predicted background yields and the numbers of observed events in the control region, for different bins of the GBDT scores. The background predictions in different bins are correlated, since the events that are used for background predictions in lower bins are also used in the background predictions in higher bins. The error bands for the predictions represent statistical uncertainties and systematic uncertainties added in quadrature. The error bars for the observed events represent statistical uncertainties, assuming Poisson statistics.

png pdf
Figure 6:
The predicted background yields and the number of observed events for the data in the signal region, with $N_{\mathrm {tracks}}^{\mathrm {3D}}$ smaller than 3 for both jets, shown for different bins of the GBDT scores. The background predictions in different bins are correlated, since the events that are used for background predictions in lower bins are also used in the background predictions in higher bins. For comparison, three benchmark signal points are also shown (dashed lines) for the jet-jet model with $m_{\mathrm {X}} = $ 300 GeV and different lifetimes. For visualization purposes, each signal process is given a cross section that yields 100 events produced in the analyzed data sample.

png pdf
Figure 7:
The 95% CL upper limits on the pair production cross section of the LLP $\mathrm {X}$, where a 100% branching fraction for $\mathrm {X}$ to decay to a quark-antiquark pair is assumed. Left : the upper limits as functions of $c\tau _{0}$ for different masses. Right : the upper limits as functions of the particle mass for different $c\tau _{0}$. The solid (dashed) curves show the observed (median expected) limits. The shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 7-a:
The 95% CL upper limits on the pair production cross section of the LLP $\mathrm {X}$, where a 100% branching fraction for $\mathrm {X}$ to decay to a quark-antiquark pair is assumed: the upper limits as functions of $c\tau _{0}$ for different masses. The solid (dashed) curves show the observed (median expected) limits. The shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 7-b:
The 95% CL upper limits on the pair production cross section of the LLP $\mathrm {X}$, where a 100% branching fraction for $\mathrm {X}$ to decay to a quark-antiquark pair is assumed: the upper limits as functions of the particle mass for different $c\tau _{0}$. The solid (dashed) curves show the observed (median expected) limits. The shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 8:
The expected and observed 95$%$ CL upper limits on the branching fraction of the SM-like Higgs boson to decay to two long-lived scalar particles, assuming the gluon-gluon fusion Higgs boson production cross section of 49 pb at 13 TeV with $m_{\mathrm{H}} = $ 125 GeV, shown at different masses and $c\tau _{0}$ for the scalar particle $\mathrm {S}$. Left : the upper limits when each scalar particle decays to a down quark-antiquark pair. Right : the upper limits when each scalar particle decays to a bottom quark-antiquark pair. The solid (dashed) curves represent the observed (median expected) limits. The shaded bands represent the regions containing 68% of the distributions of the expected limits under the background-only hypothesis.

png pdf
Figure 8-a:
The expected and observed 95$%$ CL upper limits on the branching fraction of the SM-like Higgs boson to decay to two long-lived scalar particles, assuming the gluon-gluon fusion Higgs boson production cross section of 49 pb at 13 TeV with $m_{\mathrm{H}} = $ 125 GeV, shown at different masses and $c\tau _{0}$ for the scalar particle $\mathrm {S}$: the upper limits when each scalar particle decays to a down quark-antiquark pair. The solid (dashed) curves represent the observed (median expected) limits. The shaded bands represent the regions containing 68% of the distributions of the expected limits under the background-only hypothesis.

png pdf
Figure 8-b:
The expected and observed 95$%$ CL upper limits on the branching fraction of the SM-like Higgs boson to decay to two long-lived scalar particles, assuming the gluon-gluon fusion Higgs boson production cross section of 49 pb at 13 TeV with $m_{\mathrm{H}} = $ 125 GeV, shown at different masses and $c\tau _{0}$ for the scalar particle $\mathrm {S}$: the upper limits when each scalar particle decays to a bottom quark-antiquark pair. The solid (dashed) curves represent the observed (median expected) limits. The shaded bands represent the regions containing 68% of the distributions of the expected limits under the background-only hypothesis.

png pdf
Figure 9:
Left : the 95% CL upper limits on the pair production cross section for the long-lived gluinos with $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, where a 100% branching fraction for ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ decays is assumed. The NNLO$_{approx}$+NNLL gluino pair production cross sections for $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. The observed limits from the CMS search utilizing the timing capabilities of the ECAL system [48] are also shown for comparison. Right : the 95% CL upper limits on the pair production cross section for the ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ model as a function of the mean proper decay length $c\tau _{0}$ and the gluino mass $m_{{\mathrm{\widetilde{g}}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limit on the gluino mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 9-a:
The 95% CL upper limits on the pair production cross section for the long-lived gluinos with $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, where a 100% branching fraction for ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ decays is assumed. The NNLO$_{approx}$+NNLL gluino pair production cross sections for $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. The observed limits from the CMS search utilizing the timing capabilities of the ECAL system [48] are also shown for comparison.

png pdf
Figure 9-b:
The 95% CL upper limits on the pair production cross section for the ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ model as a function of the mean proper decay length $c\tau _{0}$ and the gluino mass $m_{{\mathrm{\widetilde{g}}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limit on the gluino mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 10:
Left : the 95% CL upper limits on the pair production cross section for the long-lived gluinos with $m_{{\mathrm{\widetilde{g}}}} = $ 2400 GeV and 1600 GeV, where a 100% branching fraction for ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ decays is assumed. The NNLO$_{approx}$+NNLL gluino pair production cross sections for $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. Right : the 95% CL limits on the pair production cross section for the ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ model as a function of the mean proper decay length $c\tau _{0}$ and the gluino mass $m_{{\mathrm{\widetilde{g}}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the gluino mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 10-a:
The 95% CL upper limits on the pair production cross section for the long-lived gluinos with $m_{{\mathrm{\widetilde{g}}}} = $ 2400 GeV and 1600 GeV, where a 100% branching fraction for ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ decays is assumed. The NNLO$_{approx}$+NNLL gluino pair production cross sections for $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 10-b:
The 95% CL limits on the pair production cross section for the ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ model as a function of the mean proper decay length $c\tau _{0}$ and the gluino mass $m_{{\mathrm{\widetilde{g}}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the gluino mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 11:
Left : the 95% CL upper limits on the pair production cross section for the long-lived gluinos with $m_{{\mathrm{\widetilde{g}}}} = $ 2000 GeV and 1400 GeV, where a 100% branching fraction for ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ decays is assumed. The NNLO$_{approx}$+NNLL gluino pair production cross sections for $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. Right : the 95% CL limits on the pair production cross section for the ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ model as a function of the mean proper decay length $c\tau _{0}$ and the gluino mass $m_{{\mathrm{\widetilde{g}}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the gluino mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 11-a:
The 95% CL upper limits on the pair production cross section for the long-lived gluinos with $m_{{\mathrm{\widetilde{g}}}} = $ 2000 GeV and 1400 GeV, where a 100% branching fraction for ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ decays is assumed. The NNLO$_{approx}$+NNLL gluino pair production cross sections for $m_{{\mathrm{\widetilde{g}}}}= $ 2400 and 1600 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 11-b:
The 95% CL limits on the pair production cross section for the ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ model as a function of the mean proper decay length $c\tau _{0}$ and the gluino mass $m_{{\mathrm{\widetilde{g}}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the gluino mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 12:
Left : the 95% CL upper limits on the pair production cross section for the long-lived top squarks with $m_{\tilde{\mathrm{t}}} = $ 1600 GeV and 800 GeV, where a 100% branching fraction for $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ decays is assumed, with equal branching fractions for e, $\mu $, and $\tau $. The NNLO$_{approx}$+NNLL top squark pair production cross sections for $m_{\tilde{\mathrm{t}}}=1600$ and 1000 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. Right : the 95% CL limits on the pair production cross section for the $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ model as a function of the mean proper decay length $c\tau _{0}$ and the top squark mass $m_{\tilde{\mathrm{t}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the top squark mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 12-a:
The 95% CL upper limits on the pair production cross section for the long-lived top squarks with $m_{\tilde{\mathrm{t}}} = $ 1600 GeV and 800 GeV, where a 100% branching fraction for $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ decays is assumed, with equal branching fractions for e, $\mu $, and $\tau $. The NNLO$_{approx}$+NNLL top squark pair production cross sections for $m_{\tilde{\mathrm{t}}}=1600$ and 1000 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 12-b:
The 95% CL limits on the pair production cross section for the $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ model as a function of the mean proper decay length $c\tau _{0}$ and the top squark mass $m_{\tilde{\mathrm{t}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the top squark mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 13:
Left : the 95% CL upper limits on the pair production cross section for the long-lived top squarks with $m_{\tilde{\mathrm{t}}} = $ 1600 GeV and 800 GeV, where a 100% branching fraction for $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ decays is assumed, with equal branching fractions for e, $\mu $, and $\tau $. The NNLO$_{approx}$+NNLL top squark pair production cross sections for $m_{\tilde{\mathrm{t}}}=1600$ and 1000 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. Right : the 95% CL limits on the pair production cross section for the $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ model as a function of the mean proper decay length $c\tau _{0}$ and the top squark mass $m_{\tilde{\mathrm{t}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the top squark mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 13-a:
The 95% CL upper limits on the pair production cross section for the long-lived top squarks with $m_{\tilde{\mathrm{t}}} = $ 1600 GeV and 800 GeV, where a 100% branching fraction for $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ decays is assumed, with equal branching fractions for e, $\mu $, and $\tau $. The NNLO$_{approx}$+NNLL top squark pair production cross sections for $m_{\tilde{\mathrm{t}}}=1600$ and 1000 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 13-b:
The 95% CL limits on the pair production cross section for the $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ model as a function of the mean proper decay length $c\tau _{0}$ and the top squark mass $m_{\tilde{\mathrm{t}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the top squark mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 14:
Left : the 95% CL upper limits on the pair production cross section for the long-lived top squarks with $m_{\tilde{\mathrm{t}}} = $ 1600 GeV and 800 GeV, where a 100% branching fraction for $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ decays is assumed. The NNLO$_{approx}$+NNLL top squark pair production cross sections for $m_{\tilde{\mathrm{t}}}=1600$ and 1000 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis. Right : the 95% CL limits on the pair production cross section for the $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ model as a function of the mean proper decay length $c\tau _{0}$ and the top squark mass $m_{\tilde{\mathrm{t}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the top squark mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure 14-a:
The 95% CL upper limits on the pair production cross section for the long-lived top squarks with $m_{\tilde{\mathrm{t}}} = $ 1600 GeV and 800 GeV, where a 100% branching fraction for $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ decays is assumed. The NNLO$_{approx}$+NNLL top squark pair production cross sections for $m_{\tilde{\mathrm{t}}}=1600$ and 1000 GeV, as well as their variations due to the theoretical uncertainties, are shown as horizontal lines. The solid (dashed) curves show the observed (median expected) limits, and the shaded bands indicate the regions containing 68% of the distributions of the limits expected under the background-only hypothesis.

png pdf
Figure 14-b:
The 95% CL limits on the pair production cross section for the $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ model as a function of the mean proper decay length $c\tau _{0}$ and the top squark mass $m_{\tilde{\mathrm{t}}}$. The thick solid black (dashed red) curve shows the observed (median expected) 95% CL limits on the top squark mass as a function of $c\tau _{0}$, assuming the NNLO$_{approx}$+NNLL cross sections. The thin dashed red curves indicate the region containing 68% of the distribution of the limits expected under the background-only hypothesis. The thin solid black curves represent the change in the observed limit when the signal cross sections are varied according to their theoretical uncertainties.

png pdf
Figure A1:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ model (upper left), the ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ model (upper right), the ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ model (middle left), the $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ model (middle right), the $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ model (lower left), and the $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ model (lower right).

png pdf
Figure A1-a:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ model.

png pdf
Figure A1-b:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ model.

png pdf
Figure A1-c:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ model.

png pdf
Figure A1-d:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ model.

png pdf
Figure A1-e:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ model.

png pdf
Figure A1-f:
The signal efficiencies as functions of the long-lived particle mass and mean proper decay length in the 2017 and 2018 analysis, for the $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ model.
Tables

png pdf
Table 1:
Summary of the preselection criteria.

png pdf
Table 2:
The definitions of the different regions used in the background estimation.

png pdf
Table 3:
Summary of the systematic uncertainties in the signal yields.

png pdf
Table 4:
Event yields after different selection requirements have been applied for data collected in 2017 and 2018. Signal efficiencies for the jet-jet model with $m_{\mathrm {X}} = $ 1000 GeV and different $c\tau _{0}$ are also shown for comparison. Selection requirements are cumulative from the first row to the last.

png pdf
Table A1:
Signal efficiencies for the jet-jet model in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{\mathrm {X}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{\mathrm {X}}$. Uncertainties are statistical only.

png pdf
Table A2:
Signal efficiencies for the model where the SM-like Higgs boson decays to two long-lived scalar particles $\mathrm {S}$ in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and with $m_{\mathrm {S}} = $ 55 GeV. The long-lived scalar particle is assumed to decay to a down quark-antiquark pair ($\mathrm {S}\to \mathrm{d} \mathrm{\bar{d}} $). Selection requirements are cumulative from the first row to the last. Uncertainties are statistical only.

png pdf
Table A3:
Signal efficiencies for the model where the SM-like Higgs boson decays to two long-lived scalar particles $\mathrm {S}$ in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and with $m_{\mathrm {S}} = $ 55 GeV. The long-lived scalar particle is assumed to decay to a bottom quark-antiquark pair ($\mathrm {S}\to \mathrm{b} \mathrm{\bar{b}} $). Selection requirements are cumulative from the first row to the last. Uncertainties are statistical only.

png pdf
Table A4:
Signal efficiencies for the ${\mathrm{\widetilde{g}}} \to \mathrm{g} \tilde{\mathrm{G}} $ model in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{{\mathrm{\widetilde{g}}}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{{\mathrm{\widetilde{g}}}}$. Uncertainties are statistical only.

png pdf
Table A5:
Signal efficiencies for the ${\mathrm{\widetilde{g}}} \to \mathrm{q} \mathrm{\bar{q}} \tilde{\chi}^{0}_{1}$ model ($m_{\tilde{\chi}^{0}_{1}} = $ 100 GeV) in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{{\mathrm{\widetilde{g}}}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{{\mathrm{\widetilde{g}}}}$. Uncertainties are statistical only.

png pdf
Table A6:
Signal efficiencies for the ${\mathrm{\widetilde{g}}} \to \mathrm{t} \mathrm{b} \mathrm{s} $ model in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{{\mathrm{\widetilde{g}}}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{{\mathrm{\widetilde{g}}}}$. Uncertainties are statistical only.

png pdf
Table A7:
Signal efficiencies for the $\tilde{\mathrm{t}} \to \mathrm{b} \ell $ model in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{\tilde{\mathrm{t}}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{\tilde{\mathrm{t}}}$. Uncertainties are statistical only.

png pdf
Table A8:
Signal efficiencies for the $\tilde{\mathrm{t}} \to \mathrm{d} \ell $ model in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{\tilde{\mathrm{t}}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{\tilde{\mathrm{t}}}$. Uncertainties are statistical only.

png pdf
Table A9:
Signal efficiencies for the $\tilde{\mathrm{t}} \to \mathrm{\bar{d}} \mathrm{\bar{d}} $ model in the 2017 and 2018 analysis at different mean proper decay lengths $c\tau _{0}$ and different masses $m_{\tilde{\mathrm{t}}}$. Selection requirements are cumulative from the first row to the last for each value of $m_{\tilde{\mathrm{t}}}$. Uncertainties are statistical only.
Summary
A search has been presented for long-lived particles decaying to jets, using proton-proton collision data collected with the CMS experiment at a center-of-mass energy of 13 TeV in 2017 and 2018. The results are combined with those of a previous CMS search for displaced jets using proton-proton collision data from 2016, accumulating to a total integrated luminosity of 132 fb$^{-1}$. After all selections, one event is observed in the data collected in 2017 and 2018, which is consistent with the predicted background yield. The search is designed to be model independent, and is sensitive to a large number of models predicting displaced-jets signatures with different final-state topologies.

The best current limits are set on a variety of models that have long-lived particles with mean proper decay lengths between 1 mm and 10 m. All limits are computed at the 95% confidence level. For a simplified model where pair-produced long-lived neutral particles decay to quark-antiquark pairs, pair production cross sections larger than 0.07 fb are excluded for mean proper decay lengths between 2 and 250 mm at high mass ($m_{\mathrm{X}} > $ 500 GeV). For a model where the standard model-like Higgs boson decays to two long-lived scalar particles and each long-lived scalar particle decays to a down (bottom) quark-antiquark pair, branching fractions for the exotic Higgs boson decay larger than 1% (10%) are excluded for mean proper decay lengths between 1 and 340 mm (530 mm) when the scalar particle mass is larger than 40 GeV. For a supersymmetric (SUSY) model in the general gauge mediation scenario, where the long-lived gluino decays to a gluon and a lightest SUSY particle, gluino masses up to 2450 GeV are excluded for mean proper decay lengths between 6 and 550 mm. For another SUSY model in the mini-split scenario, where the long-lived gluino can decay to a quark-antiquark pair and the lightest neutralino, gluino masses up to 2500 GeV are excluded for mean proper decay lengths between 7 and 360 mm. An $R$-parity violating (RPV) SUSY model is also tested, where the long-lived gluino can decay to top, bottom, and strange antiquarks, and gluino masses up to 2500 GeV are excluded for mean proper decay lengths between 3 and 1000 mm. Another RPV SUSY model is studied, where the long-lived top squark can decay to a bottom quark and a charged lepton, and top squark masses up to 1600 GeV are excluded for mean proper decay lengths between 5 and 240 mm. For an RPV SUSY model, where the long-lived top squark can decay to a down quark and a charged lepton, top squark masses up to 1600 GeV are excluded for mean proper decay lengths between 3 and 360 mm. Finally, for a dynamical-RPV SUSY model, where the long-lived top squark can decay to two down antiquarks, top squark masses up to 1600 GeV are excluded for mean proper decay lengths between 2 and 1320 mm. These are the most stringent limits to date on these models.
References
1 N. Arkani-Hamed and S. Dimopoulos Supersymmetric unification without low energy supersymmetry and signatures for fine-tuning at the LHC JHEP 06 (2005) 073 hep-th/0405159
2 G. F. Giudice and A. Romanino Split supersymmetry NPB 699 (2004) 65 hep-ph/0406088
3 J. L. Hewett, B. Lillie, M. Masip, and T. G. Rizzo Signatures of long-lived gluinos in split supersymmetry JHEP 09 (2004) 070 hep-ph/0408248
4 N. Arkani-Hamed, S. Dimopoulos, G. F. Giudice, and A. Romanino Aspects of split supersymmetry NPB 709 (2005) 3 hep-ph/0409232
5 P. Gambino, G. F. Giudice, and P. Slavich Gluino decays in split supersymmetry NPB 726 (2005) 35 hep-ph/0506214
6 A. Arvanitaki, N. Craig, S. Dimopoulos, and G. Villadoro Mini-split JHEP 02 (2013) 126 1210.0555
7 N. Arkani-Hamed et al. Simply unnatural supersymmetry 1212.6971
8 P. Fayet Supergauge invariant extension of the Higgs mechanism and a model for the electron and its neutrino NPB 90 (1975) 104
9 G. R. Farrar and P. Fayet Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry PLB 76 (1978) 575
10 S. Weinberg Supersymmetry at ordinary energies. 1. masses and conservation laws PRD 26 (1982) 287
11 L. J. Hall and M. Suzuki Explicit $ R $-parity breaking in supersymmetric models NPB 231 (1984) 419
12 R. Barbier et al. $ R $-parity violating supersymmetry PR 420 (2005) 1 hep-ph/0406039
13 G. F. Giudice and R. Rattazzi Theories with gauge mediated supersymmetry breaking PR 322 (1999) 419 hep-ph/9801271
14 P. Meade, N. Seiberg, and D. Shih General gauge mediation Prog. Theor. Phys. Suppl. 177 (2009) 143 0801.3278
15 M. Buican, P. Meade, N. Seiberg, and D. Shih Exploring general gauge mediation JHEP 03 (2009) 016 0812.3668
16 J. Fan, M. Reece, and J. T. Ruderman Stealth supersymmetry JHEP 11 (2011) 012 1105.5135
17 J. Fan, M. Reece, and J. T. Ruderman A stealth supersymmetry sampler JHEP 07 (2012) 196 1201.4875
18 M. J. Strassler and K. M. Zurek Echoes of a hidden valley at hadron colliders PLB 651 (2007) 374 hep-ph/0604261
19 M. J. Strassler and K. M. Zurek Discovering the Higgs through highly-displaced vertices PLB 661 (2008) 263 hep-ph/0605193
20 T. Han, Z. Si, K. M. Zurek, and M. J. Strassler Phenomenology of hidden valleys at hadron colliders JHEP 07 (2008) 008 0712.2041
21 D. E. Kaplan, M. A. Luty, and K. M. Zurek Asymmetric dark matter PRD 79 (2009) 115016 0901.4117
22 L. J. Hall, K. Jedamzik, J. March-Russell, and S. M. West Freeze-in production of FIMP dark matter JHEP 03 (2010) 080 0911.1120
23 I.-W. Kim and K. M. Zurek Flavor and collider signatures of asymmetric dark matter PRD 89 (2014) 035008 1310.2617
24 Y. Cui and B. Shuve Probing baryogenesis with displaced vertices at the LHC JHEP 02 (2015) 049 1409.6729
25 R. T. Co, F. D'Eramo, L. J. Hall, and D. Pappadopulo Freeze-In dark matter with displaced signatures at colliders JCAP 12 (2015) 024 1506.07532
26 L. Calibbi, L. Lopez-Honorez, S. Lowette, and A. Mariotti Singlet-Doublet dark matter freeze-in: LHC displaced signatures versus cosmology JHEP 09 (2018) 037 1805.04423
27 Y. Cui, L. Randall, and B. Shuve A WIMPy baryogenesis miracle JHEP 04 (2012) 075 1112.2704
28 Y. Cui and R. Sundrum Baryogenesis for weakly interacting massive particles PRD 87 (2013) 116013 1212.2973
29 A. Atre, T. Han, S. Pascoli, and B. Zhang The search for heavy Majorana neutrinos JHEP 05 (2009) 030 0901.3589
30 M. Drewes The phenomenology of right handed neutrinos Int. J. Mod. Phys. E 22 (2013) 1330019 1303.6912
31 F. F. Deppisch, P. S. Bhupal Dev, and A. Pilaftsis Neutrinos and collider physics New J. Phys. 17 (2015) 075019 1502.06541
32 Y. Cai, T. Han, T. Li, and R. Ruiz Lepton number violation: seesaw models and their collider tests Front. Phys. 6 (2018) 40 1711.02180
33 Z. Chacko, H.-S. Goh, and R. Harnik Natural electroweak breaking from a mirror symmetry PRL 96 (2006) 231802 hep-ph/0506256
34 H. Cai, H.-C. Cheng, and J. Terning A quirky little Higgs model JHEP 05 (2009) 045 0812.0843
35 N. Craig, S. Knapen, and P. Longhi Neutral naturalness from orbifold Higgs models PRL 114 (2015) 061803 1410.6808
36 N. Craig, A. Katz, M. Strassler, and R. Sundrum Naturalness in the dark at the LHC JHEP 07 (2015) 105 1501.05310
37 D. Curtin and C. B. Verhaaren Discovering uncolored naturalness in exotic Higgs decays JHEP 12 (2015) 072 1506.06141
38 S. Alipour-Fard et al. The second Higgs at the lifetime frontier JHEP 07 (2020) 029 1812.09315
39 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
40 CMS Collaboration Search for long-lived particles decaying into displaced jets in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PRD 99 (2019) 032011 CMS-EXO-18-007
1811.07991
41 ATLAS Collaboration Search for long-lived, massive particles in events with displaced vertices and missing transverse momentum in $ \sqrt{s} = $ 13 TeV pp collisions with the ATLAS detector PRD 97 (2018) 052012 1710.04901
42 ATLAS Collaboration Search for long-lived particles produced in $ pp $ collisions at $ \sqrt{s}= $ 13 TeV that decay into displaced hadronic jets in the ATLAS muon spectrometer PRD 99 (2019) 052005 1811.07370
43 ATLAS Collaboration Search for long-lived neutral particles in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV that decay into displaced hadronic jets in the ATLAS calorimeter EPJC 79 (2019) 481 1902.03094
44 ATLAS Collaboration Search for long-lived, massive particles in events with a displaced vertex and a muon with large impact parameter in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector PRD 102 (2020) 032006 2003.11956
45 ATLAS Collaboration Search for long-lived neutral particles produced in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV decaying into displaced hadronic jets in the ATLAS inner detector and muon spectrometer PRD 101 (2020) 052013 1911.12575
46 CMS Collaboration Search for new long-lived particles at $ \sqrt{s} = $ 13 TeV PLB 780 (2018) 432 CMS-EXO-16-003
1711.09120
47 CMS Collaboration Search for long-lived particles with displaced vertices in multijet events in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PRD 98 (2018) 092011 CMS-EXO-17-018
1808.03078
48 CMS Collaboration Search for long-lived particles using nonprompt jets and missing transverse momentum with proton-proton collisions at $ \sqrt{s} = $ 13 TeV PLB 797 (2019) 134876 CMS-EXO-19-001
1906.06441
49 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
50 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
51 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ {k_{\mathrm{T}}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
52 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
53 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
54 CMS Collaboration Technical proposal for the Phase-II upgrade of the CMS detector CMS-PAS-TDR-15-002 CMS-PAS-TDR-15-002
55 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
56 T. Sjostrand et al. An introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
57 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
58 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
59 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
60 LHC Higgs Cross Section Working Group Handbook of LHC Higgs cross sections: 4. Deciphering the nature of the Higgs sector CERN-2017-002-M 1610.07922
61 G. Burdman, Z. Chacko, H.-S. Goh, and R. Harnik Folded supersymmetry and the LEP paradox JHEP 02 (2007) 009 hep-ph/0609152
62 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
63 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
64 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
65 E. Bagnaschi, G. Degrassi, P. Slavich, and A. Vicini Higgs production via gluon fusion in the POWHEG approach in the SM and in the MSSM JHEP 02 (2012) 088 1111.2854
66 Z. Liu and B. Tweedie The fate of long-lived superparticles with hadronic decays after LHC Run 1 JHEP 06 (2015) 042 1503.05923
67 C. Csaki, Y. Grossman, and B. Heidenreich Minimal flavor violation supersymmetry: a natural theory for $ R $-parity violation PRD 85 (2012) 095009 1111.1239
68 P. W. Graham, D. E. Kaplan, S. Rajendran, and P. Saraswat Displaced supersymmetry JHEP 07 (2012) 149 1204.6038
69 C. Csaki, E. Kuflik, and T. Volansky Dynamical $ R $-parity violation PRL 112 (2014) 131801 1309.5957
70 C. Csaki, E. Kuflik, O. Slone, and T. Volansky Models of dynamical $ R $-parity violation JHEP 06 (2015) 045 1502.03096
71 C. Csaki et al. Phenomenology of a long-lived LSP with $ R $-parity violation JHEP 08 (2015) 016 1505.00784
72 G. R. Farrar and P. Fayet Bounds on $ R $-hadron production from calorimetry experiments PLB 79 (1978) 442
73 M. Fairbairn et al. Stable massive particles at colliders PR 438 (2007) 1 hep-ph/0611040
74 R. Mackeprang and A. Rizzi Interactions of coloured heavy stable particles in matter EPJC 50 (2007) 353 hep-ph/0612161
75 R. Mackeprang and D. Milstead An updated description of heavy-hadron interactions in GEANT-4 EPJC 66 (2010) 493 0908.1868
76 GEANT4 Collaboration GEANT4--a simulation toolkit NIMA 506 (2003) 250
77 R. Fruhwirth, W. Waltenberger, and P. Vanlaer Adaptive vertex fitting JPG 34 (2007) N343
78 CMS Collaboration Precision measurement of the structure of the CMS inner tracking system using nuclear interactions JINST 13 (2018) P10034 CMS-TRK-17-001
1807.03289
79 CMS Collaboration CMS technical design report for the pixel detector upgrade CERN-LHCC-2012-016, CMS-TDR-011
80 S. C. Johnson Hierarchical clustering schemes Psychometrika 32 (1967) 241
81 Y. Freund and R. E. Schapire Experiments with a new boosting algorithm in Proceedings of the Thirteenth International Conference on International Conference on Machine Learning, ICML'96, San Francisco (1996)
82 J. Friedman, T. Hastie, and R. Tibshirani Additive logistic regression: a statistical view of boosting (with discussion and a rejoinder by the authors) Ann. Statist. 28 (2000) 337
83 J. H. Friedman Greedy function approximation: A gradient boosting machine. Ann. Statist. 29 (2001) 1189
84 H. Voss, A. Hocker, J. Stelzer, and F. Tegenfeldt TMVA, the toolkit for multivariate data analysis with ROOT in XIth International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT) 2007 physics/0703039
85 F. Pedregosa et al. Scikit-learn: machine learning in Python J. Mach. Learn. Res. 12 (2011) 2825 1201.0490
86 G. Punzi Sensitivity of searches for new signals and its optimization in Statistical problems in particle physics, astrophysics and cosmology. Proceedings, Conference, PHYSTAT 2003, Stanford physics/0308063
87 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s} = $ 13 TeV CMS-PAS-LUM-17-004 CMS-PAS-LUM-17-004
88 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} = $ 13 TeV CMS-PAS-LUM-18-002 CMS-PAS-LUM-18-002
89 S. Dulat et al. New parton distribution functions from a global analysis of quantum chromodynamics PRD 93 (2016) 033006 1506.07443
90 L. A. Harland-Lang, A. D. Martin, P. Motylinski, and R. S. Thorne Parton distributions in the LHC era: MMHT 2014 PDFs EPJC 75 (2015) 204 1412.3989
91 J. Butterworth et al. PDF4LHC recommendations for LHC Run II JPG 43 (2016) 023001 1510.03865
92 A. Buckley et al. LHAPDF6: parton density access in the LHC precision era EPJC 75 (2015) 132 1412.7420
93 T. Junk Confidence level computation for combining searches with small statistics NIMA 434 (1999) 435 hep-ex/9902006
94 A. L. Read Presentation of search results: The CLs technique JPG 28 (2002) 2693
95 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
96 ATLAS and CMS Collaborations, The LHC Higgs Combination group Procedure for the LHC Higgs boson search combination in Summer 2011 CMS-NOTE-2011-005
97 W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas Squark and gluino production at hadron colliders NPB 492 (1997) 51 hep-ph/9610490
98 A. Kulesza and L. Motyka Threshold resummation for squark-antisquark and gluino-pair production at the LHC PRL 102 (2009) 111802 0807.2405
99 A. Kulesza and L. Motyka Soft gluon resummation for the production of gluino-gluino and squark-antisquark pairs at the LHC PRD 80 (2009) 095004 0905.4749
100 W. Beenakker et al. Squark and gluino hadroproduction Int. J. Mod. Phys. A 26 (2011) 2637 1105.1110
101 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
102 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
Compact Muon Solenoid
LHC, CERN