Measurement of the Λ 0 b lifetime and mass in the ATLAS experiment

A measurement of the Λ 0 b lifetime and mass in the decay channel Λ 0 b → J/ψ ( µ + µ − )Λ 0 ( pπ − ) is presented. The analysis uses a signal sample of about 2200 Λ 0 b and ¯Λ 0 b decays that are reconstructed in 4.9 fb − 1 of ATLAS pp collision data collected in 2011 at the LHC center-of-mass energy of 7 TeV. A simultaneous mass and decay time maximum likelihood ﬁt is used to extract the Λ 0 b lifetime and mass. They are measured to be τ Λ b = 1 . 449 ± 0 . 036(stat) ± 0 . 017(syst) ps and m Λ b = 5619 . 7 ± 0 . 7(stat) ± 1 . 1(syst) MeV. A measurement of the Λ 0 b lifetime and mass in the decay channel Λ 0 b → J/ψ ( µ + µ − )Λ 0 ( pπ − ) is presented. The analysis uses a signal sample of about 2200 Λ 0 b and ¯Λ 0 b decays that are reconstructed in 4.9 fb − 1 of ATLAS pp collision data collected in 2011 at the LHC center-of-mass energy of 7 TeV. A simultaneous mass and decay time maximum likelihood ﬁt is used to extract the Λ 0 b lifetime and mass. They are measured to be τ Λ b = 1 . 449 ± 0 . 036(stat) ± 0 . 017(syst) ps and m Λ b = 5619 . 7 ± 0 . 7(stat) ± 1 . 1(syst) MeV.

Measurement of the Λ 0 b lifetime and mass in the ATLAS experiment The ATLAS Collaboration (Dated: April 23, 2013) A measurement of the Λ 0 b lifetime and mass in the decay channel Λ 0 b → J/ψ(µ + µ − )Λ 0 (pπ − ) is presented. The analysis uses a signal sample of about 2200 Λ 0 b andΛ 0 b decays that are reconstructed in 4.9 fb −1 of ATLAS pp collision data collected in 2011 at the LHC center-of-mass energy of 7 TeV. A simultaneous mass and decay time maximum likelihood fit is used to extract the Λ 0 b lifetime and mass. They are measured to be τΛ b = 1.449 ± 0.036(stat) ± 0.017(syst) ps and mΛ b = 5619.7 ± 0.7(stat) ± 1.1(syst) MeV.

I. INTRODUCTION
The Λ 0 b baryon (and its charge conjugateΛ 0 b ) is the lightest baryon containing a b (b) quark. With a mass of about 5620 MeV [1,2] it is not produced at B factories, where the collision center-of-mass energy is tuned to produce pairs of B mesons. Currently, hadron colliders are the only facilities where the properties of b baryons can be studied. This paper presents a measurement of the Λ 0 b mass and lifetime in the ATLAS experiment [3] using the decay channel Λ 0 b → J/ψ(µ + µ − )Λ 0 (pπ − ) (the charge conjugate mode is implied throughout the paper unless explicitly stated otherwise). The Λ 0 b lifetime, although measured by many experiments [1,[4][5][6], still suffers from a large experimental uncertainty.
The decay B 0 d → J/ψ(µ + µ − )K 0 S (π + π − ) has the same topology as the studied Λ 0 b decay. The B 0 d mass and lifetime are measured with good precision [1], and therefore this decay provides a useful tool to validate the Λ 0 b results, as both measurements are subject to similar systematic uncertainties. The lifetime ratio, τ Λ b /τ B d , can be predicted by heavy quark expansion calculations [7] and perturbative QCD [8] and is of great theoretical interest. The lifetime and mass are determined using a simultaneous unbinned maximum likelihood fit to the reconstructed mass and decay time of each selected candidate.

II. DATA SAMPLES AND TRIGGER SELECTION
The ATLAS experiment [3] is a general-purpose detector at the Large Hadron Collider (LHC). It covers nearly the entire solid angle around the interaction point with layers of tracking detectors, calorimeters, and muon chambers. The coordinate system has the z axis aligned with the beam direction. The transverse momentum p T and pseudorapidity, η, of reconstructed particles are defined with respect to that direction. This analysis uses two ATLAS subsystems: the inner detector (ID) and the muon spectrometer (MS). Both are situated in a magnetic field and serve as tracking detectors. The ID consists of three types of detector: the silicon pixel detector (Pixel), the silicon microstrip detector (SCT), and the transition radiation tracker (TRT). The MS consists of monitored drift tube chambers (MDT) and cathode strip chambers (CSC) for precision muon measurements, resistive plate chambers (RPC) and thin gap chambers (TGC) employed by the muon-trigger system. Tracks are reconstructed in the ID, and the MS is used to identify muons. Only tracks with p T above 400 MeV and pseudorapidity |η| < 2.5 are used in this analysis.
The analysis is based on data collected in 2011 using single-muon, dimuon, and J/ψ triggers. The AT-LAS trigger system [9] has three levels: the hardwarebased Level-1 trigger and the two-stage High Level Trigger (HLT). At Level-1 the muon trigger uses dedicated fast muon-trigger chambers to search for patterns of hits corresponding to muons passing different p T thresholds. Regions of interest (RoI) around these Level-1 hit patterns then serve as seeds for the HLT muon reconstruction. Since the rate from the low-p T muon triggers was too high for all accepted events to be saved, prescale factors were applied to a subset of the triggers to reduce the output rate. The muon transverse momentum thresholds for single and dimuon triggers range from 4 to 22 GeV. The J/ψ dimuon triggers require that the muons originate from a common vertex and have opposite charge, and that the dimuon mass is in the range 2.5 GeV < m µµ < 4.3 GeV. The majority of the sample was collected by the J/ψ trigger with a p T threshold of 4 GeV applied to each muon. This was the lowest-p T unprescaled trigger in the 2011 data taking; however, other complementary triggers were used, too. The p T spectrum of the selected muons peaks at 5 GeV; the lowest muon p T is above 2.5 GeV.
A Monte Carlo (MC) sample of 5 × 10 6 antibaryonΛ 0 b events is used to study systematic effects and to correct for the efficiency and acceptance of the detector. The sample is generated using the Pythia 6 MC generator [10] with the 2011 ATLAS AUET2B L0 * * tune [11], and the events are filtered so that each accepted event has a decayΛ 0 b → J/ψ(µ + µ − )Λ 0 (pπ + ) with the muons having transverse momenta of at least 2.5 GeV. The MC sample is generated with a Λ 0 b lifetime of τ MC Λ b = 1.391 ps.

III. RECONSTRUCTION AND SIGNAL SELECTION
A. J/ψ and Λ 0 preselection The decay Λ 0 b → J/ψ(µ + µ − )Λ 0 (pπ − ) has a cascade topology. The J/ψ decays instantly at the same point as the Λ 0 b (secondary vertex), while the Λ 0 lives long enough to form a displaced tertiary vertex. There are four finalstate particles: two muons from the J/ψ, and a proton and a pion from the Λ 0 decay.
The dimuon and dihadron (V 0 ) pairs are preselected by requiring that their tracks be successfully fitted to a common vertex satisfying some basic quality requirements. The J/ψ and V 0 pre-selection is very loose, so that potential candidates are not excluded at this stage. The dimuon candidates are accepted if the J/ψ vertexrefitted invariant mass lies in the range 2.8 GeV < m µµ < 3.4 GeV. The dihadron candidates are accepted if the invariant mass is in the range 1.08 GeV < m pπ < 1.15 GeV. The masses of a proton and a pion are assigned to the tracks when the invariant mass is calculated; pπ − and pπ + combinations are tested so that both Λ 0 andΛ 0 candidates are accepted.
The muon and hadron track pairs preselected with the criteria described in the previous section are then refitted with a constraint of a Λ 0 b → J/ψ(µ + µ − )Λ 0 (pπ − ) topology. The muons are constrained to intersect at a single vertex, while their invariant mass is set equal to the known mass of the J/ψ, m J/ψ = 3096.92 MeV [1]. The two hadronic tracks are constrained to a second vertex, and their invariant mass is fixed to the mass of the Λ 0 , m Λ 0 = 1115.68 MeV [1]. The combined momentum of the refitted V 0 track pair is constrained to point to the dimuon vertex in three dimensions. The fit is performed on all four tracks simultaneously, taking into account the constraints described above (cascade topology fit) and the full track error matrices. The quality of the fit is characterized by the value of χ 2 /N dof , where a global χ 2 involving all four tracks is used. The corresponding number of degrees of freedom, N dof , is six. Furthermore, for each track quadruplet, that can be successfully fitted to ) topology fit is attempted (i.e. a pion mass is assigned to the hadronic tracks and the V 0 mass is constrained to the mass of K 0 S , m K S = 497.65 MeV [1]). This is to identify possible B 0 d decays mistaken for Λ 0 b . The Λ 0 b candidates are then subjected to the following selections: • The global χ 2 /N dof < 3.
• The transverse momentum of the cascade-refitted V 0 , p T,V 0 > 3.5 GeV. • The transverse decay length of the cascade-refitted V 0 vertex measured from the Λ 0 b vertex, L xy,V 0 > 10 mm.
• The invariant mass must be in the range 5.38 GeV < m J/ψΛ 0 < 5.90 GeV.
• If the four tracks forming a Λ 0 b candidate also result in an acceptable B 0 d fit, the candidate must have a difference of cumulative χ 2 probabilities of the two fits, P Λ 0 b − P B d > 0.05. With these criteria, 4074 Λ 0 b and 4081Λ 0 b candidates (including background) are selected. No track quadruplet is successfully fitted as both a Λ 0 b and aΛ 0 b decay. The mass distributions of the selected candidates are shown in Fig. 1. In the rest of the paper the Λ 0 b andΛ 0 b samples are combined.

IV. MASS AND PROPER DECAY TIME FIT
The proper decay time of the Λ 0 b candidate is calculated from the measured decay distance and the candidate's momentum as follows: where m PDG = 5619.4 MeV [1], p T is the reconstructed Λ 0 b transverse momentum, and L xy is the Λ 0 b transverse decay distance measured from the primary vertex (PV). On average there are 6.8 collision vertices per event in the selected data resulting from multiple collisions at each LHC bunch crossing (pileup events). The collision vertex that in three-dimensional space lies closest to the trajectory of the reconstructed Λ 0 b candidate is used as the PV.
An unbinned maximum likelihood fit is used to determine the Λ 0 b mass and lifetime. The mass and proper decay time are fitted using a likelihood function defined as follows: where f sig denotes the fraction of signal candidates; m i is the invariant mass of the i-th candidate and τ i is its proper decay time. The corresponding errors, δ mi and δ τi , are estimated on a candidate-by-candidate basis by the cascade topology fit. M s and M b are probability density functions (PDFs) describing the signal and background mass dependence; T s and T b describe the dependence on the proper decay time. The invariant mass and proper decay time error distributions, , are extracted from data. It has been verified that using separate PDFs for the signal and background component produces the same result when a single PDF is used, w ≡ w s = w b . For this reason the latter case is used.
The background can be divided into two categories: prompt and nonprompt backgrounds. The prompt background consists of J/ψ candidates produced directly in the pp collision that are randomly combined with V 0 candidates, which also include fake combinatorial Λ 0 or K 0 S candidates. The prompt background decay length is due to the finite resolution of the vertex reconstruction. The nonprompt background includes events where the J/ψ candidate originates in the decay of a b hadron. This type of background has a lifetime due to its origin in long-lived b hadrons [e.g.
S meson misidentified as Λ 0 , forming a nonprompt background for Λ 0 b ]. The signal component of the mass PDF, M s , is a Gaussian function with a mean equal to m Λ b and width S m δ m . The mass error scale factor S m determines how much the errors δ mi are overestimated or underestimated. The background component is a first order polynomial with a slope b.
Using the estimated decay time error δ τ , the proper decay time resolution is modeled with a Gaussian function: where S τ denotes the proper decay time error scale factor, and τ and τ stand for the reconstructed and true proper decay times, respectively. The signal and nonprompt background proper decay time distributions are modeled as exponential functions, E(τ ; τ B ), for τ > 0, with τ B being the fitted parameter denoting either the Λ 0 b lifetime or the pseudolifetime of the long-lived background. The prompt background component is modeled by a sum of two functions: a Dirac δ function δ Dirac (τ ) and a symmetric exponential (Laplace distribution) E sym (τ ), to account for the non-Gaussian tails of the prompt background observed in data.
The functions are convolved with the resolution model (1) to obtain the PDFs of the measured proper decay time: with the nonprompt and prompt components defined as The efficiency correction function ε(τ ) in Eq. (2) accounts for the decay-time-dependent selection bias. Two sources are responsible for the selection bias in the Λ 0 b decay time: the V 0 reconstruction efficiency and the trigger selection. The V 0 reconstruction efficiency depends on the decay distance from the center of the detector, as tracks from decays further away from the center leave fewer hits in the ID. Since the Λ 0 b decay length and the distance of the Λ 0 vertex from the center of the detector are correlated (the latter includes the former), this biases the measured proper decay time toward smaller values. The other source of the bias is the muon trigger, which affects the distribution of the muon transverse impact parameter, d 0 . Applying the tag-and-probe method to J/ψ decays, the trigger efficiency as a function of d 0 is measured for a single-muon trigger in data. The simulation shows that the dimuon-trigger efficiency can be expressed as a product of single-muon efficiencies.
The MC events are reweighted to reproduce the observed trigger bias. The efficiency correction, ε(τ ), is denotes the slope of the efficiency correction. The exponential form is chosen for ε(τ ) because it describes the MC well and is particularly easy to convolve with the resolution model. The slope of the exponential, c Λ b , is extracted from a fit to the MC decay time efficiency plot shown in Fig. 2. The extracted value is c Λ b = 113±56 ps; i.e. for a decay time of 6 ps the efficiency decreases by 5%.

A. Parameters determined from the fit
The full PDF has 12 free parameters: the Λ 0 b mass and lifetime, m Λ b and τ Λ b ; the fraction of signal events, f sig ; the error scale factors, S m and S τ ; the slope of the mass dependence of the background, b; the pseudolifetimes of the long-lived background, τ bkg,1 and τ bkg,2 ; the exponential slope of the non-Gaussian prompt background, τ bkg,3 ; and the relative fractions of the various background contributions, f 1 , f 2 , and f 3 .
Other quantities are calculated from the fit parameters. The number of signal and background candidates, N sig and N bkg , are calculated as N sig = f sig N and N bkg = (1 − f sig )N , where N is the total number of candidates. The mass and proper decay time resolutions are calculated from the fit parameters, too. By analogy with a Gaussian distribution, the mass resolution σ m is defined as half of that mass range for which the integral of M s retains 68.3% of the number of signal events symmetrically around the fitted Λ 0 b mass. The proper decay time resolution, σ τ , is determined in the same fashion by integrating the prompt background PDF.

A. Results of the maximum likelihood fit
The results of the maximum likelihood fit are listed in Table I. The table shows  parameters, calculated parameters, and a χ 2 /N dof value which quantifies the fit quality. The χ 2 /N dof value is calculated from the dataset binned in mass and decay time with 61 degrees of freedom. The sizes of the bins are commensurate with the measured mass and decay time resolutions, and only bins with more than 11 entries are used for the χ 2 calculation. This requirement is imposed so that the error on the number of entries in each bin can be taken as Gaussian. The lifetime result is corrected for the selection bias (see Sec. IV); the size of the correction is +19 fs. The estimated correlation between the mass and lifetime is small, 0.002. Projections of the PDF onto the mass and proper decay time axes are shown in Fig. 3.

B. Systematic uncertainties
Systematic errors are estimated by changing various parameters of the analysis and observing the shift in the extracted mass and lifetime. The shift with respect to the baseline result is then quoted as a systematic uncertainty. The non-negligible systematic uncertainties are summarized in Table II. The individual errors are added in quadrature, yielding total systematic errors of the lifetime and mass measurements, σ syst τ = 17 fs and σ syst m = 1.1 MeV, respectively. Details of the determination of the systematic uncertainties follow.

Event selection and reconstruction bias
Two effects that lead to a selection bias for Λ 0 b candidates as a function of decay time have been identified: the dominant contribution comes from the muon trigger, which slightly biases the transverse impact parameter of muons, d 0 , toward smaller values. The second bias comes from the V 0 reconstruction. The event selection bias is corrected using the MC simulation to determine the efficiency as a function of the decay time as described in Sec. IV. The slope of the efficiency correction function, c Λ b , is determined with a statistical accuracy of 50% (56 ps). Using standard error propagation, the contribution of this uncertainty to the overall error is evaluated to be 9 fs.
Since the two bias corrections rely on MC simulation, systematic uncertainties due to the corrections are estimated. These are determined separately for the V 0 reconstruction bias and for the trigger bias. In the MC simulation the bias correction shifts the Λ 0 b lifetime by 26 fs, of which 10 fs are due to the V 0 reconstruction and 16 fs are due to the trigger requirement. The systematic error due to the V 0 bias correction is estimated by varying the Λ 0 b transverse momentum in the MC simulation (using kinematic reweighting) to probe the p T dependence of the correction. The magnitude of the variation is about 3 times the difference between the mean p T in data and the MC simulation. The corresponding error is estimated to be 4 fs. To ensure a good description of the trigger bias, the MC sample is reweighted using single-muon-trigger efficiencies expressed as a function of muon d 0 values that are extracted from data. The weighting functions are parameterized as linear functions, w(d 0 ) ∝ 1 + ad 0 , and their slope a is determined by a linear fit in bins of the measured muon p T and η. To assess the systematic error on the trigger bias correction, the weighting parameters a are varied by their errors. This produces a lifetime shift of 7 fs, which is used as a systematic error. The total systematic error calculated as a quadratic sum of the individual contributions is 12 fs.
To assess the systematic error in the mass measurement due to the event selection, the MC distribution of ∆m = m MC − m, where m MC is the generated mass, is fitted with a double Gaussian. The systematic error, given by the shift of the mean of the double Gaussian, is estimated to be 0.9 MeV. The mass shift is caused by the muon-trigger p T thresholds: muons with larger p T have higher probability of being selected than low-p T muons.
As a consequence, muons whose p T is mismeasured as larger than the true value have a higher probability of being reconstructed than muons whose p T is mismeasured as smaller, which creates a small asymmetry of the mass peak.

Background fit models
Alternative background models are used to assess the sensitivity of the results to the choice of background parametrization. A second-order polynomial and an exponential mass dependence of the M b PDF are tested. In addition the decay time dependence is modified by adding a third exponential into the nonprompt background component, T np . The alternative background PDFs fit the data well. These changes result in a lifetime shift of 2 fs and a mass shift of 0.2 MeV. In the fit model the decay time and mass are assumed to be uncorrelated. To test this assumption the fit's mass range limits, m min and m max , are varied independently by 60 MeV. This changes the relative contribution of the background from the left and right sidebands, and the mass and lifetime are extracted again for these new mass ranges. While the change of m max has a minimal impact on the extracted mass and lifetime, the change of m min produces a lifetime shift of 9 fs. This value is added to the total systematic error due to background modeling.

B 0 d contamination
The number of B 0 d candidates misidentified as Λ 0 b is estimated by a fit to the mass distribution of the candidates which fall in the Λ 0 b signal region, 5.52 GeV < m J/ψΛ 0 < 5.72 GeV, under the hypothesis that they are B 0 d → J/ψ(µ + µ − )K 0 S (π + π − ) decays. A fit to a Gaussian peak on a linear background yields 82±46 B 0 d candidates. Since these candidates are treated as Λ 0 b , their pseudolifetime is scaled up by the ratio of the Λ 0 b and B 0 d masses, = 1617 fs (the decay time change due to the difference in p T reconstructed under the two hypotheses is negligible). If all such background candidates contribute to the fitted Λ 0 b lifetime, it would cause a shift of 7 fs. This is quoted as a conservative estimate of the systematic error. The error on the mass measurement is estimated by relaxing the P Λ 0 b −P B d cut to double the estimated B 0 d background. This results in a Λ 0 b mass shift of 0.2 MeV.

Residual misalignment of the ID
The distribution of the transverse impact parameter, d 0 , of tracks originating from the PV is used to estimate the geometrical distortions due to residual misalignment. The geometry in the MC simulation is distorted by adjusting the positions of the ID modules so that the d 0 of tracks coming from the PV is biased by the same amount as observed in data. The mass-lifetime fit is performed with simulated data using the default (ideal) geometry and the sample with geometry distortions. A shift of 1 fs is observed between the two measurements and is assigned as a systematic error due to residual misalignment. No significant mass shift is observed.

Uncertainty in the amount of ID material
Inaccurate modeling of the amount of material in the ID could affect the measurement since the tracking algorithm estimates the particle energy loss using a material map. To explore this uncertainty, the MC simulation is repeated with 20% more material in the ID silicon detectors (Pixel and SCT) and their supporting services, which is large compared to the estimated uncertainty of 6%(9%) in the Pixel (SCT) detectors (see Ref. [12]). The resulting shifts of 3 fs in lifetime and 0.2 MeV in mass are conservative estimates of the systematic uncertainties from this source.

Uncertainty in the tracking momentum scale
The K 0 S mass value is used to estimate the uncertainty in the track momentum determination. The K 0 S mass extracted from a fit to the invariant mass agrees with the PDG's world average within 0.03%. Such a shift corresponds to a track momentum scale shift of 0.05%. The momentum scale can be further tested using the reconstructed J/ψ mass. The observed mass shift corresponds to a momentum scale error of −0.03%, in agreement with the assumption of ±0.05%. Shifting the momenta of all tracks in the MC simulation by this amount yields a Λ 0 b mass shift of 0.5 MeV. No significant lifetime shift is observed.

Other systematic errors
Other sources of systematic errors are investigated, such as an alternative choice of the PV (e.g. using the collision vertex whose tracks have the highest sum of p 2 T ) and the use of a PDF with the error distributions modeled separately for signal and background. These changes do not result in a significant mass or lifetime shift.
The B 0 d → J/ψ(µ + µ − )K 0 S (π + π − ) channel has the same decay topology as Λ 0 b → J/ψ(µ + µ − )Λ 0 (pπ − ) and can be used to cross-check the Λ 0 b results and to determine the ratio of the Λ 0 b and B 0 d lifetimes. The analysis and systematic error studies, described in the previous sections, are repeated for the B 0 d . The B 0 d channel is subjected to exactly the same kinematic cuts as for the Λ 0 b channel and therefore has similar systematic effects. The mass range used for the K 0 S pre-selection is 440 MeV < m π + π − < 570 MeV, and the B 0 d invariant mass must lie in the range 5.1 GeV < m J/ψK 0 S < 5.5 GeV. This selection is chosen to be as close as possible to the Λ 0 b selection rather than to maximize the B 0 d signal yield. Using the maximum likelihood fit, the B 0 d lifetime and mass are measured to be τ B d = 1.509±0.012(stat)±0.018(syst) ps and m B d = 5279.6 ± 0.2(stat) ± 1.0(syst) MeV. These values are consistent with the world averages, τ PDG
The statistical and systematic errors are propagated from the errors of the lifetime measurements. The systematic errors are conservatively assumed to be uncorrelated. This value is between the recent determination by DØ, R DØ = 0.864 ± 0.052(stat) ± 0.033(syst) [6], and the measurement by CDF, R CDF = 1.020 ± 0.030(stat) ± 0.008(syst) [5]. It agrees with the heavy quark expansion calculations which predict the value of the ratio between 0.88 and 0.97 [7] and is compatible with the next-to-leading-order QCD predictions with central values ranging between 0.86 and 0.88 (uncertainty of ±0.05) [8].