CMSPASEXO20015  
Search for longlived particles decaying in the CMS endcap muon system in protonproton collisions at $\sqrt{s}= $ 13 TeV  
CMS Collaboration  
May 2021  
Abstract: A search for longlived particles (LLPs) produced in decays of standard model (SM) Higgs bosons in 137 fb$^{1}$ of protonproton collisions at $\sqrt{s}= $ 13 TeV recorded by the CMS experiment during 20162018 is presented. The search employs a novel technique to reconstruct hadronic decays of LLPs in the endcap muon system. The search is sensitive to a broad range of LLP decay modes including $\tau^{}\tau^{+}$, LLP masses as low as a few GeV, and is largely modelindependent. No excess of events above the SM background is observed and stringent limits on the Higgs boson ($\mathrm{h}^0$) branching fraction to LLPs ($\mathrm{S}$) are obtained, particularly for proper decay lengths greater than a few meters. This search result represents the most stringent limits on the branching fraction $\mathcal{B}(\mathrm{h}^0 \rightarrow\mathrm{S}\mathrm{S})$ for proper decay lengths greater than 640 m for S masses between 740 GeV.  
Links:
CDS record (PDF) ;
inSPIRE record ;
CADI line (restricted) ;
These preliminary results are superseded in this paper, PRL 127 (2021) 261804. 
Figures  Summary  Additional Figures & Material  References  CMS Publications 

Instructions for reinterpretation can be found here 
Figures  
png pdf 
Figure 1:
The signal efficiency of the cluster reconstruction, veto, and identification selections as a function of the simulated $r$ and $z$ decay positions of S decaying to $\mathrm{b} {}\mathrm{\bar{b}} $, for a mass of 15 GeV and a uniform mixture of $c\tau $ between 1 and 10 m. The barrel and endcap muon stations are drawn as black boxes and labeled by their station names. Regions occupied by steel shielding are shaded in gray. 
png pdf 
Figure 2:
The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {\mathrm{d} \mathrm{\bar{d}}} $, and $c\tau = $ 1 m) and data distributions of $N_\mathrm {hits}$ (left) and $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\text{cluster})$ (right) in the search region. The $N_\mathrm {hits}$ plot includes events in bins C and D, while the $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\text{cluster})$ includes events in bins A and D. The last bin in the $N_\mathrm {hits}$ distributions include all events in the overflow. The bin width for the signal $N_\mathrm {hits}$ distribution is twice the size of the bin width for the data distribution. 
png pdf 
Figure 2a:
The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {\mathrm{d} \mathrm{\bar{d}}} $, and $c\tau = $ 1 m) and data distribution of $N_\mathrm {hits}$ $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\text{cluster})$ in the search region. The plot includes events in bins C and D. The last bin includes all events in the overflow. The bin width for the signal distribution is twice the size of the bin width for the data distribution. 
png pdf 
Figure 2b:
The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {\mathrm{d} \mathrm{\bar{d}}} $, and $c\tau = $ 1 m) and data distribution of $N_\mathrm {hits}$ $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\text{cluster})$ in the search region. The plot includes events in bins A and D. 
png pdf 
Figure 3:
The 95% CL expected (dotted curves) and observed (solid curves) limits on the branching fraction $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} $ as a function of $c\tau $ for the $ {\mathrm {S}} \rightarrow {\mathrm{d} \mathrm{\bar{d}}} $ (left) and $ {\mathrm {S}} \rightarrow \tau^{}\tau^{+} $ (right) decay modes. The exclusion limits are shown for four different mass hypotheses: 7 (green), 15 (black), 40 (red), and 55 (blue) GeV. 
png pdf 
Figure 3a:
The 95% CL expected (dotted curves) and observed (solid curves) limits on the branching fraction $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} $ as a function of $c\tau $ for the $ {\mathrm {S}} \rightarrow {\mathrm{d} \mathrm{\bar{d}}} $ decay mode. The exclusion limits are shown for four different mass hypotheses: 7 (green), 15 (black), 40 (red), and 55 (blue) GeV. 
png pdf 
Figure 3b:
The 95% CL expected (dotted curves) and observed (solid curves) limits on the branching fraction $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} $ as a function of $c\tau $ for the $ {\mathrm {S}} \rightarrow \tau^{}\tau^{+} $ decay mode. The exclusion limits are shown for four different mass hypotheses: 7 (green), 15 (black), 40 (red), and 55 (blue) GeV. 
Summary 
In summary, 137 fb$^{1}$ of protonproton collision data at $\sqrt{s} = $ 13 TeV recorded by the CMS experiment were used to conduct the first search for hadronically decaying LLPs using the endcap muon system. This unique detector signature makes the search sensitive to a broad range of LLP decay modes including TT, with LLP masses as low as a few GeV, and is largely modelindependent. The excellent shielding provided by the CMS detector helped suppress backgrounds to a low level, thus enabling a search for a single LLP decay. We observe no significant deviation from the SM background and set the most stringent limits on the branching fraction $\mathcal{B}(\mathrm{h^0 \to SS})$ for $c\tau$ greater than 6, 20, and 40 m, and $m_{\mathrm{S}}$ of 7, 15, and 40 GeV, respectively. For $c\tau$ above 100 m, this search outperforms the previous best limits [31] by a factor of 6 for $m_{\mathrm{S}} = $ 7 GeV and by a factor 2 for $m_{\mathrm{S}} \ge $ 15 GeV. 
Additional Figures  
png pdf 
Additional Figure 1:
The 95% CL expected (dotted curves) and observed (solid curves) upper limits on the branching fraction $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} $ as a function of $c\tau $ for the $ {\mathrm {S}} \rightarrow {{\mathrm {b}} {\overline {\mathrm {b}}}} $ decay mode. The exclusion limits are shown for three different mass hypotheses: 15 (black), 40 (red), and 55 (blue) GeV. 
png pdf 
Additional Figure 2:
The $N_\mathrm {hits}$ (left) and $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\mathrm {cluster})$ (right) distributions. The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {{\mathrm {d}} {\overline {\mathrm {d}}}} $, and $c\tau =$ 1 m) and data (black) distributions have $\Delta \phi < $ 0.75 ($N_\mathrm {hits} > $ 130) selections applied in the left (right) plot. The data distribution shown in magenta includes the inverted cut $\Delta \phi > $ 0.75 ($N_\mathrm {hits} < $ 130) applied in the left (right) plot. The last bin in the $N_\mathrm {hits}$ distributions includes all events in the overflow. 
png pdf 
Additional Figure 2a:
The $N_\mathrm {hits}$ distributions. The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {{\mathrm {d}} {\overline {\mathrm {d}}}} $, and $c\tau =$ 1 m) and data (black) distributions have $\Delta \phi < $ 0.75 selections applied. The data distribution shown in magenta includes the inverted cut $\Delta \phi > $ 0.75 applied in the left (right) plot. The last bin includes all events in the overflow. 
png pdf 
Additional Figure 2b:
The $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\mathrm {cluster})$ distribution. The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {{\mathrm {d}} {\overline {\mathrm {d}}}} $, and $c\tau =$ 1 m) and data (black) distributions have $N_\mathrm {hits} > $ 130 selection applied. The data distribution shown in magenta includes the inverted cut $N_\mathrm {hits} < $ 130 applied. 
png pdf 
Additional Figure 3:
The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {{\mathrm {d}} {\overline {\mathrm {d}}}} $, and $c\tau =$ 1 m), intime data ($N_\mathrm {hits} < $130), and outoftime data (inclusive in $N_\mathrm {hits}$) distributions of $N_\mathrm {hits}$ (left) and $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\mathrm {cluster})$ (right). The signal distributions include events from regions A, B, C, and D (as defined in the analysis summary), the intime data distributions include events from regions B and C, and the OOT data are outoftime clusters from previous bunch crossings, with the yield normalized to that of the intime distribution. The last bin in the $N_\mathrm {hits}$ distributions includes all events in the overflow. 
png pdf 
Additional Figure 3a:
The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {{\mathrm {d}} {\overline {\mathrm {d}}}} $, and $c\tau =$ 1 m), intime data ($N_\mathrm {hits} < $130), and outoftime data (inclusive in $N_\mathrm {hits}$) distributions of $N_\mathrm {hits}$ The signal distributions include events from regions A, B, C, and D (as defined in the analysis summary), the intime data distributions include events from regions B and C, and the OOT data are outoftime clusters from previous bunch crossings, with the yield normalized to that of the intime distribution. 
png pdf 
Additional Figure 3b:
The signal (assuming $ {\mathcal {B}({{\mathrm {h}^0} \rightarrow {\mathrm {S}} {\mathrm {S}}})} = $ 1%, $ {\mathrm {S}} \rightarrow {{\mathrm {d}} {\overline {\mathrm {d}}}} $, and $c\tau =$ 1 m), intime data ($N_\mathrm {hits} < $130), and outoftime data (inclusive in $N_\mathrm {hits}$) distributions of $\Delta \phi (\vec{p}_{\mathrm{T}}^{\,\text{miss}},\mathrm {cluster})$. The signal distributions include events from regions A, B, C, and D (as defined in the analysis summary), the intime data distributions include events from regions B and C, and the OOT data are outoftime clusters from previous bunch crossings, with the yield normalized to that of the intime distribution. 
png pdf 
Additional Figure 4:
The event display of a simulated signal event for $m_ {\mathrm {S}} = $ 40 GeV, $c\tau =$ 1 m, and $ {\mathrm {S}} \rightarrow {{\mathrm {b}} {\overline {\mathrm {b}}}} $ in $rz$plane (left) and $r\phi $plane (right). The LLP decayed at the position $x, y, z = $ $$364, 92, 796 cm and produced a CSC hit cluster with 711 hits (orange dots). The green lines represent the tracks. The yellow lines represent jets. The red arrow represents the MET direction. The red and blue cones represent the ECAL and HCAL energy deposits, respectively. 
png pdf 
Additional Figure 4a:
The event display of a simulated signal event for $m_ {\mathrm {S}} = $ 40 GeV, $c\tau =$ 1 m, and $ {\mathrm {S}} \rightarrow {{\mathrm {b}} {\overline {\mathrm {b}}}} $ in the $rz$plane. The LLP decayed at the position $x, y, z = $ $$364, 92, 796 cm and produced a CSC hit cluster with 711 hits (orange dots). The green lines represent the tracks. The yellow lines represent jets. The red arrow represents the MET direction. The red and blue cones represent the ECAL and HCAL energy deposits, respectively. 
png pdf 
Additional Figure 4b:
The event display of a simulated signal event for $m_ {\mathrm {S}} = $ 40 GeV, $c\tau =$ 1 m, and $ {\mathrm {S}} \rightarrow {{\mathrm {b}} {\overline {\mathrm {b}}}} $ in the $r\phi $plane. The LLP decayed at the position $x, y, z = $ $$364, 92, 796 cm and produced a CSC hit cluster with 711 hits (orange dots). The green lines represent the tracks. The yellow lines represent jets. The red arrow represents the MET direction. The red and blue cones represent the ECAL and HCAL energy deposits, respectively. 
png pdf 
Additional Figure 5:
Distributions of cluster time in data and in simulated $ {\mathrm {S}} \rightarrow {{\mathrm {b}} {\overline {\mathrm {b}}}} $ signal events. The signal distribution contains equal fractions of events with mass of 15, 40, and 55 GeV, each simulated with $c\tau $ uniformly distributed between 0.1 and 10 m. Both distributions are normalized to unit area and are required to pass the veto requirements as defined in the analysis summary. 
png pdf 
Additional Figure 6:
The clustering efficiency requiring $ > $50 hits (left) and the overall signal efficiency of the clustering requiring $ > $50 hits, veto, and identification selections (right), as defined in the analysis summary, as a function of the simulated $z$ decay positions of S decaying to $ {{\mathrm {b}} {\overline {\mathrm {b}}}} $, for a mass of 15 GeV and a uniform mixture of $c\tau $ between 1 and 10 m. The endcap muon stations are labeled by their station names. Regions occupied by steel shielding are shaded in gray. 
png pdf 
Additional Figure 6a:
The clustering efficiency requiring $ > $50 hits, as defined in the analysis summary, as a function of the simulated $z$ decay positions of S decaying to $ {{\mathrm {b}} {\overline {\mathrm {b}}}} $, for a mass of 15 GeV and a uniform mixture of $c\tau $ between 1 and 10 m. The endcap muon stations are labeled by their station names. Regions occupied by steel shielding are shaded in gray. 
png pdf 
Additional Figure 6b:
The overall signal efficiency of the clustering requiring $ > $50 hits, veto, and identification selections, as defined in the analysis summary, as a function of the simulated $z$ decay positions of S decaying to $ {{\mathrm {b}} {\overline {\mathrm {b}}}} $, for a mass of 15 GeV and a uniform mixture of $c\tau $ between 1 and 10 m. The endcap muon stations are labeled by their station names. Regions occupied by steel shielding are shaded in gray. 
png pdf 
Additional Figure 7:
The cluster efficiency in bins of hadronic and EM energy in region A (left) and B (right), estimated with LLPs decaying to ${{\tau}^{+} {\tau}^{}}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. Region A is defined as 391 cm $ < r < $ 695.5 cm and 400 cm $ < z < $ 671 cm. Region B is defined as 671 cm $ < z < $ 1100 cm, $r < $695.5 cm and $\eta  < $ 2. The cluster efficiency includes all clusterlevel selections described in the paper, except for the jet veto, time cut, and $\Delta \phi ({{\vec{p}}_{\mathrm {T}}^{\,\text {miss}}},\mathrm {cluster})$ cut. The full simulation signal yield prediction for samples with various LLP mass between 755 GeV, lifetime between 0.1100 m, and decay mode to $ {{\mathrm {d}} {\overline {\mathrm {d}}}} $ and $ {{\tau}^{+} {\tau}^{}} $ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively. The statistical uncertainty for each bin is documented in Additional Figure 7 of the HEPData record of this paper. 
png pdf 
Additional Figure 7a:
The cluster efficiency in bins of hadronic and EM energy in region A, estimated with LLPs decaying to ${{\tau}^{+} {\tau}^{}}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. Region A is defined as 391 cm $ < r < $ 695.5 cm and 400 cm $ < z < $ 671 cm. The cluster efficiency includes all clusterlevel selections described in the paper, except for the jet veto, time cut, and $\Delta \phi ({{\vec{p}}_{\mathrm {T}}^{\,\text {miss}}},\mathrm {cluster})$ cut. The full simulation signal yield prediction for samples with various LLP mass between 755 GeV, lifetime between 0.1100 m, and decay mode to $ {{\mathrm {d}} {\overline {\mathrm {d}}}} $ and $ {{\tau}^{+} {\tau}^{}} $ can be reproduced using this parameterization to within 35%. The statistical uncertainty for each bin is documented in Additional Figure 7 of the HEPData record of this paper. 
png pdf 
Additional Figure 7b:
The cluster efficiency in bins of hadronic and EM energy in region B, estimated with LLPs decaying to ${{\tau}^{+} {\tau}^{}}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. Region B is defined as 671 cm $ < z < $ 1100 cm, $r < $695.5 cm and $\eta  < $ 2. The cluster efficiency includes all clusterlevel selections described in the paper, except for the jet veto, time cut, and $\Delta \phi ({{\vec{p}}_{\mathrm {T}}^{\,\text {miss}}},\mathrm {cluster})$ cut. The full simulation signal yield prediction for samples with various LLP mass between 755 GeV, lifetime between 0.1100 m, and decay mode to $ {{\mathrm {d}} {\overline {\mathrm {d}}}} $ and $ {{\tau}^{+} {\tau}^{}} $ can be reproduced using this parameterization to within 20%. The statistical uncertainty for each bin is documented in Additional Figure 7 of the HEPData record of this paper. 
png pdf 
Additional Figure 8:
The efficiency of $N_{station} > $ 1 requirement in bins of hadronic energy in region B, estimated with LLPs decaying to ${{\tau}^{+} {\tau}^{}}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. This efficiency is independent of the EM energy. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. Region B is defined as 671 cm $ < z < $ 1100 cm, $r < $ 695.5 cm and $\eta  < $ 2. The efficiency is calculated with respect to clusters that pass all the clusterlevel cuts described in the paper, except for the jet veto, time cut, and $\Delta \phi ({{\vec{p}}_{\mathrm {T}}^{\,\text {miss}}},\mathrm {cluster})$ cut. The full simulation signal yield prediction for samples with various LLP mass between 755 GeV, lifetime between 0.1100 m, and decay mode to $ {{\mathrm {d}} {\overline {\mathrm {d}}}} $ and $ {{\tau}^{+} {\tau}^{}} $ can be reproduced using this parameterization to within 10%. 
Instructions for Reinterpretation 
We provide signal efficiency parameterization for the clusterlevel selections that allows for reproduction of the fullsimulation signal yield for various LLP masses (755 GeV), lifetimes (0.1  100 m) and decay modes (dd̅ and τ^{+}τ^{}). In order to recast this analysis, only the generator level LLP hadronic energy, EM energy, and decay position are needed. The following selection efficiencies are needed to account for all clusterlevel selections mentioned in the paper:
The LLP EM and hadronic energy is assigned by matching stable (status 1) GenParticles that are produced 0.1 m from the LLP decay vertex. Neutral pions, electrons, and photons are assigned as EM energy. All other particles, except for neutrinos and muons, are assigned as hadronic energy. Neutrinos and muons are ignored because they do not produce showers in the muon system. The LLP decay region is categorized into 2 regions. These 2 regions have qualitatively different behavior. Within each region, they have quantitatively similar behavior, so we will provide the efficiency parameterization for each region separately. Region A is defined as 391 cm < r < 695.5 cm and 400 cm < z < 671 cm. Region B is defined as 671 cm < z < 1100 cm, r < 695.5 cm, and η < 2. The fraction of LLPs that decay in each region are dependent on the LLP mass and cτ. However, the signal efficiency in B is much larger than A, so for the models considered, more than 90% of clusters passing all selections are from LLPs that decay in region B. The cluster efficiency is parameterized in bins of LLP hadronic energy and EM energy in each LLP decay region (Additional Figure 7). The efficiency includes clusterlevel selections mentioned in the paper, including the segment and rechit vetos, muon veto, time spread cut, and N_{hits} ≥ 130, except for the jet veto, time cut, and Δφ(p_{T}^{miss}, cluster) cut. These latter requirements are model dependent, but they can be calculated accurately for each model using generator level information. When recasting the analysis, these additional selections need to be implemented, to be consistent with applying all selections described in the paper. The full simulation signal yield prediction for samples with various LLP mass between 7  55 GeV, lifetime between 0.1  100 m, and decay mode to dd̅ and τ^{+}τ^{} can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively. The parameterization of the cutbased ID is provided in python file cut_based_id.py attached. The function uses the LLP decay position to predict the average station of the cluster (AvgStation function in the code) and the N_{station} > 1 efficiency parameterization (Additional Figure 8). As described in the paper, the cluster ID requirement applies different η cuts depending on the N_{station} and average station number. We need a parametrization of the efficiency of the N_{station} > 1 requirement and a transfer function that takes genlevel LLP decay position to RECOlevel cluster average station (only for clusters with N_{station} = 1 ). Since the entire region A is in an η region (η < 1.3) that passes all the η selections used in the cutbased ID (tightest cut is η < 1.6), the cutbased ID efficiency in this region is 100%. The focus of the parameterization will be on region B only. The efficiency of N_{station} > 1 in region B is provided in bins of LLP hadronic energy (Additional Figure 8) and the efficiency is independent of the LLP EM energy. The average station transfer function is provided as the AvgStation function in the code attached. The full simulation signal yield prediction for samples with various LLP mass between 7  55 GeV, lifetime between 0.1  100 m, and decay mode to dd̅ and τ^{+}τ^{} can be reproduced using this parameterization to within 10%.

References  
1  G. F. Giudice and A. Romanino  Split supersymmetry  NPB 699 (2004) 65  hepph/0406088 
2  J. L. Hewett, B. Lillie, M. Masip, and T. G. Rizzo  Signatures of longlived gluinos in split supersymmetry  JHEP 09 (2004) 070  hepph/0408248 
3  N. ArkaniHamed, S. Dimopoulos, G. F. Giudice, and A. Romanino  Aspects of split supersymmetry  NPB 709 (2005) 3  hepph/0409232 
4  P. Gambino, G. F. Giudice, and P. Slavich  Gluino decays in split supersymmetry  NPB 726 (2005) 35  hepph/0506214 
5  A. Arvanitaki, N. Craig, S. Dimopoulos, and G. Villadoro  Minisplit  JHEP 02 (2013) 126  1210.0555 
6  N. ArkaniHamed et al.  Simply unnatural supersymmetry  1212.6971  
7  P. Fayet  Supergauge invariant extension of the Higgs mechanism and a model for the electron and its neutrino  NPB 90 (1975) 104  
8  G. R. Farrar and P. Fayet  Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry  PLB 76 (1978) 575  
9  S. Weinberg  Supersymmetry at ordinary energies. Masses and conservation laws  PRD 26 (1982) 287  
10  R. Barbier et al.  $ R $parity violating supersymmetry  PR 420 (2005) 1  hepph/0406039 
11  G. F. Giudice and R. Rattazzi  Theories with gauge mediated supersymmetry breaking  PR 322 (1999) 419  hepph/9801271 
12  P. Meade, N. Seiberg, and D. Shih  General gauge mediation  Prog. Theor. Phys. Suppl. 177 (2009) 143  0801.3278 
13  M. Buican, P. Meade, N. Seiberg, and D. Shih  Exploring general gauge mediation  JHEP 03 (2009) 016  0812.3668 
14  J. Fan, M. Reece, and J. T. Ruderman  Stealth supersymmetry  JHEP 11 (2011) 012  1105.5135 
15  J. Fan, M. Reece, and J. T. Ruderman  A stealth supersymmetry sampler  JHEP 07 (2012) 196  1201.4875 
16  M. J. Strassler and K. M. Zurek  Echoes of a hidden valley at hadron colliders  PLB 651 (2007) 374379  hepph/0604261 
17  M. J. Strassler and K. M. Zurek  Discovering the Higgs through highlydisplaced vertices  PLB 661 (2008) 263  hepph/0605193 
18  T. Han, Z. Si, K. M. Zurek, and M. J. Strassler  Phenomenology of hidden valleys at hadron colliders  JHEP 07 (2008) 008  0712.2041 
19  Y. Cui, L. Randall, and B. Shuve  A WIMPy baryogenesis miracle  JHEP 04 (2012) 075  1112.2704 
20  Y. Cui and R. Sundrum  Baryogenesis for weakly interacting massive particles  PRD 87 (2013) 116013  1212.2973 
21  Y. Cui and B. Shuve  Probing baryogenesis with displaced vertices at the LHC  JHEP 02 (2015) 049  1409.6729 
22  D. Smith and N. Weiner  Inelastic dark matter  PRD 64 (2001) 043502  
23  Z. Chacko, H.S. Goh, and R. Harnik  Natural electroweak breaking from a mirror symmetry  PRL 96 (2006) 231802  hepph/0506256 
24  D. Curtin and C. B. Verhaaren  Discovering uncolored naturalness in exotic Higgs decays  JHEP 12 (2015) 072  1506.06141 
25  H.C. Cheng, S. Jung, E. Salvioni, and Y. Tsai  Exotic quarks in twin Higgs models  JHEP 03 (2016) 074  1512.02647 
26  N. Craig, A. Katz, M. Strassler, and R. Sundrum  Naturalness in the Dark at the LHC  JHEP 07 (2015) 105  1501.05310 
27  M. J. Strassler  On the Phenomenology of Hidden Valleys with Heavy Flavor  0806.2385  
28  J. E. Juknevich, D. Melnikov, and M. J. Strassler  A PureGlue Hidden Valley I. States and Decays  JHEP 07 (2009) 055  0903.0883 
29  CMS Collaboration  Search for longlived particles using displaced jets in protonproton collisions at $ \sqrt{s} = $ 13 TeV  Submitted to \emphPRD  CMSEXO19021 2012.01581 
30  ATLAS Collaboration  Search for longlived 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 
31  ATLAS Collaboration  Search for longlived 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 
32  CMS Collaboration  The CMS experiment at the CERN LHC  JINST 3 (2008) S08004  CMS00001 
33  CMS Collaboration  Particleflow reconstruction and global event description with the CMS detector  JINST 12 (2017) P10003  CMSPRF14001 1706.04965 
34  M. Cacciari and G. P. Salam  Dispelling the $ N^{3} $ myth for the $ {k_{\mathrm{T}}} $ jetfinder  PLB 641 (2006) 57  hepph/0512210 
35  M. Cacciari, G. P. Salam, and G. Soyez  The anti$ {k_{\mathrm{T}}} $ jet clustering algorithm  JHEP 04 (2008) 063  0802.1189 
36  M. Cacciari, G. P. Salam, and G. Soyez  FastJet user manual  EPJC 72 (2012) 1896  1111.6097 
37  CMS Collaboration  Technical proposal for the PhaseII upgrade of the Compact Muon Solenoid  CMSPASTDR15002  CMSPASTDR15002 
38  M. Cacciari and G. P. Salam  Pileup subtraction using jet areas  PLB 659 (2008) 119  0707.1378 
39  P. Nason  A new method for combining NLO QCD with shower Monte Carlo algorithms  JHEP 11 (2004) 040  hepph/0409146 
40  S. Frixione, P. Nason, and C. Oleari  Matching NLO QCD computations with parton shower simulations: the POWHEG method  JHEP 11 (2007) 070  0709.2092 
41  S. Alioli, P. Nason, C. Oleari, and E. Re  A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEGbox  JHEP 06 (2010) 043  1002.2581 
42  E. Re  Singletop Wt channel production matched with parton showers using the POWHEG method  EPJC 71 (2011) 1547  1009.2450 
43  T. Sjostrand et al.  An introduction to $ PYTHIA $ 8.2  CPC 191 (2015) 159  1410.3012 
44  CMS Collaboration  Event generator tunes obtained from underlying event and multiparton scattering measurements  EPJC 76 (2016) 155  CMSGEN14001 1512.00815 
45  CMS Collaboration  Extraction and validation of a new set of CMS PYTHIA8 tunes from underlyingevent measurements  EPJC 80 (2020) 4  CMSGEN17001 1903.12179 
46  NNPDF Collaboration  Parton distributions for the LHC Run II  JHEP 04 (2015) 040  1410.8849 
47  NNPDF Collaboration  Parton distributions from highprecision collider data  EPJC 77 (2017) 663  1706.00428 
48  GEANT4 Collaboration  GEANT4a simulation toolkit  NIMA 506 (2003) 250  
49  CMS Collaboration  The CMS trigger system  JINST 12 (2017) P01020  CMSTRG12001 1609.02366 
50  CMS Collaboration  Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC  Submitted to \emphJINST  CMSEGM17001 2012.06888 
51  CMS Collaboration  Performance of the CMS muon detector and muon reconstruction with protonproton collisions at $ \sqrt{s} = $ 13 TeV  JINST 13 (2018) P06015  CMSMUO16001 1804.04528 
52  CMS Collaboration  Performance of the reconstruction and identification of highmomentum muons in protonproton collisions at $ \sqrt{s} = $ 13 TeV  JINST 15 (2020) P02027  CMSMUO17001 1912.03516 
53  M. Ester, H.P. Kriegel, J. Sander, and X. Xu  A densitybased algorithm for discovering clusters in large spatial databases with noise  Data Mining and Knowledge Discovery (1996) 226  
54  CMS Collaboration  Missing transverse energy performance of the CMS detector  JINST 6 (2011) P09001  CMSJME10009 1106.5048 
55  T. Junk  Confidence level computation for combining searches with small statistics  NIMA 434 (1999) 435  
56  A. L. Read  Presentation of search results: the CLs technique  Journal of Physics G: Nuclear and Particle Physics 28 (2002) 2693  
57  The ATLAS Collaboration, The CMS Collaboration, The LHC Higgs Combination Group Collaboration  Procedure for the LHC Higgs boson search combination in Summer 2011  technical report, CERN, Geneva, Aug 
Compact Muon Solenoid LHC, CERN 