Distribution of the NN output in the $\geq$1fj$\,$SR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the NN output in the $t\bar{t}\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Total event yield in the $W\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the scalar sum of the jet transverse momenta in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark+photon system in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the scalar sum of the jet transverse momenta in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the invariant mass of the $b$-tagged jet and the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of $p_{\mathrm{T}}$ of the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the difference in $\eta$ between the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the energy of the system formed by the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark and photon system in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the parton-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the particle-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

This table lists the kinematic requirements on parton-level objects used to define of the fiducial phase space for the parton-level measurement. Frixione isolation ($\href{https://arxiv.org/abs/hep-ph/9801442}{\text{hep-ph/9801442}}$) with a chosen radius of $\Delta R = 0.2$ is applied to photons ($\gamma$). The measured fiducial parton-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b) = 688\pm 23(\text{stat.})^{+75}_{-71}(\text{syst.})\,$fb.

This table lists the kinematic requirements on particle-level objects used to define of the fiducial phase space for the particle-level measurement. The particle level objects are photons ($\gamma$) not from a hadron decay, neutrinos not from a hadron decay ($\nu$), prompt electrons and muons ($\ell$) "dressed" by adding close-by ($\Delta R < 0.1$) photons, and anti-$k_t$ $R = 0.4$ jets built from stable particles ($\tau > 30\,$ps) and tau leptons excluding neutrinos and prompt dressed muons. Jets are $b$-tagged ($b$-jet) using ghost-matched $b$-hadrons with $p_{\text{T}} > 5\,$GeV. Apart from the kinematic requirements, isolation and overlap removal criteria are applied. Jets within $\Delta R = 0.4$ of a photon are removed if the $p_{\text{T}}$ of charged particles within $\Delta R = 0.3$ of the photon is smaller than $10\,\%$ of its $p_{\text{T}}$. Jets within $\Delta R = 0.4$ of a lepton are removed. Events are removed where a photon is close ($\Delta R < 0.4$) to a lepton or a surviving jet. The measured fiducial particle-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b)+\sigma_{t(\rightarrow l\nu b\gamma)q} = 303\pm 9(\text{stat.})^{+33}_{-32}(\text{syst.})\,$fb.

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the FCNC top-quark candidate in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the SM top-quark candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the transverse momentum of the Z boson candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the third lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

$p_\mathrm{T}$-differential production yield of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 0-10% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 0-10% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 30-50% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 30-50% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{T}$-differential elliptic flow coefficient $v_2$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$.

$p_\mathrm{T}$-differential elliptic flow coefficient $v_2$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 5-10%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 5-10%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 5-10%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-20%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-20%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-20%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 20-30%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 20-30%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 20-30%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-40%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-40%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-40%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 40-50%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 40-50%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 40-50%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-60%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-60%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-60%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 60-70%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 60-70%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 60-70%.

$p_{T}$-distributions of pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

$p_{T}$-distributions of protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

$p_{T}$-distributions of phi ($\phi$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 0-10%.

$p_{T}$-distributions of phi ($\phi$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-30%.

$p_{T}$-distributions of phi ($\phi$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-50%.

$p_{T}$-distributions of phi ($\phi$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-70%.

$p_{T}$-distributions of phi ($\phi$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 0-5%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 0-5%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 5-10%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 5-10%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-20%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-20%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 20-30%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 20-30%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-40%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-40%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 40-50%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 40-50%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-60%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-60%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 60-70%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 60-70%.

$p_{T}$-distributions of kaons ($K^{+}+K^{-}$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

$p_{T}$-distributions of protons ($p+pbar$)/pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

$p_{T}$-distributions of phi ($\phi$)/protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 0-10%.

$p_{T}$-distributions of phi ($\phi$)/protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 10-30%.

$p_{T}$-distributions of phi ($\phi$)/protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 30-50%.

$p_{T}$-distributions of phi ($\phi$)/protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 50-70%.

$p_{T}$-distributions of phi ($\phi$)/protons ($p+pbar$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV. Centrality class 70-90%.

Integrated yield ratios of kaons ($K^{+}+K^{-}$) to pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

Integrated yield ratios of protons ($p+pbar$) to pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

Integrated yield ratios of phi ($\phi$) to pions ($\pi^{+}+\pi^{-}$) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

Example of $\Lambda$ signal for 0-40% most central events, over mixed-event background for the bin $-0.05 < y_{cm} < 0.05$ and reduced transverse masses between $100-150 MeV/c^{2}$.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{0}_{S}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{0}_{S}$ for the 0-40% most central events. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 2.1997626 x 10^{9}$

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\Lambda$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\Lambda$ for the 0-40% most central events. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 2.1997626 x 10^{9}$

Rapidity distribution of $K^{0}_{S}$

Rapidity distribution of $K^{0}_{S}$

Rapidity distribution of $\Lambda$

Rapidity distribution of $\Lambda$

Compilation of mid-rapidity yields for central Au+Au collisions as a function of $\sqrt{s_{NN}}$ of $K^{0}_{S}$ and $\Lambda$.

Compilation of mid-rapidity yields for central Au+Au collisions as a function of $\sqrt{s_{NN}}$ of $K^{0}_{S}$ and $\Lambda$.

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $K^{0}_{S}$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $K^{0}_{S}$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $\Lambda$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $\Lambda$

Comparison of the shape of the rapidity distribution of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the rapidity distribution of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the rapidity distribution of $\Lambda$ to various transport model versions.

Comparison of the shape of the rapidity distribution of $\Lambda$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $\Lambda$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $\Lambda$ to various transport model versions.

Differential yield of $K^{0}_{S}$ as function of rapdidity and reduced transverse mass for the four centrality classes.

Differential yield of $K^{0}_{S}$ as function of rapdidity and reduced transverse mass for the four centrality classes. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 5.64247975 x 10^{8} (0 - 10\%)$, $5.85743394 x 10^{8} (10 - 20\%)$, $5.76369451 x 10^{8} (20 - 30\%)$ or $4.73401752 x 10^{8} (30 - 40\%)$

Differential yield of $\Lambda$ as function of rapdidity and reduced transverse mass for the four centrality classes.

Differential yield of $\Lambda$ as function of rapdidity and reduced transverse mass for the four centrality classes. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 5.64247975 x 10^{8} (0 - 10\%)$, $5.85743394 x 10^{8} (10 - 20\%)$, $5.76369451 x 10^{8} (20 - 30\%)$ or $4.73401752 x 10^{8} (30 - 40\%)$

$K^{0}_{S}$ and $\Lambda$ multiplicities in full phase space and inverse slopes at mid-rapidity $T_{Eff}$ for a given centrality.

$K^{0}_{S}$ and $\Lambda$ multiplicities in full phase space and inverse slopes at mid-rapidity $T_{Eff}$ for a given centrality.

Value of the parameter $\alpha$ extracted for $K^{+}$, $K^{-}$, $K^{0}_{S}$, $\Lambda$ and $\Phi$, as displayed in Figure 5 - Global $Mult/\langle A_{Part}\rangle$ Fit

Value of the parameter $\alpha$ extracted for $K^{+}$, $K^{-}$, $K^{0}_{S}$, $\Lambda$ and $\Phi$, as displayed in Figure 5 - Global $Mult/\langle A_{Part}\rangle$ Fit

The K0S cross section for 5 jet ET ranges.

The LAMBDA cross section for 5 jet ET ranges.

The parameter Alpha as a function of x_F. Statistical and systematic uncertainties are indicated separately.

SPHERICITY distribution.

PLANARITY distribution.

OBLATENESS distribution.

HEAVY JET Mass distribution in the standard scheme definition.

HEAVY JET Mass distribution in the E scheme definition.

LIGHT JET Mass distribution in the standard scheme definition.

JET MASS DIFFERENCE distribution in the standard scheme definition.

TOTAL JET MASS distribution in the standard scheme definition.

TOTAL JET MASS distribution in the E scheme definition.

Wide Jet Broadening (BMAX) distribution.

Narrow Jet Broadening (BMIN) distribution.

Total Jet Broadening (BTOT) distribution.

Jet Broadening Difference (BDIFF) distribution.

C-PARAMETER distribution.

Integrated 2 Jet rates for YCUT=0.04 (Durham, angle ordered Durham and Cambridge alogrithms) and YCUT = 0.8 (JADE alogrithm).

Integrated 3 Jet rates for YCUT=0.04 (Durham, angle ordered Durham and Cambridge alogrithms) and YCUT = 0.8 (JADE alogrithm).

Integrated Jet Cone Energy Fractions (JCEF) for different anglular intervals.

Integrated values of the Energy Energy Correlations (EEC) for different angular intervals.

Integrated values of the Energy Energy Correlations Asymmetry (AEEC) for different angular intervals.

Mean values of Y23.

Mean values of Thrust related observables.

Mean values of the Jet Broadenings.

Mean values of jet masses in the standard definition scheme.

Mean values of jet masses in the E definition scheme.

Mean values of jet masses in the p definition scheme.

Mean values of the C- and D- Parameter variables.

Rapidity distributions (relative to the thrust axis) for c.m. energies from183 to 196 GeV.

Rapidity distributions (relative to the thrust axis) for c.m. energies from200 to 207 GeV.

LN(1/XP) distributions for c.m. energies from 183 to 196 GeV.

LN(1/XP) distributions for c.m. energies from 200 to 207 GeV.

PTIN distributions for c.m. energies from 183 to 196 GeV.

PTIN distributions for c.m. energies from 200 to 207 GeV.

PTOUT distributions for c.m. energies from 183 to 196 GeV.

PTOUT distributions for c.m. energies from 200 to 207 GeV.

1-THRUST distributions for c.m. energies from 183 to 196 GeV.

1-THRUST distributions for c.m. energies from 200 to 207 GeV.

THRUST-MAJOR distributions for c.m. energies from 183 to 196 GeV.

THRUST-MAJOR distributions for c.m. energies from 200 to 207 GeV.

THRUST-MINOR) distributions for c.m. energies from 183 to 196 GeV.

THRUST-MINOR distributions for c.m. energies from 200 to 207 GeV.

Oblateness distributions for c.m. energies from 183 to 196 GeV.

Oblateness distributions for c.m. energies from 200 to 207 GeV.

Wide Jet Broadening (BMAX) distributions for c.m. energies from 183 to 196 GeV.

Wide Jet Broadening (BMAX) distributions for c.m. energies from 200 to 207 GeV.

Total Jet Broadening (BTOT) distributions for c.m. energies from 183 to 196GeV.

Total Jet Broadening (BTOT) distributions for c.m. energies from 200 to 207GeV.

Jet Broadening Difference (BDIFF) distributions for c.m. energies from 183 to 196 GeV.

Jet Broadening Difference (BDIFF) distributions for c.m. energies from 200 to 207 GeV.

C-PARAMETER distributions for c.m. energies from 183 to 196 GeV.

C-PARAMETER distributions for c.m. energies from 200 to 207 GeV.

D-PARAMETER distributions for c.m. energies from 183 to 196 GeV.

D-PARAMETER distributions for c.m. energies from 200 to 207 GeV.

Heavy Jet Mass distributions in the standard scheme definition for c.m. energies from 183 to 196 GeV.

Heavy Jet Mass distributions in the standard scheme definition for c.m. energies from 200 to 207 GeV.

Heavy Jet Mass distributions in the p scheme definition for c.m. energies from 183 to 196 GeV.

Heavy Jet Mass distributions in the p scheme definition for c.m. energies from 200 to 207 GeV.

Heavy Jet Mass distributions in the E scheme definition for c.m. energies from 183 to 196 GeV.

Heavy Jet Mass distributions in the E scheme definition for c.m. energies from 200 to 207 GeV.

Light Jet Mass distributions in the standard definition for c.m. energies from 183 to 196 GeV.

Light Jet Mass distributions in the standard definition for c.m. energies from 200 to 207 GeV.

Jet Mass Difference distributions from c.m.energies from 183 to 196 GeV.

Jet Mass Difference distributions from c.m.energies from 200 to 207 GeV.

Sphericity distributions for c.m.energies from 183 to 196 GeV.

Sphericity distributions for c.m.energies from 200 to 207 GeV.

Planarity distributions for c.m.energies from 183 to 196 GeV.

Planarity distributions for c.m.energies from 200 to 207 GeV.

Aplanarity distributions for c.m.energies from 183 to 196 GeV.

Aplanarity distributions for c.m.energies from 200 to 207 GeV.