
Predictions for cold nuclear matter effects on charged hadrons, identified
light hadrons, quarkonium and heavy flavor hadrons, DrellYan dileptons, jets,
photons, gauge bosons and top quarks produced in $p+$Pb collisions at
$\sqrt{s_{_{NN}}} = 8.16$ TeV are compiled and, where possible, compared to
each other. Predictions of the normalized ratios of $p+$Pb to $p+p$ cross
sections are also presented for most of the observables, providing new insights
into the expected role of cold nuclear matter effects. In particular, the role
of nuclear parton distribution functions on particle production can now be
probed over a wider range of phase space than ever before.

We present numerical tests and predictions of the KrkNLO method for matching
of NLO QCD corrections to hard processes with LO parton shower Monte Carlo
generators (NLO+PS). This method was described in detail in our previous
publications, where it was also compared with other NLO+PS matching approaches
(MC@NLO and POWHEG) as well as fixedorder NLO and NNLO calculations. Here we
concentrate on presenting some numerical results (cross sections and
distributions) for $Z/\gamma^*$ (DrellYan) and Higgsboson production
processes at the LHC. The DrellYan process is used mainly to validate the
KrkNLO implementation in the Herwig 7 program with respect to the previous
implementation in Sherpa. We also show predictions for this process with the
new, complete, MCscheme parton distribution functions and compare them with
our previously published results. Then, we present the first results of the
KrkNLO method for Higgs production in gluongluon fusion at the LHC and compare
them with MC@NLO and POWHEG predictions from Herwig 7 fixedorder results from
HNNLO and a resummed calculation from HqT, as well as with experimental data
from the ATLAS collaboration.

Using the framework that interpolates between the leading power limit of the
Color Glass Condensate and the High Energy (or $k_{T}$) factorization we
calculate the direct component of the forward dijet production in
ultraperipheral $\mathrm{Pb}$$\mathrm{Pb}$ collisions at CM energy
$5.1\,\mathrm{TeV}$ per nucleon pair. The formalism is applicable when the
average transverse momentum of the dijet system $P_{T}$ is much bigger than the
saturation scale $Q_{s}$, $P_{T}\gg Q_{s}$, while the imbalance of the dijet
system can be arbitrary. The cross section is uniquely sensitive to the
Weizs\"ackerWilliams (WW) unintegrated gluon distribution, which is far less
known from experimental data than the most common dipole gluon distribution
appearing in inclusive small$x$ processes. We have calculated cross sections
and nuclear modification ratios using WW gluon distribution obtained from the
dipole gluon density through the Gaussian approximation. The dipole gluon
distribution used to get WW was fitted to the inclusive HERA data with the
nonlinear extension of unified BFKL+DGLAP evolution equation. The saturation
effects are visible but rather weak for realistic $p_{T}$ cut on the dijet
system, reaching about $20\%$ with the cut as low as $6\,\mathrm{GeV}$. We find
that the LO collinear factorization with nuclear leading twist shadowing
predicts quite similar effects.

The Future Circular Collider (FCC) Study is aimed at assessing the physics
potential and the technical feasibility of a new collider with centreofmass
energies, in the hadronhadron collision mode, seven times larger than the
nominal LHC energies. Operating such machine with heavy ions is an option that
is being considered in the accelerator design studies. It would provide, for
example, PbPb and pPb collisions at sqrt{s_NN} = 39 and 63 TeV, respectively,
per nucleonnucleon collision, with integrated luminosities above 30 nb^1 per
month for PbPb. This is a report by the working group on heavyion physics of
the FCC Study. First ideas on the physics opportunities with heavy ions at the
FCC are presented, covering the physics of the QuarkGluon Plasma, of gluon
saturation, of photoninduced collisions, as well as connections with other
fields of highenergy physics.

Next steps in development of the KrkNLO method of implementing NLO QCD
corrections to hard processes in parton shower Monte Carlo programs are
presented. This new method is a simpler alternative to other wellknown
approaches, such as MC@NLO and POWHEG. The KrkNLO method owns its simplicity to
the use of parton distribution functions (PDFs) in a new, socalled Monte Carlo
(MC), factorization scheme which was recently fully defined for the first time.
Preliminary numerical results for the Higgsboson production process are also
presented.

A next step in development of the KrkNLO method of including complete NLO QCD
corrections to hard processes in a LO partonshower Monte Carlo (PSMC) is
presented. It consists of generalisation of the method, previously used for the
DrellYan process, to Higgsboson production. This extension is accompanied
with the complete description of parton distribution functions (PDFs) in a
dedicated, Monte Carlo (MC) factorisation scheme, applicable to any process of
production of one or more colourneutral particles in hadronhadron collisions.

We study the production of forward dijets in protonlead and protonproton
collisions at the Large Hadron Collider. Such configurations, with both jets
produced in the forward direction, impose a dilutedense asymmetry which allows
to probe the gluon density of the lead or proton target at small longitudinal
momentum fractions. Even though the jet momenta are always much bigger than the
saturation scale of the target, $Q_s$, the transverse momentum imbalance of the
dijet system may be either also much larger than $Q_s$, or of the order $Q_s$,
implying that the small$x$ QCD dynamics involved is either linear or
nonlinear, respectively. The small$x$ improved TMD factorization framework
deals with both situation in the same formalism. In the latter case, which
corresponds to nearly backtoback jets, we find that saturation effects induce
a significant suppression of the forward dijet azimuthal correlations in
protonlead versus protonproton collisions.

Predictions made in Albacete {\it et al} prior to the LHC $p+$Pb run at
$\sqrt{s_{NN}} = 5$ TeV are compared to currently available data. Some
predictions shown here have been updated by including the same experimental
cuts as the data. Some additional predictions are also presented, especially
for quarkonia, that were provided to the experiments before the data were made
public but were too late for the original publication are also shown here.

A new method of including the complete NLO QCD corrections to hard processes
in the LO partonshower Monte Carlo (PSMC) is presented. This method, called
KrkNLO, requires the use of parton distribution functions in a dedicated Monte
Carlo factorization scheme, which is also discussed in this paper. In the
future, it may simplify introduction of the NNLO corrections to hard processes
and the NLO corrections to PSMC. Details of the method and numerical examples
of its practical implementation, as well as comparisons with other
calculations, such as MCFM, MC@NLO, POWHEG, for single $Z/\gamma^*$boson
production at the LHC, are presented.

In these proceedings, we review results for several diboson production
processes beyond NLO QCD at high transverse momenta using the VBFNLO
MonteCarlo program together with the LOOPSIM method. Additionally, we show for
the WZ production process how higher order QCD corrections can resemble
anomalous coupling effects.

We study forward dijet production in dilutedense hadronic collisions. By
considering the appropriate limits, we show that both the
transversemomentumdependent (TMD) and the highenergy factorization formulas
can be derived from the Color Glass Condensate framework. Respectively, this
happens when the transverse momentum imbalance of the dijet system, $k_t$, is
of the order of either the saturation scale, or the hard jet momenta, the
former being always much smaller than the latter. We propose a new formula for
forward dijets that encompasses both situations and is therefore applicable
regardless of the magnitude of $k_t$. That involves generalizing the TMD
factorization formula for dijet production to the case where the incoming
small$x$ gluon is offshell. The derivation is performed in two independent
ways, using either Feynman diagram techniques, or colorordered amplitudes.

We study saturation effects in the production of forward dijets in
protonlead collisions at the Large Hadron Collider, using the framework of
High Energy Factorization. Such configurations, with both jets produced in the
forward direction, probe the gluon density of the lead nucleus at small
longitudinal momentum fraction, and also limit the phase space for emissions of
additional jets. We find significant suppression of the forward dijet azimuthal
correlations in protonlead versus protonproton collisions, which we attribute
to stronger saturation of the gluon density in the nucleus than in the proton.
In order to minimize model dependence of our predictions, we use two different
extensions of the BalitskyKovchegov equation for evolution of the gluon
density with subleading corrections.

We provide a description, within the High Energy Factorization formalism, of
centralforward dijet decorrelation data measured by the CMS experiment and the
predictions for nuclear modification ratio~$R_{pA}$ in p+Pb collisions. In our
study, we use the unintegrated gluon density derived from the BFKL and BK
equations supplemented with subleading corrections and a hard scale dependence.
The latter is introduced at the final step of the calculation by reweighting
the Monte Carlo generated events using suitable Sudakov form factors, without
changing the total cross section. We achieve a good description of data in the
whole region of the azimuthal angle.

This is the summary report of the energy frontier QCD working group prepared
for Snowmass 2013. We review the status of tools, both theoretical and
experimental, for understanding the strong interactions at colliders. We
attempt to prioritize important directions that future developments should
take. Most of the efforts of the QCD working group concentrate on protonproton
colliders, at 14 TeV as planned for the next run of the LHC, and for 33 and 100
TeV, possible energies of the colliders that will be necessary to carry on the
physics program started at 14 TeV. We also examine QCD predictions and
measurements at leptonlepton and leptonhadron colliders, and in particular
their ability to improve our knowledge of strong coupling constant and parton
distribution functions.

Predictions for charged hadron, identified light hadron, quarkonium, photon,
jet and gauge bosons in p+Pb collisions at sqrt s_NN = 5 TeV are compiled and
compared. When test run data are available, they are compared to the model
predictions.

The 2011 Les Houches workshop was the first to confront LHC data. In the two
years since the previous workshop there have been significant advances in both
soft and hard QCD, particularly in the areas of multileg NLO calculations, the
inclusion of those NLO calculations into parton shower Monte Carlos, and the
tuning of the nonperturbative parameters of those Monte Carlos. These
proceedings describe the theoretical advances that have taken place, the impact
of the early LHC data, and the areas for future development.

The workshop on "HadronHadron and CosmicRay Interactions at multiTeV
Energies" held at the ECT* centre (Trento) in Nov.Dec. 2010 gathered together
both theorists and experimentalists to discuss issues of the physics of
highenergy hadronic interactions of common interest for the particle, nuclear
and cosmicray communities. QCD results from collider experiments  mostly
from the LHC but also from the Tevatron, RHIC and HERA  were discussed and
compared to various hadronic Monte Carlo generators, aiming at an improvement
of our theoretical understanding of soft, semihard and hard parton dynamics.
The latest cosmicray results from various groundbased observatories were also
presented with an emphasis on the phenomenological modeling of the first
hadronic interactions of the extended airshowers generated in the Earth
atmosphere. These miniproceedings consist of an introduction and short
summaries of the talks presented at the meeting.

We present the miniproceedings of the workshop on ``Parton fragmentation in
the vacuum and in the medium'' held at the European Centre for Theoretical
Studies in Nuclear Physics and Related Areas (ECT*, Trento) in February 2008.
The workshop gathered both theorists and experimentalists to discuss the
current status of investigations of quark and gluon fragmentation into hadrons
at different accelerator facilities (LEP, Bfactories, JLab, HERA, RHIC, and
Tevatron) as well as preparations for extension of these studies at the LHC.
The main physics topics covered were: (i) lightquark and gluon fragmentation
in the vacuum including theoretical (global fits analyses and MLLA) and
experimental (data from e+e, pp, ep collisions) aspects, (ii) strange and
heavyquark fragmentation, (iii) parton fragmentation in cold QCD matter
(nuclear DIS), and (iv) mediummodified fragmentation in hot and dense QCD
matter (highenergy nucleusnucleus collisions). These miniproceedings consist
of an introduction and short summaries of the talks presented at the meeting.

This writeup is a compilation of the predictions for the forthcoming Heavy
Ion Program at the Large Hadron Collider, as presented at the CERN Theory
Institute 'Heavy Ion Collisions at the LHC  Last Call for Predictions', held
from May 14th to June 10th 2007.

The saturation model with DGLAP evolution is shown to give good description
of the production of the charm and beauty quarks in deep inelastic scattering.
The modifications of saturation properties caused by the presence of heavy
quarks are also discussed.

The most recent H1 and ZEUS data for diffractive structure functions are
analyzed under three different theoretical approaches. This includes the
Pomeron Structure Function (PSF) framework, BartelsEllisKowalskiWusthoff
(BEKW) color dipole approach and the GolecBiernatWusthoff (GBW) saturation
model. The models are shown to successfully fit the set of combined data.
Conceptual differences between the models are discussed and, as an example, the
prediction for the longitudinal diffractive structure function is presented.

We derive the asymptotic travelingwave solutions of the nonlinear
1dimensional BalitskyKovchegov QCD equation for rapidity evolution in
momentumspace, with 1loop running coupling constant and equipped with the
BalitskyKovchegovKuraevLipatov kernel at nexttoleading logarithmic
accuracy, conveniently regularized by different resummation schemes. Traveling
waves allow to define "universality classes" of asymptotic solutions, i.e.
independent of initial conditions and of the nonlinear damping. A dependence on
the resummation scheme remains, which is analyzed in terms of geometric scaling
properties.

We describe the most recent data on the diffractive structure functions from
the H1 and ZEUS Collaborations at HERA using four models. First, a Pomeron
Structure Function (PSF) model, in which the Pomeron is considered as an object
with parton distribution functions. Then, the Bartels Ellis Kowalski Wusthoff
(BEKW) approach is discussed, assuming the simplest perturbative description of
the Pomeron using a twogluon ladder. A third approach, the Bialas Peschanski
(BP) model, based on the dipole formalism is then described. Finally, we
discuss the GolecBiernatW\"usthoff (GBW) saturation model which takes into
account saturation effects. The best description of all avaible measurements
can be achieved with either the PSF based model or the BEKW approach. In
particular, the BEKW prediction allows to include the highest $\beta$
measurements, which are dominated by higher twists effects and provide an
efficient and compact parametrisation of the diffractive cross section. The two
other models also give a good description of cross section measurements at
small $x$ with a small number of parameters. The comparison of all predictions
allows us to identify interesting differences in the behaviour of the effective
pomeron intercept and in the shape of the longitudinal component of the
diffractive structure functions. In this last part, we present some features
that can be discriminated by new experimental measurements, completing the HERA
program.

Predictions for the total, elastic and single diffractive cross sections
calculated for the LHC in the framework of the MiettinenPumplin model are
presented. The total cross section is expected to be 15% smaller than that
determined by Donnachie and Landshoff in the model with soft pomeron. The
diffractive cross section is almost constant in the TevatronLHC energy range.

We discuss polarized leptonproton scattering with special emphasis on the
difference between target polarization defined relative to the lepton beam or
to the virtual photon direction. In particular, this difference influences
azimuthal distributions in the final state. We provide a general framework of
analysis and apply it to the specific cases of semiinclusive deep inelastic
scattering, of exclusive meson production, and of deeply virtual Compton
scattering.