The production of $J/\psi$ mesons with rapidity $1.5<y<4.0$ or $5.0<y<2.5$ and transverse momentum $p_\mathrm{T}<14 \mathrm{GeV}/c$ is studied with the LHCb detector in protonlead collisions at a nucleonnucleon centreofmass energy $\sqrt{s_{NN}}=5 \mathrm{TeV}$. The analysis is based on a data sample corresponding to an integrated luminosity of about $1.6 \mathrm{nb}^{1}$. For the first time the nuclear modification factor and forwardbackward production ratio are determined separately for prompt $J/\psi$ mesons and $J/\psi$ from $b$hadron decays. Clear suppression of prompt $J/\psi$ production with respect to protonproton collisions at large rapidity is observed, while the production of $J/\psi$ from $b$hadron decays is less suppressed. These results show good agreement with available theoretical predictions. The measurement shows that cold nuclear matter effects are important for interpretations of the related quarkgluon plasma signatures in heavyion collisions.
Projections of the combined fit on (a, b) dimuon invariant mass and (c, d) $t_z$ in two representative bins in the (a, c) % forward and (b, d) % backward samples. For the mass projections the (red solid curve) total fitted function is shown together with the (blue dotted curve) $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ signal and (green dotted curve) background contributions. For the $t_z$ projections the total fitted function is indicated by the solid red curve, the background by the green hatched area, the prompt signal by the blue area and $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$ by the solid black curve. 
Fig1a.pdf [20 KiB] HiDef png [259 KiB] Thumbnail [91 KiB] *.C file 

Fig1b.pdf [20 KiB] HiDef png [267 KiB] Thumbnail [91 KiB] *.C file 

Fig1c.pdf [23 KiB] HiDef png [449 KiB] Thumbnail [103 KiB] *.C file 

Fig1d.pdf [23 KiB] HiDef png [472 KiB] Thumbnail [109 KiB] *.C file 

Distributions (normalised to unitary integral) of (a) track multiplicity and the $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ (b) transverse momentum $p_\mathrm{T}$, (c) momentum $p$, and (d) rapidity in laboratory frame $y_{\mathrm{lab}}$ in (black dots) forward and (red squares) backward regions of $ p\mathrm{Pb} $ collisions, and in (blue triangles) simulated $pp$ collisions. The distributions are background subtracted using the sPlot technique. 
Fig2a.pdf [15 KiB] HiDef png [143 KiB] Thumbnail [65 KiB] *.C file 

Fig2b.pdf [17 KiB] HiDef png [178 KiB] Thumbnail [70 KiB] *.C file 

Fig2c.pdf [17 KiB] HiDef png [194 KiB] Thumbnail [77 KiB] *.C file 

Fig2d.pdf [17 KiB] HiDef png [166 KiB] Thumbnail [65 KiB] *.C file 

Single differential production crosssections for (black dots) prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ and (red squares) $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$ as functions of (a, b) $p_\mathrm{T}$ and (c, d) $y$ in the (a, c) forward and (b, d) backward regions. 
Fig3a.pdf [14 KiB] HiDef png [122 KiB] Thumbnail [58 KiB] *.C file 

Fig3b.pdf [15 KiB] HiDef png [129 KiB] Thumbnail [59 KiB] *.C file 

Fig3c.pdf [14 KiB] HiDef png [136 KiB] Thumbnail [63 KiB] *.C file 

Fig3d.pdf [14 KiB] HiDef png [138 KiB] Thumbnail [65 KiB] *.C file 

Double differential production crosssections for (a) prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ mesons and (b) $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$ in the forward samples. 
Fig4a.pdf [16 KiB] HiDef png [173 KiB] Thumbnail [71 KiB] *.C file 

Fig4b.pdf [17 KiB] HiDef png [185 KiB] Thumbnail [73 KiB] *.C file 

Nuclear modification factor $R_{ p\mathrm{Pb} }$ as a function of $y$ for (a) prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ mesons and (b) $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$, together with the theoretical predictions from (yellow dashed line and brown band) Refs. \cite{Ferreiro:2013pua,delValle:2014wha}, (blue band) Ref. \cite{Albacete:2013ei}, and (green solid and blue dashdotted lines) Ref. \cite{Arleo:2012rs}. The inner error bars (delimited by the horizontal lines) show the statistical uncertainties; the outer ones show the statistical and systematic uncertainties added in quadrature. The uncertainty due to the interpolated $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ crosssection in $pp$ collisions at $\sqrt{s}=5\mathrm{ Te V} $ is $5.5\%$ ($8.4\%$) for prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ mesosns ($ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$). 
Fig5a.pdf [16 KiB] HiDef png [699 KiB] Thumbnail [111 KiB] *.C file 

Fig5b.pdf [14 KiB] HiDef png [461 KiB] Thumbnail [84 KiB] *.C file 

Forwardbackward production ratio $R_{ {FB}}$ as a function of $y$ for (a) prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ mesons and (b) $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$, together with the theoretical predictions from (yellow dashed line and brown band) Refs. \cite{Ferreiro:2013pua,delValle:2014wha}, (blue band) Ref. \cite{Albacete:2013ei}, and (green solid and blue dashdotted lines) Ref. \cite{Arleo:2012rs}. The inner error bars (delimited by the horizontal lines) show the statistical uncertainties; the outer ones show the statistical and systematic uncertainties added in quadrature. 
Fig6a.pdf [15 KiB] HiDef png [843 KiB] Thumbnail [119 KiB] *.C file 

Fig6b.pdf [14 KiB] HiDef png [737 KiB] Thumbnail [105 KiB] *.C file 

Forwardbackward production ratio $R_{ {FB}}$ for (black dot) prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ mesons and (red square) $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ from $b$ as a function of $p_\mathrm{T}$ in the rapidity range $2.5<y<4.0$. The theoretical predictions from (blue band) Ref. \cite{Albacete:2013ei} and (green solid and blue dashdotted lines) Ref. \cite{Arleo:2013zua} are for prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} $ mesons. The inner error bars (delimited by the horizontal lines) show the statistical uncertainties; the outer ones show the statistical and systematic uncertainties added in quadrature. 
Fig7.pdf [17 KiB] HiDef png [239 KiB] Thumbnail [88 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2013052.gif Thumbnail 
Relative systematic uncertainties on the differential production crosssections. The uncertainty due to the radiative tail and branching fraction cancels in both $R_{ p\mathrm{Pb} }$ and $R_{ {FB}}$. The uncertainty due to the tracking efficiency and the luminosity partially cancels for $R_{ {FB}}$. 
Table_1.pdf [67 KiB] HiDef png [113 KiB] Thumbnail [17 KiB] tex code 

Single differential production crosssections (in $ {\rm \upmu b} /( {\mathrm{ Ge V /}c} )$) for prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ mesons and $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ from $b$ as functions of transverse momentum. The first uncertainty is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins, and the third is the correlated component. 
Table_2.pdf [39 KiB] HiDef png [233 KiB] Thumbnail [44 KiB] tex code 

Single differential production crosssections (in $ {\rm \upmu b}$ ) for prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ mesons and $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ from $b$ as functions of rapidity. The first uncertainty is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins, and the third is the correlated component. 
Table_3.pdf [39 KiB] HiDef png [170 KiB] Thumbnail [31 KiB] tex code 

Nuclear modification factor $R_{ p\mathrm{Pb} }$ as a function of $y$ with $ p_{\rm T} <14 {\mathrm{ Ge V /}c} $. The first uncertainty is statistical, the second is the systematic, and the third is the uncertainty related to the interpolation error. 
Table_4.pdf [45 KiB] HiDef png [56 KiB] Thumbnail [9 KiB] tex code 

Forwardbackward production ratio $R_{ {FB}}$ as a function of $y$ with $ p_{\rm T} <14 {\mathrm{ Ge V /}c} $. The first uncertainty is statistical, the second is the uncorrelated systematic component, and the third is the systematic uncertainty correlated between bins. 
Table_5.pdf [38 KiB] HiDef png [49 KiB] Thumbnail [7 KiB] tex code 

Forwardbackward production ratio $R_{ {FB}}$ as a function of $p_\mathrm{T}$ with $2.5<y<4.0$. The first uncertainty is statistical, the second is the uncorrelated systematic component, and the third is the systematic uncertainty correlated between bins. 
Table_6.pdf [46 KiB] HiDef png [122 KiB] Thumbnail [18 KiB] tex code 
Created on 14 June 2019.Citation count from INSPIRE on 14 June 2019.