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.

Data selection efficiency as a function of nuclear recoil energy

Data points used in analysis in log_10(S2)-S1 space

Direct photon spectra for minimum bias (0-86%) Au+Au collision at $\sqrt{s_{NN}} = 62.4$ GeV (a) and $39$ GeV (b). For $62.4$ GeV also centrality bins of 0-20% (c) and 20-40% (d) are shown. Data points are shown with statistical and systematic uncertainties, unless the central value is negative (arrows) or is consistent with zero within the statistical uncertainties (arrows with data point). In these cases upper limit with CL = 95$%$ are given.

Inverse slopes, $T_{eff}$ , obtained from fitting the combined data from central collisions shown in Fig. 11 is compared to the fit results of the individual data sets at $\sqrt{s_{NN}} = 62.4$ GeV, $200$ GeV, and $2760$ GeV. Also included is the value for $\sqrt{s_{NN}} = 39$ GeV obtained from fitting the min. bias data set in the lower $p_T$ range.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated direct-photon yields from A+A collisions for $p_{T,min}$ of $5$ GeV/c (a) and $8$ GeV/c. The representation is the same as in Fig. 13. Also shown are the results from pQCD calculations scaled by $N_{coll}$

Integrated direct-photon yields from A+A collisions for $p_{T,min}$ of $5$ GeV/c (a) and $8$ GeV/c. The representation is the same as in Fig. 13. Also shown are the results from pQCD calculations scaled by $N_{coll}$

The $\alpha$ values extracted using fits to integrated direct photon yields. The dashed line gives the average $\alpha$ value for the 4 lower $p_{T,min}$ points. Also shown is a model calculation for $\alpha$ discussed in the text.

Exponent according to the Eq. 5 as a function of transverse momenta extracted from p+Au/p+Al and $^{3}$HeAu/p+Au collision systems. The uncertainties from the $N_{coll}$ calculations and from the overall normalization of p+p cancel in these ratios

Nuclear modification factors in p+Al, p+Au, d+Au, and $^{3}$He+Au in five centrality bins and for inelastic collisions at $\sqrt{s}$ = 200 GeV. The right boxes are the uncertainties of the $N_{coll}$ collisions from the Glauber model, while the left box represents the overall normalization uncertainty from p+p collisions

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

$R_{xA}$ for inelastic collisions compared to three different nuclear PDF calculations and their uncertainties. The data points include the statistical and systematical uncertainties. Relevant data has been tabulated in Figure 7.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. the number of collisions per projectile participant in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. Average $R_{xA}$ in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Panel(b) and panel(c) are tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

Integrated away-side $\gamma_{dir}$-h per-trigger $I_{AA}$ and $I_{dA}$ of Au+Au, d+Au, and p+p, as a function of $\xi$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

Ratios of $I_{AA}$ as a function of direct photon $p_T$ for three different away-side integration ranges.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Direct photon spectra(Nucl. Part. Phys. 23, A1 (1997)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra(Sov. J. Nucl. Phys. 51, 836 (1990)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 62.4 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 39.0 GeV.

Direct photon spectra(Physical Review C91, 064904 (2015)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 104, 132301 (2010)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 109, 152302 (2012)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review C98, 054902 (2018)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Cu+Cu at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physics Letters B754 235 (2016)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Pb+Pb at $\sqrt{s_{NN}}$= 2760 GeV.

Number of binary collisions N_{coll} vs. Charged particle multiplicity in Au+Au at various cm energies.

Charged particle multiplicity vs. centrality in Au+Au at various cm energies.

Number of binary collisions N_{coll} vs. centrality in Au+Au at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the near side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the away side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $Y(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: (a) the negative of the dipole coefficient, $-c_1$; (b) the quadrupole coefficient $c_2$; (c) the ratio ${-c_2}/{c_1}$.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au and $p$+$p$.

System centrality dependence of correlation coefficients $-c_1$, $c_2$, and the ratio ${-c_2}/{c_1}$ as a function of $N_{part}$.

System centrality dependence of correlation coefficients $-c_1$, $c_2$, and the ratio ${-c_2}/{c_1}$ as a function of $N_{part}$.

Lévy scale parameter $R$ versus average $m_T$ of the pair. The graphical representation of statistical and systematic uncertainties is the same as in Fig. 4.

Lévy index parameter $\alpha$ versus average $m_T$ of the pair. Statistical and systematic uncertainties are indicated similarly to Fig. 4. The horizontal line, %\aplha$ = 1.207, represents the 0%-30% centrality average value of %\aplha$.

Inverse square of the Lévy scale parameter 1/$R^2$ versus average $m_T$ of the pair. Statistical and systematic uncertainties shown as bars and boxes, respectively.

New scale parameter $\hat{R}$ versus average $m_T$ of the pair, with a linear fit. Statistical and systematic uncertainties shown as bars and boxes, respectively.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.