
Particle physics has an ambitious and broad experimental programme for the
coming decades. This programme requires large investments in detector hardware,
either to build new facilities and experiments, or to upgrade existing ones.
Similarly, it requires commensurate investment in the R&D of software to
acquire, manage, process, and analyse the shear amounts of data to be recorded.
In planning for the HLLHC in particular, it is critical that all of the
collaborating stakeholders agree on the software goals and priorities, and that
the efforts complement each other. In this spirit, this white paper describes
the R&D activities required to prepare for this software upgrade.

We present a systematic study of differences between NLL resummation and
parton showers. We first construct a Markovian MonteCarlo algorithm for
resummation of additive observables in electronpositron annihilation.
Approximations intrinsic to the pure NLL result are then removed, in order to
obtain a traditional, momentum and probability conserving parton shower based
on the coherent branching formalism. The impact of each approximation is
studied, and an overall comparison is made between the parton shower and pure
NLL resummation. Differences compared to modern partonshower algorithms
formulated in terms of color dipoles are analyzed.

The first computation of Higgs production in association with three jets at
NLO in QCD has recently been performed using the effective theory, where the
top quark is treated as an infinitely heavy particle and integrated out. This
approach is restricted to the regions in phase space where the typical scales
are not larger than the top quark mass. Here we investigate this statement at a
quantitative level by calculating the leadingorder contributions to the
production of a Standard Model Higgs boson in association with up to three jets
taking full topquark and bottomquark mass dependence into account. We find
that the transverse momentum of the hardest particle or jet plays a key role in
the breakdown of the effective theory predictions, and that discrepancies can
easily reach an order of magnitude for transverse momenta of about 1 TeV. The
impact of bottomquark loops are found to be visible in the small transverse
momentum region, leading to corrections of up to 5 percent. We further study
the impact of mass corrections when VBF selection cuts are applied and when the
centerofmass energy is increased to 100 TeV.

We present a detailed phenomenological analysis of the production of a
Standard Model Higgs boson in association with up to three jets. The Higgs is
produced via gluon fusion, which is an irreducible background to the vector
boson fusion mechanism. We calculate the nexttoleading order corrections in
QCD in the limit of an infinitely heavy top quark. Numerical results are
presented for a large variety of observables, for different selection cuts, and
for different choices of the jet tagging scheme.

Computing plays an essential role in all aspects of high energy physics. As
computational technology evolves rapidly in new directions, and data throughput
and volume continue to follow a steep trendline, it is important for the HEP
community to develop an effective response to a series of expected challenges.
In order to help shape the desired response, the HEP Forum for Computational
Excellence (HEPFCE) initiated a roadmap planning activity with two key
overlapping drivers  1) software effectiveness, and 2) infrastructure and
expertise advancement. The HEPFCE formed three working groups, 1) Applications
Software, 2) Software Libraries and Tools, and 3) Systems (including systems
software), to provide an overview of the current status of HEP computing and to
present findings and opportunities for the desired HEP computational roadmap.
The final versions of the reports are combined in this document, and are
presented along with introductory material.

In this talk, we discuss recent developments in combining parton showers and
fixedorder calculations. We focus on the UNNLOPS method for matching
nexttonexttoleading order computations to the parton shower, and we present
results from Sherpa for DrellYan leptonpair and Higgsboson production at the
LHC.

We present a detailed phenomenological analysis of the production of a
Standard Model Higgs boson in association with up to three jets. We consider
the gluon fusion channel using an effective theory in the large topquark mass
limit. Higgs boson production in gluon fusion constitutes an irreducible
background to the vector boson fusion (VBF) process; hence the precise
knowledge of its characteristics is a prerequisite for any measurement in the
VBF channel. The calculation is carried out at nexttoleading order (NLO) in
QCD in a fully automated way by combining the two programs GoSam and Sherpa. We
present numerical results for a large variety of observables for both standard
cuts and VBF selection cuts. We find that for all jet multiplicities the NLO
corrections are sizeable. This is particularly true in the presence of
kinematic selections enhancing the VBF topology, which are based on vetoing
additional jet activity. In this case, precise predictions for the background
can be made using our calculation by taking the difference between the
inclusive H+2 jets and the inclusive H+3 jets result.

We present a new framework for computing resummed and matched distributions
in processes with many hard QCD jets. The intricate color structure of soft
gluon emission at large angles renders resummed calculations highly nontrivial
in this case. We automate all ingredients necessary for the color evolution of
the soft function at nexttoleadinglogarithmic accuracy, namely the selection
of the color bases and the projections of color operators and Born amplitudes
onto those bases. Explicit results for all QCD processes with up to $2\to 5$
partons are given. We also devise a new treelevel matching scheme for resummed
calculations which exploits a quasilocal subtraction based on the
CataniSeymour dipole formalism. We implement both resummation and matching in
the Sherpa event generator. As a proof of concept, we compute the resummed and
matched transversethrust distribution for hadronic collisions.

We present a simple approach to combine NNLO QCD calculations and parton
showers, based on the UNLOPS technique. We apply the method to the computation
of DrellYan leptonpair production at the Large Hadron Collider. We comment on
possible improvements and intrinsic uncertainties.

In this talk we summarize some recent developments in perturbative QCD and
their application to particle physics phenomenology.

Uncertainties in the simulation of Higgs boson production with up to two jets
at nextto leading order accuracy are investigated. Traditional uncertainty
estimates based on scale variations are extended employing different functional
forms for the central scale, and the impact of details in the implementation of
the parton shower is discussed.

We present precise predictions for fourlepton plus jets production at the
LHC obtained within the fully automated Sherpa+OpenLoops framework. Offshell
intermediate vector bosons and related interferences are consistently included
using the complexmass scheme. Fourlepton plus 0 and 1jet final states are
described at NLO accuracy, and the precision of the simulation is further
increased by squared quarkloop NNLO contributions in the gg > 4l, gg > 4l+g,
gq > 4l+q, and qq > 4l+g channels. These NLO and NNLO contributions are
matched to the Sherpa parton shower, and the 0 and 1jet final states are
consistently merged using the MEPS@NLO technique. Thanks to Sudakov
resummation, the parton shower provides improved predictions and uncertainty
estimates for exclusive observables. This is important when jet vetoes or jet
bins are used to separate fourlepton final states arising from Higgs decays,
diboson production, and toppair production. Detailed predictions are presented
for the ATLAS and CMS H>WW analyses at 8 TeV in the 0 and 1jet bins.
Assessing renormalisation, factorisation and resummationscale uncertainties,
which reflect also unknown subleading Sudakov logarithms in jet bins, we find
that residual perturbative uncertainties are as small as a few percent.

In this talk the application of the recently introduced methods to merge NLO
calculations of successive jet multiplicities to the production of top pairs in
association with jets will be discussed, in particular a fresh look is taken at
the top quark forwardbackward asymmetries. Emphasis will be put on the
achieved theoretical accuracy and the associated perturbative and
nonperturbative error estimates.

The merging of matrix elements and parton showers is an established
calculational tool for the description of multijet final states at hadron
colliders. These methods have recently been promoted to nexttoleading order
accuracy in the description of hard well separated jets.
This talk introduces such a method and discusses its application to
phenomenologically relevant signal and background processes. The systematic
assessment of its theoretical uncertainty is a prime focus.

We present a study on highperformance computing and largescale distributed
computing for perturbative QCD calculations.

We present an analysis of the forwardbackward asymmetry in the production of
top quark pairs at the Tevatron collider. We use novel Monte Carlo methods for
merging matrix elements and parton showers to combine NLO QCD predictions for
tt and tt+jet production. Theoretical uncertainties are quantified in detail.
We find agreement with experimental data on the transverse momentum dependence
of the asymmetry.

We present a new approach to combine multiple NLO partonlevel calculations
matched to parton showers into a single inclusive event sample. The method
provides a description of hard multijet configurations at nextto leading
order in the perturbative expansion of QCD, and it is supplemented with the
allorders resummed modelling of jet fragmentation provided by the parton
shower. The formal accuracy of this technique is discussed in detail, invoking
the example of electronpositron annihilation into hadrons. We focus on the
effect of renormalisation scale variations in particular. Comparison with
experimental data from LEP underlines that this novel formalism describes data
with a theoretical accuracy that has hitherto not been achieved in standard
Monte Carlo event generators.

For the first time, differential cross sections for the production of
Wbosons in conjunction with up to three jets, computed at nextto leading
order in QCD and including parton shower corrections, are presented and
compared to recent experimental data from the Large Hadron Collider.

The MC@NLO and MEPS@NLO methods, as implemented in the MonteCarlo event
generator framework Sherpa, are used to estimate the perturbative and
nonperturbative uncertainties in various processes such as dijet production
and the production of a W boson in association with (multiple) jets.

In this contribution we review recent progress with fixedorder QCD
predictions for the production of a vector boson in association with jets at
hadron colliders, using the programs BlackHat and SHERPA. We review general
features of nexttoleadingorder (NLO) predictions for the production of a
massive vector boson in association with four jets. We also discuss how precise
descriptions of vectorboson production can be applied to the determination of
backgrounds to new physics signals. Here we focus on datadriven backgrounds to
a missingenergyplusjets search performed by CMS. Finally, we review recent
progress in developing theoretical tools for highmultiplicity loopcomputation
within the BlackHatlibrary. In particular, we discuss methods for handling the
color degrees of freedom in multijet predictions at NLO.

In this publication, uncertainties in and differences between the MC@NLO and
POWHEG methods for matching nexttoleading order QCD calculations with parton
showers are discussed. Implementations of both algorithms within the event
generator Sherpa and based on CataniSeymour subtraction are employed to assess
the impact on a representative selection of observables. In the case of MC@NLO
a substantial simplification is achieved by using dipole subtraction terms to
generate the first emission. A phase space restriction is employed, which
allows to vary in a transparent way the amount of nonsingular radiative
corrections that are exponentiated. Effects on various observables are
investigated, using the production of a Higgs boson in gluon fusion, with or
without an associated jet, as a benchmark process. The case of H+jet production
is presented for the first time in an NLO+PS matched simulation. Uncertainties
due to scale choices and nonperturbative effects are explored in the
production of W and Z bosons in association with a jet. Corresponding results
are compared to data from the Tevatron and LHC experiments.

This is a written account of the computer tutorial offered at the Sixth
MC4BSM workshop at Cornell University, March 2224, 2012. The tools covered
during the tutorial include: FeynRules, LanHEP, MadGraph, CalcHEP, Pythia 8,
Herwig++, and Sherpa. In the tutorial, we specify a simple extension of the
Standard Model, at the level of a Lagrangian. The software tools are then used
to automatically generate a set of Feynman rules, compute the invariant matrix
element for a sample process, and generate both partonlevel and fully
hadronized/showered Monte Carlo event samples. The tutorial is designed to be
selfpaced, and detailed instructions for all steps are included in this
writeup. Installation instructions for each tool on a variety of popular
platforms are also provided.

We quantify uncertainties in the MonteCarlo simulation of inclusive and
dijet final states, which arise from using the MC@NLO technique for matching
nexttoleading order parton level calculations and parton showers. We analyse
a large variety of data from early measurements at the LHC. In regions of phase
space where Sudakov logarithms dominate over highenergy effects, we observe
that the main uncertainty can be ascribed to the free parameters of the parton
shower. In complementary regions, the main uncertainty stems from the
considerable freedom in the simulation of underlying events.

We present a processindependent technique to consistently combine
nexttoleading order partonlevel calculations of varying jet multiplicity and
parton showers. Double counting is avoided by means of a modified truncated
shower scheme. This method preserves both the fixedorder accuracy of the
partonlevel result and the logarithmic accuracy of the parton shower. We
discuss the renormalisation and factorisation scale dependence of the approach
and present results from an automated implementation in the Sherpa event
generator using the test case of Wboson production at the Large Hadron
Collider. We observe a dramatic reduction of theoretical uncertainties compared
to existing methods which underlines the predictive power of our novel
technique.

The pMSSM provides a broad perspective on SUSY phenomenology. In this paper
we generate two new, very large, sets of pMSSM models with sparticle masses
extending up to 4 TeV, where the lightest supersymmetric particle (LSP) is
either a neutralino or gravitino. The existence of a gravitino LSP necessitates
a detailed study of its cosmological effects and we find that Big Bang
Nucleosynthesis places strong constraints on this scenario. Both sets are
subjected to a global set of theoretical, observational and experimental
constraints resulting in a sample of \sim 225k viable models for each LSP type.
The characteristics of these two model sets are briefly compared. We confront
the neutralino LSP model set with searches for SUSY at the 7 TeV LHC using both
the missing (MET) and nonmissing ET ATLAS analyses. In the MET case, we employ
Monte Carlo estimates of the ratios of the SM backgrounds at 7 and 8 TeV to
rescale the 7 TeV datadriven ATLAS backgrounds to 8 TeV. This allows us to
determine the pMSSM parameter space coverage for this collision energy. We find
that an integrated luminosity of \sim 520 fb^{1} at 8 TeV would yield a
substantial increase in this coverage compared to that at 7 TeV and can probe
roughly half of the model set. If the pMSSM is not discovered during the 8 TeV
run, then our model set will be essentially void of gluinos and lightest first
and second generation squarks that are \lesssim 700800 GeV, which is much less
than the analogous mSUGRA bound. Finally, we demonstrate that nonMET SUSY
searches continue to play an important role in exploring the pMSSM parameter
space. These two pMSSM model sets can be used as the basis for investigations
for years to come.