$v_2\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 20-30% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 30-40% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 40-50% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 50-60% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 0-5% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 5-10% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 10-20% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 20-30% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 30-40% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 40-50% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 50-60% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 0-5% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 5-10% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 10-20% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 20-30% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 30-40% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 40-50% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 50-60% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 0-5% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 5-10% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 10-20% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 20-30% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 30-40% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 40-50% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of ${\rm K}^{0}_{\rm{S}}$ as a function of $p_{\rm T}$ for the 50-60% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 0-5% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 5-10% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 10-20% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 20-30% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 30-40% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 40-50% centrality interval.

$v_2\{2, |\Delta\eta| > 2.0\}$ of $\Lambda$+$\overline{\Lambda}$ as a function of $p_{\rm T}$ for the 50-60% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 0-10% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 10-30% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of $\pi^{\pm}$ as a function of $p_{\rm T}$ for the 30-50% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 0-10% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 10-30% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of $\rm{K}^{\pm}$ as a function of $p_{\rm T}$ for the 30-50% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 0-10% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 10-30% centrality interval.

$v_3\{2, |\Delta\eta| > 2.0\}$ of p+$\rm\overline{p}$ as a function of $p_{\rm T}$ for the 30-50% centrality interval.

The $p_{\rm T}$ of $v_2^{\mathrm{max}}$ for $\pi^{\pm}$ as a function of centrality.

The $p_{\rm T}$ of $v_2^{\mathrm{max}}$ for p+$\rm\overline{p}$ as a function of centrality.

Pb-remnant side ZN signal normalized to MB value vs. ZN centrality percentile in p-Pb collisions at 8.16 TeV

p-remnant side ZN signal normalized to MB value vs. ZN centrality percentile in p-Pb collisions at 5.02 TeV

p-remnant side ZN signal normalized to MB value vs. ZN centrality percentile in p-Pb collisions at 8.16 TeV

Pb-remnant side ZN signal normalized to MB value vs. average Ncoll in p-Pb collisions at 5.02 TeV

Pb-remnant side ZN signal normalized to MB value vs. average Ncoll in p-Pb collisions at 8.16 TeV

p-remnant side ZN signal normalized to MB value vs. average Ncoll in p-Pb collisions at 5.02 TeV

p-remnant side ZN signal normalized to MB value vs. average Ncoll in p-Pb collisions at 8.16 TeV

ZN signal normalized to MB value vs. normalized charged-particle multiplicity in p-p collisions at 13 TeV

p-remnant side ZN signal normalized to MB value vs. normalized charged-particle multiplicity in p-Pb collisions at 8.16 TeV

ZP signal normalized to MB value vs. normalized charged-particle multiplicity in p-p collisions at 13 TeV

p-remnant side ZP signal normalized to MB value vs. average charged-particle multiplicity in p-Pb collisions at 8.16 TeV

ZN signal normalized to MB value vs. leading transverse momentum in |eta|<0.8 in p-p collisions at 13 TeV

ZP signal normalized to MB value vs. leading transverse momentum in |eta|<0.8 in p-p collisions at 13 TeV

Differential cross section as a function of $m_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $p_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $a_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\phi_{\eta}*$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\pi-\Delta\phi_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $|cos\theta*|^{(CS)}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

$v_3$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_3$ vs $p_T$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_3$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_2$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ vs $p_T$, d+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ ratio vs $p_T$, p/d/3He+Au at 200 GeV, 0-5% central

$v_3$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_3$ vs $p_T$, d+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_3$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$c_2$ vs $\eta$, p+Au at 200 GeV, 0-5% central

$c_2$ vs $\eta$, d+Au at 200 GeV, 0-5% central

$c_2$ vs $\eta$, 3He+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, p+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, d+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, 3He+Au at 200 GeV, 0-5% central

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$c_n$ vs detector combination, p+p at 200 GeV, Minimum Bias

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The geometric signal acceptance as a function of $c\tau$. The acceptance is shown for four different LLP mass hypotheses: 7, 15, 40, and 55 GeV. The acceptance is defined by requiring at least 1 LLP to decay in the region defined as 400 cm $< |z| <$ 1100 cm, $r < $ 695.5 cm, and $|\eta| <$ 2.4.

Distribution of $m_{T}(W)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $M_{3}$ in the $N_{jet}\geq 3$ signal region.

Distribution of $M_{3}$ in the $N_{jet}\geq 3$ signal region.

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

Summary of the one-dimensional intervals at 68 and 95% CL.

Summary of the one-dimensional intervals at 68 and 95% CL.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.

Distribution of the W boson transverse mass in the WZ 3l signal region. Events from conversions, and DY processes are grouped into the 'Others' category. The vertical error bars represent the statistical uncertainty in the data and the shaded band the uncertainty in the prediction. The signal contributions are scaled to the measured cross sections (postfit).

Results obtained in this analysis and other diboson production cross section measurements at different center-of-mass energies for the CMS, ATLAS, CDF, and D0 Collaborations are presented, and compared with the NNLO QCD x NLO EW and NLO predictions from MATRIX. The vertical error bars represent the uncertainty in the measured cross section.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 30-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 40-50% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 50-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 60-70% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 70-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 20-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 40-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 60-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 0-10% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 10-20% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 20-30% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 30-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 40-50% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 50-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 60-70% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 70-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 20-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 40-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 60-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Two-particle correlations $c_{1}\{2\}$ versus $Q^2$ with a rapidity separation and high-$p_{\textrm{T}}$ constraint: $\Delta \eta > 2$, $p_{\textrm{T}}$ > 0.5 GeV. Photoproduction data are shown at $Q^2$ = 0 GeV$^2$, while NC DIS is for $Q^2$ > 5 GeV$^2$.

Two-particle correlations $c_{2}\{2\}$ versus $Q^2$. Photoproduction data are shown at $Q^2$ = 0 GeV$^2$, while NC DIS is for $Q^2$ > 5 GeV$^2$.

Two-particle correlations $c_{2}\{2\}$ versus $Q^2$ with a rapidity separation: $\Delta \eta > 2$. Photoproduction data are shown at $Q^2$ = 0 GeV$^2$, while NC DIS is for $Q^2$ > 5 GeV$^2$.

Two-particle correlations $c_{2}\{2\}$ versus $Q^2$ with a high-$p_{\textrm{T}}$ constraint: $p_{\textrm{T}}$ > 0.5 GeV. Photoproduction data are shown at $Q^2$ = 0 GeV$^2$, while NC DIS is for $Q^2$ > 5 GeV$^2$.

Two-particle correlations $c_{2}\{2\}$ versus $Q^2$ with a rapidity separation and high-$p_{\textrm{T}}$ constraint: $\Delta \eta > 2$, $p_{\textrm{T}}$ > 0.5 GeV. Photoproduction data are shown at $Q^2$ = 0 GeV$^2$, while NC DIS is for $Q^2$ > 5 GeV$^2$.

Normalised charged-particle multiplicity distribution $dN/dN_{\textrm{ch}}$ in photoproduction.

Normalised charged-particle transverse momentum distribution $dN/dp_{\textrm{T}}$ in photoproduction.

Normalised charged-particle pseudorapidity distribution $dN/d\eta$ in photoproduction.

Two-particle correlations $c_1\{2\}$ versus $|\Delta \eta|$ in photoproduction.

Two-particle correlations $c_2\{2\}$ versus $|\Delta \eta|$ in photoproduction.

Two-particle correlations $c_1\{2\}$ versus $\left< p_{\textrm{T}} \right>$ in photoproduction.

Two-particle correlations $c_2\{2\}$ versus $\left< p_{\textrm{T}} \right>$ in photoproduction.

Four-particle cumulant correlations $c_1\{4\}$ versus $p_{T,1}$ in photoproduction.

Four-particle cumulant correlations $c_2\{4\}$ versus $p_{T,1}$ in photoproduction.