Source author record

Konstantin T. Matchev

Konstantin T. Matchev appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2026arXiv

Hunting for "Oddballs" with Machine Learning: Detecting Anomalous Exoplanets Using a Deep-Learned Low-Dimensional Representation of Transit Spectra with Autoencoders

This study explores the application of autoencoder-based machine learning techniques for anomaly detection to identify exoplanet atmospheres with unconventional chemical signatures using a low-dimensional data representation. We use the Atmospheric Big Challenge (ABC) database, a publicly available dataset with over 100,000 simulated exoplanet spectra, to construct an anomaly detection scenario by defining CO2-rich atmospheres as anomalies and CO2-poor atmospheres as the normal class. We benchmarked four different anomaly detection strategies: Autoencoder Reconstruction Loss, One-Class Support Vector Machine (1 class-SVM), K-means Clustering, and Local Outlier Factor (LOF). Each method was evaluated in both the original spectral space and the autoencoder's latent space using Receiver Operating Characteristic (ROC) curves and Area Under the Curve (AUC) metrics. To test the performance of the different methods under realistic conditions, we introduced Gaussian noise levels ranging from 10 to 50 ppm. Our results indicate that anomaly detection is consistently more effective when performed within the latent space across all noise levels. Specifically, K-means clustering in the latent space emerged as a stable and high-performing method. We demonstrate that this anomaly detection approach is robust to noise levels up to 30 ppm (consistent with realistic space-based observations) and remains viable even at 50 ppm when leveraging latent space representations. On the other hand, the performance of the anomaly detection methods applied directly in the raw spectral space degrades significantly with increasing the level of noise. This suggests that autoencoder-driven dimensionality reduction offers a robust methodology for flagging chemically anomalous targets in large-scale surveys where exhaustive retrievals are computationally prohibitive.

preprint2023arXiv

Deep Learning Symmetries and Their Lie Groups, Algebras, and Subalgebras from First Principles

We design a deep-learning algorithm for the discovery and identification of the continuous group of symmetries present in a labeled dataset. We use fully connected neural networks to model the symmetry transformations and the corresponding generators. We construct loss functions that ensure that the applied transformations are symmetries and that the corresponding set of generators forms a closed (sub)algebra. Our procedure is validated with several examples illustrating different types of conserved quantities preserved by symmetry. In the process of deriving the full set of symmetries, we analyze the complete subgroup structure of the rotation groups $SO(2)$, $SO(3)$, and $SO(4)$, and of the Lorentz group $SO(1,3)$. Other examples include squeeze mapping, piecewise discontinuous labels, and $SO(10)$, demonstrating that our method is completely general, with many possible applications in physics and data science. Our study also opens the door for using a machine learning approach in the mathematical study of Lie groups and their properties.

preprint2022arXiv

Kinematic Variables and Feature Engineering for Particle Phenomenology

Kinematic variables have been playing an important role in collider phenomenology, as they expedite discoveries of new particles by separating signal events from unwanted background events and allow for measurements of particle properties such as masses, couplings, spins, etc. For the past 10 years, an enormous number of kinematic variables have been designed and proposed, primarily for the experiments at the Large Hadron Collider, allowing for a drastic reduction of high-dimensional experimental data to lower-dimensional observables, from which one can readily extract underlying features of phase space and develop better-optimized data-analysis strategies. We review these recent developments in the area of phase space kinematics, summarizing the new kinematic variables with important phenomenological implications and physics applications. We also review recently proposed analysis methods and techniques specifically designed to leverage the new kinematic variables. As machine learning is nowadays percolating through many fields of particle physics including collider phenomenology, we discuss the interconnection and mutual complementarity of kinematic variables and machine learning techniques. We finally discuss how the utilization of kinematic variables originally developed for colliders can be extended to other high-energy physics experiments including neutrino experiments.

preprint2022arXiv

Uncertainties associated with GAN-generated datasets in high energy physics

Recently, Generative Adversarial Networks (GANs) trained on samples of traditionally simulated collider events have been proposed as a way of generating larger simulated datasets at a reduced computational cost. In this paper we point out that data generated by a GAN cannot statistically be better than the data it was trained on, and critically examine the applicability of GANs in various situations, including a) for replacing the entire Monte Carlo pipeline or parts of it, and b) to produce datasets for usage in highly sensitive analyses or sub-optimal ones. We present our arguments using information theoretic demonstrations, a toy example, as well as in the form of a formal statement, and identify some potential valid uses of GANs in collider simulations.

preprint2022arXiv

Unsupervised Machine Learning for Exploratory Data Analysis of Exoplanet Transmission Spectra

Transit spectroscopy is a powerful tool to decode the chemical composition of the atmospheres of extrasolar planets. In this paper we focus on unsupervised techniques for analyzing spectral data from transiting exoplanets. We demonstrate methods for i) cleaning and validating the data, ii) initial exploratory data analysis based on summary statistics (estimates of location and variability), iii) exploring and quantifying the existing correlations in the data, iv) pre-processing and linearly transforming the data to its principal components, v) dimensionality reduction and manifold learning, vi) clustering and anomaly detection, vii) visualization and interpretation of the data. To illustrate the proposed unsupervised methodology, we use a well-known public benchmark data set of synthetic transit spectra. We show that there is a high degree of correlation in the spectral data, which calls for appropriate low-dimensional representations. We explore a number of different techniques for such dimensionality reduction and identify several suitable options in terms of summary statistics, principal components, etc. We uncover interesting structures in the principal component basis, namely, well-defined branches corresponding to different chemical regimes of the underlying atmospheres. We demonstrate that those branches can be successfully recovered with a K-means clustering algorithm in fully unsupervised fashion. We advocate for a three-dimensional representation of the spectroscopic data in terms of the first three principal components, in order to reveal the existing structure in the data and quickly characterize the chemical class of a planet.

preprint2021arXiv

Finding Wombling Boundaries in LHC Data with Voronoi and Delaunay Tessellations

We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

preprint2021arXiv

Superfluid Effective Field Theory for Dark Matter Direct Detection

We develop an effective field theory (EFT) framework for superfluid ${}^4$He to model the interactions among quasiparticles, helium atoms and probe particles. Our effective field theory approach brings together symmetry arguments and power-counting and matches to classical fluid dynamics. We then present the decay and scattering rates for the relevant processes involving quasiparticles and helium atoms. The presented EFT framework and results can be used to understand the dynamics of thermalization in the superfluid, and can be further applied to sub-GeV dark matter direct detection with superfluid ${}^4$He.

preprint2020arXiv

InClass Nets: Independent Classifier Networks for Nonparametric Estimation of Conditional Independence Mixture Models and Unsupervised Classification

We introduce a new machine-learning-based approach, which we call the Independent Classifier networks (InClass nets) technique, for the nonparameteric estimation of conditional independence mixture models (CIMMs). We approach the estimation of a CIMM as a multi-class classification problem, since dividing the dataset into different categories naturally leads to the estimation of the mixture model. InClass nets consist of multiple independent classifier neural networks (NNs), each of which handles one of the variates of the CIMM. Fitting the CIMM to the data is performed by simultaneously training the individual NNs using suitable cost functions. The ability of NNs to approximate arbitrary functions makes our technique nonparametric. Further leveraging the power of NNs, we allow the conditionally independent variates of the model to be individually high-dimensional, which is the main advantage of our technique over existing non-machine-learning-based approaches. We derive some new results on the nonparametric identifiability of bivariate CIMMs, in the form of a necessary and a (different) sufficient condition for a bivariate CIMM to be identifiable. We provide a public implementation of InClass nets as a Python package called RainDancesVI and validate our InClass nets technique with several worked out examples. Our method also has applications in unsupervised and semi-supervised classification problems.

preprint2020arXiv

OASIS: Optimal Analysis-Specific Importance Sampling for event generation

We propose a technique called Optimal Analysis-Specific Importance Sampling (OASIS) to reduce the number of simulated events required for a high-energy experimental analysis to reach a target sensitivity. We provide recipes to obtain the optimal sampling distributions which preferentially focus the event generation on the regions of phase space with high utility to the experimental analyses. OASIS leads to a conservation of resources at all stages of the Monte Carlo pipeline, including full-detector simulation, and is complementary to approaches which seek to speed-up the simulation pipeline.

preprint2020arXiv

Singularity Variables for Missing Energy Event Kinematics

We discuss singularity variables which are properly suited for analyzing the kinematics of events with missing transverse energy at the LHC. We consider six of the simplest event topologies encountered in studies of leptonic W-bosons and top quarks, as well as in SUSY-like searches for new physics with dark matter particles. In each case, we illustrate the general prescription for finding the relevant singularity variable, which in turn helps delineate the visible parameter subspace on which the singularities are located. Our results can be used in two different ways - first, as a guide for targeting the signal-rich regions of parameter space during the stage of discovery, and second, as a sensitive focus point method for measuring the particle mass spectrum after the initial discovery.

preprint2019arXiv

Dreaming Awake: Disentangling the Underlying Physics in Case of a SUSY-like Discovery at the LHC

The purpose of this review is to investigate what kind of physics can be extracted at the LHC, assuming a discovery is made in events with missing transverse momentum, as generically expected in supersymmetry (SUSY) with R-parity conservation. To set the scene, we first discuss the collider phenomenology of the six possible electroweakino benchmark scenarios, as they provide valuable insight into what one might be facing at the LHC. We review the existing methods for mass reconstruction from measured kinematic endpoints in the distributions of suitable variables, e.g., the invariant masses of various sets of visible decay products, as well as the $M_{T2}$ and the $M_2$ types of variables. We propose to extend the application of these methods to the various topologies of fully hadronic final states, possibly with hadronically reconstructed massive bosons (W, Z or h). We test the idea with a simplified simulation of events in the main electroweakino benchmark scenarios. We find that the fully hadronic events allow the complete determination of the relevant mass spectrum. For comparison, we also review the potential of the standard kinematic endpoint methods for final states involving leptons from the decays of (on-shell or off-shell) sleptons. We find that with 300 $fb^{-1}$, the statistics for the leptonic events is very marginal and they look less promising than the fully hadronic channels. This corresponds to a complete reversal of the usual paradigm, where leptonic events comprised the gold-plated SUSY channels. Finally, we put together all available information and summarize what level of understanding of the underlying physics can be achieved. We show that, as a by-product of the mass reconstruction, it is also possible to determine the production cross sections and decay branching ratios, which in turn enable us to pinpoint the underlying model.

preprint2019arXiv

Kinematic Focus Point Method for Particle Mass Measurements in Missing Energy Events

We investigate the solvability of the event kinematics in missing energy events at hadron colliders, as a function of the particle mass ansatz. To be specific, we reconstruct the neutrino momenta in dilepton $t\bar{t}$-like events, without assuming any prior knowledge of the mass spectrum. We identify a class of events, which we call extreme events, with the property that the kinematic boundary of their allowed region in mass parameter space passes through the true mass point. We develop techniques for recognizing extreme events in the data and demonstrate that they are abundant in a realistic data sample, due to expected singularities in phase space. We propose a new method for mass measurement whereby we obtain the true values of the mass parameters as the focus point of the kinematic boundaries for all events in the data sample. Since the masses are determined from a relatively sharp peak structure (the density of kinematic boundary curves), the method avoids some of the systematic errors associated with other techniques. We show that this new approach is complementary to previously considered methods in the literature where one studies the solvability of the kinematic constraints throughout the mass parameter space. In particular, we identify a problematic direction in mass space of nearly 100% solvability, and then show that the focus point method is effective in lifting the degeneracy.

preprint2016arXiv

OPTIMASS: A Package for the Minimization of Kinematic Mass Functions with Constraints

Reconstructed mass variables, such as $M_2$, $M_{2C}$, $M_T^\star$, and $M_{T2}^W$, play an essential role in searches for new physics at hadron colliders. The calculation of these variables generally involves constrained minimization in a large parameter space, which is numerically challenging. We provide a C++ code, OPTIMASS, which interfaces with the MINUIT library to perform this constrained minimization using the Augmented Lagrangian Method. The code can be applied to arbitrarily general event topologies and thus allows the user to significantly extend the existing set of kinematic variables. We describe this code and its physics motivation, and demonstrate its use in the analysis of the fully leptonic decay of pair-produced top quarks using the $M_2$ variables.

preprint2016arXiv

The 750 GeV Diphoton Excess May Not Imply a 750 GeV Resonance

We discuss non-standard interpretations of the 750 GeV diphoton excess recently reported by the ATLAS and CMS Collaborations which do not involve a new, relatively broad, resonance with a mass near 750 GeV. Instead, we consider the sequential cascade decay of a much heavier, possibly quite narrow, resonance into two photons along with one or more invisible particles. The resulting diphoton invariant mass signal is generically rather broad, as suggested by the data. We examine three specific event topologies - the antler, the sandwich, and the 2-step cascade decay, and show that they all can provide a good fit to the observed published data. In each case, we delineate the preferred mass parameter space selected by the best fit. In spite of the presence of invisible particles in the final state, the measured missing transverse energy is moderate, due to its anti- correlation with the diphoton invariant mass. We comment on the future prospects of discriminating with higher statistics between our scenarios, as well as from more conventional interpretations.

preprint2015arXiv

Discovering New Physics with Voronoi Tessellations

High energy experimental data can be viewed as a sampling of the relevant phase space. We point out that one can apply Voronoi tessellations in order to understand the underlying probability distributions in this phase space. Interesting features in the data can then be discovered by studying the properties of the ensemble of Voronoi cells. For illustration, we demonstrate the detection of kinematic "edges" in two dimensions, which may signal physics beyond the standard model. We motivate the algorithm with some analytical results derived for perfect lattices, and show that the method is further improved with the addition of a few Voronoi relaxation steps via Lloyd's method.

preprint2015arXiv

Edge Detecting New Physics the Voronoi Way

We point out that interesting features in high energy physics data can be determined from properties of Voronoi tessellations of the relevant phase space. For illustration, we focus on the detection of kinematic "edges" in two dimensions, which may signal physics beyond the standard model. After deriving some useful geometric results for Voronoi tessellations on perfect grids, we propose several algorithms for tagging the Voronoi cells in the vicinity of kinematic edges in real data. We show that the efficiency is improved by the addition of a few Voronoi relaxation steps via Lloyd's method. By preserving the maximum spatial resolution of the data, Voronoi methods can be a valuable addition to the data analysis toolkit at the LHC.

preprint2015arXiv

Improving the sensitivity of stop searches with on-shell constrained invariant mass variables

The search for light stops is of paramount importance, both in general as a promising path to the discovery of beyond the standard model physics and more specifically as a way of evaluating the success of the naturalness paradigm. While the LHC experiments have ruled out much of the relevant parameter space, there are "stop gaps", i.e., values of sparticle masses for which existing LHC analyses have relatively little sensitivity to light stops. We point out that techniques involving on-shell constrained M_2 variables can do much to enhance sensitivity in this region and hence help close the stop gaps. We demonstrate the use of these variables for several benchmark points and describe the effect of realistic complications, such as detector effects and combinatorial backgrounds, in order to provide a useful toolkit for light stop searches in particular, and new physics searches at the LHC in general.

preprint2015arXiv

On-shell constrained $M_2$ variables with applications to mass measurements and topology disambiguation

We consider a class of on-shell constrained mass variables that are 3+1 dimensional generalizations of the Cambridge $M_{T2}$ variable and that automatically incorporate various assumptions about the underlying event topology. The presence of additional on-shell constraints causes their kinematic distributions to exhibit sharper endpoints than the usual $M_{T2}$ distribution. We study the mathematical properties of these new variables, e.g., the uniqueness of the solution selected by the minimization over the invisible particle 4-momenta. We then use this solution to reconstruct the masses of various particles along the decay chain. We propose several tests for validating the assumed event topology in missing energy events from new physics. The tests are able to determine: 1) whether the decays in the event are two-body or three-body, 2) if the decay is two-body, whether the intermediate resonances in the two decay chains are the same, and 3) the exact sequence in which the visible particles are emitted from each decay chain.

preprint2015arXiv

The Hierarchy Solution to the LHC Inverse Problem

Supersymmetric (SUSY) models, even those described by relatively few parameters, generically allow many possible SUSY particle (sparticle) mass hierarchies. As the sparticle mass hierarchy determines, to a great extent, the collider phenomenology of a model, the enumeration of these hierarchies is of the utmost importance. We therefore provide a readily generalizable procedure for determining the number of sparticle mass hierarchies in a given SUSY model. As an application, we analyze the gravity-mediated SUSY breaking scenario with various combinations of GUT-scale boundary conditions involving different levels of universality among the gaugino and scalar masses. For each of the eight considered models, we provide the complete list of forbidden hierarchies in a compact form. Our main result is that the complete (typically rather large) set of forbidden hierarchies among the eight sparticles considered in this analysis can be fully specified by just a few forbidden relations involving much smaller subsets of sparticles.

preprint2015arXiv

Using sorted invariant mass variables to evade combinatorial ambiguities in cascade decays

The classic method for mass determination in a SUSY-like cascade decay chain relies on measurements of the kinematic endpoints in the invariant mass distributions of suitable collections of visible decay products. However, the procedure is complicated by combinatorial ambiguities: e.g., the visible final state particles may be indistinguishable (as in the case of QCD jets), or one may not know the exact order in which they are emitted along the decay chain. In order to avoid such combinatorial ambiguities, we propose to treat the final state particles fully democratically and consider the sorted set of the invariant masses of all possible partitions of the visible particles in the decay chain. In particular, for a decay to N visible particles, one considers the sorted sets of all possible n-body invariant mass combinations (2 <= n <= N) and determines the kinematic endpoint m_(n,r)^max of the distribution of the r-th largest n-body invariant mass m_(n,r) for each possible value of n and r. For the classic example of a squark decay in supersymmetry, we provide analytical formulas for the interpretation of these endpoints in terms of the underlying physical masses. We point out that these measurements can be used to determine the structure of the decay topology, e.g., the number and position of intermediate on-shell resonances.

preprint2014arXiv

Beyond Geolocating: Constraining Higher Dimensional Operators in $H \to 4\ell$ with Off-Shell Production and More

We extend the study of Higgs boson couplings in the "golden" $gg\to H \to ZZ^\ast \to 4\ell$ channel in two important respects. First, we demonstrate the importance of off-shell Higgs boson production ($gg\to H^\ast \to ZZ \to 4\ell$) in determining which operators contribute to the $HZZ$ vertex. Second, we include the five operators of lowest non-trivial dimension, including the $Z_μZ^μ\Box H$ and $H Z_μ\Box Z^μ$ operators that are often neglected. We point out that the former operator can be severely constrained by the measurement of the off-shell $H^\ast \to ZZ$ rate and/or unitarity considerations. We provide analytic expressions for the off-peak cross-sections in the presence of these five operators. On-shell, the $Z_μZ^μ\Box H$ operator is indistinguishable from its Standard Model counterpart $H Z_μZ^μ$, while the $H Z_μ\Box Z^μ$ operator can be probed, in particular, by the $Z^\ast$ invariant mass distribution.

preprint2014arXiv

Discoveries far from the Lamppost with Matrix Elements and Ranking

The prevalence of null results in searches for new physics at the LHC motivates the effort to make these searches as model-independent as possible. We describe procedures for adapting the Matrix Element Method for situations where the signal hypothesis is not known a priori. We also present general and intuitive approaches for performing analyses and presenting results, which involve the flattening of background distributions using likelihood information. The first flattening method involves ranking events by background matrix element, the second involves quantile binning with respect to likelihood (and other) variables, and the third method involves reweighting histograms by the inverse of the background distribution.

preprint2014arXiv

Exploring Theory Space with Monte Carlo Reweighting

Theories of new physics often involve a large number of unknown parameters which need to be scanned. Additionally, a putative signal in a particular channel may be due to a variety of distinct models of new physics. This makes experimental attempts to constrain the parameter space of motivated new physics models with a high degree of generality quite challenging. We describe how the reweighting of events may allow this challenge to be met, as fully simulated Monte Carlo samples generated for arbitrary benchmark models can be effectively re-used. In particular, we suggest procedures that allow more efficient collaboration between theorists and experimentalists in exploring large theory parameter spaces in a rigorous way at the LHC.

preprint2014arXiv

Model-Independent Searches Using Matrix Element Ranking

Thus far the LHC experiments have yet to discover beyond-the-standard-model physics. This motivates efforts to search for new physics in model independent ways. In this spirit, we describe procedures for using a variant of the Matrix Element Method to search for new physics without regard to a specific signal hypothesis. To make the resulting variables more intuitive, we also describe how these variables can be "flattened", which makes the resulting distributions more visually meaningful.

preprint2014arXiv

Probing resonance decays to two visible and multiple invisible particles

We consider the decay of a generic resonance to two visible particles and any number of invisible particles. We show that the shape of the invariant mass distribution of the two visible particles is sensitive to both the mass spectrum of the new particles, as well as the decay topology. We provide the analytical formulas describing the invariant mass shapes for the nine simplest topologies (with up to two invisible particles in the final state). Any such distribution can be simply categorized by its endpoint, peak location and curvature, which are typically sufficient to discriminate among the competing topologies. In each case, we list the effective mass parameters which can be measured by experiment. In certain cases, the invariant mass shape is sufficient to completely determine the new particle mass spectrum, including the overall mass scale.

preprint2013arXiv

Geolocating the Higgs Boson Candidate at the LHC

The latest results from the ATLAS and CMS experiments at the CERN Large Hadron Collider (LHC) unequivocally confirm the existence of a resonance, $X$, with mass near 125 GeV which could be the Higgs boson of the Standard Model. Measuring the properties (quantum numbers and couplings) of this resonance is of paramount importance. Initial analyses by the LHC collaborations disfavor specific alternative benchmark hypotheses, e.g. pure pseudoscalars or gravitons. However, this is just the first step in a long-term program of detailed measurements. We consider the most general set of operators in the decay channels $X \to ZZ$, $WW$, $Zγ$, $γγ$ and derive the constraint implied by the measured rate. This allows us to provide a useful parametrization of the orthogonal independent Higgs coupling degrees of freedom as coordinates on a suitably defined sphere.

preprint2013arXiv

Kaluza-Klein Dark Matter: Direct Detection vis-a-vis LHC (2013 update)

We present updated results on the complementarity between high-energy colliders and dark matter direct detection experiments in the context of Universal Extra Dimensions (UED). In models with relatively small mass splittings between the dark matter candidate and the rest of the (colored) spectrum, the collider sensitivity is diminished, but direct detection rates are enhanced. UED provide a natural framework to study such mass degeneracies. We discuss the detection prospects for the KK photon $γ_1$ and the KK $Z$-boson $Z_1$, combining the expected LHC reach with cosmological constraints from WMAP/Planck, and the sensitivity of current or planned direct detection experiments. Allowing for general mass splittings, neither colliders, nor direct detection experiments by themselves can explore all of the relevant KK dark matter parameter space. Nevertheless, they probe different parameter space regions, and the combination of the two types of constraints can be quite powerful.

preprint2013arXiv

The Matrix Element Method: Past, Present, and Future

The increasing use of multivariate methods, and in particular the Matrix Element Method (MEM), represents a revolution in experimental particle physics. With continued exponential growth in computing capabilities, the use of sophisticated multivariate methods-- already common-- will soon become ubiquitous and ultimately almost compulsory. While the existence of sophisticated algorithms for disentangling signal and background might naively suggest a diminished role for theorists, the use of the MEM, with its inherent connection to the calculation of differential cross sections will benefit from collaboration between theorists and experimentalists. In this white paper, we will briefly describe the MEM and some of its recent uses, note some current issues and potential resolutions, and speculate about exciting future opportunities.

preprint2013arXiv

The role of interference in unraveling the ZZ-couplings of the newly discovered boson at the LHC

We present a general procedure for measuring the tensor structure of the coupling of the scalar Higgs-like boson recently discovered at the LHC to two Z bosons, including the effects of interference among different operators. To motivate our concern with this interference, we explore the parameter space of the couplings in the effective theory describing these interactions and illustrate the effects of interference on the differential dilepton mass distributions. Kinematic discriminants for performing coupling measurements that utilize the effects of interference are developed and described. We present projections for the sensitivity of coupling measurements that use these discriminants in future LHC operation in a variety of physics scenarios.

preprint2012arXiv

Focus Point Supersymmetry Redux

Recent results from Higgs boson and supersymmetry searches at the Large Hadron Collider provide strong new motivations for supersymmetric theories with heavy superpartners. We reconsider focus point supersymmetry (FP SUSY), in which all squarks and sleptons may have multi-TeV masses without introducing fine-tuning in the weak scale with respect to variations in the fundamental SUSY-breaking parameters. We examine both FP SUSY and its familiar special case, the FP region of mSUGRA/CMSSM, and show that they are beautifully consistent with all particle, astroparticle, and cosmological data, including Higgs boson mass limits, null results from SUSY searches, electric dipole moments, b -> s gamma, B_s -> mu^+ mu^-, the thermal relic density of neutralinos, and dark matter searches. The observed deviation of the muon's anomalous magnetic moment from its standard model value may also be explained in FP SUSY, although not in the FP region of mSUGRA/CMSSM. In light of recent data, we advocate refined searches for FP SUSY and related scenarios with heavy squarks and sleptons, and we present a simplified parameter space to aid such analyses.

preprint2012arXiv

Precision studies of the Higgs boson decay channel H -> ZZ -> 4l with MEKD

The importance of the H -> ZZ -> 4l "golden" channel was shown by its major role in the discovery, by the ATLAS and CMS collaborations, of a Higgs-like boson with mass near 125 GeV. We analyze the discrimination power of the matrix element method both for separating the signal from the irreducible ZZ background and for distinguishing various spin and parity hypotheses describing a signal in this channel. We show that the proper treatment of interference effects associated with permutations of identical leptons in the four electron and four muon final states plays an important role in achieving the best sensitivity in measuring the properties of the newly discovered boson. We provide a code, MEKD, that calculates kinematic discriminants based on the full leading order matrix elements and which will aid experimentalists and phenomenologists in their continuing studies of the H -> ZZ -> 4l channel.

preprint2012arXiv

Spin effects in the antler event topology at hadron colliders

We investigate spin correlation effects in the "antler" event topology pp-> A-> B1, B2 -> l^{-}, C1, l^{+}, C2 at the LHC. We study the shapes of several kinematic variables, including the relative pseudorapidity, relative azimuthal angle and the energies of the two leptons, as well as several mass variables M_{ll}, Meff, \sqrt{s}_{min}, MT2, MCT and MCTx. We focus on the two kinematic extremes of \sqrt{s} - threshold and infinity - and derive analytical expressions for the differential distributions of several variables, most notably the cos{θ_{ll}}^* variable proposed by Barr in hep-ph/0511115. For all possible spin assignments of particles A, B and C, we derive the cos{θ_{ll}}^* differential distribution at threshold, including the effects of spin correlations. Our analytical results help identify the problematic cases for spin discrimination.

preprint2010arXiv

A general method for determining the masses of semi-invisibly decaying particles at hadron colliders

We present a general solution to the long standing problem of determining the masses of pair-produced, semi-invisibly decaying particles at hadron colliders. We define two new transverse kinematic variables, $M_{CT_\perp}$ and $M_{CT_\parallel}$, which are suitable one-dimensional projections of the contransverse mass $M_{CT}$. We derive analytical formulas for the boundaries of the kinematically allowed regions in the $(M_{CT_\perp},M_{CT_\parallel})$ and $(M_{CT_\perp},M_{CT})$ parameter planes, and introduce suitable variables $D_{CT_\parallel}$ and $D_{CT}$ to measure the distance to those boundaries on an event per event basis. We show that the masses can be reliably extracted from the endpoint measurements of $M_{CT_\perp}^{max}$ and $D_{CT}^{min}$ (or $D_{CT_\parallel}^{min}$). We illustrate our method with dilepton $t\bar{t}$ events at the LHC.

preprint2010arXiv

Dark Matter Particle Spectroscopy at the LHC: Generalizing MT2 to Asymmetric Event Topologies

We consider SUSY-like missing energy events at hadron colliders and critically examine the common assumption that the missing energy is the result of two identical missing particles. In order to experimentally test this hypothesis, we generalize the subsystem MT2 variable to the case of asymmetric event topologies, where the two SUSY decay chains terminate in different "children" particles. In this more general approach, the endpoint MT2max of the MT2 distribution now gives the mass Mp(Mc(a),Mc(b)) of the parent particle as a function of two input children masses Mc(a) and Mc(b). We propose two methods for an independent determination of the individual children masses Mc(a) and Mc(b). First, in the presence of upstream transverse momentum P(UTM) the corresponding function Mp(Mc(a),Mc(b),P(UTM)) is independent of P(UTM) at precisely the right values of the children masses. Second, the previously discussed MT2 "kink" is now generalized to a "ridge" on the 2-dimensional surface Mp(Mc(a),Mc(b)). As we show in several examples, quite often there is a special point along that ridge which marks the true values of the children masses. Our results allow collider experiments to probe a multi-component dark matter sector directly and without any theoretical prejudice.

preprint2010arXiv

How to look for supersymmetry under the lamppost at the LHC

We apply a model-independent, agnostic approach to the collider phenomenology of supersymmetry (SUSY), in which all mass parameters are taken as free inputs at the weak scale. We consider the gauginos, higgsinos, and the first two generations of sleptons and squarks, and analyze all possible mass hierarchies among them ($4\times 8!=161,280$ in total) in which the lightest superpartner is neutral, leading to missing energy. In each case, we identify the full set of the dominant (i.e. least suppressed by phase space, small mixing angles or Yukawa couplings) decay chains originating from the lightest colored superpartner. Our exhaustive search reveals several quite dramatic yet unexplored multilepton signatures with up to 8 isolated leptons (plus possibly up to 2 massive gauge or Higgs bosons) in the final state. Such events are spectacular, background-free for all practical purposes, and may lead to a discovery of SUSY in the very early stage ($\sim 10\ {\rm pb}^{-1}$) of LHC operations at 7 TeV.

preprint2010arXiv

Minimal Universal Extra Dimensions in CalcHEP/CompHEP

We present an implementation of the model of minimal universal extra dimensions (MUED) in CalcHEP/CompHEP. We include all level-1 and level-2 Kaluza-Klein (KK) particles outside the Higgs sector. The mass spectrum is automatically calculated at one loop in terms of the two input parameters in MUED: the radius of the extra dimension and the cut-off scale of the model. We implement both the KK number conserving and the KK number violating interactions of the KK particles. We also account for the proper running of the gauge coupling constants above the electroweak scale. The implementation has been extensively cross-checked against known analytical results in the literature and numerical results from other programs. Our files are publicly available and can be used to perform various automated calculations within the MUED model.

preprint2010arXiv

Precision sparticle spectroscopy in the inclusive same-sign dilepton channel at LHC

The inclusive same-sign dilepton channel is already recognized as a promising discovery signature for supersymmetry in the early days of the LHC. We point out that it can also be used for precision measurements of sparticle masses after the initial discovery. As an illustration, we consider the LM6 CMS study point in minimal supergravity, where the same-sign leptons most often result from chargino decays to sneutrinos. We discuss three different techniques for determining the chargino and sneutrino masses in an inclusive manner, i.e. using only the two well measured lepton momenta, while treating all other upstream objects in the event as a single entity of total transverse momentum PT. This approach takes full advantage of the large production rates of colored superpartners, but does not rely on the poorly measured hadronic jets, and avoids any jet combinatorics problems. We discuss the anticipated precision of our methods in the early LHC data.

preprint2010arXiv

RECO level \sqrt{s}_{min} and subsystem \sqrt{s}_{min}: improved global inclusive variables for measuring the new physics mass scale in missing energy events at hadron colliders

The variable \sqrt{s}_{min} was originally proposed in arXiv:0812.1042 as a model-independent, global and fully inclusive measure of the new physics mass scale in missing energy events at hadron colliders. In the original incarnation of \sqrt{s}_{min}, however, the connection to the new physics mass scale was blurred by the effects of the underlying event, most notably initial state radiation and multiple parton interactions. In this paper we advertize two improved variants of the \sqrt{s}_{min} variable, which overcome this problem. First we show that by evaluating the \sqrt{s}_{min} variable at the RECO level, in terms of the reconstructed objects in the event, the effects from the underlying event are significantly diminished and the nice correlation between the peak in the \sqrt{s}_{min}^{(reco)} distribution and the new physics mass scale is restored. Secondly, the underlying event problem can be avoided altogether when the \sqrt{s}_{min} concept is applied to a subsystem of the event which does not involve any QCD jets. We supply an analytic formula for the resulting subsystem \sqrt{s}_{min}^{(sub)} variable and show that its peak exhibits the usual correlation with the mass scale of the particles produced in the subsystem. Finally, we contrast \sqrt{s}_{min} to other popular inclusive variables such as H_T, M_{Tgen} and M_{TTgen}. We illustrate our discussion with several examples from supersymmetry, and with dilepton events from top quark pair production.

preprint2010arXiv

Superpartner mass measurements with 1D decomposed MT2

We propose a new model-independent technique for mass measurements in missing energy events at hadron colliders. We illustrate our method with the most challenging case of a short, single-step decay chain. We consider inclusive same-sign chargino pair production in supersymmetry, followed by leptonic decays to sneutrinos. We introduce two one-dimensional decompositions of the Cambridge MT2 variable: MT2_\parallel and MT2_\perp, on the direction of the upstream transverse momentum PT and the direction orthogonal to it, respectively. We show that the sneutrino mass can be measured directly by minimizing the number of events N in which MT2 exceeds a certain threshold, conveniently measured from the endpoint MT2^max_\perp.

preprint2009arXiv

Shedding Light on the Dark Sector with Direct WIMP Production

A Weakly Interacting Massive Particle (WIMP) provides an attractive dark matter candidate, and should be within reach of the next generation of high-energy colliders. We consider the process of direct WIMP pair-production, accompanied by an initial-state radiation photon, in electron-positron collisions at the proposed International Linear Collider (ILC). We present a parametrization of the differential cross section for this process which conveniently separates the model-independent information provided by cosmology from the model-dependent inputs from particle physics. As an application, we consider two simple models, one supersymmetric, and another of the "universal extra dimensions" (UED) type. The discovery reach of the ILC and the expected precision of parameter measurements are studied in each model. In addition, for each of the two examples, we also investigate the ability of the ILC to distinguish between the two models through a shape-discrimination analysis of the photon energy spectrum. We show that with sufficient beam polarization the alternative model interpretation can be ruled out in a large part of the relevant parameter space.

preprint2008arXiv

Kaluza-Klein Dark Matter: Direct Detection vis-a-vis LHC (IDM 2008 proceeding)

We explore the phenomenology of Kaluza-Klein (KK) dark matter in very general models with universal extra dimensions (UEDs), emphasizing the complementarity between high-energy colliders and dark matter direct detection experiments. In models with relatively small mass splittings between the dark matter candidate and the rest of the (colored) spectrum, the collider sensitivity is diminished, but direct detection rates are enhanced. UEDs provide a natural framework for such mass degeneracies. We consider both 5-dimensional and 6-dimensional non-minimal UED models, and discuss the detection prospects for various KK dark matter candidates: the KK photon γ_1 (5D) the KK Z-boson Z_1 (5D) and the spinless KK photon γ_H (6D). We combine collider limits such as electroweak precision data and expected LHC reach, with cosmological constraints from WMAP and the sensitivity of current or planned direct detection experiments. Allowing for general mass splittings, we show that neither colliders, nor direct detection experiments by themselves can explore all of the relevant KK dark matter parameter space. Nevertheless, they probe different parameter space regions and the combination of the two types of constraints can be quite powerful. For example, in the case of γ_1 in 5D UEDs the relevant parameter space will be almost completely covered by the combined LHC and direct detection sensitivities expected in the near future.

preprint2005arXiv

Discrimination of Supersymmetry and Universal Extra Dimensions at Hadron Colliders

We contrast the experimental signatures of low energy supersymmetry and the model of Universal Extra Dimensions and discuss various methods for their discrimination at hadron colliders. We study the discovery reach of the Tevatron and the LHC for level 2 Kaluza-Klein modes, which would indicate the presence of extra dimensions. We find that with 100 ${\rm fb}^{-1}$ of data the LHC will be able to discover the $γ_2$ and $Z_2$ KK modes as separate resonances if their masses are below 2 TeV. We also investigate the possibility to differentiate the spins of the superpartners and KK modes by means of the asymmetry method of Barr.

preprint2005arXiv

Robust Gamma Ray Signature of WIMP Dark Matter

If dark matter consists of weakly interacting massive particles (WIMPs), annihilation of WIMPs in the galactic center may lead to an observable enhancement of high energy gamma ray fluxes. We predict the shape and normalization of the component of the flux due to final state radiation by charged particles produced in WIMP annihilation events. The prediction is made without any assumptions about the microscopic theory responsible for WIMPs, and depends only mildly on the unknown distribution of the total WIMP annihilation cross section among the possible final states. In particular, if the WIMPs annihilate into a pair of charged fermions (leptons or quarks), the photon spectrum possesses a sharp edge feature, dropping abruptly at a photon energy equal to the WIMP mass. If such a feature is observed, it would provide strong evidence for the WIMP-related nature of the flux enhancement, as well as a measurement of the WIMP mass. We discuss the prospects for observational discovery of this feature at ground-based and space-based gamma ray telescopes.