Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average core acoplanarity ($<\alpha^{core}>$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average invariant mass ($<m_{\mu\mu}>$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

The product of signal acceptance and efficiency in the 2l categories for the signal produced via vector boson fusion.

$m_{X}^{T}$ distribution in data in the 0l 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the non-VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis. The CMS search for a heavy resonance using 2016 data and the same final state [JHEP 11 (2018) 172] is shown as a comparison.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 2 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 3 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 4 TeV for the non-VBF signal.

Three-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for small-angle radiation ($\Delta R_{23}$ < 1.0)

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Three-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Three-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Z + two-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for small-angle radiation ($\Delta R_{23}$ < 1.0)

Z + two-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Z + two-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Z + two-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Distributions of $\Delta\phi_{\mathrm{trk,Z}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of $\Delta\phi_{\mathrm{trk,Z}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 70-90% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 50-70% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 70-90% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 50-70% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.

The number of background events obtained from the 1- and 2-leg predictions derived from data, together with the direct observation in data, in bins in ChF, where either the leading or subleading jet has a ChF within the bin edges, and both have a ChF below the upper bin threshold. The bottom panel shows the ratios of the observation in data to the 1-leg and the 2-leg background predictions.

The expected and observed 95% CL upper limits on the production cross section for SIMPs with masses between 1 and 1000 GeV, with the assumption that the SIMP interaction in the detector can be approximated as neutron-like. The theoretical prediction of a simplified model incorporating this approximation and including a scalar mediator with couplings g_chi = -1 and g_q = 1 is also shown (red line). For masses above 100 GeV, where the modelling of the SIMP-nucleon interaction becomes more speculative, the obtained cross section upper limits are increasingly uncertain (shaded area).

Exclusion limit for BrHZX_BrXee

Exclusion limit for BrHZX_BrXmumu

Exclusion limit for BrHZX_BrXll

Exclusion limit for cah

Exclusion limit for cZh

Exclusion limit for kappa