GLAS-PPE/28-9 24 th July 28 arxiv:88.867v1 [hep-ex] 6 Aug 28 Department of Physics and Astronomy Experimental Particle Physics Group Kelvin Building, University of Glasgow, Glasgow, G12 8QQ, Scotland Telephone: +44 ()141 33 2 Fax: +44 ()141 33 5881 Tracking and physics validation studies of the simplified geometry description with 2 hh decays M. Gersabeck 1, E. Rodrigues 1 1 University of Glasgow, Glasgow G12 8QQ, Scotland Abstract This note validates the usage of the simplified detector geometry description. A sample of2 hhdecayswasusedtoassessthetrackingandphysicsperformancewithrespecttowhat is obtained with the full detector description. No significant degradation of performance was found. LHCb Public Note, LHCb-28-3
Contents 1 Introduction 2 2 Pattern ognition 2 3 Track fitting 3 4 Physics analysis 1 4.1 Effect on the event selection............................. 1 4.2 Effect on resolutions................................. 1 5 Conclusions and Final Remarks 15
1 Introduction The onstruction of tracks is an important but time-consuming task. It is well known that track fitting contributes substantially to the onstruction time budget. Detailed studies show that a large fraction of the time spent in fitting tracks is due to the many calculations of material intersections along a particle s path. These are necessary in order to account for the detector material by means of multiple scattering and energy loss cortions. A simplified description of the detector material as seen by a particle traversing the LHCb detector has ently been implemented [1]. It replaces the full detector material description by a small set of simple modules (mostly boxes and cylinders) that model the average material properties. In this note we study the implications of using this simplified geometry when accounting for detector material during track fitting. We compare the performance of the simplified versus the full geometry with a sample of B π + π events in terms of track fit quality, quality of onstruction and event selection, and physics analysis. Note that all results are obtained starting from the same data sample generated and simulated with the full geometry in Geant4. 2 Pattern ognition LHCb pattern ognition algorithms ignore any material effects and should therefore be insensitive to whether the simplified geometry description is used (an exception is explained below). Those considered in this note are: finding of tracks in the vertex locator (VELO) in r-z and 3D-space. The algorithms are hereafter denoted by VeloR and VeloSpace, respectively [2]; finding of tracks that traverse the whole LHCb detector (called long tracks ). The two existing long tracking algorithms are hereafter denoted Forward [3] and Matching [4]. In Table 1 the efficiencies 1) for the VeloR, VeloSpace, Forward, and Matching pattern ognition algorithms are compared. All efficiencies are quoted for long tracks with no momentum cut applied. As expected, all efficiencies are identical, with the exception of the Matching efficiency, whose difference can be understood as the algorithm matches T-station seed tracks to VELO track segments by extrapolating them to the magnet region, and the extrapolation internally looks at the material along the track s trajectory. An identical conclusion can be drawn for the number of clone tracks and the ghost rate of all four algorithms. VeloR VeloSpace Forward Matching Geometry efficiency (%) efficiency (%) efficiency (%) efficiency (%) full 98.±.1 97.±.1 85.9±.2 81.1±.2 simplified 98. ±.1 97. ±.1 85.9 ±.2 81.4 ±.2 Table 1: VeloR, VeloSpace, Forward and Matching pattern ognition efficiencies for the full and the simplified geometries. 1) For more details about the definitions of the pattern ognition efficiencies see [5].
For completeness the tracks pseudorapidity, η, distributions as obtained with the Forward and the Matching pattern ognition algorithms are compared in Figure 1. No significant differences are observed, as expected, even in the very forward η region where effects of the simplified description are most likely to be evident as high-η tracks traverse more material..24.22.2.18.16.14 (a) Full: <>=3.5875 RMS=.76 Simplified: <>=3.5897 RMS=.7599.12.1.8.6.4.2 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 Pseudorapidity η.22.2.18.16.14.12.1.8.6.4.2 Full: <>=3.446 RMS=.8272 Simplified: <>=3.4395 RMS=.8274 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 Pseudorapidity η Figure 1: Distributions in tracks pseudorapidity as obtained with the (a) Forward and the Matching pattern ognition algorithms for the full and the simplified geometries. 3 Track fitting Pattern ognition tracks are fitted in order to obtain the best estimate of the track parameter values and errors. During the fitting procedure some of the hits on the tracks (called LHCbIDs) are flagged as outliers and removed from the tracks. The distributions of outliers removed by the fitter for Forward and Matching tracks are compared in Figure 2. Irrespective of the
.4.35.3.25 (a) Full: <>=.884 RMS=.8238 Simplified: <>=.8867 RMS=.824.2.15.1.5 1 2 3 4 5 Number of hits - number of measurements.6.5.4 Full: <>=.5487 RMS=.7467 Simplified: <>=.5738 RMS=.7591.3.2.1 1 2 3 4 5 Number of hits - number of measurements Figure 2: Distributions of outlier hits as obtained with the (a) Forward and the Matching pattern ognition algorithms for the full and the simplified geometries. geometry used, Forward tracks tend to have more outliers than Matching tracks. When using the simplified geometry, this tendency is less pronounced. In particular, Matching tracks fitted with the simplified geometry lose slightly more hits compared to when they are fitted with the full geometry. The quality of track fitting is straightforwardly assessed looking at the resolutions and the pull distributions of the track state parameters: positions x and y, slopes t x and t y, and charge-over-momentum ratio q/p. All the distributions shown in this section were obtained with the Forward algorithm. However, it has been checked that all the following conclusions also hold for long tracks from the Matching algorithm. The track parameter resolutions at the first track measurement point quantities dominated by the VELO measurements are collected in Figure 3. Neither the position nor the slope resolutions deteriorate when using the simplified rather than the full geometry. A slight increase in the momentum resolution (here taken as the root mean squared, RMS, rather than the σ of a single-gaussian fit) from.6% to.63% is observed. The increase originates mostly from a broadening in the left part of the distribution, where p < p.
(a).9.8.7 Full: <>=-. RMS=.14 Simplified: <>=-. RMS=.14.9.8.7 Full: <>=-. RMS=.12 Simplified: <>=-. RMS=.13.6.6.5.5.4.4.3.3.2.2.1.1 -.1 -.8 -.6 -.4 -.2.2.4.6.8.1 x - x (cm) -.1 -.8 -.6 -.4 -.2.2.4.6.8.1 y - y (cm) (c) (d).12.1.8 Full: <>=. RMS=.4 Simplified: <>=. RMS=.4.12.1.8 Full: <>=. RMS=.4 Simplified: <>=. RMS=.4.6.6.4.4.2.2 -.25-.2-.15-.1-.5.5.1.15.2.25 t x() - t x() (rad) -.25-.2-.15-.1-.5.5.1.15.2.25 t y() - t y() (rad) (e).1.8 Full: <>=-.3 RMS=.6 Simplified: <>=-.6 RMS=.63.6.4.2 -.5 -.4 -.3 -.2 -.1.1.2.3.4.5 (p - p )/p Figure 3: Resolutions on the track parameters at the first measurement point for the full and the simplified geometries. This sample of long tracks was obtained with the Forward pattern ognition algorithm.
(a).35.3.25.2.15 Full: <>=.17 RMS=1.2321 Simplified: <>=.4 RMS=1.2466.3 Full: <>=-.18 RMS=1.343.25.2.15 Simplified: <>=.7 RMS=1.326.1.5.1.5-5 -4-3 -2-1 1 2 3 4 5 x pull -5-4 -3-2 -1 1 2 3 4 5 y pull (c) (d).35.3.25 Full: <>=.88 RMS=1.385 Simplified: <>=.76 RMS=1.3165.35.3.25 Full: <>=.39 RMS=1.3281 Simplified: <>=.19 RMS=1.348.2.2.15.15.1.1.5.5-5 -4-3 -2-1 1 2 3 4 5 t x pull -5-4 -3-2 -1 1 2 3 4 5 t y pull (e).3 Full: <>=.489 RMS=1.3979.25 Simplified: <>=.1184 RMS=1.4124.2.15.1.5-5 -4-3 -2-1 1 2 3 4 5 1/p pull Figure 4: Pull distributions of the track parameters at the first measurement point for the full and the simplified geometries. This sample of long tracks was obtained with the Forward pattern ognition algorithm.
(a).1 Full: <>=-.1 RMS=.113 Simplified: <>=-. RMS=.113.8.1 Full: <>=-.8 RMS=.1275.8 Simplified: <>=-.6 RMS=.1278.6.6.4.4.2.2 -.1 -.8 -.6 -.4 -.2.2.4.6.8.1 x - x (cm) -1 -.8 -.6 -.4 -.2.2.4.6.8 1 y - y (cm) (c) (d).18.16.14.12.1 Full: <>=. RMS=.6 Simplified: <>=-. RMS=.6.1.8.6 Full: <>=-. RMS=.9 Simplified: <>=-. RMS=.9.8.6.4.4.2.2 -.5 -.4 -.3 -.2 -.1.1.2.3.4.5 t x() - t x() (rad) -.5 -.4 -.3 -.2 -.1.1.2.3.4.5 t y() - t y() (rad) (e).1.8 Full: <>=.2 RMS=.53 Simplified: <>=.2 RMS=.55.6.4.2 -.5 -.4 -.3 -.2 -.1.1.2.3.4.5 (p - p )/p Figure 5: Resolutions on the track parameters in the Outer Tracker region for the full and the simplified geometries. This sample of long tracks was obtained with the Forward pattern ognition algorithm.
(a) )/p.8 Full Simplified p.7 (p.6.5.4 2 2.5 3 3.5 4 4.5 5 η η.12 η.1.8.6.4.2 2 2.5 3 3.5 4 4.5 5 η (c).2 )/p.18 p.16 (p.14.12.1.8.6.4 2 4 6 8 1 12 14 16 18 2 p (GeV) (d) 2 Full: <>=-.1 σ=.437 Simplified: <>=-.3 σ=.45 18 16 14 12 1 8 6 4 2 -.1 -.8 -.6 -.4 -.2.2.4.6.8.1 (p - p )/p Figure 6: Resolutions in (a) momentum versus pseudorapidity, in pseudorapidity versus pseudorapidity and in (c) momentum versus momentum, as given by the σ values of single- Gaussian fits. Figure (d) shows the result of a single-gaussian fit to the momentum resolution averaged over the momentum range in (c). All the distributions were obtained for the full and the simplified geometries with the Forward pattern ognition algorithm.
Figure 4 shows the pull distributions at the first track measurement point. No differences are observed between the full and the simplified geometries apart from a slight increase in the momentum bias; it increases from.5±.6 to.12±.6. The exercise was repeated at different locations along the track s trajectory: resolutions and pull distributions were calculated at the track s origin vertex position and at positions in the various tracking detectors. One such example of resolution distributions in the region of the Outer Tracker is presented in Figure 5. In all cases the same conclusions can be drawn as for the distributions at the first measurement point. The momentum resolution was also studied as a function of the momentum and the pseudorapidity of the tracks. The distributions collected in Figure 6 profile the σ values of single- Gaussian fits to the momentum resolution. A small deterioration can be observed over the full momentum and η spectra. The momentum resolution versus momentum was projected onto the y-axis to obtain an average resolution over(most of) the spectrum; the results of single-gaussian fits to the obtained projection distributions (Figure 6(d)) show a core momentum resolution of.44% and.45% with the full and the simplified geometry, respectively. Figure 6 further shows the tracks pseudoradipity resolution as a function of pseudorapidity. No degradation of resolution was observed..16.14.12 Full: <>=.3367 RMS=.323 Simplified: <>=.3416 RMS=.355.1.8.6.4.2.1.2.3.4.5.6.7.8.9 1 Probability of χ 2 Figure 7: Track fit probability of χ 2 distribution for long tracks found by the Forward pattern ognition algorithm for the full and the simplified geometries. The probability of track fit χ 2 is another measure of the fit quality; an accurate fit model shouldgiverisetoaflatdistribution(adiscussion canbefoundin[5]). Figure7comparestheχ 2 probability distribution obtained with the full and the simplified geometries. Both distributions agree very well.
4 Physics analysis In this section the impact of using the simplified geometry for track fitting on the quality of the selection and onstruction of B decays is studied. The B π + π decay (extensively described in [6]) was used for the sake of example. 4.1 Effect on the event selection B decays are typically selected exploiting the high mass and long lifetime of B mesons. The discriminating variables used are the transverse momentum and impact parameter of the B and its daughters and the flight-distance of the B. In Table 2 the selection cuts applied to the generic 2 hh channels are shown (a more detailed explanation of all cuts can be found in [6]). 2 selection parameter cut value hh smallest p t (GeV) of the daughters > 1. largest p t (GeV) of the daughters > 3. B(s) p t(gev) > 1.2 smallest IP/σ IP of the daughters > 6 largest IP/σ IP of the daughters > 12 B(s) IP/σ IP < 2.5 B(s) vertex fit χ2 < 5 (L)/σ L > 18 m (MeV) < 5 Table 2: Selection cuts applied to the 2 hh channels [6]. The distributions of the various 2 hh selection variables are shown in Figures 8 to 1. Positively and negatively charged pions were looked at independently to track down any possible charge-induced biases. Note that all plots were obtained after applying the full selection on all the variables but the plotted one. In case of the p T and impact parameter significance cuts on the pions, where one threshold is applied to both pions and another has to be exceeded by at least one of them, these cuts have been switched off simultaneously. In addition, normalised integrals are shown to give a dit comparison of the acceptances obtained with the full and the simplified setups. The distributions obtained with the simplified geometry agree very well with those obtained with the full setup. The observed differences were at most at the percent level. As the differences in the single-cut variables partly cancel out, the overall change in the number of selected 2 hh events is 1% (Table 3). Still, as can be seen from the table, the percentage of common events selected with the full and the simplified geometries is 95-96%, whereas 4-5% of the events in each sample are only present in that particular sample. In the following all comparison distributions were obtained with all the selected events, i.e. using both common and non-common events. 4.2 Effect on resolutions Finally, the resolution on the most important physics analysis observables were compared: momentum resolutions have been studied as well as the resolutions on the primary and secondary
25 2 (a) acceptance 1.8 15.6 1.4 5.2 1 2 3 4 5 6 7 8 9 1 IP significance B 1 2 3 4 5 6 7 8 9 1 B IP significance 14 12 (c) acceptance 1.8 (d) 1 8.6 6.4 4 2.2 1 2 3 4 5 6 7 8 9 1 π + IP significance 1 2 3 4 5 6 7 8 9 1 π + IP significance 12 1 (e) acceptance 1.8 (f) 8.6 6 4.4 2.2 1 2 3 4 5 6 7 8 9 1 - π IP significance 1 2 3 4 5 6 7 8 9 1 - π IP significance Figure 8: Distributions of the event selection variables of impact parameter significances for (a) the B candidate and its daughter (c) positively and (e) negatively charged pions for the full and simplified geometries(full and dashed lines, respectively). The right-hand-side distributions correspond to the integrated left-hand-side distributions.
3 25 2 (a) acceptance 1.8 15.6 1.4 5.2 1 2 3 4 5 B p (MeV) T 1 2 3 4 5 B p (MeV) T 6 5 4 (c) acceptance 1.8 (d) 3.6 2.4 1.2 1 2 3 π + p (MeV) T 1 2 3 π + p (MeV) T 5 4 (e) acceptance 1.8 (f) 3.6 2.4 1.2 1 2 3 π - p (MeV) T 1 2 3 π - p (MeV) T Figure9: Distributions of the event selection variables of (a) B and daughter (c) positively and (e) negatively charged pions transverse momentum for the full and simplified geometries (full and dashed lines, respectively). The right-hand-side distributions correspond to the integrated left-hand-side distributions.
25 2 (a) acceptance 1.95.9.85 15.8.75 1.7 5.65.6 5 1 15 2 25 3 35 4 45 5 B decay vertex χ 2 5 1 15 2 25 3 35 4 45 5 B decay vertex χ 2 16 14 12 (c) acceptance 1.8 (d) 1 8.6 6.4 4 2.2 1 2 3 4 5 6 7 8 9 1 B flight distance significance 1 2 3 4 5 6 7 8 9 1 B flight distance significance 3 25 (e) 2 15 1 5 51 515 52 525 53 535 54 545 55 Mass (MeV) Figure 1: Distributions of B (a) decay vertex χ 2, (c) flight distance significance and (e) invariant mass for the full and simplified geometries (full and dashed lines, respectively). The right-hand-side distributions correspond to the integrated left-hand-side distributions.
Number of selected Events only In common Geometry events in the sample with other sample full 4141 162 (3.9%) 3979 (96.1%) simplified 4186 27 (4.9%) 3979 (95.1%) Table 3: Number ofselected events afterrunning the2 hhselection forthefull andthe simplified geometries. The third and the forth rows indicate, respectively, the number (and percentage) of selected events only present in the full and in the simplified samples and the number of events in common. (B decay) vertices and on the B proper time. These resolutions are shown in Figures 11 and 12 while their values (the σ s of single-gaussian fits) are summarised in Tables 4 and 5. No significant degradation on either of the quantities was observed. Comparing with the momentum resolution quoted in Section 3, the numbers in Table 4 differ by roughly 2%. This difference is understood as the latter numbers correspond to the σ values of single-gaussian fits (rather than the width of the distribution) for the sub-sample of relatively high momentum B-daughter pions. Momentum Mass Proper time Geometry resolution resolution resolution (%) (MeV) (fs) full.495(5) 22.5(3) 37.7(5) simplified.52(6) 22.9(4) 37.7(6) Table 4: Values of the resolutions of the daughter pion momenta, the B mass and the B proper time for the full and the simplified geometries. The resolutions correspond to the σ values of single-gaussian fits. The errors on the last digit are specified in parenthesis. Primary vertex B vertex Geometry resolutions (µm) resolutions (µm) x y z x y z full 9.2(1) 8.8(2) 41.4(7) 14.2(2) 14.(2) 147(3) simplified 8.9(1) 8.8(1) 41.4(7) 14.3(2) 14.3(2) 145(3) Table 5: Values of the position resolutions on the primary and the B decay vertices for the full and the simplified geometries. The resolutions correspond to the σ values of single-gaussian fits. The errors on the last digit are specified in parenthesis. Themomentumandxandy slopesofthepositivelyandnegativelychargedbdaughterpions were further investigated as a function of the track azimuthal angle φ. The detector geometry in φ is highly non-trivial, which makes the simplified description a potentially inappropriate
35 35 3 25 (a) 3 25 2 2 15 15 1 1 5 5 -.1 -.5.5.1 - p )/p (p -.1 -.5.5.1 - p )/p (p 9 8 7 (c) 6 5 4 3 2 1 -.5 -.4 -.3 -.2 -.1 -.1.2.3.4.5 t - t (ps) Figure 11: Resolutions on the (a) positively and the negatively charged daughter pion momenta and on the (c) B proper time for the full and simplified geometries (full and dashed lines, respectively). replacement. As can be seen from Figure 13, no significant disagreement was found (in spite of low statistics). Additionally, a dit comparison of the B daughter pion momenta as onstructed with the full and simplified geometries has been made. Figure 14 shows the relative difference of the onstructed momenta for positively and negatively charged pions. A single-gaussian fit gives a σ of.6% without any significant bias for both distributions. 5 Conclusions and Final Remarks The alternative simplified geometry for track fitting has been validated with respect to the full detector geometry. The implications for physics analysis in terms of tracking and physics
Figure12: Resolutionsonthe(a)primaryvertexandtheB vertexforthefullandsimplified geometries (full and dashed lines, respectively). 12 1 8 6 4 2 9 8 7 6 5 4 3 2 1 (a) -.2 -.15 -.1 -.5 -.5.1.15.2 x vtx() - x vtx() (mm) -.2 -.15 -.1 -.5 -.5.1.15.2 B decay vertex x - x (mm) 14 12 1 8 6 4 2 9 8 7 6 5 4 3 2 1 -.2 -.15 -.1 -.5 -.5 - y.1.15 (mm).2 y vtx() vtx() -.2 -.15 -.1 -.5 -.5.1.15.2 B decay vertex y - y (mm) 8 7 6 5 4 3 2 1 -.5 -.4 -.3 -.2 -.1 -.1.2.3.4.5 z vtx() - z vtx() (mm) 25 2 15 1 5 -.5 -.4 -.3 -.2 -.1 -.1.2.3.4.5 B decay vertex z - z (mm)
)/p - p (p.65.6.55.5.45 (a) -3-2 -1 1 2 3 φ )/p - p (p.58.56.54.52.5.48.46.44.42.4.38.36-3 -2-1 1 2 3 φ x() - slope.14.13-3 1 (c) y() - slope.11-3 1 (d) x().12 y().1 slope.11.1 slope.9.9.8.8.7.6.7-3 -2-1 1 2 3 φ -3-2 -1 1 2 3 φ x() - slope x() slope.13.12.11.1.9.8.7-3 1 (e).6-3 -2-1 1 2 3 φ y() - slope y() slope.13.12.11.1.9.8.7.6-3 1 (f) -3-2 -1 1 2 3 φ Figure 13: Resolutions of negatively charged (left-hand-side distributions) and positively charged (right-hand-side distributions) daughter pions in momentum (a,b) and slopes in x (c,d) and y (e,f) as function of φ for the full and simplified geometries (full and dashed lines, respectively).
N events 35 3 (a) N events 35 3 25 25 2 2 15 15 1 1 5 5 -.4 -.2.2.4 (p -p )/p simple full full -.4 -.2.2.4 (p -p )/p simple full full Figure 14: Relative difference in the B daughter pion momenta between full and simplified geometry for positive (a) and negative tracks. performance were assessed. No significant degradation of performance was found in this study of 2 hh events. With the LHC start-up date approaching, it is foreseen to onstruct 28 data with the full detector geometry description for track fitting. A new simplified description will then be re-derived at a later stage, which in turn will need to be validated again before the decision to switch to the simplified description can be taken. From this study no major problem is foreseen.
References [1] W. Hulsbergen, http://indico.cern.ch/conferenceotherviews.py?confid=1171. [2] D. Hutchcroft et al., VELO Pattern Recognition, LHCb Note 27-13, 27. [3] O. Callot, S. Hansmann-Menzemer, The Forward Tracking: Algorithm and Performance Studies, LHCb Note 27-15, 27. [4] J. van Tilburg, Matching VELO tracks with seeding tracks, LHCb Note 21-13, 21. [5] M. Needham, Performance of the LHCb track onstruction software, LHCb Note 27-144, 27. [6] A. Carbone et al., Charmless charged two-body B decays at LHCb, LHCb Note 27-56, 27; A. Carbone et al., Analysis of two-body B decays at LHCb using DC6 data, LHCb Note in preparation.