psi(2S) double-differential cross section times branching fraction and the corresponding scaling factors to obtain the cross sections for different polarization scenarios (azimuthal polarization parameter in the center of mass helicity frame lambda_theta^HX = +1, -1, +0.03) as a function of pT for |y| < 1.2.

Ratio of psi(2S) to J/psi differential cross section times branching fractions assuming unpolarized productions as a function of pT for |y| < 1.2.

Lambda-Phi in the CS frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the CS frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the CS frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the CS frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the CS frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Phi in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta-Phi in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Tilde in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Phi in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta-Phi in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Tilde in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Phi in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta-Phi in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Tilde in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

The (unpolarised) acceptance corrected J/PSI cross section times branching ratio to MU+ MU- for (PROMPT) and NON-PROMPT) production in the |rapidity| bin 1.6-2.1.

The (unpolarised) acceptance corrected J/PSI cross section times branching ratio to MU+ MU- for (PROMPT) and NON-PROMPT) production in the |rapidity| bin 2.1-2.4.

The (unpolarised) acceptance corrected PSI(2S) cross section times branching ratio to MU+ MU- for (PROMPT) and NON-PROMPT production in the |rapidity| bin 0.0-1.2.

The (unpolarised) acceptance corrected PSI(2S) cross section times branching ratio to MU+ MU- for (PROMPT) and NON-PROMPT production in the |rapidity| bin 1.2-1.6.

The (unpolarised) acceptance corrected PSI(2S) cross section times branching ratio to MU+ MU- for (PROMPT) and NON-PROMPT production in the |rapidity| bin 1.6-2.4.

The uncorrected visible J/PSI cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 0.0-0.9.

The uncorrected visible J/PSI cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 0.9-1.2.

The uncorrected visible J/PSI cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 1.2-1.6.

The uncorrected visible J/PSI cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 1.6-2.1.

The uncorrected visible J/PSI cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 2.1-2.4.

The uncorrected visible PSI(2S) cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 0.0-1.2.

The uncorrected visible PSI(2S) cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 1.2-1.6.

The uncorrected visible PSI(2S) cross section times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 1.6-2.4.

The ratio between the PSI(2S) and J/PSI cross sections times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 0.0-1.2.

The ratio between the PSI(2S) and J/PSI cross sections times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 1.2-1.6.

The ratio between the PSI(2S) and J/PSI cross sections times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the |rapidity| bin 1.6-2.4.

The ratio between the PSI(2S) and J/PSI cross sections times branching ratio to MU+ MU- for PROMPT and NON-PROMPT production in the integrated |rapidity| bin 0.0-2.4.

The non-prompt fraction (FB) of J/PSI times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 0.0-0.9.

The non-prompt fraction (FB) of J/PSI times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 0.9-1.2.

The non-prompt fraction (FB) of J/PSI times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 1.2-1.6.

The non-prompt fraction (FB) of J/PSI times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 1.6-2.1.

The non-prompt fraction (FB) of J/PSI times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 2.1-2.4.

The non-prompt fraction (FB) of PSI(2S) times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 0.0-1.2.

The non-prompt fraction (FB) of PSI(2S) times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 1.2-1.6.

The non-prompt fraction (FB) of PSI(2S) times branching ratio to MU+ MU- from B-HADRON decays in the |rapidity| bin 1.6-2.4.

The difference between the J/PSI cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the J/PSI cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the J/PSI cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the J/PSI cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the J/PSI cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the PSI(2S) cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the PSI(2S) cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

The difference between the PSI(2S) cross section times branching ratio to MU+ MU- for non-prompt production calculated for four different polarization scenarios and the zero polarization values. The notation is (-1,CS), (+1,CS), (-1,HX), (+1,HX), where CS(HX) is the Collins-Soper (Helicity) frame and -1(+1) refer to full longitudinal (transverse) polarization.

Prompt J/$\psi$ $v_{2}$ as a function of rapidity.

Prompt J/$\psi$ $v_{2}$ as a function of $p_{T}$.

Prompt J/$\psi$ $v_{2}$ as a function of $p_{T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality at high $p_{\rm T}$, 6.5$<p_{\rm T}<$30 GeV/c, in three different $|y|$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality at high $p_{\rm T}$, 6.5$<p_{\rm T}<$30 GeV/c, in three different $|y|$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality, at forward rapidity, 1.6$<|y|<$2.4, for two different $p_{\rm T}$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality, at forward rapidity, 1.6$<|y|<$2.4, for two different $p_{\rm T}$ regions.

Nonprompt J/$\psi$ $v_{2}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $v_{2}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of pT, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

95% CL interval on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Nuclear modification factor RAA of psi(2S) integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Mean PT^2 values at negative, null and positive XF for D+AU reaction The mean PT is obtained with a phenomenological fit of the J/PSI distribution in PT of the form (1/(2*PI*PT))*D2(SIG)/DPT*DY = A ( 1+(PT/B)^2)^-6.

Mean PT^2 values at negative,null and positive XF for P+P reaction The mean PT is obtained with a phenomenological fit of the J/PSI distribution in PT of the form (1/(2*PI*PT))*D2(SIG)/DPT*DY = A ( 1+(PT/B)^2)^-6.

PT dependence of the Alpha power for Xf ~ -0.08.

PT dependance of the Alpha power for Xf ~ 0.

PT dependance of the Alpha power for Xf ~ 0.09.

Centrality dependance of the Nuclear modification factor for rapidity : y=-1.7. Centrality is represented by the number of collisions.

Centrality dependance of the Nuclear modification factor for rapidity y=0.Centrality is represented by the number of collisions.

Centrality dependance of the Nuclear modification factor for rapidity y=1.8.Centrality is represented by the number of 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.

Mean PT value for J/PSI production. 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 value given is a combined value of the mean PT for the electron channel results and the muon channel results. The systematic uncertainties were estimated from the spread in mean pt from a weighted mean of the binned data, the phenomenological fit and an exponential fit.An additionnal 3% incertainty is added in the systenatic error for the muon PT due to the incertainty in momentum scale.

Value of the total cross-section for J/PSI production. The total cross section is estimated using a pythia calculation, which reproduces the sigma distribution in rapidity the best, normalized to our data.The systematic error given is obtained by setting the measured cross sections to their upper and lower systematic error limits and noting how the total cross section changed.The choice of the PDF and model was estimated to have little impact in the value of the total cross section (~3%).In addition to the systematic error given, there is a +-40 microbarn error due to normalization.

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.