The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

Differential cross-section for $\eta_c(1S)$ from inclusive charmonium production in b-hadrons decays for $p_T > 6.5$ [GeV/$c$] and $2.0 < y < 4.5$. The reported uncertainties are total errors.

Trans.min. scalar pT sum density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Away charged multiplicity density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards <mean pT(charged)> vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse <mean pT(charged)> vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Away <mean pT(charged)> vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Towards <mean pT(charged)> vs charged multiplicity, Born leptons : Statistical and systematic errors are added in quadrature.

Transverse <mean pT(charged)> vs charged multiplicity, Born leptons : Statistical and systematic errors are added in quadrature.

Away <mean pT(charged)> vs charged multiplicity, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.diff. scalar pT sum density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.diff. charged multiplicity density vs Z-boson pT, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Born leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Away scalar pT sum density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Away charged multiplicity density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. charged multiplicity distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards <mean pT(charged)> vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse <mean pT(charged)> vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Away <mean pT(charged)> vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Towards <mean pT(charged)> vs charged multiplicity, dressed leptons : Statistical and systematic errors are added in quadrature.

Transverse <mean pT(charged)> vs charged multiplicity, dressed leptons : Statistical and systematic errors are added in quadrature.

Away <mean pT(charged)> vs charged multiplicity, dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.min. scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.diff. scalar pT sum density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.diff. charged multiplicity density vs Z-boson pT, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 0 - 5 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 5 - 10 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 10 - 20 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 20 - 50 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 50 - 110 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

Trans.max. scalar pT sum distribution, Z-boson pT = 110 - 500 GeV, Dressed leptons : Statistical and systematic errors are added in quadrature.

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.