
We calculate the soft function for the global event variable 1jettiness at
nexttonexttoleading order (NNLO) in QCD. We focus specifically on the
nonAbelian contribution, which, unlike the Abelian part, is not determined by
the nexttoleading order result. The calculation uses the known general forms
for the emission of one and two soft partons and is performed using a
sectordecomposition method that is spelled out in detail. Results are
presented in the form of numerical fits to the 1jettiness soft function for
LHC kinematics (as a function of the angle between the incoming beams and the
finalstate 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 Njettiness event variable to handle infrared
singularities.

In this paper we present a calculation of the $\gamma+j$ process at
nexttonexttoleading 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
nexttonextto 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 reevaluation of any approximations inherent in the
theoretical calculations. In this paper we address one such issue, the
inclusion of mass effects in topquark loops. We ameliorate an existing
calculation of $Z+1$~jet and $Z+2$~jet production by presenting exact analytic
formulae for amplitudes containing topquark loops that enter at
nexttoleading order in QCD. Although approximations based on an expansion in
powers of $1/m_t^2$ can lead to poor highenergy behavior, an exact treatment
of topquark 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
objectoriented 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, 19731985.

We present results for the production of a pair of onshell 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 twoloop 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 NexttoNextto 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 partonlevel 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 offshell 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 finalstate phase space by studying
$pp\rightarrow V(H\rightarrow WW^*) \rightarrow$ leptons.

In this paper we present a nexttonexttoleading 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
gluoninitiated 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 datadriven fits.

We present the first complete calculation of Zboson production in
association with a jet in hadronic collisions through nexttonexttoleading
order in perturbative QCD. Our computation uses the recentlyproposed
Njettiness subtraction scheme to regulate the infrared divergences that appear
in the realemission contributions. We present phenomenological results for 13
TeV protonproton collisions with fully realistic fiducial cuts on the
finalstate particles. The remaining theoretical uncertainties after the
inclusion of our calculations are at the percentlevel, making the Z+jet
channel ready for precision studies at the LHC Run II.

We present the implementation of several colorsinglet finalstate processes
at NexttoNextto Leading Order (NNLO) accuracy in QCD to the publicly
available partonlevel 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 nonlocal $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
multiprocessor systems.

We present a nexttoleading 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
multithreading. 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 4lepton + 2jet 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 highmass behaviour of the
Higgs boson mediated diagrams, which unlike onshell production, gives
information about the Higgs couplings without assumptions on the Higgs boson
total width. We assess the sensitivity of the highmass 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 CaolaMelnikov (CM) method,
the offshell 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 bosonmediated diagrams via the decay $H
\to ZZ$ (signal diagrams), as well as from diagrams where the Zbosons 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
offresonant effect for the Higgs boson, we treat the Zbosons as onshell.
Thus our analysis is limited to the region above threshold, where the invariant
mass of the Zpair, $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 jetbinned data
become available.

We present results for the Standard model description of the fourlepton
production, mediated both by Higgs boson production and by other oneloop
standard model processes. The description of fourlepton 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 offshell $H\rightarrow ZZ$
events have shown that Run 1 LHC data can constrain the Higgs width, $\Gamma_H
< (2545) \Gamma_{H}^{\rm SM}$. Using 20 fb1 of 8 TeV ATLAS data, we estimate
a bound on the Higgs boson width from the WW channel between $\Gamma_H <
(100500) \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 onshell Higgs boson limit the potential number of offshell
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 $(2550) \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 fourlepton 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 offshell 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 Higgsmediated and continuum production of fourlepton
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 offshell region to constrain the total width of the Higgs boson. Using a
simple cutandcount method, existing LHC data should bound the width at the
level of 2545 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 nexttoleading order (NLO), including the decay of the top quark and
the Z boson. This electroweak process gives rise to the trilepton signature
l+ll'+ + 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 btag 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 partonlevel 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 multiparton
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 recentlyproposed 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 NexttoLeading 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 eventbyevent 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 reinforced the idea
that precision physics at hadron colliders is desirable and, indeed, possible.
The Tevatron data strongly suggests that oneloop computations in QCD describe
hard scattering well. Extrapolating this observation to the LHC, we conclude
that knowledge of many shortdistance processes at nexttoleading order may be
required to describe the physics of hard scattering. While the field of
oneloop 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
oneloop scattering amplitudes in quantum field theory. These methods are
especially suitable for the description of multiparticle 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
Wboson at nexttoleading order. We have implemented this process into the
partonlevel integrator MCFM including the decays of both the top quarks and
the Wbosons with full spin correlations. Although the cross section for this
process is small, it is a Standard Model source of samesign lepton events that
must be accounted for in many new physics searches. For a particular analysis
of samesign lepton events in which bquarks 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
partonlevel Monte Carlo program MCFM. By treating the top quark as being
onshell, we can factorize the amplitudes for toppair production, schannel
singletop production, and tchannel singletop 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 nexttoleading order in alpha_s. The full
dependence on the bquark mass is also kept. Phenomenological results are
presented for various kinematic distributions at the LHC and for the top quark
forwardbackward asymmetry at the Tevatron.

We present a nexttoleading 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 nexttoleading
order results.