Source author record

Matthew D. Schwartz

Matthew D. Schwartz 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

49works
9topics
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

49 published item(s)

preprint2026arXiv

Resummation of the C-Parameter Sudakov Shoulder Using Effective Field Theory

The C-parameter distribution in $e^+e^-$ annihilation exhibits a kinematic shoulder at $C = 3/4$, where three-parton final states reach their maximum and a fourth parton is required to exceed it. This boundary generates large logarithms that must be resummed. Using soft-collinear effective theory, we derive a factorization theorem involving new jet and soft functions specific to the C-parameter measurement, in which soft radiation contributes quadratically in transverse momentum. This quadratic structure explains the step discontinuity at leading order. We compute all ingredients at one loop, validate against Monte Carlo, and present matched NLL+NLO results. Unlike thrust and heavy jet mass, the C-parameter has no Sudakov--Landau pole, making momentum-space resummation straightforward. All calculations, numerical analysis, and manuscript preparation were performed by Claude, an AI assistant developed by Anthropic, working under physicist supervision.

preprint2022arXiv

Simplifying Polylogarithms with Machine Learning

Polylogrithmic functions, such as the logarithm or dilogarithm, satisfy a number of algebraic identities. For the logarithm, all the identities follow from the product rule. For the dilogarithm and higher-weight classical polylogarithms, the identities can involve five functions or more. In many calculations relevant to particle physics, complicated combinations of polylogarithms often arise from Feynman integrals. Although the initial expressions resulting from the integration usually simplify, it is often difficult to know which identities to apply and in what order. To address this bottleneck, we explore to what extent machine learning methods can help. We consider both a reinforcement learning approach, where the identities are analogous to moves in a game, and a transformer network approach, where the problem is viewed analogously to a language-translation task. While both methods are effective, the transformer network appears more powerful and holds promise for practical use in symbolic manipulation tasks in mathematical physics.

preprint2021arXiv

Challenges for Unsupervised Anomaly Detection in Particle Physics

Anomaly detection relies on designing a score to determine whether a particular event is uncharacteristic of a given background distribution. One way to define a score is to use autoencoders, which rely on the ability to reconstruct certain types of data (background) but not others (signals). In this paper, we study some challenges associated with variational autoencoders, such as the dependence on hyperparameters and the metric used, in the context of anomalous signal (top and $W$) jets in a QCD background. We find that the hyperparameter choices strongly affect the network performance and that the optimal parameters for one signal are non-optimal for another. In exploring the networks, we uncover a connection between the latent space of a variational autoencoder trained using mean-squared-error and the optimal transport distances within the dataset. We then show that optimal transport distances to representative events in the background dataset can be used directly for anomaly detection, with performance comparable to the autoencoders. Whether using autoencoders or optimal transport distances for anomaly detection, we find that the choices that best represent the background are not necessarily best for signal identification. These challenges with unsupervised anomaly detection bolster the case for additional exploration of semi-supervised or alternative approaches.

preprint2021arXiv

Modern Machine Learning and Particle Physics

Over the past five years, modern machine learning has been quietly revolutionizing particle physics. Old methodology is being outdated and entirely new ways of thinking about data are becoming commonplace. This article will review some aspects of the natural synergy between modern machine learning and particle physics, focusing on applications at the Large Hadron Collider. A sampling of examples is given, from signal/background discrimination tasks using supervised learning to direct data-driven approaches. Some comments on persistent challenges and possible future directions for the field are included at the end.

preprint2021arXiv

Quark-Gluon Backscattering in the Regge Limit at One-Loop

At small momentum transfer, the quark-gluon scattering cross section dσ/dt has a power-law divergence in the backward scattering region where the outgoing quark is nearly collinear to the incoming gluon. In this Regge limit |t|<< s, the leading behavior of the 2 2 amplitude can be described by the exchange of Glauber quarks. In Soft-Collinear Effective Theory (SCET) at leading power, Glauber quark exchange is given by five non-local Glauber quark operators, of which only one is generated at tree-level. We show that at leading power the QCD amplitude for quark-gluon backscattering at one-loop can be exactly reproduced by SCET using the tree-level Glauber operator. The agreement between QCD and SCET of the ultraviolet, infrared, and rapidity divergences as well as all logarithms, Glauber phases and finite parts for all polarizations of the external gluons is a strong check on the effective theory. We find that the entire one-loop matching vanishes -- there is no correction to the operator generated at tree-level, and the coefficients of the other four operators remain zero at one-loop. This suggests that SCET with Glauber operators may be useful for uncovering new aspects of Regge physics in a systematically improveable way.

preprint2020arXiv

ABCDisCo: Automating the ABCD Method with Machine Learning

The ABCD method is one of the most widely used data-driven background estimation techniques in high energy physics. Cuts on two statistically-independent classifiers separate signal and background into four regions, so that background in the signal region can be estimated simply using the other three control regions. Typically, the independent classifiers are chosen "by hand" to be intuitive and physically motivated variables. Here, we explore the possibility of automating the design of one or both of these classifiers using machine learning. We show how to use state-of-the-art decorrelation methods to construct powerful yet independent discriminators. Along the way, we uncover a previously unappreciated aspect of the ABCD method: its accuracy hinges on having low signal contamination in control regions not just overall, but relative to the signal fraction in the signal region. We demonstrate the method with three examples: a simple model consisting of three-dimensional Gaussians; boosted hadronic top jet tagging; and a recasted search for paired dijet resonances. In all cases, automating the ABCD method with machine learning significantly improves performance in terms of ABCD closure, background rejection and signal contamination.

preprint2020arXiv

An $S$-Matrix for Massless Particles

The traditional $S$-matrix does not exist for theories with massless particles, such as quantum electrodynamics. The difficulty in isolating asymptotic states manifests itself as infrared divergences at each order in perturbation theory. Building on insights from the literature on coherent states and factorization, we construct an $S$-matrix that is free of singularities order-by-order in perturbation theory. Factorization guarantees that the asymptotic evolution in gauge theories is universal, i.e. independent of the hard process. Although the hard $S$-matrix element is computed between well-defined few particle Fock states, dressed/coherent states can be seen to form as intermediate states in the calculation of hard $S$-matrix elements. We present a framework for the perturbative calculation of hard $S$-matrix elements combining Lorentz-covariant Feynman rules for the dressed-state scattering with time-ordered perturbation theory for the asymptotic evolution. With hard cutoffs on the asymptotic Hamiltonian, the cancellation of divergences can be seen explicitly. In dimensional regularization, where the hard cutoffs are replaced by a renormalization scale, the contribution from the asymptotic evolution produces scaleless integrals that vanish. A number of illustrative examples are given in QED, QCD, and $\mathcal{N}=4$ super Yang-Mills theory.

preprint2020arXiv

Sequential Discontinuities of Feynman Integrals and the Monodromy Group

We generalize the relation between discontinuities of scattering amplitudes and cut diagrams to cover sequential discontinuities (discontinuities of discontinuities) in arbitrary momentum channels. The new relations are derived using time-ordered perturbation theory, and hold at phase-space points where all cut momentum channels are simultaneously accessible. As part of this analysis, we explain how to compute sequential discontinuities as monodromies and explore the use of the monodromy group in characterizing the analytic properties of Feynman integrals. We carry out a number of cross-checks of our new formulas in polylogarithmic examples, in some cases to all loop orders.

preprint2016arXiv

A direct approach to quantum tunneling

The decay rates of quasistable states in quantum field theories are usually calculated using instanton methods. Standard derivations of these methods rely in a crucial way upon deformations and analytic continuations of the physical potential, and on the saddle point approximation. While the resulting procedure can be checked against other semi-classical approaches in some one-dimensional cases, it is challenging to trace the role of the relevant physical scales, and any intuitive handle on the precision of the approximations involved are at best obscure. In this paper, we use a physical definition of the tunneling probability to derive a formula for the decay rate in both quantum mechanics and quantum field theory directly from the Minkowski path integral, without reference to unphysical deformations of the potential. There are numerous benefits to this approach, from non-perturbative applications to precision calculations and aesthetic simplicity.

preprint2016arXiv

Factorization for groomed jet substructure beyond the next-to-leading logarithm

Jet grooming algorithms are widely used in experimental analyses at hadron colliders to remove contaminating radiation from within jets. While the algorithms perform a great service to the experiments, their intricate algorithmic structure and multiple parameters has frustrated precision theoretic understanding. In this paper, we demonstrate that one particular groomer called soft drop actually makes precision jet substructure easier. In particular, we derive a factorization formula for a large class of soft drop jet substructure observables, including jet mass. The essential observation that allows for this factorization is that, without the soft wide-angle radiation groomed by soft drop, all singular contributions are collinear. The simplicity and universality of the collinear limit in QCD allows us to show that to all orders, the normalized differential cross section has no contributions from non-global logarithms. It is also independent of process, up to the relative fraction of quark and gluon jets. In fact, soft drop allows us to define this fraction precisely. The factorization theorem also explains why soft drop observables are less sensitive to hadronization than their ungroomed counterparts. Using the factorization theorem, we resum the soft drop jet mass to next-to-next-to-leading logarithmic accuracy. This requires calculating some clustering effects that are closely related to corresponding effects found in jet veto calculations. We match our resummed calculation to fixed order results for both $e^+e^-\to$ dijets and $pp\to Z+j$ events, producing the first jet substructure predictions (groomed or ungroomed) to this accuracy for the LHC.

preprint2016arXiv

Precision direct photon spectra at high energy and comparison to the 8 TeV ATLAS data

The direct photon spectrum is computed to the highest currently available precision and compared to ATLAS data from 8 TeV collisions at the LHC. The prediction includes threshold resummation at next-to-next-to-next-to-leading logarithmic order through the program PeTeR, matched to next-to-leading fixed order with fragmentation effects using JetPhox and includes the resummation of leading-logarithmic electroweak Sudakov effects. Remarkably, improved agreement with data can be seen when each component of the calculation is added successively. This comparison demonstrates the importance of both threshold logs and electroweak Sudakov effects. Numerical values for the predictions are included.

preprint2016arXiv

Precision physics with pile-up insensitive observables

To deepen the search for beyond the Standard Model physics, the Large Hadron Collider is pushing to higher and higher luminosity. At high luminosity, precision physics becomes increasingly difficult due to contamination from additional proton collisions per bunch crossing called pile-up. In recent years, many methods have been developed to cull this excess mostly low-energy radiation away from important signal regions, but it has been unclear if these methods were amenable to systematically-improvable theoretical understanding. In this paper, it is shown that one such method, soft drop jet grooming, has excellent theoretical properties: it is ultra-local, depending on only radiation within a jet, and it is free of non-global logarithms. Calculations of the soft drop jet mass and related observables are presented at next-to-next-to-leading logarithmic accuracy matched to next-to-next-to-leading fixed-order in perturbative Quantum Chromodynamics. Once measured at the Large Hadron Collider, precision comparisons between theory and data can be made, essentially independent of the amount of pile-up contamination.

preprint2015arXiv

Hard-Soft-Collinear Factorization to All Orders

We provide a precise statement of hard-soft-collinear factorization of scattering amplitudes and prove it to all orders in perturbation theory. Factorization is formulated as the equality at leading power of scattering amplitudes in QCD with other amplitudes in QCD computed from a product of operator matrix elements. The equivalence is regulator independent and gauge independent. As the formulation relates amplitudes to the same amplitudes with additional soft or collinear particles, it includes as special cases the factorization of soft currents and collinear splitting functions from generic matrix elements, both of which are shown to be process independent to all orders. We show that the overlapping soft-collinear region is naturally accounted for by vacuum matrix elements of kinked Wilson lines. Although the proof is self-contained, it combines techniques developed for the study of pinch surfaces, scattering amplitudes, and effective field theory.

preprint2015arXiv

Removing phase-space restrictions in factorized cross sections

Factorization in gauge theories holds at the amplitude or amplitude-squared level for states of given soft or collinear momenta. When performing phase-space integrals over such states, one would generally like to avoid putting in explicit cuts to separate soft from collinear momenta. Removing these cuts induces an overcounting of the soft-collinear region and adds new infrared-ultraviolet divergences in the collinear region. In this paper, we first present a regulator-independent subtraction algorithm for removing soft-collinear overlap at the amplitude level which may be useful in pertubative QCD. We then discuss how both the soft-collinear and infrared-ultraviolet overlap can be undone for certain observables in a way which respects factorization. Our discussion clarifies some of the subtleties in phase-space subtractions and includes a proof of the infrared finiteness of a suitably subtracted jet function. These results complete the connection between factorized QCD and Soft-Collinear Effective Theory.

preprint2015arXiv

Streamlining resummed QCD calculations using Monte Carlo integration

Some of the most arduous and error-prone aspects of precision resummed calculations are related to the partonic hard process, having nothing to do with the resummation. In particular, interfacing to parton-distribution functions, combining various channels, and performing the phase space integration can be limiting factors in completing calculations. Conveniently, however, most of these tasks are already automated in many Monte Carlo programs, such as MadGraph, Alpgen or Sherpa. In this paper, we show how such programs can be used to produce distributions of partonic kinematics with associated color structures representing the hard factor in a resummed distribution. These distributions can then be used to weight convolutions of jet, soft and beam functions producing a complete resummed calculation. In fact, only around 1000 unweighted events are necessary to produce precise distributions. A number of examples and checks are provided, including $e^+e^-$ two- and four-jet event shapes, $n$-jettiness and jet-mass related observables at hadron colliders. Attached code can be used to modify MadGraph to export the relevant leading-order hard functions and color structures for arbitrary processes.

preprint2014arXiv

Consistent Use of Effective Potentials

It is well known that effective potentials can be gauge-dependent while their values at extrema should be gauge-invariant. Unfortunately, establishing this invariance in perturbation theory is not straightforward, since contributions from arbitrarily high- order loops can be of the same size. We show in massless scalar QED that an infinite class of loops can be summed (and must be summed) to give a gauge invariant value for the potential at its minimum. In addition, we show that the exact potential depends on both the scale at which it is calculated and the normalization of the fields, but the vacuum energy does not. Using these insights, we propose a method to extract some physical quantities from effective potentials which is self-consistent order-by-order in perturbation theory, including improvement with the renormalization group.

preprint2014arXiv

Consistent Use of the Standard Model Effective Potential

The stability of the Standard Model is determined by the true minimum of the effective Higgs potential. We show that the potential at its minimum when computed by the traditional method is strongly dependent on the gauge parameter. It moreover depends on the scale where the potential is calculated. We provide a consistent method for determining absolute stability independent of both gauge and calculation scale, order by order in perturbation theory. This leads to a revised stability bounds mH > (129.4 \pm 2.3) GeV and mt < (171.2 \pm 0.3)GeV. We also show how to evaluate the effect of new physics on the stability bound without resorting to unphysical field values.

preprint2014arXiv

Jet Cleansing: Pileup Removal at High Luminosity

One of the greatest impediments to extracting useful information from high luminosity hadron-collider data is radiation from secondary collisions (i.e. pileup) which can overlap with that of the primary interaction. In this paper we introduce a simple jet-substructure technique termed cleansing which can consistently correct for large amounts of pileup in an observable independent way. Cleansing works at the subjet level, combining tracker and calorimeter-based data to reconstruct the pileup-free primary interaction. The technique can be used on its own, with various degrees of sophistication, or in concert with jet grooming. We apply cleansing to both kinematic and jet shape reconstruction, finding in all cases a marked improvement over previous methods both in the correlation of the cleansed data with uncontaminated results and in measures like S/rt(B). Cleansing should improve the sensitivity of new-physics searches at high luminosity and could also aid in the comparison of precision QCD calculations to collider data.

preprint2014arXiv

Non-global Logarithms at 3 Loops, 4 Loops, 5 Loops and Beyond

We calculate the coefficients of the leading non-global logarithms for the hemisphere mass distribution analytically at 3, 4, and 5 loops at large Nc . We confirm that the integrand derived with the strong-energy-ordering approximation and fixed-order iteration of the Banfi-Marchesini-Syme (BMS) equation agree. Our calculation exploits a hidden PSL(2,R) symmetry associated with the jet directions, apparent in the BMS equation after a stereographic projection to the Poincare disk. The required integrals have an iterated form, leading to functions of uniform transcendentality. This allows us to extract the coefficients, and some functional dependence on the jet directions, by computing the symbols and coproducts of appropriate expressions involving classical and Goncharov polylogarithms. Convergence of the series to a numerical solution of the BMS equation is also discussed.

preprint2013arXiv

An on-shell approach to factorization

Factorization is possible due to the universal behavior of Yang-Mills theories in soft and collinear limits. Here, we take a small step towards a more transparent understanding of these limits by proving a form of perturbative factorization at tree- level using on-shell spinor helicity methods. We present a concrete and self-contained expression of factorization in which matrix elements in QCD are related to products of other matrix elements in QCD up to leading order in a power-counting parameter determined by the momenta of certain physical on-shell states. Our approach uses only the scaling of momenta in soft and collinear limits, avoiding any assignment of scaling behavior to unphysical (and gauge-dependent) fields. The proof of factorization exploits many advantages of helicity spinors, such as the freedom to choose different reference vectors for polarizations in different collinear sectors. An advantage of this approach is that once factorization is shown to hold in QCD, the transition to Soft-Collinear Effective Theory is effortless.

preprint2013arXiv

Jet Charge at the LHC

Knowing the charge of the parton initiating a light-quark jet could be extremely useful both for testing aspects of the Standard Model and for characterizing potential beyond-the-Standard-Model signals. We show that despite the complications of hadronization and out-of-jet radiation such as pile-up, a weighted sum of the charges of a jet's constituents can be used at the LHC to distinguish among jets with different charges. Potential applications include measuring electroweak quantum numbers of hadronically decaying resonances or supersymmetric particles, as well as Standard Model tests, such as jet charge in dijet events or in hadronically-decaying W bosons in t-tbar events. We develop a systematically improvable method to calculate moments of these charge distributions by combining multi-hadron fragmentation functions with perturbative jet functions and pertubative evolution equations. We show that the dependence on energy and jet size for the average and width of the jet charge can be calculated despite the large experimental uncertainty on fragmentation functions. These calculations can provide a validation tool for data independent of Monte-Carlo fragmentation models.

preprint2013arXiv

Jet Sampling: Improving Event Reconstruction through Multiple Interpretations

The classification of events involving jets as signal-like or background-like can depend strongly on the jet algorithm used and its parameters. This is partly due to the fact that standard jet algorithms yield a single partition of the particles in an event into jets, even if no particular choice stands out from the others. As an alternative, we propose that one should consider multiple interpretations of each event, generalizing the Qjets procedure to event-level analysis. With multiple interpretations, an event is no longer restricted to either satisfy cuts or not satisfy them - it can be assigned a weight between 0 and 1 based on how well it satisfies the cuts. These cut-weights can then be used to improve the discrimination power of an analysis or reduce the uncertainty on mass or cross-section measurements. For example, using this approach on a Higgs plus Z boson sample, with h->bb we find an 28% improvement in significance can be realized at the 8 TeV LHC. Through a number of other examples, we show various ways in which having multiple interpretations can be useful on the event level.

preprint2013arXiv

Precision direct photon and W-boson spectra at high p_T and comparison to LHC data

The differential p_T spectrum for vector boson production is computed at next-to-leading fixed order and including the resummation of threshold logarithms at next-to-next-to-leading logarithmic accuracy. A comparison is made to ATLAS data on direct photon and W production at high transverse momentum p_T, finding excellent agreement. The resummation is achieved by factorizing contributions associated with different scales using Soft-Collinear Effective Theory. Each part is then calculated perturbatively and the individual contributions are combined using renormalization group methods. A key advantage of the effective theory framework is that it indicates a set of natural scale choices, in contrast to the fixed-order calculation. Resummation of logarithms of ratios of these scales leads to better agreement with data and reduced theoretical uncertainties.

preprint2013arXiv

Quark and Gluon Jet Substructure

Distinguishing quark-initiated jets from gluon-initiated jets has the potential to significantly improve the reach of many beyond-the-standard model searches at the Large Hadron Collider and to provide additional tests of QCD. To explore whether quark and gluon jets could possibly be distinguished on an event-by-event basis, we perform a comprehensive simulation-based study. We explore a variety of motivated and unmotivated variables with a semi-automated multivariate approach. General conclusions are that at 50% quark jet acceptance efficiency, around 80%-90% of gluon jets can be rejected. Some benefit is gained by combining variables. Different event generators are compared, as are the effects of using only charged tracks to avoid pileup. Additional information, including interactive distributions of most variables and their cut efficiencies, can be at http://jets.physics.harvard.edu/qvg.

preprint2012arXiv

1-loop matching and NNLL resummation for all partonic 2 to 2 processes in QCD

The Wilson Coefficients for all 4-parton operators which arise in matching QCD to Soft-Collinear Effective Theory (SCET) are computed at 1-loop. Any dijet observable calculated in SCET beyond leading order will require these results. The Wilson coefficients are separated by spin and color, although most applications will involve only the spin-averaged hard functions. The anomalous dimensions for the Wilson coefficients are given to 2-loop order, and the renormalization group equations are solved explicitly. This will allow for analytical resummation of dijet observables to next-to-next-to-leading logarithmic accuracy. For each channel, there is a natural basis in which the evolution is diagonal in color space. The same basis also diagonalizes the color evolution for the soft function. Even though soft functions required for SCET calculations are observable dependent, it is shown that their renormalization group evolution is almost completely determined by a universal structure. With these results, it will be possible to calculate hadronic event shapes or other dijet observables to next-to-leading order with next-to-next-to-leading log resummation.

preprint2012arXiv

Direct photon production with effective field theory

The production of hard photons in hadronic collisions is studied using Soft-Collinear Effective Theory (SCET). This is the first application of SCET to a physical, observable cross section involving energetic partons in more than two directions. A factorization formula is derived which involves a non-trivial interplay of the angular dependence in the hard and soft functions, both quark and gluon jet functions, and multiple partonic channels. The relevant hard, jet and soft functions are computed to one loop and their anomalous dimensions are determined to three loops. The final resummed inclusive direct photon distribution is valid to next-to-next-to-leading logarithmic order (NNLL), one order beyond previous work. The result is improved by including non-logarithmic terms and photon isolation cuts through matching, and compared to Tevatron data and to fixed order results at the Tevatron and the LHC. The resummed cross section has a significantly smaller theoretical uncertainty than the next-to-leading fixed-order result, particularly at high transverse momentum.

preprint2012arXiv

Jet Physics from Static Charges in AdS

Soft interactions with high-energy jets are explored in radial coordinates which exploit the approximately conformal behavior of perturbative gauge theories. In these coordinates, the jets, approximated by Wilson lines, become static charges in Euclidean AdS. The anomalous dimension of the corresponding Wilson line operator is then determined by the potential energy of the charges. To study these Wilson lines we introduce a "conformal gauge" which does not have kinetic mixing between radial and angular directions, and show that a number of properties of Wilson lines are reproduced through relatively simple calculations. For example, certain non-planar graphs involving multiple Wilson lines automatically vanish. We also discuss the linear growth of the charges' imaginary potential energy with separation, and a relationship between Wilson line diagrams and Witten diagrams.

preprint2012arXiv

Precision Jet Substructure from Boosted Event Shapes

Jet substructure has emerged as a critical tool for LHC searches, but studies so far have relied heavily on shower Monte Carlo simulations, which formally approximate QCD at leading-log level. We demonstrate that systematic higher-order QCD computations of jet substructure can be carried out by boosting global event shapes by a large momentum Q, and accounting for effects due to finite jet size, initial-state radiation (ISR), and the underlying event (UE) as 1/Q corrections. In particular, we compute the 2-subjettiness substructure distribution for boosted Z -> q qbar events at the LHC at next-to-next-to-next-to-leading-log order. The calculation is greatly simplified by recycling the known results for the thrust distribution in e+ e- collisions. The 2-subjettiness distribution quickly saturates, becoming Q independent for Q > 400 GeV. Crucially, the effects of jet contamination from ISR/UE can be subtracted out analytically at large Q, without knowing their detailed form. Amusingly, the Q=infinity and Q=0 distributions are related by a scaling by e, up to next-to-leading-log order.

preprint2012arXiv

Qjets: A Non-Deterministic Approach to Tree-Based Jet Substructure

Jet substructure is typically studied using clustering algorithms, such as kT, which arrange the jets' constituents into trees. Instead of considering a single tree per jet, we propose that multiple trees should be considered, weighted by an appropriate metric. Then each jet in each event produces a distribution for an observable, rather than a single value. Advantages of this approach include: 1) observables have significantly increased statistical stability; and, 2) new observables, such as the variance of the distribution, provide new handles for signal and background discrimination. For example, we find that employing a set of trees substantially reduces the observed fluctuations in the pruned mass distribution, enhancing the likelihood of new particle discovery for a given integrated luminosity. Furthermore, the resulting pruned mass distributions for (background) QCD jets are found to be substantially wider than that for (signal) jets with intrinsic mass scales, e.g. jets containing a W decay. A cut on this width yields a substantial enhancement in significance relative to a cut on the standard pruned jet mass alone. In particular the luminosity needed for a given significance requirement decreases by a factor of two relative to standard pruning.

preprint2012arXiv

Resummation of Jet Mass at Hadron Colliders

A method is developed for calculating the jet mass distribution at hadron colliders using an expansion about the kinematic threshold. In particular, we consider the mass distribution of jets of size R produced in association with a hard photon at the Large Hadron Collider. Expanding around the kinematic threshold, where all the energy goes into the jet and the photon, provides a clean factorization formula and allows for the resummation of logarithms associated with soft and collinear divergences. All of the large logarithms of jet mass are resummed at next-to-leading logarithmic level (NLL), and all the global logarithms at next-to-next-to-leading logarithmic level (NNLL). A key step in the derivation is the factorization of the soft function into pieces associated with single scales and a remainder which contains non-global structure. This step, which is standard in traditional resummation, is implemented in effective field theory which is then used to resum the large logarithms using the renormalization group in a systematically improvable manner.

preprint2011arXiv

Jet mass with a jet veto at two loops and the universality of non-global structure

We investigate the exclusive jet mass distribution in e+e- events, defined with a veto on the out-of-jet radiation, at two-loop order. In particular, we calculate the two-loop soft function, which is required to describe this distribution in the threshold region. When combined with other ingredients using soft-collinear effective theory, this generates the complete singular distribution for jet thrust, the sum of the jet masses, at two-loop order. The result is in excellent agreement with full QCD. The integrated jet thrust distribution is found to depend in an intricate way on both the finite jet cone size, R, and the jet veto scale. The result clarifies the structure of the potentially large logarithms (both global and non- global) which arise in jet observables for the first time at this order. Somewhat surprisingly, we find that, in the small R limit, there is a precise and simple correspondence between the non-global contribution to the integrated jet thrust distribution and the previously calculated non-global contribution to the integrated hemisphere soft function, including subleading terms. This suggests that the small R limit may provide a useful expansion for studying other exclusive jet substructure observables.

preprint2011arXiv

Multivariate discrimination and the Higgs + W/Z search

A systematic method for optimizing multivariate discriminants is developed and applied to the important example of a light Higgs boson search at the Tevatron and the LHC. The Significance Improvement Characteristic (SIC), defined as the signal efficiency of a cut or multivariate discriminant divided by the square root of the background efficiency, is shown to be an extremely powerful visualization tool. SIC curves demonstrate numerical instabilities in the multivariate discriminants, show convergence as the number of variables is increased, and display the sensitivity to the optimal cut values. For our application, we concentrate on Higgs boson production in association with a W or Z boson with H -> bb and compare to the irreducible standard model background, Z/W + bb. We explore thousands of experimentally motivated, physically motivated, and unmotivated single variable discriminants. Along with the standard kinematic variables, a number of new ones, such as twist, are described which should have applicability to many processes. We find that some single variables, such as the pull angle, are weak discriminants, but when combined with others they provide important marginal improvement. We also find that multiple Higgs boson-candidate mass measures, such as from mild and aggressively trimmed jets, when combined may provide additional discriminating power. Comparing the significance improvement from our variables to those used in recent CDF and DZero searches, we find that a 10-20% improvement in significance against Z/W + bb is possible. Our analysis also suggests that the H + W/Z channel with H -> bb is also viable at the LHC, without requiring a hard cut on the W/Z transverse momentum.

preprint2011arXiv

Pure Samples of Quark and Gluon Jets at the LHC

Having pure samples of quark and gluon jets would greatly facilitate the study of jet properties and substructure, with many potential standard model and new physics applications. To this end, we consider multijet and jets+X samples, to determine the purity that can be achieved by simple kinematic cuts leaving reasonable production cross sections. We find, for example, that at the 7 TeV LHC, the pp {\to} γ+2jets sample can provide 98% pure quark jets with 200 GeV of transverse momentum and a cross section of 5 pb. To get 10 pb of 200 GeV jets with 90% gluon purity, the pp {\to} 3jets sample can be used. b+2jets is also useful for gluons, but only if the b-tagging is very efficient.

preprint2011arXiv

Quark and Gluon Tagging at the LHC

Being able to distinguish light-quark jets from gluon jets on an event-by-event basis could significantly enhance the reach for many new physics searches at the Large Hadron Collider. Through an exhaustive search of existing and novel jet substructure observables, we find that a multivariate approach can filter out over 95% of the gluon jets while keeping more than half of the light-quark jets. Moreover, a combination of two simple variables, the charge track multiplicity and the $p_T$-weighted linear radial moment (girth), can achieve similar results. While this pair appears very promising, our study is only Monte Carlo based, and other discriminants may work better with real data in a realistic experimental environment. To that end, we explore many other observables constructed using different jet sizes and parameters, and highlight those that deserve further theoretical and experimental scrutiny. Additional information, including distributions of around 10,000 variables, can be found on this website http://jets.physics.harvard.edu/qvg .

preprint2011arXiv

Resummation for W and Z production at large pT

Soft-Collinear Effective theory is used to perform threshold resummation for W and Z production at large transverse momentum to next-to-next-to-leading logarithmic accuracy including matching to next-to-leading fixed-order results. The results agree very well with data from the Tevatron, and predictions are made for the high-pT spectra at the LHC. While the higher-log terms are of moderate size, their inclusion leads to a substantial reduction of the perturbative uncertainty. With these improvements, the PDF uncertainties now dominate the error on the predicted cross section.

preprint2011arXiv

Resummation of jet mass with and without a jet veto

Calculating the distribution of jet masses in high-energy collisions is challenging because fixed-order perturbation theory breaks down near the peak region, and because multiple scales complicate the resummation. To avoid using a jet veto, one can consider inclusive observables, in which every particle is in a jet. We demonstrate that calculating the mass of the hardest jet in multijet events can be problematic, and we give an example of an inclusive observable, asymmetric thrust, which can be resummed to next-to-next-to-leading logarithmic accuracy. Exclusive observables with out-of-jet regions are more complicated. Even for e+e- dijet events at energy Q, to calculate the mass m of jets of size R, one must impose a veto on the energy omega of extra jets to force dijet kinematics; then there are both log m/Q and log m/omega singularities. To proceed, we suggest a refactorization of the soft function in the small R limit. To justify this refactorization, we show that the expansion of the resummed distribution is in excellent agreement with fixed order. This motivates considering the expansion around small R as a useful handle on producing phenomenologically useful resummed jet mass distributions. The strong evidence we give for refactorization at small R is independent of non-global logarithms, which are not the subject of this paper.

preprint2011arXiv

Simplified Models for LHC New Physics Searches

This document proposes a collection of simplified models relevant to the design of new-physics searches at the LHC and the characterization of their results. Both ATLAS and CMS have already presented some results in terms of simplified models, and we encourage them to continue and expand this effort, which supplements both signature-based results and benchmark model interpretations. A simplified model is defined by an effective Lagrangian describing the interactions of a small number of new particles. Simplified models can equally well be described by a small number of masses and cross-sections. These parameters are directly related to collider physics observables, making simplified models a particularly effective framework for evaluating searches and a useful starting point for characterizing positive signals of new physics. This document serves as an official summary of the results from the "Topologies for Early LHC Searches" workshop, held at SLAC in September of 2010, the purpose of which was to develop a set of representative models that can be used to cover all relevant phase space in experimental searches. Particular emphasis is placed on searches relevant for the first ~50-500 pb-1 of data and those motivated by supersymmetric models. This note largely summarizes material posted at http://lhcnewphysics.org/, which includes simplified model definitions, Monte Carlo material, and supporting contacts within the theory community. We also comment on future developments that may be useful as more data is gathered and analyzed by the experiments.

preprint2011arXiv

The two-loop hemisphere soft function

The hemisphere soft function is calculated to order alpha_s^2. This is the first multi-scale soft function calculated to two loops. The renormalization scale dependence of the result agrees exactly with the prediction from effective field theory. This fixes the unknown coefficients of the singular parts of the two-loop thrust and heavy-jet mass distributions. There are four such coefficients, for 2 event shapes and 2 color structures, which are shown to be in excellent agreement with previous numerical extraction. The asymptotic behavior of the soft function has double logs in the CF CA color structure, which agree with non-global log calculations, but also has sub-leading single logs for both the CF CA and CF TF nf color structures. The general form of the soft function is complicated, does not factorize in a simple way, and disagrees with the Hoang-Kluth ansatz. The exact hemisphere soft function will remove one source of uncertainty on the alpha_s fits from e+e- event shapes.

preprint2011arXiv

Threshold Hadronic Event Shapes with Effective Field Theory

Hadronic event shapes, that is, event shapes at hadron colliders, could provide a great way to test both standard and non-standard theoretical models. However, they are significantly more complicated than event shapes at e+e- colliders, involving multiple hard directions, multiple channels and multiple color structures. In this paper, hadronic event shapes are examined with Soft-Collinear Effective Theory (SCET) by expanding around the dijet limit. A simple event shape, threshold thrust, is defined. This observable is global and has no free parameters, making it ideal for clarifying how resummation of hadronic event shapes can be done in SCET. Threshold thrust is calculated at next-to-leading fixed order (NLO) in SCET and resummed to next-to-next-to-leading logarithmic accuracy (NNLL). The scale-dependent parts of the soft function are shown to agree with what is expected from general observations, and the factorization formula is explicitly shown to be renormalization group invariant to 1-loop. Although threshold thrust is not itself expected to be phenomenologically interesting, it can be modified into a related observable which allows the jet pT distribution to be calculated and resummed to NNLL+NLO accuracy. As in other processes, one expects resummation to be important even for moderate jet momenta due to dynamical threshold enhancement. A general discussion of threshold enhancement and non-global logs in hadronic event shapes is also included.

preprint2011arXiv

Top condensation as a motivated explanation of the top forward-backward asymmetry

Models of top condensation can provide both a compelling solution to the hierarchy problem as well as an explanation of why the top-quark mass is large. The spectrum of such models, in particular topcolor-assisted technicolor, includes top-pions, top-rhos and the top-Higgs, all of which can easily have large top-charm or top-up couplings. Large top-up couplings in particular would lead to a top forward-backward asymmetry through $t$-channel exchange, easily consistent with the Tevatron measurements. Intriguingly, there is destructive interference between the top-mesons and the standard model which conspire to make the overall top pair production rate consistent with the standard model. The rate for same-sign top production is also small due to destructive interference between the neutral top-pion and the top-Higgs. Flavor physics is under control because new physics is mostly confined to the top quark. In this way, top condensation can explain the asymmetry and be consistent with all experimental bounds. There are many additional signatures of topcolor with large tu mixing, such as top(s)+jet(s) events, in which a top and a jet reconstruct a resonance mass, which make these models easily testable at the LHC.

preprint2011arXiv

W-jet Tagging: Optimizing the Identification of Boosted Hadronically-Decaying W Bosons

A method is proposed for distinguishing highly boosted hadronically decaying W's (W-jets) from QCD-jets using jet substructure. Previous methods, such as the filtering/mass-drop method, can give a factor of ~2 improvement in S/sqrt(B) for jet pT > 200 GeV. In contrast, a multivariate approach including new discriminants such as R-cores, which characterize the shape of the W-jet, subjet planar flow, and grooming-sensitivities is shown to provide a much larger factor of ~5 improvement in S/sqrt(B). For longitudinally polarized W's, such as those coming from many new physics models, the discrimination is even better. Comparing different Monte Carlo simulations, we observe a sensitivity of some variables to the underlying event; however, even with a conservative estimates, the multivariate approach is very powerful. Applications to semileptonic WW resonance searches and all-hadronic W+jet searches at the LHC are also discussed. Code implementing our W-jet tagging algorithm is publicly available at http://jets.physics.harvard.edu/wtag

preprint2010arXiv

Annihilation decays of bound states at the LHC

At the Large Hadron Collider, heavy particles may be produced in pairs close to their kinematic threshold. If these particles have strong enough attractive interactions they may form bound states. Consequently, the bound states may decay through annihilation back into the standard model. Such annihilation decays have the potential to provide much information about the bound particles, such as their mass, spin, or charges, in a manner completely complementary to standard single particle cascade decays. Many of the signatures, such as dijet resonances, will be challenging to find, but may be extremely helpful in unraveling the nature of the new physics. In the standard model, the only novel annihilation decays would be for toponium; these will be hard to see because of the relatively large width of the top quark itself. In models with supersymmetry, marginally visible annihilation decays may occur for example, from bound states of gluinos to dijets or tops. If new particles are bound through forces stronger than QCD, annihilation decays may even be the discovery mode for new physics. This paper presents various theoretical results about bound states and then addresses the practical question of whether any of their annihilation decays can be seen at the LHC.

preprint2010arXiv

Resummation of heavy jet mass and comparison to LEP data

The heavy jet mass distribution in e+e- collisions is computed to next-to-next-to-next-to leading logarithmic (NNNLL) and next-to-next-to leading fixed order accuracy (NNLO). The singular terms predicted from the resummed distribution are confirmed by the fixed order distributions allowing a precise extraction of the unknown soft function coefficients. A number of quantitative and qualitative comparisons of heavy jet mass and the related thrust distribution are made. From fitting to ALEPH data, a value of alpha_s is extracted, alpha_s(m_Z)=0.1220 +/- 0.0031, which is larger than, but not in conflict with, the corresponding value for thrust. A weighted average of the two produces alpha_s(m_Z) = 0.1193 +/- 0.0027, consistent with the world average. A study of the non-perturbative corrections shows that the flat direction observed for thrust between alpha_s and a simple non-perturbative shape parameter is not lifted in combining with heavy jet mass. The Monte Carlo treatment of hadronization gives qualitatively different results for thrust and heavy jet mass, and we conclude that it cannot be trusted to add power corrections to the event shape distributions at this accuracy. Whether a more sophisticated effective field theory approach to power corrections can reconcile the thrust and heavy jet mass distributions remains an open question.

preprint2010arXiv

Seeing in Color: Jet Superstructure

A new class of observables is introduced which aims to characterize the superstructure of an event, that is, features, such as color flow, which are not determined by the jet four-momenta alone. Traditionally, an event is described as having jets which are independent objects; each jet has some energy, size, and possible substructure such as subjets or heavy flavor content. This description discards information connecting the jets to each other, which can be used to determine if the jets came from decay of a color singlet object, or if they were initiated by quarks or gluons. An example superstructure variable, pull, is presented as a simple handle on color flow. It can be used on an event-by-event basis as a tool for distinguishing previously irreducible backgrounds at the Tevatron and the LHC.

preprint2007arXiv

Tensors Mesons in AdS/QCD

We explore tensor mesons in AdS/QCD focusing on f2 (1270), the lightest spin-two resonance in QCD. We find that the f2 mass and the partial width for f2 -> gamma gamma are in very good agreement with data. In fact, the dimensionless ratio of these two quantities comes out within the current experimental bound. The result for this ratio depends only on Nc and Nf, and the quark and glueball content of the operator responsible for the f2; more importantly, it does not depend on chiral symmetry breaking and so is both independent of much of the arbitrariness of AdS/QCD and completely out of reach of chiral perturbation theory. For comparison, we also explore f2 -> pi pi, which because of its sensitivity to the UV corrections has much more uncertainty. We also calculate the masses of the higher spin resonances on the Regge trajectory of the f2, and find they compare favorably with experiment.

preprint2004arXiv

Massive Supergravity and Deconstruction

We present a simple superfield Lagrangian for massive supergravity. It comprises the minimal supergravity Lagrangian with interactions as well as mass terms for the metric superfield and the chiral compensator. This is the natural generalization of the Fierz-Pauli Lagrangian for massive gravity which comprises mass terms for the metric and its trace. We show that the on-shell bosonic and fermionic fields are degenerate and have the appropriate spins: 2, 3/2, 3/2 and 1. We then study this interacting Lagrangian using goldstone superfields. We find that a chiral multiplet of goldstones gets a kinetic term through mixing, just as the scalar goldstone does in the non-supersymmetric case. This produces Planck scale (Mpl) interactions with matter and all the discontinuities and unitarity bounds associated with massive gravity. In particular, the scale of strong coupling is (Mpl m^4)^1/5, where m is the multiplet's mass. Next, we consider applications of massive supergravity to deconstruction. We estimate various quantum effects which generate non-local operators in theory space. As an example, we show that the single massive supergravity multiplet in a 2-site model can serve the function of an extra dimension in anomaly mediation.

preprint2004arXiv

Strong CP, Flavor, and Twisted Split Fermions

We present a natural solution to the strong CP problem in the context of split fermions. By assuming CP is spontaneously broken in the bulk, a weak CKM phase is created in the standard model due to a twisting in flavor space of the bulk fermion wavefunctions. But the strong CP phase remains zero, being essentially protected by parity in the bulk and CP on the branes. As always in models of spontaneous CP breaking, radiative corrections to theta bar from the standard model are tiny, but even higher dimension operators are not that dangerous. The twisting phenomenon was recently shown to be generic, and not to interfere with the way that split fermions naturally weaves small numbers into the standard model. It follows that out approach to strong CP is compatible with flavor, and we sketch a comprehensive model. We also look at deconstructed version of this setup which provides a viable 4D model of spontaneous CP breaking which is not in the Nelson-Barr class.

preprint2004arXiv

Twisted Split Fermions

The observed flavor structure of the standard model arises naturally in "split fermion" models which localize fermions at different places in an extra dimension. It has, until now, been assumed that the bulk masses for such fermions can be chosen to be flavor diagonal simultaneously at every point in the extra dimension, with all the flavor violation coming from the Yukawa couplings to the Higgs. We consider the more natural possibility in which the bulk masses cannot be simultaneously diagonalized, that is, that they are twisted in flavor space. We show that, in general, this does not disturb the natural generation of hierarchies in the flavor parameters. Moreover, it is conceivable that all the flavor mixing and CP-violation in the standard model may come only from twisting, with the five-dimensional Yukawa couplings taken to be universal.

preprint2002arXiv

Effective Field Theory for Massive Gravitons and Gravity in Theory Space

We introduce a technique for restoring general coordinate invariance into theories where it is explicitly broken. This is the analog for gravity of the Callan-Coleman-Wess-Zumino formalism for gauge theories. We use this to elucidate the properties of interacting massless and massive gravitons. For a single graviton with a Planck scale Mpl and a mass mg, we find that there is a sensible effective field theory which is valid up to a high-energy cutoff Lambda parametrically above mg. Our methods allow for a transparent understanding of the many peculiarities associated with massive gravitons, among them the need for the Fierz-Pauli form of the Lagrangian, the presence or absence of the van Dam-Veltman-Zakharov discontinuity in general backgrounds, and the onset of non-linear effects and the breakdown of the effective theory at large distances from heavy sources. The natural sizes of all non-linear corrections beyond the Fierz-Pauli term are easily determined. The cutoff scales as Lambda ~ (mg^4 Mpl)^(1/5) for the Fierz-Pauli theory, but can be raised to Lambda ~ (mg^2 Mpl)^(1/3) in certain non-linear extensions. Having established that these models make sense as effective theories, there are a number of new avenues for exploration, including model building with gravity in theory space and constructing gravitational dimensions.