Study of the production of Λ 0 b and B 0 hadrons in pp collisions and ﬁrst measurement of the Λ 0 b → J/ψ pK − branching fraction

The product of the Λ 0 b ( B 0 ) diﬀerential production cross-section and the branching fraction of the decay Λ 0 b → J/ψ pK − ( B 0 → J/ψ K ∗ (892) 0 ) is measured as a function of the beauty hadron transverse momentum, p T , and rapidity, y . The kinematic region of the measurements is p T < 20 GeV /c and 2 . 0 < y < 4 . 5. The measurements use a data sample corresponding to an integrated luminosity of 3 fb − 1 collected by the LHCb detector in pp collisions at centre-of-mass energies √ s = 7 TeV in 2011 and √ s = 8 TeV in 2012. Based on previous LHCb results of the fragmentation fraction ratio, f Λ 0 b /f d , the branching fraction of the decay Λ 0 b → J/ψ pK − is measured to be where the ﬁrst uncertainty is statistical, the second is systematic, the third is due to the uncertainty on the branching fraction of the decay B 0 → J/ψ K ∗ (892) 0 , and the fourth is due to the knowledge of f Λ 0 b /f d . The sum of the asymmetries in the production and decay between Λ 0 b and Λ 0 b is also measured as a function of p T and y . The previously published branching fraction of Λ 0 b → J/ψ pπ − , relative to that of Λ 0 b → J/ψ pK − , is updated. The branching fractions of Λ 0 b → P + c ( → J/ψ p ) K − are determined.


Introduction
In quantum chromodynamics (QCD) the production process of b hadrons can be divided into two steps, assuming factorisation: a hard process for b production and a soft process to describe hadronisation. The hard process can be predicted by perturbative calculations in QCD; the soft process is parameterised by the fragmentation function, which has large uncertainties due to non-perturbative QCD contributions. The study of the production of b hadrons tests the factorisation ansatz. The ground state of the b-baryon family, Λ 0 b , has a wide range of decay modes. The study of its production and decays can offer complementary information to that obtained from the study of B mesons. The kinematic dependence of the production of Λ 0 b baryons relative to that of B mesons can test differences in the b quark hadronisation process between the two [1,2]. Furthermore, the asymmetry of heavy flavoured baryons and antibaryons produced in pp collisions is an important input for various asymmetry measurements. Leading-order QCD calculations predict equal production cross-sections for heavy baryons and heavy anti-baryons, while measurements at the ISR showed that Λ + c production is favoured in pp collisions at forward rapidity, y [3,4]. The CMS experiment measured the Λ 0 b and Λ 0 b production ratio in pp collisions at 7 TeV, and no asymmetry was observed, but the large uncertainties preclude definitive conclusions [5]. Measurements at LHCb can provide further tests of existing mechanisms, e.g., the string drag effect or the leading quark effect [6].
Measurements of Λ 0 b production to date have mostly been based on semileptonic decays and the hadronic decays Λ 0 b → J/ψ Λ and Λ 0 b → Λ + c π − (charge-conjugation is implied throughout the paper unless otherwise specified). Using semileptonic decays, the LHCb experiment measured the ratio of Λ 0 b baryon production to light B meson production, [7]. The kinematic dependence of the ratio of Λ 0 b to B 0 production, f Λ 0 b /f d , was measured using Λ 0 b → Λ + c π − and B 0 → D + π − decays, and the absolute branching fraction B(Λ 0 b → Λ + c π − ) was determined [8]. In this paper, the Λ 0 b candidates are reconstructed in the decay channel Λ 0 b → J/ψ pK − , which was first observed by LHCb in 2013 [9]. Compared with the open-charm decays of Λ 0 b baryons, this channel has higher trigger efficiencies, especially in the region of low transverse momentum, p T . Two pentaquark-charmonium states P c (4380) + and P c (4450) + were observed by LHCb [10] in the amplitude analysis of the Λ 0 b → J/ψ pK − decay. The measurement of the absolute branching fraction of Λ 0 b → J/ψ pK − in the current paper allows the pentaquark branching fractions to be determined. Other Λ 0 b decays with a charmonium meson in the final state, such as the Cabibbo-suppressed decay Λ 0 b → J/ψ pπ − [11], can use the Λ 0 b → J/ψ pK − decay as a reference to measure their absolute branching fractions.
The product of the Λ 0 b (B 0 ) differential production cross-section and the branching fraction of the Λ 0 b → J/ψ pK − (B 0 → J/ψ K * 0 ) decay is measured as a function of p T and y, where K * 0 indicates the K * (892) 0 meson throughout the text. The kinematic region of these measurements is p T < 20 GeV/c and 2.0 < y < 4.5 for the b hadron. The production ratio of the two b hadrons, defined as is determined, taking advantage of the cancellation of some uncertainties in both experimental measurements and theoretical calculations. Here, σ(Λ 0 b ) and σ(B 0 ) represent the production cross-sections of Λ 0 b and B 0 hadrons in pp collisions. The branching fraction B(Λ 0 b → J/ψ pK − ) is calculated from this result using previous measurements [7,8] and B(B 0 → J/ψ K * 0 ) [12]. The kinematic dependence of the sum of the asymmetries in the production and decay, a p+d ≡ a prod + a decay , between Λ 0 b and Λ 0 b is studied using Λ 0 b → J/ψ pK − and Λ 0 b → J/ψ pK + decays. Furthermore, using the measurement of B(Λ 0 b → J/ψ pK − ), the branching fractions of the decays Λ 0 b → J/ψ pπ − and Λ 0 b → P + c (→ J/ψ p)K − are determined. The measurements in this paper are based on a data sample corresponding to an integrated luminosity of 3 fb −1 , collected by the LHCb experiment in pp collisions at centreof-mass energies √ s = 7 TeV in 2011 and √ s = 8 TeV in 2012. Separate measurements are performed for each of the two centre-of-mass energies.

Detector and simulation
The LHCb detector [13,14] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region [15], a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes [16] placed downstream of the magnet.
The tracking system provides a measurement of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of (15 + 29/p T ) µm, where p T is the component of the momentum transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors [17]. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillatingpad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [18]. The online event selection is performed by a trigger [19], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction. In the hardware trigger, events are selected by requiring at least one high-p T track that is consistent with a muon hypothesis. In the software trigger, two well-reconstructed muons are required to form a vertex with good fit χ 2 and to have an invariant mass consistent with that of the J/ψ meson [20]. The trigger also requires a significant displacement between the J/ψ vertex and the associated PV of the pp collision.
In the simulation, pp collisions are generated using Pythia [21,22] with a specific LHCb configuration [23]. Decays of hadronic particles are described by EvtGen [24], in which final-state radiation is generated using Photos [25]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [26,27] as described in Ref. [28]. The physics models used by LHCb in Geant4 for hadronic interactions have been tested against experimental data from COMPAS [20], and good agreement was observed. 1

Event selection
where the J/ψ mesons are reconstructed in the dimuon final state, and K * 0 candidates are reconstructed from K * 0 → K − π + decays. Since the Λ 0 b → J/ψ pK − and B 0 → J/ψ K * 0 decays have the same topology, a similar event selection is adopted for both.
An offline selection is applied after the trigger and is divided into two steps: a preselection and a multivariate selection based on a boosted decision tree (BDT) [29][30][31][32]. In the preselection, each track of the Λ 0 b (B 0 ) candidate is required to be of good quality [14,[33][34][35]. Identified muons are required to have p T greater than 550 MeV/c, while hadrons are required to have p T greater than 250 MeV/c. The muons should be inconsistent with originating from any PV, as determined by their impact parameter. Each J/ψ candidate is required to have a good vertex fit χ 2 and an invariant mass within +43 −48 MeV/c 2 of the known J/ψ mass [20]. Particle identification (PID) requirements are imposed on the final-state tracks. For the kaon and proton in the Λ 0 b → J/ψ pK − decay, the sum of the kaon and proton p T should be larger than 1 GeV/c. Each K * 0 candidate is required to have a good vertex fit χ 2 and to have p T greater than 1 GeV/c. The invariant mass of the reconstructed K * 0 is required to be within ±70 MeV/c 2 of the K * 0 mass [20]. Each b hadron candidate must have a good vertex fit χ 2 , be consistent with originating from the PV, and have a decay time greater than 0.2 ps.
Some non-combinatorial backgrounds exist in the Λ 0 b → J/ψ pK − data sample, originating from B 0 → J/ψ K − π + and B 0 s → J/ψ K − K + decays with the π + and K + misidentified as a proton. In order to suppress these events, the invariant mass is recalculated by interpreting the proton candidate as a pion or a kaon, and the two relevant invariant mass regions are vetoed: m(J/ψ K − π + ) ∈ [5250, 5310] MeV/c 2 and m(J/ψ K + K − ) ∈ [5340, 5400] MeV/c 2 . After the mass vetoes these background contributions are reduced to a negligible level.
After the preselection, the Λ 0 b → J/ψ pK − (B 0 → J/ψ K * 0 ) candidates are filtered with the BDT to further suppress combinatorial background. For the decays Λ 0 b → J/ψ pK − and B 0 → J/ψ K * 0 , the same BDT classifier is applied. Independent BDT classifiers are used for the 2011 and 2012 samples. In the BDT training a simulated Λ 0 b sample is used as the signal. The background is taken from the lower, (5420, 5560) MeV/c 2 , and upper, (5680, 5820) MeV/c 2 , sidebands of the Λ 0 b invariant mass distribution in data. Events in the sidebands are randomly divided into two parts, one for the training and the other for the test. No overtraining is observed. The following information is used by the BDT classifier: the kinematic properties and the impact parameters of the tracks; and the vertex quality, the decay length and the impact parameter of the reconstructed b hadron candidate. The variables used for the training are chosen based on their power to discriminate signal from background and on the similarity of their distributions for Λ 0 b → J/ψ pK − and B 0 → J/ψ K * 0 decays. The threshold for the BDT response is chosen to maximise S/ √ S + B, where B represents the number of background events estimated from the sideband region and S the number of signal events in the mass peak.

Cross-section and branching fraction determination
The product of the differential production cross-section of each b hadron and the corresponding branching fraction is calculated as where N (p T , y) and ε(p T , y) are respectively the signal yield and the efficiency as functions of p T and y of the b hadron, ∆p T and ∆y are the bin widths, L is the integrated luminosity, B is the absolute branching fraction of the Λ 0 b → J/ψ pK − (B 0 → J/ψ K * 0 ) decay, and B inter represents the branching fractions of the intermediate decays: The luminosity is measured with van der Meer scans and a beam-gas imaging method [36]. The 2011 and 2012 data samples correspond to 1019 ± 17 pb −1 and 2056 ± 23 pb −1 , respectively. The branching fraction B(J/ψ → µ + µ − ) = (5.961 ± 0.033)% [20], while B(K * 0 → K − π + ) is taken to be 2/3 assuming isospin symmetry. The branching fraction B(B 0 → J/ψ K * 0 ) = (1.29 ± 0.05 ± 0.13) × 10 −3 as measured by Belle [12] is used in preference to the world average value, since in the Belle result the S-wave component is subtracted.

Signal determination
The signal yields of the Λ 0 b → J/ψ pK − and B 0 → J/ψ K * 0 decays are determined from unbinned extended maximum likelihood fits to the invariant mass distributions of the reconstructed b hadron candidates in each p T and y bin. In order to improve the mass resolution, the b hadron is refitted with constraints [37] that it originates from the PV and that the reconstructed J/ψ mass equals its known mass [20]. Figure 1 shows, as an example of one of the fit results, the invariant mass distributions of Λ 0 b and B 0 candidates in the kinematic region p T ∈ [6,7] GeV/c and y ∈ [3.0, 3.5] for the 2012 data sample. The signal shape in the fits is modelled by a double-sided Crystal Ball (DSCB) function, an empirical function comprising a Gaussian core together with power-law tails on both sides. The mean and the mass resolution of the DSCB function are free in the fits, while the tail parameters are determined from simulation in each kinematic bin according to the empirical function given in Ref. [38]. The combinatorial background is modelled by an exponential function whose parameters are left free in the fits.
In the fits to the Λ 0 b mass distribution, there is a contribution from the Λ 0 b → J/ψ pK − decay in which the proton is misidentified as a kaon and the kaon is misidentified as a proton. This background is denoted as the doubly misidentified background, and it is modelled by a DSCB function. All parameters of this DSCB function are fixed from the simulation study, including: the difference between the mean of this DSCB function and that of the signal shape; the ratio of the mass resolution between these two DSCB; the yield fraction relative to the Λ 0 b signal channel; and the tail parameters. In the B 0 → J/ψ K * 0 sample, in addition to the combinatorial background, there are two further sources of background. One is the decay B 0 s → J/ψ K * 0 , which populates the upper sideband of the invariant mass distribution, and is modelled with a DSCB function. The tail parameters of this DSCB function are the same as those of the B 0 signal shape and the remaining parameters are free in the fits. The other comes from partially reconstructed B mesons and is described by the tail of a Gaussian function. The associated mean and width are free parameters in the fits.
According to a previous LHCb measurement [39], the fraction of the K * 0 meson contribution in the B 0 → J/ψ K * 0 decay is calculated as (89.9 ± 0.4 ± 1.3)%, where the remainder is due to the S-wave component in the K − π + system and its interference with the K * 0 meson. The fitted B 0 yields are subtracted with this number to remove the components from S-wave and its interference.

Efficiencies
The efficiency ε Λ 0 b ,B 0 (p T , y) consists of the geometrical acceptance of the detector, the trigger efficiency, the reconstruction and preselection efficiency, the hadron PID efficiency, and the BDT selection efficiency. All the efficiencies are determined from a sample of simulated signal events, except the hadron PID efficiency, which is determined from data with tracks from the decays [10] is not modelled in the simulation. The simulated sample is weighted to reproduce the distributions of the BDT training variables and the two-dimensional distribution of m(pK − ) and m(J/ψ p) observed in the background-subtracted data sample, which has been obtained using the sPlot technique [40], with the b-hadron invariant mass as the discriminating variable. It is found that the correlations between the discriminating variable and the control variables are negligible.

Asymmetry determination
The observed (raw) asymmetry for Λ 0 b and Λ 0 b is defined as The symbol N (x) is the signal yield in the given bin of x from the fits to the invariant mass distribution of the Λ 0 b (Λ 0 b ) sample, where x denotes p T or y. The observed asymmetry is a sum of several contributions: the asymmetry between the numbers of the produced Λ 0 b and Λ 0 b baryons in pp collisions, a prod (x); the decay asymmetry between the Λ 0 b → J/ψ pK − and Λ 0 b → J/ψ pK + channels, a decay (x); the asymmetry between the p and p detection efficiencies, a p D (x); the asymmetry between the K − and K + detection efficiencies, a K D (x); and the asymmetry between the PID efficiencies for Λ 0 b and Λ 0 b baryons, a PID (x). Other possible asymmetries are neglected. Assuming that all these asymmetries are small, the asymmetries a p+d (x) of Λ 0 b and Λ 0 b baryons can be calculated as The value A raw (x) − a PID (x) can be calculated as where ε The kaon detection asymmetry a K D (x) as a function of p T and y is obtained from a previous LHCb study [41].
The proton detection asymmetry a p D (x) as a function of p T or y is estimated from simulation, which uses the Geant4 model as described in Section 2. The proton detection asymmetry as a function of p T or y is calculated with the proton and antiproton track reconstruction efficiencies in the corresponding kinematic bin. It is checked that the kinematic distributions of protons and Λ 0 b baryons in the simulation sample are consistent with those in the data sample. As a crosscheck, the proton detection asymmetries are also estimated through a simulation sample, where the Λ 0 b signals are partially reconstructed without using the proton information, and the results are consistent.

Systematic uncertainties
Several sources of systematic uncertainties are studied in the analysis and are summarised in Tables 1 and 2. For the production cross-section measurements, the uncertainties originate from the determination of the signal yields, efficiencies, branching fractions and luminosities. The total systematic uncertainties are obtained from the sum in quadrature of all components.
Imperfect knowledge of the mass distributions for the signal and backgrounds causes systematic uncertainties in the signal yield determination. For the signal shape, the Apollonios function [42] and the sum of a Gaussian function and a Crystal Ball function  The fits are repeated with a linear function substituted for the exponential model for the combinatorial background. The fits are also repeated without the double misidentified components. The maximum differences of the signal yields from the nominal results are quoted as systematic uncertainties due to the background shape.
Most efficiencies are estimated from simulation. The limited size of the simulation sample leads to systematic uncertainties on the efficiencies ranging from 1.7% to 16.5%.
The tracking efficiency is estimated from simulation and calibrated by data [43]. The uncertainty of the calibration is 0.4% per track. Additional systematic uncertainties are assigned to hadrons due to imperfect knowledge of hadron interactions in the detector, 1.1% for kaons, 1.4% for pions and 1.4% for protons.
The BDT efficiency is estimated with the weighted simulation sample to ensure that the distributions of the two training variables, the kinematic properties of the tracks and the vertex quality, agree with those in data. The uncertainties on the weights are propagated to the final results to give the corresponding systematic uncertainty.
The trigger efficiency is determined in the simulation and validated in a control sample of J/ψ → µ + µ − decays [19]. The difference of the central values of this determination in data and the simulation in each bin is taken as the systematic uncertainty. Uncertainties due to the limited sample size of the simulation are added in quadrature.
The PID efficiency is estimated with a data-driven method. A sample of J/ψ → µ + µ − , D * + → D 0 (→ K − π + )π + and Λ + c → pK − π + decays obtained without using PID information is used to evaluate the PID efficiency. The limited sample size used to calculate the PID efficiency introduces a systematic uncertainty in each kinematic bin. To study the bin-by-bin migration effect, the number of the bins is doubled or halved and the PID efficiency is recalculated. The largest deviation from the nominal result is taken as the uncertainty.
To account for the rich and complex structure of multiple intermediate resonances in the Λ 0 b → J/ψ pK − decay, the simulation sample is weighted in two-dimensional bins of m(K − p) and m(J/ψ p) to match the data. Pseudoexperiments are performed to estimate the systematic uncertainties due to the weights. The weight in each bin is varied according to its uncertainty and the total efficiency is recalculated. The RMS of the distribution obtained from the pseudoexperiments is taken as the systematic uncertainty. As mentioned in Section 3, the preselection includes mass vetoes. The preselection efficiencies are estimated from the simulation sample. A fit to the Λ 0 b invariant mass distribution in the vetoed data sample is performed, which gives the number of Λ 0 b → J/ψ pK − signal events rejected by the vetoes. The fraction of the vetoed signal events in the data sample is compared with that in the simulation sample. A difference of 1.3% (1.9%) is observed for the 2011 (2012) sample, and this is taken as the systematic uncertainty.
The uncertainty in the determination of the integrated luminosity is 1.7% (1.2%) for the 2011 (2012) data sample [36]. An uncertainty of 0.6% is taken on B(J/ψ → µ + µ − ) [20]. The fractions of the S-wave component in the K − π + system and their interference were determined by a previous LHCb measurement, and their 1.4% uncertainty [39] is taken as a systematic uncertainty for the B 0 → J/ψ K * 0 decay.
In the Λ 0 b and Λ 0 b asymmetry measurement, all of the uncertainties mentioned above cancel in the ratio, except for those due to the signal shape, the background shape, the limited sample size and the PID efficiency. Since a data-driven determination of proton detection asymmetries is not available, the difference in the determination of the kaon detection asymmetries in data and simulation is taken as a systematic uncertainty for the proton detection asymmetry. The uncertainties vary from 0.1% to 1.9% in kinematic bins, with large values occurring in bins of low p T or low signal yields. In the LHCb Geant4 physics models, the cross-sections of interactions between particles and the material are checked with test beam data as discussed in Section 2. There are more data for protons than for kaons. Therefore, these uncertainties can be considered to be conservative.

Cross-section results
The product of the Λ 0 b (B 0 ) double-differential cross-section and the branching fraction of the decay Λ 0 b → J/ψ pK − (B 0 → J/ψ K * 0 ) is shown in Fig. 2, and the values are listed in Tables 3, 4, 5 and 6 in the Appendix. By integrating over y or p T , the single differential production cross-sections, shown in Fig. 3, are obtained. Figure 4 shows the p T distribution of the Λ 0 b production, fitted by a power-law function with the Tsallis parameterisation [44,45]: dσ p T dp T ∝ 1 where T is a temperature-like parameter, N determines the power-law behaviour at large E k⊥ , and E k⊥ ≡ p  other and with the values found by CMS [5]. Other functions suggested in Ref. [46] do not give acceptable fits to the data. In Fig. 4 the data points are placed in the bin according to the prescription of Ref. [47].
The integrated cross-sections of the b hadrons with 0 < p T < 20 GeV/c and 2.0 < y < 4.5 are measured to be Taking the branching fraction B(B 0 → J/ψ K * 0 ) from Belle [12], the measured B 0 production cross-section at 7 TeV is consistent with the previous LHCb measurement [38]. The where the first uncertainties are statistical and the second systematic. Many systematic uncertainties cancel totally or partially in these ratios: the ratio of the luminosities is known with a precision of 1.44% [36]; the tracking efficiency is considered to be fully correlated, due to the fact that the associated systematic uncertainty is dominated by hadronic interactions of the tracks in the detector; the mass veto efficiency, the branching fraction of the J/ψ → µ + µ − decay, and the S-wave contribution in the K − π + system are also fully correlated. All other sources are considered uncorrelated. The ratio of the integrated production cross-sections agrees with FONLL predictions [48][49][50]. Figure 5 shows the p T and y dependence of the ratios for Λ 0 b and B 0 production cross-sections at 8 TeV with respect to those at 7 TeV, together with linear fits to the distributions: where p T = 6.7 (6.9) GeV/c is the mean p T of Λ 0 b (B 0 ) hadrons in the data sample, y = 3.1 is the mean y, and ndf is the number of degrees of freedom. The p T dependence of the ratio agrees with FONLL predictions, while the y dependence does not.
The measured values of the ratio R Λ 0 b /B 0 , defined in Eq. 1, as a function of p T and y are shown in Fig. 6. In the region p T < 5 GeV/c, no p T dependence of the ratio R Λ 0 b /B 0 is observed, while the ratio decreases for p T > 5 GeV/c. No dependence with rapidity is observed. In Fig. 6 the p T dependence of the ratio R Λ 0 b /B 0 is fitted with the fragmentation function ratio f Λ 0 b /f d (p T ) given in Ref. [8], which is only defined in the range p T > 3 GeV/c. The asymmetry a p+d between Λ 0 b and Λ 0 b is shown in Fig. 7. The values are listed in Table 7 in the Appendix. The data points are fitted with linear functions. The slope fitted to the asymmetry as a function of p T is consistent with zero, (2.3 ± 3.0) × 10 −3 /(GeV/c) for 7 TeV and (3.5 ± 2.0) × 10 −3 /(GeV/c) for 8 TeV. The fit to a p+d (y) gives a non-zero slope, and a combination of the results for 7 TeV and 8 TeV gives a p+d (y) = (−0.001 ± 0.007) + (0.058 ± 0.014)(y − y ), where y = 3.1 is the average rapidity of Λ 0 b hadrons in the data sample. The non-zero slope suggests some baryon number transport from the beam particles to the less centrally produced Λ 0 b , which leads to a Λ 0 b /Λ 0 b cross-section ratio that increases with rapidity and which can be interpreted as, for example, a string drag effect or leading quark effect [6,51].

Branching fraction results
The ratio R Λ 0 b /B 0 can be calculated in bins of p T as: It is related to the fragmentation fraction ratio where S ≡ B(Λ 0 b → J/ψ pK − )/B(B 0 → J/ψ K * 0 ) is a constant factor, which can be determined from the fit in Fig. 6. The absolute branching fraction of the decay Λ 0 b → J/ψ pK − can then be measured as  The first uncertainty is statistical, the second is systematic, the third is due to the uncertainty on the branching fraction of the B 0 → J/ψ K * 0 decay, and the fourth is due to the knowledge of f Λ 0 b /f d . In Ref. [11] the ratio B(Λ 0 b → J/ψ pπ − )/B(Λ 0 b → J/ψ pK − ) was reported. Combining this with the value of B(Λ 0 b → J/ψ pK − ) above, the branching fraction of Λ 0 b → J/ψ pπ − is determined as where the first uncertainty is statistical, the second is due to the systematic uncertainty on B(Λ 0 b → J/ψ pπ − )/B(Λ 0 b → J/ψ pK − ), and the third is due to systematic uncertainty on B(Λ 0 b → J/ψ pK − ). Two pentaquark-charmonium states, P c (4380) + and P c (4450) + , were observed by LHCb in the amplitude analysis of the Λ 0 b → J/ψ pK − decay [10], and the fractions f (P + c ) of the two pentaquark-charmonium states in the Λ 0 b → J/ψ pK − decay were measured. Using these fractions and the value of B(Λ 0 b → J/ψ pK − ) obtained in this analysis, the branching fractions where the first uncertainty is statistical, the second is due to the systematic uncertainty on f (P + c ), and the third is due to the systematic uncertainty on B(Λ 0 b → J/ψ pK − ).

Conclusion
Using a data sample corresponding to an integrated luminosity of 3 fb −1 collected by the LHCb detector in 2011 and 2012, the product of the Λ 0 b differential production cross-section and the branching fraction of the decay Λ 0 b → J/ψ pK − is measured as a function of the Λ 0 b baryon's transverse momentum and rapidity. The product of the B 0 differential production cross-section and the branching fraction of the decay B 0 → J/ψ K * 0 is also measured. The kinematic region of the measurements is p T < 20 GeV/c and 2.0 < y < 4.5.