Exclusive Reconstruction of B-Decays with Missing Neutrals

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 4momentum conservations in vertices with available geometric information in the event. The paper gives details of such methods on the B0 s → D− sK , D− s → 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.


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 Loc atio n fro m 4-m om ent um con ser vat ion In underground experiments the geometric equivalent to accelerator B-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 mode 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 2fold 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 fieldsuch 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 O III helical correction, giving a solution also for other detectors with full-strength magnetic field in the vertexing region .
Vertexing methods reject reasonably well QCD backgrounds originating at the interaction point IP , with some resolution there around ; however, other B-modes can constitute significant backgrounds.One such mode is the sister mode with pion in the Bvertex, roughly 4-5 times more abundant.Kinematically and topologically similar modes from B 0 d also exist-the latter roughly 3 times more abundant than B 0 s .Although there is no exact replica of the mode in B 0 d 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 B 0 d being less "kaon-productive" than B 0 s and that there are two kaons present in the D-vertex.All other kinematic parameters offer little protection against B 0 d 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 B 0 d originating modes ca.45% rejection .

Dynamic Vertexing
The 4-momentum conservations in the B-and D-vertices are 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
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 E D with some coefficients: , resulting in a second order equation in E D :

2.4
For convenience Q has the direction of Q in the figure and the absolute value of 1 GeV/c.

Advances in High Energy Physics
The ambiguity raised by the second order equation is solved in favor of the " " sign.What actually happens is that β Q b p k is rather small and the equation looks like The solution with "−" sign is therefore negative in value, hence unphysical.Once E D 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 − λ K n K where n K p K /p K , 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 This implies minimising the above 3 constraints simultaneously: 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
Minimising with respect to the B-vertex B, the following is obtained: 2.9 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 errorellipsoid operator reduces to the perpendicular to track projector.Therefore in q only the perpendicular to track error of r 0K will contribute to the B-vertex error.
Similarly, for spherical error-ellipsoids, the D-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 energy bb 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/c 2 , 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 d flight M B /cp B .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.
The obtained resolution is similar to usual resolutions in B-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 of J/ψ 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 B 0 d → 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 M lν 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., −M lν l and so the rejection power should be better especially for the electron, where the mass is also much lighter .

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.

Figure 1 :
Figure 1: Topology of the mode considered in the absence of magnetic field: B 0 s → D − s K K − π − π 0 K .The line Q is the perpendicular from the D-vertex onto the kaon track.

Figure 2 :
Figure 2: Residuals of the reconstructed B-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.