First measurement of $\Xi_{\rm c}^0$ production in pp collisions at $\mathbf{\sqrt{s}}$ = 7 TeV

The production of the charm-strange baryon $\Xi_{\rm c}^0$ is measured for the first time at the LHC via its semileptonic decay into e$^+\Xi^-\nu_{\rm e}$ in pp collisions at $\sqrt{s}=7$ TeV with the ALICE detector. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval 1 $<$ $p_{\rm T}$ $<$ 8 GeV/$c$ at mid-rapidity, $|y|$ $<$ 0.5. The transverse momentum dependence of the $\Xi_{\rm c}^0$ baryon production relative to the D$^0$ meson production is compared to predictions of event generators with various tunes of the hadronisation mechanism, which are found to underestimate the measured cross-section ratio.

In this paper, we report the first measurement of the p T -differential production cross section of Ξ 0 c multiplied by the branching ratio (BR) into the semileptonic decay mode, Ξ 0 c → e + Ξ − ν e , and its ratio to the measured production cross section of D 0 mesons [21] as a function of p T , up to 8 GeV/c.The absolute branching ratio of this Ξ 0 c decay is currently unknown [34].Using a data sample of pp collisions at √ s = 7 TeV recorded with the ALICE detector in 2010, the measurement is performed by analyzing e + Ξ − pairs formed by combining positrons and Ξ − baryons reconstructed with the detectors of the ALICE central barrel, covering the pseudorapidity interval |η| < 0.9.The missing momentum of the neutrino is corrected using unfolding techniques.Charge conjugate modes are implied everywhere, unless otherwise stated.Only the sub-detectors relevant for this data analysis are described below.A more complete and detailed description of the ALICE detector and its performance can be found in Refs.[35,36].
The detectors used in this analysis include the Inner Tracking System (ITS), the Time Projection Chamber (TPC) and the Time-Of-Flight detector (TOF).These detectors are located in a large solenoid magnet producing a magnetic field of 0.5 T parallel to the LHC beam axis.The ITS consists of six cylindrical layers of silicon detectors, placed at radial distances ranging from 3.9 cm to 43 cm from the nominal beam axis and covering the full azimuth.The two innermost layers consist of Silicon Pixel Detectors (SPD), the two intermediate layers of Silicon Drift Detectors (SDD) and the two outermost layers of Silicon Strip Detectors (SSD).The total material budget of the ITS is on average 7.7% of a radiation length, for particles with η = 0 [37].The ITS spatial resolution enables the measurement of the distance of closest approach (d 0 ) of tracks to the primary vertex with a resolution better than 75 µm in the Ξ 0 c production in pp collisions at √ s = 7 TeV ALICE Collaboration transverse plane for p T > 1 GeV/c in pp collisions [38].The TPC is a cylindrical gaseous detector with a volume of about 90 m 3 .The TPC provides track reconstruction with up to 159 space points at radial distances from the beam axis ranging between 85 cm and 247 cm, within the full azimuth.The TPC cluster-position resolution is about 500 µm along the beam direction and in the transverse direction for tracks with η = 0 [39].The TPC also provides particle identification capabilities via the measurement of the specific ionisation energy loss, dE/dx, with a resolution of approximately 5.2% in pp collisions [36].
The TOF detector consists of multi-gap resistive plate chambers placed at a radial distance of 3.7 m from the beam axis and also covers the full azimuth.The TOF detector, with a timing resolution of about 80 ps, measures the time-of-flight of particles relative to the time of the collision, which is determined by the arrival time of the particles at the TOF detector and by the T0 detector, an array of Cherenkov counters placed at +370 cm and −70 cm from the nominal interaction point along the beam axis [40].
The analysed data sample consists of pp collisions at √ s = 7 TeV recorded during the 2010 LHC data taking period with a minimum bias trigger that requires at least one hit in either the SPD or the V0 detectors.The two layers of the SPD detector cover |η| < 2.0.The two V0 detectors, each comprising 32 scintillator tiles, are installed on both sides of the interaction point and cover −3.7 < η < −1.7 and 2.8 < η < 5.1.The trigger condition captures 87% of the pp inelastic cross section [41].The collision vertex is reconstructed with an efficiency of 88% and only events with a reconstructed vertex within 10 cm from the nominal interaction point along the beam direction are used in this analysis.Pile-up events are identified by searching for a second interaction vertex, reconstructed with at least three SPD tracklets (that are two-point track segments connecting hits in the two SPD layers) pointing to a common vertex, which is separated from the first vertex by at least 8 mm.After the selections, the analysed sample corresponds to an integrated luminosity L int = 5.9 ± 0.2 nb −1 .
The Ξ 0 c candidates are defined from e + Ξ − pairs by combining a track originating from the primary vertex (denoted by "electron track" in the following) and a reconstructed Ξ − baryon.Electron tracks satisfying |η| < 0.8 and p T > 0.5 GeV/c are required to have at least 100 associated clusters in the TPC (out of which at least 80 are used for the calculation of the dE/dx signal), a χ 2 normalised to the number of TPC clusters smaller than 4 and at least 4 hits in the ITS.It is also required that the electron track has associated hits in the two innermost layers of the ITS, in order to reject electrons from photon conversions occurring in the detector material outside the innermost SPD layer [13].Electrons are identified using the dE/dx measurement in the TPC and the time-of-flight measurement of the TOF detector.In both cases, the selection is applied on the n TPC σ and n TOF σ variables defined as the difference between the measured dE/dx or time-of-flight values and the one expected for electrons, divided by the corresponding detector resolution.The following selection criteria are applied: |n TOF σ | < 3 and −3.9 + 1.2p T − 0.094p 2 T < n TPC σ (p T ) < 3. The p T -dependent lower limit on n TPC σ was optimised to reject hadrons.Thus, an electron purity of 98% is achieved over the whole p T range.
The background from "photonic" electrons (originating from Dalitz decays of neutral mesons and photon conversions in the detector material) remaining in the electron sample are identified using a technique based on the invariant mass of e + e − pairs [42].The electron tracks are paired with opposite-sign tracks from the same event passing loose selection criteria (|n TPC σ | < 5 without TOF requirement) and are identified as photonic electrons if there is at least one pair with an invariant mass smaller than 50 MeV/c 2 .Setting such loose electron identification criteria is meant to increase the efficiency of finding the partners.This improves the signal-to-background ratio for Ξ 0 c by about 50%, while the fraction of the signal lost due to misidentifications is less than 2 %.The Ξ − baryons are reconstructed from the decay chain Ξ − → π − Λ, followed by Λ → pπ − .Tracks used to define Ξ − candidates are required to have at least 80 clusters in the TPC and a dE/dx signal in the TPC consistent with the expected values for protons (pions) within 4σ .The Ξ − and Λ baryons have long lifetimes (cτ of about 4.91 cm and 7.89 cm, respectively [34]), and thus they can be identified using their characteristic cascade-like or V-shaped decay topologies [43][44][45].Pions originating directly from Ξ − Ξ 0 c production in pp collisions at √ s = 7 TeV ALICE Collaboration decays are selected by requiring d 0 > 0.02 cm; protons and pions originating from Λ decays are required to have d 0 > 0.07 cm.The d 0 of the Λ trajectory to the primary vertex is required to be larger than 0.03 cm, while its cosine of the pointing angle, which is the angle between the reconstructed Λ momentum and the line connecting the Λ and Ξ − decay vertices, is required to be larger than 0.98.The distances of the Ξ − and Λ decay vertices from the beam line are required to be larger than 0.4 and 2.7 cm, respectively.These selection criteria are tuned to reduce the background, while keeping a high efficiency for the signal.Figure 1 shows the Ξ − peak in the π − Λ invariant-mass distribution integrated over p T .Only Ξ − candidates with invariant masses within 8 MeV/c 2 from the Ξ − mass (1321.71± 0.07 MeV/c 2 [34]) indicated by an arrow in Figure 1 are kept for further analysis.In this interval, the signal-to-background ratio is about 8.
The e + Ξ − pairs are formed from selected positrons and Ξ − candidates.Only pairs with an opening angle smaller than 90 degrees are used for the analysis.The background in the e + Ξ − pair distribution is estimated by exploiting the fact that Ξ 0 c baryons decay into e + Ξ − ν e (right-sign, RS), but not into e − Ξ − ν e (wrong-sign, WS), while most of the background sources contribute equally to RS and WS pairs.The yield of WS pairs is therefore used to estimate the background and is subtracted from the yield of RS pairs to obtain the Ξ 0 c raw yield.The procedure is verified with PYTHIA 6.4.21 [46] simulations using the Perugia-0 tune [47] and the GEANT3 transport code [48], including a realistic description of the detector response and alignment during the data taking period.A similar procedure was adopted by the ARGUS and CLEO collaborations studying e + e − collisions [24, 25].In order to obtain the p T -differential production cross section of Ξ 0 c baryons, the background-subtracted (WS-subtracted) yield needs to be corrected for: the signal loss due to misidentification of photonic electrons, the Ξ b contribution in the WS pairs, the missing neutrino momentum, the detector acceptance and the track-reconstruction and the candidate-selection efficiencies.No correction is applied for possible c mass [34].
differences in the acceptance of RS and WS pairs, which are found to be negligible for the current analysis based on a study with the mixed-event technique (i.e. by pairing electrons and Ξ − from different events).
The first correction accounts for the signal loss caused by the misidentification of photonic electrons.The misidentification occurs when electrons from Ξ 0 c decays accidentally have opposite-sign partners giving rise to a very small invariant mass of the e + e − pair.The misidentification probability is estimated to be less than 2% by applying the tagging algorithm to e + e + and e − e − pairs.The correction is applied as a function of the p T of the e + Ξ − pair.
The second correction accounts for the overestimation of the background caused by Ξ b → e − Ξ − ν e X decays, which produce WS pairs.Since the branching ratio of Ξ b into e − Ξ − ν e X and the Ξ b cross section in pp collisions at LHC energies have not been measured yet, two assumptions are made to estimate this contribution.First, the shape of the transverse momentum distribution of the Ξ b baryon is assumed to be the same as that of Λ 0 b , which was measured for p T > 10 GeV/c and |y| < 2 by the CMS collaboration [49].This measurement is extrapolated to p T = 0 using the Tsallis function, whose parameters were also determined by the CMS collaboration by fitting the measured distribution.
The fit parameters are consistent with those determined by the LHCb collaboration for the measurement of Λ 0 b down to p T = 0 at forward rapidity (2 < y < 4.5) [50].The second hypothesis is made for the total yield of Ξ b → e − Ξ − ν e X , which is determined by using the measurements of BR(b [52] in e + e − collisions and by assuming that the fraction of beauty quarks that hadronise into Λ 0 b and Ξ b baryons are the same as those in e + e − collisions.This assumption is supported by B-meson measurements, which show that the yield of B 0 s mesons relative to non-strange B mesons is consistent in e + e − and pp collisions [53].The Ξ b distribution obtained with these assumptions is further processed to take into account the detector acceptance, efficiency and the momentum carried by non-reconstructed decay particles.This is done with the PYTHIA 6 simulation using GEANT3 for particle transport through the detector.The correction increases with p T and reaches 2% at the highest p T interval.
The transverse momentum distribution of e + Ξ − pairs is corrected for the missing momentum of the neutrino using unfolding techniques.The response matrix to correct for the missing neutrino momentum is generated based on the correlation between the p T of the Ξ 0 c baryon and that of the reconstructed e + Ξ − pair, which is obtained from the simulation described above and is shown in Figure 3.The response matrix includes both the decay kinematics and the instrumental effects, such as energy loss and bremsstrahlung in the detector material.The response matrix needs to be determined using a realistic Ξ 0 cbaryon p T distribution.However, the distribution is not known a priori.Therefore, the response matrix is prepared in two steps.In the first step, the response matrix is obtained with the p T distribution generated with PYTHIA 6.The resulting Ξ 0 c momentum distribution is used to produce the response matrix for the second iteration.The unfolding is performed with the RooUnfold [54] implementation of the Bayesian unfolding technique [55], which is an iterative method based on Bayes' theorem.Convergence of the Bayesian method is achieved after three iterations.The p T -differential production cross section of Ξ 0 c baryons multiplied by the branching ratio into the considered semileptonic decay channel is calculated from the yields obtained by the unfolding approach as follows: where N Ξ 0 c is the yield in a given p T interval with width ∆p T .The yield is divided by the integrated luminosity L int of the analysed sample and by the product of the branching ratios of the decays Ξ − → π − Λ (99.887 ± 0.035% [34]) and Λ → pπ − (63.9 ± 0.5% [34]), which is indicated as BR Ξ − .The factor 1/2 is needed because the cross section is computed for the average of Ξ 0 c and Ξ 0 c , while the raw yield includes both contributions.The factor (A × ε) is the product of the geometrical acceptance (A) and the reconstruction and selection efficiency (ε) for Ξ 0 c → e + Ξ − ν e decays determined for Ξ 0 c generated in |y| < 0.8.Finally, the yield is normalized to one unit of rapidity by dividing it by ∆y = 1.6 under the assumption that the rapidity distribution of Ξ 0 c is uniform in the range |y| < 0.8.This assumption is verified with an accuracy of 1% using PYTHIA 6.Note that the flatness of the rapidity distribution in |y| < 0.8 is also relevant for the comparison to the D 0 meson cross section, which was determined in |y| < 0.5 [21].
The acceptance and the efficiency are calculated from the simulations with an additional correction to take into account the fact that the elastic cross section of anti-protons is not accurate in GEANT3 [56].The correction is calculated using the GEANT4 transport code [57], which has a more accurate description of the cross section, and found to be less than 2%.Since the acceptance and the efficiency depend on the Ξ 0 c -baryon p T , the Ξ 0 c should be generated with a realistic momentum distribution.This was obtained via a two-step procedure similar to that used for the response matrix.The systematic uncertainty on the Ξ 0 c cross section has different contributions, which are the uncertainties on the raw yield (owing to the procedure of background estimation), on the (A × ε) factor (due to imperfections in the simulated samples), on the correction of the missing neutrino momentum (related to the unfolding procedure) and on the normalisation.Table 1 summarises the estimated systematic uncertainties, reporting their values in all the p T intervals.The total systematic uncertainty is determined by adding the individual contributions in quadrature in each p T interval.
The systematic uncertainty on the raw yield includes the uncertainties due to the WS subtraction procedure and to the estimation of the Ξ b contribution.In the WS subtraction procedure described above, it was assumed that all the background sources contribute equally to RS and WS pairs.This is true as long as the background comprises uncorrelated pairs of electrons and Ξ − .A systematic uncertainty of 4% on the Ξ 0 c signal yield due to possible differences between RS and WS is estimated from simulations with the PYTHIA 6 event generator by checking the remaining contamination of background pairs in the RS yield after the subtraction of the WS pairs.The WS subtraction could also be affected by the amount of hadron contamination in the electron sample and the signal-to-background ratio of the Ξ 0 c signal.This effect is studied by repeating the analysis with different electron identification criteria.The results obtained with these modified criteria are found to be consistent with the ones from the default selections and therefore no systematic uncertainty is assigned.The systematic uncertainty due to the Ξ b contribution to the WS pairs is estimated by varying the Ξ b momentum distribution within the quoted uncertainty of about 50% on the cross section of Λ 0 b in pp collisions [49] and the quoted uncertainty of about 50% on the ratio of the fragmentation fractions of beauty quarks into Λ 0 b and Ξ b in e + e − collisions [51,52].The effect on the final results is found to be about 1% because the contribution from Ξ b is small.These systematic uncertainties add up to a total uncertainty of 5% for the raw yield extraction.The systematic uncertainties arising from the reconstruction and selection efficiencies are estimated by repeating the analysis with different selection criteria for electrons, Ξ − and e + Ξ − pairs and by comparing the corrected yields.Due to the statistical limitations of the Ξ 0 c sample, the electron efficiencies are studied via variations of the track-quality criteria and of the n σ values for the electron identification with TPC and TOF in the Λ + c → e + Λν e decays, which are analysed with the same procedure and have higher statistical significance.The RMS of the deviations of the corrected yields relative to the value obtained with the standard selection criteria, which amounts to 4% and 3%, is then assigned as a systematic uncertainty on the reconstruction and selection efficiency.Similarly, a systematic uncertainty of 1% on both the Ξ − reconstruction and selection efficiency is estimated from the RMS deviation of the inclusive Ξ − corrected yield against variations of the criteria applied to select the Ξ − decay tracks and its cascade 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 [44] is considered and summed in quadrature with that estimated from the variation of the selection criteria.The uncertainties on the electron and Ξ − track-quality criteria are considered as correlated and combined linearly.The uncertainty on the e + Ξ − pair selection efficiency is estimated by varying the selection criteria on the opening angle and the invariant mass of the pair and a systematic uncertainty of 3-27% is assigned depending on p T .Finally, a systematic uncertainty may also arise from an imperfect description of the acceptance of e + Ξ − pairs in the simulation.It is estimated to be 11% by comparing the azimuthal distributions of inclusive electrons and Ξ − baryons in the data and in the simulation.The uncertainty on the e + Ξ − pair acceptance is summed in quadrature with that on the electron and Ξ − selection efficiencies, resulting in a systematic uncertainty on the (A × ε) correction factor ranging from 13% to 30% depending on p T .
The systematic uncertainty on the missing neutrino momentum correction with the unfolding procedure is evaluated by varying the prior distribution to the Bayesian unfolding and by using different unfolding techniques, such as the χ 2 minimisation method [58, 59] and the Singular Value Decomposition (SVD) method [60].The RMS deviation of the results, ranging between 4% and 29% depending on p T , is assigned as a systematic uncertainty.A systematic uncertainty of 3% is also assigned due to the imperfect knowledge of the Ξ 0 c -baryon p T distributions used as input for the efficiency calculation and the unfolding procedure from the simulation.It is estimated from the difference induced in the result by adding an additional step in the iterative procedure described above to obtain the input p T distributions.These systematic uncertainties add up to an uncertainty ranging between 6% and 29% depending on p T .
Finally, the results have a 3.5% normalisation systematic uncertainty arising from the uncertainty in the determination of the minimum-bias trigger cross section in pp collisions at √ s = 7 TeV [41].
The p T -differential cross section of Ξ 0 c baryons multiplied by the branching ratio into e + Ξ − ν e is shown in Figure 5 for the p T interval 1 < p T < 8 GeV/c at mid-rapidity, |y| < 0.5.The error bars and boxes represent the statistical and systematic uncertainties, respectively.The feed down contribution from Ξ b , e. g.
, is not subtracted due to the lack of knowledge of the absolute branching ratios of Ξ b → Ξ 0 c + X .The ratio of the p T -differential cross section of Ξ 0 c baryons to that of D 0 mesons [21] is shown in Figure 6.The p T intervals of the cross-section measurements are combined to have the same p T bin boundaries for Ξ 0 c and D 0 .The systematic uncertainty in a merged p T interval is defined by propagating the yield extraction uncertainties of the D 0 measurement as uncorrelated among p T intervals and all the other uncertainties of the D 0 and Ξ 0 c measurements as correlated.The systematic uncertainty on the Ξ 0 c /D 0 ratio is calculated treating all the uncertainties on the Ξ 0 c and D 0 cross sections as uncorrelated, except for the normalisation uncertainty that cancels out in the ratio.The ratio integrated in the transverse momentum interval 1 < p T < 8 GeV/c is found to be (7.0 ± 1.5(stat) ± 2.6(syst)) × 10 −3 .
In Figure 6(a), the measured transverse momentum dependence of the Ξ 0 c /D 0 ratio is compared with predictions from the PYTHIA 8.211 event generator [46,65].PYTHIA 8 uses 2 → 2 processes followed by a leading-logarithmic p T -ordered parton shower for the charm quark pair production and the hadronisation is treated with the Lund string model [66].The figure shows the results obtained with different tunes of hadronisation: the Monash 2013 tune [62] and the Mode 0 tune from [28].The latter is based on a model for the hadronisation of multi-parton systems, which includes string formation beyond the leading-colour approximation and is implemented in PYTHIA 8 with specific tuning of the colour reconnection parameters.As compared to the Monash 2013 tune, this model provides a better description of the measured baryon-to-meson ratios in the light-flavour sector.Two other tunes (Mode 2 and Mode 3) provided in Ref. [28] give similar Ξ 0 c /D 0 ratios as Mode 0. In Figure 6(b), the measured ratio is also compared to other models implementing different hadronisation mechanisms: DIPSY [63] with the rope hadronisation [67] and HERWIG 7.0.4[64] with the cluster hadronisation [68].To compare the data with these models, theoretical calculations of the branching ratio, which range between 0.83% and 4.2% [69][70][71], are used.This range defines the width of the bands shown for the model calculations represented in Figure 6.Although the predictions of the Mode 0 tune of PYTHIA 8 are the closest to the data compared to the other models, all calculations underestimate the measured ratio significantly.Thus, this new measurement can provide an important constraint to the models of charm quark hadronisation in pp collisions, once a measurement of the absolute branching ratio of the Ξ 0 c will become available.In summary, we reported on the first LHC measurement of the inclusive p T -differential production cross section of the charm-strange baryon Ξ 0 c multiplied by the branching ratio into e + Ξ − ν e in pp collisions at  [63] and HERWIG 7 [64], are shown as shaded bands representing the range of the currently available theoretical predictions for the branching ratio of the considered Ξ 0 c decay mode.
√ s = 7 TeV.The ratio of this measurement integrated over 1 < p T < 8 GeV/c to the production cross section of the D 0 meson integrated over the same p T interval was found to be (7.0 ± 1.5(stat) ± 2.6(syst)) × 10 −3 .Several event generators with various models and tunes for the hadronisation mechanism underestimate the measured ratio.

Figure 2 (
Figure 2(a) shows the invariant-mass distributions of RS and WS pairs, integrated over the whole p T interval.The invariant-mass distribution of Ξ 0 c candidates obtained by subtracting the WS pair yield from the RS one is shown in Figure 2(b) together with the signal distribution from the simulation, which is normalised to the measured RS−WS yield.The shapes of the two distributions are found to be consistent with each other.Due to the missing momentum of the neutrino, the invariant-mass distribution of the e + Ξ − pair does not peak at the Ξ 0 c mass (2470.85+0.28 −0.40 MeV/c 2 [34]) indicated by an arrow in Figure 2(b).The invariant mass of e + Ξ − pairs from Ξ 0 c decays is bounded by the Ξ 0 c mass due to the missing momentum of the neutrino.Thus only e + Ξ − pairs satisfying m eΞ < 2.5 GeV/c 2 are selected for further analysis.

Fig. 1 :
Fig. 1: Invariant-mass distribution of Ξ − → π − Λ (and charge conjugate) candidates integrated over p T .The arrow indicates the world average Ξ − mass from Ref. [34] and the dashed lines indicate the selected interval for the Ξ − candidates.

Fig. 2 :
Fig. 2: (a) Invariant-mass distributions of right-sign and wrong-sign (and charge conjugate) pairs integrated over the whole p T interval.(b) Invariant-mass distribution of Ξ 0 c candidates obtained by subtracting the wrong-sign pair yield from the right-sign one compared with the signal distribution from the simulation, which is normalised to the measured RS−WS yield.The arrow indicates the Ξ 0 c mass [34].

Fig. 3 :
Fig. 3: Correlation between the generated Ξ 0 c -baryon p T and the reconstructed e + Ξ − pair p T , obtained from the simulation based on PYTHIA 6 described in the text.

Figure 4 Fig. 4 :
Fig. 4: Product of acceptance and efficiency (A × ε) of Ξ 0 c baryons generated in |y| < 0.8 decaying into e + Ξ − ν e as a function of p T , determined from simulations PYTHIA 6 (see text).

Fig. 5 :
Fig.5: Inclusive Ξ 0 c -baryon p T -differential production cross section multiplied by the branching ratio into e + Ξ − ν e , as a function of p T for |y| < 0.5, in pp collisions at √ s = 7 TeV.The error bars and boxes represent the statistical and systematic uncertainties, respectively.The contribution from Ξ b decays is not subtracted.

Fig. 6 :
Fig.6: Ratio of the p T -differential cross sections of Ξ 0 c baryons (multiplied by the branching ratio into e + Ξ − ν e ) and D 0 mesons [21] as a function of p T for |y| < 0.5, in pp collisions at √ s = 7 TeV.The error bars and boxes represent the statistical and systematic uncertainties, respectively.Predictions from theoretical models, (a) PYTHIA 8 with different tunes[28,62] (b) DIPSY[63] and HERWIG 7[64], are shown as shaded bands representing the range of the currently available theoretical predictions for the branching ratio of the considered Ξ 0 c decay mode.

Table 1 :
Summary of systematic uncertainties on the p T -differential cross section of Ξ 0 c → e + Ξ − ν e for 5 p T intervals.The uncertainty on the missing neutrino momentum is denoted as p ν T in the table.