The observed limit contour in the SQUARK-NEUTRALINO plane.

The expected limit contour in the SQUARK-NEUTRALINO plane.

Cross section x BR limit for the UED model.

The transverse momentum distribution of the 4th-leading lepton for events with at least 4 leptons and no Z-boson candidate.

Distribution of missing transverse momentum for events with at least 4 leptons and no Z-boson candidate.

Distribution of effective mass for events with at least 4 leptons and no Z-boson candidate.

Vertex-level efficiency if no re-tracking is performed, as a function of the vertex radial position for an RPV SUSY model of squark production with $\tilde{q}\to q[\tilde{\chi}_1^0\to\mu qq]$, where $m(\tilde{q}) = 700$ GeV, $m(\tilde{\chi}_1^0) = 494$ GeV and $c\tau(\tilde{\chi}_1^0)$ = 175 mm. The result with re-tracking is tabulated in http://hepdata.cedar.ac.uk/view/ins1362183/d1.

Upper limits (95% CL) on the number of DV+muon displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q[\tilde{\chi}_1^0\to\mu qq]$ or $\tilde{g}\to qq[\tilde{\chi}_1^0\to\mu qq]$ decays, respectively.

Upper limits (95% CL) on the number of DV+muon displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q\tilde{\chi}_1^0$ or $\tilde{g}\to qq\tilde{\chi}_1^0$ decays, respectively.

Upper limits (95% CL) on the number of DV+electron displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7c) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q[\tilde{\chi}_1^0\to e qq]$ or $\tilde{g}\to qq[\tilde{\chi}_1^0\to e qq]$ decays, respectively.

Upper limits (95% CL) on the number of DV+electron displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7d) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q\tilde{\chi}_1^0$ or $\tilde{g}\to qq\tilde{\chi}_1^0$ decays, respectively.

Upper limits (95% CL) on the number of dilepton $e^+e^-$ displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to ee\nu]$ decays and masses in GeV as indicated. A dash indicates results omitted due to limited statistical precision.

Upper limits (95% CL) on the number of dilepton $e^{\pm}\mu^{\mp}$ displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu]$ decays and masses in GeV as indicated.

Upper limits (95% CL) on the number of $\mu^+\mu^-$ displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5c) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to \mu\mu\nu]$ decays and masses in GeV as indicated. A dash indicates results omitted due to limited statistical precision.

Upper limits (95% CL) on the number of dilepton ($e^+e^-+e^{\pm}\mu^{\mp}+\mu^+\mu^-$) displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5d) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and $m(\tilde{\chi}_1^0)$ = 400 GeV.

Upper limits (95% CL) on the number of dilepton ($e^+e^-+e^{\pm}\mu^{\mp}+\mu^+\mu^-$) displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5d) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and $m(\tilde{\chi}_1^0)$ = 1000 GeV.

Upper limits (95% CL) from the dilepton $e^+e^-+e^{\pm}\mu^{\mp}$ channels on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 6a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu/ee\nu]$ decays and masses in GeV as indicated. For comparison, the production cross-sections for $m(\tilde{g})$ = 600 GeV and 1300 GeV are $1280\pm230$ fb and $1.7\pm0.7$ fb, respectively.

Upper limits (95% CL) from the dilepton $e^{\pm}\mu^{\mp}+\mu^+\mu^-$ channels on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 6b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu/\mu\mu\nu]$ decays and masses in GeV as indicated. For comparison, the production cross-sections for $m(\tilde{g})$ = 600 GeV and 1300 GeV are $1280\pm230$ fb and $1.7\pm0.7$ fb, respectively. A dash indicates results omitted due to limited statistical precision.

Upper limits (95% CL) from the dilepton $e^+e^-+e^{\pm}\mu^{\mp}+\mu^+\mu^-$ channels on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 6c) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and $m(\tilde{\chi}_1^0)$ = 400 GeV. For comparison, the production cross-section for $m(\tilde{g})$ = 1100 GeV is $7.6\pm2.8$ fb.

Upper limits (95% CL) from the dilepton $e^+e^-+e^{\pm}\mu^{\mp}+\mu^+\mu^-$ channels on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 6c) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and $m(\tilde{\chi}_1^0)$ = 1000 GeV. For comparison, the production cross-section for $m(\tilde{g})$ = 1100 GeV is $7.6\pm2.8$ fb.

Upper limits (95% CL) from the DV+jets channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 8a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark pair production, with $\tilde{q}\to q\tilde{\chi}_1^0$ decays, $m(\tilde{q})$ = 700 GeV and $m(\tilde{\chi}_1^0)$ = 494 GeV. For comparison, the production cross-section for $m(\tilde{q})$ = 700 GeV is $124\pm17$ fb. A dash indicates where the limit-setting procedure did not converge.

Upper limits (95% CL) from the DV+jets channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 8b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark pair production, with $\tilde{q}\to q\tilde{\chi}_1^0$ decays, $m(\tilde{q})$ = 1000 GeV and $m(\tilde{\chi}_1^0)$ = 108 GeV. For comparison, the production cross-section for $m(\tilde{q})$ = 1000 GeV is $11.9\pm1.5$ fb. A dash indicates where the limit-setting procedure did not converge.

Upper limits (95% CL) from the DV+$E_T^{miss}$ channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 9a) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and a $\tilde{\chi}_1^0$ mass in GeV as indicated. For comparison, the production cross-section for $m(\tilde{g})$ = 1100 GeV is $7.6\pm2.8$ fb.

Upper limits (95% CL) from the DV+jets channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 9b) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and a $\tilde{\chi}_1^0$ mass in GeV as indicated. For comparison, the production cross-section for $m(\tilde{g})$ = 1100 GeV is $7.6\pm2.8$ fb. A dash indicates where the limit-setting procedure did not converge.

Upper limits (95% CL) from the DV+$E_T^{miss}$ channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 10a) for four split-SUSY models as a function of the $\tilde{g}$ proper decay distance $c\tau$. The models consider gluino pair production, with $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays, $\tilde{g}$ masses in GeV as indicated and $m(\tilde{\chi}_1^0)$ = 100 GeV. For comparison, the production cross-sections for $m(\tilde{g})$ = 400 GeV, 800 GeV and 1400 GeV are $18700\pm2800$ fb, $149\pm31$ fb and $0.71\pm0.32$ fb, respectively.

Upper limits (95% CL) from the DV+jets channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 10b) for four split-SUSY models as a function of the $\tilde{g}$ proper decay distance $c\tau$. The models consider gluino pair production, with $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays, $\tilde{g}$ masses in GeV as indicated and $m(\tilde{\chi}_1^0)$ = 100 GeV. For comparison, the production cross-sections for $m(\tilde{g})$ = 400 GeV, 800 GeV and 1400 GeV are $18700\pm2800$ fb, $149\pm31$ fb and $0.71\pm0.32$ fb, respectively. A dash indicates where the limit-setting procedure did not converge or the limit is greater than $10^8$ fb.

Upper limits (95% CL) from the DV+$E_T^{miss}$ channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 10c) for three split-SUSY models as a function of the $\tilde{g}$ proper decay distance $c\tau$. The models consider gluino pair production, with $[\tilde{g}\to tt \tilde{\chi}_1^0]$ decays, $\tilde{g}$ masses in GeV as indicated and $m(\tilde{\chi}_1^0)$ = 100 GeV. For comparison, the production cross-sections for $m(\tilde{g})$ = 600 GeV, 1000 GeV and 1400 GeV are $1270\pm230$ fb, $22\pm6$ fb and $0.71\pm0.32$ fb, respectively.

Upper limits (95% CL) from the DV+jets channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 10d) for three split-SUSY models as a function of the $\tilde{g}$ proper decay distance $c\tau$. The models consider gluino pair production, with $[\tilde{g}\to tt \tilde{\chi}_1^0]$ decays, $\tilde{g}$ masses in GeV as indicated and $m(\tilde{\chi}_1^0)$ = 100 GeV. For comparison, the production cross-sections for $m(\tilde{g})$ = 600 GeV, 1000 GeV and 1400 GeV are $1270\pm230$ fb, $22\pm6$ fb and $0.71\pm0.32$ fb, respectively. A dash indicates where the limit-setting procedure did not converge.

Upper limits (95% CL) from the DV+$E_T^{miss}$ channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 10e) for four split-SUSY models as a function of the $\tilde{g}$ proper decay distance $c\tau$. The models consider gluino pair production, with $[\tilde{g}\to tt \tilde{\chi}_1^0]$ decays, $\tilde{g}$ masses in GeV as indicated and $m(\tilde{\chi}_1^0)$ = $m(\tilde{g}) - 480$ GeV. For comparison, the production cross-sections for $m(\tilde{g})$ = 500 GeV, 800 GeV and 1400 GeV are $4450\pm700$ fb, $149\pm31$ fb and $0.71\pm0.32$ fb, respectively.

Upper limits (95% CL) from the DV+jets channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 10f) for four split-SUSY models as a function of the $\tilde{g}$ proper decay distance $c\tau$. The models consider gluino pair production, with $[\tilde{g}\to tt \tilde{\chi}_1^0]$ decays, $\tilde{g}$ masses in GeV as indicated and $m(\tilde{\chi}_1^0)$ = $m(\tilde{g}) - 480$ GeV. For comparison, the production cross-sections for $m(\tilde{g})$ = 500 GeV, 800 GeV and 1400 GeV are $4450\pm700$ fb, $149\pm31$ fb and $0.71\pm0.32$ fb, respectively. A dash indicates where the limit-setting procedure did not converge.

Observed 95% CL exclusion contour in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY model with gluino production, $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays and $m(\tilde{\chi}_1^0)$ = 100 GeV. These results are obtained from the DV+$E_T^{miss}$ and DV+jets channels, see http://hepdata.cedar.ac.uk/view/ins1362183/d34 for details.

Expected 95% CL exclusion contour in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY model with gluino production, $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays and $m(\tilde{\chi}_1^0)$ = 100 GeV. These results are obtained from the DV+$E_T^{miss}$ and DV+jets channels, see http://hepdata.cedar.ac.uk/view/ins1362183/d34 for details.

Observed 95% CL exclusion contour in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY model with gluino production, $[\tilde{g}\to tt \tilde{\chi}_1^0]$ decays and $m(\tilde{\chi}_1^0)$ = 100 GeV. These results are obtained from the DV+$E_T^{miss}$ and DV+jets channels, see http://hepdata.cedar.ac.uk/view/ins1362183/d34 for details.

Expected 95% CL exclusion contour in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY model with gluino production, $[\tilde{g}\to tt \tilde{\chi}_1^0]$ decays and $m(\tilde{\chi}_1^0)$ = 100 GeV. These results are obtained from the DV+$E_T^{miss}$ and DV+jets channels, see http://hepdata.cedar.ac.uk/view/ins1362183/d34 for details.

Observed 95% CL exclusion contour in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY model with gluino production, $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays and $m(\tilde{\chi}_1^0)$ = $m(\tilde{g})$ - 480 GeV. These results are obtained from the DV+$E_T^{miss}$ and DV+jets channels, see http://hepdata.cedar.ac.uk/view/ins1362183/d34 for details.

Expected 95% CL exclusion contour in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY model with gluino production, $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays and $m(\tilde{\chi}_1^0)$ = $m(\tilde{g})$ - 480 GeV. These results are obtained from the DV+$E_T^{miss}$ and DV+jets channels, see http://hepdata.cedar.ac.uk/view/ins1362183/d34 for details.

Choice of signal region in the $m(\tilde{g})$-$c\tau$ plane for the split-SUSY models with gluino decays and neutralino masses in GeV as indicated. A value of 1 indicates that the DV+$E_T^{miss}$ channel provided the best expected constraint for that point, while a value of -1 indicates the DV+jets channel. A dash indicates where results were not calculated for a particular mass/lifetime combination.

Upper limits (95% CL) from the DV+$E_T^{miss}$ channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 11a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q[\tilde{\chi}_1^0\to \nu qq]$ or $\tilde{g}\to qq[\tilde{\chi}_1^0\to \nu qq]$ decays, respectively. For comparison, the production cross-sections for $m(\tilde{g})$ = 700 GeV, $m(\tilde{q})$ = 700 GeV and $m(\tilde{q})$ = 1000 GeV are $430\pm80$ fb, $124\pm17$ fb and $11.9\pm1.5$ fb, respectively.

Upper limits (95% CL) from the DV+$E_T^{miss}$ channel on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 11b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark pair production, with $\tilde{q}\to q\tilde{\chi}_1^0$ decays and masses in GeV as indicated. For comparison, the production cross-sections for $m(\tilde{q})$ = 700 GeV and 1000 GeV are $124\pm17$ fb and $11.9\pm1.5$ fb, respectively.

Vertex mass distributions in data and one RPV SUSY signal model, for vertices with 3, 4 and $\geq$5 tracks in the DV+muon channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for squark pair production, with $\tilde{q}\to q[\tilde{\chi}_1^0\to\mu qq]$ decays, $m(\tilde{q})$ = 700 GeV, $m(\tilde{\chi}_1^0)$ = 494 GeV and $c\tau(\tilde{\chi}_1^0)$ = 175 mm.

Vertex mass distributions in data and one GGM SUSY signal model, for vertices with 3, 4 and $\geq$5 tracks in the DV+jets channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV, $m(\tilde{\chi}_1^0)$ = 400 GeV and $c\tau(\tilde{\chi}_1^0)$ = 175 mm.

Vertex mass distributions in data and one RPV SUSY signal model, for vertices with 3, 4 and $\geq$5 tracks in the DV+electron channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for squark pair production, with $\tilde{q}\to q[\tilde{\chi}_1^0\to e qq]$ decays, $m(\tilde{q})$ = 700 GeV, $m(\tilde{\chi}_1^0)$ = 494 GeV and $c\tau(\tilde{\chi}_1^0)$ = 175 mm.

Vertex mass distributions in data and one split-SUSY signal model, for vertices with 3, 4 and $\geq$5 tracks in the DV+$E_T^{miss}$ channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for gluino pair production, with $[\tilde{g}\to g/qq \tilde{\chi}_1^0]$ decays, $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi}_1^0)$ = 100 GeV and $c\tau(\tilde{g})$ = 175 mm.

Vertex mass distributions in data and one RPV SUSY signal model, for vertices with 0, 1 and $\geq$2 charged leptons in the dilepton $e^+e^-$ channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu/ee\nu]$ decays, $m(\tilde{g})$ = 1300 GeV, $m(\tilde{\chi}_1^0)$ = 50 GeV and $c\tau(\tilde{\chi}_1^0)$ = 30 mm.

Vertex mass distributions in data and one RPV SUSY signal model, for vertices with 0, 1 and $\geq$2 charged leptons in the dilepton $e^{\pm}\mu^{\mp}$ channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu/ee\nu]$ decays, $m(\tilde{g})$ = 1300 GeV, $m(\tilde{\chi}_1^0)$ = 50 GeV and $c\tau(\tilde{\chi}_1^0)$ = 30 mm.

Vertex mass distributions in data and one RPV SUSY signal model, for vertices with 0, 1 and $\geq$2 charged leptons in the dilepton $\mu^+\mu^-$ channel. The entries are the number of DVs in each bin, not normalized according to bin size. The signal model is for gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu/\mu\mu\nu]$ decays, $m(\tilde{g})$ = 1300 GeV, $m(\tilde{\chi}_1^0)$ = 50 GeV and $c\tau(\tilde{\chi}_1^0)$ = 30 mm.

Details of sparticle decay modes, model parameters, approximate number of generated events $N_{\rm evts}$ and total production cross-section $\sigma$ for selected models considered in this analysis. The long-lived particle (LLP) is the $\tilde{\chi}_1^0$ in all cases except for the model used for DV+$E_T^{miss}$, where it is the $\tilde{g}$. Events in some models are filtered to make better use of computing resources, details are given in the original table and the associated efficiency is already taken into account in the quoted cross-sections.

Numbers of events satisfying each selection criterion in turn for the DV+muon channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Filter and trigger'' line includes the requirements placed on the primary vertex.

Numbers of events satisfying each selection criterion in turn for the DV+electron channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Filter and trigger'' line includes the requirements placed on the primary vertex.

Numbers of events satisfying each selection criterion in turn for the DV+jets channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Filter and trigger'' line includes the requirements placed on the primary vertex.

Numbers of events satisfying each selection criterion in turn for the DV+$E_T^{miss}$ channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Filter and trigger'' line includes the requirements placed on the primary vertex.

Numbers of events satisfying each selection criterion in turn for the dilepton $e^+e^-$ channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Trigger'' line includes the event-level cosmic-muon veto and the requirements placed on the primary vertex, but not the offline filter selection. This is included in the ''Lepton reconstruction and DV match'' line, along with all other requirements placed on offline-reconstructed leptons. About 20,000 events of the full sample are analyzed, which contain an equal mixture of $\tilde{\chi}_1^0 \to e^{\pm}\mu^{\mp}\nu$ and $\tilde{\chi}_1^0 \to e^+e^-\nu$ decays.

Numbers of events satisfying each selection criterion in turn for the dilepton $e^{\pm}\mu^{\mp}$ channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Trigger'' line includes the event-level cosmic-muon veto and the requirements placed on the primary vertex, but not the offline filter selection. This is included in the ''Lepton reconstruction and DV match'' line, along with all other requirements placed on offline-reconstructed leptons. The sample contains about 50% $\tilde{\chi}_1^0 \to e^{\pm}\mu^{\mp}\nu$ decays, 25% $\tilde{\chi}_1^0 \to e^+e^-\nu$ and 25% $\tilde{\chi}_1^0 \to \mu^+\mu^-\nu$ decays.

Numbers of events satisfying each selection criterion in turn for the dilepton $\mu^+\mu^-$ channel, using the model listed in http://hepdata.cedar.ac.uk/view/ins1362183/d44. The ''Trigger'' line includes the event-level cosmic-muon veto and the requirements placed on the primary vertex, but not the offline filter selection. This is included in the ''Lepton reconstruction and DV match'' line, along with all other requirements placed on offline-reconstructed leptons. About 20,000 events of the full sample are analyzed, which contain an equal mixture of $\tilde{\chi}_1^0 \to e^{\pm}\mu^{\mp}\nu$ and $\tilde{\chi}_1^0 \to \mu^+\mu^-\nu$ decays.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Despite listing this as an exclusive final state (as there must be at least one b-jets), there is no explicit selection on the presence of additional light-flavour jets. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 300, 500 and 800 GeV and $\tan\beta$ = 10 in the hMSSM scenario are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. In the paper, the first bin is cut off at 60 GeV for aesthetics but contains underflows down to 50 GeV as in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the Drell Yan production cross section times ditau branching fraction as a function of the Zprime boson mass.

Observed and expected 95% CL upper limits on the Higgs boson production cross section times ditau branching fraction as a function of the boson mass and the relative strength of the b-associated production.

Ratio of the 95% CL upper limits on the production cross section times branching fraction for alternate Zprime models with respect to the SSM, both observed and expected are shown.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance and acceptance times efficiency for a heavy gauge boson produced by Drell Yan as a function of the gauge boson mass.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV and $m_{\tilde g}=$1.6 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$1.6 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV with $m_{\tilde g}=$2.2 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$2.2 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV and $m_{\tilde g}=$3.0 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$3.0 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed and expected limits on the cross-section times branching ratio in fb at 95% CL for the gluino-mediated stop model with off-shell stops, for mass points with four- or five-body gluino decays, obtained by the SS/3L and 0/1L3B analyses.

Data (bold boxes) and background estimates (colour fill) for m_beta vs. m_betagamma for the gluino R-hadron search (1000 GeV). The blue thin-line boxes illustrate the expected signal (on top of background) for the given R-hadron mass hypothesis. The black dashed vertical/horizontal lines at 500 GeV show the mass selection (signal region in the top-right). Two events pass this selection.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of bottom-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of top-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Final selection requirements as a function of the simulated R-hadron mass.

Summary of all studied systematic uncertainties. Ranges indicate a dependency on the R-hadron mass hypothesis (from low to high masses).

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range after the final selection using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Distribution of the truth-level beta for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Distribution of the truth-level betagamma for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 50 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 30 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 10 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the gluino R-hadron search.

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the squark R-hadron search.

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range in the meta-stable gluino R-hadron search using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

The observed 95\% C.L. limit for $pp\rightarrow G^{*}_{KK}\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ in the bulk RS model with $k/\bar{M}_{Pl} = 2$, as a function of resonance mass.

The observed 95\% C.L. limit for $pp\rightarrow H\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ with fixed $\Gamma_{H} = 1$ GeV, as a function of resonance mass.

Excluded range of lifetimes as a function of gluino R-hadron mass. The central value of the lower-limit of the observed and expected excluded masses at the 95% CL is given.

Cross-section as a function of mass for stable gluino R-hadrons. The central value of the observed and expected 95% upper limit on excluded cross-section is given. The 1 sigma band on the expected UL is given as an uncertainty.