# Exclusive Reconstruction ofB-Decays with Missing Neutrals

Authors: M. Dima
Journal: Advances in High Energy Physics
 (2012)

---

## Abstract

Often decay channels that are of theoretical interest cannot be reconstructed exclusively due to missing neutrals (such as neutrinos), or due to single-track vertices. This situation appears both in underground astrophysics experiments as well as in conventional accelerator experiments. A method to “recover” such missing particles from their kinematics and reconstruct “exclusively” the modes would benefit both domains in a number of ways. The main idea is to combine 4-momentum conservations in vertices with available geometric information in the event. The paper gives details of such methods on theBs0→Ds−K+, Ds−→K+K−π−(π0) prototype decay, which also encounters 2-fold ambiguities in its solutions. Such ambiguities can be lifted and the paper shows how, while also addressing the potential the method has in physics analyses and detector studies.


---

## Body

## 1. Introduction

The kinematics of decay modes with a missing particle, usually an elusive neutrino, are of interest in both underground astrophysics experiments as well as in conventional accelerator experiments. Exclusive reconstruction in today’s new generation experiments is important for precision physics, performant background rejections needed in rare decays, and so forth. Some modes that are of considerable theoretical interest, from the experimental point of view, pose inherent difficulties with missing neutrals and single-track vertices. While such is not an impossible situation (since the vertex must be somewhere on the track), the precise location eludes in general the reconstruction of the mode. Nonetheless, in certain cases, information from the rest of the event provides this capability—in general in the form of 4-momentum conservation in subsequent vertices. Work on recovery using dynamic information has been explored before [1, 2], not necessarily in association with vertexing. Code addressing such methods stems from the days of the CLEO collaboration, KWFIT [3, 4] and of BaBar Collaboration, KinFitter [5].In underground experiments the geometric equivalent to acceleratorB-physics programs are strange-particle decays [6, 7], where kinematic imprecisions associated with the missing neutrino and a traditional lack of some form of vertexing have prevented exclusive reconstruction of the modes. At the new underground astrophysics experiments (such as LAGUNA [8–12]) tracking-vertexing facilities are being discussed in some of the studies. The methods studied could in principle be applied for all “stealth” particles (neutrino, π0, other neutrals) decaying further downstream.The objective in such cases is to recover if possible exclusively the event, by combining conservation equations and geometric aspects in the decays. The case made in this paper is for tracking-vertexing assemblies in underground astrophysics experiments and for more precise vertexing in accelerator-based experiments.Consider the prototype modeBs0→Ds-K+, Ds-→K+K-π- (π0), which is of interest for γCKM measurement, however, with a somewhat fair branching ratio. With π0 in the final state the branching ratio is roughly equal to that without the π0 [13], thus putting the two together essentially doubles the statistics. Reconstruction of the π0 would in principle be possible in the calorimeter, however, with less resolution than using tracking.Fortunately there is sufficient information in the event for reconstructing the mode without explicit calorimeter information on theπ0—in principle with better resolution also. This information is from 4-momentum conservation in the B- and D-vertices. A particular 2-fold ambiguity from the quadratic equations arises in this (and other) modes. It is important to be able to lift this ambiguity without imposing cuts (which due to detector resolution can eliminate signal). Figure 1 shows the topology of the event (in the absence of magnetic field—such as for the LHCb experiment at the LHC. A small track bending over the vertexing region can also be included, in the form of an 𝒪III helical correction, giving a solution also for other detectors with full-strength magnetic field in the vertexing region).Figure 1
Topology of the mode considered in the absence of magnetic field:Bs0→Ds-(K+K-π-π0) K+. The line Q→ is the perpendicular from the D-vertex onto the kaon track.Vertexing methods reject reasonably well QCD backgrounds (originating at the interaction point (IP), with some resolution there around); however, otherB-modes can constitute significant backgrounds. One such mode is the sister mode with pion in the B-vertex, roughly 4-5 times more abundant. Kinematically and topologically similar modes from Bd0 also exist—the latter roughly 3 times more abundant than Bs0. Although there is no exact replica of the mode in Bd0 version, there are many modes with nonnegligible branching ratios that “losing” one track can mimic this mode. Due to the similarities in topology (displaced vertices) and kinematics (B-masses, center of mass decay energies), the method cannot distinguish between such cases—not with current detector resolution. The only remaining solution is performant Particle ID, making sure that the B-vertex track is a kaon (Bd0 being less “kaon-productive” than Bs0) and that there are two kaons present in the D-vertex. All other kinematic parameters offer little protection against Bd0 modes—due to the high boost factors, mass similarities, and center of mass decay energies. A pleiad of such parameters has been studied: the discriminant of the second order equation, pointback, reconstructed invariant masses, angles between track-combinations or decay-planes, 2-3 parameters combinations of the aforementioned, and so forth. For all the detector resolution was insufficient to tell apart the current mode from its sister mode (best achievable ca. 41% rejection), or from other Bd0 originating modes (ca. 45% rejection).


## 2. Dynamic Vertexing

The 4-momentum conservations in theB- and D-vertices are
(2.1)(EK+ED)2=MB2+(pD2+pK2+2p→Dp→K),(ED-Ev)2=mπ2+(pD2+pv2-2p→Dp→v),
where v refers to the “virtual” particle having 4-momentum equal to the sum of the charged tracks in the D-vertex. These can be rewritten as
(2.2)p→Dp→K=MK2+MD2-MB22+EKED,p→Dp→v=-MD2+Mv2-mπ22+EvED.
Since p→B=p→D+p→K, the three vectors lie in the same plane. Also, the B-vertex is on the kaon track, hence p→D=λQ→+μp→k, the two constants λ and μ being determined from (2.1) as linear functions of ED with some coefficients:
(2.3)λ=αED+β,μ=aED+b.
However, pD2=ED2-MD2=(λQ→+μp→K)2, resulting in a second order equation in ED:
(2.4)ED2-MD2=ED2(αQ→+ap→K)2+ED·2(αQ→+ap→K)(βQ→+bp→K)+(βQ→+bp→K)2.
For convenience Q→ has the direction of Q→ in the figure and the absolute value of 1 GeV/c.The ambiguity raised by the second order equation is solved in favor of the “+” sign. What actually happens is that(βQ→+bp→k) is rather small (ϵ→) and the equation looks like
(2.5)ED2-MD2=χ→2ED2+ϵ→(2χ→ED+ϵ→).
The solution with “−” sign is therefore negative in value, hence unphysical.OnceED is known, the coefficients λ and μ are known and so are p→D and p→B.Having the D-vertex and-p→D pointing back at the kaon track, a B-vertex can be now inferred. The problem is formulated as follows: the B-vertex is a point inside the 1σK kaon track error-tube  with the property that aiming along the p→D direction, the D-vertex is missed by about 1σD, while aiming along the -p→B direction, the IP is missed by about 1σIP. The kaon track is parametrised as r→0K-λKn→K where n→K=p→K/pK, r→0K, some origin on the track, and λK, a running parameter describing the track. Therefore, if B→ is the B-vertex, the above requirements translate into(1)
kaon track tube,σK, from -r→0K+λKn→K+B→,(2)
D-vertex ±σD, from B→+λDp→D,(3)
IP±σIP, from B→+λBp→B.This implies minimising the above 3 constraints simultaneously:
(2.6)〈B→+λKn→K-r→0K|σK-2|B→+λKn→K-r→0K〉+〈B→+λDp→D-v→D|σD-2|B→+λDp→D-v→D〉+〈B→+λBp→B-v→IP|σIP-2|B→+λBp→B-v→IP〉
with respect to the three running parameters λK (units [m]), λD, and λB (both in units of [m]/[GeV/c]) and to B→. Annulling the respective derivatives:
(2.7)λK=〈B→-r→0K|σK-2|n→K〉〈n→K|σK-2|n→K〉,λD=〈B→-v→D|σD-2|n→D〉〈n→D|σD-2|n→D〉,λB=〈B→-v→IP|σIP-2|n→B〉〈n→B|σIP-2|n→B〉.Minimising with respect to theB-vertex B→, the following is obtained:(2.8)|B→〉=[∑x=K,D,B‍σx-2(1-|n→x〉〈n→x|σx-2〈n→x|σx-2|n→x〉)]-1|q→〉
with
(2.9)|q→〉=∑x=K,D,B‍σx-2(1-|n→x〉〈n→x|σx-2〈n→x|σx-2|n→x〉)|v→x〉.For the kaon track error tube an extended acception of error ellipsoid was used: one with infinite long-axis. In the case of a round cross-section (equal axes) the error-ellipsoid operator reduces to the perpendicular to track projector. Therefore (inq→) only the perpendicular to track error of r→0K will contribute to the B-vertex error.Similarly, for spherical error-ellipsoids, theD-vertex and the IP would contribute only with their respective errors perpendicular to n→D and n→B. This is not exactly the case, both ellipsoids being elongated, the IP along z and the D-vertex along its flight direction. This extra elongation simply translates into the addition of a parallel projector along said directions to the respective error-ellipsoid operators. For the D-vertex the elongation being exactly along n→D and greater than the spherical part of the operator, the combined effect in q→ is the same as having only the spherical part. For the IP the elongation is along z, which is not n→B, but almost. So again, to first approximation, only the waist of the IP will play a determinant part in the B-vertex errors. The IP being very accurate waist-wise (for the LHCb-IP: σxy=20 μm versus σz=60 μm [14]), the main contribution will come from the D-vertex—set by smearing to a typical resolution of this detector [15].The code was tested using a PYTHIA simulation of the LHC environment for the LHCb experiment (LHC energybb- events, “HardQCD:gg2bbbar”), for which the tracks were smeared both in momentum value and direction (like a needle deviating in direction within an error-tube). The momentum value smearing was chosen such as to give 2-body mass resolutions as those reported by this experiment in literature [16, 17] (ca. 12 MeV/c2), while the error-tube diameter (ca. 15 μm) was chosen such as to give the vertex locator resolutions reported in the literature [17] (ca. 150 μm longitudinally and 25 μm transversally). Once the particle tracks and vertices are obtained in a form simulating detector response, they are fed to the code for reconstruction. Figure 2 shows the reconstructed B-vertex resolution in terms of proper time, determined from τproper=dflightMB/cpB. The physically relevant directions are parallel to flight direction and transverse to this. The method uses 3 constraints combined, hence a rudimentary approximation is to allot a gaussian in the resolution fit for each. The top part of the figure shows the residuals with respect to flight-direction (which is almost z-aligned, fitted with 3 Gaussians). The combined resolution is on the order of 32 fs. The bottom part shows residuals transverse to the flight direction (approximately the xy-plane) fitted with 3 half-Gaussians centered at zero, resolution on the order of 5.4 fs. It is evident that the resolution is dominated by the boost direction (flight direction elongated error-ellipsoid). This is to be expected due to the high boosts, where the shallow opening angles among tracks translate pixel detector errors into large longitudinal ones.Residuals of the reconstructedB-vertex with respect to the Monte Carlo position: (a) projected onto the B-flight direction, (b) transverse to the B-flight direction. The dimensions are proper time in fs. Three Gaussians were used in the fits (in the case of the transverse centered at zero), the global resolutions being ca. 32 fs (longitudinally) and 5.4 fs (transversally)—details are in the text.
(a)
(b)The obtained resolution is similar to usual resolutions inB-reconstructions achieved in 4 track decays [18] (40 fs). The obvious advantage here is that although particles are missing, their kinematics helps “recover” them as if they were present and this is an exclusive reconstruction.The method can also be used for detector studies, such as calibrating Particle ID, or calorimetry, in the lepton (electron/muon) sector. This implies in general using semileptonic modes. Although the traditional leptonic decay ofJ/ψ is standard, it yields two hard tracks, quite close to each other (which in the electron case makes difficult estimating the individual radiative losses). Also Particle ID performs lesser on tracks very close to one another, the real performance (used in B-physics) being actually better. Therefore, the larger variety given by semileptonic modes would be preferable, if only a pure enough (better than Particle ID miss-ID rates) sample could be selected. Consider, for instance, the mode Bd0→D*- (K-π+l-ν-l) π+π-π+. The B-vertex provides 3 “clean" pion tracks useful for Particle-ID studies in the hadron sector. The D-vertex would also be good, if only labels could be attached to the tracks. Owing to the large mass difference between kaon/pion and electron, the equations would prevent a vertexable solution if the masses were wrongly attached. However, in the case l=μ, the pion and lepton tracks are similar in mass. Still, they are quite different in dynamics, one coming from a lepton line, united to the baryon half of the reaction via a W-, while the other from a continuous quark line. The imbalance in this case would come from the different momentum properties (putting the wrong momentum in Mlνl, for instance, versus other dynamic quantities in the reaction). In this sense the method can act as a track-labeler in the D-vertex, offering a quite clean sample of B-origin leptons for Particle-ID studies—which is rare. A study was performed for kaon/pion separation power of the method. The mode with no π0 in the final state was chosen. Of the existing tracks in the D-vertex one was purposely reported “lost” to the program and its mass fed to the code. If the kaon track is reported “lost,” but pion mass is given for it to the code, the inner kinematics of the reactions will work out wrong, further influencing the (highly nonlinear) vertexing procedure described aforehand. The vertexing produces now new momenta for all particles implied and an invariant mass can be recomputed for the “lost” particle’s mass (different, but close, to the one inputted). This mass error versus pointback is a criterion that gives a ca. 45% kaon-pion separation. As pointed above, for leptons additional handles (due to inner dynamics) exist (i.e., -Mlνl) and so the rejection power should be better (especially for the electron, where the mass is also much lighter).


## 3. Conclusions

The method aims to recover modes otherwise nonusable, with missing neutrals lost (neutrino), or less precise (calorimetry). This is very important to both underground based experiments as well as accelerator-based experiments, where background rejection or measurement precision depends on the exclusive reconstruction of the event. The paper lifts the two-fold ambiguity arising in the second order equation involved, which otherwise relied on other criteria (pointback for instance) which cut on signal.The potential of the method for Particle-ID source studies (lepton sector) is also discussed.The methods studied could in principle be applied for all “stealth” particles (neutrino,π0, other neutrals) decaying further downstream. Plans for extending the method could involve these options together with the Particle-ID studies.



---

*Source: 123083-2012-12-19.xml*