• We calculate the soft function for the global event variable 1-jettiness at next-to-next-to-leading order (NNLO) in QCD. We focus specifically on the non-Abelian contribution, which, unlike the Abelian part, is not determined by the next-to-leading order result. The calculation uses the known general forms for the emission of one and two soft partons and is performed using a sector-decomposition method that is spelled out in detail. Results are presented in the form of numerical fits to the 1-jettiness soft function for LHC kinematics (as a function of the angle between the incoming beams and the final-state jet) and for generic kinematics (as a function of three independent angles). These fits represent one of the needed ingredients for NNLO calculations that use the N-jettiness event variable to handle infrared singularities.
  • In this paper we present a calculation of the $\gamma+j$ process at next-to-next-to-leading order (NNLO) in QCD and compare the resulting predictions to 8 TeV CMS data. We find good agreement with the shape of the photon $p_T$ spectrum, particularly after the inclusion of additional electroweak corrections, but there is a tension between the overall normalization of the theoretical prediction and the measurement. We use our results to compute the ratio of $Z(\to \ell^+\ell^-)+j$ to $\gamma+j$ events as a function of the vector boson transverse momentum at NNLO, a quantity that is used to normalize $Z(\rightarrow\nu\overline{\nu}) +j$ backgrounds in searches for dark matter and supersymmetry. Our NNLO calculation significantly reduces the theoretical uncertainty on this ratio, thus boosting its power for future searches of new physics.
  • We present the first calculation of direct photon production at next-to-next-to leading order (NNLO) accuracy in QCD. For this process, although the final state cuts mandate only the presence of a single electroweak boson, the underlying kinematics resembles that of a generic vector boson plus jet topology. In order to regulate the infrared singularities present at this order we use the $N$-jettiness slicing procedure, applied for the first time to a final state that at Born level includes colored partons but no required jet. We compare our predictions to ATLAS 8 TeV data and find that the inclusion of the NNLO terms in the perturbative expansion, supplemented by electroweak corrections, provides an excellent description of the data with greatly reduced theoretical uncertainties.
  • The sophistication of current predictions for $Z+$jet production at hadron colliders necessitates a re-evaluation of any approximations inherent in the theoretical calculations. In this paper we address one such issue, the inclusion of mass effects in top-quark loops. We ameliorate an existing calculation of $Z+1$~jet and $Z+2$~jet production by presenting exact analytic formulae for amplitudes containing top-quark loops that enter at next-to-leading order in QCD. Although approximations based on an expansion in powers of $1/m_t^2$ can lead to poor high-energy behavior, an exact treatment of top-quark loops demonstrates that their effect is small and has limited phenomenological interest.
  • We present a new release of the QCDLoop library based on a modern object-oriented framework. We discuss the available new features such as the extension to the complex masses, the possibility to perform computations in double and quadruple precision simultaneously, and useful caching mechanisms to improve the computational speed. We benchmark the performance of the new library, and provide practical examples of phenomenological implementations by interfacing this new library to Monte Carlo programs.
  • I describe the contributions of Guido Altarelli to the development of Quantum Chromodynamics from the discovery of asymptotic freedom until the end of the S$p\bar{p}$S collider era, 1973-1985.
  • We present results for the production of a pair of on-shell Z bosons via gluon fusion. This process occurs both through the production and decay of the Higgs boson, and through continuum production where the Z boson couples to a loop of massless quarks or to a massive quark. We calculate the interference of the two processes and its contribution to the cross section up to and including order O(alpha_s^3). The two-loop contributions to the amplitude are all known analytically, except for the continuum production through loops of top quarks of mass m. The latter contribution is important for the invariant mass of the two Z bosons, (as measured by the mass of their leptonic decay products, m_4l), in a regime where m_4l > 2m because of the contributions of longitudinal bosons. We examine all the contributions to the virtual amplitude involving top quarks, as expansions about the heavy top quark limit. Comparison with the analytic results, where known, allows us to assess the validity of the heavy quark expansion, and it extensions. We give results for the NLO corrections to this interference, including both real and virtual radiation.
  • In this paper we present a Next-to-Next-to Leading Order (NNLO) calculation of the production of a Higgs boson in association with a massive vector boson. We include the decays of the unstable Higgs and vector bosons, resulting in a fully flexible parton-level Monte Carlo implementation. We also include all $\mathcal{O}(\alpha_s^2)$ contributions that occur in production for these processes: those mediated by the exchange of a single off-shell vector boson in the $s$-channel, and those which arise from the coupling of the Higgs boson to a closed loop of fermions. We study final states of interest for Run II phenomenology, namely $H\rightarrow b\bar{b}$, $\gamma\gamma$ and $WW^*$. The treatment of the $H\rightarrow b\bar{b}$ decay includes QCD corrections at NLO. We use the recently developed $N$-jettiness regularization procedure, and study its viability in the presence of a large final-state phase space by studying $pp\rightarrow V(H\rightarrow WW^*) \rightarrow$ leptons.
  • In this paper we present a next-to-next-to-leading order (NNLO) calculation of the process $pp\rightarrow \gamma\gamma$ that we have implemented into the parton level Monte Carlo code MCFM. We do not find agreement with the previous calculation of this process in the literature. In addition to the $\mathcal{O}(\alpha_s^2)$ corrections present at NNLO, we include some effects arising at $\mathcal{O}(\alpha_s^3)$, namely those associated with gluon-initiated closed fermion loops. We investigate the role of this process in the context of studies of QCD at colliders and as a background for searches for new physics, paying particular attention to the diphoton invariant mass spectrum. We demonstrate that the NNLO QCD prediction for the shape of this spectrum agrees well with functional forms used in recent data-driven fits.
  • We present the first complete calculation of Z-boson production in association with a jet in hadronic collisions through next-to-next-to-leading order in perturbative QCD. Our computation uses the recently-proposed N-jettiness subtraction scheme to regulate the infrared divergences that appear in the real-emission contributions. We present phenomenological results for 13 TeV proton-proton collisions with fully realistic fiducial cuts on the final-state particles. The remaining theoretical uncertainties after the inclusion of our calculations are at the percent-level, making the Z+jet channel ready for precision studies at the LHC Run II.
  • We present the implementation of several color-singlet final-state processes at Next-to-Next-to Leading Order (NNLO) accuracy in QCD to the publicly available parton-level Monte Carlo program MCFM. Specifically we discuss the processes $pp\rightarrow H$, $pp\rightarrow Z$, $pp\rightarrow W$, $pp\rightarrow HZ$, $pp\rightarrow HW$ and $pp\rightarrow\gamma\gamma$. Decays of the unstable bosons are fully included, resulting in a flexible fully differential Monte Carlo code. The NNLO corrections have been calculated using the non-local $N$-jettiness subtraction approach. Special attention is given to the numerical aspects of running MCFM for these processes at this order. We pay particular attention to the systematic uncertainties due to the power corrections induced by the $N$-jettiness regularization scheme and the evaluation time needed to run the hybrid openMP/MPI version of MCFM at NNLO on multi-processor systems.
  • We present a next-to-leading order (NLO) calculation of $t\bar{t}$ production in hadronic collisions interfaced to shower generators according to the POWHEG method. We start from an NLO result from previous work, obtained in the zero width limit, where radiative corrections to both production and decays are included. The POWHEG interface required an extension of the POWHEG BOX framework, in order to deal with radiation from the decay of resonances. This extension is fully general (i.e. it can be applied in principle to any process considered in the zero width limit), and is here applied for the first time. In order to perform a realistic simulation, we introduce finite width effects using different approximations, that we validated by comparing with published exact NLO results. We have interfaced our POWHEG code to the PYTHIA8 shower Monte Carlo generator. At this stage, we dealt with novel issues related to the treatment of resonances, especially with regard to the initial scale for the shower that needs to be set appropriately. This procedure affects, for example, the fragmentation function of the b quark, that we have studied with particular attention. We believe that the tool presented here improves over previous generators for all aspects that have to do with top decays, and especially for the study of issues related to top mass measurements that involve B hadrons or b jets. The work presented here also constitutes a first step towards a fully consistent matching of NLO calculations involving intermediate resonances decaying into coloured particles, with parton showers.
  • We report on our findings modifying MCFM using OpenMP to implement multi-threading. By using OpenMP, the modified MCFM will execute on any processor, automatically adjusting to the number of available threads. We modified the integration routine VEGAS to distribute the event evaluation over the threads, while combining all events at the end of every iteration to optimize the numerical integration. Special care has been taken that the results of the Monte Carlo integration are independent of the number of threads used, to facilitate the validation of the OpenMP version of MCFM.
  • We present results on 4-lepton + 2-jet production, the partonic processes most commonly described as vector boson pair production in the Vector Boson Fusion (VBF) mode. This final state contains diagrams that are mediated by Higgs boson exchange. We focus particularly on the high-mass behaviour of the Higgs boson mediated diagrams, which unlike on-shell production, gives information about the Higgs couplings without assumptions on the Higgs boson total width. We assess the sensitivity of the high-mass region to Higgs coupling strengths, considering all vector boson pair channels, $W^-W^+$, $W^\pm W^\pm$, $W^\pm Z$ and $ZZ$. Because of the small background, the most promising mode is $W^+ W^+$ which has sensitivity to Higgs couplings because of Higgs boson exchange in the $t$-channel. Using the Caola-Melnikov (CM) method, the off-shell couplings can be interpreted as bounds on the Higgs boson total width. We estimate the bound that can be obtained with current data, as well as the bounds that could be obtained at $\sqrt{s}=13$~TeV in the VBF channel for data samples of 100 and 300 fb$^{-1}$. The CM method has already been successfully applied in the gluon fusion (GGF) production channel. The VBF production channel gives important complementary information, because both production and decay of the Higgs boson occur already at tree graph level.
  • We study interference effects in the production channel ZZ + jet, in particular focusing on the role of the Higgs boson. This production channel receives contributions both from Higgs boson-mediated diagrams via the decay $H \to ZZ$ (signal diagrams), as well as from diagrams where the Z-bosons couple directly to a quark loop (background diagrams). We consider the partonic processes $gggZZ$ and $gq \bar{q} ZZ$ in which interference between signal and background diagrams first occurs. Since interference is primarily an off-resonant effect for the Higgs boson, we treat the Z-bosons as on-shell. Thus our analysis is limited to the region above threshold, where the invariant mass of the Z-pair, $m_{ZZ}$, satisfies the condition $m_{ZZ}>2m_Z$. In the region $m_{ZZ} > 300$ GeV we find that the interference in the ZZ + jet channel is qualitatively similar to interference in the inclusive ZZ channel. Moreover, the rates are sufficient to study these effects at the LHC once jet-binned data become available.
  • We present results for the Standard model description of the four-lepton production, mediated both by Higgs boson production and by other one-loop standard model processes. The description of four-lepton final states in MCFM v6.8 is reviewed, with special reference to the interference effects that can occur for identical species of leptons. We present results both for interference in the $l^+l^- l^- l^+$ and in the $l^-l^+ \nu_l \bar{\nu}_l$ final state. Prospects for further improvement in the theoretical description of four lepton production are also reviewed.
  • We investigate the potential of the process $gg \to H \to WW$ to provide bounds on the Higgs width. Recent studies using off-shell $H\rightarrow ZZ$ events have shown that Run 1 LHC data can constrain the Higgs width, $\Gamma_H < (25-45) \Gamma_{H}^{\rm SM}$. Using 20 fb-1 of 8 TeV ATLAS data, we estimate a bound on the Higgs boson width from the WW channel between $\Gamma_H < (100-500) \Gamma_H^{SM}$. The large spread in limits is due to the range of cuts applied in the existing experimental analysis. The stricter cuts designed to search for the on-shell Higgs boson limit the potential number of off-shell events, weakening the constraints. As some of the cuts are lifted the bounds improve. We show that there is potential in the high transverse mass region to produce upper bounds of the order of $(25-50) \Gamma_H^{SM}$, depending strongly on the level of systematic uncertainty that can be obtained. Thus, if these systematics can be controlled, a constraint on the Higgs boson width from the $H \to WW$ decay mode can complement a corresponding limit from $H \to ZZ$.
  • We revisit the hadronic production of the four-lepton final state, e^- e^+ \mu^- \mu^+, through the fusion of initial state gluons. This process is mediated by loops of quarks and we provide first full analytic results for helicity amplitudes that account for both the effects of the quark mass in the loop and off-shell vector bosons. The analytic results have been implemented in the Monte Carlo program MCFM and are both fast, and numerically stable in the region of low Z transverse momentum. We use our results to study the interference between Higgs-mediated and continuum production of four-lepton final states, which is necessary in order to obtain accurate theoretical predictions outside the Higgs resonance region. We have confirmed and extended a recent analysis of Caola and Melnikov that proposes to use a measurement of the off-shell region to constrain the total width of the Higgs boson. Using a simple cut-and-count method, existing LHC data should bound the width at the level of 25-45 times the Standard Model expectation. We investigate the power of using a matrix element method to construct a kinematic discriminant to sharpen the constraint. In our analysis the bound on the Higgs width is improved by a factor of about 1.6 using a simple cut on the MEM discriminant, compared to an invariant mass cut m_{4l} > 300 GeV.
  • We present results for the production of a Z boson in association with single top at next-to-leading order (NLO), including the decay of the top quark and the Z boson. This electroweak process gives rise to the trilepton signature l+l-l'+- + jets + missing energy. We present results for this signature and show that the rate is competitive with the contribution of the mixed strong and electroweak production process, ttZ. As such it should be observable in the full data sample from LHC running at $\sqrt{s}=8$ TeV. The single top + Z process is a hitherto unconsidered irreducible background in searches for flavour changing neutral current decays of the top quark in tt production. For a selection of cuts used at the LHC involving a b-tag it is the dominant background. In an appendix we also briefly discuss the impact of NLO corrections on the related tH process.
  • In this work we present the implementation of generators for W and Z bosons in association with two jets interfaced to parton showers using the POWHEG BOX. We incorporate matrix elements from the parton-level Monte Carlo program MCFM in the POWHEG BOX, allowing for a considerable improvement in speed compared to previous implementations. We address certain problems that arise when processes that are singular at the Born level are implemented in a shower framework using either a generation cut or a Born suppression factor to yield weighted events. In such a case, events with very large weights can be generated after the shower through a number of mechanisms. Events with very small transverse momentum at the Born level can develop large transverse momentum either after the hardest emission, after the shower, or after the inclusion of multi-parton interactions. We present a solution to this problem that can be easily implemented in the POWHEG BOX. We also show that a full solution to this problem can only be achieved if the generator maintains physical validity also when the transverse momentum of the emitted partons becomes unresolved. One such scheme is the recently-proposed MiNLO method for the choice of scale and the exponentiation of Sudakov form factors in NLO computations. We present a validation study of our generators, by comparing their output to available LHC data.
  • We illustrate how the Matrix Element Method at Next-to-Leading Order (MEM@NLO) can be used to discriminate between events arising from the production of a Higgs boson, which subsequently decays to a final state consisting of ell^+ell^-gamma, and the background production of the same final state. We illustrate how the method could be used in an experimental analysis by devising cuts on the signal (P_S) and background (P_B) weights that are computed event-by-event in this approach. We find that we can increase the S/sqrt(B) ratio by around 50\% compared to an invariant mass fit on its own. Considering only statistical uncertainty, this is equivalent to recording a factor of around two times more integrated luminosity.
  • The success of the experimental program at the Tevatron re-inforced the idea that precision physics at hadron colliders is desirable and, indeed, possible. The Tevatron data strongly suggests that one-loop computations in QCD describe hard scattering well. Extrapolating this observation to the LHC, we conclude that knowledge of many short-distance processes at next-to-leading order may be required to describe the physics of hard scattering. While the field of one-loop computations is quite mature, parton multiplicities in hard LHC events are so high that traditional computational techniques become inefficient. Recently new approaches based on unitarity have been developed for calculating one-loop scattering amplitudes in quantum field theory. These methods are especially suitable for the description of multi-particle processes in QCD and are amenable to numerical implementations. We present a systematic pedagogical description of both conceptual and technical aspects of the new methods.
  • We present results for the production of a top pair in association with a W-boson at next-to-leading order. We have implemented this process into the parton-level integrator MCFM including the decays of both the top quarks and the W-bosons with full spin correlations. Although the cross section for this process is small, it is a Standard Model source of same-sign lepton events that must be accounted for in many new physics searches. For a particular analysis of same-sign lepton events in which b-quarks are also present, we investigate the effect of the NLO corrections as a function of the signal region cuts.
  • We describe the implementation of top production and decay processes in the parton-level Monte Carlo program MCFM. By treating the top quark as being on-shell, we can factorize the amplitudes for top-pair production, s-channel single-top production, and t-channel single-top production into the product of an amplitude for production and an amplitude for decay. In this way we can retain all spin correlations. Both the production and the decay amplitudes are calculated consistently at next-to-leading order in alpha_s. The full dependence on the b-quark mass is also kept. Phenomenological results are presented for various kinematic distributions at the LHC and for the top quark forward-backward asymmetry at the Tevatron.
  • We present a next-to-leading order calculation of Higgs boson production plus one and two jets via gluon fusion interfaced to shower Monte Carlo programs, implemented according to the POWHEG method. For this implementation we have used a new interface of the POWHEG BOX with MadGraph4, that generates the codes for generic Born and real processes automatically. The virtual corrections have been taken from the MCFM code. We carry out a simple phenomenological study of our generators, comparing them among each other and with fixed next-to-leading order results.