Exponent according to the Eq. 5 as a function of transverse momenta extracted from p+Au/p+Al and $^{3}$HeAu/p+Au collision systems. The uncertainties from the $N_{coll}$ calculations and from the overall normalization of p+p cancel in these ratios

Nuclear modification factors in p+Al, p+Au, d+Au, and $^{3}$He+Au in five centrality bins and for inelastic collisions at $\sqrt{s}$ = 200 GeV. The right boxes are the uncertainties of the $N_{coll}$ collisions from the Glauber model, while the left box represents the overall normalization uncertainty from p+p collisions

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

$R_{xA}$ for inelastic collisions compared to three different nuclear PDF calculations and their uncertainties. The data points include the statistical and systematical uncertainties. Relevant data has been tabulated in Figure 7.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. the number of collisions per projectile participant in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. Average $R_{xA}$ in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Panel(b) and panel(c) are tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 9 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 11 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 12 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 13 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 14 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for at least 15 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 4 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 5 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 6 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 9 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 15 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 40 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 12 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 60 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 11 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 80 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 100 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 40 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 60 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 80 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 100 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 6 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category, as well as the observed and expected 95% CL model-independent upper limits on product of cross-section, acceptance and efficiency (in fb). The parameters of the model are determined in a fit to a reduced set of bins, corresponding to the model-independent fit discussed in the text.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category, as well as the observed and expected 95% CL model-independent upper limits on product of cross-section, acceptance and efficiency (in fb). The parameters of the model are determined in a fit to a reduced set of bins, corresponding to the model-independent fit discussed in the text.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Expected 95% CL exclusion contour for the stop pair production model with bino LSP.

Observed 95% CL exclusion contour for the stop pair production model with bino LSP.

Expected 95% CL exclusion contour for the stop pair production model with higgsino LSP.

Observed 95% CL exclusion contour for the stop pair production model with higgsino LSP.

Expected 95% CL exclusion contour for the stop pair production model with wino LSP.

Observed 95% CL exclusion contour for the stop pair production model with wino LSP.

Expected 95% CL excluded cross section for the RPV model with electroweakino prodiction with higgsino LSP hypothesis.

Observed 95% CL excluded cross section for the RPV model with electroweakino prodiction with higgsino LSP hypothesis.

Expected 95% CL excluded cross section for the RPV model with electroweakino prodiction with wino LSP hypothesis.

Observed 95% CL excluded cross section for the RPV model with electroweakino prodiction with wino LSP hypothesis.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category with 20 and 40 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category with 60, 80 and 100 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category with 20 and 40 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category with 60, 80 and 100 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Expected 95% CL excluded cross section for the stop pair production model with bino LSP.

Observed 95% CL excluded cross section for the stop pair production model with bino LSP.

Expected 95% CL excluded cross section for the stop pair production model with wino LSP.

Observed 95% CL excluded cross section for the stop pair production model with wino LSP.

Expected 95% CL excluded cross section for the stop pair production model with higgsino LSP.

Observed 95% CL excluded cross section for the stop pair production model with higgsino LSP.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the stop model with higgsino LSP in the $1\ell$ category.

Acceptance and efficiency for the stop model with higgsino LSP in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 4 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 5 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 6 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 7 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 8 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $2\ell^{\mathrm{sc}}$ category, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 4 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 5 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 6 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 7 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 8 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $2\ell^{\mathrm{sc}}$ category, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Cut flow for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model with $m_{\tilde{g}} = 2$ TeV and $m_{\tilde{\chi}^0_1} = 1$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model with with $m_{\tilde{g}} = 1.6$ TeV and $m_{\tilde{t}} = 1$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with $m_{\tilde{g}} = 2.2$ TeV and $m_{\tilde{\chi}^0_1} = 1.05$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the stop model with $m_{\tilde{t}} = 1.175$ TeV and $m_{\tilde{\chi}^0_1} = 0.7$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the electroweakino production model, considering only the production of $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$, with $m(\tilde{\chi}^{\pm}_1,\tilde{\chi}^0_1)= 250$ GeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. In the $2\ell^{\mathrm{sc}}$ category no events are expected, as only one lepton is expected to be produced in the decay.Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the electroweakino production model, considering only the production of $\tilde{\chi}^0_1 \tilde{\chi}^0_2$, with $m(\tilde{\chi}^0_1,\tilde{\chi}^0_2)= 250$ GeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

The best-fit predicted and observed distributions of the combined NN output in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ with leptonically decaying $\tau$ at 95% confidence level.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ at 95% confidence level, combination of hadronically and leptonically decaying $\tau$ channels.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(e)$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(\mu)$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

Inclusive differential cross section for two electrons and two muons (Max = 1710~GeV).

Mass of the first dilepton pair in the $Z$ mass region (Max = 88.1~GeV).

Mass of the first dilepton pair in the $H$ mass region (Max = 107~GeV).

Mass of the first dilepton pair in the offshell region (Max = 146~GeV).

Mass of the first dilepton pair in the $ZZ$ mass region (Max = 232~GeV).

Mass of the second dilepton pair in the $Z$ mass region (Max = 44.2~GeV).

Mass of the second dilepton pair in the $H$ mass region (Max = 59.6~GeV).

Mass of the second dilepton pair in the offhell region (Max = 88.3~GeV).

Mass of the second dilepton pair in the $ZZ$ mass region (Max = 1000~GeV).

$p_T$ of the first dilepton pair in the $Z$ mass region (Max = 291~GeV).

$p_T$ of the first dilepton pair in the $H$ mass region (Max = 245~GeV).

$p_T$ of the first dilepton pair in the offshell region (Max = 216~GeV).

$p_T$ of the first dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$p_T$ of the second dilepton pair in the $Z$ mass region (Max = 117~GeV).

$p_T$ of the second dilepton pair in the $H$ mass region (Max = 124~GeV).

$p_T$ of the second dilepton pair in the offhell region (Max = 192~GeV).

$p_T$ of the second dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$\cos\theta^*$ between first dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between first dilepton pair in the offshell region.

$\cos\theta^*$ between first dilepton pair in the $H$ mass region.

$\cos\theta^*$ between first dilepton pair in the $ZZ$ mass region.

$\cos\theta^*$ between second dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between second dilepton pair in the $H$ mass region.

$\cos\theta^*$ between second dilepton pair in the offshell mass region.

$\cos\theta^*$ between second dilepton pair in the $ZZ$ mass region.

$\Delta y$ between the dilepton pairs in the $Z$ mass region.

$\Delta y$ between the dilepton pairs in the $H$ mass region.

$\Delta y$ between the dilepton pairs in the offshell region.

$\Delta y$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the dilepton pairs in the $Z$ mass region.

$\Delta \phi$ between the dilepton pairs in the $H$ mass region.

$\Delta \phi$ between the dilepton pairs in the offshell region.

$\Delta \phi$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $Z$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $H$ mass region.

$\Delta \phi$ between the leding lepton pairs in the offshell region.

$\Delta \phi$ between the leading lepton pairs in the $ZZ$ mass region.

Four lepton mass in first $p_T^{4l}$ slice (Max = 726~GeV)

Four lepton mass in second $p_T^{4l}$ slice (Max = 1079~GeV)

Four lepton mass in third $p_T^{4l}$ slice (Max = 1153~GeV)

Four lepton mass in fourth $p_T^{4l}$ slice (Max = 1320~GeV)

Four lepton mass in fifth $p_T^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in first $y^{4l}$ slice (Max = 1710~GeV)

Four lepton mass in second $y^{4l}$ slice (Max = 1116~GeV)

Four lepton mass in third $y^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in fourth $y^{4l}$ slice (Max = 1044~GeV)

Four lepton mass in fifth $y^{4l}$ slice (Max = 918~GeV)

Covariance matrix for M4l

Covariance matrix for M4lvChannel

Covariance matrix for m12_vs_M4l

Covariance matrix for m34_vs_M4l

Covariance matrix for pt12_vs_M4l

Covariance matrix for pt34_vs_M4l

Covariance matrix for CTS12_vs_M4l

Covariance matrix for CTS34_vs_M4l

Covariance matrix for dYpairs_vs_M4l

Covariance matrix for dPhiPairs_vs_M4l

Covariance matrix for dPhiLeadLep_vs_M4l

Covariance matrix for M4lvPt4lbin

Covariance matrix for M4lvRapiditybin

EFT results (linear)

EFT results (full)

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\chi})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\tau})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}$(alignment) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_{h_1}^{125}$(CPV) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Event yields of the SM background processes and the 800 GeV $H^{+}$ sample in the four analysis regions before the fit to the data. Uncertainties include both statistical and systematic uncertainties. The yields of the $H^{+}$ signal sample correspond to a cross-section times branching fraction of 10 pb.

Event acceptance for the different $H^+$ mass signal samples.

$D_s^{\pm}$ invariant yield as a function of $p_{T}$ in 40-80% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 0-10% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 10-40% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.0 < $p_T$ < 2.0 GeV/c, 2.0 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 40-80% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 10-20% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 20-40% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

The integrated $D_{s}/D^{0}$ yield ratio within 1.5 < $p_{T}$ < 5.0 GeV/c as a function of collision centrality. The centrality bins are 60-80%, 40-60%, 20-40%, 10-20% and 0-10%.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

The best-fit expected and observed distributions of the combined NN output in the VRSS for the $e\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $e\tau$ channel for events with 1-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $e\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $\mu\tau$ channel for events with 1-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $\mu\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ at 95% confidence level.

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for diphoton $|cos(\theta*)|$ are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line).

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for average photon transverse momentum are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line) with bands denoting the theoretical uncertainty.

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for average photon transverse momentum are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line).

Measured differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for diphoton rapidity are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line) with bands denoting the theoretical uncertainty.

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for diphoton rapidity are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line).

The 95% CL upper limit on the ALP cross section $\sigma_{\gamma\gamma\rightarrow a \rightarrow\gamma\gamma}$ for the $\gamma\gamma\rightarrow a \rightarrow \gamma\gamma$ process as a function of ALP mass m$_{a}$. The observed upper limit is shown as a solid black line and the expected upper limit is shown by the dashed black line with its $\pm1$ and $\pm2$ standard deviation bands. The discontinuity at $m_{a}=70 GeV$ is caused by the increase of the mass-bin width which brings an increase in signal acceptance.

The 95% CL upper limit on the ALP coupling $1/\Lambda_{a}$ for the $\gamma\gamma\rightarrow a \rightarrow \gamma\gamma$ process as a function of ALP mass m$_{a}$. The observed upper limit is shown as a solid black line and the expected upper limit is shown by the dashed black line with its $\pm1$ and $\pm2$ standard deviation bands. The discontinuity at $m_{a} = 70 GeV$ is caused by the increase of the mass-bin width which brings an increase in signal acceptance.

Fiducial cross section for light-by-light scattering

The $p_{T}$ dependance of $v_{4}^{Inclusive}$, $v_{4}^{Linear}$, $v_{4}^{Mode-coupled}$, $\chi_{4,22}$ and $\rho_{4,22}$

The $p_{T}$ dependance of $v_{5}^{Inclusive}$, $v_{5}^{Linear}$, $v_{5}^{Mode-coupled}$, $\chi_{5,23}$ and $\rho_{5,23}$