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.

Distribution of $m_{\mathrm{JJ}}$ in the SR (black points) compared to the background-only expectation, where the shape is derived from the CR and the normalisation is fitted in the SR.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model A, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model B, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d} = 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model C, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d}= 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model D, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Relative efficiency (in %) of each stage of the analysis selection with respect to the previous one for models A through D and $m_{Z^\prime}=2.5$ TeV. Approximately 90000 events were generated per signal mass hypothesis.

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.10 < $\alpha$ < 0.12.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.12 < $\alpha$ < 0.14.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.14 < $\alpha$ < 0.16.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The product of the analysis acceptance and selection efficiency for all analysis selection criteria is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

The cross-section of the signal samples is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^&beta;=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01 and m_&chi;=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

The Bose-Einstein correlation (BEC) parameter R as a function of n_ch for HMT events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameter R as a function of k_T for MB events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameter &lambda; as a function of k_T for MB events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The two-particle double-ratio correlation function, R₂(Q), for pp collisions for track p_T &gt;100&nbsp;MeV at &radic;s=13&nbsp;TeV in the multiplicity interval 71 &le; n_ch &lt; 80 for the minimum-bias (MB) events. The blue dashed and red solid lines show the results of the exponential and Gaussian fits, respectively. The region excluded from the fits is shown. The statistical uncertainty and the systematic uncertainty for imperfections in the data reconstruction procedure are added in quadrature.

The two-particle double-ratio correlation function, R₂(Q), for pp collisions for track p_T &gt;100&nbsp;MeV at &radic;s=13&nbsp;TeV in the multiplicity interval 231 &le; n_ch &lt; 300 for the high-multiplicity track (HMT) events. The blue dashed and red solid lines show the results of the exponential and Gaussian fits, respectively. The region excluded from the fits is shown. The statistical uncertainty and the systematic uncertainty for imperfections in the data reconstruction procedure are added in quadrature.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

Comparison of single-ratio two-particle correlation functions, using the unlike-charge particle (UCP) pair reference sample, for minimum-bias (MB) events, showing C₂^data(Q) (top panel) at 13&nbsp;TeV (black circles) and 7&nbsp;TeV (open blue circles), and the ratio of C₂^7&nbsp;TeV (Q) to C₂^13&nbsp;TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative multiplicity region 3.09 &lt; m_ch &le; 3.86. The statistical and systematic uncertainties, combined in quadrature, are presented. The systematic uncertainties include track efficiency, Coulomb correction, non-closure and multiplicity-unfolding uncertainties.

Comparison of single-ratio two-particle correlation functions, using the unlike-charge particle (UCP) pair reference sample, for minimum-bias (MB) events, showing C₂^data(Q) (top panel) at 13 &nbsp;TeV (black circles) and 7 &nbsp;TeV (open blue circles), and the ratio of C₂^7&nbsp;TeV (Q) to C₂^13&nbsp;TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative k_T region 400 &lt; k_T &le;500 &nbsp;MeV. The statistical and systematic uncertainties, combined in quadrature, are presented. The systematic uncertainties include track efficiency, Coulomb correction, non-closure and multiplicity-unfolding uncertainties.

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The two-dimensional dependence on m_ch and k_T for p_T &gt; 100&nbsp;MeV for the correlation strength, &lambda;, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The two-dimensional dependence on m_ch and k_T for p_T &gt; 100&nbsp;MeV for the source radius, R, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The two-dimensional dependence on m_ch and k_T for p_T &gt; 500&nbsp;MeV for the correlation strength, &lambda;, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The two-dimensional dependence on m_ch and k_T for p_T &gt; 500&nbsp;MeV for the source radius, R, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of &#8731;n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of &#8731;n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of &#8731;n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

Systematic uncertainties (in percent) in the correlation strength, &lambda;, and source radius, R, for the exponential fit of the two-particle double-ratio correlation functions, R₂(Q), for p_T &gt; 100&nbsp;MeV at &radic;s= 13&nbsp;TeV for the MB and HMT events. The choice of MC generator gives rise to asymmetric uncertainties, denoted by uparrow and downarrow. This asymmetry propagates through to the cumulative uncertainty. The columns under ‘Uncertainty range’ show the range of systematic uncertainty from the fits in the various n_ch intervals.

The results of the fits to the dependencies of the correlation strength, &lambda;, and source radius, R, on the average rescaled charged-particle multiplicity, m_ch, for |&eta;| &lt; 2.5 and both p_T &gt; 100&nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV for the minimum-bias (MB) and the high-multiplicity track (HMT) events. The parameters &gamma; and &delta; resulting from a joint fit to the MB and HMT data are presented. The total uncertainties are shown.

The results of the fits to the dependencies of the correlation strength, &lambda;, and source radius, R, on the pair average transverse momentum, k_T, for various functional forms and for minimum-bias (MB) and high-multiplicity track (HMT) events for p_T &gt; 100&nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV. The total uncertainties are shown.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the multiplicity, m_ch, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 100&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the multiplicity, m_ch, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the pair transverse momentum, k_T, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 100&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the pair transverse momentum, k_T, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 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.

Distributions of $s_{\text{T}}$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $m_{\text{T}}(\tau)$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $\Sigma m_{\text{T}}(b_{1,2})$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $p_{\text{T}}(\tau)$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Observed event yields in data ('Observed') and expected event yields for SM background processes obtained from the background-only fit ('Total bkg.' and rows below) in the signal regions of the di-tau and single-tau channels. The quoted uncertainties include both the statistical and systematic uncertainties and are truncated at zero yield. By construction, no $t\bar{t}$ (2 real $\tau$) events can pass the selections in the single-tau channel. As the individual uncertainties are correlated, they do not add in quadrature to equal the total background uncertainty.

From left to right: upper limits at the 95% confidence level (CL) on the visible cross section ($\sigma_\text{vis}$) and on the number of signal events ($S_{\text{obs}}^{95}$). The third column ($S_{\text{exp}}^{95}$) shows the upper limit at the 95% CL on the number of signal events, given the expected number (and $\pm 1\,\sigma$ excursions on the expectation) of background events. The last two columns indicate the confidence level observed for the background-only hypothesis ($\text{CL}_{b}$), the discovery $p$-value ($p(s=0)$) and the significance ($Z$). In the di-tau SR, where fewer events are observed than predicted by the fitted background estimate, the $p$-value is capped at 0.5.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Upper limits on the signal cross section at the 95 % confidence level for the stop-stau signal model.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with up-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with down-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with minimal coupling (MC).

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with additional gauge couplings (YM).

Acceptance of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the di-tau SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau one-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau multi-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

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

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

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

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

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

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