JHEP03(2017)157
Published for SISSA by SpringerReceived: January 20, 2017 Accepted: March 16, 2017 Published: March 29, 2017
Measurement of charged-particle distributions
sensitive to the underlying event in
√
s = 13 TeV
proton-proton collisions with the ATLAS detector at
the LHC
The ATLAS collaboration
E-mail: atlas.publications@cern.ch
Abstract: We present charged-particle distributions sensitive to the underlying event, measured by the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV, in low-luminosity Large Hadron Collider fills corresponding to an integrated luminosity of 1.6 nb−1. The distributions were constructed using charged particles with ab-solute pseudorapidity less than 2.5 and with transverse momentum greater than 500 MeV, in events with at least one such charged particle with transverse momentum above 1 GeV. These distributions characterise the angular distribution of energy and particle flows with respect to the charged particle with highest transverse momentum, as a function of both that momentum and of charged-particle multiplicity. The results have been corrected for detector effects and are compared to the predictions of various Monte Carlo event gen-erators, experimentally establishing the level of underlying-event activity at LHC Run 2 energies and providing inputs for the development of event generator modelling. The cur-rent models in use for UE modelling typically describe this data to 5% accuracy, compared with data uncertainties of less than 1%.
Keywords: Hadron-Hadron scattering (experiments)
JHEP03(2017)157
Contents1 Introduction 1
2 Underlying-event observables 2
3 Monte Carlo modelling of the underlying event 4
4 The ATLAS detector 7
5 Event and object selection 8
6 Correction to particle level 9
6.1 Event and track weighting 9
6.2 Re-orientation correction 11
7 Systematic uncertainties 12
8 Results 13
9 Conclusions 20
The ATLAS collaboration 25
1 Introduction
To perform precise Standard Model measurements or to search for new physics phenomena at hadron colliders, it is important to have a good understanding not only of the primary short-distance hard scattering process, but also of the accompanying interactions of the rest of the proton-proton collision — collectively termed the underlying event (UE). As the UE is an intrinsic part of the same proton-proton collision as any “signal” partonic interaction, accurate description of its properties by Monte Carlo (MC) event generators is important for the LHC physics programme.
In modelling terms, the UE can receive contributions from initial- and final-state radia-tion (ISR, FSR), from the QCD evoluradia-tion of colour connecradia-tions between the hard scattering and the beam-proton remnants, and from additional hard scatters in the same p-p collision, termed multiple partonic interactions (MPI). As it is significantly influenced by physics not currently calculable from first principles, the measurement of the UE’s properties is crucial not only for better understanding of the mechanisms involved but also to provide input for empirical tuning of the free parameters of phenomenological UE models in MC event generators.
JHEP03(2017)157
It is impossible to uniquely separate the UE from the hard scattering process on an event-by-event basis, but observables can be defined which are particularly sensitive to the properties of the UE. Measurements of such observables have been performed in pp collisions between √s = 900 GeV and 7 TeV in ATLAS [1–5], ALICE [6] and CMS [7–9]. UE observables were also previously measured in p¯p collisions in dijet and Drell-Yan events at CDF, with centre-of-mass energies of√s = 1.8 TeV [10, 11] and 1.96 TeV [12].
In this paper we report the measurement of UE observables with the ATLAS detec-tor [13] at the LHC, using charged particles in 1.6 nb−1 of proton-proton collisions at a centre-of-mass energy of 13 TeV.
2 Underlying-event observables
The UE observables in this study are constructed from “primary” charged particles in the pseudorapidity range |η| < 2.5,1 whose transverse momentum (pT) is required to be
greater than 500 MeV. Primary charged particles are defined as those with a mean lifetime τ > 300 ps, which are either directly produced in pp interactions or from subsequent decays
of particles with a lifetime τ < 30 ps. This measurement follows the fiducial particle
definition used in the ATLAS 13 TeV minimum-bias measurement [14] excluding particle
species with τ in the range 30 ps to 300 ps and their decay products. The charged particles falling in this range are strange baryons with a very low reconstruction probability, whose decays are inconsistently modelled in MC generators. Their exclusion from the fiducial acceptance definition hence avoids the need to apply large and poorly defined corrections to particle level for these species, making the measurement more accurate.
This measurement uses the established form of UE observables [1–7,9,10,12], in which the azimuthal plane of the event is segmented into several distinct regions with differing sen-sitivities to the UE. As illustrated in figure1, the azimuthal angular difference with respect to the leading (highest-pT) charged particle, |∆φ| = |φ−φlead|, is used to define the regions:
• |∆φ| < 60◦, the “towards region”;
• 60◦ < |∆φ| < 120◦, the “transverse region”; and • |∆φ| > 120◦, the “away region”.
As the scale of the hard scattering increases, the leading charged particle acts as a conve-nient indicator of the main flow of hard-process energy. The towards and away regions are dominated by particle production from the hard process and are hence relatively insensi-tive to the softer UE. In contrast, the transverse region is more sensiinsensi-tive to the UE, and observables defined inside it are the primary focus of UE measurements.
1
ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).
JHEP03(2017)157
∆φ−∆φ
leading charged particle
towards |∆φ| < 60◦ away |∆φ| > 120◦ transverse (max) 60◦< |∆φ| < 120◦ transverse (min) 60◦< |∆φ| < 120◦
Figure 1. Definition of regions in the azimuthal angle with respect to the leading (highest-pT) charged particle, with arrows representing particles associated with the hard scattering process and the leading charged particle highlighted in red. Conceptually, the presence of a hard-scatter particle on the right-hand side of the transverse region, increasing itsPpT, typically leads to that side being identified as the “trans-max” and hence the left-hand side as the “trans-min”, with maximum sensitivity to the UE.
A further refinement is to distinguish on a per-event basis between the more and the less active sides of the transverse region [15, 16], defined in terms of their relative scalar
sums of primary charged-particle pTand termed “trans-max” and “trans-min” respectively.
The trans-min region is relatively insensitive to wide-angle emissions from the hard process, and the difference between trans-max and trans-min observables (termed the “trans-diff”) hence represents the effects of hard-process contamination. In this analysis, an event must have a non-zero primary charged-particle multiplicity in the trans-min region in order to be included in either the trans-min, -max, or -diff observables.
The variables measured in this analysis, constructed using charged particles with pT>
0.5 GeV and |η| < 2.5, with a higher-pT requirement of pleadT > 1 GeV placed on the
leading charged particle, are described in table 1. These variables are divided into two
groups: first the “binned” per-event or per-particle quantities used to define the horizontal
axes in the unfolded observables of section 8; and secondly the “averaged” mean values
of distributions of per-event quantities to be studied as functions of the binned variables — a construction known as a “profile”. The second-group variables are defined for each
bin and (except for hmean pTi, in which the η–φ area factors cancel) are scaled by the
JHEP03(2017)157
Symbol Description
Binned variables
pleadT Transverse momentum of the leading charged particle Nch(transverse) Number of charged particles in the transverse region
|∆φ| Absolute difference in particle azimuthal angle from the leading charged particle Averaged variables
hNch/δηδφi Mean number of charged particles per unit η–φ (in radians) hPpT/δηδφi Mean scalar pTsum of charged particles per unit η–φ (in radians)
hmean pTi Mean per-event average pTof charged particles (≥ 1 charged particle required)
Table 1. Definitions of the measured observables in terms of primary charged particles. The upper group of observables are used to define the x-axes of the plots to be shown in section8, and the lower group are the mean values of distributions constructed in each x-bin and plotted on the y-axes. The δηδφ scale factors convert the raw measurements of regional Nch andP pTinto densities per unit η–φ and their values change depending on the region/bin-sizes being considered, so that the density variables are everywhere directly comparable.
including between various experiments and collider energies. The area factor δφ is 2π/3 for the toward, transverse & away regions, whereas it is π/3 for the single-sided trans-min & trans-max regions, and 2π/nbins for each of the nbins equally sized bins in distributions
plotted against |∆φ|. Due to the |η| < 2.5 fiducial acceptance and the η-independence of the region definitions, δη = 5 in all cases. The profile observables are implemented as profile histograms, presenting the mean values of the “averaged” variables as measured in each bin of another observable. These hence measure the degree of correlation between two event features, either between the UE and hard scattering, or between different UE
aspects. The mean charged-particle momentum hpTi is constructed on an event-by-event
basis and then averaged over all events to give hmean pTi.
The majority of underlying-event observables study the dependences of the averaged quantities on the transverse momentum of the leading object — here the leading charged particle. The development of this from low to high pTcorresponds to the smooth transition
from “minimum bias” interactions to the hard-scattering regime focused on by most LHC analyses, and the correlation distributions characterise how soft QCD effects co-evolve with the hard process through this transition. This analysis also studies the dependence of the observables on the azimuthal angle with respect to the leading particle and each region’s charged-particle multiplicity. For the observables studied as a function of relative azimuthal angle, the leading particle is excluded from the spectrum.
3 Monte Carlo modelling of the underlying event
As a physics process related to the bulk structure of protons and not calculable from first-principles perturbative QCD, the underlying event is modelled in Monte Carlo event generator programs by various phenomenological approaches. The scattering subprocess type which contributes most to UE observables is non-diffractive inelastic scattering. In event generator implementations, model-specific choices are made to regulate these
pro-JHEP03(2017)157
cesses’ QCD divergences at low scales. In this analysis the observable definitions restrict the effect of diffractive scattering, i.e. colour-singlet exchange, to play a minor role: the hpTi profile observables have been found to be completely insensitive to diffraction, while the Nch and PpT profiles exhibit a 2% effect at low pleadT and the whole of the Nch vs. ∆φ
distribution is affected by diffraction at a 1–2% level.
This section reviews relevant features of the Pythia 8 [17, 18], Herwig 7 [19] and
Epos [20, 21] MC event generator models, which are used in this study either for data
correction or for comparison to the final corrected data distributions. A summary of the Monte Carlo generator configurations used is given in table2.
Pythia 8.18 / 8.21: the Pythia generator family is very widely used at the LHC and elsewhere for event modelling, and Pythia 8 is its most recent release series. Its key features are leading-logarithmic initial- and final-state parton showers and a Lund string hadronisation model, in addition to particle decays and soft-QCD modelling. Hard-process calculations are performed either via an internal set of leading-order matrix elements, or by reading externally generated hard scattering events.
The Pythia 8 approach to soft-QCD modelling uses a parameterised total pp cross-section, which is split into elastic, diffractive, and non-diffractive (ND) inelastic scat-tering subprocesses. The first two of these again have parameterised cross-sections, and the hadronic ND inelastic cross-section is set by subtracting their contributions from the total.
The ND contribution is modelled using perturbative 2 → 2 QCD matrix elements, dominated by t-channel gluon exchange. As the partonic cross-section for this exceeds
the hadronic ND cross-section at low pT, the existence of MPI is implied and may
be used to regularise the cross-section growth. Pythia 8 uses an evolved form of the Sj¨ostrand-van Zijl MPI model [22] in which the eikonal formalism is applied to give a model with a Poisson distribution of multiple perturbative scatterings whose mean rate depends on the scale of the hard process (interpreted as the reciprocal of the pp impact parameter), the proton form factor, and the ratio of hadronic to partonic ND cross-sections. An ansatz is used to regularise cross-section growth at
low pT, with a weak power-law evolution of the regularisation parameter with
√ s. The current form of the model also interleaves MPI emissions with parton shower evolution, allows several forms of matter overlap parameterisation, and interacts with a non-perturbative annealing procedure for reconfiguration of colour strings during hadronisation (the “colour reconnection” mechanism).
As MPI and hadronisation modelling are phenomenological and even the perturbative parton shower formalism has some configurational freedom, many parameter optimi-sations (“tunes”) of Pythia 8 have been performed. The following are considered in this study:
• ATLAS’s dedicated underlying-event tune is “A14” [23]. Its configuration is
based on the NNPDF2.3 LO [24] parton density function (PDF), and was
observ-JHEP03(2017)157
Generator Version Tune PDF Focus From
Pythia 8 8.185 A2 MSTW2008 LO MB ATLAS
Pythia 8 8.185 A14 NNPDF2.3 LO UE ATLAS
Pythia 8 8.186 Monash NNPDF2.3 LO MB/UE Authors
Herwig 7 7.0.1 UE-MMHT MMHT2014 LO UE/DPS Authors
Epos 3.4 LHC — MB Authors
Table 2. Details of the MC models used. Some tunes are focused on describing the minimum-bias (MB) distributions better, while the rest are tuned to describe underlying event (UE) or double parton scattering (DPS) distributions.
ables, including jet, Drell-Yan, and t¯t data, with an emphasis on high-scale
events.
• An alternative tune, “Monash” [25], is used for comparison. Like A14, it was
constructed using Drell-Yan and underlying-event data from ATLAS, but also data from CMS, from the SPS, and from the Tevatron in order to constrain energy scaling. It also uses the NNPDF 2.3 LO PDF, and has more of a
general-purpose / low-pT emphasis than A14. This tune gives an excellent description
of the ATLAS 7, 8 and 13 TeV minimum-bias pT spectra [14,26,27].
• The ATLAS minimum-bias tune “A2” [28] was used for deriving detector
cor-rections. This is based on the MSTW2008 LO PDF [29], and was tuned using
ATLAS minimum-bias data at 7 TeV for the MPI parameters, in addition to the older Pythia 8 tune “4C” values for fragmentation parameters. It provides a good description of minimum-bias and transverse energy flow data [30].
The Pythia 8 predictions shown in this paper use large MC samples from versions 8.185 and 8.186, but checks against the newer 8.2xx release series (specifically, version 8.210) found no distinguishable difference.
Herwig 7.0.1: the Herwig family of MC generators has also been heavily used in many collider physics studies, and Herwig 7 is the most recent major series [19,31]. Like Pythia, it is a fully exclusive hadron-level generator, containing leading-logarithmic parton showers, hadronisation and decays, an MPI mechanism, and capabilities for parton-shower matching to higher-order hard processes. It uses a cluster hadronisa-tion scheme, and angular-ordered and dipole parton showers.
The soft-QCD modelling in Herwig 7 (and Herwig++ before it) uses an eikonal
model similar to Pythia’s, but with some distinctions. The same treatment with Poisson-distributed simulation of many independent perturbative QCD scatters is then used, but with a simpler MPI parameterisation than in Pythia: the functional form of the proton electromagnetic form factor is used to parameterise the hadronic matter distribution rather than Pythia’s very flexible parameterisation, and a
JHEP03(2017)157
Pythia’s phenomenological regularisation of eikonal scattering below the pT cutoff.
The Herwig soft eikonal scattering model uses distinct matter distributions for soft
and hard MPI scattering, and introduces a √s-dependence of MPI parameters
sim-ilar to that found in Pythia. The parameters of this model are highly constrained by fits to total cross-section and elastic scattering data [31,32]. A colour-disruption mechanism is used in hadronisation, as a cluster-oriented analogy to the Pythia colour reconnection, to improve the quality of minimum-bias observable description.
The Herwig 7 default tune, H7-UE-MMHT with the MMHT2014 LO PDF [33], has
been used. This tune, like its Herwig++predecessors, is based on LHC and Tevatron
underlying-event measurements, as well as double parton scattering data [34]. It
provides a good description of all these observables for √s from Tevatron 300 GeV
to LHC 7 TeV.
Epos 3.4: an alternative approach is taken by Epos, a specialist soft-QCD/cosmic-ray air-shower MC generator based on an implementation of parton-based Gribov-Regge theory [20]. This is a QCD-inspired effective-field theory describing the hard and soft scattering simultaneously. It incorporates elements of collective flow modelling from nuclear and heavy-ion physics, using hydrodynamic flow modelling in high-density regions and a string-based hadronisation model elsewhere. As a result the Epos cal-culations do not make use of standard parton density functions. Using this version, equivalent to the so-called “LHC tune” of version 1.99 [21], Epos gives a very good description of ATLAS’s 13 TeV minimum-bias data, including the tails of distribu-tions where UE physics should be involved, but as it lacks a dedicated hard scattering component it is unclear how accurate its description of UE correlations can be. Detector-level simulation
The Pythia 8 A2 sample and an MC simulation of single particles distributed to populate
the high-pT region were used to derive the detector corrections for these measurements.
Smaller samples produced with the Epos generator were used to cross-check the validity of the correction procedure.
These events were processed through the ATLAS detector simulation framework [35],
which is based on Geant 4 [36]. They were then reconstructed and analysed using the
same processing chain as for data. In all Monte Carlo samples the distribution of the primary vertex position was reweighted to match the distribution in data.
4 The ATLAS detector
The ATLAS detector is described in detail in ref. [13]. In this analysis, the trigger system and inner tracking detectors are of particular relevance.
The inner tracking detector is immersed in the 2 T axial magnetic field of a supercon-ducting solenoid, and measures the trajectories of charged particles in the pseudorapidity range |η| < 2.5 with full azimuthal coverage. It consists of a silicon pixel detector (pixel), a silicon microstrip detector (SCT) and a straw-tube transition radiation tracker (TRT),
JHEP03(2017)157
split into a barrel and two endcap components. The barrel consists of 4 pixel layers, 4 layers of SCT modules with back-to-back silicon strip sensors, and 73 layers of TRT straws; each endcap has 3 pixel layers, 9 SCT layers, and 160 TRT straw layers. The pixel, SCT and
TRT have r-φ position resolutions of 10µm, 17 µm, and 130 µm respectively, and the pixel
(apart from the innermost barrel layer) and SCT have resolutions of 115µm and 580 µm
respectively in the z-direction for the barrel modules and in the r-direction for the endcaps. The innermost pixel layer, the “insertable B-layer” (IBL), was added between LHC
Runs 1 and 2, around a new narrower (radius of 25 mm) and thinner beam pipe [37]. It is
composed of 14 lightweight staves arranged in a cylindrical geometry, each made of 12 silicon planar sensors in its central region and 2×4 three-dimensional sensors at the ends. The IBL
pixel dimensions are 50 × 250µm in the φ and z directions (compared with 50 × 400 µm for
other pixel layers). The smaller radius and the reduced pixel size result in improvements of both the transverse and longitudinal impact parameter resolutions. In addition, the services for the existing pixel detector were upgraded, significantly reducing the material in |η| > 1.5, in particular at the boundaries of the active tracking volume. A track from a charged particle traversing the barrel detector typically has 12 silicon measurement points (hits), of which 4 are pixel and 8 are SCT, and more than 30 TRT straw hits.
The ATLAS detector has a two-level trigger system: level 1 and the high-level trig-ger [38]. Events used in this analysis were required to satisfy level-1 triggers using the minimum-bias trigger scintillators (MBTS). These are mounted at each end of the detector in front of the liquid-argon endcap-calorimeter cryostats at z = ±3.56 m and are segmented into two rings in pseudorapidity (2.07 < |η| < 2.76 and 2.76 < |η| < 3.86), with 8 azimuthal sectors in the inner ring and 4 in the outer. The MBTS scintillator was replaced for Run 2 due to radiation damage incurred in Run 1. The MBTS triggers fired if at least one MBTS hit from either side of the detector was recorded above threshold.
5 Event and object selection
This measurement uses data taken in a special configuration of the LHC, with low beam currents and reduced beam focusing, producing a low mean number of interactions per bunch-crossing, hµi, between 0.003 and 0.03. This configuration and event selection were
earlier used and documented in detail in the ATLAS 13 TeV minimum-bias analysis [14],
and so only a summary is given here.
Trigger: Events were selected from colliding proton bunches using a trigger which
re-quired one or more MBTS counters above threshold on either side of the detector. The efficiency of this trigger was measured to be 99% for low-multiplicity events and to rapidly rise to 100% for events with higher track multiplicities. The trigger requirement does not bias the pT and η distributions of selected tracks in this analysis due to the pleadT > 1 GeV
requirement.
Primary vertex: Each event was required to contain a primary vertex reconstructed
from at least two tracks with pT> 100 MeV and selection requirements specific to
vertex-ing [39]. The canonical primary vertex was identified as that with the highest P p2 T of its
JHEP03(2017)157
associated tracks. To reduce contamination from events with more than one interaction in a bunch crossing (“pile-up”), events containing more than one primary vertex with four or more associated tracks were removed. The contributions from non-collision backgrounds and events where two interactions are reconstructed as a single vertex were studied in data and found to be negligible.
A total of 66 million data events passed the trigger and vertex selection requirements for this analysis.
Tracks: Tracks were reconstructed starting from hits in the silicon detectors and then
extrapolated to include information from the TRT. Each track required hits in both the pixel system and the SCT, including a requirement of a hit in the innermost expected pixel layer to reject secondary particles.
All tracks were reconstructed within the geometric acceptance of the inner detector, |η| < 2.5. In this analysis, all selected tracks were additionally required to have a transverse momentum above 500 MeV, and both the transverse impact parameter and the projected longitudinal impact parameter with respect to the primary vertex were required to be less
than 1.5 mm. For high-pT tracks, a further requirement was placed on the χ2 probability
of the track fit, to suppress mismeasured tracks. The selected events were additionally required to contain at least one track with a transverse momentum above 1 GeV.
Backgrounds: The contributions from non-collision background events, events with
more than one interaction, and fake tracks (those formed by a random combination of hits or from a combination of hits from several particles) were found to be negligible. The contribution from secondary particles was estimated as in ref. [14].
6 Correction to particle level
In order to compare the measured underlying-event distributions with model predictions in the particle-level fiducial phase space, the observables have been corrected for detector effects. These corrections include explicit accounting for inefficiencies due to the trigger selection, vertexing, and track reconstruction, by weighting the events and tracks by in-verse efficiencies. A further correction has been applied to account for non-linear effects, particularly azimuthal re-orientation of the event which occurs should the leading charged particle not be reconstructed. In this situation the identification of the towards, transverse, and away regions can differ from that at particle level, leading to “wrong” track-region as-sociations. Re-orientation mainly affects events with a low number of charged particles, and has been corrected using a dedicated method. The track-weight correction is described first, followed by the method to account for event re-orientation.
6.1 Event and track weighting
The weighting procedure for correcting measured distributions to particle level is affected by the primary vertex reconstruction efficiency, the tracking efficiency, the rate of non-primary particles, and the rate of charged strange baryons. Before defining the weights, we summarise these aspects of the vertexing and tracking performance:
JHEP03(2017)157
Vertexing efficiency: The vertex reconstruction efficiency, vtx, was determined from
data by taking the ratio of the number of selected events with a reconstructed vertex to the total number of events with the requirement of a primary vertex removed. The efficiency was close to 100% except for events in which exactly one pT> 500 MeV analysis track was
found in the inner detector — for these the vertexing efficiency was found to be 90%. Since the leading-track pT > 1 GeV requirement in the analysis is correlated with higher track
multiplicities, the effect of the vertex reconstruction on the analysis is small.
Tracking efficiency: The efficiency to reconstruct primary charged particles was
de-termined from simulation by matching MC-generator primary charged particles to tracks reconstructed from simulated inner-detector hits. For particles at the analysis
transverse-momentum cutoff of pT = 500 MeV, the track reconstruction is ∼ 85% efficient for central
pseudorapidity, and decreases to 65% for |η| ∼ 2.5 at the edge of the inner tracking de-tector; these efficiencies are higher for tracks with higher pT, for example 90% and 80%
respectively for pT= 10 GeV.
Most of the loss in tracking efficiency is due to particle interactions with the detector material and support structures, and hence a good description of the detector material is needed. The ATLAS inner tracker upgrades for LHC Run 2 necessitated a new study of the detector material uncertainties, which are detailed in ref. [14]. The track reconstruction efficiency for |η| > 1.5 is corrected using a method that compares the efficiency to extend a track reconstructed in the pixel detector into the SCT in data and simulation.
Rate of non-primary tracks: There are several sources of non-primary tracks: 1) tracks
from hadronic interaction of particles with the detector material, 2) decay products from
particles with strange quark content, mostly K0 and Λ0 decays, 3) photon conversions.
The last of these is negligible for the tracks with pT > 500 MeV used in this analysis. The
rate of non-primary tracks was determined by comparing side-band fits of MC simulation impact parameter distributions to the data, leading to an estimated 2.3% of non-primary tracks in data for the pT > 500 MeV track selection, and a smaller effect for higher-pT
tracks.
Fraction of charged strange baryons: Charged strange baryons and their decay
prod-ucts have low reconstruction efficiency, unless they have large transverse momentum. With the fiducial primary particle definition used here, the fraction of tracks due to strange
baryons with pT = 10 GeV is 1%, while it is much smaller than 1% at lower pT. With
our fiducial primary particle definition, tracks originating from charged strange baryons are classified as a background. An estimate of their contribution was made using the Epos generator, a choice motivated both by Epos providing the best current
descrip-tion of charged strange baryon rates as measured by the ALICE experiment [40], and by
substantial variation of strange baryon modelling between different MC generators.
Weighting: The effect of events lost due to the trigger and vertex requirements was
corrected by an event-by-event weight, wev nBLsel, η = 1 trig nBLsel · 1 vtx nBLsel, η , (6.1)
JHEP03(2017)157
where trigis the trigger efficiency and nBLsel is the multiplicity of “beam line” selected tracks,
which have the same selection requirements as analysis tracks except that no vertexing requirement is applied and and hence there is no restriction on the longitudinal impact parameter. To correct for inefficiencies in the track reconstruction, the distributions of the selected analysis tracks were corrected with track-by-track weights,
wtrk(pT, η) =
1 trk(pT, η)
· (1 − fnonp(pT, η) − fokr(pT, η) − fsb(pT)) , (6.2)
where trk is the tracking efficiency, and fnonp, fokr and fsb are respectively the fractions
of non-primary tracks, of out-of-kinematic-range tracks (i.e. tracks mis-reconstructed as within the fiducial pT and η acceptances, but which actually originated from outside those
ranges), and of weakly decaying charged strange baryons.
This track weight was used in the construction of the Nch, P pT, and ∆φ
distri-butions, and the mean pT in the event, determined using the weighted average hpTi =
P
i∈ trackspTiwtrki /
P
i∈ trackswtrki .
For the distributions binned in charged multiplicity, this weight was not used to correct the binning variable, since this induces migrations subject to fluctuations not accounted for by track weighting. This correction was instead handled by the residual correction method described in the next section.
6.2 Re-orientation correction
The re-orientation correction was based on the HBOM method [41], in which the effect of
the detector and reconstruction algorithms on an observable (i.e. a histogram bin value) is treated as an operator A. The observed reconstruction-level value of an observable, Xobs, is hence related to its true value Xtrue as Xobs = AXtrue, but one can continue to apply A to the modified event, so that the kth additional application gives Xkobs = AkXobs = Ak+1Xtrue. If the evolution of the observable’s value is a smoothly varying function under such iterations, an nth-order polynomial function Pn(k) can be fitted to Xkobsfor k ≥ 0, and
then extrapolated back to k = −1, i.e. Xtrue. This procedure is carried out independently for each bin of each distribution. It is distinct from the unfolding approach used in refs. [4] and [14], chosen because it is fully data-driven, i.e. does not rely on simulation performance or require reweighting of MC events to data.
In this analysis, A encodes the effects of track reconstruction, i.e. each application of A in principle smears the track kinematics and considers the possibility that some tracks would not have been reconstructed — directly inducing re-orientation should the leading track experience such a fate. These effects are described by the measured tracking res-olutions and efficiencies, as functions of track η and pT, but the efficiency was found to
dominate and the resolution to be negligible, so A reduces to random discarding of tracks according to their reconstruction efficiency. For the profile distributions, a weight propor-tional to the inverse of tracking efficiency (as above) was applied in each HBOM iteration step, k, to make the effect on the observable less dramatic and more easily parameterisable. To avoid correlations between observables with different k values, each Xkobs was calculated independently, using a distinct random seed and k iterations of track discarding starting
JHEP03(2017)157
from the set of reconstructed (k = 0) tracks. Second-order polynomials and a maximum of k = 4 HBOM iterations were used as the nominal configuration, with variations being used to derive an extra systematic uncertainty, described in the following section.
7 Systematic uncertainties
Several categories of systematic uncertainties that may influence the distributions after corrections and unfolding were quantified, and their magnitudes are summarised in table3. The sources of uncertainty and the methods used to estimate them were the following: Trigger and vertexing: the systematic uncertainties were found to be negligible. Track reconstruction: the uncertainties in track reconstruction efficiency principally
arise from imperfect knowledge of the material in the inner detector. The new in-sertable B-layer and changes to pixel detector services in the |η| > 1.5 region add un-certainties not included in the inner-detector material assessments performed during Run 1: these have been evaluated by comparison of data to simulation, and by com-parisons of simulations with different compositions of interacting particles. The result is an estimate of reconstruction efficiency uncertainty between 0.4% and 1.5% [14]. An additional systematic uncertainty arises from possible biases and degradation in the plead
T measurement: these effects were determined in ref. [14] and estimated to
affect the tail of the pleadT distribution by 4–5%.
Non-primary particles: the systematic uncertainties were propagated by modification of track weights. The systematic uncertainty on the selected non-primary track fraction was estimated to be 24%, using variations of the fit range in the tail of the impact parameter distributions, and different MC generators and different shape assumptions for the extrapolation of the fraction from the side-bands to the signal region. The MC variation is responsible for the most significant contribution to this uncertainty. The systematic uncertainty due to the fraction of charged strange baryons was derived using the deviations of generator predictions from Epos [14].
Unfolding: the systematic uncertainties associated with the HBOM unfolding have two distinct sources:
Non-closure: the HBOM correction procedure is in principle independent of the Monte Carlo generator modelling. However, the method shows deficiencies in some regions and distributions. The relative size of the observed closure (i.e. non-reproduction of a known input) derived using Pythia 8 A2 was therefore applied as a correction. The size of the correction is included as a systematic uncertainty, everywhere less than 2.5%. The measured transverse momentum distribution of the leading charged particle was found to be consistent with the result obtained using Bayes-inspired iterative unfolding [42] within the estimated bias.
Parameterisation: the statistical uncertainty of the HBOM method was derived for each bin using Monte Carlo sampling: for each HBOM(k) iteration, samples were
JHEP03(2017)157
Range of values
Observable Material Non-primaries Non-closure Parameterisation
Nch orP pT vs. ∆φ 0.9% 0.6% 0–0.6% 0–0.4%
Nch orP pT vs. pleadT 0.5–1.0% 0.3–0.6% 0–2.5% 0–0.4%
hmean pTi vs. Nch 0–0.5% 0–0.5% — 0.5% (combined) —
hmean pTi vs. plead
T 0–0.4% 0–0.3% — 0.5% (combined) —
Table 3. Summary of systematic uncertainties for each class of UE observable, broken down by origin.
drawn from a Gaussian distribution corresponding to the point uncertainty, and 1000 independent replica HBOM fits were performed. The resulting uncertainty contri-bution is derived from the replicas’ central 68% confidence interval. An additional systematic uncertainty, estimating the stability of the fit method, was derived by using different numbers of HBOM iterations k and polynomial degrees n.
For the results presented in the following section, these independent sources of systematic uncertainty have been combined in quadrature to form single representative systematic uncertainty estimates.
8 Results
The unfolded distributions and their main features are discussed here, and in figures2,3,5
and6are compared to model predictions from Pythia 8 (the A14, A2 and Monash tunes),
Herwig 7, and Epos. These model predictions were obtained using the analysis’ associated Rivet routine [43].
Leading charged-particle pT: Figure 2 shows the unit-normalised distribution of
events with respect to the transverse momentum of the leading charged particle, pleadT . This is a steeply falling distribution, with a change of slope for pleadT & 5 GeV: a form which is broadly modelled by all generators. The Pythia 8 A14 and Monash tunes, as well as Epos, model the distribution within 15% out to peaks strongly at the lowest pleadT
and alternates between under- and over-shooting the data at higher scales, finally producing a softer tail than seen in data.
Angular distributions versus leading charged particle: Figure 3 shows the mean
multiplicity and P pT distributions as a function of azimuthal angle with respect to the
leading particle for different pleadT requirements. Two event selections are shown here: the plead
T > 1 GeV cut common to all observables, and a harder pleadT > 10 GeV requirement.
The difference between these two selections illustrates the transition from relatively isotropic minimum-bias scattering to the emergence of hard partonic scattering structure and hence a dominant axis of energy flow. This event structure with least activity per-pendicular to the leading-object axis, i.e. away from ∆φ = 0◦ and 180◦, is seen for both
JHEP03(2017)157
] -1 [GeV lead T / d p ev dN ev 1/N 5 − 10 4 − 10 3 − 10 2 − 10 1 − 10 1 10 2 10 3 10 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 1 1.2 1.4Figure 2. Unit-normalised distribution of the transverse momentum of the leading charged particle, plead
T > 1 GeV, compared to various generator models. The error bars on data points represent statistical uncertainty and the blue band the total combined statistical and systematic uncertainty.
〉 φ δ η δ / ch N〈 1 10 ATLAS | < 2.5 η > 0.5 GeV, | T p -1 = 13 TeV, 1.6 nb s > 10 GeV lead T p > 1 GeV lead T p PYTHIA 8 A14 PYTHIA 8 Monash Herwig7 Epos MC / Data 0.8 0.9 1 1.1 lead> 10 GeV T p | [degrees] φ ∆ | 0 20 40 60 80 100 120 140 160 180 MC / Data 0.8 1 1.2 pleadT > 1 GeV [GeV]〉 φ δ η δ / T p Σ 〈 1 10 ATLAS | < 2.5 η > 0.5 GeV, | T p -1 = 13 TeV, 1.6 nb s > 10 GeV lead T p > 1 GeV lead T p PYTHIA 8 A14 PYTHIA 8 Monash Herwig7 Epos MC / Data 0.8 0.9 1 1.1 lead> 10 GeV T p | [degrees] φ ∆ | 0 20 40 60 80 100 120 140 160 180 MC / Data 0.8 1 1.2 pleadT > 1 GeV
Figure 3. Distributions of mean densities of charged-particle multiplicity Nch (left), and P pT (right) as a function of |∆φ| (with respect to the leading charged particle) for plead
T > 1 GeV and plead
T > 10 GeV separately, with comparisons to MC generator models. The error bars on data points represent statistical uncertainty and the blue band the total combined statistical and systematic uncertainty.
JHEP03(2017)157
[GeV] lead T p 5 10 15 20 25 30 〉 φ δ η δ / ch N〈 0 0.5 1 1.5 2 2.5 ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p Towards region Transverse region Away region [GeV] lead T p 5 10 15 20 25 30 [GeV]〉 φ δ η δ / T p Σ 〈 0 1 2 3 4 5 6 7 8 ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p Towards region Transverse region Away regionFigure 4. Mean η–φ densities of charged-particle multiplicities (left) andPpT(right) as a function of the transverse momentum of the leading charged particle in the transverse, towards, and away azimuthal regions. The error bars, which are mostly hidden by the data markers, represent combined statistical and systematic uncertainty.
selections and both observables but is much stronger for the event subset with the higher pleadT > 10 GeV requirement: this demonstrates the evolution of event shape as a hard scattering component develops.
There is no clear “best” MC model for these observables. Epos performs best in the more inclusive pleadT > 1 GeV selection, followed by Pythia 8 A14; Herwig 7 significantly undershoots while Pythia 8 Monash is everywhere above the data. But in the
hard-scattering pleadT > 10 GeV event selection, Herwig 7 and Monash perform best, with a
slight undershoot from Pythia 8 A14 and a large one from Epos. These orderings apply to both the Nch andPpT variables, although to different extents.
Nch and PpT densities in transverse/towards/away regions: Figure 4shows the
evolution of the mean charged-particle multiplicity and PpT densities with the pT of the
leading charged particle, from 1 to 30 GeV. For both observables, the towards, transverse, and away regions are shown overlaid for ease of comparison.
The primary feature is the general shape seen in all curves: starting close to zero at low plead
T , there is first a very rapid rise in activity in which the three regions are not
strongly distinguished, then an abrupt transition at pleadT ≈ 5 GeV above which the three regions have quite distinct behaviours. The initial rapid rise to a roughly stable value of ∼1 charged particle or ∼1 GeV per unit η–φ area is known as the “pedestal effect”. In modelling terms this reflects a reduction of the pp impact parameter with increasing pleadT , and hence the transition between the minimum-bias and hard-scattering regimes.
Secondly, the shape of the transverse region is different from the other two in both
observables — it almost completely plateaus after pleadT ≈ 5 GeV, while the towards and
away regions continue to rise nearly linearly with pleadT . This characteristic feature of UE profile observables is the empirical demonstration of the azimuthal-region paradigm for UE analysis: the hard process dominates the towards and away regions, which continue to increase in activity as the hard-process scale grows, but the transverse region is rela-tively unaffected. This is consistent with the pedestal effect where the overlap between
JHEP03(2017)157
colliding protons is complete and hence any further growth is due to connections to, or contaminations from, the hard process rather than more MPI scattering.
An interesting feature is that for pleadT & 7 GeV the away region actually becomes the one with highest charged-particle multiplicity, despite not containing the highest-pT
charged particle. By the PpT density measure, however, the towards region is
unambigu-ously the most active region for all pleadT values.
Nch and PpT densities in trans-min/max/diff regions: Figure 5 focuses on the
UE-dominated transverse region, and its per-event specialisations trans-min, -max, and -diff. Between these, the trans-min is the most sensitive to MPI effects, i.e. the pedestal, while the trans-max includes both MPI and hard-process contaminations: the trans-diff is hence the clearest measure of those contaminations. The comparisons to MC models are again made in these plots.
There is significant variation in performance between the models, with Pythia 8’s Monash tune and Herwig 7 giving the best description of data in the plateau region of trans-min, followed within 10% by the other Pythia 8 tunes. Epos, however, slightly overestimates in the “ramp” region to the pedestal effect plateau, and on the plateau it underestimates the pedestal height by around 20%. Herwig 7 and Pythia 8 A2 both
mismodel the transition, with a severe undershoot for Herwig below pleadT of 5 GeV, and
a milder but broader undershoot from A2 which extends up to pleadT ≈ 20 GeV.
The predictions cluster together more tightly in the process-inclusive trans-max region, with all generators other than Epos providing a description of the Nch density data within
a few percent for pleadT ≥ 10 GeV; Epos continues to undershoot the data significantly.
Pythia 8 A14 also significantly undershoots the PpT density data, by around 10% as
compared to Epos’s 20%. Looking in trans-diff gives a clearer view of how non-MPI contributions to the UE are modelled, with mostly flat ∼10% overshoots from all models other than Epos in Nchdensity, and a spread of performance in describing thePpTdensity
evolution. In the latter the best performance comes from the Pythia 8 Monash and A2 tunes, with Herwig 7 and Pythia 8 A14 both wrong by 5–10% but in opposite directions — Epos’s prediction is again separated from the ATLAS data by more than 20%.
There is no obvious best model for all observables, but the Pythia 8 Monash tune
agrees well with the data in all observables other than trans-diff Nch density, and
Her-wig 7 has comparable performance for hard-scattering events away from the “minimum-bias region” of pleadT < 5 GeV. The Pythia 8 A14 tune, used for much of the hard-process simulation in ATLAS, predicts activity 5% to 10% below the data, indicating that some re-tuning for 13 TeV event modelling may yield performance benefits. Epos is not able to model the level of underlying-event activity well in events with higher pleadT .
Mean transverse momentum in transverse & trans/min/max regions: The
per-event mean transverse momentum of charged particles in the transverse azimuthal regions is of interest since it illustrates the balance in UE physics between the P pT and
multi-plicity observables. This balance is affected in some MC models by colour-reconnection or -disruption mechanisms, which stochastically reconfigure the colour structures in the hadro-nising system into energetically favourable states and typically increase the pTper particle.
JHEP03(2017)157
〉 φ δ η δ / ch N〈 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-min region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 1 [GeV]〉 φ δ η δ / T p Σ 〈 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-min region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 1 〉 φ δ η δ / ch N〈 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-max region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 1 [GeV]〉 φ δ η δ / T p Σ 〈 0.5 1 1.5 2 2.5 3 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-max region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 1 〉 φ δ η δ / ch N〈 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-diff region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 0.9 1 1.1 1.2 [GeV]〉 φ δ η δ / T p Σ 〈 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-diff region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.8 1Figure 5. Mean densities of charged-particle multiplicity Nch(left) andPpT(right) as a function of leading charged-particle pT, in the trans-min (top), trans-max (middle) and trans-diff (bottom) azimuthal regions. The error bars on data points represent statistical uncertainty and the blue band the total combined statistical and systematic uncertainty.
JHEP03(2017)157
Measurements of the correlations of mean pTwith other event features are an important
in-put to constrain such ad hoc models. Figure6shows, for the transverse and trans-min/max
azimuthal regions, the correlations of the average per-event mean transverse momentum, hmean pTi, with the transverse region’s charged-particle multiplicity and the event’s pleadT .
The correlation with transverse charged-particle multiplicity Nch is a correlation
be-tween two different “soft” properties. This distribution is modelled to within 5% for all
Nch by most of the generators and in all transverse region variants, but in all cases an
underestimation of hmean pTi is visible until Nch & 15. The best modelling is from Epos,
whose maximum undershoot is ∼3% at low Nch but which follows the data closely in
all region definitions for higher transverse multiplicities. The Monash tune is the best-performing Pythia 8 configuration, with a performance similar to that of Epos, while Pythia 8 A14 undershoots in the low-multiplicity events and Pythia 8 A2 overshoots at high-multiplicities — notably the regions not included in each tune’s construction.
Her-wig 7 shows the largest variations, from a ∼7% undershoot at low Nch to a 5% overshoot
at Nch ≈ 30.
However, Herwig 7 performs better than all the other generators when considering the correlation between transverse region hmean pTi and pleadT . As seen in figure 5, the
transverse P pT does not reach as flat a plateau as does Nch with increasing pleadT , and
hence the event-wise hmean pTi increases with pleadT . Herwig’s behaviour in this observable
is within 1% of the data except in the “minimum bias” pleadT . 5 GeV phase space (where
the Herwig model is not expected to work), while all Pythia 8 tunes undershoot the data by between 5% and 10%, and Epos even more so.
JHEP03(2017)157
[GeV]〉 T mean p〈 0.7 0.8 0.9 1 1.1 1.2 1.3 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Transverse region (Transverse) ch N 5 10 15 20 25 Model / Data 0.95 1 1.05 [GeV]〉 T mean p〈 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Transverse region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.9 0.95 1 [GeV]〉 T mean p〈 0.7 0.8 0.9 1 1.1 1.2 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-min region (Transverse) ch N 5 10 15 20 25 Model / Data 0.95 1 1.05 [GeV]〉 T mean p〈 0.7 0.8 0.9 1 1.1 1.2 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-min region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.95 1 [GeV]〉 T mean p〈 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-max region (Transverse) ch N 5 10 15 20 25 Model / Data 0.95 1 1.05 [GeV]〉 T mean p〈 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 Data PYTHIA 8 A14 PYTHIA 8 A2 PYTHIA 8 Monash Herwig7 Epos ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p > 1 GeV lead T p Trans-max region [GeV] lead T p 5 10 15 20 25 30 Model / Data 0.9 0.95 1Figure 6. Mean charged-particle average transverse momentum as a function of charged-particle multiplicity in transverse region Nch(Transverse) (left) and as a function of pleadT (right), for each of the transverse (top), trans-min (middle) and trans-max (bottom) azimuthal regions. The error bars on data points represent statistical uncertainty and the blue band the total combined statistical and systematic uncertainty.
JHEP03(2017)157
[GeV] lead T p 5 10 15 20 25 30 〉 φ δ η δ / ch N〈 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p Transverse region = 13 TeV s = 7 TeV s = 0.9 TeV s [GeV] lead T p 5 10 15 20 25 30 [GeV]〉 φ δ η δ / T p Σ 〈 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 ATLAS -1 = 13 TeV, 1.6 nb s | < 2.5 η > 0.5 GeV, | T p Transverse region = 13 TeV s = 7 TeV s = 0.9 TeV sFigure 7. Mean charged-particle multiplicity (left) and PpT (right) densities as a function of transverse momentum of the leading charged particle measured for√s = 0.9, 7 TeV [1] and 13 TeV centre-of-mass energies. The fiducial acceptance definitions of the √s = 0.9 and 7 TeV measure-ments did not exclude charged strange baryons, but this effect is limited to a few percent at most.
9 Conclusions
We have presented several distributions sensitive to properties of the underlying event,
measured by the ATLAS experiment at the LHC using 1.6 nb−1of low-luminosity 13 TeV pp
collision events. The measured observables are defined using charged tracks and have been corrected to the level of primary charged particles. Correction of the results for detector effects was performed via a combination of tracking and vertexing efficiency weighting and HBOM unfolding, with the latter accounting for re-orientation and similar indirect effects not accounted for by the reweighting. Systematic uncertainties are included for the effects of all data processing steps, including the unfolding.
The analysis observables follow the established strategy in which each event is az-imuthally segmented into “towards”, “transverse”, and “away” regions with respect to the highest-pT(“leading”) charged particle, with a further per-event distinction into “min” and
“max” sides within the transverse region. This construction produces a set of observables collectively sensitive to each of the components of the underlying event. The observables
in this study characterise the correlations of a) each azimuthal region’s mean PpT and
mean charged-particle multiplicity densities with the transverse momentum, pleadT , of the event’s leading charged particle; and b) the per-event mean charged-particle pT with both
pleadT and the transverse azimuthal region’s charged-particle multiplicity. The correlations
with pleadT characterise the degree of connection between the underlying event and the
hardest partonic scattering in the pp collision, while correlations with the transverse re-gions’ charged-particle multiplicities probe the particle production mechanisms between underlying-event components.
The presented measurement improves upon previous ATLAS measurements of the underlying event using leading-track alignment, both in the reach in pleadT and the precision
achieved. Figure 7 shows the transverse region’s mean Nch and PpT densities in bins of
pleadT as measured in leading-charged-particle UE studies at√s = 0.9 and 7 TeV, compared with the new 13 TeV data. An increase in UE activity of approximately 20% is observed when going from 7 TeV to 13 TeV pp collisions.
JHEP03(2017)157
These results are a crucial input to physics studies in LHC pp collisions throughout the 13 TeV run, since the underlying event is present in all inelastic collision processes. Comparisons against predictions from several commonly used MC generator configura-tions indicate that for most observables the models describe the underlying-event data to better than 5% accuracy. But this level of variation is far greater than the experimental uncertainty of the measurements, and there is evidence of systematic mismodelling. The model defects are particularly acute in the transition from generic soft inelastic interac-tions — so-called “minimum bias” scattering — to secondary interacinterac-tions in the presence of hard partonic scattering. Some improvement to MC tunes in light of these data will hence be of benefit to physics studies in LHC Run 2. The Epos MC generator, specialised for simulation of inclusive soft QCD processes, displays particularly discrepant features as the pleadT scale increases, casting doubt on its suitability for modelling LHC multiple pp interactions despite currently providing the best description of minimum-bias data.
This study hence provides important data from the new collider energy frontier for the tuning of MC generators in LHC Run 2, and further incentives for the development of a nucleon collision model capable of describing collider data in both hard and soft scatterings and at multiple centre-of-mass energies.
Acknowledgments
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently.
We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Aus-tralia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DSM/IRFU, France; GNSF, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; FOM and NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian
Fed-eration; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZˇS, Slovenia; DST/NRF,
South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, the Canada Council, CANARIE, CRC, Compute Canada, FQRNT, and the Ontario Innovation Trust, Canada; EPLANET, ERC, ERDF, FP7, Horizon 2020 and Marie Sk lodowska-Curie Actions, European Union;
In-vestissements d’Avenir Labex and Idex, ANR, R´egion Auvergne and Fondation Partager
le Savoir, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF; BSF, GIF and Minerva, Israel; BRF, Norway; CERCA Programme Generalitat de Catalunya, Generalitat Valenciana, Spain; the Royal Society and Leverhulme Trust, United Kingdom.
JHEP03(2017)157
The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (U.K.) and BNL (U.S.A.), the Tier-2 facilities worldwide and large non-WLCG resource providers. Ma-jor contributors of computing resources are listed in ref. [44].
Open Access. This article is distributed under the terms of the Creative Commons
Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.
References
[1] ATLAS collaboration, Measurement of underlying event characteristics using charged particles in pp collisions at√s = 900 GeV and 7 TeV with the ATLAS detector,Phys. Rev. D
83 (2011) 112001[arXiv:1012.0791] [INSPIRE].
[2] ATLAS collaboration, Measurements of underlying-event properties using neutral and charged particles in pp collisions at 900 GeV and 7 TeV with the ATLAS detector at the
LHC,Eur. Phys. J. C 71 (2011) 1636[arXiv:1103.1816] [INSPIRE].
[3] ATLAS collaboration, Underlying event characteristics and their dependence on jet size of charged-particle jet events in pp collisions at√s = 7 TeV with the ATLAS detector,Phys.
Rev. D 86 (2012) 072004[arXiv:1208.0563] [INSPIRE].
[4] ATLAS collaboration, Measurement of the underlying event in jet events from 7 TeV proton-proton collisions with the ATLAS detector,Eur. Phys. J. C 74 (2014) 2965
[arXiv:1406.0392] [INSPIRE].
[5] ATLAS collaboration, Measurement of distributions sensitive to the underlying event in inclusive Z-boson production in pp collisions at√s = 7 TeV with the ATLAS detector,Eur.
Phys. J. C 74 (2014) 3195[arXiv:1409.3433] [INSPIRE].
[6] ALICE collaboration, Underlying event measurements in pp collisions at√s = 0.9 and 7 TeV with the ALICE experiment at the LHC,JHEP 07 (2012) 116[arXiv:1112.2082] [INSPIRE].
[7] CMS collaboration, First measurement of the underlying event activity at the LHC with√ s = 0.9 TeV,Eur. Phys. J. C 70 (2010) 555[arXiv:1006.2083] [INSPIRE].
[8] CMS collaboration, Measurement of the underlying event activity using charged-particle jets in proton-proton collisions at√s = 2.76 TeV,JHEP 09 (2015) 137[arXiv:1507.07229] [INSPIRE].
[9] CMS collaboration, Measurement of the underlying event in the Drell-Yan process in proton-proton collisions at√s = 7 TeV,Eur. Phys. J. C 72 (2012) 2080[arXiv:1204.1411] [INSPIRE].
[10] CDF collaboration, T. Affolder et al., Charged jet evolution and the underlying event in p¯p collisions at 1.8 TeV,Phys. Rev. D 65 (2002) 092002[INSPIRE].
[11] CDF collaboration, D. Acosta et al., The underlying event in hard interactions at the Tevatron ¯pp collider, Phys. Rev. D 70 (2004) 072002[hep-ex/0404004] [INSPIRE].
JHEP03(2017)157
[12] CDF collaboration, T. Aaltonen et al., Studying the underlying event in Drell-Yan and hightransverse momentum jet production at the Tevatron,Phys. Rev. D 82 (2010) 034001
[arXiv:1003.3146] [INSPIRE].
[13] ATLAS collaboration, The ATLAS experiment at the CERN Large Hadron Collider,2008
JINST 3 S08003[INSPIRE].
[14] ATLAS collaboration, Charged-particle distributions in √s = 13 TeV pp interactions measured with the ATLAS detector at the LHC,Phys. Lett. B 758 (2016) 67
[arXiv:1602.01633] [INSPIRE].
[15] G. Marchesini and B.R. Webber, Associated transverse energy in hadronic jet production,
Phys. Rev. D 38 (1988) 3419[INSPIRE].
[16] J. Pumplin, Hard underlying event correction to inclusive jet cross-sections,Phys. Rev. D 57
(1998) 5787[hep-ph/9708464] [INSPIRE].
[17] T. Sj¨ostrand, S. Mrenna and P.Z. Skands, A brief introduction to PYTHIA 8.1, Comput.
Phys. Commun. 178 (2008) 852[arXiv:0710.3820] [INSPIRE].
[18] T. Sj¨ostrand, S. Mrenna and P.Z. Skands, PYTHIA 6.4 physics and manual,JHEP 05
(2006) 026[hep-ph/0603175] [INSPIRE].
[19] J. Bellm et al., HERWIG 7.0/HERWIG++ 3.0 release note, Eur. Phys. J. C 76 (2016) 196
[arXiv:1512.01178] [INSPIRE].
[20] T. Pierog and K. Werner, EPOS model and ultra high energy cosmic rays, [Nucl. Phys. Proc.
Suppl. B 196 (2009) 102] [arXiv:0905.1198] [INSPIRE].
[21] T. Pierog, I. Karpenko, J.M. Katzy, E. Yatsenko and K. Werner, EPOS LHC: test of collective hadronization with data measured at the CERN Large Hadron Collider,Phys. Rev.
C 92 (2015) 034906[arXiv:1306.0121] [INSPIRE].
[22] T. Sj¨ostrand and M. van Zijl, A multiple interaction model for the event structure in hadron collisions,Phys. Rev. D 36 (1987) 2019[INSPIRE].
[23] ATLAS collaboration, ATLAS Run 1 PYTHIA8 tunes,ATL-PHYS-PUB-2014-021, CERN, Geneva Switzerland, (2014).
[24] R.D. Ball et al., Parton distributions with LHC data, Nucl. Phys. B 867 (2013) 244
[arXiv:1207.1303] [INSPIRE].
[25] P. Skands, S. Carrazza and J. Rojo, Tuning PYTHIA 8.1: the Monash 2013 tune,Eur. Phys.
J. C 74 (2014) 3024[arXiv:1404.5630] [INSPIRE].
[26] ATLAS collaboration, Charged-particle multiplicities in pp interactions measured with the ATLAS detector at the LHC,New J. Phys. 13 (2011) 053033 [arXiv:1012.5104] [INSPIRE].
[27] ATLAS collaboration, Charged-particle distributions in pp interactions at√s = 8 TeV measured with the ATLAS detector,Eur. Phys. J. C 76 (2016) 403[arXiv:1603.02439] [INSPIRE].
[28] ATLAS collaboration, Further ATLAS tunes of PYTHIA6 and PYTHIA8,
ATL-PHYS-PUB-2011-014, CERN, Geneva Switzerland, (2011).
[29] A.D. Martin, W.J. Stirling, R.S. Thorne and G. Watt, Parton distributions for the LHC,
JHEP03(2017)157
[30] ATLAS collaboration, Measurements of the pseudorapidity dependence of the totaltransverse energy in proton-proton collisions at√s = 7 TeV with ATLAS,JHEP 11 (2012)
033[arXiv:1208.6256] [INSPIRE].
[31] M. Bahr et al., HERWIG++ physics and manual, Eur. Phys. J. C 58 (2008) 639
[arXiv:0803.0883] [INSPIRE].
[32] M. Bahr, J.M. Butterworth and M.H. Seymour, The underlying event and the total cross section from Tevatron to the LHC,JHEP 01 (2009) 065[arXiv:0806.2949] [INSPIRE].
[33] L.A. Harland-Lang, A.D. Martin, P. Motylinski and R.S. Thorne, Parton distributions in the LHC era: MMHT 2014 PDFs,Eur. Phys. J. C 75 (2015) 204[arXiv:1412.3989] [INSPIRE].
[34] M.H. Seymour and A. Siodmok, Constraining MPI models using σeff and recent Tevatron and LHC underlying event data,JHEP 10 (2013) 113[arXiv:1307.5015] [INSPIRE].
[35] ATLAS collaboration, The ATLAS simulation infrastructure,Eur. Phys. J. C 70 (2010) 823
[arXiv:1005.4568] [INSPIRE].
[36] GEANT4 collaboration, S. Agostinelli et al., GEANT4: a simulation toolkit,Nucl. Instrum.
Meth. A 506 (2003) 250[INSPIRE].
[37] ATLAS collaboration, ATLAS insertable B-layer technical design report,
CERN-LHCC-2010-013, CERN, Geneva Switzerland, (2010) [CERN-LHCC-2010-013]
[ATLAS-TDR-19].
[38] ATLAS collaboration, 2015 start-up trigger menu and initial performance assessment of the ATLAS trigger using Run-2 data,ATL-DAQ-PUB-2016-001, CERN, Geneva Switzerland, (2016).
[39] G. Piacquadio, K. Prokofiev and A. Wildauer, Primary vertex reconstruction in the ATLAS experiment at LHC,J. Phys. Conf. Ser. 119 (2008) 032033[INSPIRE].
[40] ALICE collaboration, Strange particle production in proton-proton collisions at√
s = 0.9 TeV with ALICE at the LHC, Eur. Phys. J. C 71 (2011) 1594[arXiv:1012.3257] [INSPIRE].
[41] J.W. Monk and C. Oropeza-Barrera, The HBOM method for unfolding detector effects, Nucl.
Instrum. Meth. A 701 (2013) 17[arXiv:1111.4896] [INSPIRE].
[42] G. D’Agostini, A multidimensional unfolding method based on Bayes’ theorem, Nucl.
Instrum. Meth. A 362 (1995) 487[INSPIRE].
[43] A. Buckley et al., Rivet user manual, Comput. Phys. Commun. 184 (2013) 2803
[arXiv:1003.0694] [INSPIRE].
[44] ATLAS collaboration, ATLAS computing acknowledgements 2016–2017,
