Differential cross section times dimuon branching fraction of Y(1S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Nuclear modification factor of Y(1S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

RFB of Y(1S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(2S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(3S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(1S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(2S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(3S) versus $E^{|\eta_{lab}|>4}_{T}$.

Nuclear modification factor of Y(1S), Y(2S), and Y(3S) integrated over pT and $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

The W' boson 95% CL limits on the product of cross section and branching fraction. The expected and observed limits are shown for the high VLQ mass range.

The W' boson 95% CL limits on the product of cross section and branching fraction. The expected and observed limits are shown for the low VLQ mass range.

Expected 95% CL limits for generalized hypotheses varying the fraction of tB (F(VLQ=B)) and bT (F(VLQ=T)) from the W′ decay

Expected 95% CL limits for generalized hypotheses varying the VLQ branching fraction to qH and qZ

Observed 95% CL limits for generalized hypotheses varying the fraction of tB (F(VLQ=B)) and bT (F(VLQ=T)) from the W′ decay

Observed 95% CL limits for generalized hypotheses varying the VLQ branching fraction to qH and qZ

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 40 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 45 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 50 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 55 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 60 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 65 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 70 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 75 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 80 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 85 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 90 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 95 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

The observed 95% CL upper limits on the product of the cross section, and branching fraction for the production of X$\to\phi\phi\to$4b, with $m_\phi$ = 100 GeV. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to various extended higgs sector model cross sections as a function of the parameter $m_XN/f$.

Data event yields and predicted signal and background yields (including subdominant $t\bar{t}$) in bins of dijet and average jet mass.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Efficiency corrected $\phi$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.1 in the method section.

Efficiency and acceptance corrected $K^{*0}$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.4 in the method section.

$\phi$-meson $\rho_{00}$ obtained from 1st- and 2nd-order event planes. The red stars (gray squares) show the $\phi$-meson $\rho_{00}$ as a function of beam energy, obtained with the 2nd-order (1st-order) EP.

$\phi$-meson $\rho_{00}$ with respect to different quantization axes. $\phi$-meson $\rho_{00}$ as a function of beam energy, for the out-of-plane direction (stars) and the in-plane direction (diamonds). Curves are fits based on theoretical calculations with a $\phi$-meson field. The corresponding $G_s^{(y)}$ values obtained from the fits are shown in the legend.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

The $M_{JJ}$ distribution for the number of observed events (black markers) compared with the estimated backgrounds (filled histograms) and their uncertainties (hatched areas) in the SR1. The distributions expected from the signal under three $M_X$ and $M_Y$ hypotheses and assuming a cross section of 1 fb are also shown. The lower panels show the ''Pulls'' defined as (observed events - expected events)/$\sqrt{\smash[b]{\sigma_{obs}^{2} - \sigma_{exp}^{2}}}$, where $\sigma_{obs}$ and $\sigma_{exp}$ are the statistical and total uncertainties in the observation and the background estimation, respectively. The minus sign accounts for the correlation between data and the data-driven estimation.

The 95% confidence level expected and observed upper limits on $\sigma$(pp->X->YH->4b) for different values of $M_X$ and $M_Y$. The areas within the red and black contours represent the regions where the cross sections predicted by NMSSM and TRSM, respectively, are larger than the experimental limits. The areas within the dashed and dotted contours show the excluded masses at -1 standard deviation from the expected limits.

The soft-drop mass distribution of the top quark candidate jets in the 2018 jets+lepton category, in the tight ParticleNet region, after the joint fit in all-jets and jets+lepton categories. Observed data (black markers) and the postfit estimate (filled histograms) are shown for the three jet categories. The lower panel shows the ''Pulls'' defined as (observed events - expected events)/$\sqrt{\smash[b]{\sigma_{obs}^{2} + \sigma_{exp}^{2}}}$, where $\sigma_{obs}$ and $\sigma_{exp}$ are the statistical and total uncertainties in the observation and the background estimation, respectively.

The 95% confidence level expected and observed upper limits on $\sigma$(pp->X->YH->4b) for different values of $M_X$ and $M_Y$. The areas within the red and black contours represent the regions where the cross sections predicted by NMSSM and TRSM, respectively, are larger than the experimental limits. The areas within the dashed and dotted contours on the left show the excluded masses at -1 standard deviation from the expected limits.

Measurements of the signal strengths $\mu_{\text{s}}$ in the STXS stage-1.2 bins for the ggH, qqH, and VH processes, for the combination of the CB- and VH-analyses. Central values and a split of uncertainties are also provided with each result.

Correlation matrix of the POIs of the measured STXS stage-1.2 signal strengths for the combination of the NN- an VH-analyses.

Correlation matrix of the POIs of the measured STXS stage-1.2 signal strengths for the combination of the CB- an VH-analyses.

Cross section measurements in the STXS stage-0 bins for the combination of the NN- and VH-analyses. Central values and combined statistical and systematic uncertainties are given for each measurement.

Cross section measurements in the STXS stage-1.2 bins for the combination of the NN- and VH-analyses. Central values and combined statistical and systematic uncertainties are given for each measurement.

Cross section measurements in the STXS stage-0 bins for the combination of the CB- and VH-analyses. Central values and combined statistical and systematic uncertainties are given for each measurement.

Cross section measurements in the STXS stage-1.2 bins for the combination of the CB- and VH-analyses. Central values and combined statistical and systematic uncertainties are given for each measurement.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 4 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Observed signal strength ($\mu$) for a SM Higgs boson at 125.38 GeV. The labels untagged combined,dijet combined and combined represent the results obtained from simultaneous fits of the untagged categories, dijet categories, and full set of categories, respectively. The black solid line shows $\mu=1$, and the red dashed line shows the best fit value $\hat{\mu}=2.4$ of all categories combined.

Upper limit (95% CL) on the signal strength ($\mu$) relative to the SM prediction, as a function of the assumed value of the Higgs boson mass used in the fit.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=1600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2000\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2200\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2400\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2800\;\mathrm{GeV}$.

$V_{n\Delta}$ coefficients for minimum bias events as a function of N$_{\text{trk}}$ for $ 1.0 < p_\mathrm{T} < 3.0 GeV$ in pPb collisions at 8.16 TeV.

Single-particle azimuthal anisotropy $v_2$ for $\gamma$p events as a function of N$_{\text{trk}}$ for $0.3 < p_\mathrm{T} < 3.0 GeV$ in pPb collisions at 8.16 TeV.

Single-particle azimuthal anisotropy $v_2$ for minimum bias events as a function of N$_{\text{trk}}$ for $0.3 < p_\mathrm{T} < 3.0 GeV$ in pPb collisions at 8.16 TeV.

Single-particle azimuthal anisotropy $v_2$ for $\gamma$p events as a function of N$_{\text{trk}}$ for $1.0 < p_\mathrm{T} < 3.0 GeV$ in pPb collisions at 8.16 TeV.

Single-particle azimuthal anisotropy $v_2$ for minimum bias events as a function of N$_{\text{trk}}$ for $1.0 < p_\mathrm{T} < 3.0 GeV$ in pPb collisions at 8.16 TeV.