A search is performed for massive long-lived particles (LLPs) decaying semileptonically into a muon and two quarks. Two kinds of LLP production processes were considered. In the first, a Higgs-like boson with mass from 30 to 200\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document} is produced by gluon fusion and decays into two LLPs. The analysis covers LLP mass values from 10\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document} up to about one half the Higgs-like boson mass. The second LLP production mode is directly from quark interactions, with LLP masses from 10 to 90\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}. The LLP lifetimes considered range from 5 to 200 ps. This study uses LHCb data collected from proton-proton collisions at s=13\,TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s} = 13\text {\,TeV} $$\end{document}, corresponding to an integrated luminosity of 5.4\,fb-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,fb} ^{-1}$$\end{document}. No evidence of these long-lived states has been observed, and upper limits on the production cross-section times branching ratio have been set for each model considered.
publisher-imprint-nameSpringervolume-issue-count12issue-article-count94issue-toc-levels0issue-pricelist-year2022issue-copyright-holderThe Author(s)issue-copyright-year2022article-contains-esmNoarticle-numbering-styleContentOnlyarticle-registration-date-year2022article-registration-date-month3article-registration-date-day7article-toc-levels0toc-levels0volume-typeRegularjournal-productNonStandardArchiveJournalnumbering-styleContentOnlyarticle-grants-typeOpenChoicemetadata-grantOpenAccessabstract-grantOpenAccessbodypdf-grantOpenAccessbodyhtml-grantOpenAccessbibliography-grantOpenAccessesm-grantOpenAccessonline-firstfalsepdf-file-referenceBodyRef/PDF/10052_2022_Article_10186.pdfpdf-typeTypesettarget-typeOnlinePDFissue-typeRegulararticle-typeOriginalPaperjournal-subject-primaryPhysicsjournal-subject-secondaryElementary Particles, Quantum Field Theoryjournal-subject-secondaryNuclear Physics, Heavy Ions, Hadronsjournal-subject-secondaryQuantum Field Theories, String Theoryjournal-subject-secondaryMeasurement Science and Instrumentationjournal-subject-secondaryAstronomy, Astrophysics and Cosmologyjournal-subject-secondaryNuclear Energyjournal-subject-collectionPhysics and Astronomyopen-accesstrueIntroduction
Supersymmetry (SUSY) is one of the most popular extensions of the Standard Model (SM), which can solve the hierarchy problem, can unify the gauge couplings at the Planck scale and proposes dark matter candidates. The minimal supersymmetric extension of the Standard Model (MSSM) is the simplest phenomenologically viable realisation of SUSY [1, 2]. The present study addresses a subset of models featuring massive long-lived particles (LLPs) with a measurable flight distance [3, 4], decaying semileptonically. Long-lived particles decaying semileptonically with displaced jets composed of SM particles have been studied by the experiments at the LHC [5–9] Additional information on searches for LLPs at collider experiments can be found in Refs. [10–12].
This analysis uses proton-proton (pp) collision data at a centre-of-mass energy s\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s}$$\end{document} =13\,TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,TeV}$$\end{document} collected by the LHCb experiment at the LHC, corresponding to a total integrated luminosity of 5.4\,fb-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,fb} ^{-1}$$\end{document}. It extends the analysis of Ref. [9] on data collected at s=7\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s} =7$$\end{document} and 8\,TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,TeV}$$\end{document}. The adopted theoretical framework is inspired by the SUper GRAvity (mSUGRA) with R-parity violation (RPV) [13], in which the neutralino can decay into a muon and two quarks: χ~10→μ+qiqj(μ-qi¯qj¯)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1} \rightarrow \mu ^{+} q_{i} q_{j}( \mu ^{-} \bar{q_{i}} \bar{q_{j}})$$\end{document}. Neutralinos can be produced by a variety of processes. In this paper the analysis has been performed assuming the two mechanisms depicted in Fig. 1. In the first process, a Higgs-like particle, h0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$h^0$$\end{document}, is produced by gluon fusion and decays into two LLPs. The analysis covers h0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$h^0$$\end{document} masses from 30 to 200\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, LLP lifetimes from 5 to 200 ps and LLP mass values from 10\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document} up to about one half the h0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$h^0$$\end{document} mass. The second mode is a direct LLP production from quark interactions. The LLP lifetime range considered is from 5 to 200 ps and the mass range from 10 to 90\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}. The LLP lifetime range begins at 5 ps, well above the typical b-hadron lifetime, and extends up to 200 ps, where most of the vertices are still within the LHCb vertex locator (VELO). The mass range avoids the region of the SM b-quark states, but also takes into account the forward acceptance of the LHCb detector within which the decay products of relatively light LLPs can be efficiently detected.
The LLP signature is a displaced vertex made of charged particle tracks accompanied by an isolated muon with high transverse momentum with respect to the proton beam direction, pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document}. This study benefits from the excellent vertex reconstruction provided by the VELO, and by the low pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} threshold of the muon trigger, compared to the other LHC experiments. In addition, the LHCb experiment is probing a rapidity region only partially accessible by other LHC experiments. These properties allow the LHCb experiment to be complementary to similar analyses performed by the two central detectors at the LHC and even explore regions of the theoretical parameter space where these experiments are limited by their low efficiency to reconstruct highly boosted LLPs.
Detector description and simulation
The LHCb detector [14, 15] is a single-arm forward spectrometer covering the pseudorapidity range 2<η<5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2<\eta <5$$\end{document}, designed for the study of particles containing b\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$b $$\end{document} or c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c $$\end{document} quarks. The detector includes a high-precision tracking system consisting of the VELO which is a silicon-strip detector surrounding the pp interaction region [16], a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4Tm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$4{\mathrm {\,Tm}}$$\end{document}, and three stations of silicon-strip detectors and straw drift tubes [17, 18] placed downstream of the magnet. The tracking system provides a measurement of the momentum, p\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p$$\end{document}, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200\,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c$$\end{document}. The minimum distance of a track to a primary pp collision vertex (PV), the impact parameter, is measured with a resolution of (15+29/pT)μm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(15+29/p_{\mathrm {T}})\,\upmu \text {m} $$\end{document}, where pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} is in \,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c$$\end{document}. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors [19]. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic (ECAL) and a hadronic calorimeter (HCAL) [20]. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [21]. The online event selection is performed by a trigger [22], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction. During data taking an alignment and calibration of the detector is performed in near real-time and used in the software trigger [23]. The same alignment and calibration information is propagated to the offline reconstruction.
LLP production processes considered in this paper, where the χ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1}$$\end{document} represents the LLP: a di-LLP production via a scalar particle h0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$h^0$$\end{document}; b non-resonant, direct LLP production from quark interactions, where X is a stable particle, with mass identical to the LLP. The LLP decays into a muon and two quarks: χ~10→μ+qiqj(μ-qi¯qj¯)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1} \rightarrow \mu ^{+} q_{i} q_{j}( \mu ^{-} \bar{q_{i}} \bar{q_{j}})$$\end{document}
Simulation is used to model the effects of the detector acceptance and the imposed selection requirements. In the simulation, pp collisions are generated using Pythia 8 [24, 25] with a specific LHCb configuration [26] and with parton density functions taken from CTEQ6L [27]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [28, 29] as described in Ref. [30]. The simulation includes pileup events with an average of 1.1 pp visible interactions per bunch crossing.
Several sets of signal events have been produced assuming the processes illustrated in Fig. 1, where the χ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1}$$\end{document} plays the role of a long-lived particle. For the first process considered, two χ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1}$$\end{document} particles are obtained from the decay of the Higgs-like boson produced by gluon fusion, gg→h0→χ~10χ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ gg \rightarrow h^0 \rightarrow \tilde{\chi }^{0}_{1} \tilde{\chi }^{0}_{1} $$\end{document}. For the second process, the LLP is produced in a non-resonant mode, qq¯→χ~10X\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ q\bar{q} \rightarrow \tilde{\chi }^{0}_{1} X $$\end{document}. Here X is a stable neutral particle with the same mass as that of the χ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1}$$\end{document} state. This production of a LLP in association with a stable particle X is included, which enables probing the sensitivity to this topology, with the signal LLP recoiling against such a particle.
The LLP decays into a muon and two quarks; the branching ratio of χ~10→μ+qiqj(μ-qi¯qj¯)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tilde{\chi }^{0}_{1} \rightarrow \mu ^{+} q_{i} q_{j}( \mu ^{-} \bar{q_{i}} \bar{q_{j}})$$\end{document} is set to be equal for each quark combination (qi=u,c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$q_{i} = u,c$$\end{document} and qj=d¯,s¯,b¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$q_{j} = \bar{d},\bar{s},\bar{b}$$\end{document}), with an equal proportion of μ+\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu ^{+}$$\end{document} and μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu ^{-}$$\end{document}.
In the following, the model name is indicated by the values of mh0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0}$$\end{document}, mχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{\tilde{\chi }^{0}_{1}}$$\end{document} and τχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau _{\tilde{\chi }^{0}_{1}}$$\end{document}; h125-chi40-10ps, for example, corresponds to mh0=125\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0} =125\text {\,GeV\!/}c^2 $$\end{document}, mχ~10=40\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{\tilde{\chi }^{0}_{1}}=40\text {\,GeV\!/}c^2 $$\end{document}, τχ~10=10\,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau _{\tilde{\chi }^{0}_{1}}=10\text {\,ps} $$\end{document}. For the direct production, the Higgs mass is omitted from this notation, such as for example in chi30-10ps.
The most relevant background in this analysis is from events containing heavy quarks. The background from heavy quarks directly produced in pp collisions, as well as from W\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${W} $$\end{document}, Z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${Z} $$\end{document}, Higgs boson and top quark decays, is studied using the simulation. The simulation of inclusive bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} and cc¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${c} {\overline{{c}}} $$\end{document} events is not efficient to produce a large enough sample to cover the relevant high-pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} muon kinematic region. Hence, a dedicated sample of 20×106\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$20 \times 10^6$$\end{document} (1×106)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1 \times 10^6)$$\end{document} simulated bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} (cc¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${c} {\overline{{c}}} $$\end{document}) events has been produced with a minimum parton p^T\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hat{p}_\mathrm{T}$$\end{document} of 20\,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c$$\end{document} and requiring a muon with pT>12\,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}} >12\text {\,GeV\!/}c $$\end{document} and 1.5<η<5.0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.5<\eta <5.0$$\end{document}. All the simulated background species are suppressed by the multivariate analysis presented in the next section. Therefore, a data-driven approach is employed for the final background estimation.
Signal selection
Signal events are selected by requiring a vertex displaced from any PV in the event and containing one isolated, high-pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} muon. Due to the relatively high LLP mass, the muons from the LLP decay are expected to be more isolated than muons from hadron decays. The events from pp collisions are selected online by a trigger requiring muons with pT>10\,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}} >10\text {\,GeV\!/}c $$\end{document}. The offline analysis requires that the triggering muon has an impact parameter, IPμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathrm{IP^{\mu }}$$\end{document}, with respect to any PV, larger than 0.25\,mm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0.25\text {\,mm} $$\end{document} and a transverse momentum, pTμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${p^{\mu }_{\mathrm {T}}}$$\end{document}, larger than 12\,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$12\text {\,GeV\!/}c $$\end{document}. Primary and displaced vertices are reconstructed offline from charged particle tracks [31]. Genuine PVs are identified by a small radial distance from the beam axis, Rxy<0.3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy} <0.3$$\end{document} mm. Once the set of PVs is identified, all the other vertices are candidates for the decay position of LLPs. An LLP candidate is formed by requiring three or more tracks including the muon and having an invariant mass above 4.5\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}. There is no requirement for the reconstructed momentum to point to a specific PV. Particles interacting with the detector material are an important source of background. Therefore, a geometric veto is used to reject candidates with vertices in regions occupied by detector material [32]. The event preselection requires at least one PV in the event and at least one LLP candidate.
Figure 2 compares the distributions from data and from the simulated bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events for the relevant observables, after preselection. For illustration the shapes of simulated h125-chi40-10ps events are also superimposed. The effect of the geometric veto is visible in the Rxy\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy}$$\end{document} distribution, for candidates with Rxy\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy}$$\end{document} above 5\,mm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,mm}$$\end{document}. From simulation, the veto introduces a loss of efficiency of 3% (27%) for the detection of LLPs with a 50\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document} mass and a 10 ps (200 ps) lifetime, mh0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0}$$\end{document} =125\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}. The muon-isolation variable is defined as the sum of the energy of tracks surrounding the muon direction, including the muon itself, in a cone of radius Rηϕ=0.3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_{\eta \phi }= 0.3$$\end{document} in the pseudorapidity-azimuthal (η,ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(\eta , \phi )$$\end{document} space, divided by the energy of the muon track. The radius is reduced to Rηϕ=0.2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_{\eta \phi }= 0.2$$\end{document} when the theoretical hypothesis assumes a LLP mass of 10\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, to account for the reduced aperture of the jet of particles produced by the LLP decay. A muon-isolation value of unity denotes a fully isolated muon. In simulation the muon from the signal is found to be more isolated than the hadronic background. The variables σR\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma _R$$\end{document} and σZ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma _Z$$\end{document} are the vertex uncertainties in the radial direction and in the z direction respectively.
Distributions from data compared to simulated bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events (blue) and the simulated signal h125-chi40-10p (red), after preselection. From a to j: muon transverse momentum; muon impact parameter; muon isolation; the calorimetric energy, Ecalorimeters\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {E}_{\mathrm{calorimeters}}$$\end{document}, associated with the muon normalised by the muon energy, Emuon\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hbox {E}_{\mathrm{muon}}$$\end{document}; the number of tracks used to reconstruct the LLP vertex including the muon; the radial distance to the beam line of the reconstructed vertex; longitudinal and radial vertex fit errors, σZ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma _Z$$\end{document} and σR\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma _R$$\end{document}; reconstructed transverse momentum and mass of the LLP candidate. The distributions from simulated events are normalised to the data
The reconstructed vertex mass is very broad and does not peak at the neutralino mass values, because it misses some charged particle tracks, and any neutrals produced in the LLP decay.
The shapes of the distributions in Fig. 2 are all consistent with a dominant bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} composition of the background. This is confirmed by comparing the yields in data and simulation: after preselection and requiring the isolation parameter below 1.2, the total number of LLP candidates in data is 148×103\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$148\times 10^3$$\end{document}. The predicted background yields from bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} and cc¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${c} {\overline{{c}}} $$\end{document} events are (120±20)×103\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(120 \pm 20)\times 10^3$$\end{document} and (14±4)×103\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(14 \pm 4)\times 10^3$$\end{document}, respectively. Small contributions are expected from processes with W\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$W $$\end{document}, Z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$Z $$\end{document} bosons plus jets, top and Standard Model Higgs events: 260, 20, 2, and 1 candidates, respectively. The bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} and cc¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${c} {\overline{{c}}} $$\end{document} prediction uses the cross-sections measured by the LHCb experiment at 13 TeV [33, 34]. The acceptance of this analysis is computed with MadGraph5-aMC@NLO [35] and the detection efficiency is obtained from simulated events. As already stated, these background estimations are only used for cross-checks.
A multivariate analysis (MVA) based on a boosted decision tree [36, 37] is used to further purify the data sample. Ten MVA input variables are selected to optimise the signal–background separation. They are: pTμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${p^{\mu }_{\mathrm {T}}}$$\end{document} and IPμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathrm{IP^{\mu }}$$\end{document}, the ratio of the energies associated with the muon measured in ECAL and HCAL normalised to the muon energy, the LLP candidate pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document}, its pseudorapidity, the number of tracks forming the LLP, the vertex uncertainties σR\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma _R$$\end{document} and σZ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma _Z$$\end{document}, and the vertex Rxy\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy}$$\end{document} distance.
Larger vertex uncertainties are expected on the vertices of candidates from bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events compared to signal LLPs. The former are more boosted and produce more collimated tracks, while the relatively heavier signal LLPs decay into more divergent tracks. This effect decreases when the mass of the LLP approaches the mass of b\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$b $$\end{document}-quark hadrons. The selection based on the energy deposit in the calorimeters is efficient to suppress the background due to kaons or pions punching through the calorimeters and being misidentified as muons. The muon-isolation variable and the reconstructed mass of the long-lived particles are not included in the classifier; the discrimination power of these two variables is subsequently exploited for the signal determination.
The signal MVA training samples are provided by simulation. The background training sample is obtained from data, based on the hypothesis that the fraction of signal in the data after preselection is small. This automatically includes all possible background sources, with the correct relative abundance.
The training is performed independently for each simulated model. The MVA classifier is subsequently applied to the data and to the simulated signal. For each model, the optimal MVA cut value is chosen by an iterative minimization procedure to give the best expected cross-section upper limit, but keeping at least ten candidates to allow the invariant-mass fit to work properly.
The classifier can be biased by the presence of signal in data used as background training set. To quantify the potential bias, the MVA training is performed adding a fraction of simulated signal events (up to 5%) to the background set. This test demonstrates a negligible effect on the MVA performance for all the signal models.
Reconstructed invariant mass of the LLP candidates. Subfigures a, c, and e correspond to the signal selections which assume the models h80-chi30-10ps, h200-chi20-10ps, and the non-resonant model chi30-10ps, respectively. Subfigures b, d, and f are the corresponding distributions for candidates selected in the background region. The results of the fits are superimposed
Determination of the signal yield
The signal yield is determined with an unbinned extended maximum-likelihood fit to the distribution of the reconstructed LLP mass. The shape of the signal component is taken from the simulated models, and a background component is added. After the MVA selection, no simulated background survives, therefore the background shape is determined by a data-driven method, which also avoids potential simulation mismodeling of the reconstructed mass. The data candidates are separated into a signal region with muon isolation below 1.2 and a background region with isolation values from 1.4 to 2.0. The signal-region selection accepts more than 80% of the signal for all the models considered (see e.g. Fig. 2). Any potential signal yield in the background region is considered negligible. The reconstructed mass distribution obtained from the background candidates is used to constrain an empirical probability density function (PDF) consisting of the sum of two negative-slope exponential functions, one of them convolved with a Gaussian function. Shape parameters and amplitudes are left to vary in the fit. It is possible that the mass distribution obtained after selection of the background region does not represent exactly the background component in the signal region. Hence, a correction is applied before performing the fit: the mass distribution selected in the background region is weighted with weights deduced from the comparison of the candidate mass distributions of signal and background regions obtained from data with a relaxed MVA selection. This relaxed selection is required to have sufficiently populated samples and to minimise the correlation with the final distributions from which signal yields are obtained. The consistency of this procedure is tested on bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} simulated events.
Signal detection efficiency, in percent, after preselection, εpresel\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\varepsilon } _{presel}$$\end{document}, including the geometrical acceptance, and after MVA selection, ε\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\varepsilon } $$\end{document}, and numbers of fitted signal and background events in the signal region, Ns\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N_\mathrm{s}$$\end{document} and Nb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N_\mathrm{b}$$\end{document}, for the different signal hypotheses with resonant LLP production (masses units are \,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, lifetimes in \,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,ps}$$\end{document}). The last column gives the value of χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^2$$\end{document} per degree of freedom, ndf, from the fit
Signal detection efficiency, in percent, after preselection, εpresel\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\varepsilon } _{presel}$$\end{document}, including the geometrical acceptance, and after MVA selection, ε\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\varepsilon } $$\end{document}, and numbers of fitted signal and background events in the signal region, Ns\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N_\mathrm{s}$$\end{document} and Nb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N_\mathrm{b}$$\end{document}, for the different signal hypotheses for the LLP non-resonant production models (masses units are \,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, lifetimes in \,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,ps}$$\end{document}). The last column gives the value of χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^2$$\end{document} per degree of freedom, ndf, from the fit
Examples of the invariant mass of the selected LLP candidates are shown in Fig. 3 for the signal and background regions. The invariant-mass fit is performed simultaneously on LLP candidates from the signal and from the background regions. In the former, the numbers of signal and background events are free parameters of the fit. The results of the fit are shown in the figure. The sensitivity of the fit procedure is studied by adding a small number of simulated signal events to the data according to a given signal model. The fitted yields are on average consistent with the numbers of added events. The fitted signal yields, given in Tables 1 and 2 are compatible with the background-only hypothesis for all the theoretical models.
Detection efficiency and systematic uncertainties
The detection efficiency required in the calculation of the signal yield is estimated from the simulated signal events. The efficiencies after preselection and after MVA selection are shown in Tables 1 and 2, for the considered models of resonant and non-resonant LLP productions, respectively. The values include the geometrical acceptance. Several phenomena compete to determine the detection efficiency. In general the efficiency after preselection increases with the LLP mass because more particles are produced in the decay of heavier LLPs. There is a loss of particles outside the spectrometer acceptance, especially when the LLPs are produced from the decay of heavier states, such as the Higgs-like particle. In addition, the lower boost of heavier LLPs results in a shorter average flight length, which is disfavoured by the requirement of a minimum Rxy\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy}$$\end{document} value. With increasing LLP lifetimes a larger portion of the decays falls into the material region and is vetoed. Finally, a drop of sensitivity is expected for LLPs with a lifetime close to the b\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$b $$\end{document}-hadron lifetimes, where the contamination from bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events becomes even more important, especially for low-mass LLPs. The detection efficiency is reduced by up to one order of magnitude after the optimised MVA selection while the background is reduced by 3–4 orders of magnitude.
Contributions to the relative systematic uncertainties. The indicated ranges cover the theoretical models considered. The contributions from the signal and background models used in the LLP mass fit are treated separately
A breakdown of the relative systematic uncertainties is shown in Table 3. The uncertainties of the partonic luminosity depend upon the process considered; they are estimated following the procedure explained in Refs. [38, 39] and vary from 3% up to 6%, which is found for the gluon fusion process. The integrated luminosity [40] contributes with an uncertainty of 2%. The statistical precision of the efficiencies determined from simulation is in the range 2–4% for the different models. Different sources of systematic uncertainty arising from discrepancies between data and simulation have been considered. The size of those discrepancies for the relevant observables are inferred from a comparison of the distributions obtained from data and from bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} simulated events, which describes the data quite completely, or from other calibration processes.
Upper limits at 95% CL on the production cross-section times branching ratio for signal models with a resonant production via an Higgs-like boson. Masses are given in \,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, lifetimes in \,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,ps}$$\end{document}, cross-sections in \,pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,pb}$$\end{document}
The muon detection efficiency, including trigger, tracking, and muon identification efficiencies, is studied by a tag-and-probe technique applied to muons from J/ψ→μ+μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J/\psi \rightarrow {\mu ^+\mu ^-} $$\end{document}, Υ(1S)→μ+μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\varUpsilon } (1S) \rightarrow {\mu ^+\mu ^-} $$\end{document} and Z→μ+μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${Z} \rightarrow {\mu ^+\mu ^-} $$\end{document} decays. The corresponding systematic effects due to differences between data and simulation are estimated to be between 2 and 3.7%, depending on the theoretical model considered.
A comparison of the simulated and observed pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} distributions of muons from Z→μ+μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${Z} \rightarrow {\mu ^+\mu ^-} $$\end{document} decays shows a maximum difference of 0.2\,GeV\!/c\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c$$\end{document} in the selected region; this difference is propagated to the LLP analysis by shifting the muon pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} threshold by the same amount. The corresponding systematic uncertainty is below 1% for all models under consideration.
Upper limits at 95% CL on the production cross-section times branching ratio for signal models with a non-resonant production. Masses are given in \,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, lifetimes in \,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,ps}$$\end{document}, cross-sections in \,pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,pb}$$\end{document}
Expected (open dots and 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} and 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\sigma $$\end{document} bands) and observed (full dots) cross-section times branching fraction upper limits (95% CL) as a function of mχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hbox {m}_{\tilde{\chi }^{0}_{1}}$$\end{document} for the resonant production processes with mh0=125\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0} =125\text {\,GeV\!/}c^2 $$\end{document}, and, from a to g, τχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau _{\tilde{\chi }^{0}_{1}}$$\end{document} of 5, 10, 20, 30, 50, 100, and 200\,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,ps}$$\end{document}
Expected (open dots and 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} and 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} bands) and observed (full dots) cross-section times branching fraction upper limits (95% CL) as a function of τχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau _{\tilde{\chi }^{0}_{1}}$$\end{document} for the resonant production with mh0=125\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0} =125\text {\,GeV\!/}c^2 $$\end{document}, and, from a to e, mχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{\tilde{\chi }^{0}_{1}}$$\end{document} of 20, 30, 40, 50, and 60\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}
Expected (open dots and 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} and 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} bands) and observed (full dots) cross-section times branching fraction upper limits (95% CL) as a function of mh0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0}$$\end{document} and, from a to f, mχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{\tilde{\chi }^{0}_{1}}$$\end{document} of 10, 20, 30, 40, 50, and 60\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}
Expected (open dots and 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} and 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} bands) and observed (full dots) cross-section times branching fraction upper limits (95% CL), a: as a function of τχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau _{\tilde{\chi }^{0}_{1}}$$\end{document} with mχ~10=30\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{\tilde{\chi }^{0}_{1}}=30\text {\,GeV\!/}c^2 $$\end{document}, b: as a function of mχ~10\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{\tilde{\chi }^{0}_{1}}$$\end{document} with τχ~10=10\,ps\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau _{\tilde{\chi }^{0}_{1}}=10\text {\,ps} $$\end{document}. The processes are from direct, non-resonant, LLP production
The muon impact-parameter distribution is also studied from Z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${Z} $$\end{document} decays and shows a discrepancy between data and simulation of about 10μm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\,\upmu \text {m}$$\end{document} close to the pTμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${p^{\mu }_{\mathrm {T}}}$$\end{document} threshold. By changing the minimum IPμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathrm{IP^{\mu }}$$\end{document} requirement by this amount, the change in the detection efficiency is below 1% for all the models.
The vertex reconstruction efficiency has a complicated spatial structure due to the geometry of the VELO and the material veto. Uncertainties in the estimated vertex-finding efficiency are due to the per-track efficiency, track resolution, and differences in the contribution from background tracks due to the underlying interaction and pile-up. In the material-free region, Rxy<4.5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy} < 4.5$$\end{document} mm, the efficiency as a function of the flight distance has been studied in the context of lifetime measurements [41], showing that the simulation reproduces the data within 1%. In the region Rxy>4.5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy} > 4.5$$\end{document} mm a deviation of less than 6% is inferred from the study of inclusive bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events in data and simulation. By altering the efficiency in the simulation program as a function of the true vertex position, the effect on the LLP detection efficiency is estimated to be 1–2%. A second method to determine this contribution uses vertices from B0→J/ψK∗0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B^0 \rightarrow {{J /\psi }} {{K} ^{*0}} $$\end{document} decays with J/ψ→μ+μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{J /\psi }} \!\rightarrow {\mu ^+\mu ^-} $$\end{document} and K∗0→K+π-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{K} ^{*0}} \rightarrow K^+ \pi ^-$$\end{document}. For this process the vertex detection efficiencies in data and simulation agree within 10%. This result, obtained from a process with four final-state particles, is propagated to the LLP decay into a larger number of charged particle tracks and a detection threshold of three tracks. A discrepancy of at most 2% between the LLP efficiency in data and simulation is found, which is adopted as a contribution to the systematic detection uncertainty.
The uncertainty on the position of the beam line in the transverse plane is less than 20μm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\,\upmu \text {m}$$\end{document} [16]. It can affect the secondary-vertex selection, mainly via the requirement on Rxy\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R_\mathrm{xy}$$\end{document}. By altering the PV position in simulated signal events, the effect is estimated to be below 1%.
The effect of the imperfect modelling on the observables used in the MVA training is estimated with pseudoexperiments. As previously stated, the bias on each input variable is determined by comparing simulated and experimental distributions of muons and LLP candidates from Z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$Z $$\end{document} and W\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$W $$\end{document} events, as well as from bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events. At the MVA test stage, each input variable is modified by a scale factor randomly selected from a Gaussian distribution of width equal to the corresponding bias. The standard deviation of the signal efficiency distribution is taken as a systematic uncertainty.
The signal and background samples are obtained through a selection on the muon isolation parameter. By a comparison of data and muons from simulated bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events, the maximum uncertainty on this variable is estimated to be ±0.015\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\pm 0.015$$\end{document} in the proximity of the thresholds, with a maximal effect on the efficiency of 1.7%.
Comparing the mass distributions of bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} and Z→bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$Z\rightarrow {{b} {\overline{{b}}}} $$\end{document} events, a maximum mass-scale discrepancy between data and simulated events of 10% is estimated in the proximity of the threshold, which translates into a 1.4% contribution to the detection efficiency uncertainty.
Finally, the total systematic uncertainty is obtained as the sum in quadrature of all contributions, where the different components of the detection efficiency are assumed to be fully correlated.
The choice of the signal and background invariant-mass templates can affect the results of the LLP mass fits. The uncertainty due to the signal model accounts for the mass scale and the mass resolution. The mass scale and resolution discrepancies between data and simulation are below 1% and 1.5% respectively, as obtained from bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} and Z→bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$Z \rightarrow {{b} {\overline{{b}}}} $$\end{document} events. Pseudoexperiments are used to estimate the effect on the cross-section calculation. For each theoretical model, ten simulated signal events are added to the selected data after a Gaussian smearing or after changing the mass scale. The average deviation of the observed upper limits with respect to the one obtained from the default signal and background distributions is below 2%.
The background shape is deduced from data selected in the poorly isolated region after reweighting, with weights inferred from the data distributions obtained with relaxed selection criteria. The overall uncertainty is estimated by reducing by half the weights and running pseudoexperiments as before. The average deviation of the observed upper limits is below 14%.
Results
The 95% confidence level (CL) upper limits, expected and observed, on the production cross-sections times branching fraction are computed for each model using the CLs approach [42]. Statistical and systematic uncertainties on the signal efficiencies are included as nuisance parameters of the likelihood function, assuming Gaussian distributions. Finally, the upper limit values are corrected by the factors which account for the imperfect modelling of signal and background templates.
The numerical results for all the models are given in Tables 4 and 5. Figures 4, 5, 6 and 7 show the measured cross-section times branching ratio upper limits, for different theoretical models. The decrease of sensitivity for relatively low LLP mass value is explained by the above-mentioned effects on the detection efficiency. The upper limits for the processes with mh0=125\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$m_{h^0} =125\text {\,GeV\!/}c^2 $$\end{document} can be compared to the prediction of the Standard Model Higgs production cross-section from gluon fusion of about 46\,pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,pb}$$\end{document} at s=13\,TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s} =13\text {\,TeV} $$\end{document} [43].
Conclusion
Long-lived massive particles decaying into a muon and two quarks have been searched for using proton-proton collision data collected by the LHCb experiment at s=13\,TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s} =13\text {\,TeV} $$\end{document}, corresponding to an integrated luminosity of 5.4 \,fb-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,fb} ^{-1}$$\end{document}. The LLP lifetime range considered is from 5 to 200 ps. The background is dominated by bb¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b} {\overline{{b}}} $$\end{document} events and is reduced by tight selection requirements, including a dedicated multivariate classifier. The signal yield is determined by a fit to the LLP reconstructed mass with a signal shape inferred from the theoretical models.
The forward acceptance of the LHCb experiment makes it complementary to other LHC experiments, while its low trigger pT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$p_{\mathrm {T}}$$\end{document} threshold allows exploring relatively small LLP masses. Two types of LLP productions have been assumed. In the first a Higgs-like particle is produced by gluon fusion and decays into two LLPs. The analysis covers Higgs-like boson masses from 30 to 200\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}, and LLP mass range from 10 \,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document} up to about one half of the mass of the parent boson. The second mode is a direct LLP production from quark interactions, covering the LLP mass range from 10 up to 90 \,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}.
The results for all theoretical models considered are compatible with the background-only hypothesis. The upper limits at 95% CL set on the cross-section times branching fractions are mostly of O(0.1\,pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,pb}$$\end{document}), but the sensitivity is limited to O(10\,pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,pb}$$\end{document}) for the lowest LLP mass value considered of 10\,GeV\!/c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {\,GeV\!/}c^2$$\end{document}.
Acknowledgements
We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); MOST and NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); NWO (The Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MSHE (Russia); MICINN (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (UK); DOE NP and NSF (USA). We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (The Netherlands), PIC (Spain), GridPP (UK), RRCKI and Yandex LLC (Russia), CSCS (Switzerland), IFIN-HH (Romania), CBPF (Brazil), PL-GRID (Poland) and NERSC (USA). We are indebted to the communities behind the multiple open-source software packages on which we depend. Individual groups or members have received support from ARC and ARDC (Australia); AvH Foundation (Germany); EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union); A*MIDEX, ANR, IPhU and Labex P2IO, and Région Auvergne-Rhône-Alpes (France); Key Research Program of Frontier Sciences of CAS, CAS PIFI, CAS CCEPP, Fundamental Research Funds for the Central Universities, and Sci. & Tech. Program of Guangzhou (China); RFBR, RSF and Yandex LLC (Russia); GVA, XuntaGal and GENCAT (Spain); the Leverhulme Trust, the Royal Society and UKRI (UK).
Data Availability
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: All LHCb scientific output is published in journals, with preliminary results made available in Conference Reports. All are Open Access, without restriction on use beyond the standard conditions agreed by CERN. Data associated to the plots in this publication are made available on the CERN document server at http://cdsweb.cern.ch/record/2706539.]
ReferencesDimopoulosSRabySWilczekFSupersymmetry and the scale of unification1981D2416811981PhRvD..24.1681D10.1103/PhysRevD.24.1681S.P. Martin, A supersymmetry primer. Adv. Ser. Direct. High Energy Phys. 21, 1 (2010). https://doi.org/10.1142/9789812839657_0001. arXiv:hep-ph/9709356P.W. Graham, D.E. Kaplan, S. Rajendran, P. Saraswat, Displaced supersymmetry. JHEP 07, 149 (2012). https://doi.org/10.1007/JHEP07(2012)149. arXiv:1204.6038M.J. Strassler, K.M. Zurek, Discovering the Higgs through highly-displaced vertices. Phys. Lett. B 661, 263 (2008). https://doi.org/10.1016/j.physletb.2008.02.008. arXiv:hep-ph/0605193ATLAS Collaboration, G. Aad et al., Search for displaced vertices arising from decays of new heavy particles in 7 TeV pp collisions at ATLAS. Phys. Lett. B 707, 478 (2012). https://doi.org/10.1016/j.physletb.2011.12.057. arXiv:1109.2242ATLAS Collaboration, G. Aad et al., Search for massive, long-lived particles using multitrack displaced vertices or displaced lepton pairs in pp collisions at s=8\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s} = 8$$\end{document} TeV with the ATLAS detector. Phys. Rev. D 92, 072004 (2015). https://doi.org/10.1103/PhysRevD.92.072004. arXiv:1504.05162CMS Collaboration, A.M. Sirunyan et al., Search for new long-lived particles at s=13TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s}=13~{{\rm TeV}}$$\end{document}. Phys. Lett. B 780, 432 (2018). https://doi.org/10.1016/j.physletb.2018.03.019. arXiv:1711.09120ATLAS Collaboration, G. Aad et al., Search for long-lived, massive particles in events with a displaced vertex and a muon with large impact parameter in pp collisions at s=13TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s}=13~{{\rm TeV}}$$\end{document} with the atlas detector. Phys. Rev. D 102, 032006 (2020). https://doi.org/10.1103/physrevd.102.032006. arXiv:2003.11956LHCb Collaboration, R. Aaij et al., Search for massive long-lived particles decaying semileptonically in the LHCb detector. Eur. Phys. J. C 77, 224 (2017). https://doi.org/10.1140/epjc/s10052-017-4744-6. arXiv:1612.00945P.W. Graham, S. Rajendran, P. Saraswat, Supersymmetric crevices: missing signatures of R-parity violation at the LHC. Phys. Rev. D 90, 075005 (2014). https://doi.org/10.1103/PhysRevD.90.075005. arXiv:1403.7197L. Lee, C. Ohm, A. Soffer, T.-T. Yu, Collider searches for long-lived particles beyond the Standard Model. Prog. Part. Nucl. Phys. 106, 210 (2019). https://doi.org/10.1016/j.ppnp.2019.02.006. arXiv:1810.12602J. Alimena et al., Searching for long-lived particles beyond the standard model at the large hadron collider. J. Phys. G Nucl. Part. Phys. 47, 090501 (2020). https://doi.org/10.1088/1361-6471/ab4574. arXiv:1903.04497B.C. Allanach, A. Dedes, H.K. Dreiner, R parity violating minimal supergravity model. Phys. Rev. D 69, 115002 (2004). https://doi.org/10.1103/PhysRevD.69.115002, https://doi.org/10.1103/PhysRevD.72.079902. arXiv:hep-ph/0309196LHCb Collaboration, A.A. Alves Jr. et al., The LHCb detector at the LHC. JINST3, S08005 (2008). https://doi.org/10.1088/1748-0221/3/08/S08005LHCb Collaboration, R. Aaij et al., LHCb detector performance. Int. J. Mod. Phys. A 30, 1530022 (2015). https://doi.org/10.1142/S0217751X15300227. arXiv:1412.6352R. Aaij et al., Performance of the LHCb vertex locator. JINST 9, P09007 (2014). https://doi.org/10.1088/1748-0221/9/09/P09007. arXiv:1405.7808R. Arink et al., Performance of the LHCb outer tracker. JINST 9, P01002 (2014). https://doi.org/10.1088/1748-0221/9/01/P01002. arXiv:1311.3893P. d’Argent et al., Improved performance of the LHCb outer tracker in LHC Run 2. JINST 12, P11016 (2017). https://doi.org/10.1088/1748-0221/12/11/P11016. arXiv:1708.00819M. Adinolfi et al., Performance of the LHCb RICH detector at the LHC. Eur. Phys. J. C 73, 2431 (2013). https://doi.org/10.1140/epjc/s10052-013-2431-9. arXiv:1211.6759C. Abellan Beteta et al., Calibration and performance of the LHCb calorimeters in Run 1 and 2 at the LHC. Submitted to JINST. arXiv:2008.11556A.A. Alves Jr. et al., Performance of the LHCb muon system. JINST 8, P02022 (2013). https://doi.org/10.1088/1748-0221/8/02/P02022. arXiv:1211.1346R. Aaij et al., The LHCb trigger and its performance in 2011. JINST 8, P04022 (2013). https://doi.org/10.1088/1748-0221/8/04/P04022. arXiv:1211.3055R. Aaij et al., Performance of the LHCb trigger and full real-time reconstruction in Run 2 of the LHC. JINST 14, P04013 (2019). https://doi.org/10.1088/1748-0221/14/04/P04013. arXiv:1812.10790T. Sjöstrand, S. Mrenna, P. Skands, A brief introduction to PYTHIA 8.1. Comput. Phys. Commun. 178, 852 (2008). https://doi.org/10.1016/j.cpc.2008.01.036. arXiv:0710.3820T. Sjöstrand, S. Mrenna, P. Skands, PYTHIA 6.4 physics and manual. JHEP 05, 026 (2006). https://doi.org/10.1088/1126-6708/2006/05/026. arXiv:hep-ph/0603175BelyaevIHandling of the generation of primary events in Gauss, the LHCb simulation framework201133103204710.1088/1742-6596/331/3/032047J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD analysis. JHEP 07, 012 (2002). https://doi.org/10.1088/1126-6708/2002/07/012. arXiv:hep-ph/0201195Geant4 Collaboration, J. Allison et al., Geant4 developments and applications. IEEE Trans. Nucl. Sci. 53, 270 (2006). https://doi.org/10.1109/TNS.2006.869826Geant4 Collaboration, S. Agostinelli et al., Geant4: a simulation toolkit. Nucl. Instrum. Methods A 506, 250 (2003). https://doi.org/10.1016/S0168-9002(03)01368-8M. Clemencic et al., The LHCb simulation application, Gauss: design, evolution and experience. J. Phys. Conf. Ser. 331, 032023 (2011). https://doi.org/10.1088/1742-6596/331/3/032023M. Kucharczyk, P. Morawski, M. Witek, Primary vertex reconstruction at LHCb. CERN-LHCb-PUB-2014-044, LHCb-PUB-2014-044 (2014)LHCb Collaboration, R. Aaij et al., Search for the rare decay, KS0→μ+μ-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$K_{{\rm S}}^{0}\rightarrow \mu ^{+}\,\mu ^{-}$$\end{document}. JHEP01, 090 (2013). https://doi.org/10.1007/JHEP01(2013)090. arXiv:1209.4029LHCb Collaboration, R. Aaij et al., Measurement of the b-quark production cross-section in 7 and 13 TeV pp collisions. Phys. Rev. Lett. 118, 052002 (2017). https://doi.org/10.1103/PhysRevLett.118.052002. arXiv:1612.05140 [Erratum ibid. 119, 169901 (2017). https://doi.org/10.1103/PhysRevLett.119.169901]A. Pearce, LHCb early Run 2 measurements of B and charm production, LHCb-PROC-2016-001J. Alwall et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations. J. High Energ. Phys. 2014, 79 (2014). https://doi.org/10.1007/jhep07(2014)079. arXiv:1405.0301BreimanLFriedmanJHOlshenRAStoneCJ1984BelmontWadsworth International Group0541.62042FreundYSchapireREA decision-theoretic generalization of on-line learning and an application to boosting199755119147305510.1006/jcss.1997.15040880.68103M. Botje et al., The PDF4LHC working group interim recommendations. arXiv:1101.0538J. Butterworth et al., PDF4LHC recommendations for LHC Run II. J. Phys. G43, 023001 (2016). https://doi.org/10.1088/0954-3899/43/2/023001. arXiv:1510.03865LHCb Collaboration, R. Aaij et al., Precision luminosity measurements at LHCb. JINST 9, P12005 (2014). https://doi.org/10.1088/1748-0221/9/12/P12005. arXiv:1410.0149LHCb Collaboration, R. Aaij et al., Measurements of the B+\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B^{+}$$\end{document}, B0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B^{0}$$\end{document}, Bs0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B^{0}_{s}$$\end{document} meson and Λb0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\varLambda ^{0}_{b}$$\end{document} baryon lifetimes. JHEP 04, 114 (2014). https://doi.org/10.1007/JHEP04(2014)114. arXiv:1402.2554ReadALPresentation of search results: the CLs\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text{CL}_s$$\end{document} technique2002G2826932002JPhG...28.2693R10.1088/0954-3899/28/10/313ATLAS Collaboration, G. Aad et al., Combined measurements of Higgs boson production and decay using up to 80fb-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$80~{\rm fb}^{-1}$$\end{document} of proton–proton collision data at s=13TeV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sqrt{s}=13~{\rm TeV}$$\end{document} collected with the ATLAS experiment. Phys. Rev. D 101, 012002 (2020). https://doi.org/10.1103/PhysRevD.101.012002. arXiv:1909.02845