Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e u}|$, versus $S_{e u}$ mass for scalar LQs coupled to electrons.

Observed and Expected UL exclusions on the $BR(H\to S)$ of generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $S_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Observed and Expected UL exclusions on the number of 'signal' yields in sideband A from the model agnostic fit, assuming no 'SUEP-like' signal contamination in the CRs.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e d}|$, versus $S_{e d}$ mass for scalar LQs coupled to electrons.

Postfit values for SR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for PR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu u}|$, versus $S_{\mu u}$ mass for scalar LQs coupled to muons.

Postfit values for CRDY in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for CRTT in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu d}|$, versus $S_{\mu d}$ mass for scalar LQs coupled to muons.

Per-selection efficiencies on Generic - $m_{A'} = 1.0\;GeV$ $BR(A' \rightarrow \pi\pi) = 1.0$ Hadronic - $m_{A'} = 0.7\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.15$ $BR(A' \rightarrow \pi\pi) = 0.7$ Leptonic - $m_{A'} = 0.5\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.2$ and $BR(A' \rightarrow \pi\pi) = 0.6$ samples. Selections: 2 OSSF leptons, one SUEP cluster, $60\leq Z_{m}\leq 120$ GeV, $p_{T}^Z\geq 25$ GeV, veto events with b-tagged jets, $p_{T}^{SUEP} > 60\;GeV$, and $dR^{Ak4}_{Ak15} < 1.5$. Efficiencies are the ratio of the histogram size at the selection step to the previous step (except twoLepton vs. initial). The final column is the total cumulative efficiency (dR vs. initial). Note that these selection numbers were all calculated using samples generated with 600,000 events for each signal point considered.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Interpolated observed and expected UL exclusions on $BR(H\to S) = 1$ for generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$. The discreate heatmap points can be found in the UL Generic (Observed & Expected) figure in this sandbox

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e u}|$, versus $V_{e u}$ mass for vector LQs coupled to electrons.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e d}|$, versus $V_{e d}$ mass for vector LQs coupled to electrons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu u}|$, versus $V_{\mu u}$ mass for vector LQs coupled to muons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu d}|$, versus $V_{\mu d}$ mass for vector LQs coupled to muons.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the SR (t-pass) and SR (t-fail) after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. A representative mono-top signal (vector coupling scenario) with a mediator mass of 1 TeV, a DM candidate mass of 150 GeV, and a cross section of 1 pb is overlaid as an orange line. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the SR (t-pass) and SR (t-fail) after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. A representative mono-top signal (vector coupling scenario) with a mediator mass of 1 TeV, a DM candidate mass of 150 GeV, and a cross section of 1 pb is overlaid as an orange line. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W(e$\nu$) (t-pass) and W(e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W(e$\nu$) (t-pass) and W(e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W($\mu\nu$) (t-pass) and W($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W($\mu\nu$) (t-pass) and W($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z(ee) (t-pass) and Z(ee) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z(ee) (t-pass) and Z(ee) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z($\mu\mu$) (t-pass) and Z($\mu\mu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z($\mu\mu$) (t-pass) and Z($\mu\mu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the t$\bar{t}$ (e$\nu$) (t-pass) and t$\bar{t}$ (e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the t$\bar{t}$ (e$\nu$) (t-pass) and t$\bar{t}$ (e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the tt ($\mu\nu$) (t-pass) and tt ($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the tt ($\mu\nu$) (t-pass) and tt ($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the $\gamma$ (t-pass) and $\gamma$ (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the $\gamma$ (t-pass) and $\gamma$ (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Upper limits at 95% CL on $\sigma\mathcal{B}$ of mono-top production presented in the two- dimensional plane spanned by the mediator and DM candidate masses for a mediator mass between 200 and 2250 GeV and a DM candidate mass between 1 and 1125 GeV only consider- ing on-shell decays of the mediator to the DM candidates. The mediator has vector couplings to quarks and DM candidates in the upper plot and axial-vector couplings in the lower plot. The median expected exclusion range is indicated by a black solid line, demonstrating the search sensitivity of the analysis. The 68% probability interval of the expected exclusion is shown in black dashed lines. Contours of theory predictions for constant values of $\sigma\mathcal{B}$ are shown in gray dashed lines. The observed exclusion contour of mediator and DM candidate masses is rep- resented by the red solid line. The exclusion contour obtained from measurements of the DM relic density $\Omega_{\mathrm{nbm}}h^2$ by the Planck Collaboration is shown in the gray solid line.

Upper limits at 95% CL on $\sigma\mathcal{B}$ of mono-top production presented in the two- dimensional plane spanned by the mediator and DM candidate masses for a mediator mass between 200 and 2250 GeV and a DM candidate mass between 1 and 1125 GeV only consider- ing on-shell decays of the mediator to the DM candidates. The mediator has vector couplings to quarks and DM candidates in the upper plot and axial-vector couplings in the lower plot. The median expected exclusion range is indicated by a black solid line, demonstrating the search sensitivity of the analysis. The 68% probability interval of the expected exclusion is shown in black dashed lines. Contours of theory predictions for constant values of $\sigma\mathcal{B}$ are shown in gray dashed lines. The observed exclusion contour of mediator and DM candidate masses is rep- resented by the red solid line. The exclusion contour obtained from measurements of the DM relic density $\Omega_{\mathrm{nbm}}h^2$ by the Planck Collaboration is shown in the gray solid line.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}} > 2$ GeV and $p_{\mathrm{T}}^{\mathrm{ch}} > 1$ GeV.

$\kappa^{2}_\mathrm{K}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{2}_\mathrm{p}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{2}_\mathrm{p}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\pi,\mathrm{K}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\pi,\mathrm{K}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\pi,\mathrm{p}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\pi,\mathrm{p}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\mathrm{p},\mathrm{K}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\mathrm{p},\mathrm{K}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{2}_{\mathrm{Q}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{2}_{\mathrm{Q}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\mathrm{Q},\mathrm{K}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\mathrm{Q},\mathrm{K}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\mathrm{Q},\mathrm{p}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_{\mathrm{Q},\mathrm{p}}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$C_\mathrm{Q,p}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$C_\mathrm{Q,p}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$C_\mathrm{Q,K}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$C_\mathrm{Q,K}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$C_\mathrm{p,K}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$C_\mathrm{p,K}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_\mathrm{Q,p}/\kappa^{2}_\mathrm{Q}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_\mathrm{Q,K}/\kappa^{2}_\mathrm{Q}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$(2\kappa^{11}_\mathrm{Q,K}-\kappa^{11}_\mathrm{p,K})/\kappa^{2}_\mathrm{K}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$(2\kappa^{11}_\mathrm{Q,p}-\kappa^{11}_\mathrm{p,K})/\kappa^{2}_\mathrm{p}$ as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_\mathrm{Q,p}/\kappa^{2}_\mathrm{Q}$ normalized by its value at 0$-$5% centrality as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 1 $p_\mathrm{T}$ acceptance.

$\kappa^{11}_\mathrm{Q,p}/\kappa^{2}_\mathrm{Q}$ normalized by its value at 0$-$5% centrality as a function of centrality (%) in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV for Set 2 $p_\mathrm{T}$ acceptance.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the signal region (SR). The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the signal region (SR). The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 300 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 500 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W+ signal region.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W- signal region.

Post-fit distributions of the transverse mass in the W+ signal region.

Post-fit distributions of the transverse mass in the W- signal region.

Post-fit distributions of the dimuon mass in the Z boson signal region.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 12 < mT < 18 GeV. The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 12 < mT < 18 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (larger value) and quadratic (smaller value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

The table compares theoretical and measured values of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured values of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured ratios of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

The table compares theoretical and measured ratios of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

Comparison of the measured and theoretical cross sections for Z, W+, W-, and W production at different center-of-mass energies. The theoretical predictions are generated with DYTurbo in combination with the NNPDF 3.1 PDF set.

The nuclear suppression factor $S^{\mathrm{J}/\psi}$ of incoherent $\mathrm{J}/\psi$ meson photoproduction as a function of Bjorken $x$ in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux.

The ratio between incoherent and coherent $\mathrm{J}/\psi$ meson photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The total covariance matrix of the total incoherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total incoherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total incoherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.