
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.

This document collects the proceedings of the "Parton Radiation and
Fragmentation from LHC to FCCee" workshop (http://indico.cern.ch/e/ee\_jets16)
held at CERN in Nov. 2016. The writeup reviews the latest theoretical and
experimental developments on parton radiation and partonhadron fragmentation
studies including analyses of LEP, Bfactories, and LHC data with a focus
on the future perspectives reacheable in $e^+e^$ measurements at the Future
Circular Collider (FCCee), with multiab$^{1}$ integrated luminosities
yielding 10$^{12}$ and 10$^{8}$ jets from Z and W bosons decays as well as
10$^5$ gluon jets from Higgs boson decays. The main topics discussed are: (i)
parton radiation and partontohadron fragmentation functions (splitting
functions at NNLO, small$z$ NNLL resummations, global FF fits including Monte
Carlo (MC) and neuralnetwork analyses of the latest Belle/BaBar highprecision
data, parton shower MC generators), (ii) jet properties (quarkgluon
discrimination, $e^+e^$ event shapes and multijet rates at NNLO+N$^{n}$LL,
jet broadening and angularities, jet substructure at smallradius, jet charge
determination, $e^+e^$ jet reconstruction algorithms), (iii) heavyquark jets
(dead cone effect, charmbottom separation, gluonto$b\bar{b}$ splitting), and
(iv) nonperturbative QCD phenomena (colour reconnection, baryon and
strangeness production, BoseEinstein and FermiDirac finalstate correlations,
colour string dynamics: spin effects, helix hadronization).

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.

This report was prepared in the context of the LPCC "Electroweak Precision
Measurements at the LHC WG" and summarizes the activity of a subgroup dedicated
to the systematic comparison of public Monte Carlo codes, which describe the
DrellYan processes at hadron colliders, in particular at the CERN Large Hadron
Collider (LHC). This work represents an important step towards the definition
of an accurate simulation framework necessary for very highprecision
measurements of electroweak (EW) observables such as the $W$ boson mass and the
weak mixing angle. All the codes considered in this report share at least
nexttoleadingorder (NLO) accuracy in the prediction of the total cross
sections in an expansion either in the strong or in the EW coupling constant.
The NLO fixedorder predictions have been scrutinized at the technical level,
using exactly the same inputs, setup and perturbative accuracy, in order to
quantify the level of agreement of different implementations of the same
calculation. A dedicated comparison, again at the technical level, of three
codes that reach nexttonexttoleadingorder (NNLO) accuracy in quantum
chromodynamics (QCD) for the total cross section has also been performed. These
fixedorder results are a welldefined reference that allows a classification
of the impact of higherorder sets of radiative corrections. Several examples
of higherorder effects due to the strong or the EW interaction are discussed
in this common framework. Also the combination of QCD and EW corrections is
discussed, together with the ambiguities that affect the final result, due to
the choice of a specific combination recipe.

We show that already at the NLO level the DGLAP evolution kernel Pqq starts
to depend on the choice of the evolution variable. We give an explicit example
of such a variable, namely the maximum of transverse momenta of emitted partons
and we identify a class of evolution variables that leave the NLO Pqq kernel
unchanged with respect to the known standard MSbar results. The kernels are
calculated using a modified CurciFurmanskiPetronzio method which is based on
a direct Feynmangraphs calculation.

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.

The decadeold technique of combining NLOcorrected hard process with
LOlevel parton shower Monte Carlo is now mature and used in practice of the
QCD calculations in the LHC data analysis. The next step, its extension to an
NNLOcorrected hard process combined with the NLOlevel parton shower Monte
Carlo, will require development of the latter component. It does not exist yet
in a complete form. In this note we describe recent progress in developing the
NLO parton shower for the initialstate hadron beams. The technique of adding
NLO corrections in the fully exclusive form (defined in recent years) is now
simplified and tested numerically, albeit for a limited set of NLO diagrams in
the evolution kernels.

We present the Monte Carlo event generator WINHAC for DrellYan processes in
protonproton, protonantiproton, protonion and ionion collisions. It
features multiphoton radiation within the YennieFrautschiSuura exclusive
exponentiation scheme with O(alpha) electroweak corrections for the
chargedcurrent (W+/W) processes and multiphoton radiation generated by PHOTOS
for neutralcurrent (Z+gamma) ones. For the initialstate QCD/QED parton shower
and hadronisation it is interfaced with PYTHIA. It includes several options,
e.g. for the polarized Wboson production, generation of weighted/unweighted
events, etc. WINHAC was crosschecked numerically at the permille level with
independent Monte Carlo programs, such as HORACE and SANC. It has been used as
a basic tool for developing and testing some new methods of precise
measurements of the Standard Model parameters at the LHC, in particular the
Wboson mass. Recently, it has been applied to simulations of double DrellYan
processes resulting from doubleparton scattering, in order to assess their
influence on the Higgsboson detection at the LHC in its ZZ and W+W decay
channels.

The concept of new methodology of adding QCD NLO corrections in the initial
state Monte Carlo parton shower (hard process part) is tested numerically
using, as an example, the process of the heavy boson production at
hadronhadron colliders such as LHC. In spite of the use of a simplified model
of the process, all presented numerical results prove convincingly that the
basic concept of the new methodology works correctly in practice, that is in
the numerical environment of the Monte Carlo parton shower event generator. The
differences with the other well established methods, like MC@NLO and POWHEG,
are briefly discussed and future refinements of the implementation of the new
method are also outlined.

Charged lepton transverse momenta in the DrellYan processes play an
important role at the LHC in precision measurements of the Standard Model
parameters, such as the Wboson mass and width, their charge asymmetries and
sin^2(theta_W). Therefore, their distributions should be described as accurate
as possible by the Monte Carlo event generators. In this paper we discuss the
problem of matching the hardprocess kinematics of the Monte Carlo generator
WINHAC with the partonshower kinematics of the PYTHIA 6.4 generator while
interfacing these two programs. We show that improper assignment of the quark
and antiquark effective momenta in the LO matrix element computations may
affect considerably the predicted lepton transverse momenta and even completely
reverse their charge asymmetries at the LHC. We propose two matching schemes in
which the NLO QCD distributions of the leptonic kinematical variables can be
well reproduced by the LO WINHAC generator.

The project of constructing a complete NLOlevel Parton Shower Monte Carlo
for the QCD processes developed in IFJ PAN in Krakow is reviewed. Four issues
are discussed: (1) the extension of the standard inclusive collinear
factorization into a new, fully exclusive scheme; (2) reconstruction of the LO
Parton Shower in the new scheme; (3) inclusion of the exclusive NLO corrections
into the hard process and (4) inclusion of the exclusive NLO corrections into
the evolution (ladder) part.

At the LHC, the measurement of the W mass with a precision of \cal{O}(10)
MeV/c^2 is both mandatory and difficult. In the analysis strategies proposed so
far, shortcuts have been made that are justified for protonantiproton
collisions at the Tevatron, but not for protonproton collisions at the LHC.
The root of the problem lies in the inadequate knowledge of parton density
functions of the proton. It is argued that in order to reach a 10 MeV/c^2
precision for the W mass, more precise parton density functions of the proton
are needed, and an LHCspecific analysis strategy ought to be pursued.
Proposals are made on both issues.

This paper is the second of the series of papers proposing dedicated
strategies for precision measurements of the Standard Model parameters at the
LHC. The common feature of these strategies is their robustness with respect to
the systematic measurement and modeling error sources. Their impact on the
precision of the measured parameters is reduced using dedicated observables and
dedicated measurement procedures which exploit flexibilities of the collider
and detector running modes. In the present paper we focus our attention on the
measurement of the charge asymmetry of the Wboson mass. This measurement is of
primordial importance for the LHC experimental program, both as a direct test
of the chargesignindependent coupling of the Wbosons to the matter particles
and as a necessary first step towards the precision measurement of the
chargeaveraged Wboson mass. We propose and evaluate the LHCspecific strategy
to measure the mass difference between the positively and negatively charged
Wbosons, MW+  MW. We show that its present precision can be improved at the
LHC by a factor of 20. We argue that such a precision is beyond the reach of
the standard measurement and calibration methods imported to the LHC from the
Tevatron program.

This work covers methodology of solving QCD evolution equation of the parton
distribution using Markovian Monte Carlo (MMC) algorithms in a class of models
ranging from DGLAP to CCFM. One of the purposes of the above MMCs is to test
the other more sophisticated Monte Carlo programs, the socalled Constrained
Monte Carlo (CMC) programs, which will be used as a building block in the
parton shower MC. This is why the mapping of the evolution variables (eikonal
variable and evolution time) into fourmomenta is also defined and tested. The
evolution time is identified with the rapidity variable of the emitted parton.
The presented MMCs are tested independently, with ~0.1% precision, against the
nonMC program APCheb especially devised for this purpose.

In this paper we describe a check of the implementation of SANC system
generated modules into the framework of the WINHAC Monte Carlo event generator.
At this stage of work we limit ourselves to inclusion of complete oneloop
electroweak corrections to the chargedcurrent DrellYan process. We perform
tuned comparisons of the results derived with the aid of two codes: 1) the
standard SANC integrator with YFSinspired treatment of the ISR QED corrections
and 2) the WINHAC generator, upgraded with the SANC electroweak modules and
downgraded to the O(\alpha) QED corrections. The aim of these comparisons is to
prove the correctness of implementation of the SANC electroweak modules into
the WINHAC generator. This is achieved through the presented tuned comparisons.

With the imminent start of LHC experiments, development of phenomenological
tools, and in particular the Monte Carlo programs and algorithms, becomes
urgent. A new algorithm for the generation of a parton shower initiated by the
single initial hadron beam is presented. The new algorithm is of the class of
the so called ``constrained MC'' type algorithm (an alternative to the backward
evolution MC algorithm), in which the energy and the type of the parton at the
end of the parton shower are constrained (predefined). The complete kinematics
configurations with explicitly constructed four momenta are generated and
tested. Evolution time is identical with rapidity and minimum transverse
momentum is used as an infrared cutoff. All terms of the leadinglogarithmic
approximation in the DGLAP evolution are properly accounted for. In addition,
the essential improvements towards the socalled CCFM/BFKL models are also
properly implemented. The resulting parton distributions are crosschecked up
to the 0.1% precision level with the help of a multitude of comparisons with
other MC and nonMC programs. We regard these tests as an important asset to be
exploited at the time when the presented MC will enter as a building block in a
larger MC program for W/Z production process at LHC.

We discuss precision Monte Carlo (MC) calculations for solving the QCD
evolution equations up to the nexttoleadingorder (NLO) level. They employ
forward Markovian Monte Carlo algorithms, which provide rigorous solutions of
the above equations. These algorithms are implemented in the form of the Monte
Carlo program EvolFMC. This program has been crosschecked with independent,
nonMC, programs (QCDNum16 and APCheb33) and the numerical agreement at the
level of 0.1% has been found.

A systematic extension of the Monte Carlo (MC) algorithm, that solves the
DGLAP equation, into the socalled the oneloop CCFM evolution is presented.
Modifications are related to a zdependent coupling constant; transverse
momentum dependence is added to the xdependence of the parton distributions.
The presented Markovian algorithm for oneloop CCFM evolution is the first step
in extending it to other more sophisticated schemes beyond DGLAP. In
particular, implementing the complete CCFM will be the next step. The presently
implemented oneloop CCFM option will be a useful tool in testing the
forthcoming MC solutions. Numerical results of the new MC are confronted with
other nonMC numerical solutions. The agreement within the MC statistical error
of ~0.1% is found. Also, numerical results for kTdependent structure functions
are presented.

In this paper we propose a strategy for measuring the inclusive Wboson
production processes at LHC. This strategy exploits simultaneously the unique
flexibility of the LHC collider in running variable beam particle species at
variable beam energies, and the configuration flexibility of the LHC detectors.
We propose their concrete settings for a precision measurement of the Standard
Model parameters. These dedicated settings optimise the use of the Z boson and
DrellYan pair production processes as ``the standard reference candles''. The
presented strategy allows to factorise and to directly measure those of the QCD
effects which affect differently the W and Z production processes. It reduces
to a level of 10^{4} the impact of uncertainties in the partonic distribution
functions (PDFs) and in the transverse momentum of the quarks on the
measurement precision. Last but not the least, it reduces by a factor of 10 an
impact of systematic measurement errors, such as the energy scale and the
measurement resolution, on the Wboson production observables.

We present precision Monte Carlo calculations solving the QCD evolution
equations up to the nexttoleadingorder (NLO) level. They employ forward
Markovian Monte Carlo (FMC) algorithms, which provide the rigorous solutions of
the QCD evolution equations. Appropriate Monte Carlo algorithms are described
in detail. They are implemented in the form of the Monte Carlo program EvolFMC,
which features the NLO kernels for the QCD evolution. The presented numerical
results agree with those from independent, nonMC, programs (QCDNum16,
APCheb33) at the level of 0.1%. In this way we have demonstrated the
feasibility of the precision MC calculations for the QCD evolution and provided
very useful numerical tests (benchmarks) for other, nonMarkovian, MC
algorithms developed recently.

The LHC has been designed as a collider of proton and ion beams. However, in
its experimental program, which is focused mainly on studies of high energy
transfer collisions of Standard Model pointlike particles, protons and ions
will play a backstage role. For the majority of the LHC experimentalists, their
role will be confined to providing standardized, accelerationprocessstable
envelopes for tunabledensity and tunableisospin bunches of Standard Model
constituents: quarks and gluons. The interbunch environment of collisions of
these Standard Model particles specific to hadronic colliders and absent in the
leptonic ones, has always been considered as a burden  an annoying but
unavoidable price to pay for increasing the collisionenergy of pointlike
particles in the storage rings. In this paper we shall argue that such a burden
can be converted into an important merit of the highenergy hadronic colliders
 a cornerstone for a fermilengthscale ``collisionexperiment'' employing
the bunches of spectator quarks and gluons as tunable ``femtoscopic'' targets
for the beams of shortliving electroweak bosons.

We present the comparisons of two independent Monte Carlo event generators,
HORACE and WINHAC, for singleWboson production in hadronic collisions with
multiphoton effects in leptonic W decays. These comparisons were performed
first at the parton level with fixed quarkbeams energy, and then at the hadron
level for protonproton collisions at the LHC. In general, a good agreement
between the two programs has been found. Possible sources of differences in
some of the presented results are discussed. We also present and discuss the
effects of including nonzero quark masses for the main singleWboson
observables at the LHC.

We present the calculation of multiphoton radiation effects in leptonic
Wboson decays in the framework of the YennieFrautschiSuura exclusive
exponentiation. This calculation is implemented in the Monte Carlo event
generator WINHAC for single Wboson production in hadronic collisions at the
parton level. Some numerical results obtained with the help of this program are
also presented.

We present the recent developments in the precision studies of Wpair and
singleW processes in e+e collisions achieved with the help of the KoralW and
YFSWW3 Monte Carlo generators. We focus on the theoretical precision of the
measurements of M_W and anomalous couplings on the example of lambda coupling.
We present the mechanism of running these two independent codes in the form of
one Concurrent Monte Carlo code. We describe also the extensions of KoralW
necessary to emulate the kinematical region of singleW process.