Inclusive J/$\psi$ $v_2$ as a function of $p_{T}$ in the centrality interval 0$-$50 %

Inclusive J/$\psi$ $v_2$ as a function of $p_{T}$ in the centrality interval 0$-$10 %

Inclusive J/$\psi$ $v_2$ as a function of $p_{T}$ in the centrality interval 30$-$50 %

Inclusive J/$\psi$ $v_3$ as a function of $p_{T}$ in the centrality interval 0$-$10 %

Inclusive J/$\psi$ $v_3$ as a function of $p_{T}$ in the centrality interval 10$-$30 %

Inclusive J/$\psi$ $v_3$ as a function of $p_{T}$ in the centrality interval 30$-$50 %

Inclusive J/$\psi$ $v_3$ as a function of $p_{T}$ in the centrality interval 0$-$50 %

Inclusive J/$\psi$ $v_3/v_2$ as a function of $p_{T}$ in the centrality interval 0$-$50 %

Inclusive J/$\psi$ $v_3/v_2^{3/2}$ as a function of $p_{T}$ in the centrality interval 0$-$50 %

Inclusive J/$\psi$ $v_2$ as a function of $p_{T}$ in the centrality interval 20$-$40 %

Inclusive J/$\psi$ $v_2$ as a function of centrality in the $p_{T}$ range 0$-$5 GeV/$c$

Inclusive J/$\psi$ $v_2$ as a function of centrality in the $p_{T}$ range 5$-$20.0 GeV/$c$

Inclusive J/$\psi$ $v_2^{\pi}/v_2^{J/\psi}$ as a function of centrality in the $p_{T}$ range 0$-$5 GeV/$c$

Inclusive J/$\psi$ $v_2^{\pi}/v_2^{J/\psi}$ as a function of centrality in the $p_{T}$ range 5$-$20.0 GeV/$c$

Inclusive J/$\psi$ $v_3$ as a function of centrality in the $p_{T}$ range 0$-$5 GeV/$c$

Inclusive J/$\psi$ $v_3$ as a function of centrality in the $p_{T}$ range 5$-$20.0 GeV/$c$

Average $v_3${SP} vs. $p_\mathrm{T}$ of prompt D$^0$, D$^+$, and D$^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% in the rapidity interval $|y|<0.8$. The first (sys,data) error is the systematic uncertainty from all the other sources except for the B feed-down. The second (sys,FD) error is the systematic uncertainty from the B feed-down contribution.

Average $v_2${SP} vs. $p_\mathrm{T}$ of prompt D$^0$, D$^+$, and D$^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 0-10% and 20% small-$q_2$ sample, in the rapidity interval $|y|<0.8$. The first (sys,data) error is the systematic uncertainty from all the other sources except for the B feed-down. The second (sys,FD) error is the systematic uncertainty from the B feed-down contribution.

Average $v_2${SP} vs. $p_\mathrm{T}$ of prompt D$^0$, D$^+$, and D$^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 0-10% and 20% large-$q_2$ sample, in the rapidity interval $|y|<0.8$. The first (sys,data) error is the systematic uncertainty from all the other sources except for the B feed-down. The second (sys,FD) error is the systematic uncertainty from the B feed-down contribution.

Average $v_2${SP} vs. $p_\mathrm{T}$ of prompt D$^0$, D$^+$, and D$^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% and 20% small-$q_2$ sample, in the rapidity interval $|y|<0.8$. The first (sys,data) error is the systematic uncertainty from all the other sources except for the B feed-down. The second (sys,FD) error is the systematic uncertainty from the B feed-down contribution.

Average $v_2${SP} vs. $p_\mathrm{T}$ of prompt D$^0$, D$^+$, and D$^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% and 20% large-$q_2$ sample, in the rapidity interval $|y|<0.8$. The first (sys,data) error is the systematic uncertainty from all the other sources except for the B feed-down. The second (sys,FD) error is the systematic uncertainty from the B feed-down contribution.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $v_2${SP} in the 20% small-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 0-10% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $v_2${SP} in the 20% large-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 0-10% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $v_2${SP} in the 20% small-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $v_2${SP} in the 20% large-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $\mathrm{d}N/\mathrm{d}p_\mathrm{T}$ in the 20% small-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 0-10% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $\mathrm{d}N/\mathrm{d}p_\mathrm{T}$ in the 20% large-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 0-10% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $\mathrm{d}N/\mathrm{d}p_\mathrm{T}$ in the 20% small-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% in the rapidity interval $|y|<0.8$.

Ratio of average prompt D$^0$, D$^+$, and D$^{*+}$ meson $\mathrm{d}N/\mathrm{d}p_\mathrm{T}$ in the 20% large-$q_2$ sample and the unbiased sample vs. $p_\mathrm{T}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in the centrality class 30-50% in the rapidity interval $|y|<0.8$.

$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.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic measured with the two-particle cumulants as a function of transverse momentum in centrality bin 60-80%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 0-2%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 60-80%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 60-80%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic measured with the six-particle cumulats as a function of transverse momentum in centrality bin 60-80%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic measured with the eight-particle cumulats as a function of transverse momentum in centrality bin 60-80%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 40-50%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 10-20%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 20-30%.

The second flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 30-40%.

The triangular flow harmonic measured with the two-particle cumulats as a function of transverse momentum in centrality bin 0-25%.

The triangular flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 0-25%.

The triangular flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 0-25%.

The triangular flow harmonic measured with the two-particle cumulats as a function of transverse momentum in centrality bin 25-60%.

The triangular flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 25-60%.

The triangular flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 25-60%.

The quadrangular flow harmonic measured with the two-particle cumulats as a function of transverse momentum in centrality bin 0-25%.

The quadrangular flow harmonic measured with the Event Plane method as a function of transverse momentum in centrality bin 0-25%.

The quadrangular flow harmonic measured with the four-particle cumulats as a function of transverse momentum in centrality bin 0-25%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 0-2%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 2-5%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 5-10%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 10-15%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 15-20%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 20-25%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 25-30%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 30-35%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 35-40%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 40-45%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 45-50%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 50-55%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 55-60%.

The second flow harmonic measured with the two-particle cumulants as a function of pseudorapidity in centrality bin 60-80%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 0-2%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 2-5%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 5-10%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 10-15%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 15-20%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 20-25%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 25-30%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 30-35%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 35-40%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 40-45%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 45-50%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 50-55%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 55-60%.

The second flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 60-80%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 2-5%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 5-10%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 10-15%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 15-20%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 20-25%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 25-30%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 30-35%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 35-40%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 40-45%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 45-50%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 50-55%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 55-60%.

The second flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 60-80%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 2-5%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 5-10%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 10-15%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 15-20%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 20-25%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 25-30%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 30-35%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 35-40%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 40-45%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 45-50%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 50-55%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 55-60%.

The second flow harmonic measured with the six-particle cumulats as a function of pseudorapidity in centrality bin 60-80%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 2-5%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 5-10%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 10-15%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 15-20%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 20-25%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 25-30%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 30-35%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 35-40%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 40-45%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 45-50%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 50-55%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 55-60%.

The second flow harmonic measured with the eight-particle cumulats as a function of pseudorapidity in centrality bin 60-80%.

The triangular flow harmonic measured with the two-particle cumulats as a function of pseudorapidity in centrality bin 0-60%.

The triangular flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 0-60%.

The triangular flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 0-60%.

The quadrangular flow harmonic measured with the two-particle cumulats as a function of pseudorapidity in centrality bin 0-25%.

The quadrangular flow harmonic measured with the Event Plane method as a function of pseudorapidity in centrality bin 0-25%.

The quadrangular flow harmonic measured with the four-particle cumulats as a function of pseudorapidity in centrality bin 0-25%.

The second flow harmonic measured with the two-particle cumulats as a function of <Npart>.

The second flow harmonic measured with the four-particle cumulats as a function of <Npart>.

The second flow harmonic measured with the six-particle cumulats as a function of <Npart>.

The second flow harmonic measured with the eight-particle cumulats as a function of <Npart>.

The ratio of second flow harmonics measured with the six- and four-particle cumulants as a function of <Npart>.

The ratio of second flow harmonics measured with the eight- and four-particle cumulants as a function of <Npart>.

The second flow harmonic measured with the Event Plane method as a function of <Npart>.

The triangular flow harmonic measured with the Event Plane method as a function of <Npart>.

The triangular flow harmonic measured with the two-particle cumulants as a function of <Npart>.

The triangular flow harmonic measured with the two-particle cumulants as a function of <Npart>.

The quadrangular flow harmonic measured with the Event Plane method as a function of <Npart>.

The quadrangular flow harmonic measured with the two-particle cumulants as a function of <Npart>.

The quadrangular flow harmonic measured with the two-particle cumulants as a function of <Npart>.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic fluctiuations, F(v2), as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic fluctuations, F(v2), as a function of <Npart>.

The triangular flow harmonic fluctuations, F(v3), as a function of <Npart>.

The triangular flow harmonic fluctuations, F(v4), as a function of <Npart>.

The second flow harmonic measured with the two-particle cumulats as a function of <Npart>.

The second flow harmonic measured with the four-particle cumulats as a function of <Npart>.

The second flow harmonic measured with the six-particle cumulats as a function of <Npart>.

The second flow harmonic measured with the eight-particle cumulats as a function of <Npart>.

The ratio of second flow harmonics measured with the six- and four-particle cumulants as a function of <Npart>.

The ratio of second flow harmonics measured with the eight- and four-particle cumulants as a function of <Npart>.

The triangular flow harmonic measured with the two-particle cumulants as a function of <Npart>.

The quadrangular flow harmonic measured with the Event Plane method as a function of <Npart>.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{EP} and v2{4}, as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 2-5%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 5-10%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 10-15%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 15-20%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 20-25%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 25-30%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 30-35%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 35-40%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 40-45%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 45-50%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 50-55%.

The second flow harmonic fluctiuations, F(v2), calculated from v2{2} and v2{4}, as a function of transverse momentum in centrality bin 55-60%.

The second flow harmonic fluctuations, F(v2), as a function of <Npart>.

The triangular flow harmonic fluctuations, F(v3), as a function of <Npart>.

The triangular flow harmonic fluctuations, F(v4), as a function of <Npart>.

The transverse momentum, $p_{T}$, dependence of the pixel track (PXT) reconstruction efficiency for three pseudorapidity ranges in 0-10% centrality interval.

The transverse momentum, $p_{T}$, dependence of the pixel track (PXT) reconstruction efficiency for three pseudorapidity ranges in 40-50% centrality interval.

The transverse momentum, $p_{T}$, dependence of the pixel track (PXT) reconstruction efficiency for three pseudorapidity ranges in 70-80% centrality interval.

The transverse momentum, $p_{T}$, dependence of the pixel track (PXT) reconstruction fake rate for three pseudorapidity ranges in 0-10% centrality interval.

The transverse momentum, $p_{T}$, dependence of the pixel track (PXT) reconstruction fake rate for three pseudorapidity ranges in 40-50% centrality interval.

The transverse momentum, $p_{T}$, dependence of the pixel track (PXT) reconstruction fake rate for three pseudorapidity ranges in 70-80% centrality interval.

The transverse momentum, $p_{T}$, dependence of the inner detector track (IDT) reconstruction efficiency for three pseudorapidity ranges in 0-10% centrality interval.

The transverse momentum, $p_{T}$, dependence of the inner detector track (IDT) reconstruction efficiency for three pseudorapidity ranges in 40-50% centrality interval.

The transverse momentum, $p_{T}$, dependence of the inner detector track (IDT) reconstruction efficiency for three pseudorapidity ranges in 70-80% centrality interval.

The transverse momentum, $p_{T}$, dependence of the inner detector track (IDT) reconstruction fake rate for three pseudorapidity ranges in 0-10% centrality interval.

The transverse momentum, $p_{T}$, dependence of the inner detector track (IDT) reconstruction fake rate for three pseudorapidity ranges in 40-50% centrality interval.

The transverse momentum, $p_{T}$, dependence of the inner detector track (IDT) reconstruction fake rate for three pseudorapidity ranges in 70-80% centrality interval.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 0-10% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 10-20% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 20-30% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 30-40% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 40-50% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 50-60% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 60-70% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Elliptic flow $v_{2}$ integrated over transverse momentum $p_{T}>p_{T,0}$ as a function of $p_{T,0}$ for 70-80% centrality interval, obtained with different charged-particle reconstruction methods: the tracklet (TKT) method with $p_{T,0}=0.07$ GeV, the pixel track (PXT) method with $p_{T,0} \geq 0.1$ GeV and the ID track (IDT) method with $p_{T,0}=0.5$ GeV. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 0-10% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 10-20% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 20-30% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 30-40% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 40-50% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 50-60% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 60-70% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Pseudorapidity dependence of elliptic flow, $v_{2}$, integrated over transverse momentum, $p_{T}$, for different charged particle reconstruction methods and different low-$p_{T}$ thresholds for the 70-80% centrality interval. Error bars indicate statistical and systematic uncertainties added in quadrature.

Integrated elliptic flow, $v_{2}$, as a function of $|\eta| - y_{beam}$ for three centrality intervals Error bars indicate statistical and systematic uncertainties added in quadrature.

The transverse momentum, $p_{T}$, dependence of the TKT track reconstruction efficiency for $\pi^{\pm}$, $K^{\pm}$ and $p^{\pm}$ in the pseudorapidity range $|\eta| < 1$ for 0-10% centrality interval.

The transverse momentum, $p_{T}$, dependence of the TKT track reconstruction efficiency for $\pi^{\pm}$, $K^{\pm}$ and $p^{\pm}$ in the pseudorapidity range $|\eta| < 1$ for 40-50% centrality interval.

The transverse momentum, $p_{T}$, dependence of the TKT track reconstruction efficiency for $\pi^{\pm}$, $K^{\pm}$ and $p^{\pm}$ in the pseudorapidity range $|\eta| < 1$ for 70-80% centrality interval.

The transverse momentum, $p_{T}$, dependence of the PXT track reconstruction efficiency for $\pi^{\pm}$, $K^{\pm}$ and $p^{\pm}$ in the pseudorapidity range $|\eta| < 1$ for 0-10% centrality interval.

The transverse momentum, $p_{T}$, dependence of the PXT track reconstruction efficiency for $\pi^{\pm}$, $K^{\pm}$ and $p^{\pm}$ in the pseudorapidity range $|\eta| < 1$ for 40-50% centrality interval.

The transverse momentum, $p_{T}$, dependence of the PXT track reconstruction efficiency for $\pi^{\pm}$, $K^{\pm}$ and $p^{\pm}$ in the pseudorapidity range $|\eta| < 1$ for 70-80% centrality interval.

The MEAN(Npart) dependence of SIGMA(V2)/MEAN(V2) for three pT ranges together with the total systematic uncertainties.

The MEAN(Npart) dependence of MEAN(V3) for three pT ranges together with the total systematic uncertainties.

The MEAN(Npart) dependence of SIGMA(V3) for three pT ranges together with the total systematic uncertainties.

The MEAN(Npart) dependence of SIGMA(V3)/MEAN(V3) for three pT ranges together with the total systematic uncertainties.

The MEAN(Npart) dependence of MEAN(V4) for three pT ranges together with the total systematic uncertainties.

The MEAN(Npart) dependence of SIGMA(V4) for three pT ranges together with the total systematic uncertainties.

The MEAN(Npart) dependence of SIGMA(V4)/MEAN(V4) for three pT ranges together with the total systematic uncertainties.

Eccentricity curves for EPSILON2 in Figure 12.

Eccentricity curves for EPSILON3 in Figure 12.

Eccentricity curves for EPSILON4 in Figure 12.

Comparison of MEAN(V2) and SQRT(MEAN(V2**2)), derived from the EbyE V2 distributions, with the V2(EP), for charged particles in the pT > 0.5 GeV range.

The ratios of SQRT(MEAN(V2**2)) and V2(EP) to MEAN(V2), for charged particles in the pT > 0.5 GeV range.

Comparison of MEAN(V3) and SQRT(MEAN(V3**2)), derived from the EbyE V3 distributions, with the V3(EP), for charged particles in the pT > 0.5 GeV range.

The ratios of SQRT(MEAN(V3**2)) and V3(EP) to MEAN(V3), for charged particles in the pT > 0.5 GeV range.

Comparison of MEAN(V4) and SQRT(MEAN(V4**2)), derived from the EbyE V4 distributions, with the V4(EP), for charged particles in the pT > 0.5 GeV range.

The ratios of SQRT(MEAN(V4**2)) and V4(EP) to MEAN(V4), for charged particles in the pT > 0.5 GeV range.

Comparison of MEAN(V2) and SQRT(MEAN(V2**2)), derived from the EbyE V2 distributions, with the V2(EP), for charged particles in the 0.5 < pT < 1 GeV range.

The ratios of SQRT(MEAN(V2**2)) and V2(EP) to MEAN(V2), for charged particles in the 0.5 < pT < 1 GeV range.

Comparison of MEAN(V3) and SQRT(MEAN(V3**2)), derived from the EbyE V3 distributions, with the V3(EP), for charged particles in the 0.5 < pT < 1 GeV range.

The ratios of SQRT(MEAN(V3**2)) and V3(EP) to MEAN(V3), for charged particles in the 0.5 < pT < 1 GeV range.

Comparison of MEAN(V4) and SQRT(MEAN(V4**2)), derived from the EbyE V4 distributions, with the V4(EP), for charged particles in the 0.5 < pT < 1 GeV range.

The ratios of SQRT(MEAN(V4**2)) and V4(EP) to MEAN(V4), for charged particles in the 0.5 < pT < 1 GeV range.

Comparison of MEAN(V2) and SQRT(MEAN(V2**2)), derived from the EbyE V2 distributions, with the V2(EP), for charged particles in the pT > 1 GeV range.

The ratios of SQRT(MEAN(V2**2)) and V2(EP) to MEAN(V2), for charged particles in the pT > 1 GeV range.

Comparison of MEAN(V3) and SQRT(MEAN(V3**2)), derived from the EbyE V3 distributions, with the V3(EP), for charged particles in the pT > 1 GeV range.

The ratios of SQRT(MEAN(V3**2)) and V3(EP) to MEAN(V3), for charged particles in the pT > 1 GeV range.

Comparison of MEAN(V4) and SQRT(MEAN(V4**2)), derived from the EbyE V4 distributions, with the V4(EP), for charged particles in the pT > 1 GeV range.

The ratios of SQRT(MEAN(V4**2)) and V4(EP) to MEAN(V4), for charged particles in the pT > 1 GeV range.

Bessel-Gaussian fit parameters from Eq. (1.4) and total errors.

The dependence of MEAN(V2) and V2(RP) on MEAN(Npart).

The dependence of SIGMA(V2) and DELTA(V2) on MEAN(Npart).

The dependence of SIGMA(V2) / MEAN(V2) and DELTA(V2) / V2(RP) on MEAN(Npart).

Comparison of the V2(RP) obtained from the Bessel-Gaussian fit of the V2 distributions with the values for two-particle (V2(calc){2}), four-particle (V2(calc){4}), six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants calculated directly from the unfolded V2 distributions.

The ratios of the four-particle (V2(calc){4}), six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants to the fit results (V2(RP)), with the total uncertainties.

The ratios of the six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants to the four-particle (V2(calc){4}) cumulants, with the total uncertainties.

Comparison of the V3(RP) obtained from the Bessel-Gaussian fit of the V3 distributions with the values for two-particle (V3(calc){2}), four-particle (V3(calc){4}), six-particle (V3(calc){6}) and eight-particle (V3(calc){8}) cumulants calculated directly from the unfolded V3 distributions.

The ratios of the four-particle (V3(calc){4}), six-particle (V3(calc){6}) and eight-particle (V3(calc){8}) cumulants to the fit results (V3(RP)), with the total uncertainties.

The ratios of the six-particle (V3(calc){6}) and eight-particle (V3(calc){8}) cumulants to the four-particle (V3(calc){4}) cumulants, with the total uncertainties.

The standard deviation (SIGMA(V2)), the width obtained from Bessel-Gaussian function (DELTA(V2)), the width F1 = SQRT( ( V2(calc){2}**2 - V2(calc){4}**2 ) / 2 ) estimated from the two-particle cumulant (V2(calc){2}) and four-particle cumulant (V2(calc){4}), where these cumulants are calculated analytically via Eq. (5.3) from the V2 distribution.

Various estimates of the relative fluctuations given as SIGMA(V2) / MEAN(V2), DELTA(V2) / V2(RP), F2 = SQRT( ( V2(calc){2}**2 - V2(calc){4}**2) / ( 2*V2(calc){4}**2 ) ) and F3 = SQRT( ( V2(calc){2}**2 - V2(calc){4}**2) / ( V2(calc){2}**2 + V2(calc){4}**2 ) ).

Comparison in 0.5 < pT < 1 GeV of the V2(RP) obtained from the Bessel-Gaussian fit of the V2 distributions with the values for two-particle (V2(calc){2}), four-particle (V2(calc){4}), six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants calculated directly from the unfolded V2 distributions.

The ratios for 0.5 < pT < 1 GeV of the four-particle (V2(calc){4}), six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants to the fit results (V2(RP)), with the total uncertainties.

The ratios for 0.5 < pT < 1 GeV of the six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants to the four-particle (V2(calc){4}) cumulants, with the total uncertainties.

Comparison in pT > 1 GeV of the V2(RP) obtained from the Bessel-Gaussian fit of the V2 distributions with the values for two-particle (V2(calc){2}), four-particle (V2(calc){4}), six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants calculated directly from the unfolded V2 distributions.

The ratios for pT > 1 GeV of the four-particle (V2(calc){4}), six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants to the fit results (V2(RP)), with the total uncertainties.

The ratios for pT > 1 GeV of the six-particle (V2(calc){6}) and eight-particle (V2(calc){8}) cumulants to the four-particle (V2(calc){4}) cumulants, with the total uncertainties.

The values of V2(RP) and V2(RP,obs) obtained from the Bessel-Gaussian fits to the V2 and V2(obs) distributions, with the statistical uncertainties.

The values of DELTA(V2) and DELTA(V2,obs) obtained from the Bessel-Gaussian fits to the V2 and V2(obs) distributions, with the statistical uncertainties.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.