Hunting the Higgs Boson Using the Cholesky Decomposition of an Indefinite Matrix

Abstract

Linear models have found widespread use in statistical investigations. For every linear model there exists a matrix representation for which the ReML (Restricted Maximum Likelihood) can be constructed from the elements of the corresponding matrix. This method works in the standard manner when the covariance structure is non-singular. It can also be used in the case where the covariance structure is singular, because the method identifies particular non-stochastic linear combinations of the observations which must be constrained to zero. In order to use this method, the Cholesky decomposition has to be generalized to symmetric and indefinite matrices using complex arithmetic methods. This method is applied to the problem of determining the spatial size (vertex) for the Higgs Boson decay in the Higgs à 4 lepton channel. A comparison based on the X² variable from the vertex fit for Higgs signal and tt background is presented and shows that the background can be greatly suppressed using the X² variable. One of the major advantages of this method over the currently adopted technique of b-tagging (Tomalin 2008) is that it is not affected by multiple interactions (pile up).