Updated search for long-lived particles decaying to jet pairs

Eur. Phys. J. C (2017) 77:812 https://doi.org/10.1140/epjc/s10052-017-5178-x Regular Article - Experimental Physics Updated search for long-lived particles decaying to jet pairs LHCb Collaboration CERN, 1211 Geneva 23, Switzerland Received: 23 May 2017 / Accepted: 31 August 2017 / Published online: 29 November 2017 © CERN for the benefit of the LHCb collaboration 2017. This article is an open access publication Abstract A search is presented for long-lived particles with a mass between 25 and 50 GeV/c2and a lifetime between 2 and 500 ps, using proton–proton collision data corresponding to an integrated luminosity of 2.0 fb−1, collected by the LHCb detector at centre-of-mass energies of 7 and 8 TeV. The particles are assumed to be pair-produced in the decay of a 125 GeV/c2Standard-Model-like Higgs boson. The experimental signature is a single long-lived particle, identified by a displaced vertex with two associated jets. No excess above background is observed and limits are set on the production cross-section as a function of the mass and lifetime of the long-lived particle. 1 Introduction Various extensions of the Standard Model (SM) feature new particles whose couplings to lighter states are sufficiently small to result in detectable lifetimes. In this paper we report on a search for such long-lived particles, which are assumed to be pair-produced in the decay of a Standard-Model-like Higgs boson, and subsequently decay into a quark–antiquark pair. Such a signature is present in models with a hiddensector non-Abelian gauge group, where the Standard Model Higgs boson acts as a portal [1–5]. The new scalar particle represents the lightest state in the hidden sector and is called a hidden-valley pion (πv) throughout this paper. Experimental constraints on the properties of the Higgs boson of mass 125 GeV/c2observed by the ATLAS and CMS collaborations [6,7] still allow for branching fractions of non-SM decay modes of up to 30% [8]. Data collected with the LHCb experiment in 2011 and 2012 are used for this analysis, restricted to periods in which suitable triggers were available. The data sample analysed corresponds to 0.62 fb−1at a centre-of-mass energy of √s=7 TeV and 1.38 fb−1at √s=8 TeV. In simulated events with πvpairs originating from a Higgs boson decay it e-mail: pieter.da[email protected] is found that in most cases no more than one of the two πv decays occurs inside the LHCb acceptance. Consequently, the experimental signature is a single πvparticle. The candidate is identified by its decay to two hadronic jets originating from a displaced vertex, with a transverse distance to the proton-proton collision axis (Rxy)ofatleast0.4mm.The vertex is required to have at least five tracks reconstructed in the LHCb vertex detector. The analysis is sensitive to πv particles with a mass between 25 and 50 GeV/c2and a lifetime between 2 and 500 ps. The lifetime range is limited due to the presence of large prompt backgrounds at short decay times and the acceptance of the vertex detector for long decay times. The lower boundary on the mass range arises from the requirement to identify two hadronic jets while the upper boundary is driven by the geometric acceptance of the detector. This paper presents an update of an earlier analysis, which considered only the data set corresponding to an integrated luminosity of 0.62 fb−1collected at √s=7TeV[9]. Similar searches for hidden-valley particles decaying to jet pairs were performed by the D0 [10], CDF [11], ATLAS [12–14] and CMS [15] collaborations. Compared to these analyses, this search is sensitive to πvparticles with relatively low mass and lifetime. The LHCb collaboration has also performed a search for events with two displaced high-multiplicity vertices [16] and a search for events with a lepton from a highmultiplicity displaced vertex [17] in the context of SUSY models, and several searches for so far unknown long-lived particles in B-meson decays [18–21]. 2 Detector and event simulation The LHCb detector [22,23] is a single-arm forward spectrometer covering the pseudorapidity range 2 <η<5, designed for the study of particles containing bor cquarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector (VELO) surrounding the pp interaction region, a large-area silicon-strip detec- 123 812 Page 2 of 14 Eur. Phys. J. C (2017) 77 :812 tor located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The tracking system provides a measurement of the momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of (15+(29 GeV/c)/pT)μm, where pTis the component of the momentum transverse to the collision axis. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad (SPD) and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers. The model for the production of πvparticles through the Higgs portal is fully specified by three parameters: the mass of the Higgs boson and the mass and lifetime of the πv.The Higgs boson mass is taken to be 125 GeV/c2, and its production through the gluon-gluon fusion process is simulated with the Pythia8 generator [24], with a specific LHCb configuration [25] and using the CTEQ6 leading-order set of parton density functions [26]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [27,28] as described in Ref. [29]. Signal samples with πvmasses of 25, 35, 43 and 50 GeV/c2 and lifetimes of 10 and 100 ps are generated. In the simulated events the long-lived particles decay exclusively as πv→b¯ b, since this decay mode is generally preferred in the Higgs portal model. Samples with decays to c- and s-quark pairs are generated as well, but only in the scenario with a mass of 35 GeV/c2and a lifetime of 10 ps. 3 Event selection The experimental signature for this analysis is a single displaced vertex with two associated jets. Only decays that produce a sufficient number of tracks in the VELO for a vertex to be reconstructed are considered. Due to the geometry of the vertex detector, this restricts the sample to decay points up to about 200 mm from the nominal interaction point along the beam direction, and up to about 30 mm in the transverse direction, thereby limiting the decay time acceptance. The selection strategy is the same as used in the analysis of Ref. [9]. Reconstructed tracks are used to find the decay vertex, and jets are built out of reconstructed particles compatible with originating from that vertex. Constraints on the signal yield are determined from a fit to the dijet invariant mass distribution. The main source of background is displaced vertices from heavy-flavour decays or interactions of particles with detector material. To take into account the strong dependence of the background level on the separation from the beam axis, different selection criteria are used in different bins of Rxy, and the final fit is performed in bins of this variable. The selection consists of online (trigger) and offline parts. The trigger [30] is divided into a hardware (L0) and a software (HLT) stage. The L0 requires a muon with high pTor a hadron, photon or electron with high transverse energy in the calorimeters. In order to reduce the processing time of the subsequent trigger stages, events with a large hit multiplicity in the SPD are discarded. The software stage is divided into two parts, which for this analysis differ between the 2011 and 2012 data. In the 2011 sample, the first software stage (HLT1) requires a single high-pTtrack with a large impact parameter. The HLT1 selection for the 2012 sample was complemented with a two-track vertex signature with looser track quality criteria, in order to improve the efficiency at large displacements. At the second stage of the software trigger (HLT2), events are required to pass either a dedicated inclusive displaced-vertex selection or a standard topological B decay selection, which requires a two-, three- or four-track vertex with a significant displacement from all PVs [30]. The inclusive displaced-vertex selection uses an algorithm similar to that used for the LHCb primary vertex reconstruction [31]. A combination of requirements on the minimum number of tracks in the vertex (at least four), the distance Rxy of the vertex to the beam axis (at least 0.4 mm), the invariant mass of the particles associated with the vertex (at least 2 GeV/c2) and the scalar sum pTof the tracks that form the vertex (at least 3 GeV/c), is used to define a set of trigger selections with sufficiently low rate. Before the offline selection can be applied, the displaced vertex corresponding to the decay of the πvcandidate must be reconstructed. For those events in which the HLT2 inclusive displaced-vertex selection was successful, the same vertex candidate found in the trigger is used; this approach differs from that used in the previous LHCb analysis [9] and simplifies the evaluation of systematic uncertainties. For events selected only by the topological Btrigger, a modified version of the algorithm is run on the output of the offline reconstruction with the following criteria: vertices with 0.4<Rxy <1 mm must have at least eight tracks and the invariant mass of the system must exceed 10 GeV/c2, vertices with 1 <Rxy <5 mm must have at least six tracks, and those with Rxy >5 mm must have at least five tracks. To exclude background due to interactions with the detector material, vertices inside a veto region around the VELO detector elements are discarded. Events with many parallel displaced tracks, which can arise from machine background, are identified by the azimuthal distribution of hits in the VELO and are also discarded. Next, jets are reconstructed following a particle flow approach. The same set of inputs as in Ref. [32]isused, 123 Eur. Phys. J. C (2017) 77 :812 Page 3 of 14 812 namely tracks of charged particles and calorimeter energy deposits, after subtraction of the energy associated with charged particles. To remove background, tracks that are compatible with coming from a PV, tracks with a smaller impact parameter to any primary vertex than to the displaced vertex, and tracks that have an impact parameter to the displaced vertex larger than 2 mm are all discarded. The anti-kTjet clustering algorithm is used [33], with a distance parameter of R=0.7. The jet momentum and jet mass are calculated from the four-vectors of all constituents of the jet. In simulated events the jet energy response is found to be close to unity except for the lowest jet momenta, near the minimally required transverse momentum of 5 GeV/c. Therefore, no jet energy correction was applied for this search. To enhance the jet purity the fraction of the jet energy carried by charged particles should be at least 0.1, there should be at least one track with transverse momentum above 0.9 GeV/c, no pair of constituents should carry 90% of the jet energy, and no single charged or neutral constituent should contribute more than 70 or 50% of the total energy, respectively. To ensure that they can reliably be associated to a vertex, the jets are also required to have at least two constituents with track segments in the VELO. To account for differences in trigger and background conditions, for the 2012 data this requirement was tightened to at least four segments for Rxy <1 mm, and at least three segments for 1<Rxy <2 mm. For each jet an origin point is reconstructed from the jet constituents with VELO information. The jet trajectory is defined based on this origin point and the momentum of the jet. Any jet whose trajectory does not point back to the candidate vertex within 2 mm, or points more closely to a primary vertex, is removed. Only candidates with at least two jets passing these criteria are retained. Two final criteria are applied to the dijet candidates. The first is that the momentum vector of the dijet candidate should be aligned with the displacement vector from a PV to the reconstructed vertex position. This is implemented as a requirement on the dijet invariant mass divided by the corrected mass, m/mcorr >0.7. The corrected mass is computed as mcorr =m2+(psin θ)2+psin θ[34], where mand pare the reconstructed mass and momentum of the dijet, and θis the minimum angle between the momentum vector and the displacement vectors to the vertex from any PV in the event. A requirement on m/mcorr is preferred over one on the angle θitself, since its efficiency depends less strongly on the boost and the mass of the candidate [35]. The second criterion is that the kinematic separation of the jets should satisfy R=(η)2+(φ)2<2.2, where η and φ are the pseudorapidity and azimuthal angle differences between the two jets, respectively. This reduces the tail in the dijet invariant mass distribution by suppressing the remaining back-to-back dijet background. Table 1 Number of selected candidates per generated H0→πvπv event (efficiency) in percent for different πv→qq,q=b,c,smodels for 2011 and 2012 data taking conditions, as derived from simulation. The relative statistical uncertainty on the efficiency due to the limited size of the simulated sample is less than a few percent πvmass (GeV/c2) 2011 2012 10 ps 100 ps 10 ps 100 ps πv→b¯ b25 0.45 0.097 0.46 0.111 πv→b¯ b35 0.80 0.176 0.83 0.224 πv→b¯ b43 0.73 0.190 0.77 0.222 πv→b¯ b50 0.49 0.141 0.54 0.171 πv→c¯c35 1.35 1.35 πv→s¯s35 1.30 1.19 The overall efficiency to reconstruct and select displaced πvdecays in the simulated samples is summarized in Table 1for the 2011 and 2012 data taking conditions. A large part of the inefficiency is due to the detector acceptance, which is about 13% (8%) and 6.5% (5.5%) for πvparticles with a lifetime of 10 ps (100 ps) and masses of 25 and 50 GeV/c2, respectively. Other important contributions are due to the selection on the displacement from the beamline, requirements on the minimum number of tracks forming the vertex, the material interaction veto, the reduction in VELO tracking efficiency at large displacements, and the jet selection [36]. The efficiency for long-lived particles decaying to s- and cquark pairs is higher than for decays to b-quark pairs due to the larger number of tracks originating directly from the πv decay vertex. 4 Systematic uncertainties Systematic uncertainties on the efficiency are obtained from studies of data-simulation differences in control samples. They are reported in Tables 2and 3, for the 2011 and 2012 conditions, respectively, and discussed in more detail below. Uncertainties on the signal efficiency due to parton-density distributions, the simulation of fragmentation and hadronization, and the Higgs boson production cross-section and kinematics are not taken into account. The vertex reconstruction efficiency can be split into two parts, namely the track reconstruction efficiency and the vertex finding efficiency. The track reconstruction efficiency is described by the simulation to within a few percent, including for highly displaced and low-momentum tracks [37–39]. The effect of a systematic change in this efficiency is studied by randomly removing 2% of the signal tracks and reapplying all selection criteria. The vertex finding algorithm is not fully efficient even if all tracks are reconstructed. In particular, the efficiency to 123 812 Page 4 of 14 Eur. Phys. J. C (2017) 77 :812 Table 2 Overview of the contributions to the relative systematic uncertainty on the signal efficiency and luminosity (in percent) for different signal samples in 2011 conditions. The uncertainty on the total efficiency is obtained by summing the individual contributions in quadrature πvmass ( GeV/c2)2535435035,c¯c35, s¯s πvlifetime (ps) 10 100 10 100 10 100 10 100 10 10 Tracking efficiency 4.24.13.33.23.02.83.02.71.81.7 Vertex finding 3.84.23.33.92.83.73.72.62.92.8 Jet reconstruction 3.13.11.61.60.70.70.50.50.91.0 Jet identification 3.03.03.03.03.03.03.03.03.03.0 Jet direction 7.07.06.06.07.47.48.58.55 .95.7 L0 4.04.03.03.03.03.02.02.01.82.1 NSPD 1.71.72.02.01.61.62.32.31.71.6 HLT1 2.02.02.02.02.02.02.02.02.02.0 HLT2 3.03.03.03.03.03.03.03.03.03.0 Total efficiency 11.511.69.810.010.310.511.210.98.78.6 Luminosity 1.71.71.71.71.71.71.71.71.71.7 Table 3 Overview of the contributions to the relative systematic uncertainty on the signal efficiency and luminosity (in percent) for different signal samples in 2012 conditions. The uncertainty on the total efficiency is obtained by summing the individual contributions in quadrature πvmass ( GeV/c2)2535435035,c¯c35, s¯s πvlifetime (ps) 10 100 10 100 10 100 10 100 10 10 Tracking efficiency 3.12.82.42.42.22.12.01.71.21.1 Vertex finding 4.24.53.84.43.44.13.13.93.43.5 Jet reconstruction 2.72.71.11.10.70.70.30.30.91.0 Jet identification 3.03.03.03.03.03.03.03.03.03.0 Jet direction 5.85.85.35.36.16.17.97.95 .35.8 L0 4.04.02.52.52.02.02.02.02.02.0 NSPD 2.22.22.52.52.52.52.52.52.42.1 HLT1 2.02.02.02.02.02.02.02.02.02.0 HLT2 3.03.03.03.03.03.03.03.03.03.0 Total efficiency 10.510.69.29.49.19.510.410.68.68.9 Luminosity 1.21.21.21.21.21.21.21.21.21.2 find a low-multiplicity secondary vertex is reduced in the proximity of a high-multiplicity PV. The effect is studied in data and simulation using exclusively reconstructed B0→ J/Ψ K∗0decays, which can be selected with high purity without tight requirements on the vertex. The efficiency for the displaced vertex reconstruction algorithm to find the B0 candidate is measured as a function of the displacement Rxy in data and simulation [36]. The difference, weighted by the Rxy distribution of the signal candidates, is used to derive a systematic uncertainty. Systematic uncertainties related to the jet reconstruction can be introduced in two ways: through differences between data and simulation in the jet reconstruction efficiency and through differences between data and simulation in the resolution on the jet energy and direction, which enter the dijet candidate kinematic and m/mcorr selection and the dijet invariant mass shape. The jet reconstruction efficiency has been studied previously in measurements of the Z+jet and Z+b-jet cross-sections and was found to be consistent between data and simulation [32,40]. The Z→μ+μ−+jet sample is used to study jet-related systematic effects for this analysis as well. To mimic the selection of the particle-flow inputs, the PV associated to the Zis used as a proxy for the displaced vertex. The difference between data and simulation with the largest impact on the jet reconstruction efficiency is the energy response to low-pTjets, close to the threshold of 5GeV/c. The sensitivity to a different energy response in data and simulation is evaluated by increasing the minimum jet pTfor candidates passing the full offline selection by 10%, which is the uncertainty on the jet energy scale. The change in the overall selection efficiency is assigned as a systematic uncertainty. By replacing the jet identification criteria with a requirement on the pTbalance between the leading jet and the Zboson, the Z→μ+μ−sample can also be used to study the difference in jet identification efficiency between data and simulation. No difference larger than 3% relative is seen, which is assigned as a systematic uncertainty. To validate the simulation of the jet-direction resolution the jet-direction is estimated separately with the charged and 123 Eur. Phys. J. C (2017) 77 :812 Page 5 of 14 812 Fig. 1 Dijet invariant mass distribution in the different Rxy bins, for the 2011 data sample. For illustration, the best fit with a signal πvmodel with mass 35 GeV/c2and lifetime 10 ps is overlaid. The solid blue line indicates the total background model, the short-dashed green line indicates the signal model for signal strength μ=1, and the long-dashed red line indicates the best-fit signal strength R R R R R R neutral components of the jet in Z+jet events. The distribution of the charged-neutral difference in the estimated direction is found to be consistent between data and simulation for both the ηand the φprojection, and across the full range of pT. To quantify the effect on the πvsignal efficiency, an additional smearing to the jet-direction is applied to jets of selected candidates in the simulation. The jet angles with respect to the beam direction are smeared independently in the horizontal and vertical planes by about one third of the resolution, which is the largest value compatible with the comparison of data and simulation in Z+jet events. The systematic uncertainty related to the L0 trigger selection consists of two parts, due to differences in the L0 calorimeter trigger response between data and simulation, and due to the difference between data and simulation in the distribution of the SPD hit multiplicity NSPD. The first is evaluated by studying the L0 calorimeter trigger response on jets reconstructed in Z+jet events, where the trigger decision is made based on the Z→μ+μ−decay products, and is independent of the jet. The observed data-simulation differences are propagated to the πvreconstruction efficiency and correspond to systematic uncertainties of 2–4%, depending on the πvmass. Jets in Z+jet events are mostly light-quark jets, while our benchmark signal decays to bquarks. It is found in simulated events that the efficiency of the L0 calorimeter trigger is practically independent of jet flavour. A small fraction of b-quark jets is triggered exclusively by the L0 muon trigger, which is well modelled in the simulation. The second part of the L0 systematic uncertainty arises because the SPD multiplicity is not well described in the simulation. This effect is studied with a Z→μ+μ−sample triggered by the dimuon L0 selection, which applies only a loose selection on this quantity. An efficiency correction is derived, which is about 90% for 2011 data, and about 85% for 2012 data, with an uncertainty of 2–3%. The difference in the correction between the different πvmodels is smaller 123 812 Page 6 of 14 Eur. Phys. J. C (2017) 77 :812 than the systematic variation. This correction is applied to the overall detection efficiency derived from the simulation and the uncertainty is taken as a systematic uncertainty. The differences between data and simulation in the HLT1 selection are dominated by the track reconstruction efficiency, which was discussed above, and additional track quality criteria. One such difference is due to a requirement on the number of VELO hits for displaced tracks. It is characterized using B0→J/Ψ K∗0decays selected with triggers that do not apply such a requirement. For this sample the selection efficiency was found to be 2% higher in data than in simulated events, which is assigned as a systematic uncertainty. For πvdecays the final-state track multiplicity is larger, which dilutes effects due to a mismodelling of the single-track efficiency. The main source of systematic uncertainty in the HLT2 selection is the vertex reconstruction efficiency, which was discussed above. The efficiency of the topological Btrigger, which is relevant for a subset of the candidates, is accurately described in simulation. It is measured as a function of Rxy in data and simulation using B0→J/Ψ K∗0candidates that are selected by a different, dimuon-based, trigger criterion. A maximum difference of 2–3% is observed, which is assigned as a systematic uncertainty. 5 Results Constraints on the presence of a signal are derived from a fit to the dijet invariant mass distributions, shown in Figs. 1 and 2. To take advantage of the difference in the Rxy distribution for background and signal, the data are divided into six Rxy bins. The data are further split according to data taking year to account for differences in running conditions and Higgs boson production cross-section. The signal efficiency for each Rxy bin is obtained from the simulated samples with πvlifetimes of 10 and 100 ps, with the decay time dis- Fig. 2 Dijet invariant mass distribution in the different Rxy bins, for the 2012 data sample. For illustration, the best fit with a signal πvmodel with mass 35 GeV/c2and lifetime 10 ps is overlaid. The solid blue line indicates the total background model, the short-dashed green line indicates the signal model for signal strength μ=1, and the long-dashed red line indicates the best-fit signal strength R R R R R R 123 Eur. Phys. J. C (2017) 77 :812 Page 7 of 14 812 Fig. 3 Expected (open circles and dotted line) and observed (filled circles and solid line) upper limit versus lifetime for different πvmasses and decay modes. The green (dark) and yellow (light) bands indicate the quantiles of the expected upper limit corresponding to ±1σand ±2σfor a Gaussian distribution. The decay πv→b¯ bis assumed, unless specified otherwise tributions reweighted to mimic other lifetime hypotheses as needed. Results are presented as upper limits on the signal strength μ≡(σ/σ SM gg→H0)·B(H0→πvπv), where σis the excluded signal cross-section, σSM gg→H0is the SM Higgs boson production cross-section via the gluon fusion process and B(H0→πvπv)is the branching fraction of the Higgs boson decay to πvparticles. The branching fraction Bqq of the πvparticle to the q¯qfinal state (with q¯q=b¯ b,c¯cor s¯s depending on the final state under study) is assumed to be 100%. If the decay width of the πvparticle is dominated by other decays than that under study, the limits scale as 1/(Bqq (2−Bqq )). The Higgs boson production cross-section is assumed to be 15.11 pb at 7 TeV and 19.24 pb at 8 TeV [41]. The CLsmethod [42] is used to determine upper limits. The profile likelihood ratio qμ PLL =L(μ, ˆ θ(μ))/L(ˆμ, ˆ θ) is chosen as a test statistic, where L(μ, θ) denotes the likelihood as a function of μand a set of nuisance parameters θ, which are also extracted from the data; L(μ, ˆ θ(μ)) is the maximum likelihood for a hypothesized value of μand L(ˆμ, ˆ θ) is the global maximum likelihood. To estimate the sensitivity of the analysis and the significance of a potential signal, the expected upper limit quantiles in the case of zero signal are also evaluated. For each value of μand θthe likelihood is evaluated as L(μ, θ) =iP(xi;μ, θ), where Pis the probability density for event iand the product runs over all selected events. The observables xifor each candidate include the dijet mass, Rxy bin and data taking year. For each Rxy bin and data taking year, the invariant mass distribution is modelled by the sum of background and signal components. The distribution for the signal is modelled as a Gaussian distribution whose parameters are obtained from fully simulated signal events. For the background distribution an empirical model, outlined below, is adopted. Background candidates can be categorized into two contributions. The first category is mostly due to the combination 123 812 Page 8 of 14 Eur. Phys. J. C (2017) 77 :812 Table 4 Observed 95% CL signal strength (μ) upper limits for different πvmodels πvmass πvlifetime (ps) 2 5 10 20 50 100 200 500 25 GeV/c21.64 0.83 1.12 1.22 2.84 4.37 9.28 22.82 35 GeV/c20.63 0.35 0.32 0.41 0.76 1.37 2.56 5.86 43 GeV/c20.52 0.21 0.16 0.21 0.35 0.63 1.12 2.77 50 GeV/c20.50 0.17 0.14 0.15 0.25 0.41 0.76 1.72 35 GeV/c2,πv→c¯c0.33 0.17 0.16 0.20 0.39 0.64 1.19 2.90 35 GeV/c2,πv→s¯s0.40 0.20 0.19 0.24 0.42 0.77 1.41 3.51 of a heavy-flavour decay vertex or an interaction with detector material with particles from a primary interaction. This contribution has a steeply decreasing invariant mass spectrum. Following the approach in Ref. [9], the distribution is modelled by the convolution of a falling exponential distribution with a bifurcated Gaussian. All parameters of this background model are free to vary in the fit. The second category is due to Standard Model dijet events. These events have candidates with jets that are approximately back-to-back in the transverse plane. It is suppressed by the selection on the dijet opening angle R. Its remaining contribution has a less steeply falling mass spectrum. It is described in the fit with a similar functional shape as for the first category, but with the parameters and the relative yields in the different bins fixed from a fit to the invariant mass distribution of candidates that fail the Rrequirement. In the final fit only the total normalization of this component is varied. The second component is new compared to the model used for the previous analysis [9]. It leads to a better description of the high-mass tail, at the expense of one extra fit parameter for each data taking year. It was found that the result of the fit is not sensitive to the exact Rrequirement used to select the events for this component. All parameters of the fit to the invariant mass distribution are allowed to float independently in each bin, except for the following nuisance parameters: the dijet invariant mass scale, the overall signal efficiency, and the normalization for the second background contribution. All relevant systematic uncertainties are incorporated in the fit model: the overall uncertainty on the efficiency, as described in Sect. 4,the uncertainty on the dijet invariant mass scale, and the uncertainties on the shape parameters and relative normalisation arising from the finite size of the simulated signal samples. Gaussian constraints on these parameters are added to the likelihood. Alternatives have been considered for the background mass model, in particular with an additional less steeply falling exponential to describe the tail. With these models the estimated background yield at higher mass is similar or Fig. 4 Observed upper limit versus lifetime for different πvmasses and decay modes. The decay πv→b¯ bis assumed, unless specified otherwise larger than with the nominal background model, leading to tighter limits on the signal. As the nominal model gives the most conservative limit, no additional systematic uncertainty is assigned for background modeling. There is no significant excess of signal in the data. Upper limits at 95% confidence level (CL) as a function of lifetime for hidden-valley models with different πvmass and decay mode are shown in Fig. 3and summarized in Table 4and Fig. 4. The best sensitivity is obtained for a mass of about 50 GeV/c2and a lifetime of about 10 ps. The main improvements with respect to the previous result [9] are due to the enlarged data sample, the improved trigger selections, and the addition of the Rxy bin above 5 mm, which contributes to the increased sensitivity at larger lifetimes. 6 Conclusion Results have been presented from a search for long-lived particles with a mass in the range 25–50 GeV/c2and a lifetime between 2 and 500 ps. The particles are assumed to be pair- 123 Eur. Phys. J. C (2017) 77 :812 Page 9 of 14 812 produced in the decay of a 125 GeV/c2Standard-Model-like Higgs boson and to decay into two jets. Besides decays to b¯ b, which are the best motivated in the context of hidden-valley models [1,2], also decays to c¯cand s¯squark pairs are considered. No evidence for so far unknown long-lived particles is observed and limits are set as a function of mass and lifetime. These measurements complement other constraints on this production model at the LHC [13,15] by placing stronger constraints at small masses and lifetimes. 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); MinES and FASO (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); 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 (United Kingdom), RRCKI and Yandex LLC (Russia), CSCS (Switzerland), IFINHH (Romania), CBPF (Brazil), PL-GRID (Poland) and OSC (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 AvH Foundation (Germany), EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union), Conseil Général de Haute-Savoie, Labex ENIGMASS and OCEVU, Région Auvergne (France), RFBR and Yandex LLC (Russia), GVA, Xunta- Gal and GENCAT (Spain), Herchel Smith Fund, The Royal Society, Royal Commission for the Exhibition of 1851 and the Leverhulme Trust (United Kingdom). Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. References 1. M.J. Strassler, K.M. Zurek, Echoes of a hidden valley at hadron colliders. Phys. Lett. B 651, 374 (2007). arXiv:hep-ph/0604261 2. M.J. Strassler, K.M. Zurek, Discovering the Higgs through highly-displaced vertices. Phys. Lett. B 661, 263 (2008). arXiv:hep-ph/0605193 3. S. Chang, R. Dermisek, J.F. Gunion, N. Weiner, Nonstandard Higgs boson decays. Annu. Rev. Nucl. Part. Sci. 58, 75 (2008). arXiv:0801.4554 4. N. Craig, A. Katz, M. Strassler, R. Sundrum, Naturalness in the dark at the LHC. JHEP 07, 105 (2015). arXiv:1501.05310 5. D. Curtin, C.B. Verhaaren, Discovering uncolored naturalness in exotic Higgs decays. JHEP 12, 072 (2015). arXiv:1506.06141 6. ATLAS collaboration, G. Aad et al., Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. Phys. Lett. B 716, 1 (2012). arXiv:1207.7214 7. CMS collaboration, S. Chatrchyan et al., Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC. Phys. Lett. B 716, 30 (2012). arXiv:1207.7235 8. ATLAS and CMS collaborations, G. Aad et al., Measurements of the Higgs boson production and decay rates and constraints on its couplings from a combined ATLAS and CMS analysis of the LHC pp collision data at √s=7and8TeV.JHEP08, 045 (2016). arXiv:1606.02266 9. LHCb collaboration, R. Aaij et al., Search for long-lived particles decaying to jet pairs. Eur. Phys. J. C 75, 152 (2015). arXiv:1412.3021 10. D0 collaboration, V.M. Abazov et al., Search for resonant pair production of long-lived particles decaying to b¯ bin p¯pcollisions at √s=1.96 TeV. Phys. Rev. Lett. 103, 071801 (2009). arXiv:0906.1787 11. CDF collaboration, T. Aaltonen et al., Search for heavy metastable particles decaying to jet pairs in pp collisions at √s=1.96 TeV. Phys.Rev.D85, 012007 (2012). arXiv:1109.3136 12. ATLAS collaboration, G. Aad et al., Search for a light Higgs boson decaying to long-lived weakly-interacting particles in proton– proton collisions at √s=7 TeV with the ATLAS detector. Phys. Rev. Lett. 108, 251801 (2012). arXiv:1203.1303 13. ATLAS collaboration, G. Aad et al., Search for long-lived, weakly interacting particles that decay to displaced hadronic jets in proton– proton collisions at √s=8 TeV with the ATLAS detector. Phys. Rev. D 92, 012010 (2015). arXiv:1504.03634 14. ATLAS collaboration, G. Aad et al., Search for pair-produced longlived neutral particles decaying in the ATLAS hadronic calorimeter in pp collisions at √s=8 TeV. Phys. Lett. B 743, 15 (2015). arXiv:1501.04020 15. CMS collaboration, V. Khachatryan et al., Search for long-lived neutral particles decaying to quark–antiquark pairs in proton– proton collisions at √s=8 TeV. Phys. Rev. D 91, 012007 (2015). arXiv:1411.6530 16. LHCb collaboration, R. Aaij et al., Search for Higgs-like bosons decaying into long-lived exotic particles. Eur. Phys. J. C 76, 664 (2016). arXiv:1609.03124 17. LHCb collaboration, R. Aaij et al., Search for massive long-lived particles decaying semileptonically in the LHCb detector. Eur. Phys.J.C77, 224 (2017). arXiv:1612.00945 18. LHCb collaboration, R. Aaij et al., Search for Majorana neutrinos in B−→π+μ−μ−decays. Phys. Rev. Lett. 112, 131802 (2014). arXiv:1401.5361 19. LHCb collaboration, R. Aaij et al., Searches for Majorana neutrinos in B−decays. Phys. Rev. D 85, 112004 (2012). arXiv:1201.5600 20. LHCb collaboration, R. Aaij et al., Search for hidden-sector bosons in B0→K∗0χ(μ+μ−)decays. Phys. Rev. Lett. 115, 161802 (2015). arXiv:1508.04094 21. LHCb collaboration, R. Aaij et al., Search for long-lived scalar particles in B+→K+χ(μ+μ−)decay. Phys. Rev. D 95, 071101 (2017). arXiv:1612.07818 22. LHCb collaboration, A.A. Alves Jr. et al., The LHCb detector at the LHC. JINST 3, S08005 (2008) 23. LHCb collaboration, R. Aaij et al., LHCb detector performance. Int. J. Mod. Phys. A 30, 1530022 (2015). arXiv:1412.6352 24. T. Sjöstrand, S. Mrenna, P. Skands, A brief introduction to PYTHIA 8.1. Comput. Phys. Commun. 178, 852 (2008). arXiv:0710.3820 25. I. Belyaev et al., Handling of the generation of primary events in Gauss, the LHCb simulation framework. J. Phys. Conf. Ser. 331, 032047 (2011) 26. J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD analysis. JHEP 07, 012 (2002). arXiv:hep-ph/0201195 27. Geant4 collaboration, J. Allison et al., Geant4 developments and applications. IEEE Trans. Nucl. Sci. 53, 270 (2006) 123