The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $55 < p_{\mathrm{T}} < 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the FAIL region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

Mass spectrum predicted in the signal region defined by $G_{\text{i}}^{\text{Strips}} > 0.22$ and $p_{\mathrm{T}} > 70$ GeV.

Cross section limits for gluino and supersymmetric top R-hadrons for both background prediction methods.

Cross section limits for supersymmetric tau models for both background prediction methods.

Cross section limits for DY-produced tau prime models (single and multicharged) for both background prediction methods.

Cross section limits for Z prime to multicharged tau prime model for both background prediction methods. All Z prime models assume a narrow width, and a 100% branching fraction to 600 GeV tau primes.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the ionization method.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the mass method.

Cumulative selection efficiency for the data and for two signal hypotheses.

Expected and observed mass limits obtained using 2017-2018 data for various HSCP candidate models,for the two background estimate methods.

Mass windows used in the mass method as a function of the signal target mass for the signal samples assuming a charge of $1e$ and $2e$, as appropriate to the signal model.

Mass spectrum predicted in the validation region defined by $0.018 < G_{\text{i}}^{\text{Strips}} < 0.057$ and $p_{\mathrm{T}} > 70$ GeV. The thresholds used in the $G_{\text{i}}^{\text{Strips}}$ requirement represent the 50% and 90% quantile of the distribution.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency as a function of $\beta$, for the $\tilde{g}$ signals.

Trigger efficiency as a function of $\beta$, for the $\tilde{\tau}$ signals.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{\tau}$ obtained with the ionization method.

Cross section limits for $\tilde{\tau}$ obtained with the mass method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 modelobtained with the ionization method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 model obtained with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the mass method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the ionization method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the mass method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the ionization method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the mass method.

Expected background yield, expected signal yield, and observed data for the mass method for 2017.

Expected background yield, expected signal yield, and observed data for the mass method for 2018.

Selection efficieny for GMSB supersymmetric $\tau$.

Selection efficiency for gluino R-hardon samples.

Selection efficiency for supersymmetric top R-hadrons.

Selection efficiency for pair-produced supersymmetric $\tau$.

Selection efficiency for $\tau'$ with $|Q| = 1e$.

Selection efficiecny for $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass = 3 TeV) going to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass 5 TeV) goes to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 600 GeV) with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 1 TeV) with charge $|Q| = 2e$.

Observed profile likelihood projection on mH, for the 2e2mu final state, using the N-2D′ VXBS approach. Both statistical and systematic uncertainties have been considered.

Observed profile likelihood projection on mH, for the 2mu2e final state, using the N-2D′ VXBS approach. Both statistical and systematic uncertainties have been considered.

Result of the Higgs boson width measurement using the on-shell production only, shown as 1 − CL curve extracted using the Feldman–Cousins approach (black dots). The error bars represent the spread of the simulated experiments for each point.

Observed profile likelihood projections from the Higgs boson width fit using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Expected profile likelihood projections from the Higgs boson width fit using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Observed profile likelihood projections from the Higgs boson width fit using the on- and off-shell H → ZZ → 4l channel.

Expected (dashed) profile likelihood projections from the Higgs boson width fit using the on- and off-shell H → ZZ → 4l channel.

Observed 2D profile likelihood projection of the off-shell signal strength parameters (mu off-shell, mu off-shell) from the fit to the combined off-shell H → ZZ → 4l and 2l2nu channels. [Nature Phys. 18 (2022) 1329] The best fit value is shown by the black cross and the SM prediction by the red x. The 68% and 95%CL contours are given by the dashed and solid curves, respectively. The color scale to the right of the plot relates the quantitative values.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2016 pre-VPF, for positive muons.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2016 pre-VPF, for negative muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2016 pre-VPF, for positive muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2016 pre-VPF, for negative muons.

Difference between data and simulation in Z → 2e events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) ranges in 2016 pre-VPF.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2016 post-VPF, for positive muons.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2016 post-VPF, for negative muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2016 post-VPF, for positive muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2016 post-VPF, for negative muons.

Difference between data and simulation in Z → 2e events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) ranges in 2016 post-VPF.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2017, for positive muons.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2017, for negative muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2017, for positive muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2017, for negative muons.

Difference between data and simulation in Z → 2e events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) ranges in 2017.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2018, for positive muons.

Difference between data and simulation in Z → 2mu events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) range, in 2018, for negative muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2018, for positive muons.

Difference between data and simulation in J/Psi → 2mu events, normalized to simulation, as a function of the transverse momentum (pt) for two different pseudorapidity (|eta|) ranges, in 2018, for negative muons.

Difference between data and simulation in Z → 2e events, normalized to simulation, as a function of transverse momentum (pt) for three different pseudorapidity (|eta|) ranges in 2018.

The observed likelihood scan as a function of mH using the N-2D'_VXBS$ approach.

The expected likelihood scan as a function of mH using the N-2D'_VXBS$ approach. Expected likelihood has been obtained with the hypothesis mH = 125.38 GeV

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2e events in data in 2016 pre-VPF.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2mu events in data 2016 pre-VPF.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2e events in data in 2016 post-VPF.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2mu events in data 2016 post-VPF.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2e events in data in 2017.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2mu events in data 2017.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2e events in data in 2018.

Comparison of measured with the predicted dilepton mass resolution using the event-by-event mass uncertainty for Z → 2mu events in data 2018.

Observed profile likelihood projections from the Higgs boson width fit using the on- and off-shell H → ZZ → 4l channel with kappa Q constrained to 0.

Expected profile likelihood projections from the Higgs boson width fit using the on- and off-shell H → ZZ → 4l channel with kappa Q constrained to 0.

Observed profile likelihood projection on the Higgs boson width using on- and off-shell production. The analysis of H → ZZ → 4l, considering an unconstrained coupling strength κQ is presented.

Expected profile likelihood projection on the Higgs boson width using on- and off-shell production. The analysis of H → ZZ → 4l, considering an unconstrained coupling strength κQ is presented.

Expected profile likelihood projection on muF off-shell using on- and off-shell production for H → ZZ → 4l.

Observed profile likelihood projection on muF off-shell using on- and off-shell production for H → ZZ → 4l.

Expected profile likelihood projection on muF off-shell using using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Observed profile likelihood projection on muF off-shell using using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Expected profile likelihood projection on muV off-shell using on- and off-shell production for H → ZZ → 4l.

Observed profile likelihood projection on muV off-shell using on- and off-shell production for H → ZZ → 4l.

Expected profile likelihood projection on muV off-shell using using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu.

Observed profile likelihood projection on muV off-shell using using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Expected profile likelihood projection on mu off-shell using on- and off-shell production for H → ZZ → 4l.

Observed profile likelihood projection on mu off-shell using on- and off-shell production for H → ZZ → 4l.

Expected profile likelihood projection on mu off-shell using using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Observed profile likelihood projection on mu off-shell using using the full combination of H → ZZ → 4l with the off-shell H → ZZ → 2l2nu [Nature Phys. 18 (2022) 1329].

Covariance for full matrix measurement all bins from $m(t\bar{t})$ fit

Results for full matrix measurement inclusive from $p_{T}(t)$

Covariance for full matrix measurement inclusive from $p_{T}(t)$

Results for full matrix measurement all bins from $p_{T}(t)$ fit

Covariance for full matrix measurement all bins from $p_{T}(t)$ fit

Results for full matrix measurement differential in $m(t\bar{t})$

Covariance for full matrix measurement differential in $m(t\bar{t})$

Results for full matrix measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Covariance for full matrix measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Results for full matrix measurement differential in $p_{T}(t)$

Covariance for full matrix measurement differential in $p_{T}(t)$

Results for $D$ measurement inclusive from $m(t\bar{t})$

Covariance for $D$ measurement inclusive from $m(t\bar{t})$

Results for $D$ measurement inclusive from $p_{T}(t)$

Covariance for $D$ measurement inclusive from $p_{T}(t)$

Results for $\tilde{D}$ measurement inclusive in $m(t\bar{t})$

Covariance for $\tilde{D}$ measurement inclusive in $m(t\bar{t})$

Results for $\tilde{D}$ measurement inclusive in $p_{T}(t)$

Covariance for $\tilde{D}$ measurement inclusive in $p_{T}(t)$

Results for $D$ measurement all bins from from $m(t\bar{t})$ fit

Covariance for $D$ measurement all bins from from $m(t\bar{t})$ fit

Results for $D$ measurement all bins from from $p_{T}(t)$ fit

Covariance for $D$ measurement all bins from from $p_{T}(t)$ fit

Results for $\tilde{D}$ measurement all bins from from $m(t\bar{t})$ fit

Covariance for $\tilde{D}$ measurement all bins from from $m(t\bar{t})$ fit

Results for $\tilde{D}$ measurement all bins from from $p_{T}(t)$ fit

Covariance for $\tilde{D}$ measurement all bins from from $p_{T}(t)$ fit

Results for $D$ measurement differential in $m(t\bar{t})$

Covariance for $D$ measurement differential in $m(t\bar{t})$

Results for $\tilde{D}$ measurement differential in $m(t\bar{t})$

Covariance for $\tilde{D}$ measurement differential in $m(t\bar{t})$

Results for $\tilde{D}$ measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Covariance for $\tilde{D}$ measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Results for $D$ measurement differential in $p_{T}(t)$

Covariance for $D$ measurement differential in $p_{T}(t)$

Results for $\tilde{D}$ measurement differential in $p_{T}(t)$

Covariance for $\tilde{D}$ measurement differential in $p_{T}(t)$

The 95% CL upper limits on the branching fraction $\mathcal{B}(\mathrm{H \to SS})$ with $\mathrm{S}\to\tau\tau$ decay, for different LLP masses $m_{\mathrm{S}}$ and proper decay lengths $c\tau_{0}$.

The 95% CL limits on the LLP mass $m_{\mathrm{S}}$ for different proper decay lengths $c\tau_{0}$ assuming a branching fraction of 1% for the $\mathrm{H \to SS}$ decay, and with subsequent $\mathrm{S \to b\overline{b}}$ decay.

The 95% CL limits on the LLP mass $m_{\mathrm{S}}$ for different proper decay lengths $c\tau_{0}$ assuming a branching fraction of 1% for the $\mathrm{H \to SS}$ decay, and with subsequent $\mathrm{S \to d \overline{d}}$ decay.

The 95% CL limits on the dark-sector top quark partner mass $m_{T}$ for different hidden glueball masses $m_{0}$ in the fraternal twin Higgs model.

The 95% CL limits on the dark-sector top quark partner mass $m_{T}$ for different hidden glueball masses $m_{0}$ in the folded SUSY model.

Observed and theoretical cross sections. In the observed, the first uncertainty is statistical, the second is the systematic, and the third the luminosity. In the theoretical, the first uncertainty is due to scale variations, the second due to the choice of PDF. The theoretical cross section has been computed at approximate third-order in QCD, with the addition of third-order corrections of soft gluon emission terms, assuming a top quark mass of 172.5 GeV and using the PDF4LHC21 PDF set.

The number of estimated signal and background events after the fit in the 1j1b, 2j1b, and 2j2b regions compared to the observed number of events. The total uncertainties in the estimated events after the fit are given.

Normalised fiducial differential tW production cross section as a function of the $Leading$ $lepton$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $Jet$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $m(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The p-values from the goodness-of-fit tests comparing the six differential cross section measurements with the predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5 aMC@NLO (aMC) DR, DR2, DS, and DS with a dynamic factor + PYTHIA 8. The complete covariance matrix from the results and the statistical uncertainties in the predictions are taken into account.

Distribution of the RF output for events in the 1j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Response matrix between detector and particle level for the $Leading$ $lepton$ $p_{T}$.

Response matrix between detector and particle level for the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $Jet$ $p_{T}$.

Response matrix between detector and particle level for the $m(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$.

Response matrix between detector and particle level for the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Leading$ $lepton$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $p_{Z}(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Jet$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$ in units of 1.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$ in units of $1/GeV^{2}$.

The model-independent 95\% \CL expected and observed upper limits set on ${\sigma(\PP\to 2\Pa+\PX)\mathcal{B}^2(\Pa\to 2\PGm)\alphaGen}$ over the range $0.21 < \MPa < 9\GeV$ for the combined 2016, 2017, and 2018 analyses. Mass ranges that overlap with \JPsi and \PgU resonances are excluded from the search

The 95\% \CL observed upper limits on the effective coupling $\CahLambda$ of the ALP to the SM Higgs assuming the branching fraction of the ALP to muons is 1 (blue) and 0.1 (orange), for both the expected (dashed) and observed (solid) limits

The 95\% \CL observed upper limits on the effective coupling $\cll/\Lambda$ of the ALP to the SM leptons. The orange shaded region represents the parameter space excluded by this search under three choices of $\CahLambda$ $1\TeV^{-2}$ (solid), $0.1\TeV^{-2}$ (dashed), and $0.01\TeV^{-2}$ (dotted)

The 95\% \CL observed upper limits on $\varepsilon^{2} \mathcal{B}(\ZDark \to \sDark \overline{\sDark}) \mathcal{B}^{2}(\sDark \rightarrow 2\PGm)$ for the vector portal model as a function of the dark scalar mass $\MsDark$ and dark vector boson mass $\MZDark$. The yellow and green bands indicate one and two standard deviation values of the limit, respectively. Because the model-independent limits are calculated only up to a dimuon mass of 60\GeV, the \sDark mass considered for these limits is below 60\GeV

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 90~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 125~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 150~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits (black solid curves) from this search as interpreted in the dark SUSY scenario for the process $\PP \to \Ph \to 2\nOne \to 2\gammaDark + 2\nDark \to 4\PGm + \PX$, with \mbox{$\MnOne = 60\GeV$} and $\MnDark = 1\GeV$. The limits are presented in the plane of $\varepsilon$ and $\MgammaDark$. The color gradient represents different branching fraction assumptions for \mbox{$\mathcal{B}(\Ph \to 2\nOne \to 2\gammaDark + 2\nDark)$}, ranging from dark orange (0.05\%) to light orange (10\%). The degradation of the limit around 1\GeV is attributed to the drop in the dimuon branching fraction $\mathcal{B}(\gammaDark \to 2\mu)$ due to the dimuon resonance of the $\phi$ meson~\cite{Batell:2009yf}

$\chi_{c0}(4500) \to J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

$\chi_{c1}(4685) \ \mathrm{or} \ \chi_{c0}(4700) \to J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

Nonresonant $J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.