The 95% CL limits on the effective coupling constants. The scale of new physics is set to $\Lambda = 1$ TeV.

The 95% CL limits on the branching fractions. The scale of new physics is set to $\Lambda = 1$ TeV.

The pre-fit distribution of the NN discriminant for the FCNC $tu\gamma$ LH signal scaled to 10 times the observed limit and all SM backgrounds as well as measured data in the SR. The last bin of the distribution contains the overflow and the first bin also contains the underflow. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction.

The pre-fit distribution of the photon transverse momentum for SM backgrounds as well as measured data in the $t\bar{t}\gamma$ CR. The last bin of the distribution contains the overflow. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction.

The pre-fit distribution of the photon transverse momentum for the all SM backgrounds as well as measured data in the $W\gamma$+jets CR. The last bin of the distribution contains the overflow. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $q\bar{q}\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $gg\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the WMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-0 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $q\bar{q} \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The jet mass distribution of large-$R$ jets originating from the hadronic decay of $W$ and $Z$ bosons produced from the decay of BSM bosons with mass $m_X = 2000$ GeV. The decays simulated are for the production models $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ with a spin-1 resonance $X^{\pm}$ and $gg\to X^0 \to Z\gamma$ with a spin-0 resonance $X^{0}$.

The photon $E_{T}$ selection criteria as a function of resonance mass for D2 categories. Selections for four different signal channels are shown together.

The photon $E_{T}$ selection criteria as a function of resonance mass for ZMASS or WMASS categories. Selections for four different signal channels are shown together.

Di-muon mass distribution for the $\mu\mu\mu$ signal region for data and expected background. The expected signal distribution for $m_a = 35$ GeV is shown assuming $\sigma(t\bar{t}a)\times \text{Br}(a\rightarrow\mu\mu) = 4$ fb. Rare background processes include $ZZ+$jets, $WWZ$, $WZZ$, $ZZZ$, and $t\bar{t}t\bar{t}$.

Fit of the Double Crystal Ball width as a function of the $a$-boson mass in the $e\mu\mu$ channel.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $t\bar{t}a$.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 120$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 140$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 160$ GeV.

Comparison of the expected background with data in a region dominated by non-prompt background. The region is defined exactly like the $e\mu\mu$ signal region but requiring the two selected muons to have the same electrical charge.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $e\mu\mu$ channel.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $\mu\mu\mu$ channel.

Acceptance times efficiency for the $t\bar{t}a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ hypothesis.

Acceptance times efficiency for the $t \rightarrow bH^+$, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $e\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Acceptance times efficiency for the $t \rightarrow bH^+ $, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $\mu\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Expected lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Observed lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Cutflow for two signal mass points for $t\bar{t}a$, $m_{a} = 20$ GeV and $m_{a} = 60$ GeV, and one signal mass point for $H^+ \rightarrow Wa$ for the electron channel $e\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the electron channel $e\mu\mu$. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 20$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 60$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for $H^+ \rightarrow W a$ with $m_{H^{+}} = 120$ GeV, $m_{a} = 30$ GeV for the muon channel $\mu\mu\mu$. A cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Total yields in the $e\mu\mu$ and $\mu\mu\mu$ signal regions for the expected background and data after the profile likelihood fit for a $m_a = 30$ GeV signal mass hypothesis.

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Fitted event-level b-tagging selection efficiency as a function of the dijet invariant mass. The uncertainties are evaluated by propagating the uncertainties on the fit parameters to the fitted efficiency.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 1.5, 1.0, and 0.75 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 0.5, 0.35, and 0.25 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Expected and observed 95% CL for the 1SFH mu Majorana model.

Expected and observed 95% CL for the 2QDH NH Dirac model.

Expected and observed 95% CL for the 2QDH NH Majorana model.

Expected and observed 95% CL for the 2QDH IH Dirac model.

Expected and observed 95% CL for the 2QDH IH Majorana model.

Cutflow for six simulated signal channels showing the weighted number of expected events based on the single-flavour mixing model in the Majorana limit. Each column uses the generated signal sample with the mass hypothesis $m_N = 10$ GeV and proper decay length $c\tau_N = 10$ mm.

Cutflow for the six channels in data showing the number of events passing each successive signal selection for Majorana HNLs.

The event selection efficiency for each mass-lifetime point in all six studied channels. Shown is the fraction of the produced MC simulation events that pass all signal region selections. An entry of 0 indicates no events were selected.

The dominant signal uncertainty is due to differences in reconstruction of displaced vertices and tracks between data and MC. This is evaluated by comparing $K^{0}_{S} \rightarrow \pi^+\pi^-$ event yields in the VR and in MC produced with Pythia8.186 in bins of $p_\mathrm{T}$ and $r_\mathrm{DV}$. The data/MC ratio is normalized to the bin nearest the IP where the tracking and vertexing reconstruction algorithms are expected to be most robust. The symmetrized difference from 1.0 is applied to each signal vertex as a per-event systematic variation.

Expected and observed yields in the different analysis regions (prefit) for the 1SFH e Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH e Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH e Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH e Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH u Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH u Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH u Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH u Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (NH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (NH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (NH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (NH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (IH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (IH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (IH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (IH) Majorana model (10 GeV, 10mm).

The total displaced vertexing efficiency as a function of $r_{DV}$ for the custom configuration used in this analysis. The definition of the secondary vertex efficiency can be found in defined in \cite{ATL-PHYS-PUB-2019-013}. The efficiency is shown for $\mu-\mu\mu$, $\mu-\mu e$ and $\mu-ee$ signals with $m_N=10$~GeV and $c\tau_N=10$~mm.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.3< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.3< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.3< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$Cov_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$Cov_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$Cov_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$Cov_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ for peripheral events, Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality,

$\rho_{2}$ for peripheral events, Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality,

$\rho_{2}$ for peripheral events, Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for peripheral events, Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for peripheral events, Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality,

$\rho_{2}$ for peripheral events, Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality,

$\rho_{2}$ for peripheral events, Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for peripheral events, Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for peripheral events, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality,

$\rho_{2}$ for peripheral events, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality,

$\rho_{2}$ for peripheral events, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for peripheral events, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$, Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$, Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$, Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$, Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$, Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$, Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for central events, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for central events, Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for central events, Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ for central events, Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Three_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Three_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\Sigma E_{T}$ vs $N^{rec}_{ch}$ for Pb+Pb 5.02 TeV

$\Sigma E_{T}$ vs $N^{rec}_{ch}$ for Xe+Xe 5.44 TeV

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Standard method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Standard method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Combined_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Combined_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Three_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Three_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Combined_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality,

$\rho_{2}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Combined_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality,

$\rho_{3}$ for central events, Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ for central events, Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ for central events, Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ for central events, Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Standard method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Standard method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Combined_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{3}$ ratio between Xe+Xe 5.44 TeV and Pb+Pb 5.02 TeV for central events, Combined_subevent method, for , $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality,

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$\rho_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{2}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{2}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{3}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{3}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{3}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{4}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{4}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{4}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{3}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Two_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$Cov_{4}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$ based Centrality.

$c_{k}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$c_{k}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{2})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{2})$ Combined subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{3})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{3})$ Combined subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{4})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{4})$ Combined subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$c_{k}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$c_{k}$ Standard method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{2})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{2})$ Combined subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{3})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{3})$ Combined subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{4})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{4})$ Combined subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$c_{k}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$c_{k}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{2})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{2})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{3})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{3})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{4})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$var(v^{2}_{4})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$.

$c_{k}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$c_{k}$ Standard method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{2})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{2})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{3})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{3})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{4})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$var(v^{2}_{4})$ Combined subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N^{rec}_{ch}$.

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Pb+Pb 5.02 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<2.5, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <2.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{2}$ Three_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{3}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $\Sigma E_{T}$ based Centrality

$\rho_{4}$ Combined_subevent method, for Xe+Xe 5.44 TeV, $|\eta|$<1.0, 0.5< $p_{T}$ <5.0 GeV vs $N_{ch}^{rec}$ based Centrality