Abstract
A search is presented for massive long-lived particles, in the 20–60 \({\mathrm {Ge V\!/}c^2}\) mass range with lifetimes between 5 and 100 \(\mathrm{ps}\). The dataset used corresponds to 0.62\( \text{ fb }^{-1}\) of proton-proton collision data collected by the LHCb detector at \(\sqrt{s} =7\mathrm {\,Te V} \). The particles are assumed to be pair-produced by the decay of a Higgs-like boson with mass between 80 and 140 \({\mathrm {Ge V\!/}c^2}\). No excess above the background expectation is observed and limits are set on the production cross-section as a function of the long-lived particle mass and lifetime and of the Higgs-like boson mass.
Similar content being viewed by others
1 Introduction
The standard model of particle physics (SM) has shown great success in describing physics processes at very short distances. Nevertheless, open questions remain, such as the hierarchy problem, the imprecise unification of gauge couplings, and the absence of candidates for dark matter. Considerable efforts have been made to address these issues, resulting in a large variety of models. Supersymmetry (SUSY), in which the strong and electroweak forces are unified at a renormalisation scale near the Planck scale, provides a possible solution for the hierarchy problem; the minimal supersymmetric standard model (MSSM) is the simplest, phenomenologically viable realisation of SUSY [1, 2].
The present study focuses on a subset of models featuring massive long-lived particles (LLP) with a measurable flight distance. We concentrate on scenarios in which the LLP decays hadronically in the LHCb vertex detector, travelling distances which can be larger than those of typical \(b \) hadrons.
A large number of LLP searches have been performed by the experiments at the LHC and Tevatron, mainly using the Hidden Valley framework [3] as a benchmark model [4,5,6,7,8]. Hidden Valley processes have also been sought by LHCb [9], which is able to explore the forward rapidity region only partially covered by other LHC experiments. In addition, it is able to trigger on particles with low transverse momenta, allowing the experiment to probe relatively small LLP masses.
The event topology considered in this study is quite different from that of Hidden Valley models. The minimal supergravity model (mSUGRA) realisation of the MSSM is used as a benchmark model with baryon number violation [10], as suggested in Refs. [11, 12]. Here a Higgs-like boson produced in pp collisions decays into two LLPs (neutralinos), subsequently decaying into three quarks each. The Higgs-like particle mass ranges from 80 up to 140 \({\mathrm {Ge V\!/}c^2}\), covering the mass of the scalar boson discovered by the ATLAS and CMS experiments [13, 14]. The explored LLP lifetime range of 5–100 \(\mathrm{ps}\) is higher than the typical \(b \) hadron lifetime, and corresponds to an average flight distance of up to 30 cm, which is inside the LHCb vertex detector region. The LLP mass range considered is between 20 and 60 \({\mathrm {Ge V\!/}c^2}\).
2 Detector description
The LHCb detector [15, 16] is a single-arm forward spectrometer covering the pseudorapidity range \(2<\eta <5\), designed for the study of particles containing \(b \) or \(c \) quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region (VELO), a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about \(4{\mathrm {\,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 \({\mathrm {Ge V\!/}c}\). The minimum distance of a track to a primary vertex (PV), the impact parameter, is measured with a resolution of (15 \(+\) 29/\(p_\mathrm{T}\)) \(\upmu \mathrm{m}\), where \(p_\mathrm{T}\)is the component of the momentum transverse to the beam, in \({\mathrm {Ge V\!/}c}\). 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 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 online event selection is performed by a trigger [17], which consists of a hardware stage, L0, based on information from the calorimeter and muon systems, followed by two software stages, HLT1 and HLT2, which run a simplified version of the offline event reconstruction.
3 Event generation and detector simulation
Various simulated event samples are used in this analysis. The pp collisions are generated with Pythia 6 [18]. The process simulated is \({\mathrm {h}^0} \rightarrow {\tilde{\chi }^{0}_{1}}{\tilde{\chi }^{0}_{1}}\), where the Higgs-like boson of mass \(m_\mathrm{h^0}\) is produced via gluon-gluon fusion, with the parton density function taken from CTEQ6L [19]. The neutralino \(\tilde{\chi }^{0}_{1}\) is an LLP of mass \(m_\mathrm{{LLP}}\) and lifetime \(\tau _\mathrm{{LLP}}\), which decays into three quarks via the mSUGRA baryon number violating process available in Pythia. The corresponding decay flavour structure for the neutralino with a mass of \(48\,{\mathrm {Ge V\!/}c^2} \) is 18.5% for each of the combinations with a \(b \) quark (udb, usb, cdb, csb), and 13% for each udq and cdq, where q is not a \(b \) quark, i.e. about 75% of LLPs have a \(b \) quark in the decay. This fraction becomes 70% for \(m_\mathrm{{LLP}} =20\,{\mathrm {Ge V\!/}c^2} \).
Two separate detector simulations are used, a full simulation where the interaction of the generated particles with the detector is based on Geant4 [20, 21], and a fast simulation. In Geant4, the detector and its response are implemented as described in Ref. [22]. Signal models for a representative set of theoretical parameters have been generated and fully simulated (Appendix A, Table 5). In the remainder of this paper, the following nomenclature is chosen: a prefix “BV”, indicating baryon number violation, is followed by the LLP mass in \({\mathrm {Ge V\!/}c^2}\) and lifetime, and the prefix “mH” followed by the \(m_\mathrm{h^0}\) value in \({\mathrm {Ge V\!/}c^2}\). Most of the fully simulated models have \(m_\mathrm{h^0} \)=114 \({\mathrm {Ge V\!/}c^2}\), which is in the middle of the chosen Higgs-like particle mass range. Only events with at least one \(\tilde{\chi }^{0}_{1}\) in the pseudorapidity region \(1.8<\eta <5.0\) are processed by Geant4, corresponding to about 30% of the generated events.
The fast simulation is used to cover a broader parameter space of the theoretical models. Here the charged particles from the \({\mathrm {h}^0} \rightarrow {\tilde{\chi }^{0}_{1}} {\tilde{\chi }^{0}_{1}}\) process falling in the geometrical acceptance of the detector are processed by the vertex reconstruction algorithm. The fast simulation is validated by comparison with the full simulation. The detection efficiencies predicted by the full and the fast simulation differ by less than 5% for all the signal models. The distributions for mass, momentum and transverse momentum of the reconstructed LLP, and for the reconstructed vertex position coincide.
Events with direct production of charm, bottom and top quarks are considered as sources of background. Samples of such events were produced and fully simulated. In particular, \(17 \times 10^6\) inclusive \({b} {\overline{{b}}} \) events (\(9 \times 10^6\) inclusive \({c} {\overline{{c}}} \) events) were produced with at least two \(b \) hadrons (\(c \) hadrons) in \(1.5<\eta <5.0\), and half a million \({t} {\overline{{t}}} \) events with at least one muon in the acceptance.
4 Event selection and signal determination
This analysis searches for events with pairs of displaced high-multiplicity vertices. The main background is due to secondary interactions of particles with the detector material. These events are discarded by a material veto, which rejects vertices in regions occupied by detector material [23]. The remaining candidates are found to be compatible with \({b} {\overline{{b}}} \) events.
From simulation, LLP candidates within the detector acceptance are selected by the L0 and HLT1 triggers with an efficiency of more than 85%. The simulation indicates that the trigger activity is dominated by the hadronic component of the signal expected from high multiplicity events. In HLT2, primary vertices and displaced vertices are reconstructed from charged tracks [24]. Genuine PVs are identified by a small radial distance from the beam axis, \(R_\mathrm{xy} <0.3\) mm, and must have at least 10 tracks, including at least one forward track (i.e. in the direction of the spectrometer) and one backward track. Once the set of PVs is identified, all other reconstructed vertices are candidates for the decay position of LLPs. The preselection requires at least one PV in the event and two LLP candidates. The LLP candidates must have at least four forward tracks, no backward tracks, and a minimum invariant mass reconstructed from charged tracks larger than 3.5 \({\mathrm {Ge V\!/}c^2}\) for one candidate, and larger than 4.5 \({\mathrm {Ge V\!/}c^2}\) for the other. In addition, the two secondary vertices must have \(R_\mathrm{xy} > 0.4\mathrm \,mm \) and pass the material veto.
The preselection criteria drastically suppress the hadronic background. Only 37 events (74 LLP candidates) survive from the simulated set of \(17.1 \times 10^6\) \({b} {\overline{{b}}} \) events generated in the LHCb acceptance, corresponding to an integrated luminosity of 0.3\( \text{ pb }^{-1}\). Three simulated \({c} {\overline{{c}}} \) events pass the selection. They contain \(b \) hadrons and hence belong to the category of inclusive \({b} {\overline{{b}}} \), which is also the case of the two surviving \({t} {\overline{{t}}} \) events. From the 0.62\( \text{ fb }^{-1}\) data sample, \(42.9 \times 10^3\) events are selected. The \({b} {\overline{{b}}} \) cross-section value measured by LHCb, \( 288 \pm 4 \pm 48\) \(\upmu \)b [25, 26], predicts \((76 \pm 22)\times 10^3\) events, \(1.8\pm 0.5\) times the yield observed in data. The estimate uses the next-to-leading-order POWHEG calculation [27] to correct Pythia, and the detection efficiency obtained from the simulated events. The measured yield has also been compared to the rate observed in LHCb by a dedicated inclusive \({b} {\overline{{b}}} \) analysis, based on a topological trigger [28]. The consistency with the \({b} {\overline{{b}}} \) background is verified within a statistical precision of 10%.
The shapes of the distributions of the relevant observables are compatible with the \({b} {\overline{{b}}} \) background. Figure 1 compares the distributions for the LLP candidates taken from data and from simulated \({b} {\overline{{b}}} \) events. The distributions for three fully simulated signal models are also shown. The mass and the \(p_\mathrm{T}\)values are calculated assuming the pion mass for each charged track. Figure 1d presents the radial distribution of the displaced vertices; the drop in the number of candidates with a vertex above \(R_\mathrm{xy} \sim 5\mathrm \,mm \) is due to the material veto. The variables \(\sigma _\mathrm{R}\) and \(\sigma _\mathrm{Z}\) shown in Fig. 1e, f are the position uncertainties provided by the vertex fit in the transverse distance \(R_\mathrm{xy}\) and along the z axis, parallel to the beam. The values of \(\sigma _\mathrm{R}\) and \(\sigma _\mathrm{Z}\) are larger for the candidates from \({b} {\overline{{b}}} \) background than for the signal because light boosted particles produce close parallel tracks, with the consequence that the vertex fit has larger uncertainties than for the decay of heavier particles producing more diverging tracks. Figure 2 presents the LLP distance of flight and \(R_\mathrm{xy}\) distributions compared to three fully simulated signal models, corresponding to \(\tau _\mathrm{{LLP}}\) values of 5, 10, and 50 \(\mathrm{ps}\).
The reconstructed four-vectors of the two LLPs in the event are added to form the Higgs-like candidate (di-LLP), the corresponding invariant mass and \(p_\mathrm{T}\)distributions are given in Fig. 3.
Further cuts are applied to the preselected data, to increase the statistical sensitivity. The figure of merit used is given by \(\epsilon / \sqrt{N_d + 1}\), where \(\epsilon \) is the signal efficiency from simulation for a given selection, and \(N_d\) the corresponding number of candidates found in the data. The baseline selection (\(\mathrm Sel_1\)) is defined by a minimum number of charged tracks on each vertex \(N^\mathrm{track}_\mathrm{min} = 6\), a minimum reconstructed mass \(m^\mathrm{LLP}_\mathrm{min} =6\,{\mathrm {Ge V\!/}c^2} \), and maximum uncertainties from the vertex fit \(\sigma ^\mathrm{R}_\mathrm{max} =0.05 \mathrm \,mm \), and \(\sigma ^\mathrm{Z}_\mathrm{max} =0.25\mathrm \,mm \). All the selections used in this analysis are described in Table 1, with the indication of the number of data events selected for a di-LLP reconstructed mass above 19 \({\mathrm {Ge V\!/}c^2}\). Selection \(\mathrm Bkg_1\) is used to model the background in the fit procedure described in Sect. 5, selections \(\mathrm Sel_2\) and \(\mathrm Bkg_2\) are used to study systematic effects.
5 Determination of the di-LLP signal
The signal yield is determined by a fit of the di-LLP invariant mass, assuming that the two LLPs are the decay products of a narrow resonance. This technique is hampered by the difficulty in producing a reliable background model from simulation, despite the fact that it is reasonable to believe that only \({b} {\overline{{b}}} \) events are the surviving SM component. Therefore, in this analysis the alternative is chosen to infer the background model from data by relaxing the selection requirements, as given by lines \(\mathrm Bkg_1\) and \(\mathrm Bkg_2\) of Table 1. The comparison of the results obtained with the different signal and background selections is subsequently used to estimate the systematic effects.
The signal template is the histogram built from BV simulated events selected under the same conditions as data, i.e. \(\mathrm Sel_1\). The background template is the histogram obtained from data events selected by the \(\mathrm Bkg_1\) conditions. The number of signal (background) candidates \(N_s\) (\(N_b\)) is determined by an extended maximum likelihood fit. The results are given in Fig. 4 for the BV48 10 ps mH114 signal. The fit \(\chi ^2/\mathrm {ndf}\) is 0.6. Note that only the portion of the di-LLP mass spectrum above 19 \({\mathrm {Ge V\!/}c^2}\) is used, in order to be sufficiently above the mass threshold set by the selections. Alternatively, \(\mathrm Sel_2\) and \(\mathrm Bkg_2\) are used to assess systematic effects. The fit results for the selections (\(\mathrm Sel_1\),\(\mathrm Bkg_2\)), (\(\mathrm Sel_2\),\(\mathrm Bkg_1\)) are shown in Fig. 5. The corresponding fit \(\chi ^2/\mathrm {ndf}\) values are 0.6 and 1.0. The results are given in Table 2 for all fully simulated signal models. All fits give a negative number of signal candidates, compatible with zero. These results are correlated because the data sample is in common and the di-LLP mass shapes are almost identical for the different fully simulated models as depicted in Fig. 3. A check is performed on 142 di-LLP candidates selected from simulated \({b} {\overline{{b}}} \) background without the requirement on \(R_\mathrm{xy}\) and with \(m^\mathrm{LLP}_\mathrm{min} =4\,{\mathrm {Ge V\!/}c^2} \) for both LLPs. The fitted number of signal events is \(-0.8\pm 3.5\).
The behaviour and sensitivity of the procedure is further studied by adding a small number of signal events to the data according to a given signal model. Figure 6 shows the results for two models with 10 signal events added to the data. The fitted \(N_s\) corresponds well to the number of injected signal events.
An alternative fit procedure has been applied, using parameterised signal and background templates. The sum of two exponential functions is used for the background, and an exponential convolved with a Gaussian function for the signal. The results are consistent with a null signal for all the models.
As a final check a two-dimensional sideband subtraction method (“ABCD method” [29]) has been applied in the reconstructed mass of one LLP and the number of tracks of the other LLP, also giving results consistent with zero signal.
6 Detection efficiency and systematic uncertainties
The determination of the detection efficiency is based on simulated events. The geometrical acceptance for the detection of one \(\tilde{\chi }^{0}_{1}\) in LHCb is, depending on the model, between 20 and 30%. After selection \(\mathrm Sel_1\) the predicted total di-LLP detection efficiency is in the range 0.1–1% for most of the models. Potential discrepancies between simulation and data are considered as sources of systematic uncertainties. Table 3 summarises the contributions of the systematic uncertainties, which are valid for all fully simulated models, dominated by the 15% contribution from the trigger.
The consistency between the trigger efficiency in data and simulation is checked by selecting LLP events with an independent trigger, designed for the detection of \({J /\psi }\) events. Comparing the fraction of the data that also passes the double-LLP selection with the corresponding fraction in simulated inclusive \({J /\psi }\) events, consistent efficiencies are found within a statistical uncertainty of 30%. A more precise result is obtained when requiring only a single LLP candidate [9] and assuming uncorrelated contributions from the two LLPs to determine the efficiency for detecting two LLPs in coincidence. A maximum discrepancy between data and simulation of 15% is inferred, which is the value adopted.
The consistency between the track reconstruction efficiency in data and simulation is studied by a comparison of the number of tracks selected in displaced vertices from \({b} {\overline{{b}}} \) events. The average number of tracks per LLP in data is higher than in simulated events by about 0.07 tracks. Assuming that this small effect is entirely due to a difference in tracking efficiency, the overall di-LLP detection efficiency changes by at most 5%.
The vertex reconstruction efficiency is affected by the tracking efficiency and resolution. A study of vertices from \(B^0 \rightarrow {{J /\psi }} {{K} ^{*0}} \) with \({{J /\psi }} \!\rightarrow {\mu ^{+}\mu ^{-}} \) and \({{K} ^{*0}} \rightarrow K^{+} \pi ^{-}\) has shown that the data and simulation detection efficiencies for this four-prong process agree within 7.5% [9]. This has been evaluated to correspond at most to a 4% discrepancy between the di-LLP efficiency in data and simulation.
A maximum mismatch of 10% on both the transverse momentum and mass scales is inferred from the comparison of data and simulated \({b} {\overline{{b}}} \) distributions, which propagates to a 6% contribution to the systematic uncertainty.
The effect of the material veto corresponds to a reduction of the geometrical acceptance and depends mainly on the LLP lifetime. An analysis with the requirement of \(R_\mathrm{xy} < 4\mathrm \,mm \) allows to infer a maximum systematic uncertainty of 4%.
A small contribution to the systematic uncertainty of 0.1% is determined by reweighting the simulated events to match the PV multiplicity in the data.
The uncertainty on the position of the beam line is less than 20\({\,\upmu \mathrm m}\) [30]. It can affect the secondary vertex selection, mainly via the requirement on \(R_\mathrm{xy}\). By altering the PV position in simulated signal events, the maximum effect on the di-LLP selection efficiency is 0.7%.
The Higgs-like particle production model is mainly affected by the uncertainty on the parton luminosity. A maximum variation of the detection efficiency of 9.5% is obtained following the prescriptions given in [31]. A second contribution of 3% is obtained by reweighting the Pythia generated events to match a recent calculation of the \(p_\mathrm{T}\)distributions [32]. The total theoretical uncertainty is 9.9%, obtained by summing in quadrature the mentioned contributions.
In addition to the systematic uncertainty on the detection efficiency, the following contributions have been considered. The uncertainty on the integrated luminosity is 1.7% [33]. As previously stated, the uncertainty on the momentum scale and the invariant mass scale is smaller than 10%. This value is also assumed for the di-LLP mass calibration. To assess the impact on the signal measurement, pseudoexperiments are produced with 10 events of simulated signal added to the background following the nominal signal distribution but with the di-LLP mass value scaled by \(\pm 10\%\). The subsequent maximum variation of the fitted number of events is \(\pm 1.6\), for all the signal hypotheses. The uncertainty due to the shape of the background template is obtained by comparing the number of fitted events obtained with the \(\mathrm B\,kg_1\) and \(\mathrm B\,kg_2\) selections. The change is less than one event, for all the signal models. The difference in data and simulation in the di-LLP mass resolution and the statistical precision of the signal templates used in the fit have a negligible effect. Hence, a fit uncertainty of ±2 events is considered in the calculation of the cross-section upper limits.
For the analysis based on the fast simulation, a 5% uncertainty is added to account for the relative imprecision of the fast simulation with respect to the full simulation, as explained in Sect. 3.
7 Results
The 95% confidence level (CL) upper limits on the production cross-section times branching ratio are presented in Table 4, for the fully simulated models, based on the CLs approach [34]. The fast simulation allows the exploration of a larger region of parameter space. The cross-section times branching fraction upper limits at 95% CL for benchmark theoretical models are shown in Fig. 7 (the corresponding tables are given in Appendix C).
The estimated detection efficiencies can be found in Appendix B, Tables 6 and 7. The efficiency increases with \(m_\mathrm{{LLP}}\) because more particles are produced in the decay of heavier LLPs. This effect is only partially counteracted by the loss of particles outside of the spectrometer acceptance, which is especially the case with heavier Higgs-like particles. Another competing phenomenon is that the lower boost of heavier LLPs results in a shorter average flight length, i.e. the requirement of a minimum \(R_\mathrm{xy}\) disfavours heavy LLPs. The cut on \(R_\mathrm{xy}\) is more efficient at selecting LLPs with large lifetimes, but for lifetimes larger than \(\sim \) \(50\,\mathrm{ps} \) a portion of the decays falls into the material region and is discarded. Finally, a drop of sensitivity is expected for LLPs with a lifetime close to the \(b \) hadron lifetimes, where the contamination from \({b} {\overline{{b}}} \) events becomes important, especially for low mass LLPs.
8 Conclusion
A search for Higgs-like bosons decaying into two long-lived particles decaying hadronically has been carried out using data from pp collisions at 7\(\mathrm {\,Te V}\) collected with the LHCb detector, corresponding to a total integrated luminosity of \(0.62\) \( \text{ fb }^{-1}\).
The model used to describe the LLP decay is an mSUGRA process in which the lightest neutralino \(\tilde{\chi }^{0}_{1} \) decays through a baryon number violating coupling to three quarks. Upper limits have been placed on the production cross-section for Higgs-like boson masses from 80 to 140 \({\mathrm {Ge V\!/}c^2}\), LLP masses in the range 20–60 \({\mathrm {Ge V\!/}c^2}\), and LLP lifetimes in the range of 5–100 \(\mathrm{ps}\). The number of candidates is determined by the di-LLP invariant mass fit with signal templates inferred from simulation, and background estimates from data. For the explored parameter space of the theory all results, which are correlated, are consistent with zero. Upper limits at 95% CL for cross-section times branching ratio of 1 to 5\(\mathrm pb\) are inferred for most of the considered parameter range. They are below 2\(\mathrm pb\) for the decay of a 125 \({\mathrm {Ge V\!/}c^2}\) Higgs-like particle in two LLPs with mass in the 48–60 \({\mathrm {Ge V\!/}c^2}\) range and 10 \(\mathrm{ps}\) lifetime.
References
S. Dimopoulos, S. Raby, F. Wilczek, Supersymmetry and the scale of unification. Phys. Rev. D 24, 1681 (1981). doi:10.1103/PhysRevD.24.1681
S.P. Martin, A supersymmetry primer (2016). arXiv:hep-ph/9709356. doi:10.1142/9789812839657_0001
M.J. Strassler, K.M. Zurek, Discovering the Higgs through highly-displaced vertices. Phys. Lett. B 661, 263 (2008). doi:10.1016/j.physletb.2008.02.008. arXiv:hep-ph/0605193
D0 Collaboration, V.M. Abazov et al., Search for resonant pair production of neutral long-lived particles decaying to \(b\overline{b}\) in \(p\overline{p}\) collisions at \(\sqrt{s}=\)1.96 TeV. Phys. Rev. Lett. 103, 071801 (2009). arXiv:0906.1787. doi:10.1103/PhysRevLett.103.071801
CDF Collaboration, T. Aaltonen et al., Search for heavy metastable particles decaying to jet pairs in \(p\overline{p}\) collisions at \(\sqrt{s}=\)1.96 TeV, Phys. Rev. D 85, 012007 (2012). arXiv:1109.3136. doi:10.1103/PhysRevD.85.012007
CMS Collaboration, V. Khachatryan et al., Search for long-lived neutral particles decaying to quark-antiquark pairs in proton-proton collisions at \(\sqrt{s}=8\rm TeV\). Phys. Rev. D 91, 012007 (2015). arXiv:1411.6530. doi:10.1103/PhysRevD.91.012007
ATLAS Collaboration, G. Aad et al., Search for pair-produced long-lived neutral particles decaying in the ATLAS hadronic calorimeter in \(pp\) collisions at \(\sqrt{s}\) = 8 TeV. Phys. Lett. B 743, 15 (2015). arXiv:1501.04020. doi:10.1016/j.physletb.2015.02.015
ATLAS 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). arXiv:1109.2242. doi:10.1016/j.physletb.2011.12.057
LHCb Collaboration, R. Aaij et al., Search for long-lived particles decaying to jet pairs. Eur. Phys. J. C 75, 152 (2014). arXiv:1412.3021. doi:10.1140/epjc/s10052-015-3344-6
L.M. Carpenter, D.E. Kaplan, E.-J. Rhee, Six-quark decays of the Higgs boson in supersymmetry with R-parity violation. Phys. Rev. Lett. 99, 211801 (2007). doi:10.1103/PhysRevLett.99.211801. arXiv:hep-ph/0607204
D.E. Kaplan, K. Rehermann, Proposal for Higgs and superpartner searches at the LHCb experiment. JHEP 10, 056 (2007). doi:10.1088/1126-6708/2007/10/056. arXiv:0705.3426
P.W. Graham, D.E. Kaplan, S. Rajendran, P. Saraswat, Displaced supersymmetry. JHEP 07, 149 (2012). doi:10.1007/JHEP07(2012)149. arXiv:1204.6038
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. doi:10.1016/j.physletb.2012.08.020
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. doi:10.1016/j.physletb.2012.08.021
LHCb Collaboration, A.A. Alves Jr. et al., The LHCb detector at the LHC. JINST 3, S08005 (2008). doi:10.1088/1748-0221/3/08/S08005
LHCb Collaboration, R. Aaij et al., LHCb detector performance. Int. J. Mod. Phys. A 30, 1530022 (2015). arXiv:1412.6352. doi:10.1142/S0217751X15300227
R. Aaij et al., The LHCb trigger and its performance in 2011. JINST 8, P04022 (2013). doi:10.1088/1748-0221/8/04/P04022. arXiv:1211.3055
T. Sjöstrand, S. Mrenna, P. Skands, PYTHIA 6.4 physics and manual. JHEP 05, 026 (2006). doi:10.1088/1126-6708/2006/05/026. arXiv:hep-ph/0603175
J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD analysis. JHEP 07, 012 (2002). doi:10.1088/1126-6708/2002/07/012. arXiv:hep-ph/0201195
Geant4 Collaboration, J. Allison et al., Geant4 developments and applications. IEEE Trans. Nucl. Sci. 53, 270 (2006). doi:10.1109/TNS.2006.869826
Geant4 Collaboration, S. Agostinelli et al., Geant4: a simulation toolkit. Nucl. Instrum. Methods A 506, 250 (2003). doi:10.1016/S0168-9002(03)01368-8
M. Clemencic et al., The LHCb simulation application, Gauss: design, evolution and experience. J. Phys. Conf. Ser. 331, 032023 (2011). doi:10.1088/1742-6596/331/3/032023
LHCb Collaboration, R. Aaij et al., Search for the rare decay \(K^0_S \rightarrow \mu ^{+} \mu ^{-}\). JHEP 01, 090 (2013). arXiv:1209.4029. doi:10.1007/JHEP01(2013)090
M. Kucharczyk, P. Morawski, M. Witek, Primary vertex reconstruction at LHCb, CERN-LHCb-PUB-2014- 044, LHCb-PUB-2014-044 (2014). http://cdsweb.cern.ch/search?p=LHCb-PUB-2014-044&f=reportnumber&action_search=Search&c=LHCb+Notes
LHCb Collaboration, R. Aaij et al., Measurement of \(J/\psi \) production in \(pp\) collisions at \(\sqrt{s}=7\) TeV. Eur. Phys. J. C 71, 1645 (2011). arXiv:1103.0423. doi:10.1140/epjc/s10052-011-1645-y
LHCb Collaboration, R. Aaij et al., Measurement of \(\sigma (pp \rightarrow b\overline{b} X)\) at \(\sqrt{s}=7\) TeV in the forward region. Phys. Lett. B 694, 209 (2010). arXiv:1009.2731. doi:10.1016/j.physletb.2010.10.010
C. Oleari, The POWHEG BOX. Nucl. Phys. B-Proc. Suppl. 205–206, 36 (2010). doi:10.1016/j.nuclphysbps.2010.08.016
LHCb Collaboration, M. Williams et al., The HLT2 topological lines, CERN-LHCb-PUB-2011-002, LHCb-PUB-2011-002 (2011). http://cdsweb.cern.ch/search?p=LHCb-PUB-2011-002&f=reportnumber&action_search=Search&c=LHCb+Notes
ATLAS Collaboration, G. Aad et al., Measurement of the inclusive isolated prompt photon cross section in \(pp\) collisions at \(\sqrt{s}=7\rm TeV\) with the ATLAS detector. Phys. Rev. D 83, 052005 (2011). arXiv:1012.4389. doi:10.1103/PhysRevD.83.052005
R. Aaij et al., Performance of the LHCb vertex locator. JINST 9, P09007 (2014). doi:10.1088/1748-0221/9/09/P09007. arXiv:1405.7808
M. Botje et al., The PDF4LHC working group interim recommendations (2011). arXiv:1101.0538
D. Neill, I.Z. Rothstein, V. Vaidya, The Higgs transverse momentum distribution at NNLL and its theoretical errors. JHEP 12, 097 (2015). doi:10.1007/JHEP12(2015)097. arXiv:1503.00005
LHCb Collaboration, R. Aaij et al., Precision luminosity measurements at LHCb. JINST 9, P12005 (2014). arXiv:1410.0149. doi:10.1088/1748-0221/9/12/P12005
A.L. Read, Presentation of search results: the CL\(_s\) technique. J. Phys. G 28, 2693 (2002). doi:10.1088/0954-3899/28/10/313
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); NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); FOM and NWO (The Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MinES and FANO (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), IFIN-HH (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, XuntaGal and GENCAT (Spain), Herchel Smith Fund, The Royal Society, Royal Commission for the Exhibition of 1851 and the Leverhulme Trust (United Kingdom).
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendices
Fully simulated signal datasets
Table 5 shows the parameters used to generate the 11 fully simulated signal models with Pythia 6. The Higgs-like boson is produced by gluon-gluon fusion. In the table M\(_1\) corresponds to the Pythia parameter RMSS(1), and \(\tan {\beta }\) to RMSS(5). In addition, M\(_2\) (RMSS(2)) is set at 250 \({\mathrm {Ge V\!/}c^2}\) and \(\mu \) (RMSS(4)) has the value 140. A \(m_\mathrm{h^0}\) value of 125 \({\mathrm {Ge V\!/}c^2}\) requires RMSS(16) = 2300.
Detection efficiencies
Table 6 gives the detection efficiency as a function of \(m_\mathrm{h^0}\) and \(m_\mathrm{{LLP}}\), the LLP lifetime is 10 \(\mathrm{ps}\). Table 7 gives the efficiency as a function of \(m_\mathrm{{LLP}}\) and \(\tau _\mathrm{{LLP}}\), assuming \(m_\mathrm{h^0} =114\,{\mathrm {Ge V\!/}c^2} \).
Cross-section upper limits tables
Expected and observed 95% CL cross-section times branching ratio upper limits for benchmark models, from the fast simulation. Tables 8 and 9 give the limits as a function of \(m_\mathrm{h^0} \), covering LLP masses from 35 to 60 \({\mathrm {Ge V\!/}c^2}\), \(\tau _\mathrm{{LLP}} =10\,\mathrm{ps} \). Table 10: limits as a function of the LLP lifetime for \(m_\mathrm{h^0} =100\,{\mathrm {Ge V\!/}c^2} \) and \(m_\mathrm{{LLP}} =40\,{\mathrm {Ge V\!/}c^2} \), and for \(m_\mathrm{h^0} =125\,{\mathrm {Ge V\!/}c^2} \) and \(m_\mathrm{{LLP}} =48\,{\mathrm {Ge V\!/}c^2} \). Table 11: limits as a function of the LLP mass, for \(m_\mathrm{h^0} =125\,{\mathrm {Ge V\!/}c^2} \), \(\tau _\mathrm{{LLP}} =10\,\mathrm{ps} \).
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.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
About this article
Cite this article
Aaij, R., Adeva, B., Adinolfi, M. et al. Search for Higgs-like bosons decaying into long-lived exotic particles. Eur. Phys. J. C 76, 664 (2016). https://doi.org/10.1140/epjc/s10052-016-4489-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-016-4489-7