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 HL-LHC 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 Monte-Carlo algorithm for resummation of additive observables in electron-positron 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 parton-shower 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 leading-order contributions to the production of a Standard Model Higgs boson in association with up to three jets taking full top-quark and bottom-quark 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 bottom-quark 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 center-of-mass 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 next-to-leading 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 trend-line, 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 (HEP-FCE) initiated a roadmap planning activity with two key overlapping drivers -- 1) software effectiveness, and 2) infrastructure and expertise advancement. The HEP-FCE 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 fixed-order calculations. We focus on the UNNLOPS method for matching next-to-next-to-leading order computations to the parton shower, and we present results from Sherpa for Drell-Yan lepton-pair and Higgs-boson 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 top-quark 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 next-to-leading 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 non-trivial in this case. We automate all ingredients necessary for the color evolution of the soft function at next-to-leading-logarithmic 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 tree-level matching scheme for resummed calculations which exploits a quasi-local subtraction based on the Catani-Seymour dipole formalism. We implement both resummation and matching in the Sherpa event generator. As a proof of concept, we compute the resummed and matched transverse-thrust 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 Drell-Yan lepton-pair 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 next-to 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 four-lepton plus jets production at the LHC obtained within the fully automated Sherpa+OpenLoops framework. Off-shell intermediate vector bosons and related interferences are consistently included using the complex-mass scheme. Four-lepton plus 0- and 1-jet final states are described at NLO accuracy, and the precision of the simulation is further increased by squared quark-loop 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 1-jet 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 four-lepton final states arising from Higgs decays, diboson production, and top-pair production. Detailed predictions are presented for the ATLAS and CMS H->WW analyses at 8 TeV in the 0- and 1-jet bins. Assessing renormalisation-, factorisation- and resummation-scale 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 forward-backward asymmetries. Emphasis will be put on the achieved theoretical accuracy and the associated perturbative and non-perturbative error estimates.

The merging of matrix elements and parton showers is an established calculational tool for the description of multi-jet final states at hadron colliders. These methods have recently been promoted to next-to-leading 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 high-performance computing and large-scale distributed computing for perturbative QCD calculations.

We present an analysis of the forward-backward 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 parton-level calculations matched to parton showers into a single inclusive event sample. The method provides a description of hard multi-jet configurations at next-to leading order in the perturbative expansion of QCD, and it is supplemented with the all-orders resummed modelling of jet fragmentation provided by the parton shower. The formal accuracy of this technique is discussed in detail, invoking the example of electron-positron 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 W-bosons in conjunction with up to three jets, computed at next-to 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 Monte-Carlo event generator framework Sherpa, are used to estimate the perturbative and non-perturbative 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 fixed-order 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 next-to-leading-order (NLO) predictions for the production of a massive vector boson in association with four jets. We also discuss how precise descriptions of vector-boson production can be applied to the determination of backgrounds to new physics signals. Here we focus on data-driven backgrounds to a missing-energy-plus-jets search performed by CMS. Finally, we review recent progress in developing theoretical tools for high-multiplicity loop-computation within the BlackHat-library. In particular, we discuss methods for handling the color degrees of freedom in multi-jet predictions at NLO.

In this publication, uncertainties in and differences between the MC@NLO and POWHEG methods for matching next-to-leading order QCD calculations with parton showers are discussed. Implementations of both algorithms within the event generator Sherpa and based on Catani-Seymour 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 non-singular 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 non-perturbative 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 22-24, 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 parton-level and fully hadronized/showered Monte Carlo event samples. The tutorial is designed to be self-paced, and detailed instructions for all steps are included in this write-up. Installation instructions for each tool on a variety of popular platforms are also provided.

We quantify uncertainties in the Monte-Carlo simulation of inclusive and dijet final states, which arise from using the MC@NLO technique for matching next-to-leading 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 high-energy 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 process-independent technique to consistently combine next-to-leading order parton-level 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 fixed-order accuracy of the parton-level 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 W-boson 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 non-missing 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 data-driven 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 5-20 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 700-800 GeV, which is much less than the analogous mSUGRA bound. Finally, we demonstrate that non-MET 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.