The invariant mass spectra of dimuon events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

Observed and expected background yields for different mass ranges in the dielectron channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the GSM model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the LR model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the E6 model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

The observed local p-value for the dielectron channel as a function of the dilepton invariant mass.

The observed local p-value for the dimuon channel as a function of the dilepton invariant mass.

The observed local p-value for the combined channel as a function of the dilepton invariant mass.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dielectron channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the combination of the dielectron and dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dielectron channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dimuon channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the combination of the dielectron and dimuon channels. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with central leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with forward leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Invariant mass resolution for dielectron pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dimuon pairs as obtained from simulated DY samples as a function of generated mass.

Observed data and total background estimate in all bins used in the limit calculation for the CI limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 6 bins in dilepton invariant mass: [400,500,700,1100,1900,3500,10000] in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the ADD limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 5 bins in dilepton invariant mass: [1800,2200,2600,3000,3400,10000] ([1900,2200,2600,3000,3400,10000] for 2016) in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the CI limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the ADD limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Higgs fiducial cross section in bins of rapidity for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of Jet Multiplicity The first uncertainty is statistical, the second is systematic uncertainty.

Higgs fiducial cross section in bins of pT of leading jet. The first uncertainty is statistical, the second is the systematic uncertainty. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb. The 0-jet bin is also included here, which is not shown in the Figure.

Results of likelihood scans for the signal strength modifiers corresponding to the five main SM H boson production mechanisms. The uncertainties, corresponding to one standard deviation confidence intervals, include both statistical and systematic sources. The additional breakdown of the uncertainties into their separate statistical and systematic contributions is also shown.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 68% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 95% CL region.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the Stage 0 STXS production bins and the inclusive measurement at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the merged Stage 1.2 STXS production bins at mH= 125.38 GeV.

Observed correlation matrix between the measured cross sections for the Stage 0 for HZZ decay.

Observed correlation matrix between the measured cross sections for the Stage 1.2 for HZZ decay.

Results of the maximal merging scheme STXS fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Observed and expected correlations between the parameters in the STXS stage 1.2 maximal merging fit.

Results of the minimal merging scheme STXS fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Observed and expected correlations between the parameters in the STXS stage 1.2 minimal merging fit.

Observed likelihood surface.

Expected likelihood surface.

Observed likelihood surface.

Expected likelihood surface.

Results of the STXS stage 0 fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Observed and expected correlations between the parameters in the STXS stage 0 fit.

Observed likelihood surface.

Expected likelihood surface.

Observed likelihood surface.

Expected likelihood surface.

The product of acceptance and selection efficiency within the fiducial region for the VBF $\mathrm{H}^{\pm\pm}\rightarrow\mathrm{W}^{\pm}\mathrm{W}^{\pm}\rightarrow 2\ell 2v$ and $\mathrm{H}^{\pm}\rightarrow\mathrm{W}^{\pm}\mathrm{Z}\rightarrow 3\ell v$ processes, as a function of $m_{\mathrm{H}_{5}}$. The combination of the statistical and systematic uncertainties is shown. The theoretical uncertainties in the acceptance are also included.

Expected and observed exclusion limits at 95\% confidence level for $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm\pm}) \mathcal{B}(\mathrm{H}^{\pm\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{W}^{\pm})$ as functions of $m_{\mathrm{H}^{\pm\pm}}$. The contribution of the $\mathrm{H}^{\pm}$ boson signal is set to zero for the derivation of the exclusion limits on the $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm\pm}) \mathcal{B}(\mathrm{H}^{\pm\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{W}^{\pm})$. Values above the curves are excluded.

Expected and observed exclusion limits at 95\% confidence level for $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm}) \mathcal{B}(\mathrm{H}^{\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{Z})$ as functions of $m_{\mathrm{H}^{\pm}}$. The contribution of the $\mathrm{H}^{\pm\pm}$ boson signal is set to zero for the derivation of the exclusion limits on the $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm}) \mathcal{B}(\mathrm{H}^{\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{Z})$. Values above the curves are excluded.

Expected and observed exclusion limits at 95\% confidence level for $s_{\mathrm{H}}$ as functions of $m_{\mathrm{H}_{5}}$ in the Georgi-Machacek model. Values above the curves are excluded. The exclusion limits for $s_{\mathrm{H}}$ are shown up to $m_{\mathrm{H}_{5}}=2000$ GeV, given the low sensitivity in the Georgi-Machacek model for values above that mass.

Covariance matrix for all the bins used in Figure 4. The covariance matrix from the background-only fit is shown.

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

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.

95% CL observed upper limits on the Yukawa couplings

Distribution of the transverse momentum of the diphoton system for the $\mathrm{Z}\gamma\gamma$ muon channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{W}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{Z}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

The measured $\mathrm{W}\gamma\gamma$ and $\mathrm{Z}\gamma\gamma$ fiducial cross sections and corresponding theoretical predictions from MadGraph5_aMC@NLO. The measured yields in the electron and muon channels are extrapolated to a common fiducial phase space determined from simulated signal events at the generated particle level. Generated particles are considered stable if their mean decay length is larger than $1\,\mathrm{cm}$. Generated leptons are required to have a $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. The momenta of photons in a cone of $\Delta R=0.1$, the same cone size as the one applied to reconstructed data, are added to the charged lepton momentum to correct for final-state radiation. Generated photons are required to have $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. Additionally, the candidate photons are required to have no selected leptons or photons in a cone of radius $\Delta R=0.4$ and no other stable particles, apart from photons and neutrinos, in a cone of radius $\Delta R=0.1$. Events are then selected in the $\mathrm{W}\gamma\gamma$ channel by requiring exactly one electron (muon) with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$ and at least two photons with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Events are selected in the $\mathrm{Z}\gamma\gamma$ channel by requiring two electrons (muons), at least one of them with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$, and not less than two photons, each of them with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Additionally, the invariant mass of the dilepton system is required to be $m_{\ell\ell}>55\,\mathrm{GeV}$.

Expected and observed 95% confidence level intervals for the different anomalous couplings parameters in both the $\ell\nu\gamma\gamma$ and $\ell\ell\gamma\gamma$ channels. All parameters are fixed to their SM values except the one that is fitted. No unitarity regularisation scheme is applied. Intervals with a NAN value are not measured in the corresponding channels.

Differential cross sections as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Differential cross sections as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events.

Area normalized distribution as function of Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events

Response matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events.

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).