$P_z$ sine coefficients of $\Lambda+\bar{\Lambda}$ as a function of transverse momentum in isobar collisions at 200 GeV.

$P_z$ sine coefficients of $\Lambda+\bar{\Lambda}$ as a function of the number of participants in isobar collisions at 200 GeV. The number of participants for isobar collisions were taken from Table III of Phys.Rev.C105, 014901 (2022) and averaged. Same data as the centrality dependence in Fig. 3.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) invariant yields in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. The vertical bars and the boxes represent statistical and systematic uncertainties, respectively. The horizontal bars indicate the bin width.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $p_{\rm T}$ in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with STAR (yellow stars) and PHENIX (green squares) published results, and Duke ((modified Langevin transport model, blue line) and PHSD (parton-hadron-string dynamics model, orange line) model calculations. Vertical bars and boxes around data points represent combined statistical and systematic uncertainties from both Au+Au and $p$+$p$ measurements, respectively. Boxes at unity show the global uncertainties, which for this analysis include the 8% global uncertainty on $p$+$p$ reference and the $N_{\rm coll}$ uncertainties. The left box is for PHENIX and the right one for STAR.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $N_{\rm part}$ in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with PHENIX measurements (green squares), and Duke (blue line) and PHSD (orange line) model calculations. Vertical bars and boxes around data points represent statistical and systematic uncertainties from Au+Au measurements, respectively. The gray band represents the $N_{\rm coll}$ uncertainties. The boxes at unity show the global uncertainties including the total uncertainties of the $p$+$p$ reference.

Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for single(+z)-single(−z) proton reconstruction categories.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for different proton reconstruction categories in the $Z \rightarrow \mu\mu $ analysis, and for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the positive CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the negative CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The di-proton invariant mass is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The $m_{miss}$ is shown.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow ee$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow \mu\mu$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $\gamma$ analysis.

The distribution of the BDT discriminants for the 2016--2018 data set for three different categories in the combined single-electron and single-muon channels. The three categories are defined by the number of resolved t tags ($N_\textrm{RT}$), b tags ($N_\textrm{b}$), and jets ($N_\textrm{j}$), selected as representative based on their sensitivity to signal. Here, $\textrm{t}\bar{\textrm{t}} + \geq 1 \textrm{b}$ refers to $\textrm{t}\bar{\textrm{t}}$ events with at least one additional b jet, while $\textrm{t}\bar{\textrm{t}} + 0 \textrm{b}$ includes all other $\textrm{t}\bar{\textrm{t}}$ events not produced in association with a boson. The TOP grouping contains single top quark production along with the other $\textrm{t}\bar{\textrm{t}}$ processes not explicitly shown, and EW refers to events that contain W and Z bosons but no top quarks.

Measured signal strength ($\mu = \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}}/ \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}} ^{\textrm{SM}}$), corresponding cross section (in fb), and the expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from all analysis channels. Uncertainties shown in this table are statistical and systematic uncertainties respectively. In this table, SSDL\&ML and OSDL (2016) represent results from EPJC 80 (2020) 75 and JHEP 11 (2019) 082 respectively.

The distribution of the BDT discriminants for the 2016--2018 data set for three different categories in the combined single-electron and single-muon channels. The three categories are defined by the number of resolved t tags ($N_\textrm{RT}$), b tags ($N_\textrm{b}$), and jets ($N_\textrm{j}$), selected as representative based on their sensitivity to signal. Here, $\textrm{t}\bar{\textrm{t}} + \geq 1 \textrm{b}$ refers to $\textrm{t}\bar{\textrm{t}}$ events with at least one additional b jet, while $\textrm{t}\bar{\textrm{t}} + 0 \textrm{b}$ includes all other $\textrm{t}\bar{\textrm{t}}$ events not produced in association with a boson. The TOP grouping contains single top quark production along with the other $\textrm{t}\bar{\textrm{t}}$ processes not explicitly shown, and EW refers to events that contain W and Z bosons but no top quarks. The backgrounds and $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal (derived from the fit) are shown as a stacked histogram. The hatched bands correspond to the estimated total uncertainty after the fit. While the bins are shown to be equal width, they do not correspond to equal width in BDT value.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, HT > 1400 GeV. The background from QCD multijet and \ttbar production is derived from control regions in the data. Estimates for the \tttt signal and other backgrounds are shown using simulated samples.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The sample VR category shown is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, $H_T > 1400$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{BT} \geq 2$, HT > 1100 GeV. The background from QCD multijet and \ttbar production is derived from control regions in the data. Estimates for the \tttt signal and other backgrounds are shown using simulated samples.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The sample VR category shown is defined by $N_\textrm{RT} \geq 2$, $H_T > 1100$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty.

Expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from each final state and the combination with previous CMS results (EPJC 80 (2020) 75 and JHEP 11 (2019) 082). The same-sign dilepton and multilepton (SSDL\&ML) final state results are from EPJC 80 (2020) 75.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, $H_T > 1400$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty after the fit.

Measured signal strength ($\mu = \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}}/ \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}} ^{\textrm{SM}}$), corresponding cross section (in fb), and the expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from all analysis channels. Uncertainties shown in this table are statistical and systematic uncertainties respectively. In this table, SSDL\&ML and OSDL (2016) represent results from EPJC 80 (2020) 75 and JHEP 11 (2019) 082 respectively.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{RT} = 2$, $H_T > 1100$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty after the fit.

Expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from each final state and the combination with previous CMS results (EPJC 80 (2020) 75 and JHEP 11 (2019) 082). The same-sign dilepton and multilepton (SSDL\&ML) final state results are from EPJC 80 (2020) 75.

Measured signal strength ($\mu = \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}}/ \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}} ^{\textrm{SM}}$), corresponding cross section (in fb), and the expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from all analysis channels. This table shows production from each analysis channel in this Letter, the combination of those channels, the results from previously published results, and the full combination of all CMS 2016-2018 results.

Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the combined VBF, ttH, VH and ggH channels using all available CMS data for 2016--2018, Run1, 2015, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.

Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the VBF, ttH(lep), ttH(had), V(qq)H, Z(ll)H, Z(bb)H and ggH channels using all available CMS data, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.

The 90% CL upper limits on the spin-independent DM-nucleon scattering cross section in Higgs-portal models, assuming a scalar, fermionic or vectorial DM candidate (with a UV-complete vector model), compared to direct detection experimets CRESST-III, DarkSide-50, LUX_ZEPLIN and PandaX-4T (2022 results).

Postfit event yields in each category and hadronic recoil bin of the ttH category SR for the different background and signal processes summed across 2016, 2017 and 2018 samples. The background and signal yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only. The Postfit signal contributions for a Higgs boson, produced in ttH, VH (split into WH and ZH), ggH and VBF, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are reported.

Postfit event yields in each category and hadronic recoil bin of the VH category SR for the different background and signal processes summed across 2016, 2017 and 2018 samples. The background and signal yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only. The Postfit signal contributions for a Higgs boson, produced in ttH, VH (split into WH and ZH), ggH and VBF, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are reported.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category e$+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category e$+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category $\mu+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category $\mu+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category ee$+$jets, $\mu\mu+$jets, and $\ell\ell+$jets CRs for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category ee$+$jets, $\mu\mu+$jets, and $\gamma+$jets CRs for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

The observed and expected 95% CL upper limits on the minimal coupling vector LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.

The observed and expected 95% CL upper limits on the Yang-Mills coupling vector LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the scalar LQ.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the vector LQ minimal coupling scenario.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the vector LQ Yang-Mills coupling scenario.

Baseline subtracted $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles for $4 < p_{\rm T}^{\rm e} < 12$ ${\rm GeV}/c$ and $1 < p_{\rm T}^{\rm assoc} < 7$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Baseline subtracted $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles for $4 < p_{\rm T}^{\rm e} < 12$ ${\rm GeV}/c$ and $1 < p_{\rm T}^{\rm assoc} < 7$ ${\rm GeV}/c$ in p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Associated yield of near-side (NS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 12$ ${\rm GeV}/c$ in pp and p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Associated yield of away-side (AS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 12$ ${\rm GeV}/c$ in pp and p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Width of near-side (NS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 12$ ${\rm GeV}/c$ in pp and p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Width of away-side (AS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 12$ ${\rm GeV}/c$ in pp and p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Associated yield of near-side (NS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Associated yield of away-side (AS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Width of near-side (NS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Width of away-side (AS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Associated yield of near-side (NS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Associated yield of away-side (AS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in p--Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Width of near-side (NS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Width of away-side (AS) $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles as a function of $p_{\rm T}^{\rm{assoc}}$ for $4 < p_{\rm T}^{\rm e} < 7$ and $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ in p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Baseline subtracted $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles for $4 < p_{\rm T}^{\rm e} < 7$ ${\rm GeV}/c$ and $1 < p_{\rm T}^{\rm assoc} < 7$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Baseline subtracted $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles for $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ and $1 < p_{\rm T}^{\rm assoc} < 7$ ${\rm GeV}/c$ in pp collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Baseline subtracted $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles for $4 < p_{\rm T}^{\rm e} < 7$ ${\rm GeV}/c$ and $1 < p_{\rm T}^{\rm assoc} < 7$ ${\rm GeV}/c$ in p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Baseline subtracted $\Delta\varphi$ distribution between heavy-flavor decay electrons and associated charged particles for $7 < p_{\rm T}^{\rm e} < 16$ ${\rm GeV}/c$ and $1 < p_{\rm T}^{\rm assoc} < 7$ ${\rm GeV}/c$ in p$\textendash$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at $\sqrt{s} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in Pb-p collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at various LHC energies.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of collision energy in pp collisions.

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - CL1).

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - ZNA).

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Lambda lifetime measured in different centrality intervals. Only statistical uncertainties are shown.