$\Lambda_{\rm c}^+$ production in pp collisions at $\sqrt{s} = 7$ TeV and in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV

The $p_{\rm T}$-differential production cross section of prompt $\Lambda_{\rm c}^+$ charmed baryons was measured with the ALICE detector at the Large Hadron Collider (LHC) in pp collisions at $\sqrt{s} = 7$ TeV and in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV at midrapidity. The $\Lambda_{\rm c}^+$ and ${\overline{\Lambda}}_{\rm c}^-$ were reconstructed in the hadronic decay modes $\Lambda_{\rm c}^+\rightarrow {\rm p}{\rm K^-}\pi^+$, $\Lambda_{\rm c}^+\rightarrow {\rm p}{\rm K_{\rm S}^0}$ and in the semileptonic channel $\Lambda_{\rm c}^+\rightarrow {\rm e^+}\nu_{\rm e}\Lambda$ (and charge conjugates). The measured values of the $\Lambda_{\rm c}^+/{\rm D_0}$ ratio, which is sensitive to the c-quark hadronisation mechanism, and in particular to the production of baryons, are presented and are larger than those measured previously in different colliding systems, centre-of-mass energies, rapidity and $p_{\rm T}$ intervals, where the $\Lambda_{\rm c}^+$ production process may differ. The results are compared with the expectations obtained from perturbative Quantum Chromodynamics calculations and Monte Carlo event generators. Neither perturbative QCD calculations nor Monte Carlo models reproduce the data, indicating that the fragmentation of heavy-flavour baryons is not well understood. The first measurement at the LHC of the $\Lambda_{\rm c}^+$ nuclear modification factor, $R_{\rm pPb}$, is also presented. The $R_{\rm pPb}$ is found to be consistent with unity and with that of D mesons within the uncertainties, and consistent with a theoretical calculation that includes cold nuclear matter effects and a calculation that includes charm quark interactions with a deconfined medium.

Λ + c production in pp collisions at √ s = 7 TeV and in p-Pb collisions at √ s NN = 5.02 TeV ALICE Collaboration *

Abstract
The p T -differential production cross section of prompt Λ + c charmed baryons was measured with the ALICE detector at the Large Hadron Collider (LHC) in pp collisions at √ s = 7 TeV and in p-Pb collisions at √ s NN = 5.02 TeV at midrapidity. The Λ + c and Λ − c were reconstructed in the hadronic decay modes Λ + c → pK − π + , Λ + c → pK 0 S and in the semileptonic channel Λ + c → e + ν e Λ (and charge conjugates). The measured values of the Λ + c /D 0 ratio, which is sensitive to the c-quark hadronisation mechanism, and in particular to the production of baryons, are presented and are larger than those measured previously in different colliding systems, centre-of-mass energies, rapidity and p T intervals, where the Λ + c production process may differ. The results are compared with the expectations obtained from perturbative Quantum Chromodynamics calculations and Monte Carlo event generators. Neither perturbative QCD calculations nor Monte Carlo models reproduce the data, indicating that the fragmentation of heavy-flavour baryons is not well understood. The first measurement at the LHC of the Λ + c nuclear modification factor, R pPb , is also presented. The R pPb is found to be consistent with unity and with that of D mesons within the uncertainties, and consistent with a theoretical calculation that includes cold nuclear matter effects and a calculation that includes charm quark interactions with a deconfined medium.

Introduction
The study of charm production at the Large Hadron Collider (LHC) is an important tool to test predictions obtained from perturbative Quantum Chromodynamics (pQCD) calculations for proton-proton (pp) collisions. These calculations are based on the factorisation approach that describes heavy-flavour production as a convolution of the parton distribution functions, the parton hard-scattering cross section and the fragmentation function. The cross section for heavy-flavour hadron production can be obtained from perturbative calculations at next-to-leading order with next-to-leading-log resummation, like the General-Mass Variable-Flavour-Number Scheme (GM-VFNS [1,2]) and Fixed-Order Next-to-Leading-Log (FONLL [3,4]) approaches. No predictions are, however, available for baryons in the latter approach due to lack of knowledge of the fragmentation function of charm quarks into baryons. Cross section calculations are available also with the k T factorisation framework [5]. These theoretical calculations generally describe within uncertainties the measurements at the LHC, with the central predictions for beauty production lying closer to data than the central predictions for charm production [6]. The measured transverse momentum differential cross section of charm mesons lies in the upper part of the FONLL uncertainty band and is systematically below the central value of GM-VFNS predictions [7]. Cross sections for charm production are also available in general-purpose Monte Carlo generators such as PYTHIA [8]. The hard process amplitude is calculated with leading order (LO) accuracy and, via parton showers, effective LO+LL accuracy is provided. Next-to-Leading-Order (NLO) Monte Carlo generators were developed by matching event generators, calculating the hard scattering with NLO accuracy, as in POWHEG [9], with parton showers as in PYTHIA. In pQCD calculations, the hadronisation process is modeled via a fragmentation function, which parametrises the fraction of the quark energy transferred to the produced hadron, and by the fragmentation fractions, which account for the probability of a heavy quark to hadronise into a particular hadron species. Fragmentation functions are tuned on electron-positron data under the assumption that they are universal. Similarly, the fragmentation fractions were usually assumed to be the same in different collision systems. Among other observables, the relative production of baryons and mesons ("baryon-to-meson ratio") is particularly sensitive to the fragmentation process. A study of the Λ 0 b baryon to B − and B 0 meson production by LHCb [10] reported a transverse momentum (p T ) dependence of that ratio, interpreted as evidence of non-universality of fragmentation fractions in the beauty sector [11,12]. In Monte Carlo generators, hadronisation is implemented via formation of strings as in PYTHIA, via ropes [13] as in DIPSY [14] or via clusters as in HERWIG [15]. In hadron-hadron collisions at LHC energies, multi-parton interactions and coherence effects between multiple partonic interactions may affect the hadronisation processes. Within the existing PYTHIA8 framework a better agreement with measurements by CMS [16] of the Λ/K 0 S ratio was obtained in [17] introducing additional colour reconnection mechanisms that play a role in pp collisions and are instead expected to be highly suppressed in electron-positron collisions at LEP. For the DIPSY event generator in [18] an approach was tested where strings from independent interactions can be close in space and form colour ropes, expected to yield more baryons than independent strings. Therefore, the measurement of the Λ + c production cross section in pp collisions allows one to test these expectations at LHC energies with charmed baryons and mesons.
Furthermore, the study of charmed-baryon production could play an important role in the investigation of the state of strongly-interacting matter at very high temperatures and densities realised in heavy-ion collisions, known as the Quark-Gluon Plasma (QGP) [19]. Measurements of open heavy-flavour production in this environment allow for the study of the interaction of heavy quarks with the medium constituents and the characterisation of the properties of the plasma state [20]. The interaction with the medium constituents could modify the hadronisation: a significant fraction of low and intermediate-momentum charm and beauty quarks could hadronise via recombination (coalescence) with other quarks from the medium [21,22]. Models including coalescence predict an enhanced baryon-to-meson ratio at low and intermediate p T relative to that observed in pp collisions where hadronisation can be described by string-Λ + c production in ALICE ALICE Collaboration fragmentation models [8]. In addition, the possible existence of light diquark bound states in the QGP could further enhance the Λ + c /D 0 ratio in coalescence models [23]. An enhancement of the p T -integrated Λ + c /D 0 ratio in presence of a QGP is also predicted by statistical hadronisation models [24], where the relative abundance of hadrons depends only on their masses and on the freeze-out temperature of the medium created in the collision. Recently, such an enhancement of the Λ + c /D 0 ratio was preliminarily reported by STAR in Au-Au collisions at √ s NN = 200 GeV in the 3 < p T < 6 GeV/c interval [25]. A measurement of prompt Λ + c production at the LHC in pp collisions is needed as a baseline reference for these studies.
For the intepretation of the results in nucleus-nucleus collisions, the measurement in proton-nucleus collisions is also crucial. In such a system cold-nuclear-matter (CNM) effects can affect the production of charm hadrons: their assessment is needed to disentangle them from the effects related to the formation of the QGP (hot-medium effects). In the initial state, the Parton Distribution Functions (PDFs) are modified in bound nucleons compared to free nucleons. The nuclear shadowing at low transverse momentum can decrease, among other effects, the production cross section of open charm [26]. Moreover, the multiple scattering of partons in the nucleus before or after the hard scattering can affect the momentum distributions of the produced hadrons, especially at low p T (p T < 2 GeV/c). In addition to initial-state effects, final-state effects may also be responsible for the modification of particle yields and transverse-momentum distributions in proton-nucleus collisions as compared to pp interactions. Nuclear effects can be investigated measuring the nuclear modification factor R pPb , defined as the ratio of the cross section in p-Pb collisions to that in pp interactions scaled by the mass number of the Pb nucleus. A recent measurement [27,28] of D-meson production in p-Pb collisions showed that, within uncertainties, R pPb is compatible with unity, indicating that initial and final-state effects are either small or that they compensate each other. Several other observations in p-Pb collisions, such as the presence of di-hadron azimuthal correlations at large rapidity differences [29][30][31][32][33], the evolution of the average p T at central rapidity of identified hadrons with multiplicity [34,35] and the increased strangeness yield with increasing multiplicity [36] qualitatively resemble observations in Pb-Pb collisions. This suggests the possible formation of a hot deconfined medium also in p-Pb collisions that, in turn, can affect the propagation and hadronisation of heavy quarks, modifying the momentum distribution of the charmed hadrons with respect to that expected from pp collisions, hence inducing a deviation of R pPb from unity [37, 38].
At high energies, Λ + c production has been studied at electron-positron colliders (at the Z-resonance with LEP [39][40][41], and at B factories [42-45]), in several fixed target experiments including neutrinoproton [46], hadron-nucleon [47] and photon-nucleon [48] interactions and at electron-proton colliders (in photoproduction [49,50], and via deep inelastic scattering [51]). At the LHC, a measurement of Λ + c -baryon production at forward rapidity was reported by the LHCb Collaboration [52] in pp collisions at √ s = 7 TeV in the rapidity (y) range 2.0 < y < 4.5. Here and in the following, y is defined in the centre-of-mass system of the collision. A preliminary result in p-Pb collisions at √ s NN = 5.02 TeV has also been presented recently by LHCb [53]. Previous measurements at hadron-hadron colliders [54][55][56] are at a much lower centre-of-mass energy ( √ s = O(100) GeV).
In this paper, we present the measurement of the production cross section of the prompt charmed baryon Λ + c (udc) and its charge conjugate (c.c.). Hereafter with Λ c we will refer indistinctly to both, and all mentioned decay channels refer also to their charge conjugate. The contribution from beauty feed-down to the measured Λ c yields was subtracted by using pQCD calculations of the beauty-hadron cross section together with the acceptance and efficiency values extracted from simulation. The cross section was measured with the ALICE detector [57] in pp collisions at √ s = 7 TeV in the transverse momentum and rapidity intervals 1 < p T < 8 GeV/c and |y| < 0.5 and in p-Pb collisions at √ s NN = 5.02 TeV in 2 < p T < 12 GeV/c and −0.96 < y < 0.04.
Due to the short lifetime of the Λ c baryons (cτ = 60 µm [12]) and the statistical limitation of the data sample considered, the reconstruction of Λ c decays was particularly challenging. Three decay chan-Λ + c production in ALICE ALICE Collaboration nels of the Λ c were therefore studied, two hadronic channels (Λ + c → pK − π + and Λ + c → pK 0 S ), and a semileptonic one (Λ + c → e + ν e Λ). Furthermore, several different independent analysis strategies were developed, including the use of a Bayesian approach for particle identification [58] and a Multivariate Analysis (MVA) [59]. These developments build on top of the tools and strategies used in previous ALICE analyses of D-meson hadronic decays [7,27,28,[60][61][62] and of the Ξ c -baryon semileptonic decay [63]. After a description of the detector and the data samples in Sec. 2, we detail the different analyses and methods used for the various decay channels and collision systems in Sec. 3. The efficiency corrections applied and the treatment of the feed-down correction are described in Sec. 4. The evaluation of the systematic uncertainties is presented in Sec. 5. Finally, the results are presented and discussed in Sec. 6. Here, the cross section measured in pp collisions and the Λ + c /D 0 production ratio are compared with pQCD calculations and predictions from event generators as well as with existing measurements in different collision systems and rapidity intervals. The cross section obtained in p-Pb collisions is compared with the pp results, and the first measurement of the Λ + c nuclear modification factor in p-Pb collisions, R pPb , is presented.

Experimental setup and data samples
A comprehensive description of the ALICE apparatus and its performance can be found in [57,64]. In this section, the detectors used for the analyses discussed in this paper are described. Λ c baryons were measured by reconstructing their decay products in the pseudorapidity interval |η| < 0.8 relying on the tracking and particle identification (PID) capabilities of the central-barrel detectors, which are located in a solenoid magnet providing a B = 0.5 T field, parallel to the beam direction (z-axis in the ALICE reference frame). In particular, the Inner Tracking System (ITS) and the Time Projection Chamber (TPC) were utilised for track reconstruction, while PID was performed based on the information from the TPC and the Time-Of-Flight detector (TOF).
From the innermost radius of 3.9 cm (distance from the centre of the beam vacuum tube) to the outermost radius of 43.0 cm, the ITS cylinder includes two layers of Silicon Pixel Detector (SPD), two Silicon Drift Detector layers, and two Silicon Strip Detector layers. The different ITS detectors have full azimuth but different pseudorapidity coverage, with a common |η| < 0.9 acceptance. The spatial precision of the ITS detector, its vicinity to the beam pipe, and its very low material budget [65] allow for a precise determination of the track impact parameter (i.e. the distance of closest approach of the track to the primary vertex) in the transverse plane, for which a resolution better than 75 µm is achieved for tracks with p T > 1 GeV/c [65].
The TPC is the main tracking detector of the experiment and surrounds the ITS with an active radial range from 85 cm to 250 cm and with full azimuthal coverage in the pseudorapidity interval |η| < 0.9. It provides up to 159 space points to reconstruct the particle trajectory and determine its momentum. Additionally, it provides particle identification via the measurement of the specific energy loss, dE/dx. The TOF (an array of 1593 Multi-gap Resistive Plate Chambers) completes the set of detectors used for PID in the analyses presented in this paper. It is located at a radial distance of about 3.8 m, covering full azimuth in the pseudorapidity interval |η| < 0.9. The particle arrival time at the detector is determined with a resolution of about 80 ps. The T0 consists of two arrays of Cherenkov counters, located on both sides of the interaction point at +350 cm and −70 cm from the nominal vertex position along the beam line. The time resolution of the T0 in pp and p-Pb collisions is about 50 ps for the events in which the measurement is made on both sides [66]. The event time of the collision is obtained on an event-by-event basis either using the TOF detector, or the T0 detector, or a combination of the two [66].
The results presented in this paper were obtained from the analysis of the RUN 1 data collected by ALICE in pp collisions at √ s = 7 TeV in 2010 and in p-Pb collisions at √ s NN = 5.02 TeV from the 2013 data taking campaign. During the p-Pb run, the beam energies were 4 TeV for protons and 1.59 TeV per Λ + c production in ALICE ALICE Collaboration nucleon for lead nuclei. With this beam configuration, the proton-nucleon centre-of-mass system moves in rapidity by ∆y = 0.465 in the direction of the proton beam.
The V0 detector, used for trigger and event selection, consists of two scintillator arrays, called V0A and V0C, covering the full azimuth in the pseudorapidity intervals 2.8 < η < 5.1 and −3.7 < η < −1.7, respectively. The analyses used events recorded with a minimum bias (MB) trigger, which was based on the signals from the V0 and SPD detectors. At least one hit in either of the two scintillator arrays of the V0, or at least one hit in the SPD (pseudorapidity coverage of |η| < 2 and |η| < 1.4 for the inner and the outer layers, respectively) was required by the MB-trigger condition during the pp data taking, while in p-Pb the requirement was based on coincident hits in both V0A and V0C. The events were further selected offline using the SPD, V0 and Zero Degree Calorimeter (ZDC) information in order to remove background from beam-gas collisions, and from the machine as described in [67,68]. In the analysed sample, events with more than one interaction (pile-up) were removed according to the vertex information reconstructed from the hits in the SPD detector. To maximise the ITS acceptance, only events with a z-coordinate of the reconstructed vertex position within 10 cm from the nominal interaction point were used. With these requirements, approximately 300 and 370 million MB triggered events were analysed for the pp hadronic and semileptonic channels, respectively, corresponding to an integrated luminosity of L int = 4.8 and 5.9 nb −1 with an uncertainty of ± 3.5% [69], while approximately 100 million MB triggered events were selected for the p-Pb analyses, corresponding to L int = 47.8 µb −1 (± 3.7% [70]).

Λ c analysis overview and methods
The measurement of Λ c production was performed by reconstructing three decay modes: Λ + c → pK − π + with branching ratio (BR) equal to (6.35 ± 0.33)%, Λ + c → pK 0 S with BR = (1.58 ± 0.08)% and K 0 S → π + π − with BR = (69.20 ± 0.05)%, and Λ + c → e + ν e Λ with BR = (3.6 ± 0.4)% and Λ → pπ − with BR = (63.9 ± 0.5)% [12]. The hadronic decays were fully reconstructed while the semileptonic decay was partially reconstructed because the neutrino is not detectable with the ALICE setup. The analysis strategy for the extraction of the Λ c signals from the large combinatorial background was based on the reconstruction of charged tracks with the central-barrel detectors, on the V-shaped neutral decay topology reconstruction (V 0 ) of K 0 S and Λ, on kinematical and geometrical selections, and on the use of PID on the decay tracks.
These analyses cannot fully benefit of the reconstruction and selection of secondary vertex topologies due to the comparable resolution of the ITS on the track impact parameter and the mean decay length of the Λ c . The use of PID techniques is therefore fundamental to reduce the large combinatorial background. The identification of pions, kaons, protons, and electrons used for the Λ c analyses in all the considered decay channels and for both colliding systems was based on the information from the specific energy loss dE/dx in the TPC detector and on the time of flight measured with the TOF detector. For some of the results presented here, MVA techniques were applied additionally to the selection procedure based on classical cuts and called "standard" (STD) in the following. Finally, the Λ c raw yield was extracted with an invariant mass analysis for the hadronic decay modes or, in the semileptonic analysis, by counting the candidates with the correct combination of particle species and charge sign (i.e. e + Λ and e − Λ), indicated as "right sign" in the following, after subtracting the background estimated from "wrong sign" pairs (i.e. e − Λ and e + Λ). Table 1 summarises the various analysis methods.
Simulations were used in the analyses to determine the geometrical acceptance, the efficiencies of track reconstruction and Λ c selection, and the line shape of the Λ c invariant-mass peak. The event generator used to simulate pp collisions was PYTHIA 6.4.21 [71] with the Perugia-0 tune [72]. For p-Pb collisions, PYTHIA events containing a cc or bb pair were merged with events simulated with the HIJING 1.36 event generator [73] to obtain a better description of the multiplicity distribution observed in data. The n σ ,Bayes Table 1: Λ c decay channels studied and analysis methods presented in this paper. generated particles were transported through the ALICE detector using the GEANT3 package [74].
For all the analyses, the lower limit of the Λ c p T interval in which the signal could be extracted was imposed by the large combinatorial background, which could not be reduced enough with the applied selections. The upper limit was imposed by the limited size of the analysed data sample. This section gives an overview of the analysis methods, with Sec. 3.1 dedicated to the Λ c hadronic decay modes and Sec. 3.2 to the semileptonic channel.

Hadronic decay modes
The Λ + c → pK − π + candidates were built from triplets of reconstructed tracks with proper charge-sign combination. The Λ + c → pK 0 S candidates were constructed by combining a reconstructed track (the bachelor) with a K 0 S candidate. The Λ c and K 0 S candidates were formed by combining reconstructed tracks having |η|<0.8 and at least 70 associated space points in the TPC. Additionally, for the bachelor and the tracks used to form Λ + c → pK − π + candidates, at least one hit in either of the two SPD layers was required. The K 0 S candidates were identified by applying selections on characteristics of their decay tracks (p T > 0.1 GeV/c, a minimum transverse impact parameter to the primary vertex, d 0 , of 0.05 cm and a maximum distance of closest approach between the daughters tracks of 1.5 cm) and of their weak decay topology (a minimum transverse decay radius of 0.2 cm and a minimum cosine of the V 0 pointing angle to the primary vertex of 0.99). The invariant mass of the π + π − pair was required to be compatible with the PDG mass of the K 0 S within 1 or 2 σ depending on the p T interval and the collision system. To further improve the K 0 S signal purity, especially at lower p T , veto selections on Λ, Λ and γ PDG masses were applied to the invariant masses calculated with the pπ − , pπ + and e + e − hypotheses for the daughter tracks, respectively.
For both decay channels, cuts on kinematical and geometrical variables were also applied after a tuning procedure in each p T interval. The kinematical variables include the p T of the daughter tracks and the p T of the K 0 S in the Λ + c → pK 0 S analysis. In the Λ + c → pK − π + analysis, the geometrical variables include the separation between the interaction point and the points of closest approach of the opposite-sign track pairs, the separation between the reconstructed Λ c -decay vertex and the interaction point (decay length), the distance of closest approach of the three pairs of tracks, the quadratic sum of the minimum distances of the tracks from the reconstructed Λ c -decay vertex, and the Λ c pointing angle to the primary vertex. In the Λ + c → pK 0 S analysis, the geometrical variables include the upper cuts on the d 0 of the bachelor and K 0 S (applied to remove secondary tracks originating very far from the interaction point). For both decay channels the cuts were tuned on Monte Carlo samples for each analysis to achieve a high statistical significance in each p T interval. After the selection, the acceptance in rapidity for Λ c baryons drops steeply to zero for |y lab | > 0.5 at low p T and for |y lab |> 0.8 at p T > 5 GeV/c, where y lab is the rapidity in the laboratory frame. A p T -dependent fiducial acceptance cut was therefore applied on the Λ c rapidity, |y lab | < y fid (p T ) with y fid (p T ) increasing from 0.5 to 0.8 in the interval 0 < p T < 5 GeV/c, and y fid = 0.8 for p T > 5 GeV/c, as described in [7].
The identification of the proton in the Λ + c → pK 0 S analysis was based on the dE/dx and time-of-flight information, using as a PID-discriminating variable the difference between the measured signal and that expected under the proton mass hypothesis divided by the detector resolution (n σ ), as detailed in [58]. Figure 1 shows an example of the n σ distributions relative to the proton hypothesis as a function of momentum for TOF and TPC signals in p-Pb collisions. To reduce the pion and kaon contamination, for tracks with momentum p < 1 GeV/c, a |n σ | < 2 selection with respect to the proton hypothesis was applied on the TPC dE/dx. For p > 1 GeV/c, in order to improve the signal over background ratio, the presence of the TOF signal was requested and the track was required to be within |n σ | < 3 of the expected proton TOF signal, without any further selections based on TPC information. In this momentum region, tracks missing the TOF information were discarded. In the p-Pb analysis it was further required that the track should be within |n σ | < 3 of the expected TPC signal.
In the Λ + c → pK − π + analysis, where a larger combinatorial background is present, the Bayesian PID method [58] was adopted to increase the purity of the signal. In this method, the signals from the TOF and TPC are combined constructing a conditional probability that a given track corresponds to a given hadron species (p, K or π) based on a set of measurements in the two detectors. The computation of the Bayesian probability entails the use of priors, that are evaluated with data-driven techniques. This approach provides a smoother increase of the PID efficiency with p T than the one observed with the n σ -cut approach and it makes the best possible use of the combined information coming from the two detectors. To each of the three Λ c decay tracks, a single mass hypothesis was assigned, corresponding to the hadron species (p, K, π) for which the Bayesian probability was found to be the maximum. Candidates were rejected if the daughter-track species and charge sign did not match with a pK − π + (or charge conjugate) final state. This corresponds to the "maximum probability" strategy discussed in [58] that was, for example, successfully validated in reproducing the published results [60] for the D 0 → K − π + production cross section, which were obtained with a PID strategy based on a |n σ | < 3 selection.
Λ + c production in ALICE ALICE Collaboration In addition to the STD analyses for the study of the hadronic decay modes in p-Pb collisions, a further analysis was carried out that relies on a multivariate selection to separate the background from the signal, based on Boosted Decision Trees (BDT) [59]. This approach will be indicated as MVA in the following.
To train the algorithm, the signal sample was built using Λ c particles from Monte Carlo simulations.
For the background sample, as detailed later, both real and simulated events were used. This training sample was used to determine a mapping function, which describes a decision boundary, optimised in order to maximise signal/background separation. The learned mapping function was then applied to a real data sample, in which the type of candidate is unknown. A cut on this decision boundary aims to reject background candidates while keeping signal candidates.
Prior to the BDT decision, for both decay channels, PID selections were applied. For the Λ + c → pK − π + analysis a |n σ | < 3 cut was applied on the compatibility with the expected dE/dx and time-of-flight values. For proton and kaon identification, tracks without a TOF signal were identified using only the TPC, and tracks with incompatible TPC and TOF identifications were assigned the identity given by the TOF. For pion identification only the TPC was used. In the case of the Λ + c → pK 0 S analysis, a |n σ | < 3 compatibility cut was applied on the TPC and TOF, when available, for the bachelor track. For this analysis, an additional cut in the Armenteros-Podolanski space [75] was also applied in order to reject Λ decays.
Independent BDTs were trained per p T interval and applied on the p-Pb data sample. The BDTs were trained using signal samples consisting of Λ c decays from simulated events, required to have at least one Λ c per event decaying to either a pKπ or pK 0 S final state, and including a detailed description of the detector response, the geometry of the apparatus and the conditions of the luminous region. The background sample was taken from the sidebands of the candidate invariant-mass distribution in the data (pKπ analysis), or from the simulated events (pK 0 S analysis), and it was verified that swapping the simulated/real background sample does not change the result of the trained BDT.
For the Λ + c → pK 0 S analysis the variables related to the decay topology that were used in the multivariate analysis include the p T of the bachelor track, the d 0 of the bachelor track, the V 0 invariant mass under the hypothesis that the daughters are a π + π − pair, the d 0 and the lifetime of the V 0 . For the Λ + c → pK − π + analysis the variables related to the decay topology that were used in the multivariate analysis include all variables used in the Λ + c → pK − π + STD analysis, as well as the projection of the decay length in the transverse plane normalised by its error. PID variables were also used in both analyses, namely the Bayesian probabilities that each track is correctly identified as either a proton, a kaon, or a pion for the pKπ analysis, and the Bayesian probability that the bachelor track is a proton for the pK 0 S analysis. Figure 2 shows examples of the BDT response in the two lowest p T intervals for the analysis of the Λ + c → pK − π + decay channel. The raw signal yields were extracted by fitting the invariant mass distributions of the Λ c candidates passing the selections outlined above, for every p T interval under study. The fitting function consists of a Gaussian describing the signal, whose width was fixed to the value obtained in the simulation, and a polynomial of second order or a linear function (with the choice depending on the p T interval) to describe the background.

Semileptonic decay mode
The Λ c production cross section in pp collisions at √ s = 7 TeV was also measured from its semileptonic decay Λ + c → e + ν e Λ, based on reconstructed e + Λ pairs. This analysis follows a procedure similar to the one presented in our recent work on the measurement of Ξ 0 c via its semileptonic decay, Ξ 0 c → e + Ξ − ν e [63]. Here, we briefly describe the analysis approach for the Λ c with an emphasis on the differ- Figure 2: Normalised distribution of the BDT responses of the Λ c candidates for Monte Carlo signal (blue area) and background (red shaded area) in two p T intervals for the Λ + c → pK − π + decay channel in p-Pb collisions where the MVA method was used. The arrows correspond to the applied cuts. ences from that analysis. Λ + c candidates are defined from e + Λ pairs by combining a track originating from the primary vertex, denoted electron track in the following, and a Λ baryon reconstructed through the decay Λ → pπ − , by exploiting the fact that its V 0 -shaped decay topology is significantly displaced from the interaction point, given the additional lifetime of Λ hyperons, cτ = 7.89 cm [12]. The V 0 candidates are built from pairs of tracks with |η| < 0.8 reconstructed in the TPC and the ITS provided that they pass reconstruction quality criteria in a similar way as done for the hadronic decay channels. Additional cuts were applied to select the V 0 -shaped decays: distance of closest approach between the daughter tracks smaller than 1 cm, |d 0 | of the daughter tracks larger than 0.06 cm, and cosine of the V 0 pointing angle to the primary vertex larger than 0.99. The compatibility of the pπ − invariant mass with the Λ-baryon mass within 8 MeV/c 2 was required in the analysis. The Λ sample obtained with these selections is characterised by a signalto-background ratio of about 20 for p T > 0. Electron tracks were required to satisfy the reconstruction quality criteria described in [63]. The PID selection was based, with respect to the electron hypothesis, on a |n σ | < 3 cut on the TOF signal and a p T dependent n σ cut on the TPC signal: (−3.9 + 1.2p T − 0.094p 2 T ) < n σ < 3, with p T expressed in GeV/c. The p T -dependent lower limit for the TPC n σ is defined to have a constant purity over the measured p T interval. Reconstructed e + Λ pairs were further required to have an opening angle smaller than 90 degrees and an invariant mass smaller than the Λ c mass.
Due to the missing neutrino, the invariant-mass distribution of eΛ pairs does not show a peak at the Λ c mass and the raw yield cannot be extracted via a fit to the invariant-mass distribution with signal and background components as done for the hadronic decay channels. Here, similarly to [63], the background contributions were estimated using the fact that Λ + c baryons decay only into e + Λ pairs, denoted as rightsign (RS), and not into e − Λ pairs, denoted as wrong-sign (WS), while background candidates contribute to both RS and WS pairs. The Λ c raw yield distribution was obtained by subtracting the WS contribution from the RS yields. Other contributions to eΛ pairs, such as the contributions of Λ 0 b semileptonic decays to WS pairs and of Ξ 0,+ c decays to RS pairs, are corrected after the subtraction. The obtained Λ c raw yield in the intervals of eΛ-pair momentum are further corrected for the missing momentum of the neutrino, as discussed below. Figure 5 shows the uncorrected eΛ invariant-mass distributions for WS and RS pairs for the interval 3 < p eΛ T < 4 GeV/c.
The Ξ 0,+ c baryons contribute to RS pairs through the decay chain Ξ 0,+ c → e + Ξ −,0 ν e → e + Λπ −,0 ν e . This contribution was estimated and subtracted from the RS yield to extract the yield of eΛ pairs originating from Λ + c decays. First, the ratio of eΛ pairs from Ξ 0 c and Ξ + c was determined. Assuming that the production of Ξ 0 c and Ξ + c is the same, the difference in the eΛ pair yields arises from their different branching ratios into the relevant decay modes. The ratio BR(Ξ + c → e + Ξ 0 ν e )/BR(Ξ 0 c → e + Ξ − ν e ) was measured by CLEO in e + e − collisions below ϒ(4S) energies and found to be 2.46 ± 0.7 +0.33 −0.23 [76]. Then, the relative contribution of Ξ 0,+ c decays to the total yield of RS pairs was calculated. This was done using two different methods. In the first method, the Ξ 0,+ c contribution in the p eΛ T distribution was calculated as where N i is the yield of eΛ pairs in i-th p eΛ T bin, M j is the number of Ξ 0 c in j-th p Ξ 0 c T bin, which is computed from the measured Ξ 0 c cross section [63] as detailed below, and F Ξ 0,+ c i j are the matrices taking into account the reconstruction and selection efficiencies and the decay kinematics to convert p Ξ 0,+ c T into p eΛ T . The Ξ 0 c cross section in the p T range 1 < p T < 8 GeV/c was taken from the measurement reported in [63] and the cross section outside the measured p T range was estimated using the Tsallis function, Λ + c production in ALICE ALICE Collaboration where C is a normalisation constant, m is the Ξ 0 c baryon mass, and the parameters n and T were extracted from a fit to the data in the measured p T range. The ratio between the yield of eΛ pairs from Ξ c decays and that of inclusive eΛ pairs was found to be independent of p eΛ T in the measured interval, with an average value of 0.38 ± 0.10, where the uncertainty also includes the contribution from the branching ratios measured by CLEO.
The second approach exploits the fact that the distance between the interaction point and the decay vertex of Λ baryons originating from Λ c decays is on average smaller than that of Λ baryons from Ξ c decays, mediated by Ξ hyperons (cτ ∼ 4.91 cm [12]). In detail, for each p eΛ T interval, the Ξ 0,+ c fraction was determined by fitting the measured distribution of the distance of the baryon decay point from the interaction point with the two contributions of Λ baryons originating from Λ + c and Ξ 0,+ c decays generated with PYTHIA6.4.21 (Perugia-0 tune) [72]. Also in this case, no p T dependence of the Ξ 0,+ c relative contribution in the yield of eΛ pairs was observed, and the average was found to be 0.52 ± 0.09, consistent with the result from the first approach. By taking the weighted average of the values obtained with the two methods, we obtained 0.46 ± 0.06 as the relative contribution of Ξ 0,+ c decays.
2)%, followed by the subsequent decay Λ + c → Λ + X, with BR (35 ± 11)% [12]. This contribution was estimated using the Λ 0 b measurement at central rapidity by CMS [77], which covers the transverse momentum interval p T > 10 GeV/c. The cross section for p T < 10 GeV/c was estimated using the Tsallis parameterisation reported in [77]. The Λ 0 b distribution was further converted into an eΛ distribution via simulations, taking into account the detector acceptance, the reconstruction and selection efficiency, and the decay kinematics to determine the fraction of Λ 0 b momentum carried by eΛ pairs. The obtained yield of eΛ pairs originating from Λ 0 b decays was added to the measured eΛ yield after the WS pairs were subtracted. The correction is found to increase with p eΛ T reaching about 10% in the highest p eΛ T interval. The correction for the missing momentum of the neutrino was performed by using the response matrix determined with the full detector simulation of PYTHIA events containing Λ c baryons and using the Bayesian unfolding technique [78] implemented in the ROOUNFOLD package [79]. The number of iterations, which is a regularisation parameter of the Bayes unfolding, was chosen to be 3 in this analysis. It was verified that the final result is not sensitive to this choice.

Corrections
The p T -differential cross section of prompt Λ + c baryon production was obtained for each decay channel as: where N Λ c is the raw yield (sum of particles and antiparticles) in a given p T interval with width ∆p T , f prompt is the fraction of the raw yield from prompt Λ c , (A × ε) is the product of acceptance and efficiency for prompt Λ c baryons, BR is the branching ratio for the considered decay mode and L int is the integrated luminosity. The correction factor for the rapidity coverage c ∆y was computed, for the hadronic decay modes, as the ratio between the generated Λ c -baryon yield in |y lab | < y fid (p T ) and that in |y lab | < 0.5. For the semileptonic decay analysis, the rapidity of the Λ c candidate cannot be calculated due to the missing neutrino momentum, and the y fid cut cannot be applied. A factor c ∆y = 1.6 was used in this case assuming a flat distribution of the Λ c candidates in |y lab | < 0.8, which was verified with an accuracy of 1% using pure Monte Carlo information from PYTHIA. The factor 2 in the denominator of Eq. 3 takes into account that the raw yield is the sum of particles and antiparticles, while the cross section is given for particles only and is computed as the average of Λ + c and Λ − c .
For hadronic decays the solid lines correspond to the prompt Λ c , while the dotted lines represent (A × ε) for Λ c baryons originating from beauty-hadron decays. The efficiency for semi-leptonic decays (same for both prompt and non-prompt Λ c ) is represented with one solid line. The statistical uncertainties are smaller than the marker size.
The correction for the detector acceptance and reconstruction efficiency (A × ε) was obtained following the same approach as discussed in [60]. The correction factors were obtained from Monte Carlo simulations where the detector and data taking conditions of the corresponding data samples were reproduced.
Contrary to the case of pp collisions, for which the simulation describes in a satisfactory way the chargedparticle multiplicity in data, in p-Pb collisions a weighting procedure based on the event multiplicity was applied in the calculation of the efficiency from the simulated events. This approach accounts for the dependence of the efficiency on the event multiplicity, which is due to the fact that the resolutions of the primary vertex position and of the variables used in the geometrical selections of displaced decay vertices improve with increasing multiplicity.
The efficiency was computed separately for prompt and non-prompt Λ c (originating from Λ b -baryon decays). The Λ + c → pK − π + decay channel includes not only the direct (non-resonant) decay mode, but also three resonant channels, namely pK * (892) 0 , ∆(1232) ++ K − and Λ(1520)π + . The kinematical properties of these decays are different, resulting in different acceptances and efficiencies for each case. The final correction was determined as a weighted average of the (A × ε) values of the four decay channels, using the relative branching ratios as weights. Figure 6 shows the product of acceptance times efficiency (A × ε) for Λ c baryons with |y| < y fid (p T ) in pp collisions at √ s = 7 TeV, as a function of transverse momentum, for Λ + c → pK − π + (left panel), Λ + c → pK 0 S (middle panel), and Λ + c → e + ν e Λ (right panel). The higher efficiency for Λ c from beautyhadron decays in the Λ + c → pK − π + decay channel is due to the geometrical selections on the displaced decay-vertex topology, which enhance the non-prompt component because of the additional lifetime of the beauty hadrons. In the case of the Λ + c → pK 0 S decay, for p T < 4 GeV/c the efficiency for prompt Λ c is slightly higher because the upper cut applied on the bachelor d 0 to remove secondary tracks rejects preferentially Λ c from beauty-hadron decays. In the semileptonic analysis no selection is made on variables related to the displacement of the Λ c decay vertex from the primary vertex, and therefore the efficiency is the same for both prompt and non-prompt Λ c .
When using the Multivariate Analysis approach, a further correction factor (ε BDT ) was required. This additional ingredient corresponds to the BDT cut efficiency, quantifying the fraction of true Λ c candidates accepted by the selection on the classifier output. Since the BDT analysis employed a different set of preselections, a specific correction factor ε presel for those was also taken into account. The final efficiency Λ + c production in ALICE ALICE Collaboration  correction is ε = ε BDT × ε presel .
The BDT cut efficiency was determined from the simulations with the PYTHIA and HIJING event generators described above by applying the classification algorithm resulting from the training of the BDT on the simulated sample enriched with Λ c described in Sec. 3.1.
The efficiency and acceptance corrections for prompt and non-prompt Λ c in p-Pb collisions are reported in Fig. 7 as a function of p T in the rapidity range |y lab | < y fid (p T ) for the decay channels Λ + c → pK − π + (left panel) and Λ + c → pK 0 S (right panel) for the MVA technique. The non-monotonic trend seen in the efficiencies for both channels is a result of the non-monotonic tightness of the BDT cut chosen as a function of p T , and it was verified that these choices do not have a significant systematic effect on the result.
To obtain the factor f prompt , i.e. the fraction of prompt Λ c in the raw yield, the production cross section of Λ c from Λ b decays was estimated using the beauty hadron p T shape from FONLL [3,4] as described in detail in [60] (the contribution from B-meson decays to Λ c was checked and found negligible [12]). The fraction of beauty quarks that fragment to beauty hadrons and subsequently decay into Λ c baryons f (b → Λ c ) = 0.073 was taken from [80] and the Λ b → Λ c + X decay kinematics were modelled using the EVTGEN [81] package. The production cross section of Λ c from Λ b was then multiplied for each decay channel in each p T interval by (A × ε) feed-down , the factor c ∆y , the branching ratio BR and the integrated luminosity L int . The correction factor f prompt was calculated in pp collisions as: where N Λ c /2 is the raw yield, which was divided by a factor of two to account for particles and antiparticles.
For p-Pb collisions, a hypothesis on the nuclear modification factor R feed-down pPb of Λ c from beauty-hadron decays was added as an additional factor in the last term of Eq. 4. As in the D-meson analyses [27], it was assumed that the R pPb of prompt and feed-down Λ c were equal and their ratio was varied in the range 0.9 < R feed-down pPb /R prompt pPb < 1.3 to evaluate the systematic uncertainties. The values of f prompt range Λ + c production in ALICE ALICE Collaboration between 95% and 99% depending on the decay channel and p T .

Evaluation of systematic uncertainties
This section is dedicated to the description of the various sources of systematic uncertainties for each analysis presented here. First, the systematic uncertainties for the Λ c hadronic decay modes in both pp and p-Pb collisions will be discussed. Then, the systematic uncertainties studied for the Λ c semileptonic decay mode will be presented. For each analysis, the different sources of systematic uncertainties were assumed to be uncorrelated among each other and the total systematic uncertainty was determined in each p T interval as the quadratic sum of the different contributions.
A summary of the systematic uncertainties is shown in Tabs. 2-4, for the hadronic analyses in pp collisions, the hadronic analyses in p-Pb collisions, and the semileptonic analysis in pp collisions, respectively. These include the uncertainties specific to each analysis as well as the uncertainties associated to the branching ratios of the Λ c decay modes [12]. The measured cross sections are also affected by a global normalisation uncertainty related to the determination of the integrated luminosity of 3.5% [69] and 3.7% [70] in pp and p-Pb collisions, respectively.

Systematic uncertainties for the hadronic channels
The systematic uncertainty on the raw-yield extraction was estimated for each decay mode and in each p T interval by repeating the fit to the invariant-mass distributions under different approaches. The following variations to the fit procedure were considered: (i) the background function, for which three different functions were tested (parabolic, linear and exponential), and (ii) the lower and upper limit of the fit range of the invariant-mass distributions. For each combination of the aforementioned variations, the fit was performed under different assumptions on the width and position of the Gaussian function modelling the Λ c peak in the invariant-mass distributions, namely: (a) fixing the Gaussian width to the value obtained from simulation (used as default); (b) fixing the peak position to the value obtained from simulations; (c) leaving the peak width and position as free parameters of the fit; (d) fixing both the peak width and position. Only those cases satisfying quality criteria on the resulting fits were considered to assess the final systematic uncertainty, which was defined as the RMS of the distribution of the signal yields obtained from the different trials.
The contribution to the systematic uncertainty due to the tracking efficiency was evaluated as discussed in [7] for the D-meson analysis, i.e. by comparing the probability of matching TPC tracks to ITS points in data and simulation and by varying the quality cuts to select the tracks used in the analysis. The uncertainty on the matching efficiency was defined as the relative difference of the matching efficiencies in data and simulations. The matching efficiency for primary tracks is higher than that for secondary tracks produced far from the interaction point in strange particle decays (such as those coming from the K 0 S decay in the Λ + c → pK 0 S channel) or in interactions with the detector material. Different fractions of primary and secondary tracks, in data and simulations, could lead to a wrong estimation of the systematic uncertainty for the matching. For this reason, the comparison of the matching efficiency in data and simulations was done after weighting the relative abundances of particles in Monte Carlo to match those observed in data. The uncertainty resulting from these studies was added in quadrature with the uncertainty on the track selection for the final uncertainty on the tracking efficiency.
Systematic uncertainties on the efficiency can also arise from possible differences in the distributions and resolution of selection variables between data and the simulation. The systematic effect induced by these imperfections was estimated by repeating the analysis with several sets of selection criteria for the Λ c candidates. Each selection was varied with respect to the central value, obtaining a relative variation of the efficiency between 5% and 40%. The uncertainty due to these selections was then estimated from the RMS of the cross sections resulting from all the variations and it ranges from 4% to 10% depending Λ + c production in ALICE ALICE Collaboration on the analysis and the decay channel.
The results presented in this paper rely on an extensive use of the PID capabilities of the TPC and TOF detectors. The uncertainties arising from discrepancies in the PID efficiency in data and simulation were estimated by varying the PID strategy (with tighter or looser n σ cuts, or with different configurations for the Bayesian PID, for the Λ + c → pK 0 S and Λ + c → pK − π + analyses, respectively), and estimating the uncertainty from the RMS of the resulting corrected yields obtained from the tests.
The efficiencies determined from the simulations depend on the generated p T distribution of Λ c baryons. The central values of the correction factors were obtained by re-weighting the Λ c distribution generated by PYTHIA according to the ratio of the p T distribution of D 0 mesons from FONLL calculations and from PYTHIA simulations. A systematic uncertainty was defined by considering the RMS variation of the efficiencies determined with different generated p T shapes, namely: (i) c-quark p T distributions from FONLL, (ii) Λ c p T shapes from PYTHIA. It was found to be 3% at most, depending on the analysis.
As discussed in Sec. 4, the efficiency for Λ c reconstruction and selection depends on the multiplicity of particles produced in the collision, since the resolution on the primary vertex improves with increasing multiplicity. For p-Pb collisions, a systematic uncertainty was assigned to account for the accuracy of the multiplicity weighting procedure applied in the efficiency calculation. It amounts to 1% for the analysis using MVA, while it is negligible for the STD analysis, for which the efficiency shows a less pronounced dependence on multiplicity.
The contribution to the uncertainties coming from the subtraction of Λ c baryons from Λ b decays was calculated as the envelope of the uncertainty bands obtained (i) by varying the p T -differential cross section of beauty hadrons within the theoretical uncertainties of the FONLL calculation, and (ii) with the same method but after scaling by a factor of two the fraction f (b → Λ c ), which is used together with FONLL cross sections to determine the yield of Λ c from Λ b decays. The uncertainty in the FONLL calculation of (i) was determined by changing the b-quark mass and the perturbative scales, as explained in [4], also including the uncertainty on f (c → Λ b ) from [80], and finally adding in quadrature the uncertainty estimated for the used PDF set. The variation by a factor of two of the fraction f (b → Λ c ) in (ii) was motivated by the observation that FONLL calculations describe the available Λ b cross section measurements in pp collisions at √ s =7 TeV once the value of 0.197 measured at CDF [82] is taken. As noted in Sec. 1 the different values of this fragmentation fraction measured in hadron-hadron collisions with respect to e + e − interactions has been interpreted as a violation of its universality [12]. If the value f (b → Λ b ) = 0.088, derived from LEP measurements in electron-positron collisions [80], is used for the fragmentation fraction, the FONLL calculations underestimate by a factor of about two the Λ b measurement by LHCb at forward rapidity in the same p T region of this analysis [83] and by a factor of about 1.6 the CMS measurements at mid-rapidity in their lowest reported p T interval (10 < p T < 13 GeV/c) [77], corresponding to the high-p T region of this analysis.
Additional possible sources of systematic uncertainties were checked. The difference between the resolution on the K 0 S mass in data and simulation, the difference in the reconstruction efficiencies for Λ + c and Λ − c , and the possible contamination in the Λ c invariant-mass distribution coming from D + → π + K 0 S and D + s → K + K 0 S decays were all checked and proved to give a negligible contribution to the final uncertainties. These decays enter the candidate Λ c sample only if the kaon or the pion passes the proton PID selection.
For the analyses using MVA, specific sources of systematic uncertainty were additionally considered. The uncertainty associated to the selection on the MVA classifier output was estimated by repeating the analysis with different cutting points after verifying that these variations induce a significant modification of the efficiency, between 10 and 40%. The RMS of the distribution of the corrected yields was then used to assign the systematic uncertainty (reported under cut efficiency in Tab. 3). Branching ratio (%) 5.1 5.0 Luminosity (%) 3.7 Table 2: Summary of relative systematic uncertainties for the lowest and highest p T intervals considered in the analysis, for the two Λ c hadronic decay modes in pp collisions. When the uncertainty was found to be < 1%, it was considered negligible ("neg." in the table).
Branching ratio (%) 5.1 5.0 Luminosity (%) 3.5 Table 3: Summary of relative systematic uncertainties for the lowest and highest p T intervals considered in the analysis for the two Λ c hadronic decay modes and the two analysis techniques in p-Pb collisions. When the uncertainty was found to be < 1%, it was considered negligible ("neg." in the table).
A possible systematic effect of the specific multivariate algorithm chosen (BDT) [59] was checked by changing the configuration of the MVA method. These changes included the number of trees used to construct the forest, the maximum depth of the trees constructed, the boosting algorithm, the application of data preprocessing such as the transformation of input variables to reduce correlation or the transformation of the variable shapes into more appropriate forms, the metric defining the separation criterion in the node and the number of input variables. The effects of such modifications in the corrected yields were found to be negligible.
The PID-related variables play an important role in the multivariate selection, since they offer the largest discrimination power. As a further cross-check, the systematic uncertainty associated with the inclusion Λ + c production in ALICE ALICE Collaboration of these variables in the multivariate selection was estimated. For the Λ + c → pK − π + analysis, the kaon and pion priors used in the calculation of the Bayesian probability were modified conservatively based on the maximum mismatch between the default priors determined through an iterative procedure and the measured particle abundances [58], and the BDT efficiency was determined for each modification. For the Λ + c → pK 0 S analysis, the Bayesian probability for protons in simulation and data was compared using the daughter particles of V0 decays in order to select a pure proton sample. The resulting variations were found to be consistent within 2-4%; this effect was not included as an additional uncertainty source, since it should be accounted for in the BDT cut variation and its magnitude is smaller than the assigned systematic uncertainty. Moreover, to assess whether the Bayesian approach used in the MVA might lead to biased results, the Λ + c → pK 0 S analysis was repeated using an n σ approach for the bachelor PID and not considering any PID in the BDT. The results for the three cases were found to be compatible, and therefore no systematic uncertainty was assigned. As reported in Sec. 3.1 and Sec. 4, a loose particle identification, based on rectangular n σ -compatibility cuts on the TPC and TOF PID response for pion, kaon and proton tracks is applied prior to the BDT. The systematic uncertainty associated with this cut was studied by comparing the corrected Λ c yield obtained with and without this cut (Λ + c → pK 0 S ) and without the TOF selection (Λ + c → pK − π + ) and was found to be negligible in the p T range considered here.
The contribution of the uncertainty related to the imperfect description of the impact parameter resolution in the simulation, which could affect the input variables related to vertex reconstruction, was checked in the Λ + c → pK − π + analysis. For this check, the distribution of the input variables was altered by smearing the reconstructed track parameters to match the impact parameter resolution observed in data, and the BDT cut efficiency was recalculated. The change in efficiency was 2% at low p T , and less than 1% at high p T , consistent with being a contribution to the systematic uncertainty estimated with the cut-variation procedure.

Systematic uncertainties for the semileptonic channel
The following contributions to the systematic uncertainty on the Λ c cross section measurements through the Λ + c → eΛν e decay channel were considered: raw-yield extraction, (A × ε) correction factor, correction for the missing neutrino momentum and for feed-down from beauty-hadron decays. These contributions were added in quadrature to obtain the total systematic uncertainty in each p T interval and they are summarised in Tab. 4.
The systematic uncertainty due to the raw-yield extraction includes the uncertainties in the WS subtraction procedure, the estimation of the Ξ 0,+ c contribution to RS pairs and the Λ 0 b contribution in WS pairs. The WS pair subtraction described in Sec. 3.2 was based on the assumption that there were no charge asymmetric background sources and that the acceptance of RS and WS pairs were the same. The influence of the charge asymmetric background sources was evaluated using PYTHIA events with full detector simulation, as done in the Ξ 0 c analysis [63], and found to be about 2%. The difference in the acceptance of RS and WS pairs was estimated using a mixed-event technique and found to be negligible for this analysis. In addition, the impact on the background subtraction of the hadron contamination in the electron sample and the signal-to-background ratio was studied varying the electron identification criteria. The corrected spectra were all found to be consistent with the one obtained with the default selections and no systematic uncertainty was assigned.
The Ξ 0,+ c contribution to the RS pairs calculated as described in Sec. 3.2 also contributes to the systematic uncertainty on the raw-yield extraction. An additional uncertainty of 10% estimated from PYTHIA simulations was assigned to take into account the possible p T dependence of the fraction of Λ + c from Ξ 0 c decays and summed in quadrature to the value reported in Sec. 3.2. The systematic uncertainty on the Λ 0 b contribution in WS pairs was estimated by taking into account the uncertainty on the Λ 0 b cross section measured by CMS [77] and the uncertainty on the relevant branching ratios. The uncertainty was found Branching ratio (%) 11 Luminosity (%) 3.7 Table 4: Summary of relative systematic uncertainties for the Λ + c → e + ν e Λ analysis in pp collisions. The uncertainties smaller than 1% are considered negligible ("neg." in the table).
to increase with p T reaching about 5% in the highest p T interval.
The (A × ε) factor could be affected by imperfections in the description of the detector alignment and response in the simulation. The systematic uncertainties due to the reconstruction and selection efficiency were estimated by repeating the analysis with different selection criteria for electrons, Λ baryons, eΛ pairs and by comparing the corrected yields. The systematic uncertainty on the electron reconstruction and selection efficiency was estimated via variations of the track-quality criteria and the PID selections for electron identification. The RMS of the Λ c corrected yields, which amounted to 4% (track-quality) and 3% (PID), was assigned as a systematic uncertainty. Similarly, a systematic uncertainty of 1% on the Λ reconstruction and selection efficiency, was estimated from the RMS of the inclusive Λ corrected yield against variations of the criteria applied to select the Λ decay tracks and its V 0 decay topology. In addition, a systematic uncertainty of 4% on the Λ efficiency due to possible imperfections in the description of the detector material in the simulations was considered and summed in quadrature to the one estimated from the variation of the selection criteria. The uncertainties on the electron and Λ reconstruction efficiency were considered as correlated and combined linearly. The uncertainty on the eΛ pair selection efficiency was estimated by varying the selection criteria on the opening angle and the invariant mass of the pair and a systematic uncertainty ranging from 1 to 25% was assigned depending on p T . The systematic uncertainty due to a possible imperfect description of the acceptance of eΛ pairs in the simulation was estimated to be 11% by comparing the azimuthal distribution of inclusive electrons and Λ baryons in data and in the simulation. The uncertainty on the eΛ pair acceptance was summed in quadrature to the one on the electron, Λ and eΛ-pair selection efficiencies.
The dependence of the corrected results on the unfolding procedure was tested by (i) using as prior for the Bayesian unfolding the p T distribution from PYTHIA Monte Carlo simulations, and (ii) adopting different unfolding methods (χ 2 minimisation with regularisation [84,85] and Singular Value Decomposition [86]). The RMS of the corrected yields was used to estimate the resulting uncertainty, which increases from 3% to 11% towards higher p T .
The uncertainty arising from the subtraction of the feed-down from beauty-hadron decays was calculated in the same way as for the hadronic decays.

Results
In this section, results are first presented in Sec. 6.1 for the prompt Λ + c production cross sections in pp and p-Pb collisions obtained using the procedure discussed in Secs. 3-5. In the decay modes under study in pp collisions, it was possible to extract a stable signal in the lowest p T interval (1 < p T < 2 GeV/c) only via the semileptonic decay. In the highest p T interval (6 < p T < 8 GeV/c) it was not possible to extract a signal from the Λ + c → pK − π + invariant mass distribution. For p-Pb collisions in the two hadronic Λ + c production in ALICE ALICE Collaboration decay modes under study with two different analysis methods (standard cuts and MVA) it was possible to extract a signal in four p T intervals from 2 to 12 GeV/c.
The results from each decay mode and analysis method agree within statistical and systematic uncertainties. After averaging the results obtained from the different decay modes under study, the final result is compared with pQCD calculations and with the outcome of event generators. The Λ + c /D 0 baryon-tomeson ratio is discussed in Sec. 6.2, and the results in pp and p-Pb collisions are compared with previous measurements in different collision systems and at different centre-of-mass energies, and compared with expectations from Monte Carlo pp event generators. Finally, in Sec. 6.3 the nuclear modification factor R pPb is computed and compared with the results for D mesons and the predictions from models including cold-nuclear-matter and hot-medium effects.
6.1 Prompt Λ + c production cross section Figure 8 (left) shows the p T -differential cross section of prompt Λ + c baryons in |y| < 0.5 in pp collisions at √ s = 7 TeV as measured in the decay channels Λ + c → pK − π + , Λ + c → pK 0 S and Λ + c → e + ν e Λ (averaged with the corresponding charge conjugates). Figure 8 (right) shows the p T -differential cross section of prompt Λ + c in −0.96 < y < 0.04 in p-Pb collisions at √ s NN = 5.02 TeV in the decay channels Λ + c → pK − π + and Λ + c → pK 0 S . In this and following figures the marker is placed at the centre of the p T interval unless differently specified, the horizontal bar spans the width of the p T interval, the vertical error bar is the statistical uncertainty and the box is the systematic uncertainty.
For both collision systems, the cross sections measured from the different decay channels and analysis methods are compatible within statistical and uncorrelated systematic uncertainties, which include the uncertainty on the respective branching ratios. The largest discrepancy is observed in the p T interval 6 < p T < 8 GeV/c in pp collisions between the Λ + c → pK 0 S decay and the semileptonic decay channel, which differ by 1.7σ after adding in quadrature statistical and uncorrelated systematic uncertainties.
To obtain a more precise determination of the cross section in each collision system, these results were averaged together, taking into account the correlation between the statistical and systematic uncertainties. In the hadronic analyses (Λ + c → pK − π + and Λ + c → pK 0 S ), the sources of systematic uncertainty assumed to be uncorrelated between different decay channels are those due to the raw-yield extraction, the Λ c selection, and the PID efficiency. The sources assumed to be correlated are those due to the tracking efficiency, the generated p T shape of the Λ c in simulation, the beauty feed-down, and the luminosity. The branching ratio uncertainties were treated as partially correlated among the hadronic decay modes, as indicated in [12].
For the semileptonic analysis there are sources of systematic uncertainties that are correlated with other sources in the hadronic decay channel. In these cases, the systematic uncertainties were assumed to be fully correlated. The uncertainties due to the reconstruction of the electron and the Λ and the acceptance of the eΛ pair are assumed to be correlated with the tracking efficiency contribution in the hadronic decay modes. The uncertainties due to the generated p T shape of the Λ c in simulation are assumed to be correlated, as well as the contribution from the Λ b feed-down. Other sources, including the uncertainty coming from the cuts on the eΛ pairs, from the wrong-sign subtraction, the Ξ 0,+ c feed-down, the unfolding, the selections on the Λ decay topology, the electron identification, and the branching ratio are assumed to be fully uncorrelated between the results from the semileptonic and hadronic decay modes.
To average the different decay channels in pp collisions, where all measurements are statistically uncorrelated, the cross section from each decay channel was given a weight corresponding to the inverse of the quadratic sum of the relative statistical and uncorrelated systematic uncertainties, also taking into account the partial correlation in the branching ratios, following the approach in [87].
In the case of the analyses in p-Pb collisions, the cross sections in the two hadronic decay channels were measured with two different approaches, namely the standard cut method and the MVA method. A high degree of correlation exists between the analyses within the same decay channel, so the statistical uncertainty between analyses within the same decay channel was treated as fully correlated. The systematic uncertainty due to the yield extraction was assumed to be uncorrelated among different analyses, while all other sources of systematic uncertainty were treated as correlated. The statistically-correlated analyses are averaged using the relative uncorrelated systematic uncertainties as weights. Figure 9 shows the results of the p T -differential production cross section of prompt Λ + c baryons in pp and in p-Pb collisions obtained with the averaging procedure described above. In Fig. 9 (left) our measurement in pp collisions is compared with GM-VFNS perturbative QCD calculations [1,2] and with the results of the POWHEG event generator [9]. GM-VFNS has predictions for the Λ c baryon for p T > 3 GeV/c and the calculations were performed using CTEQ 6.6 [88] parameterisations of the PDFs, assuming the charm-quark mass m c = 1.5 GeV/c 2 , and with the fragmentation function and fractions tuned on e + e − collision data, which results in a fragmentation fraction value f (c → Λ c ) = 0.061 [89]. For the POWHEG calculation starting at p T = 1 GeV/c, the POWHEG-BOX package [90] was used for the NLO calculations and interfaced with PYTHIA 6.4.25 for the parton shower simulation and hadronisation. The POWHEG calculations were performed using CT10nlo [91] parameterisations of the PDF and m c = 1.5 GeV/c 2 . The uncertainties shown are the envelope of the predictions obtained varying the factorisation and renormalisation scales as proposed in [4]. The GM-VFNS predictions underestimate the measured cross section, which is on average higher by a factor 2.5 than the central value of the perturbative QCD calculation, as it can be seen in the bottom panel of Fig. 9. Moreover, POWHEG underpredicts the measured cross section by a factor of 18 (4) at low (high) p T . However both GM-VFNS and POWHEG describe the measured D-meson cross sections at central rapidities [7,92] and GM-VFNS describes the Λ c cross section at forward rapidities [52]. It is noted that the fragmentation functions used in these calculations were derived from e + e − collision data, and thus the underestimation of the data by GM-VFNS and POWHEG might hint at a violation of the universality of the fragmentation functions. This possibility is for example discussed in [93] considering data in the light flavour sector.
In Fig. 9 (right) the Λ c cross section in p-Pb collisions is compared with the cross section obtained with a calculation based on POWHEG using CT10nlo PDF with nuclear modification from EPS09NLO, scaled by the mass number of lead (A = 208). This calculation for p-Pb collisions underpredicts the measured Λ + c production in ALICE ALICE Collaboration values by a similar amount as observed in pp collisions. The Λ c cross section is also compared with a calculation [94], based on a data-driven modelling of the scattering at the partonic level, specifically designed to evaluate the impact of the nuclear modification of the gluon density on heavy-flavor hadrons. The tool is based on the HELAC-Onia package [95,96], originally developed for heavy-quarkonium studies, recently extended to heavy-flavor mesons and baryons. Differently from other calculations shown in Fig. 9, this is therefore a prediction for p-Pb collisions based on pp data. Specifically the authors constrained their parameterisation of the cross section to the LHCb measurements of Λ c production in pp collisions at √ s = 7 TeV and 2 < y < 4.5 [52] and they folded it with the nuclear modification of the PDFs from EPS09NLO. This model underpredicts our measurement by a factor two.

Λ +
c /D 0 baryon-to-meson ratio The Λ + c /D 0 production ratio is sensitive to hadronisation mechanisms in the charm sector. For the D 0 cross section we use the ALICE measurements [7,27]. The Λ + c /D 0 ratio is computed by integrating the p T -differential cross sections of Λ c and D 0 (both obtained as an average of particles and anti-particles) over their common p T interval, namely 1 < p T < 8 GeV/c for pp collisions and 2 < p T < 12 GeV/c for p-Pb collisions. In the integration, the systematic uncertainty due to the raw-yield extraction in the hadronic decay analyses Λ + c → pK − π + and Λ + c → pK 0 S were assumed to be fully uncorrelated between p T intervals, and the rest of the uncertainty sources were assumed to be fully correlated between p T intervals. In the Λ + c /D 0 ratio, the uncertainties due to the tracking efficiency, the beauty feed-down, and the luminosity were assumed to be fully correlated between the Λ + c and D 0 cross sections, and all other sources were assumed to be fully uncorrelated. The resulting baryon-to-meson ratio Λ + c /D 0 measured in Λ + c production in ALICE ALICE Collaboration pp collisions at √ s =7 TeV, |y| < 0.5, and 1 < p T < 8 GeV/c is In p-Pb collisions at √ s NN =5.02 TeV, −0.96 < y < 0.04, and 2 < p T < 12 GeV/c the measured baryonto-meson ratio is and is compatible within uncertainty with that measured in pp collisions. A list of existing measurements of the Λ + c /D 0 ratio in different collision systems and kinematic ranges is reported in Tab. 5. In Fig. 10, the measured Λ + c /D 0 ratio in pp and p-Pb collisions is presented as a function of p T (left panel) and rapidity (right panel) and compared with the LHCb measurement in pp collisions, with values derived by the LHCb Collaboration [97] from their published result [52].  For the measurements in e + e − and ep collisions reported in Tab. 5 and for the LHCb results reported in Fig. 10 the central values were multiplied by a correction factor that takes into account the most recent values of the BR of the Λ + c → pK − π + and D 0 → K − π + decays [12]. Wherever the systematic uncertainties for the branching ratios were quoted separately, they were updated according to the most recent values. Luminosity systematic uncertainties that cancel out in the ratio were not considered. The Λ + c /D 0 ratio was obtained, when available, from the ratio of the measured fragmentation fractions f (c → Λ c ) and f (c → D 0 ).
As shown in the table, a comparison is not straightforward given the different energy scales, the different collision systems and the fact that the extrapolation in the phase space down to p T = 0 was done for only a fraction of all measurements. The ratio Λ + c /D 0 can depend on the p T interval in which it is evaluated because of the possible differences in the fragmentation functions of charm quarks into baryons and mesons, which would result in different momentum distributions of Λ c baryons as compared to D 0 Λ + c production in ALICE ALICE Collaboration mesons. The results reported in this paper for the Λ + c /D 0 ratio are higher than previous measurements carried out in e + e − and ep collisions, and at lower centre-of-mass energies, where proposed mechanisms expected to enhance baryon production should play a negligible role as discussed in Sec. 1. In the beauty sector a difference in the fragmentation fraction f (b → Λ b ) has been reported, with larger values observed in pp and pp collisions, respectively at Tevatron [82] and at the LHC [83], with respect to e + e − collisions at LEP [80].
As shown in Fig. 10 the ratios measured by ALICE in pp and p-Pb collisions at mid-rapidity are compatible, both as a function of p T and p T -integrated, within uncertainties. The LHCb result in rapidity intervals suggests a decreasing trend towards mid-rapidity (influencing in turn the rapidity-averaged values reported in Fig. 10 (left)) that is not consistent with the ALICE result despite the large uncertainties. Such a trend is not reported by LHCb in their recent preliminary result in p-Pb collisions [53]. Although the ALICE result seems to decrease with increasing transverse momentum, a firm conclusion cannot be drawn as to whether the observed difference between the Λ + c /D 0 ratios at forward and mid-rapidity is significantly p T -dependent. Figure 11 compares the Λ + c /D 0 ratio as a function of p T (left panel) and rapidity (right panel) in pp and p-Pb collisions with predictions obtained from Monte Carlo pp event generators, namely PYTHIA8 with Monash tune and with another tune [17] that includes a model of string formation beyond the leadingcolour approximation, DIPSY with rope parameters taken from [18], and HERWIG7 which uses a cluster hadronisation mechanism. As for the cross section calculations described in Sec. 6.1, fragmentation parameters for these predictions are derived from e + e − collision data. The enhanced colour reconnection mechanisms enabled in PYTHIA8 increase the baryon-to-meson ratio in the charm sector, bringing the prediction closer to the data at mid-rapidity. The DIPSY generator with a rope configuration, which is expected to increase the baryon-to-meson ratio, instead predicts values similar to those from PYTHIA8 with Monash tune, which are lower than the values in e + e − and ep collisions as reported in Tab. 5. Similar predictions were obtained with HERWIG7. The p-Pb measurement is compared then in Fig. 11 (left) with the calculations from Lansberg and Shao [94] for p-Pb, with Λ + c and D 0 cross section obtained through a parameterisation of pp data and using EPS09NLO nuclear modification factors. Among the different predictions this calculation is the closest to data. Finally, all models predict a flat rapidity dependence which does not describe the trend observed at forward rapidity. The preliminary result from LHCb in p-Pb collisions [53] also shows a flat rapidity dependence in the 1.5 < y < 4 interval.  Figure 11: The Λ + c /D 0 ratio measured in pp and p-Pb collisions by ALICE as a function of p T (left) and as a function of y for 2 < p T < 8 GeV/c (right). The measurements from pp collisions are compared with different event generators (quoted tunes for PYTHIA and DIPSY taken respectively from [17] and [18]). The p-Pb measurement as a function of p T is also compared with calculations from Lansberg and Shao [94]. The predictions from event generators as a function of y are also compared with the LHCb measurement [52,97].
6.3 Λ c -baryon nuclear modification factor in p-Pb collisions at √ s NN = 5.02 TeV The nuclear modification factor R pPb of Λ c baryons was calculated from the results presented in Sec. 6.1 by dividing the p T -differential prompt production cross section in p-Pb collisions at √ s NN = 5.02 TeV by that in pp collisions corrected for the different centre-of-mass energy and rapidity coverage of the pp and p-Pb measurements and multiplied by the mass number A = 208.
In particular, the cross section in pp collisions measured at √ s = 7 TeV and |y| < 0.5 was scaled in each p T interval to √ s = 5.02 TeV and −0.96 < y < 0.04 using a factor f √ s,y FONLL calculated with FONLL perturbative QCD calculations [4], following a similar procedure to the D-meson R pPb measurement [27]: with FONLL cross sections calculated at 7 TeV in |y| < 0.5, and at 5.02 TeV in −0.96 < y < 0.04. The uncertainties on the scaling factor are calculated by consistently varying the charm-quark mass, the PDF, and the factorisation and renormalisation scales in the calculations at the two energies.
The fragmentation function of charm quarks into Λ c baryons is not well known. However, it has been verified that changing the fragmentation function does not change the scaling factor significantly: the f √ s,y FONLL values obtained from FONLL calculations for D 0 , D + and D * + vary by less than 1%. For this reason, the D + production cross section ratio from FONLL was chosen for the central values of f √ s,y FONLL (p T ), and the uncertainty was estimated by varying the fragmentation function. The bare c-quark cross section from FONLL defines the upper uncertainty of the scaling factor, as the "hardest" fragmentation case, where it is assumed that all the momentum of the c quark is carried by the Λ c . The c-quark cross section from FONLL, convolved with a fragmentation function modelled using the Peterson parameterisation [99] with ε = 0.1, defines the lower uncertainty of the scaling factor as the "softest" case. For both limits, the associated uncertainties from FONLL were included. These two scenarios were chosen to encompass the values reported by the PDG review for charm-and beauty-quark fragmentation for different models of hard radiation [12]. It has also been verified that the Λ + c / D 0 ratio obtained using Λ + c production in ALICE ALICE Collaboration c R pPb from the POWHEG event generator [9] with EPS09NLO parameterisation of the nuclear modification of the PDFs [100] and the charm-hadron R pPb from the POWLANG transport model [37] assuming a QGP is formed in p-Pb collisions. these fragmentation scenarios for the Λ + c and the D 0 cross section from FONLL is compatible with the measured Λ + c /D 0 ratio. The uncertainty on the scaling factor varies from +13 −5 % in the p T interval 2-4 GeV/c to +6 −4 % in the p T interval 6-8 GeV/c. For the propagation of the uncertainties in the R pPb computation, the beauty feed-down uncertainties are considered fully correlated between the pp and p-Pb cross sections and the branching ratio uncertainties are considered partially correlated due to the different decay modes considered in the two collision systems, while all the other systematic uncertainties are treated as uncorrelated. The uncertainty due to the √ s and rapidity scaling of the pp reference was added in quadrature to the aforementioned sources. The luminosity uncertainties were treated as fully uncorrelated. Figure 12 (left) shows the Λ c -baryon nuclear modification factor R pPb in the range 2 < p T < 8 GeV/c. The result is compatible with unity within the large statistical and systematic uncertainties, and is consistent with the D-meson R pPb [27], which is shown in the same figure. Predictions for the R pPb for Λ c baryons from the POWHEG event generator with PYTHIA parton shower [9] and EPS09NLO parameterisation of nuclear modification of the PDFs [100] are presented in the right panel of Fig. 12. In the same panel, the calculations for the charmed-hadron nuclear modification factor from the POWLANG model [37], which assumes that also in p-Pb collisions at LHC energies a hot deconfined medium is formed, is superimposed. The POWLANG model utilises the Langevin approach to compute the transport of heavy quarks through an expanding QGP described by relativistic viscous hydrodynamics, but it does not include any specific mechanism to modify hadronisation, such as coalescence, that could lead to a baryon enhancement. This transport model predicts a deviation of R pPb from unity which is about 20-40% at low and intermediate momentum (p T < 5 GeV/c). The precision achieved with the current measurement does not allow us to distinguish between calculations with and without hot medium effects. Λ + c production in ALICE ALICE Collaboration

Conclusions
We measured the Λ c baryon production in pp and p-Pb collisions with ALICE at the LHC using different decay channels and different analysis methods. In pp collisions, we reported the production cross section measurement at mid-rapidity (|y| < 0.5) and √ s = 7 TeV for this baryon, while in p-Pb collisions the Λ c production cross section was measured at the centre-of-mass energy per nucleon-nucleon pair √ s NN = 5.02 TeV in the centre-of-mass rapidity interval −0.96 < y < 0.04. The results were reported for pp collisions in the transverse-momentum interval 1 < p T < 8 GeV/c and for p-Pb collisions in 2 < p T < 12 GeV/c.
The measurement of the Λ c baryon, due to its short lifetime, is challenging: the p T -differential production cross sections were therefore obtained averaging the results from different decay channels (purely hadronic and semileptonic) and with different analysis approaches, using standard cuts, Multivariate Analysis techniques and a dedicated procedure to subtract background pairs for the semileptonic channel. Different PID-discriminating variables were also used. The results of all the analyses were found to be mutually consistent within uncertainties.
In the p T interval where calculations from the GM-VFNS perturbative QCD framework are available (3 < p T < 8 GeV/c), the predictions underestimate the measured cross section on average by a factor of 2.5. The comparison is, however, affected by large uncertainties, in particular from the theoretical estimates. Calculations based on POWHEG (available for p T > 1 GeV/c) with hadronisation from the PYTHIA parton shower, underpredict the measured values by a factor of 18 (4) at low (high) p T . A similar pattern is observed comparing cross section predictions obtained with POWHEG with measured values in p-Pb collisions. Calculations for this collision system based on a parameterisation of existing pp measurements for Λ c are closer to the data, even if they are still underpredicting measured values by a factor of 2.
The baryon-to-meson ratios Λ + c /D 0 measured in pp and p-Pb collisions are compatible within their statistical and systematic uncertainties. Our result in pp collisions (Λ + c /D 0 = 0.543 ± 0.061 ± 0.160 for 1 < p T < 8 GeV/c at mid-rapidity) is larger than previous measurements obtained at lower centre-ofmass energies and in different collision systems, and also higher than the results reported by LHCb at 2 < y < 4.5 for 2 < p T < 8 GeV/c in pp collisions at the same centre-of-mass energy.
We also compared the measured Λ + c /D 0 ratio to pp event generators that implement different hadronisation schemes. All underpredict the measured values: a better qualitative agreement with our results is obtained with PYTHIA tunes that include string formation beyond the leading-colour approximation, while significantly lower values are obtained with DIPSY and HERWIG7.
Finally, a first measurement of the nuclear modification factor R pPb of Λ c baryons was obtained and it was found to be compatible with unity in the transverse-momentum interval 2 < p T < 8 GeV/c, as well as with the R pPb of D mesons. The current precision of the measurement cannot constrain existing models.
When considered in their entirety, these results provide input for theoretical models based on pQCD calculations, event generators applying different hadronisation approaches and models describing CNM and/or hot-medium effects in proton-nucleus collisions. A better precision is expected to be reached with data presently collected during the LHC RUN 2, reducing in particular the statistical uncertainties, and, in the future, during the LHC RUN 3 and 4 following a major upgrade of the ALICE apparatus [101]. This set of measurements provides an initial reference for future investigation of Λ c production in Pb-Pb collisions where the interaction of charm quarks with the hot medium may affect its production.
Λ + c production in ALICE ALICE Collaboration

Acknowledgements
The ALICE Collaboration would like to thank M. Corradi and L. Gladilin from ZEUS Collaboration for providing guidance on how to treat the different sources of uncertainties in the computation of the Λ + c /D 0 ratio for the ZEUS data, P. Spradlin and V. Vagnoni from LHCb Collaboration for providing Λ + c /D 0 baryon-to-meson ratio shown in Fig. 10 in this paper, and P. Robbe for his help to configure the EVTGEN package.
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: