arxiv: v1 [hep-ph] 8 May 2013

Transcription

1 Analtic solutions for neutrino momenta in deca of top quarks at hadron colliders Burton A. Betchart, Regina Demina, Amnon Harel Universit of Rochester arxiv: v1 [hep-ph] 8 Ma 013 Ma 9, 013 Abstract We emplo a geometric approach to analticall solving equations of constraint on the deca of top quarks involving leptons at hadron colliders. Neutrino momentum is found as a function of the 4-vectors of the associated b quark and charged lepton, the masses of the top quark and W boson, and a single parameter, which constrains it to an ellipse. The measured imbalance of transverse momenta in the event helps to reduce the solution for neutrino momenta to a discrete set for the cases of one or two top quarks decaing to leptons. The solutions can be implemented concisel with common linear algebra routines. 1 Introduction Top quark reconstruction at hadron colliders from channels containing one or more top quark decas to leptons presents a challenge since the longitudinal momentum of the colliding sstem is unconstrained and the neutrinos are not directl observed. The sum of neutrino momenta in the plane transverse to the beam ais is inferred from the total momentum imbalance. For a single neutrino final state, constraining the invariant mass of the neutrino and associated charged lepton to the W boson mass provides a quadratic equation for the unmeasured longitudinal component of neutrino momentum, with zero solutions, or two solutions which can be further resolved b application of additional invariant mass constraints. Since the total momentum imbalance frequentl has the worst resolution of all constraints on the top quark deca, we suggest an alternative approach to analtic top quark reconstruction in which the invariant mass constraints from top quark and W boson are both eact, and in which the solution set for each neutrino momentum is an ellipse. For events with a single neutrino, the ellipse is analticall reduced to a unique solution b application of the transverse momentum imbalance constraint, taking its uncertaint into account. The approach etends naturall to the case of two neutrinos in the final state, allowing an alternative method for calculating the solution pairs previousl described b [1,, 3], and suggests a most likel pair in the case of no eact solution. The methods described perform well for simulated Tevatron and LHC events. These results ma also be useful in the case of top quark deca to three quarks, where one of the quarks from the deca of the intermediate W boson falls outside eperimental acceptance. Derivation The kinematics of top quark deca constrain the W boson momentum vector to an ellipsoidal surface of revolution about an ais coincident with bottom quark momentum. Simultaneousl, the kinematics of W boson deca constrain the W boson momentum to a paraboloidal surface of revolution about an ais coincident with the momentum of the resulting charged lepton. The intersection of the two surfaces is an ellipse. The neutrino momentum is consequentl constrained to a translation of the ellipse, for which a parametric epression in the laborator frame is given. Measurement of the transverse momentum imbalance further constrains the cases of one and two neutrinos from a pair of top quarks) to a discrete solution set..1 Definitions A particle p is described b its mass P m, energ P e, and momentum 3-vector P p, with the dispersion relation P m = P e P p. The azimuthal and polar angles of p in the laborator frame are denoted P φ and P θ. For definiteness we consider the deca chain t bw bµν; a top quark T ) decas to a bottom quark B) and a W boson W ), with 1

2 z, z W θ BW µ θ µw b θ Bµ Figure 1: The momenta of observed particles µ and b define the coordinate sstems F and F. A possible momentum of the W boson is drawn to show the angles θ BW and θ µw. subsequent deca of the W boson to a muon µ) and a neutrino ν). We assume established masses T m, B m, W m, for eample from world average measurements[4], and negligible mass for the neutrino and muon. Energ and momentum conservation for this sstem impl T e = + W e = + µ e + ν e T p = B p + W p = B p + µ p + ν p. Momentum coordinate sstems F =, ỹ, z) and F =, ỹ, z ) are defined to share a common ais z = z, but F is rotated relative to F b the angle θ Bµ between B p and µ p, with µ p along the -ais, and B p along the -ais. The relationships between Cartesian and polar coordinates of W p in F ) are ) = W pc ), ỹ ) + z ) = W ps ), 1) where C ) is the cosine and S ) is the sine of the angle θ µw the coordinate sstems. θ BW ) between W p and µ p B p). Figure 1 shows. Two surfaces for W momentum Energ and momentum are conserved in the deca t bw, hence Tm = Te Tp = + W e) B p + W p) = B m + W m + W e B pw pc. For compactness, define Then 1 T m Wm Bm ). ) = W e B pw pc [ + B pw pc ] = Be W m + Wp ) 0 = B ew m 4) W p B pc ) + W p 0 = B ew m 4) W p B pc ) + W p B e B pc ) B mc + B es ).

3 Incorporating relations 1) for the F coordinates, it is clear that W p is constrained to the surface B m) + ỹ ) + z ) B p + B ew m 4) = 0, 3) which is an ellipsoid of revolution about the -ais. Particle W subsequentl decas to particles µ and ν. This deca has the same kinematics as the deca of T, with the substitutions B µ, W ν, T W. Using F coordinates, ν p is constrained to the surface ỹµ p) + zµ p) W mµ p W 4 m 4 which is a paraboloid of revolution about the -ais. A congruent surface of solutions for W p is shifted from the neutrino solutions b µ p in, ỹ + z W m 0 = 0, 4) where for convenience we divided b µ p and defined = 0. 0 W m µ p. 5) Surfaces 3) and 4) are simultaneous constraints on W p, so the solution set is their intersection. Figure shows eamples..3 Etended Matri Representation with Homogeneous Coordinates The use of homogeneous coordinates r = 1 ) T or s = z 1 ) T allows etended matri representation of various - or 3-dimensional geometric objects. The 1 3 row matri L is the etended representation of the line Lr = 0. A 3 3 smmetric square matri M is the etended representation of the conic section r T Mr = 0. A quadric surface, like a paraboloid or ellipsoid, can be represented as a 4 4 smmetric square matri A, with s T As = 0. Etended representations are unique up to a multiplicative factor, and allow both rotations and translations to be epressed b matri multiplication[5, 6, 7]. The paraboloid 4) is represented for homogeneous F coordinates with the matri à µ = Wm 0. 6) The ellipsoid 3), divided b B e, is represented for homogeneous F coordinates with the matri à b = ) B m Be 0 0 B p Be B p B e 0 0 W m ). 7) The transformation from homogeneous coordinates F to F is given b a rotation of θ Bµ around the z ais, c s 0 0 K = s c , c = cos θ Bµ 8) s = sin θ Bµ The ellipsoid 3) is represented in homogeneous F coordinates with the matri Ãb = Kà bk T, ) ) c B m Be + s Bp cs 0 c B p B ) e ) Bp cs à b = s B m Be + c 0 s B p Be Wm c B p B e s B p B e ). 9) 3

4 F coordinates: slice of z = z = 0 +P 1, 1,0) z, z B Lab coordinates: projections on z=0 µ P 0, 0) χ = 4.6 F coordinates: slice of z = z = 0 +P B 1, 1,0) z, z Lab coordinates: projections on z=0 µ P 0, 0) χ = 6.3 Figure : Reconstruction of neutrino momentum from the deca t bµν for two events, with θ Bµ large top) and small bottom). Left panel arrows indicate momenta of µ black) and b gra). Corresponding constraints on the intermediate W boson momentum are indicated dashed) b their z = 0 cross sections. The degenerate quadric surface P dotted) contains the intersection of the constraints. Upper panel ellipses show the intersection in the plane + P. Lower panel filled arrows indicate laborator transverse momenta of µ black), b gra), and ν dashed). The projection of the intersection ellipse from + P is traced in light gra dashes, while the corresponding ellipse of neutrino transverse momenta solutions is solid black. The measured imbalance in total transverse momentum and its uncertaint are shown b and the encircling dots. The neutrino transverse momentum solution most compatible with the measurement is shown b the unfilled arrow. All panels are marked in units of 100 GeV/c. 4

5 Case Factorization horizontal & vertical, G 11 = 0 = G L + = G 1 0 G 3 ), L 0 G 1 G 13 G 3 ) parallel, g 33 = 0 G L ± = G 1 G G 3 ± ) ) ) g 11 intersecting, G 0 0 = g3 g 33 G + g13 g 33 G 1 ± g 33 ) Table 1: Factorizations of the degenerate conic section r T Gr = 0 into two lines L ± r = 0, for various cases. The elements of g are the cofactors of G. Numerical stabilit can be achieved b swapping indices 1 when necessar to enforce G 11 < G, which also covers the cases G = 0..4 Intersections The pencil P λ) of quadrics A and A is defined as P λ) = A λa = A A 1 A λi ). Levin provides two theorems relevant to intersections of quadrics and of conics)[8]. First, if A and A intersect, then all quadrics on their pencil share the intersection. Second, in the pencil of an two intersecting quadrics there eists a ruled quadric, parametrized b a famil of lines, which can be used to parametrize the intersection curve. Good candidates for the ruled quadric are singular quadrics, which occur for λ equal to an real eigenvalue of A 1 A..4.1 Two ellipses For later reference, we document the general solution for the intersection of two ellipses A and A. The degenerate conic G = A λa, where λ is a real eigenvalue of A 1 A, can be factored as a smmetrized product of two lines ) G = L T +L + L T L + /. 10) The factorizations for various conditions are listed in Table 1. Points of intersection between a line Lr = 0 and a conic r T Ar = 0 are eigenvectors of their cross product L A) im = ɛ ijk L ja km. 11) Since the homogeneous coordinates are defined to have a value of 1 in the third component, the eigenvectors must be scaled appropriatel..4. à µ with Ãb B inspection, the quadric P = õ Ãb is singular, since it has no dependence on z. The same degenerate ). quadric is described b P = P B p Consider the translation transformation S given b S = ỹ , ỹ0 = 1 0c + s B p The translation of P is S T P S = B m Bp ) + s cs 0 cs s 0 0, ). 1) which can be solved with the quadratic formula and translated back, showing that P is the pair of intersecting planes ỹ ỹ 0 = m 0), m = 1 ) ± Be c. 13) s B p Onl +P, the plane with positive slope, ever intersects õ. In the z = 0 plane the two points of intersection are ỹ ± = ỹ 1 ± Z, ỹ 1 = 0/m, 1 = 0 + ỹ 1 ỹ 0)/m, ± = 1 ± Z/m, Z = 14) ỹ1 ỹ 1ỹ 0 0 Wm) 1. For Z = 0 the two points coincide, and +P, à b and õ are tangent, while for Z < 0 the constraints õ and à b are not consistent. Since an cross section of an ellipsoid is an ellipse, +P cuts Ãb in an ellipse where it 5

6 intersects õ. Since both Ãb and +P are smmetric in ± z, the ellipse of intersection is smmetric in ± z, so its major ais is in the plane z = 0, between points ±, ỹ ±, 0), and its minor ais is perpendicular to z = 0, between points 1, ỹ 1, ±Z). This ellipse has F coordinates parametrized b t, W p = 1 + Z m cost) ỹ 1 + Z cost) Z sint). 15) The corresponding ellipse of solutions for the neutrino momentum has the F coordinates Z/m 0 1 µ p cos t ν p = Ẽt, Ẽ = Z 0 ỹ 1, t = sin t 16) 0 Z Laborator coordinates The transformation of the laborator frame to the F frame can be accomplished b the following series of rotations: rotate the lab frame around its z-ais b µ φ, so that µ p is in the z plane, and z = z, so all polar angles are unchanged; rotate the primed frame around its -ais b ) µ θ π, so that µp coincides with ; rotate the double primed frame around its -ais so that B p is in the ỹ plane with Bỹ >= 0. The angle of the latter rotation, α, is equal to the principal value of the argument of B + B z 1). Noting that the rotation of the frame is equivalent to the opposite rotation of the vectors, the transformation from F coordinates to laborator coordinates is the rotation R = R zµ φ )R µ θ π ) R α). 17) The set of neutrino momentum solutions are given in the laborator coordinates b the parametric form ν p = Et, E = RẼ. 18) For homogeneous coordinates in the transverse plane, ν = ν ν 1 ) T, the solutions are E 11 E 1 E 13 ν = E t, E = E 1 E E 3 19) The etended representation of the solution ellipse in the transverse plane is ) 1 N = E T UE 1 0) which follows triviall given that t parametrizes the solution set of the unit circle U = diag ). Given a solution ν, the full neutrino momentum is.6 Transverse momentum constraint ν p = E E ) 1 ν. 1) The solution sets of neutrinos from decaing top quarks can be further constrained b the measured imbalance in transverse momentum of the colliding sstem ), which has a diagonalizable uncertaint matri σ. We treat events with one or two top quarks decaing to leptons..6.1 Single neutrino in final state The displacement between the measurement ) and the solution ν p = Et is 0 0 Λ = V 0 E) t, V 0 = 0 0. ) The weighted square of the displacement is χ = t T Xt, X = V 0 E) T Σ V 0 E), Σ = σ ) ) 6

7 Note that Σ weights the components of the displacement according to the measurement uncertaint. A stud of kinematic distributions of simulated neutrinos could motivate the use of a finite longitudinal uncertaint, rather than no longitudinal constraint. We seek the neutrino solution with the minimum value of χ. Differentiation of t with respect to t can be epressed as matri multiplication, t t = Dt, D = The etrema of χ occur at values of t such that χ t ) = t T Mt = 0, M = XD) T + XD. 5) Since M is smmetric, t is a point on the conic described b M. Since t is also a point on the unit circle, solutions must be points on the intersection of M and U = diag ). Refer to Section.4.1. We epect at least one minimum of χ. Since χ is cclic, no more than of at most 4 intersections can be minima, one of which is the global minimum we seek..6. Two neutrinos in final state Given two top quarks decaing to leptons, t bνµ + and t b νµ, the respective elliptical solution sets for neutrino transverse momenta are given b Equation 0), Since the measured components ν components, the are related b ν T N ν = 0, ν T N ν = 0. ) of transverse momentum imbalance are ideall just the sum of ν and ν = Γν, Γ = Note that Γ = Γ 1. We can rewrite the ν ellipse in terms of ν,. 6) ν T N ν = 0, N = Γ T N Γ 7) Solutions for ν are on both ellipses N and N simultaneousl. Refer to Section.4.1. There can be zero, two, or four intersections, discounting cases of tangenc. Each intersection implies a solution pair ν p, ν p), which can be found using Equations 1) and 6). Eamples are shown in Figure 3. If N and N do not intersect, Equation 7) is never satisfied, so the points of closest approach constitute the single solution pair most likel to result in the observed transverse momentum imbalance. The closest approach of ν to the ellipse N is often well-approimated b the etremum closest to zero of the function f = ν T N ν = t T X t, X = E T N E, 8) where we have used the parametrization 19). Etrema of f occur for f t = tt M t = 0, M = X D) T + X D, 9) which as in Section.6.1 is solved for t on intersections of the unit circle U with the conic described b M. Once ν is known, the closest approach on N is at one of its intersections with the line perpendicular to N at ν, given b ν D N ν ). Alternativel, least squares minimization of the distance between points on the ellipses is straightforward to implement using parametrization 19). 3 Summar A top quark s momentum is constrained to an ellipsoidal surface defined b its invariant mass and the observable kinematic properties of the b quark produced in its deca. Its momentum is further constrained to a slice of that surface b the W boson invariant mass and the observed kinematic properties of either of the products from the deca of the intermediate W boson. The other deca product of the W is likewise constrained to an ellipse, for which we have given a parametrization. In the event that a single neutrino is produced from top quark deca, we have calculated the momentum on its elliptical solution set which is most likel to have produced the observed transverse momentum imbalance. In the event that two neutrinos are produced from distinct top 7

8 Figure 3: Arrows indicate the neutrino black) and antineutrino gra) laborator transverse momenta, which are constrained to respective ellipses solid) b the associated products of top quark decas, in three events. The dashed ellipse is an additional constraint on the neutrino momentum from the antineutrino constraints, under the hpothesis that the measured imbalance in transverse momenta ) is equal to the sum of neutrino transverse momenta ). Eact solutions occur at intersections of the black and dashed ellipses, and can number zero left), two center), or four right). With no intersection left), the single best solution pair is for points of closest approach on the black and dashed ellipses. Solution pairs are marked in matching shaped points on the respective ellipses. Panels are marked in units of 100 GeV/c. quark decas, we have shown a method to calculate the discrete solution set of momentum pairs which eactl produce the observed transverse momentum imbalance, or the solution pair most likel to have produced it if eact solutions do not eist. Due to the geometric nature of the constraints, their description is facilitated b the use of homogeneous coordinates and etended matri representations. Algorithms for computing these solutions can be concisel implemented using common linear algebra routines. Such algorithms have been implemented b the authors in Pthon and C++ independentl, and tested for LHC and Tevatron use cases. A reference implementation is available[9]. References [1] R. H. Dalitz, G. R. Goldstein, Deca and polarization properties of the top quark, Phs. Rev. D ) [] L. Sonnenschein, Algebraic approach to solve tt dilepton equations, Phs. Rev. D 7 005) [3] L. Sonnenschein, Analtical solution of tt dilepton equations, Phs. Rev. D ) [4] J. Beringer et al. Particle Data Group), Review of particle phsics, Phs. Rev. D 86 01) [5] K. Chan, A simple mathematical approach for determining intersection of quadratic surfaces, in: W. W. Hager, S.-J. Huang, P. M. Pardalos, O. A. Prokopev Eds.), Multiscale Optimization Methods and Applications, volume 8 of Nonconve Optimization and Its Applications, Springer US, 006, pp doi: / x_13. [6] M. E. Mortenson, Mathematics for Computer Graphics Applications, Industrial Press Inc., [7] J. J. McConnell, Computer Graphics: Theor into Practice, Jones & Bartlett Learning, 006. [8] J. Levin, A parametric algorithm for drawing pictures of solid objects composed of quadric surfaces, Commun. ACM ) [9] B. A. Betchart, Reference implementation of analtic solutions for neutrino momenta in association with deca of top quarks, 013. URL: 8