The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2^{1/2}$ difference between quadruplets at ($\eta_{\alpha}$, $\eta_{\alpha}$, $\eta_{\beta}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, $\eta_{\alpha}$, -$\eta_{\beta}$, -$\eta_{\beta}$).

The two-particle cumulants, V2{2} as a function of $\eta$ for one particle while averaged over $\eta$ of the partner particle.

The two-particle cumulants, V2{2} as a function of $\eta$.

The two-particle cumulants, V2{2} as a function of $\eta$ for one particle while averaged over $\eta$ of the partner particle.

The difference in two-particle cumulants, V2{2} as a function of $\eta$ for one particle while averaged over $\eta$ of the partner particle.

No description provided.

V3{2} as a function of $\eta$.

<$V_3^2$> as a function of $\eta$.

The $\Delta\eta$-dependent component of the two-particle cumulant with $\Delta\eta$-gap, $D^-$ in Eq. (11), of the second harmonic is shown as a function of ??-gap |$\Delta\eta$| > x.

The $\Delta\eta$-dependent component of the two-particle cumulant with $\Delta\eta$-gap, $D^-$ in Eq. (11), of the third harmonic is shown as a function of ??-gap |$\Delta\eta$| > x.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $\sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the second harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $\sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the second harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the second harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the third harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $sqrt{\delta_2}$, $\sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the third harmonic.

The relative elliptic flow fluctuation $\sigma_2/<v_2>$ centrality dependence in $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions.

$v_{2,2}^{unsub}$ and $v_{2,2}$ as a function of $\Delta\eta$ calculated from the 2-D per-trigger yields in figure 4(a) and 4(b), respectively.

$v_{3,3}^{unsub}$ and $v_{3,3}$ as a function of $\Delta\eta$ calculated from the 2-D per-trigger yields in figure 4(a) and 4(b), respectively.

$v_{4,4}^{unsub}$ and $v_{4,4}$ as a function of $\Delta\eta$ calculated from the 2-D per-trigger yields in figure 4(a) and 4(b), respectively.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

The integrated per-trigger yield, Y_{int}, on the near-side, the away-side and their difference and Y_{int} from the recoil as a function of event activity. Errors are statistical.

The integrated per-trigger yield, Y_{int}, on the near-side, the away-side and their difference and Y_{int} from the recoil as a function of event activity. Errors are statistical.

The Fourier coefficients $v_{n}$ as a function of $p_{T}^{a}$ extracted from the correlation functions, before and after the subtraction of the recoil component.

The Fourier coefficients $v_{n}$ as a function of $p_{T}^{a}$ extracted from the correlation functions, before and after the subtraction of the recoil component.

The Fourier coefficients $v_{n}$ as a function of $p_{T}^{a}$ extracted from the correlation functions, before and after the subtraction of the recoil component.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The centrality dependence of $v_{2}$ as a function of $N_{ch}^{rec}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{3}$ as a function of $N_{ch}^{rec}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{4}$ as a function of $N_{ch}^{rec}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{2}$ as a function of $E_{T}^{Pb}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{3}$ as a function of $E_{T}^{Pb}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{4}$ as a function of $E_{T}^{Pb}$. Values from before and after the recoil subtraction are included.

The $v_{2}$ as a function of $E_{T}^{Pb}$ obtained indirectly by mapping from the $N_{ch}^{rec}-dependence of $v_{2}$ using the correlation data shown in Fig. 2(b).

The $v_{3}$ as a function of $E_{T}^{Pb}$ obtained indirectly by mapping from the $N_{ch}^{rec}-dependence of $v_{3}$ using the correlation data shown in Fig. 2(b).

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic $v_1$ obtained using factorization from $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic $v_1$ obtained using factorization from $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic $v_1$ obtained using factorization from $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

$v_{2}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method.

$v_{2}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method, after the scaling.

$v_{3}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method.

$v_{3}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method, after the scaling.

$v_{4}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method.

$v_{4}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method, after the scaling.

Correlation between $E_{T}^{FCal}$ and $N_{ch}^{rec}$ for MB events (without weighting) and MB+HMT events (with weighting), errors are statistical.

Reconstruction efficiency of TYPE0 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE0 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE0 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $p_{\mathrm{T}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.05 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $p_{\mathrm{T}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.15 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $p_{\mathrm{T}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $p_{\mathrm{T}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $p_{\mathrm{T}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.05 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.15 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of $L_{\mathrm{xy}}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. The uncertainties are statistical only.

Reconstruction efficiency of TYPE0 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.4 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE0 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.9 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE1 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.4 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE1 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.9 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.05 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.15 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.4 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Reconstruction efficiency of TYPE2 LJs as a function of the $p_{\mathrm{T}}$ of the $s_{d_{1}}$ for LJs with two $\gamma_{d}$'s for an \scalar mass of 2 GeV. For the $\gamma_{d}$, the kinematically allowed mass of 0.9 GeV is considered. The distributions for the other $s_{d_{1}}$ masses are very similar. The uncertainties are statistical only.

Muon trigger efficiency, $\rm \varepsilon$(2mu6), as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. The uncertainties are statistical only.

Muon trigger efficiency, $\rm \varepsilon$(2mu6), as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. The uncertainties are statistical only.

Muon trigger efficiency, $\rm \varepsilon$(2mu6), as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. The uncertainties are statistical only.

Muon trigger efficiency, $\rm \varepsilon$(2mu6), as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. The uncertainties are statistical only.

Muon trigger efficiency, $\rm \varepsilon$(2mu6), as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. The uncertainties are statistical only.

Muon trigger efficiency, $\rm \varepsilon$(2mu6), as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to \mu\mu$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.05 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.15 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $p_{T}$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.05 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.15 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.4 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 0.9 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

Calorimetric trigger efficiency as a function of $\eta$ of the $\gamma_{d}$ for $\gamma_{d} \to ee,\pi\pi$ obtained from the LJ gun MC samples with $\gamma_{d}$ masses 1.5 GeV. Similar distributions are obtained for LJs containing two dark photons. The uncertainties are statistical only.

The experimental sensitivity and observed limit on the number of dark photon candidates as a function of the assumed dark photon mass.

Limits on the $U-\gamma$ mixing parameters from PHENIX at 90%, 85% and 80% confidence levels.

Distribution of exclusive jet multiplicity.

Distribution of pT (leading jet) [GeV] with at least one jet in the event.

Distribution of pT (leading jet) [GeV] with exactly one jet in the event.

Distribution of pT (leading jet) [GeV] with at least two jets in the event.

Distribution of pT (leading jet) [GeV] with at least three jets in the event.

Distribution of pT (2nd jet) [GeV] with at least two jets in the event.

Distribution of pT (3rd jet) [GeV] with at least three jets in the event.

Distribution of ST [GeV] with at least two jets in the event.

Distribution of ST [GeV] with exactly two jets in the event.

Distribution of ST [GeV] with at least three jets in the event.

Distribution of ST [GeV] with exactly three jets in the event.

Distribution of DR(jj) [GeV] with at least two jets in the event.

Distribution of DPhi(jj) [GeV] with at least two jets in the event.

Distribution of m(jj) [GeV] with at least two jets in the event.

Distribution of leading jet rapidity with at least one jet in the event.

Distribution of 2nd jet rapidity with at least two jets in the event.

Distribution of 3rd jet rapidity with at least three jets in the event.

The total inelastic cross section, the first systematic error accounts for all experimental uncertainties and the second error for the extrapolation t-->0.

The measured differential elastic cross section. In addition to the statistical and total systematic uncertainties, the following 24 systematic shifts are given, which are included in the profile fit with their signs: -- Constraints: Beam optics uncertainty obtained by varying the ALFA constraints in the optics fit -- QScan: Variation by +/- 0.1 % of the quadrupole strength -- Q2: Fit of the strength of Q2 using the best value for the strength of Q1 and Q3 -- MadX: Uncertainty related to the beam transport replacing matrix transport by MadX PTC tracking -- Q5Q6: Variation of the strength of Q5 and Q6 by -0.2% as indicated by machine constraints -- Qmisal: Uncertainty due to the mis-alignment of the quadrupoles in the beam line -- Q1Q3: Propagation of the optics fit uncertainty in the strenght of Q1 and Q3 on the differential elastic cross section -- Stat2: Alignment uncertainty from the choice of a reference station -- Dist: Alignment uncertainty related to the distance calibration between the upper and lower detectors -- Leff: Alignment uncertainty related to effective lever arm used in the alignment optimization procedure -- Offv: Alignment uncertainty related to the vertical beam center offset -- Offh: Alignment uncertainty related to the horizontal beam center offset -- Ang: Alignment uncertainty related to the detector rotation in the x-y plane -- BGn: Uncertainty from the background normalization -- BGs: Uncertainty from the background shape -- MCres: Error from modelling of the detector response -- Slope: Residual dependence on the physics model estimated by varying the nuclear slope in the simulation by +/- 1 GeV^-2 -- Emit: Uncertainty from the emittance used to calculate beam divergence in the simulation -- Unf: Unfolding uncertainty from the data-driven closure test -- Trac: Uncertainty from the variation of the track reconstruction selection cuts -- Xing: Uncertainty from residual crossing angle in the horizontal plane -- Eff: Uncertainty from the reconstruction efficiency -- Lumi: Luminosity uncertainty (+/- 2.3%) -- Ebeam: Uncertainty from the nominal beam energy (+/- 0.65%) A small difference in the statistical uncertainties give here compared to the published version is related to insignificant rounding issues.

The statistical covariance matrix for the differential elastic cross section.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $Corr_{\gamma\,-ch}^{real}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $Corr_{\gamma\,-ch}^{mixed}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $\omega_{\gamma}^{real}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $\omega_{\gamma}^{mixed}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $v_{dyn}^{real}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $v_{dyn}^{mixed}$ is plotted.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the same side.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the same side.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the away-side.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the away-side.

The correlation between positive and negative charged particles measured by the FTPC and photons measured by the PMD using $v_{dyn}$ for real events.

$r_{m,1}$ vs multiplicity for first order of m.

$r_{m,1}$ vs multiplicity for first order of m.

$r_{m,1}$ vs multiplicity for second order of m.

$r_{m,1}$ vs multiplicity for second order of m.

$r_{m,1}$ vs multiplicity for third order of m.

$r_{m,1}$ vs multiplicity for third order of m.

The moments $r_{m,1}$(m=1-3) for real and mixed events as a function of order m in a fixed multiplicity bin of 47$\langle \sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }\langle 54$.

Differential production yield per binary collision for $W^{-}$ bosons as a function of $|\eta_\ell|$.

The lepton charge asymmetry $A_{\ell}$ from $W^\pm$ bosons as a function of absolute pseudorapidity.