We report on the inclusive production cross sections of J/$\psi$, $\psi$(2S), $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S), measured at forward rapidity with the ALICE detector in pp collisions at a center-of-mass energy $\sqrt{s}=8$ TeV. The analysis is based on data collected at the LHC and corresponds to an integrated luminosity of 1.28 pb$^{-1}$. Quarkonia are reconstructed in the dimuon-decay channel. The differential production cross sections are measured as a function of the transverse momentum $p_{\rm T}$ and rapidity $y$, over the $p_{\rm T}$ ranges $0<p_{\rm T}<20$ GeV/$c$ for J/$\psi$, $0<p_{\rm T}<12$ GeV/$c$ for all other resonances, and for $2.5<y<4$. The cross sections, integrated over $p_{\rm T}$ and $y$, and assuming unpolarized quarkonia, are $\sigma_{{\rm J}/\psi} = 8.98\pm0.04\pm0.82$ $\mu$b, $\sigma_{\psi{\rm (2S)}} = 1.23\pm0.08\pm0.22$ $\mu$b, $\sigma_{\Upsilon{\rm(1S)}} = 71\pm6\pm7$ nb, $\sigma_{\Upsilon{\rm(2S)}} = 26\pm5\pm4$ nb and $\sigma_{\Upsilon{\rm(3S)}} = 9\pm4\pm1$ nb, where the first uncertainty is statistical and the second one is systematic. These values agree, within at most $1.4\sigma$, with measurements performed by the LHCb collaboration in the same rapidity range.

The pseudorapidity ($\eta$) and transverse-momentum ($p_{\rm T}$) distributions of charged particles produced in proton-proton collisions are measured at the centre-of-mass energy $\sqrt{s}$ = 13 TeV. The pseudorapidity distribution in $|\eta|<$ 1.8 is reported for inelastic events and for events with at least one charged particle in $|\eta|<$ 1. The pseudorapidity density of charged particles produced in the pseudorapidity region $|\eta|<$ 0.5 is 5.31 $\pm$ 0.18 and 6.46 $\pm$ 0.19 for the two event classes, respectively. The transverse-momentum distribution of charged particles is measured in the range 0.15 $<$ $p_{\rm T}$ $<$ 20 GeV/c and $|\eta|<$ 0.8 for events with at least one charged particle in $|\eta|<$ 1. The correlation between transverse momentum and particle multiplicity is also investigated by studying the evolution of the spectra with event multiplicity. The results are compared with calculations from PYTHIA and EPOS Monte Carlo generators.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonances has been measured in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV using the ALICE detector. Resonances are reconstructed via their hadronic decay channels in the rapidity interval $-$0.5 $<$$y$$<$ 0 and the transverse momentum spectra are measured for various multiplicity classes up to $p_{\rm T}$ = 20 GeV/$c$ for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $p_{\rm T}$ = 16 GeV/$c$ for $\mathrm{\phi(1020)}$. The $p_{\rm T}$ -integrated yields and mean transverse momenta are reported and compared with previous results in pp, p-Pb and Pb-Pb collisions. The $x_{\mathrm{T}}$ scaling for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonance production is newly tested in p-Pb collisions and found to hold in the high-$p_{\rm T}$ region at LHC energies. The nuclear modification factors ($R_{\rm pPb}$) as a function of $p_{\rm T}$ for $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ at $\sqrt{s_{NN}}$ = 8.16 TeV are presented along with the new $R_{\rm pPb}$ measurements of $\mathrm{K}^{*0}$, $\mathrm{\phi}$ , $\Xi$, and $\Omega$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. At intermediate $p_{\rm T}$ (2-8 GeV/$c$), $R_{\rm pPb}$ of $\Xi$, $\Omega$ show a Cronin-like enhancement, while $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ show no or little nuclear modification. At high $p_{\rm T}$ ($>$ 8 GeV/$c$), the $R_{\rm pPb}$ values of all hadrons are consistent with unity within uncertainties. The $R_{\rm pPb}$ of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ at $\sqrt{s_{\rm NN}}$ = 8.16 and 5.02 TeV show no significant energy dependence.

The study of the azimuthal anisotropy of inclusive muons produced in p-Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV, using the ALICE detector at the LHC is reported. The measurement of the second-order Fourier coefficient of the particle azimuthal distribution, $v_2$, is performed as a function of transverse momentum $p_{\rm T}$ in the 0-20% high-multiplicity interval at both forward ($2.03 < y_{\rm CMS} < 3.53$) and backward ($-4.46 < y_{\rm CMS} < -2.96$) rapidities over a wide $p_{\rm T}$ range, $0.5 < p_{\rm T} < 10$ GeV/$c$, in which a dominant contribution of muons from heavy-flavour hadron decays is expected at $p_{\rm T} > 2$ GeV/$c$. The $v_2$ coefficient of inclusive muons is extracted using two different techniques, namely two-particle cumulants, used for the first time for heavy-flavour measurements, and forward-central two-particle correlations. Both techniques give compatible results. A positive $v_2$ is measured at both forward and backward rapidities with a significance larger than $4.7\sigma$ and $7.6\sigma$, respectively, in the interval $2 < p_{\rm T} < 6$ GeV/$c$. Comparisons with previous measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, and with AMPT and CGC-based theoretical calculations are discussed. The findings impose new constraints on the theoretical interpretations of the origin of the collective behaviour in small collision systems.

Measurements of the production of electrons from heavy-flavour hadron decays in pp collisions at $\sqrt{s} = 13$ TeV at midrapidity with the ALICE detector are presented down to a transverse momentum ($p_{\rm T}$) of 0.2 GeV$/c$ and up to $p_{\rm T} = 35$ GeV$/c$, which is the largest momentum range probed for inclusive electron measurements in ALICE. In p$-$Pb collisions, the production cross section and the nuclear modification factor of electrons from heavy-flavour hadron decays are measured in the $p_{\rm T}$ range $0.5 < p_{\rm T} < 26$ GeV$/c$ at $\sqrt{s_{\rm NN}} = 8.16$ TeV. The nuclear modification factor is found to be consistent with unity within the statistical and systematic uncertainties. In both collision systems, first measurements of the yields of electrons from heavy-flavour hadron decays in different multiplicity intervals normalised to the multiplicity-integrated yield (self-normalised yield) at midrapidity are reported as a function of the self-normalised charged-particle multiplicity estimated at midrapidity. The self-normalised yields in pp and p$-$Pb collisions grow faster than linear with the self-normalised multiplicity. A strong $p_{\rm T}$ dependence is observed in pp collisions, where the yield of high-$p_{\rm T}$ electrons increases faster as a function of multiplicity than the one of low-$p_{\rm T}$ electrons. The measurement in p$-$Pb collisions shows no $p_{\rm T}$ dependence within uncertainties. The self-normalised yields in pp and p$-$Pb collisions are compared with measurements of other heavy-flavour, light-flavour, and strange particles, and with Monte Carlo simulations.

The production of $\Upsilon$ mesons in $pp$ collisions at $\sqrt=7$ and $8\,\mathrm{TeV}$ is studied with the LHCb detector using data samples corresponding to an integrated luminosity of $1\,\mathrm{fb}^{-1}$ and $2\,\mathrm{fb}^{-1}$ respectively. The production cross-sections and ratios of cross-sections are measured as functions of the meson transverse momentum $p_T$ and rapidity $y$, for $p_T<30\,\mathrm{GeV}/c$} and $2.0<y<4.5$.

The polarization of the $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S) $mesons, produced in $pp$ collisions at centre-of-mass energies $\sqrt{s}$=7 and 8TeV, is measured using data samples collected by the LHCb experiment, corresponding to integrated luminosities of 1 and 2fb$^{-1}$, respectively. The measurements are performed in three polarization frames, using $\Upsilon\to\mu^+\mu^-$ decays in the kinematic region of the transverse momentum $p_{T}(\Upsilon)<30GeV/c$, and rapidity $2.2<y(\Upsilon)<4.5$. No large polarization is observed.

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.

The production of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV is studied with the \lhcb detector. Cross-section measurements are performed as a function of the transverse momentum $p_\mathrm{T}$ and the rapidity $y$ of the $J/\psi$ meson in the region $p_\mathrm{T}<14\mathrm{GeV}/c$ and $2.0<y<4.5$, for both prompt $J/\psi$ mesons and $J/\psi$ mesons from $b$-hadron decays. The production cross-sections integrated over the kinematic coverage are $15.03\pm 0.03\pm 0.94$ $\mu$b for prompt $J/\psi$ and $2.25\pm 0.01\pm 0.14$ $\mu$b for $J/\psi$ from $b$-hadron decays, assuming zero polarization of the $J/\psi$ meson. The first uncertainties are statistical and the second systematic. The cross-section reported for $J/\psi$ mesons from $b$-hadron decays is used to extrapolate to a total $b\bar{b}$ cross-section. The ratios of the cross-sections with respect to $\sqrt{s}=8$ TeV are also determined.

The inclusive $D_s^{\pm}$ production asymmetry is measured in $pp$ collisions collected by the LHCb experiment at centre-of-mass energies of $\sqrt{s} =7$ and 8 TeV. Promptly produced $D_s^{\pm}$ mesons are used, which decay as $D_s^{\pm}\to\phi\pi^{\pm}$, with $\phi\to K^+K^-$. The measurement is performed in bins of transverse momentum, $p_{\rm T}$, and rapidity, $y$, covering the range $2.5