The production of $\Upsilon (nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{-1}$. The $\Upsilon (nS)$ mesons are reconstructed through their decays into two opposite-sign muons. The measurements comprise the differential production cross-sections of the $\Upsilon (1S)$ and $\Upsilon (2S)$ states, their forward-to-backward ratios and nuclear modification factors, performed as a function of the transverse momentum $p_{\mathrm{T}}$ and rapidity in the nucleon-nucleon centre-of-mass frame $y^*$ of the $\Upsilon (nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($-5.0<y^*<-2.5$) for $p$Pb (Pb$p$) collisions. In addition, production cross-sections for $\Upsilon (3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon (nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.

Production cross-sections of prompt charm mesons are measured with the first data from $pp$ collisions at the LHC at a centre-of-mass energy of $13\,\mathrm{TeV}$. The data sample corresponds to an integrated luminosity of $4.98 \pm 0.19\,\mathrm{pb}^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^{0}$, $D^{+}$, $D_{s}^{+}$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\mathrm{T}}$, and rapidity, $y$, and cover the range $0 < p_{\mathrm{T}} < 15\,\mathrm{GeV}/c$ and $2.0 < y < 4.5$. The inclusive cross-sections for the four mesons, including charge conjugation, within the range of $1 < p_{\mathrm{T}} < 8\,\mathrm{GeV}/c$ are found to be \begin{equation} \sigma(pp \to D^{0} X) = 2072 \pm 2 \pm 124\,\mu\mathrm{b}\\ \sigma(pp \to D^{+} X) = 834 \pm 2 \pm \phantom{1}78\,\mu\mathrm{b}\\ \sigma(pp \to D_{s}^{+} X) = 353 \pm 9 \pm \phantom{1}76\,\mu\mathrm{b}\\ \sigma(pp \to D^{*+} X) = 784 \pm 4 \pm \phantom{1}87\,\mu\mathrm{b} \end{equation} where the uncertainties are due to statistical and systematic uncertainties, respectively.

Measurements are presented of electroweak boson production using data from $pp$ collisions at a centre-of-mass energy of $\sqrt{s} = 8\mathrm{\,Te\kern -0.1em V}$. The analysis is based on an integrated luminosity of $2.0\mathrm{\,fb}^{-1}$ recorded with the LHCb detector. The bosons are identified in the $W\rightarrow\mu\nu$ and $Z\rightarrow\mu^{+}\mu^{-}$ decay channels. The cross-sections are measured for muons in the pseudorapidity range $2.0 < \eta < 4.5$, with transverse momenta $p_{\rm T} > 20{\mathrm{\,Ge\kern -0.1em V\!/}c}$ and, in the case of the $Z$ boson, a dimuon mass within $60 < M_{\mu^{+}\mu^{-}} < 120{\mathrm{\,Ge\kern -0.1em V\!/}c^{2}}$. The results are \begin{align*} \sigma_{W^{+}\rightarrow\mu^{+}\nu} &= 1093.6 \pm 2.1 \pm 7.2 \pm 10.9 \pm 12.7{\rm \,pb} \, , \sigma_{W^{-}\rightarrow\mu^{-}\bar{\nu}} &= \phantom{0}818.4 \pm 1.9 \pm 5.0 \pm \phantom{0}7.0 \pm \phantom{0}9.5{\rm \,pb} \, , \sigma_{Z\rightarrow\mu^{+}\mu^{-}} &= \phantom{00}95.0 \pm 0.3 \pm 0.7 \pm \phantom{0}1.1 \pm \phantom{0}1.1{\rm \,pb} \, , \end{align*} where the first uncertainties are statistical, the second are systematic, the third are due to the knowledge of the LHC beam energy and the fourth are due to the luminosity determination. The evolution of the $W$ and $Z$ boson cross-sections with centre-of-mass energy is studied using previously reported measurements with $1.0\mathrm{\,fb}^{-1}$ of data at $7\mathrm{\,Te\kern -0.1em V}$. Differential distributions are also presented. Results are in good agreement with theoretical predictions at next-to-next-to-leading order in perturbative quantum chromodynamics.

The production of the $\eta_c (1S)$ state in proton-proton collisions is probed via its decay to the $p \bar{p}$ final state with the LHCb detector, in the rapidity range $2.0 < y < 4.5$ and in the meson transverse-momentum range $p_T > 6.5$ GeV/c. The cross-section for prompt production of $\eta_c (1S)$ mesons relative to the prompt $J/\psi$ cross-section is measured, for the first time, to be $\sigma_{\eta_c (1S)}/\sigma_{J/\psi} = 1.74 \pm 0.29 \pm 0.28 \pm 0.18 _{B}$ at a centre-of-mass energy $\sqrt{s} = 7$ TeV using data corresponding to an integrated luminosity of 0.7 fb$^{-1}$, and $\sigma_{\eta_c (1S)}/\sigma_{J/\psi} = 1.60 \pm 0.29 \pm 0.25 \pm 0.17 _{B}$ at $\sqrt{s} = 8$ TeV using 2.0 fb$^{-1}$. The uncertainties quoted are, in order, statistical, systematic, and that on the ratio of branching fractions of the $\eta_c (1S)$ and $J/\psi$ decays to the $p \bar{p}$ final state. In addition, the inclusive branching fraction of $b$-hadron decays into $\eta_c (1S)$ mesons is measured, for the first time, to be $B ( b \rightarrow \eta_c X ) = (4.88 \pm 0.64 \pm 0.29 \pm 0.67 _{B}) \times 10^{-3}$, where the third uncertainty includes also the uncertainty on the $J/\psi$ inclusive branching fraction from $b$-hadron decays. The difference between the $J/\psi$ and $\eta_c (1S)$ meson masses is determined to be $114.7 \pm 1.5 \pm 0.1$ MeV/c$^2$.

Associated production of bottomonia and open charm hadrons in $pp$ collisions at $\sqrt{s}=7$ and $8$TeV is observed using data corresponding to an integrated luminosity of 3$fb^{-1}$ accumulated with the LHCb detector. The observation of five combinations, $\Upsilon(1S)D^0$, $\Upsilon(2S)D^0$, $\Upsilon(1S)D^+$, $\Upsilon(2S)D^+$ and $\Upsilon(1S)D^+_{s}$, is reported. Production cross-sections are measured for $\Upsilon(1S)D^0$ and $\Upsilon(1S)D^+$ pairs in the forward region. The measured cross-sections and the differential distributions indicate the dominance of double parton scattering as the main production mechanism. This allows a precise measurement of the effective cross-section for double parton scattering.

Using a data sample corresponding to an integrated luminosity of $2.0\,fb^{-1}$, collected by the LHCb experiment, the production of the $\eta_c(1S)$ state in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13 \text{ TeV}$ is studied in the rapidity range ${2.0 < y < 4.5}$ and in the transverse momentum range ${6.5 < p_{T} < 14.0\text{ GeV}}$. The cross-section for prompt production of $\eta_c(1S)$ mesons relative to that of the $J/\psi$ meson is measured using the ${p\bar{p}}$ decay mode and is found to be ${\sigma_{\eta_c(1S)}/\sigma_{J/\psi} = 1.69 \pm 0.15 \pm 0.10 \pm 0.18}$. The quoted uncertainties are, in order, statistical, systematic and due to uncertainties on the branching fractions of the ${J/\psi\to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The prompt $\eta_c(1S)$ production cross-section is determined to be ${\sigma_{\eta_c(1S)} = 1.26 \pm 0.11\pm 0.08 \pm 0.14 \,\mu b}$, where the last uncertainty includes that on the ${J/\psi}$ meson cross-section. The ratio of the branching fractions of $b$-hadron decays to the $\eta_c(1S)$ and ${J/\psi}$ states is measured to be ${\mathcal{B}_{b\to\eta_c X}/\mathcal{B}_{b\to J/\psi X} = 0.48 \pm 0.03 \pm 0.03 \pm 0.05}$, where the last uncertainty is due to those on the branching fractions of the ${J/\psi \to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The difference between the ${J/\psi}$ and $\eta_c(1S)$ masses is also determined to be ${113.0 \pm 0.7 \pm 0.1\text{ MeV}}$, which is the most precise single measurement of this quantity to date.

The production of prompt charged particles in proton-lead collisions and in proton-proton collisions at the nucleon-nucleon centre-of-mass energy ${\sqrt{s_{\scriptscriptstyle\mathrm{NN}}}=5\,\mathrm{TeV}}$ is studied at LHCb as a function of pseudorapidity ($\eta$) and transverse momentum ($p_{\mathrm{T}}$) with respect to the proton beam direction. The nuclear modification factor for charged particles is determined as a function of $\eta$ between ${-4.8<\eta<-2.5}$ (backward region) and ${2.0<\eta<4.8}$ (forward region), and $p_{\mathrm{T}}$ between ${0.2<p_{\mathrm{T}}<8.0\,\mathrm{GeV}/c}$. The results show a suppression of charged particle production in proton-lead collisions relative to proton-proton collisions in the forward region and an enhancement in the backward region for $p_{\mathrm{T}}$ larger than $1.5\,\mathrm{GeV}/c$. This measurement constrains nuclear PDFs and saturation models at previously unexplored values of the parton momentum fraction down to $10^{-6}$.

The production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons in proton-proton collisions at the centre-of-mass energy of sqrt(s)=7 TeV is studied with the LHCb detector. The analysis is based on a data sample of 25 pb-1 collected at the Large Hadron Collider. The Upsilon mesons are reconstructed in the decay mode Upsilon -&gt; mu+ mu- and the signal yields are extracted from a fit to the mu+ mu- invariant mass distributions. The differential production cross-sections times dimuon branching fractions are measured as a function of the Upsilon transverse momentum pT and rapidity y, over the range pT &lt; 15 GeV/c and 2.0 &lt; y &lt; 4.5. The cross-sections times branching fractions, integrated over these kinematic ranges, are measured to be sigma(pp -&gt; Upsilon(1S) X) x B(Upsilon(1S)-&gt;mu+ mu-) = 2.29 {\pm} 0.01 {\pm} 0.10 -0.37 +0.19 nb, sigma(pp -&gt; Upsilon(2S) X) x B(Upsilon(2S)-&gt;mu+ mu-) = 0.562 {\pm} 0.007 {\pm} 0.023 -0.092 +0.048 nb, sigma(pp -&gt; Upsilon(3S) X) x B(Upsilon(3S)-&gt;mu+ mu-) = 0.283 {\pm} 0.005 {\pm} 0.012 -0.048 +0.025 nb, where the first uncertainty is statistical, the second systematic and the third is due to the unknown polarisation of the three Upsilon states.

Charm production at the LHC in pp collisions at sqrt(s)=7 TeV is studied with the LHCb detector. The decays D0 -> K- pi+, D+ -> K- pi+ pi+, D*+ -> D0(K- pi+) pi+, D_s+ -> phi(K- K+) pi+, Lambda_c+ -> p K- pi+, and their charge conjugates are analysed in a data set corresponding to an integrated luminosity of 15 nb^{-1}. Differential cross-sections dsigma/dp_T are measured for prompt production of the five charmed hadron species in bins of transverse momentum and rapidity in the region 0 < p_T < 8 GeV/c and 2.0 < y < 4.5. Theoretical predictions are compared to the measured differential cross-sections. The integrated cross-sections of the charm hadrons are computed in the above p_T-y range, and their ratios are reported. A combination of the five integrated cross-section measurements gives sigma(c\bar{c})_{p_T < 8 GeV/c, 2.0 < y < 4.5} = 1419 +/- 12 (stat) +/- 116 (syst) +/- 65 (frag) microbarn, where the uncertainties are statistical, systematic, and due to the fragmentation functions.

We have measured antiproton production cross sections as functions of centrality in collisions of 14.6 GeV/c per nucleon Si28 ions with targets of Al, Cu, and Pb. For all targets, the antiproton yields increase linearly with the number of projectile nucleons that have interacted, and show little target dependence. We discuss the implications of this result on the production and absorption of antiprotons within the nuclear medium.