The impact of the eight most important systematic uncertainties on \(R_t\), in %.

The post-fit yields in the two SRs. All uncertainties applied in the analysis are included.

Limits on EFT parameters and Vtb extracted from the analysis. The upper and lower limits correspond to 95% CL.

Results of the main analysis.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-plus CR.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-minus CR.

Pre-fit distribution of the invariant mass of the untagged jet and the b-tagged jet, \(m(jb)\), in SR plus.

Pre-fit distribution of the absolute value of the pseudorapidity of the untagged jet, \(|\eta(j)|\), in SR plus.

Pre-fit distribution of the absolute value of the difference in \(p_{\text{T}}\) between the reconstructed W boson and the jet pair, \(|\Delta p_{\text{T}}(W, jb)|\), in SR plus.

Pre-fit distribution of the difference in azimuth angle between the reconstructed W boson and the jet pair, \(|\Delta\phi(W, jb)|\), in SR plus.

Pre-fit distribution of the invariant mass of the reconstructed top quark, \(m(t)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the charged lepton and the untagged jet, \(|\Delta\eta(\ell ,j)|\), SR plus.

Pre-fit distribution of the angular distance of the charged lepton and the untagged jet, \(\Delta R(\ell ,j)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the b-tagged jet and the charged lepton, \(|\Delta\eta(b,\ell )|\), in SR plus.

Suppl. Fig. 7. The reconstructed shower angle for all single photon events. The individual signal and background event type categories added together form the unconstrained prediction.

Suppl. Fig. 7. The constrained covariance matrix for the reconstructed shower angle for all single-photon events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 7. The unconstrained covariance matrix for reconstructed shower angle for all single-photon events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 3. The reconstructed number of protons, with a proton KE energy threshold of 35 MeV, for all single-photon events. The individual signal and background event type categories added together form the unconstrained prediction.

Fig. 4, Suppl. Fig. 8. The reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The individual signal and background event type categories added together form the unconstrained prediction.

Fig. 4. The constrained covariance matrix for the reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 4. The unconstrained covariance matrix for reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 6. The reconstructed shower energy for zero proton events binned from 150 to 1250 MeV, with LEE prediction. The individual signal and background event type categories added together form the unconstrained prediction. The "LEE" category is the unfolded MiniBooNE LEE model prediction scaled under the neutrino-target nucleon interaction assumption.

Fig. 6. The constrained covariance matrix for the reconstructed shower energy for zero proton events binned from 150 to 1250 MeV. Both with and without the LEE model added. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 6. The unconstrained covariance matrix for reconstructed shower energy for zero proton events binned from 150 to 1250 MeV. Both with and without the LEE model added. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 10. The reconstructed shower backwards projected distance for zero proton events. The individual signal and background event type categories added together form the unconstrained prediction.

Suppl. Fig. 9. The reconstructed shower angle for zero proton events. Includes LEE model prediction. The individual signal and background event type categories added together form the unconstrained prediction. The "LEE" category is the unfolded MiniBooNE LEE model prediction scaled under the neutrino-target nucleon interaction assumption.

Suppl. Fig. 9. The constrained covariance matrix for the reconstructed shower angle for zero proton events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 9. The unconstrained covariance matrix for reconstructed shower angle for zero proton events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

1$\gamma$Xp shower energy and angle selection efficiencies for true 1$\gamma$Xp events. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. Energy and Angle values represent bin lower edges.

1$\gamma$0p shower energy and angle selection efficiencies for true single-photon events with no other true particles of any energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. Energy and Angle values represent bin lower edges.

Covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Pearson data statistical uncertainties have been included, and include small correlations due to events which can be selected by both WC and Pandora. The first four bins are the WC $1\gamma Np$, WC $1\gamma 0p$, Pandora $1\gamma 1p$, and Pandora $1\gamma 0p$ channels, and the next four sets of 16 bins are the NC $\pi^0 Np$, NC $\pi^0 0p$, $\nu_\mu$CC $Np$, and $\nu_\mu$CC $0p$ channels. This corresponds to Fig. 1 of the paper and Fig. 6 of the Supplemental Material.

Constrained WC $1\gamma Np$ energy distribution. There are 6 bins from 0 to 600 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained WC $1\gamma 0p$ energy distribution. There are 14 bins from 0 to 700 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained Pandora $1\gamma 1p$ energy distribution. There are 6 bins from 0 to 600 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained Pandora $1\gamma 0p$ energy distribution. There are 14 bins from 0 to 700 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained covariance matrix showing constrained uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Pearson data statistical uncertainties have been included, and include small correlations due to events which can be selected by both WC and Pandora. The first seven bins are the WC $1\gamma Np$ reconstructed shower energy from 0-600 MeV with an overflow bin, the next 15 bins are the WC $1\gamma 0p$ reconstructed shower energy from 0-700 MeV with an overflow bin, the next seven bins are the Pandora $1\gamma 1p$ reconstructed shower energy from 0-600 MeV with an overflow bin, and the last 15 bins are the Pandora $1\gamma 0p$ reconstructed shower energy from 0-700 MeV with an overflow bin. This corresponds to Fig. 2 of the paper.

2D efficiency matrix showing the efficiency of the WC $1\gamma Np$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 15 bins from 0 to 1500 MeV in true photon energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the WC $1\gamma 0p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 1p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 15 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 0p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the WC $1\gamma 0p$ selection with respect to all true zero-primary-final-state-proton NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 0p$ selection with respect to all true zero-primary-final-state-proton NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied.

Total invariant mass distribution $m_{4\ell}$ in Signal Region 2.

Selected events in the ($m_{4\ell}$, $\left\langle m_{\ell\ell}\right\rangle$) plane, SR1

Selected events in the ($m_{4\ell}$, $\left\langle m_{\ell\ell}\right\rangle$) plane, SR2

Expected background events in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR1

Expected signal shape for $(m_{S}, m_{Z_d}) = (110~\text{GeV},30~\text{GeV}$) in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR1. The signal distributions are normalized to unit volume.

Expected background events in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR2. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to $\langle m_{\ell\ell}\rangle$ values between 83 and 99 GeV, are deleted.

Expected signal shape for $(m_{S}, m_{Z_{d}}) = (750~\text{GeV},150~\text{GeV}$) in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR2. The signal distributions are normalized to unit volume. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to $\langle m_{\ell\ell}\rangle$ values between 83 and 99 GeV, are deleted.

Expected 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$, SR1

Observed 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$, SR1

Observed / Expected 95% CL limit ratio in SR1

Expected 95% CL limit on $\mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$ with Z-veto, SR2

Observed 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$ with Z-veto, SR2

Observed / Expected 95% CL limit ratio with Z-veto, SR2

Local $p_0$-values in the ($m_{S}, m_{Z_{d}}$) plane, SR1

Local $p_0$-values in the ($m_{S}, m_{Z_{d}}$) plane, SR2 with Z-veto applied

Observed data in Signal Region 1 in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane.

Observed data in Signal Region 2 in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to <m_ll> values between 83 and 99 GeV, are deleted.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $4e$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $2e2\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $4\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $4e$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $2e2\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $4\mu$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $4e$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $2e2\mu$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $4\mu$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $4e$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $2e2\mu$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $4\mu$ final state.

Isolated-photon cross section ratio of pp collisions at $\sqrt{s}=13~\mathrm{TeV}$ over $\sqrt{s}=5.02~\mathrm{TeV}$ in data (left column) and NLO pQCD calculation from JETPHOX (right column) for $R=0.4$. Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided. The data normalisation uncertainty is provided in the paper.

Isolated-photon cross section ratio of pp collisions at $\sqrt{s}=7~\mathrm{TeV}$ over $\sqrt{s}=5.02~\mathrm{TeV}$ in data (left column) and NLO pQCD calculation from JETPHOX (right column) for $R=0.4$. Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided. The data normalisation uncertainty is provided in the paper.

Nuclear modification factor ($R_{\rm AA}$) of isolated photons for Pb$-$Pb and pp collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm AA}$ obtained with a cone size value of $R=0.2$. Each column for each Pb$-$Pb collisions centrality class. The last column for the NLO pQCD JETPHOX calculations, Pb$-$Pb (nPDF) over pp (PDF). Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided.

Nuclear modification factor ($R_{\rm AA}$) of isolated photons for Pb$-$Pb and pp collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm AA}$ obtained with a cone size value of $R=0.4$. Each column for each Pb$-$Pb collisions centrality class. The last column for the NLO pQCD JETPHOX calculations, Pb$-$Pb (nPDF) over pp (PDF). Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided.

Ratios $R_{\rm CSP}$ and $R_{\rm CSC}$ of isolated photons cross sections in data for Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm CSP}$, ratio central over semi-peripheral collisions 50$-$70%. $R_{\rm CSC}$, ratio central over semi-central collisions 30$-$50%. Ratios obtained with a cone size value of $R=0.2$. Data statistical and systematic uncertainties are provided.

Ratios $R_{\rm CSP}$ and $R_{\rm CSC}$ of isolated photons cross sections in data for Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm CSP}$, ratio central over semi-peripheral collisions 50$-$70%. $R_{\rm CSC}$, ratio central over semi-central collisions 30$-$50%. Ratios obtained with a cone size value of $R=0.4$. Data statistical and systematic uncertainties are provided.

Self-normalised ZN energy as a function of the self-normalised charged-particle-multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised ZN energy as a function of the self-normalised charged-particle-multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$\Lambda+\overline{\Lambda}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Lambda+\overline{\Lambda}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$K^{0}_{S}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$\Lambda+\overline{\Lambda}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Lambda+\overline{\Lambda}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the standalone classification (V0M event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the standalone classification (V0M event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the low-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The multiplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The multiplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $1.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 2.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, trig} < 3.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 3.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $1.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, trig} < 3.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $1.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, trig} < 3.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

The isolated-photon cross section ratio, 13 TeV to 7 TeV, from pQCD NLO calculations with JETPHOX. The calculations were obtained by choosing factorisation, normalisation, and fragmentation scales. The parton distribution function used in the calculations is NNPDF4.0 for the $\sqrt{s}=$13 TeV measurement and CT14 for the $\sqrt{s}=$7 TeV measurement. The fragmentation function is BFG II for both measurements.

Differential cross section of isolated photons as a function of $x_\mathrm{T}^{\gamma}$ measured in pp collisions at $\sqrt{s}=$13 TeV.

The $p_\mathrm{T}$-integrated nuclear modification factor $R_\mathrm{pPb}$ of charm production at midrapidity in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV.