J/psi-->e+e- invariant yield in Cu+Cu collisions as a function of p_T at mid-rapidity for the 20-40 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI yield versus transverse momentum PT, at mid rapidity : -0.35<y<0.35, for a centrality range of 40-60%.

J/psi-->e+e- invariant yield in Cu+Cu collisions as a function of p_T at mid-rapidity for the 40-60 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI yield versus transverse momentum PT, at mid rapidity : -0.35<y<0.35, for a centrality range of 60-94%.

J/psi-->e+e- invariant yield in Cu+Cu collisions as a function of p_T at mid-rapidity for the 60-94 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI yield versus transverse momentum PT, at forward rapidity : absolute value of y belongs to [1.22.2], for a centrality range of 0-20%.

J/psi-->mu+mu- invariant yield in Cu+Cu collisions as a function of p_T at forward rapidity for the 0-20, 20-40, 40-60 and 60-94 centrality ranges. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI yield versus transverse momentum PT, at forward rapidity : absolute value of y belongs to [1.22.2], for a centrality range of 20-40%.

J/psi-->e+e- average momentum squared in Cu+Cu collisions as a function of N_part at forward rapidity (p_T integrated). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value.

J/PSI yield versus transverse momentum PT, at forward rapidity : absolute value of y belongs to [1.22.2], for a centrality range of 40-60%.

J/psi-->mu+mu- average momentum squared in Cu+Cu collisions as a function of N_part at forward rapidity (p_T integrated). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value.

J/PSI yield versus transverse momentum PT, at forward rapidity : absolute value of y belongs to [1.22.2], for a centrality range of 60-94%.

J/psi-->mu+mu- and J/psi-->e+e- nuclear modification in Cu+Cu collisions as a function of p_T at forward and mid-rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Mean PT^2 versus Npart, number of participants, at mid rapidity : -0.35<y<0.35.

J/psi-->mu+mu- and J/psi-->e+e- nuclear modification in Cu+Cu collisions as a function of y at forward and mid-rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Mean PT^2 versus Npart, number of participants, at forward rapidity : absolute value of y belongs to [1.22.2].

J/psi-->e+e- nuclear modification in Cu+Cu collisions as a function of N_part at mid-rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Nuclear modification factor RAA versus transverse momentum PT at mid rapidity : -0.35<y<0.35, for the 0-20% most central events.

J/psi-->mu+mu- nuclear modification in Cu+Cu collisions as a function of N_part at forward rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Nuclear modification factor RAA versus transverse momentum PT at forward rapidity : absolute value of y belongs to [1.22.2], for the 0-20% most central events.

J/psi-->mu+mu- in Cu+Cu collisions as a function of $N_{part}$ at forward rapidity expressed as a ratio with respect to J/psi-->e+e- nuclear modification at mid-rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Nuclear modification factor RAA versus rapidity, for the 0-20% most central events.

Nuclear modification factor RAA versus number of participants, at mid rapidity : -0.35<y<0.35.

Nuclear modification factor RAA versus number of participants, at forward rapidity : absolute value of y belongs to [1.22.2].

$\Upsilon(1S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $p_{T}$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

Scaled $pp$ differential cross-section in $p_{T}$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

Scaled $pp$ differential cross-section in $y$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $p_{T}$ integrated over $\vert y^* \vert$ in the range $2.5 < \vert y^* \vert < 4.0$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $\vert y^* \vert$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ to non-prompt $J/\psi$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

Invariant yield of electrons from heavy-flavor decays for 40-60% central collisions, versus PT.

Invariant yield of electrons from heavy-flavor decays for 60-92% central collisions, versus PT.

Invariant yield of electrons from heavy-flavor decays for minimum bias collisions, versus PT.

Nuclear modification factor of heavy-flavor electrons for PT>0.3 Gev/c as a function of Centrality (represented by the number of participants) An additional systematic error of +0.2325,-0.1587 from the measurements in P+P collisions is not included.

Nuclear modification factor of heavy-flavor electrons for PT>3.0 Gev/c as a function of Centrality (represented by the number of participants) An additional systematic error of +0.1103,-0.09038 from the measurements in P+P collisions is not included.

Nuclear modification factor of heavy-flavor electrons versus PT for 0-10% central collisions An additional systematic error of +8%,-6% due to the uncertainty on the overlap integral TAA is to add. This error has been determined using g3data, a software allowing toextract data graphically from plots.

Azimuthal amisotropy (v2) versus PT for minimum bias collisions.

Mean PT^2 value at forward rapidities : absolute value of y belongs to [1.2;2.2] The mean PT is obtained with a phenomonological fit of the J/PSI distribution in PT of the form (1/(2*PI*PT))*D(SIG)/DPT = A ( 1+(PT/B)^2)^-6 .The systematic error includes the incertainty from the maximum shape deviation permitted by the point-to-point correlated errors and from allowing the exponent of the fit fonctionto be a free parameter.

J/PSI differential cross section times dilepton branching ratio versus rapidity.

Total cross-section for J/PSI production. The value of the total cross section has been determined by fitting the rapidity distribution with many theoretical and phenomenological shapes. The systematic error quotedwas estimated from the maximum variation allowed by shifting the mid and forward rapidity data independantely by their point-to-point correlated systematic errors.There is a +- 18 nanobarn normalization error in addition to the statistical and systematical errors shown.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 0-10% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 10-20% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 20-40% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 40-60% centrality.

Spectra in transverse momentum of electrons created in open heavy flavor decays, for 60-92% centrality.

Non photonic electrons differential yield scaled by the number of collisions, as a function of centrality. PT belongs to 0.8-4.0 GeV/c.

Non photonic electrons differential yield scaled by the number of collisions, for minimum bias collisions. PT belongs to 0.8-4.0 GeV/c.

Number of collisions as a fonction of centrality There are only systematic errors due to the calculation of Ncoll.

Number of collisions for minimum bias events There are only systematic errors due to the calculation of Ncoll.

Nuclear overlap function TAA as a function of centrality The axis is labelled with 1/SIG, because the unit of the nuclear overlap function is barn-1 There are only systematic errors due to the calculation of TAA.

Nuclear overlap function TAA for minimum bias events The axis is labelled with 1/SIG, as the unit of the nuclear overlap function is barn-1 There are only systematic errors due to the calculation of TAA.

Differential charm cross section scaled by the number of collisions, at mid rapidity, as a function of centrality The differential cross section per number of collisions D(SIG)/Ncoll is calculated as following : D(SIG)/Ncoll = D(N(cc_bar))/TAA.

Differential charm cross section scaled by the number of collisions, at mid rapidity, for minimum bias events The differential cross section per number of collisions D(SIG)/Ncoll is calculated as following : D(SIG)/Ncoll = D(N(cc_bar))/TAA.

Total charm cross section scaled by the number of collisions, as a function of centrality The total cross section per number of collisions SIG/Ncoll is calculated as following : SIG/Ncoll = N(cc_bar)/TAA.

Total charm cross section scaled by the number of collisions, for minimum bias events The total cross section per number of collisions SIG/Ncoll is calculated as following : SIG/Ncoll = N(cc_bar)/TAA.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus PT at mid rapidity (-0.35<y<0.35) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus PT at forward rapidity (1.2<y<2.2) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus centrality at backward rapidity (-2.2<y<-1.2) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus centrality at mid rapidity (-0.35<y<0.35) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus centrality at forward rapidity (1.2<y<2.2) in D+AU collisions.

Nuclear modification factor versus rapidity in D+AU collisions, over 3 bins of rapidity.

Nuclear modification factor versus rapidity in D+AU collisions, over 5 bins of rapidity.

Nuclear modification factor versus PT at backward rapidity (-2.2<y<-1.2) in D+AU collisions.

Nuclear modification factor versus PT at mid rapidity (-0.35<y<0.35) in D+AU collisions.

Nuclear modification factor versus PT at forward rapidity (1.2<y<2.2) in D+AU collisions.

Nuclear modification factor versus centrality at backward rapidity (-2.2<y<-1.2) in D+AU collisions.

Nuclear modification factor versus centrality at mid rapidity (-0.35<y<0.35) in D+AU collisions.

Nuclear modification factor versus centrality at forward rapidity (1.2<y<2.2) in D+AU collisions.

Mean PT^2 versus rapidity, extracted from a fit to the data, for DEUT + Au reactions. The mean value of PT^2 is calculated from the data for PT < 5 Gev/c.

Mean PT^2 versus rapidity, extracted from a fit to the data, for P + P reactions. The mean value of PT^2 is calculated from the data for PT < 5 Gev/c.

Breakup cross section of c-c_bar pairs inside cold nuclear matter,integrated over the whole domain of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Breakup cross section of c-c_bar pairs inside cold nuclear matter for different ranges of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Breakup cross section of c-c_bar pairs inside cold nuclear matter, integrated over the whole domain of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Breakup cross section of c-c_bar pairs inside cold nuclear matter for different ranges of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Nuclear modification factor as a function of PT, for 10-20% central reactions Note that the systematic error given is related to the the uncertainties in the p+p measurements.An additional systematic error, symmetrical on the + and - side, related to the uncertainties in the Au+Au measurement, is given in the second column. Another, PT-independant, 13%systematic error due to the uncertainty on the overlap function and the Pi0 yield normalization is to add.

Nuclear modification factor as a function of PT, for 20-40% central reactions Note that the systematic error given is related to the the uncertainties in the p+p measurements.An additional systematic error, symmetrical on the + and - side, related to the uncertainties in the Au+Au measurement, is given in the second column. Another, PT-independant, 13%systematic error due to the uncertainty on the overlap function and the Pi0 yield normalization is to add.

Nuclear modification factor as a function of PT, for 40-60% central reactions Note that the systematic error given is related to the the uncertainties in the p+p measurements.An additional systematic error, symmetrical on the + and - side, related to the uncertainties in the Au+Au measurement, is given in the second column. Another, PT-independant, 13%systematic error due to the uncertainty on the overlap function and the Pi0 yield normalization is to add.