Expected CLs contour in the GMSB (LAMBDA, TAU(NEUTRALINO) ) signal space.

Expected CLs contour with minus 1-sigma in the GMSB (LAMBDA, TAU(NEUTRALINO) ) signal space.

Expected CLs contour with plus 1-sigma in the GMSB (LAMBDA, TAU(NEUTRALINO) ) signal space.

The total expected and observed limits on SPS8 cross-section [fb] for each signal point.

The total expected and observed limits on number of accepted signal events for each signal point.

The total normalization systematic for each signal point.

Cutflow showing the number of events passing various cuts for the $\Lambda =$ 300 TeV, $\tau =$ 2 ns GMSB SPS8 signal point.

The number of events in the signal region sub-divided into the different bins used in the limit setting procedure for the $\Lambda =$ 300 TeV, $\tau =$ 2 ns GMSB SPS8 signal point. Maximum values of $|t_\gamma| <$ 4 ns and $\Delta z_\gamma<$ 2000 mm have been applied to both photons in the event.

Expected and observed 95% CL cross section exclusion contours for top squark pair production in the plane of top squark lifetime ($c\tau$) and top squark mass. These limits assume a branching fraction of 100\% through the RPV vertex $\tilde{t}$ $\to$ b l, where the branching fraction to any lepton flavor is equal to 1/3. As indicated in the plot, the region to the left of the contours is excluded by this search.

Electron reconstruction efficiency as function of its tranverse impact parameter, $d_0$.

Muon reconstruction efficiency as function of its tranverse impact parameter, $d_0$.

Electron selection efficiency as function of its tranverse momentum, $p_T$.

Muon selection efficiency as function of its tranverse momentum, $p_T$.

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.

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.

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

Distribution of exclusive jet multiplicity.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The statistical covariance matrix for the differential elastic cross section.