Distributions of $p_{T}^{Z_1}$. Small background contributions are denoted as "other backgrounds", including 4$\mu$ events containing non-prompt muons estimated from data and from $ttV$, $VVV$, and Higgs boson production processes.

Distributions of $p_{T}^{Z_2}$. Small background contributions are denoted as "other backgrounds", including 4$\mu$ events containing non-prompt muons estimated from data and from $ttV$, $VVV$, and Higgs boson production processes.

Distributions of the mass difference of the $Z_1$ and $Z_2$ candidates. Small background contributions are denoted as "other backgrounds", including 4$\mu$ events containing non-prompt muons estimated from data and from $ttV$, $VVV$, and Higgs boson production processes.

The pDNN output discriminant variable distributions for low mass with a signal sample at 35 GeV and 51 GeV, respectively.

The pDNN output discriminant variable distributions for high mass with a signal sample at 35 GeV and 51 GeV, respectively.

Mass spectra of $m_{Z2}$ for the pDNN-selected events with a signal sample at 15 GeV.

Mass spectra of $m_{Z1}$ for the pDNN-selected events with a signal sample at 51 GeV.

The $p_0$-value scan across the Z' mass signal regions.

95% CL upper limits (expected and observed) on the cross-sections times branching fraction. The discontinuity at 42 GeV represents the border of the low/high mass classifiers.

95% CL upper limits (expected and observed) on the coupling parameter. The discontinuity at 42 GeV represents the border of the low/high mass classifiers.

The parameterized mass resolution $\sigma_{m_{\mu\mu}}$ as a function of $m_{Z'}$ using fully simulated $Z' \rightarrow \mu^+\mu^-$ events.

The transverse momentum $p_{T}$ distributions for each muon (ordered with $p_{T}$ from (a) to (d)). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

The transverse momentum $p_{T}$ distributions for each muon (ordered with $p_{T}$ from (a) to (d)). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

The transverse momentum $p_{T}$ distributions for each muon (ordered with $p_{T}$ from (a) to (d)). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

The transverse momentum $p_{T}$ distributions for each muon (ordered with $p_{T}$ from (a) to (d)). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. The plots (a) to (d) are the $\eta$ distributions of the 4 muons ($p_{T}$ ordered). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. The plots (a) to (d) are the $\eta$ distributions of the 4 muons ($p_{T}$ ordered). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. The plots (a) to (d) are the $\eta$ distributions of the 4 muons ($p_{T}$ ordered). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. The plots (a) to (d) are the $\eta$ distributions of the 4 muons ($p_{T}$ ordered). In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. Plots (a) to (d) are the angular separations, $\Delta R$ of the two muons that formed $Z_1$ and $Z_2$, and the mass distributions of $Z_1$ and $Z_2$. In addition to the major backgrounds from the SM $Z(Z^*)\rightarrow 4\mu$ production, other background, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. Plots (a) to (d) are the angular separations, $\Delta R$ of the two muons that formed $Z_1$ and $Z_2$, and the mass distributions of $Z_1$ and $Z_2$. In addition to the major backgrounds from the SM $Z(Z^*)\rightarrow 4\mu$ production, other background, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. Plots (a) to (d) are the angular separations, $\Delta R$ of the two muons that formed $Z_1$ and $Z_2$, and the mass distributions of $Z_1$ and $Z_2$. In addition to the major backgrounds from the SM $Z(Z^*)\rightarrow 4\mu$ production, other background, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Kinematic distributions of the pre-selected $4\mu$ events. Plots (a) to (d) are the angular separations, $\Delta R$ of the two muons that formed $Z_1$ and $Z_2$, and the mass distributions of $Z_1$ and $Z_2$. In addition to the major backgrounds from the SM $Z(Z^*)\rightarrow 4\mu$ production, other background, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Data and MC comparison of the $Z \to 4\mu$ invariant mass using pre-selected 4$\mu$ events. Figure (b) shows the distribution between 80 GeV~and 110 GeV~in Figure (a). Figure (c) shows the distribution of transverse momentum of $4\mu$ system. In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Data and MC comparison of the $Z \to 4\mu$ invariant mass using pre-selected 4$\mu$ events. Figure (b) shows the distribution between 80 GeV~and 110 GeV~in Figure (a). Figure (c) shows the distribution of transverse momentum of $4\mu$ system. In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

Data and MC comparison of the $Z \to 4\mu$ invariant mass using pre-selected 4$\mu$ events. Figure (b) shows the distribution between 80 GeV~and 110 GeV~in Figure (a). Figure (c) shows the distribution of transverse momentum of $4\mu$ system. In addition to the major background from the SM $Z(Z^*)\rightarrow 4\mu$ production, other backgrounds, including 4$\mu$ events containing non-prompt muons estimated from data, and from $ttV$, $VVV$, and Higgs boson production processes, are included in the plots. Examples of the Z' signal from $pp\rightarrow Z'\mu^+\mu^- \rightarrow 4\mu$ process with masses of 15 and 51 GeV are also shown in the plots.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $p_{T}^{miss}$ after the preselection from 2017 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $p_{T}^{miss}$ after the preselection from 2018 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $N_{Jet}$ after the preselection from 2017 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $N_{Jet}$ after the preselection from 2018 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Simulated distribution of $p_{T}(l)$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are given for a top squark mass of 300 GeV and $\Delta {m}=10$, 30,50, and 80 GeV.

Simulated distribution of $p_{T}(l)$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are shown for four different $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ values, all corresponding to the same $\Delta {m} = 30$ GeV.

Simulated distribution of $p_{T}^{miss}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are given for a top squark mass of 300 GeV and $\Delta {m}=10$, 30,50, and 80 GeV.

Simulated distribution of $p_{T}^{miss}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are shown for four different $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ values, all corresponding to the same $\Delta {m} = 30$ GeV.

Simulated distribution of $N_{Jet}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are given for a top squark mass of 300 GeV and $\Delta {m}=10$, 30,50, and 80 GeV.

Simulated distribution of $N_{Jet}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are shown for four different $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ values, all corresponding to the same $\Delta {m} = 30$ GeV.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 10$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 20$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 30$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 40$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 50$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 60$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 70$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 80$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 10$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 20$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 30$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 40$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 50$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 60$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 70$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 80$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}(l)$ for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}(l)$ for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}^{miss}$ for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}^{miss}$ for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Data yield and prediction of the total background and its composition in the eight SRs as defined in for the 2017 analysis. The predictions of W + jets, ttbar, and nonprompt lepton processes are based on data, while those of rare backgrounds are based on simulation. The predictions and their associated uncertainties are pre-fit. The uncertainties are the quadratic sum of the statistical and systematic uncertainties for background processes predicted from data, and the uncertainties on the cross section for the backgrounds predicted from simulation. The yields for two signal points, with $\Delta {m} = 10$ GeV and $\Delta {m} = 80$ GeV, are also shown. The bins corresponding to different $\Delta {m}$ trainings are not statistically independent.

The ratio of the number of observed events over the predicted total background per SR.

Data yield and prediction of the total background and its composition in the eight SRs as defined in for the 2018 analysis. The predictions of W + jets, ttbar, and nonprompt lepton processes are based on data, while those of rare backgrounds are based on simulation. The predictions and their associated uncertainties are pre-fit. The uncertainties are the quadratic sum of the statistical and systematic uncertainties for background processes predicted from data, and the uncertainties on the cross section for the backgrounds predicted from simulation. The yields for two signal points, with $\Delta {m} = 10$ GeV and $\Delta {m} = 80$ GeV, are also shown. The bins corresponding to different $\Delta {m}$ trainings are not statistically independent.

The ratio of the number of observed events over the predicted total background per SR.

Exclusion limit at 95% CL for the four-body decay of the top squark as a function of $m(\mathrm{\widetilde{t}}_{1})$ and $\Delta {m}$, for the combined data of the years 2016, 2017, and 2018. The lines represent the observed limits. These limits are derived using the expected top squark pair production cross section. The green line represents the central values, and the other lines the variations due to the theoretical or experimental uncertainties.

Exclusion limit at 95% CL for the four-body decay of the top squark as a function of $m(\mathrm{\widetilde{t}}_{1})$ and $\Delta {m}$, for the combined data of the years 2016, 2017, and 2018. The lines represent the expected limits. These limits are derived using the expected top squark pair production cross section. The green line represents the central values, and the other lines the variations due to the theoretical or experimental uncertainties.

Exclusion limit at 95% CL for the four-body decay of the top squark as a function of $m(\mathrm{\widetilde{t}}_{1})$ and $\Delta {m}$, for the combined data of the years 2016, 2017, and 2018.

Reconstruction level differential cross section spectla, obtained for the central acceptance, $|\eta| < 3$, for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology compared to the to the EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The number of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The transeverse momentum of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The total energy of all Particle Flow candidatesin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

Unfolded diffraction enhanced differential cross section for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffraction enhanced differential cross section for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Reconstruction level diffraction enhanced differential cross section spectrum for Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap, compared with the spectrum obtained with events satisfying the ZDC veto requirement $E_{ZDC−} < 1 \mathrm{~TeV}$ which selects only the events without lead nuclear break up (without correction for HF undetectable energy). Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $10-20\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $20-30\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $30-40\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $40-50\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $50-60\%$ centrality.

$\Lambda K^{0}_{S} \oplus \overline{\Lambda} K^{0}_{S}$ correlation meassurement in $0-80\%$ centrality.

$\Lambda\Lambda\oplus\overline{\Lambda}\overline{\Lambda}$ correlation meassurement in $0-80\%$ centrality.

The extracted $R_{inv}$ from $K^{0}_{S} K^{0}_{S}$ correlations as a function of centrality.

The extracted $\lambda$ parameters from $K^{0}_{S} K^{0}_{S}$ correlations as a function of centrality.

Extracted effective range ($d_{0}$) and real components of the complex scattering length ($\Re f_{0}$) from the $\Lambda\Lambda\oplus\overline{\Lambda}\overline{\Lambda}$ and $\Lambda K^{0}_{S} \oplus \overline{\Lambda} K^{0}_{S}$ correlations.

Extracted imaginary components of the complex scattering length ($\Im f_{0}$) and real components of the complex scattering length ($\Re f_{0}$) from the $\Lambda K^{0}_{S} \oplus \overline{\Lambda} K^{0}_{S}$ correlation.

Event yields of the signal and SM background processes in the six analysis regions after the fit to the data under the $t \rightarrow uX$ hypothesis assuming $m_X = 30$ GeV. Total includes signal and background.The quoted uncertainties take into account correlations and constraints of the nuisance parameters and include both the statistical and systematic uncertainties. Negative correlations between the $t\bar{t} +$ light, $t\bar{t} + \geq1b$ and $t\bar{t} + \geq1c$ modelling uncertainties can make the uncertainty in the total yields smaller than in the individual components.

Event yields of the signal and SM background processes in the six analysis regions after the fit to the data under the $t \rightarrow cX$ hypothesis assuming $m_X = 30$ GeV. Total includes signal and background.The quoted uncertainties take into account correlations and constraints of the nuisance parameters and include both the statistical and systematic uncertainties. Negative correlations between the $t\bar{t} +$ light, $t\bar{t} + \geq1b$ and $t\bar{t} + \geq1c$ modelling uncertainties can make the uncertainty in the total yields smaller than in the individual components.

Event acceptance times efficiency in percent for every $t\rightarrow uX$ and $t\rightarrow cX$ mass signal sample.

Cut flow for the scalar signal in the $t\rightarrow uX$ decay combining both quark and anti-quark samples. Shown for each signal are the corresponding mass, the number of generated events, the number of reconstructed ("Reco") events, the events that pass the lepton triggers, the events that have only one electron ("el") or only one muon ("mu") with p$_{\text{T}}$ larger than 27 GeV, the number of events with at least four jets with p$_{\text{T}}$ larger than 25 GeV, the number of events with at least three $b$-tagged jets at the 70% efficiency working point, and the number of events with at least two $b$-tagged jets at the 60% efficiency working point and at least another one at the 70%. The quoted yields do not include reweighting.

Cut flow for the scalar signal in the $t\rightarrow cX$ decay combining both quark and anti-quark samples. Shown for each signal are the corresponding mass, the number of generated events, the number of reconstructed ("Reco") events, the events that pass the lepton triggers, the events that have only one electron ("el") or only one muon ("mu") with p$_{\text{T}}$ larger than 27 GeV, the number of events with at least four jets with p$_{\text{T}}$ larger than 25 GeV, the number of events with at least three $b$-tagged jets at the 70% efficiency working point, and the number of events with at least two $b$-tagged jets at the 60% efficiency working point and at least another one at the 70%. The quoted yields do not include reweighting.

Electron $p_{T}$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_{T}^{miss}$ distribution in the electron channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Muon $p_T$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_T^{miss}$ distribution in the muon channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the E+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the MU+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Theoretical prediction for the production and decay cross-section of SSM WPRIME bosons (W --> LEPTON + NU) as a function of the WPRIME mass, for the range 200-6400 GeV and inclusive in phase-space (no cuts). The uncertainty shown covers PDF and ALPHAS uncertainties. These theoretical predictions are obtained with PYTHIA (LO) and corrected to NNLO in QCD using FEWZ program.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW as a function of the mass of the WPRIME boson, for the electron channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime, being equal to gW, the SM coupling.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW, as a function of the mass of the WPRIME boson, for the muon channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime , being equal to gW, the SM coupling.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the electron channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the muon channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the combined electron and muon channels for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{231}$ coupling in the RPV SUSY model with a STAU mediator, as a function of the mass of STAU, for the electron channel and for the production coupling, ${\lambda}^{'}_{3ij}$=0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying the ones shown by the ratio ${\lambda}^{'}_{3ij}$/0.5

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{132}$ coupling in the RPV SUSY model with a stau mediator, as a function of the mass of the stau, for the muon channel and for the production coupling, ${\lambda}^{'}_{3ij}$ = 0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying these ones by the ratio ${\lambda}^{'}_{3ij}$/0.5

Region in the oblique (W,Y) parameter phase space allowed by the current analysis at 95% CL. The combination of the electron and muon channel distributions from the 2017 and 2018 data sets at sqrt(s) = 13 TeV is used, having 101 fb-1 of luminosity.

Negative Log-Likelihood for the determination of the oblique W parameter, for the combination of the electron and muon channel and for 2017 and 2018 data sets (L=101 fb-1) at sqrt(s) = 13 TeV.

Total background estimate and observed data in all bins used in the limit calculation for the SSM WPRIME limits. The total event sample is split into the three years of data taking and the two lepton flavors.The distributions being shown are $M_T$ ones. The associated integrated luminosities for each year are of 35.9, 41.5, and 59.4 fb$̣^{-1}. The mass bins are numbered in ascending order.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW +QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

The $m_{W\gamma}$ distribution for muon events (electron events are also used and the muon events are shown here as an example.) satisfying the aQGC region selection and used to set constraints on the anomalous gauge coupling parameters (left). The gray line represents a nonzero $f_{M,2}/Λ^4$ setting. Events with $m_{W\gamma} > 1500 GeV$ are included in the last bin.

Exclusion limits at the 95% CL for each aQGC coefficient, derived from the mW$\gamma$ distribution, assuming all other coefficients are set to zero. Unitarity bounds corresponding to each operator are also listed. All coupling parameter limits are in TeV^{−4}, while Ubound values are in TeV.

Results summary.

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity (0$-$10%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 1.5 GeV/$\it{c}$

Inclusive charged-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in 20$-$60% semi-central Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity (20$-$60%) Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 2 GeV/$\it{c}$

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity (20$-$60%) Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 3 GeV/$\it{c}$

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in semi-central (20$-$60%) Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 5 GeV/$\it{c}$

The profile likelihood scan as a function of Bs to mu+mu- and B0 to mu+mu- decay branching fractions in 2D.

Correlation matrix of all the parameters including nuisance parameters used in the profile likelihood unfolding fit on $A_{C}$.