# Inclusive photon production at forward rapidities in pp and p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

S Acharya, D Adamová, A Adler, G Aglieri Rinella, M Agnello, N Agrawal, Z Ahammed, S Ahmad, S.U Ahn, I Ahuja, et al.

► To cite this version:

S Acharya, D Adamová, A Adler, G Aglieri Rinella, M Agnello, et al.. Inclusive photon production at forward rapidities in pp and p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Eur.Phys.J.C, 2023, 83 (7), pp.661. 10.1140/epjc/s10052-023-11729-y. hal-04174608

HAL Id: hal-04174608

https://hal.science/hal-04174608v1

Submitted on 14 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Inclusive photon production at forward rapidities in pp and p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

ALICE Collaboration*

CERN, 1211 Geneva 23, Switzerland

Received: 21 March 2023 / Accepted: 18 June 2023 / Published online: 25 July 2023

© CERN for the benefit of the ALICE collaboration 2023

**Abstract** A study of multiplicity and pseudorapidity distributions of inclusive photons measured in pp and p-Pb collisions at a center-of-mass energy per nucleon-nucleon collision of $\sqrt{s_{NN}}$ = 5.02 TeV using the ALICE detector in the forward pseudorapidity region $2.3 < \eta_{lab} < 3.9$ is presented. Measurements in p-Pb collisions are reported for two beam configurations in which the directions of the proton and lead ion beam were reversed. The pseudorapidity distributions in p-Pb collisions are obtained for seven centrality classes which are defined based on different event activity estimators, i.e., the charged-particle multiplicity measured at midrapidity as well as the energy deposited in a calorimeter at beam rapidity. The inclusive photon multiplicity distributions for both pp and p-Pb collisions are described by double negative binomial distributions. The pseudorapidity distributions of inclusive photons are compared to those of charged particles at midrapidity in pp collisions and for different centrality classes in p-Pb collisions. The results are compared to predictions from various Monte Carlo event generators. None of the generators considered in this paper reproduces the inclusive photon multiplicity distributions in the reported multiplicity range. The pseudorapidity distributions are, however, better described by the same generators.

## 1 Introduction

The primary goal of high-energy heavy-ion physics is the study of a new state of nuclear matter, the quark-gluon plasma (QGP), a thermalized system of partons (quarks and gluons) [1–3]. The study of proton-proton (pp) and proton-nucleus (p-A) collisions provides the baseline for the interpretation of the effects of the QGP formation that are observed in heavy-ion collisions. In addition, the study of p-A collisions helps to understand the effects of cold nuclear matter on the production of final-state particles. Among these effects, an important role is played by transverse momentum ($k_T$) broadening of initial- and final-state partons [4], and

by the modification of the parton distributions functions in bound nucleons compared to those of free nucleons, leading in particular to a reduction of parton densities (shadowing) at small parton fractional momentum $x$ [5]. Interestingly, recent experimental results in pp and p-A collisions particularly at high multiplicities have shown features such as collective flow and strangeness enhancement which are usually attributed to the formation of a QGP in heavy-ion collisions [6–13]. The origin of these phenomena in pp and p-A systems is not yet fully understood and under active scrutiny by the community. It is therefore important to understand the global properties of the system produced in pp and p-A collisions. Multiplicity and pseudorapidity distributions of produced particles are some examples of global observables, providing important information about the particle production mechanisms in these collisions. The integrated yield of particle production is mostly dominated by soft quantum chromodynamics (QCD) interactions, i.e. small momentum transfer ($Q^2$) processes, and is described by non-perturbative phenomenological models. On the other hand, the hard particle production (large $Q^2$) can be described by the well-established theory of perturbative quantum chromodynamics (pQCD). Measurements of multiplicity and pseudorapidity distributions provide important constraints on these models. Moreover, the variation of pseudorapidity ($\eta$) density of produced particles ($dN/d\eta$) with collision centrality can be parametrized by using a two-component model to extract the relative contributions to particle production from hard scatterings and soft processes [14,15].

The measurement of the multiplicity of inclusive photons provides complementary information with respect to those of charged particles as the inclusive photons are mostly produced in the decay of neutral pions ($\pi^0$) [16–18]. A comparative study of charged particles and inclusive photons can shed light on the possible similarities and/or differences in the underlying mechanisms of charged and neutral particle production.

At the Large Hadron Collider (LHC) energies, the underlying mechanisms of particle production could be different

<footnote>* e-mail: alice-publications@cern.ch</footnote>
at central and forward rapidities. Charged-particle multiplic-
ity measurements at the LHC were previously performed by
ALICE, ATLAS, CMS, and LHCb experiments in pp, p–Pb,
and Pb–Pb collisions [19–43]. Inclusive photon production
at forward rapidity was studied in pp collisions at $\sqrt{s} = 0.9$,
2.76, and 7 TeV by ALICE [16]. At midrapidity, the mean
charged-particle multiplicity is found to follow a power law
dependence on $\sqrt{s}$ [20–22] whereas the mean multiplicity
of inclusive photons at forward rapidity can be described by
both a logarithmic and a power law dependence with $\sqrt{s}$ at
LHC energies [16].

In the present study, measurements of multiplicity and pseudorapidity distributions of inclusive photons at forward rapidity are reported for pp, p–Pb, and Pb–p collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with ALICE. The dependence of the inclusive photon production on the centrality of the collision, which is related to the impact parameter of the p–A collision, and its scaling behavior with the number of participating nucleons ($N_{part}$) at forward rapidity is also studied for the first time in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. This measurement enables an extension of our understanding of particle production at forward rapidity.

This paper is organized as follows. The ALICE sub-
detectors relevant for the measurement of inclusive photons
are described in Sect. 2. The data samples used in this anal-
ysis and the selection of events are discussed in Sect. 3. Sec-
tion 4 discusses the details of the used centrality estimators.
A discussion on Monte Carlo (MC) event generators and the
simulation framework is given in Sect. 5. The reconstruction
of inclusive photons is presented in Sect. 6. The correction for
instrumental effects using the unfolding method is described
in Sect. 7. Section 8 discusses the estimation of systematic
uncertainties from various sources. Section 9 presents the
results of inclusive photon multiplicity and pseudorapidity
distributions obtained in this analysis and the outcome of
this study is summarized in Sect. 10.

## 2 Experimental setup

A comprehensive description of the ALICE detectors and
their performances can be found in Refs. [44,45]. The pho-
ton multiplicity detector (PMD) is used for the detection of
inclusive photons at forward rapidity whereas for the purpose
of event selection and centrality determination, the Silicon
Pixel Detector (SPD), the V0 detector and the Zero Degree
Calorimeter (ZDC) are used in this analysis. The ALICE ref-
erence frame is defined with the *z* axis directed along the
beam line, the nominal interaction point (IP) at *z* = 0, and
the positive *z* direction pointing towards the PMD.

The PMD is a preshower detector to measure inclusive
photon multiplicity and its spatial distribution in the forward
rapidity region. It is a gaseous detector placed at *z* = 367 cm

covering the pseudorapidity region of $2.3 < \eta_{\text{lab}} < 3.9$ with
full azimuth coverage. The PMD consists of two fine granu-
lar planes: the charged particle veto plane and the preshower
(PRE) plane. A lead converter with a thickness of three radi-
ation lengths ($3X_0$) is placed between these two planes. The
PMD consists of 184,320 hexagonal cells of size $0.22 \text{ cm}^2$
and depth of $0.5 \text{ cm}$ with a copper honey-comb shaped cath-
ode extended towards a $20 \ \mu\text{m}$ thick gold-plated tungsten
wire at ground potential at the center of each cell. The cells are
arranged in 40 modules in two planes. Each cell is filled with a
gas mixture of Ar and $CO_2$ with a 70:30 ratio. Due to the pres-
ence of the lead converter, the incident photons produce elec-
tromagnetic showers by pair production and bremsstrahlung
radiation and hit several cells in the PRE plane. On the other
hand, charged hadrons generally fire only one or two cells
in the PRE plane and produce a signal typical of a minimum
ionizing particle. After receiving signals in the PRE plane a
photon reconstruction algorithm is applied and a photon-rich
sample of clusters is obtained by applying suitable rejection
criteria to remove clusters from hadrons and from secondary
particles as discussed in Sect. 6. The performance of the PMD
is described in Refs. [16,46,47].

The SPD makes up the two innermost cylindrical layers
of the ALICE Inner Tracking System (ITS) [48] surround-
ing the beam pipe. The two layers cover the pseudorapidity
ranges $|\eta_{\text{lab}}| < 2$ (inner layer) and $|\eta_{\text{lab}}| < 1.4$ (outer layer).
For this analysis, the SPD is mainly used in the determina-
tion of the position of the interaction vertex, the selection
of the minimum bias (MB) trigger (combined with V0), and
the centrality estimation [25–27]. The V0 detector [49,50]
is made of two scintillator arrays placed on either side of
the IP at $z = 330 \text{ cm}$ (VOA, covering the pseudorapidity
interval $2.8 < \eta_{\text{lab}} < 5.1$), and $z = -90 \text{ cm}$ (VOC, covering
$-3.7 < \eta_{\text{lab}} < -1.7$). The V0 is used for event selection and
the determination of collision centrality [25,26,29,32]. The
ZDC [51] consists of two sets of neutron (ZNA and ZNC) and
proton (ZPA and ZPC) calorimeters positioned at $\pm 112.5 \text{ m}$
from the IP, on both sides. The ZPA and ZPC calorimeters
cover the pseudorapidity range $6.5 < |\eta_{\text{lab}}| < 7.4$, while
the ZNA and ZNC calorimeters have a geometric coverage
$|\eta_{\text{lab}}| > 8.8$. The energy detected by the ZN calorimeter on
the Pb-remnant side (ZNA) is used to determine the centrality
in p–Pb collisions [25,26].

## 3 Data sample and event selection

The p–Pb and pp data samples used in this analysis were col-
lected in 2013 during LHC Run 1 and in 2015 during LHC
Run 2, respectively. The p–Pb collisions were recorded for
two beam configurations: in one (denoted as p–Pb), the pro-
ton beam at 4 TeV energy was pointing towards the negative *z*
direction in the ALICE reference system [52], while the lead
ions at 1.58 TeV per nucleon energy moved in the opposite direction; in the other configuration (denoted as Pb–p), the directions of both proton and lead ion beams were reversed. The nucleon–nucleon center-of-mass energy for both configurations resulted in $\sqrt{s_{NN}}$ = 5.02 TeV. Due to the asymmetry of the beam energies, the nucleon–nucleon center-of-mass system in p–Pb collisions is shifted in rapidity by $\Delta y$ = 0.465 with respect to symmetric pp collisions in the direction of the proton beam. In the following, the variable $\eta_{lab}$ is used to indicate the pseudorapidity in the laboratory reference frame. For pp collisions, the center-of-mass frame is same as the laboratory frame.

In pp collisions, a sample of inelastic (INEL) events was selected using a MB trigger condition, which required the detection of at least one particle in either of the two V0 scintillator arrays or in the SPD. For p–Pb and Pb–p data samples, non-single diffractive (NSD) events were selected using a MB trigger condition, which required signals in both V0A and V0C (V0-AND requirement). The background events such as beam–gas or beam–halo interactions occurring outside the interaction region were rejected using the timing information from the V0 detector. Pile-up events, in which two or more collisions are detected as single events, were minimized using the procedure outlined in Ref. [22]. Events were further selected by restricting the reconstructed primary vertex position along the beam axis within $\pm$ 10 cm from the nominal IP. After applying all selection criteria, about 23 M INEL pp collisions, 90 M NSD p–Pb events, and 1.2 M NSD Pb–p collisions are considered for this analysis.

## 4 Centrality determination in p–Pb collisions

In p–Pb collisions, the centrality is determined using two centrality estimators. The first is based on the number of clusters in the outer layer of the SPD (CL1) at midrapidity and the second one makes use of the energy deposited in the neutron calorimeter in the Pb-remnant side at large rapidity (ZNA for the p–Pb beam configuration). Multiplicity-based centrality selection is performed by fitting the multiplicity distribution of the number of SPD clusters with a Glauber MC model [53,54] combined with a simple model for particle production, which assumes that the multiplicity at midrapidity has an average value proportional to $N_{part}$ and a probability distribution described by a negative binomial distribution (NBD). Details about this method can be found in Ref. [25]. However, in p–Pb collisions, multiplicity fluctuations have a large influence on the centrality determination compared to the same effect in heavy-ion collisions as these fluctuations are sizable compared to the width of the $N_{part}$ distribution. Other effects, originating, e.g., from the fragmentation of partons produced in hard scatterings, were found to bias the estimation of the

centrality based on measurements of particle multiplicity at midrapidity [25].

In contrast, the centrality estimation using the energy measured at very large rapidity by the ZNA neutron calorimeter in the Pb-going direction introduces minimal bias. The ZNA detects slow neutrons emitted from nuclear deexcitation processes or knocked out by wounded nucleons [55,56]. The value of $N_{part}$ is calculated using two procedures. The first one is that the energy spectrum obtained from the calorimeter is fitted by a Slow-Nucleon Emission model (which is developed based on the parametrization of results at lower energy p–A experiments [25,55]) coupled to the Glauber MC model. However, the estimation of $N_{part}$ with this method is model-dependent. The second procedure is based on a hybrid method as discussed in detail in Ref. [25]. In this method, the values of $N_{part}$ for each centrality class are obtained under the following assumptions:

$N_{part}^{high-p_T}$: The yield of high-$p_T$ charged particles at midrapidity is proportional to the number of inelastic nucleon-nucleon collisions ($N_{coll}$);

$N_{part}^{Pb-side}$: The charged-particle multiplicity in the Pb-going direction is proportional to the number of wounded target nucleons.

The detailed calculations are discussed in Ref. [25]. The centrality classes used in this analysis are 0–5% (most central (higher multiplicity) collisions), 5–10%, 10–20%, 20–40%, 40–60%, 60–80% and 80–100% (most peripheral (lower multiplicity) collisions) for both CL1 and ZNA estimators.

## 5 Event generators and simulation framework

In the measurements of inclusive photon multiplicity, the corrections for instrumental effects and the evaluation of the systematic uncertainties are performed with the help of MC event generators. For p–Pb collisions, HIJING (v1.36) [57] and DPMJET (v3.0-5) [58] are used whereas PYTHIA 8 (v8.243) with the Monash 2013 tune [59] and EPOS LHC [60] are used for pp collisions. The generated particles from these event generators are transported through the experimental setup using the GEANT 3 software package [61], which includes a detailed description of the apparatus geometry and of the detector response via the AliRoot software framework [62]. The response of the PMD (in terms of the number of fired cells and energy deposition) to the incident particles in simulation is found to reproduce fairly well the results from the test beam data as discussed in Ref. [16].

For pp collisions 1 M (0.5 M) minimum bias INEL events were generated using the PYTHIA 8 (EPOS LHC) event generator whereas for p–Pb and Pb–p collisions 13.7 M (16 M)
and 3.4 M (1 M) NSD events were produced using HIJING (DPMJET).

The PYTHIA 8 [63] event generator is a standard tool for the generation of high energy physics collisions dominantly based on $2 \to 2$ hard scattering processes. In pp interactions, PYTHIA 8 introduces an impact-parameter dependent MultiParton Interaction (MPI) activity [64] to model the soft underlying event (UE). PYTHIA 8 uses the Lund string fragmentation model for the hadronization of partons. An MPI-based model is implemented in PYTHIA 8 to introduce beam remnants and color reconnection wherein all the gluons of lower-$p_T$ interactions are merged with the color-flow dipoles of the highest-$p_T$ one, in such a way that the total string length is minimized [65,66]. In the Monash tune of PYTHIA 8, the parameters relevant for initial-state radiation and MPI are tuned by using the MB, Drell-Yan, and UE data from the Tevatron, SPS and LHC. The EPOS event generator is a parton-based MC model with flux tube initial condition for hadron–hadron collisions. It uses Gribov–Regge theory for describing soft interactions. The EPOS LHC [60] generator is tuned to LHC data to describe the results from various collision systems at different center-of-mass energies, particularly the observed collective behavior in pp and p–Pb collisions at LHC energies.

The HIJING event generator is a pQCD inspired MC model aimed particularly at the study of jet and minijet production and jet-medium interactions in heavy-ion collisions. It also allows the study of particle production in high energy hadron–hadron and hadron–nucleus collisions. In p–A and heavy-ion collisions, multiple nucleon–nucleon interactions are simulated using binary approximation and the Glauber model. Nuclear shadowing effects are taken into account using parameterized parton distribution functions inside the nucleus. The soft interactions in HIJING are described by Lund FRITIOF [67] and the dual parton model (DPM) [68]. For jet fragmentation and hadronization, HIJING uses the prescription described in the Lund fragmentation model. The DPMJET event generator is a multi-purpose MC model based on the DPM capable of simulating hadron–hadron, hadron–nucleus, nucleus–nucleus, photon–hadron, photon–photon and photon–nucleus interactions from a few GeV up to the highest cosmic-ray energies. It uses the Glauber–Gribov multiple scattering formalisms to calculate nuclear cross sections and utilizes the Reggeon theory and pQCD to describe soft interactions and hard interactions respectively.

## 6 Photon reconstruction

A procedure similar to that described in Ref. [16] is adopted to reconstruct photons incident on the preshower plane of the PMD. An incident photon is expected to deposit a larger amount of energy (measured in ADC units) and produce

a signal in a larger number of cells compared to an incident charged hadron on the PRE plane of the PMD. Using a nearest-neighbor clustering algorithm, contiguous cells having non-zero energy deposition are grouped together forming clusters. Each cluster carries information, such as the number of fired cells contained in the cluster, the position of the centroid of the cluster, and the total energy deposition in the cluster.

Suitable thresholds on the number of fired cells and the energy deposited in clusters are applied to discriminate between photon and hadron clusters. The number of clusters that satisfy the discrimination threshold conditions is termed as $N_{\gamma}$-like clusters. A similar condition in discrimination threshold as used in Ref. [16], namely number of cells greater than 2 and energy deposition greater than 9 times the MPV (Most Probable Value of the charged pion ADC distribution from a Landau distribution function fitted to the data), is applied to obtain the distributions of $N_{\gamma}$-like clusters. The $N_{\gamma}$-like clusters provide a photon-rich sample, however, there are contaminations due to residual background from hadron clusters and from secondary particles (i.e. particles produced in interactions with the upstream material in front of the PMD) that reduce the purity of the sample. The application of the rejection criteria also removes some of the actual photons from the sample. The distributions of $N_{\gamma}$-like clusters are therefore required to be corrected for these effects as discussed in the next section.

## 7 Correction procedure

The measured distributions of $N_{\gamma}$-like clusters are distorted by several instrumental effects such as photon detection inefficiency, limited acceptance, contaminations from hadrons and secondary particles etc. A set of correction procedures is adopted to obtain the corrected $N_{\gamma}$ distributions.

These instrumental effects are corrected for via a response matrix $R_{mt}$ extracted from MC simulations using PYTHIA 8 with the Monash 2013 tune to generate pp collisions and the HIJING event generator for p–Pb collisions. These simulations include a detailed description of the experimental conditions and the detector settings and the simulated events are reconstructed and analyzed with the same algorithms used for the experimental data. Inelastic pp collisions can be classified into two processes: diffractive and non-diffractive (ND). The ND interactions are the dominant processes in pp collisions. According to Regge theory [69], diffractive scattering occurs via the exchange of Pomeron. In single diffractive (SD) processes, the exchanged Pomeron interacts and one of the protons breaks up, producing particles of the diffractive system of mass $M_X$. In double diffractive (DD) processes, both protons break up. It was found that the $M_X$ distribution for the SD processes in experimental data [70,71] differs from the
<figure><img src="image_2.png" /><figcaption>**Fig. 1** Graphical representation of the detector response matrices for three different collision systems: pp (a), p–Pb (b) and Pb–p (c). The reconstructed and true photon multiplicities are denoted as $N_{\gamma}$-like and</figcaption></figure>

$N_{\gamma}$-true, respectively, obtained from simulations using PYTHIA 8 with the Monash 2013 tune for pp collisions and HIJING for p–Pb and Pb–p collisions

one predicted by PYTHIA 8 simulations [72]. To account for this effect in simulation, diffraction-tuned event generators (which can reproduce diffraction cross sections and the shapes of the $M_X$ distribution obtained from data) were used in previous charged-particle measurements by ALICE [22,23]. In this work, the $M_X$ distribution of SD events in PYTHIA 8 with the Monash 2013 tune is reweighted to match the diffraction-tuned PYTHIA simulation used in previous 7 TeV analysis [22]. The response matrix in pp collisions is constructed using this reweighted PYTHIA 8 simulation.

Figure 1 shows the graphical representation of the detector response matrices for minimum bias pp, p–Pb, and Pb–p collisions in terms of the correlation between the true photon multiplicity ($N_{\gamma}$-true) from the event generator and the reconstructed photon multiplicity ($N_{\gamma}$-like) in the simulated events. The mean ($(N_{\gamma}$-like)) and width ($\sigma_{N_{\gamma}$-like})$ of the distribution of $N_{\gamma}$-like are presented in Figs. 12 and 13 as a function of $N_{\gamma}$-true for pp and p–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The resolution of photon multiplicity reconstruction in pp (p–Pb) collisions at $N_{\gamma}$-true = 44(140) is found to be ~12% (~6%). Simulation studies have shown that the PMD is sensitive to transverse momenta ($p_T$) as low as ~50 MeV/c [16]. In the

correction procedure, true photons of all $p_T$ are considered, which makes the present photon measurement inclusive. In simulation, the fraction of inclusive photons having $p_T$ below 50 MeV/c is estimated to be about 16% for both pp and p–Pb collisions. The detection efficiency of inclusive photons as a function of $\eta$ varies from 27% to 52%. The matrix element $R_{mt}$ represents the conditional probability that an event with true multiplicity $t$ is measured as an event with multiplicity $m$. The distorted measured distribution ($M$) can be expressed as

$$M = R_{mt}T. \qquad (1)$$

One can therefore recover the true distribution ($T$) for given $M$ as

$$T = R_{mt}^{-1} M. \qquad (2)$$

However, the matrix $R_{mt}$ may be singular and one can calculate inverse matrix $R_{mt}^{-1}$ only if $R_{mt}$ is not singular. Furthermore, even if $R_{mt}$ can be inverted, the results obtained with Eq. 2 contain oscillations mainly because of finite statistics in the response matrix. To overcome this problem, a regularized unfolding method based on Bayes' theorem [73] using the
<figure><img src="image_3.png" /><figcaption>**Fig. 2** MC closure test of the unfolding method for pp (a), p–Pb (b), and Pb–p (c) collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. The reconstructed, unfolded and true photon multiplicity distributions are presented in the top panels. The ratios of unfolded to true multiplicity distributions are shown in the bottom panels</figcaption></figure>

RooUnfold software package [74] is used to correct the measured $N_{\gamma}$-like distributions following the description outlined in Refs. [16,20–23].

In order to optimize the parameters (number of iterations) of the Bayesian unfolding method, a MC closure test was performed using simulations. The simulated data were divided

into two statistically independent samples for the closure test. The response matrix was built using one set whereas the other set was used to obtain the reconstructed ($N_{\gamma}$-like) distributions, correct them using Eq. 2, and compare the result to the true distributions. The sensitivity of the results was checked by varying the number of iterations and a value corresponding to the best performance of the closure test was chosen as the optimized number of iterations. Figure 2 shows the performance of the closure tests for pp (Fig. 2(a)), p–Pb (Fig. 2(b)), and Pb–p (Fig. 2(c)) collisions. The upper panels show the true, reconstructed ($N_{\gamma}$-like), and unfolded photon multiplicity distributions. In the bottom panels, the ratios of the unfolded to the true photon multiplicity distributions are plotted for the optimized number of iterations (1 for pp, 4 for p–Pb, and 2 for Pb–p collisions). The ratio plots indicate that there is good agreement between the unfolded and true distributions except at high multiplicities where fluc-

<table><caption>**Table 1** Trigger and vertex reconstruction efficiencies for INEL events in pp collisions and for NSD events in p–Pb and Pb–p collisions. Uncertainties listed are total uncertainties. Statistical uncertainties are negligible</caption><thead><tr><th>System</th><th>Trigger efficiency (%)</th><th>Vertex reconstruction efficiency (%)</th></tr></thead><tbody><tr><td>pp</td><td>94.5 ± 1.2</td><td>95.0 ± 0.6</td></tr><tr><td>p–Pb</td><td>98.5 ± 0.9</td><td>99.8 ± 0.1</td></tr><tr><td>Pb–p</td><td>98.5 ± 0.9</td><td>99.8 ± 0.1</td></tr></tbody></table>
tuations arise due to the limited size of the simulated sample. The measured photon distributions in experimental data were unfolded using the optimized number of iterations. An alternate method of unfolding, singular value decomposition (SVD) [75], was also used and consistent results were obtained. In this method, the response matrix can be factorized in the form $R_{mt} = USV^T$, where U and V are orthogonal matrices and S is a diagonal matrix with non-negative diagonal elements. These diagonal elements are called singular values of the matrix $R_{mt}$. Regularization is applied by using a smooth cut-off on small singular-value contributions suppressing large fluctuations in the unfolded distribution. The same procedure was applied to extract the pseudorapidity distributions, but the unfolding of photon multiplicity distributions was performed separately in eight pseudorapidity intervals of a width of 0.2 units.

The unfolded results in MB collisions are then corrected for the trigger ($\epsilon_{Trig}$) and vertex reconstruction ($\epsilon_{Vtx}$) efficiencies. A similar procedure as discussed in [16] is adopted to estimate $\epsilon_{Trig}$ and $\epsilon_{Vtx}$ using MC simulations. These efficiencies can be defined as

$$ \epsilon_{\text{Trig}} = \frac{N_{\text{Trig}}}{N_{\text{All}}}, \quad \epsilon_{\text{Vtx}} = \frac{N_{\text{TrigVtx}}}{N_{\text{Trig}}} \qquad (3) $$

where $N_{All}$ and $N_{Trig}$ correspond to the number of all simulated events and of triggered events, respectively, and $N_{TrigVtx}$ corresponds to triggered events with reconstructed vertex. The values of $\epsilon_{Trig}$ and $\epsilon_{Vtx}$ for pp, p–Pb, and Pb–p collisions are reported in Table 1. The unfolded pseudorapidity distributions are corrected for these efficiencies using an overall correction factor whereas multiplicity dependent correction factors are used to correct the multiplicity distributions.

The V0-AND trigger used to select p–Pb events is not fully efficient for NSD events. The correction for the V0-AND trigger inefficiency and the imperfect description in the MC of the vertex reconstruction efficiency would mainly concern the most peripheral event class (80–100%). In this work, centrality classes have been defined as percentiles of the visible cross section and the centrality dependent pseudorapidity distributions are not corrected for trigger inefficiencies as was also reported in Ref. [25].

## 8 Estimation of systematic uncertainties

The contributions of the different sources of systematic uncertainties associated with the measurements of multiplicity and pseudorapidity distributions of inclusive photons are summarized in Table 2 for the measurements in MB events of pp, p–Pb, and Pb–p collisions and in Tables 3 and 4 for the measurements in centrality classes of p–Pb collisions. The results of p–Pb collisions in centrality classes are not affected by the uncertainties on the event selection and diffraction mass. The total systematic uncertainties are calculated by adding the contributions from the individual sources in quadrature. The methods used for the estimation of the systematic uncertainties are discussed in the following subsections.

### 8.1 Effect of upstream material

One of the major contributions to the systematic uncertainties comes from the uncertainty in the implementation of upstream material in front of the PMD (i.e., the material between the nominal IP and the PMD) in the GEANT 3

<caption>**Table 2** Contributions to systematic uncertainties (in percent) in the measurements of pseudorapidity and multiplicity distributions of inclusive photons in pp, p–Pb, and Pb–p collisions. For multiplicity distributions,</caption>

numbers are given at multiplicity values of 0, the mean $\langle m \rangle$ and when $P(N_{\gamma}) = 10^{-3}$

<table><thead><tr><th rowspan="3">Sources</th><th colspan="3">dN<sub>γ</sub>/dη<sub>lab</sub> analysis</th><th colspan="9">P(N<sub>γ</sub>) analysis</th></tr><tr><th rowspan="2">pp</th><th rowspan="2">p–Pb</th><th rowspan="2">Pb–p</th><th colspan="3">pp</th><th colspan="3">p–Pb</th><th colspan="3">Pb–p</th></tr><tr><th>N<sub>γ</sub><br>0</th><th>N<sub>γ</sub><br>⟨m⟩</th><th>P(N<sub>γ</sub>)<br>10<sup>-3</sup></th><th>N<sub>γ</sub><br>0</th><th>N<sub>γ</sub><br>⟨m⟩</th><th>P(N<sub>γ</sub>)<br>10<sup>-3</sup></th><th>N<sub>γ</sub><br>0</th><th>N<sub>γ</sub><br>⟨m⟩</th><th>P(N<sub>γ</sub>)<br>10<sup>-3</sup></th></tr></thead><tbody><tr><td>Upstream material</td><td>7.0</td><td>8.2–10.1</td><td>8.7–9.9</td><td>3.0</td><td>2.5</td><td>20.0</td><td>15.2</td><td>3.4</td><td>34.0</td><td>18.5</td><td>3.7</td><td>37.2</td></tr><tr><td>Hadron and secondary contamination</td><td>1.9</td><td>1.9</td><td>2.0</td><td>1.4</td><td>1.4</td><td>7.5</td><td>6.7</td><td>1.2</td><td>1.2</td><td>0.7</td><td>1.1</td><td>7.2</td></tr><tr><td>Event generator for response matrix</td><td>2.5</td><td>3.2</td><td>1.0</td><td>2.6</td><td>3.0</td><td>5.0</td><td>3.2</td><td>3.6</td><td>16.0</td><td>3.0</td><td>3.5</td><td>14.4</td></tr><tr><td>Unfolding method</td><td>3.1</td><td>0.6</td><td>1.1</td><td>25.0</td><td>3.7</td><td>23.0</td><td>11.0</td><td>1.2</td><td>1.2</td><td>11.0</td><td>3.0</td><td>3.0</td></tr><tr><td>Event selection efficiency</td><td>1.2</td><td>0.9</td><td>0.9</td><td>19.0</td><td>1.6</td><td>1.6</td><td>33.5</td><td>1.2</td><td>0.7</td><td>34.0</td><td>1.0</td><td>0.8</td></tr><tr><td>Diffraction shape</td><td>2.1</td><td>n/a</td><td>n/a</td><td>1.3</td><td>4.0</td><td>5.3</td><td>n/a</td><td>n/a</td><td>n/a</td><td>n/a</td><td>n/a</td><td>n/a</td></tr><tr><td>Diffraction ratio</td><td>2.3</td><td>n/a</td><td>n/a</td><td>7.4</td><td>4.6</td><td>0.2</td><td>n/a</td><td>n/a</td><td>n/a</td><td>n/a</td><td>n/a</td><td>n/a</td></tr><tr><td>Total</td><td>8.9</td><td>9.1–10.8</td><td>9.1–10.2</td><td>32.8</td><td>8.4</td><td>32.0</td><td>39.3</td><td>5.4</td><td>37.6</td><td>40.3</td><td>6.0</td><td>40.7</td></tr></tbody></table>
<table><caption>**Table 3** Contributions to systematic uncertainties (in percent) in the measurements of centrality dependent pseudorapidity distributions of inclusive photons using the CL1 estimator in p–Pb collisions</caption><thead><tr><th>Sources</th><th>0–5%</th><th>5–10%</th><th>10–20%</th><th>20–40%</th><th>40–60%</th><th>60–80%</th><th>80–100%</th></tr></thead><tbody><tr><td>Upstream material</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td></tr><tr><td>Hadron and secondary contamination</td><td>2.1</td><td>2.3</td><td>2.1</td><td>2.3</td><td>2.4</td><td>2.1</td><td>3.1</td></tr><tr><td>Event generator for response matrix</td><td>3.0</td><td>3.0</td><td>3.1</td><td>3.0</td><td>2.9</td><td>3.1</td><td>3.0</td></tr><tr><td>Unfolding method</td><td>1.7</td><td>2.9</td><td>2.6</td><td>2.9</td><td>3.1</td><td>2.8</td><td>2.8</td></tr><tr><td>Total</td><td>9.2–10.9</td><td>9.5–11.2</td><td>9.4–11.0</td><td>9.5–11.1</td><td>9.6–11.2</td><td>9.5–11.1</td><td>9.7–11.3</td></tr></tbody></table>

<table><caption>**Table 4** Contributions to systematic uncertainties (in percent) in the measurements of centrality dependent pseudorapidity distributions of inclusive photons using the ZNA estimator in p–Pb collisions</caption><thead><tr><th>Sources</th><th>0–5%</th><th>5–10%</th><th>10–20%</th><th>20–40%</th><th>40–60%</th><th>60–80%</th><th>80–100%</th></tr></thead><tbody><tr><td>Upstream material</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td><td>8.2–10.1</td></tr><tr><td>Hadron and secondary contamination</td><td>2.1</td><td>1.3</td><td>1.6</td><td>1.6</td><td>2.1</td><td>1.7</td><td>2.6</td></tr><tr><td>Event generator for response matrix</td><td>1.4</td><td>1.8</td><td>2.2</td><td>3.2</td><td>3.8</td><td>3.3</td><td>2.7</td></tr><tr><td>Unfolding method</td><td>2.8</td><td>1.9</td><td>2.0</td><td>1.7</td><td>2.7</td><td>2.8</td><td>3.0</td></tr><tr><td>Total</td><td>9.1–10.8</td><td>8.7–10.5</td><td>9.0–10.7</td><td>9.2–10.9</td><td>9.7–11.3</td><td>9.5–11.1</td><td>9.6–11.2</td></tr></tbody></table>

description of the ALICE apparatus. At forward rapidity, the material budget was studied in detail for the charged-particle multiplicity measurements [23,30,31]. Based on this study, an upper limit of the material budget uncertainty of 10% was considered in this analysis. To study the effect of this on photon counting in p–Pb and Pb–p collisions, two response matrices using the HIJING event generator were considered, one with the default material description and the other with a 10% increase in the overall material of the ALICE apparatus. The measured distribution in data is then unfolded using these two different matrices. The difference between the two unfolded multiplicities determines the systematic uncertainty due to the material in front of the PMD. A similar procedure is applied for each $\eta_{lab}$ interval to obtain the systematic uncertainties for the pseudorapidity distributions of inclusive photons in MB events and in all centrality classes. The systematic uncertainty due to material effects for pp collisions is taken from the previously published paper on inclusive photon measurement in pp collisions at $\sqrt{s} = 0.9, 2.76$, and 7 TeV [16].

## 8.2 Contamination in the photon sample

The photon–hadron discrimination conditions in terms of energy deposition and number of cells in clusters were optimized to minimize contaminations from hadronic clusters and secondary particles. The purity of the photon sample with the default threshold condition (energy deposition greater than 9 times the MPV and number of cells greater than 2) is found to be 65%. To estimate the systematic uncertainty

due to the correction for the contamination, the analysis was repeated with a different threshold condition on energy deposition (greater than 6 times the MPV). The unfolded distributions are obtained for both cases and the difference between these two distributions is quoted as the systematic uncertainty. By applying a higher threshold condition on energy deposition (greater than 12 times the MPV), the estimated systematic uncertainty with respect to the default threshold condition is found to be of the same order as the quoted one.

## 8.3 Unfolding procedure

As discussed in Sect. 7, we rely on MC models to unfold the measured photon multiplicity and pseudorapidity distributions. The sensitivity of the unfolded results to the choice of the MC models used in the unfolding procedure is estimated by comparing the results obtained using the response matrices from two different event generators. For pp collisions, two separate response matrices are built using PYTHIA 8 with the Monash 2013 tune and EPOS LHC, whereas for p–Pb and Pb–p collisions those are constructed using HIJING and DPMJET. These matrices are used to unfold the measured distributions in data. The difference in the unfolded distributions is considered as the systematic uncertainty. The effect of the choice of the unfolding method was determined by using an alternative unfolding method, SVD, in addition to the default Bayesian method. The difference in the unfolded results is considered as the systematic uncertainty.
## 8.4 Event selection efficiency

The systematic uncertainty due to possible imperfections in the correction for the event selection efficiency is determined by estimating the trigger and vertex reconstruction efficiencies using two different event generators. For pp collisions, the efficiency values are computed from events simulated with the PYTHIA 8 with Monash 2013 tune and EPOS LHC generators, while HIJING and DPMJET are used to obtain these efficiencies in p–Pb and Pb–p collisions. This uncertainty mostly influences the zero multiplicity bin in the multiplicity distributions and it reduces significantly at large multiplicity because single and double diffraction contributions become smaller when going to higher multiplicity.

## 8.5 Diffraction mass distributions

As described in Sect. 7, the central values of the final results in pp collisions are determined using simulations with the PYTHIA 8 event generator with the Monash 2013 tune reweighted to match the measured $M_X$ distribution in SD events. The uncertainty associated with this procedure is estimated by taking the difference between the results obtained with the reweighted mass distribution and the default one from PYTHIA 8 with Monash 2013 tune. The uncertainty related to the fractions of SD and DD events produced by PYTHIA 8 is also evaluated by varying those fractions by ± 30% of their nominal values.

# 9 Results and discussions

The multiplicity and pseudorapidity distributions of inclusive photons in the forward pseudorapidity interval $2.3 < \eta_{\text{lab}} < 3.9$ measured using the PMD for pp collisions at $\sqrt{s} = 5.02$ TeV and p–Pb and Pb–p collisions at $\sqrt{s_{NN}} = 5.02$ TeV are presented and discussed in this section. The distributions are obtained without any selection on the transverse momentum of the photons. The obtained results are compared with predictions from MC models and available experimental measurements of charged particles.

## 9.1 Multiplicity distributions

The measured multiplicity distributions of inclusive photons are presented as probability distributions ($P(N_\gamma)$) as a function of $N_\gamma$ in Fig. 3. The systematic uncertainty is represented by the gray bands and the statistical uncertainty is smaller than the marker size. The measurements are obtained for pp and for both configurations of p–Pb collisions in the pseudorapidity interval $2.3 < \eta_{\text{lab}} < 3.9$. The average photon multiplicity for pp collisions is $6.44 \pm 0.36$ (sys), for p–Pb collisions is $29.38 \pm 1.88$ (sys), and for Pb–p collisions is

$24.59 \pm 1.52$ (sys). The multiplicity range in p–Pb collisions reaches up to 140 which is about 5 times larger than its mean multiplicity. The average multiplicity in pp collisions is about 4.5 times smaller than that in p–Pb collisions. The different multiplicity distributions and average values in p–Pb and Pb–p collisions arise from the different particle production in the p-going and Pb-going hemispheres as well as from the difference in the rapidity coverage of the PMD in the center-of-mass frame because of the $\Delta\gamma$ shift of $\pm 0.465$ units between p–Pb and Pb–p collisions.

The results in pp collisions (Fig. 3(a)) are compared to the predictions from PYTHIA 8 with the Monash 2013 tune [59], PYTHIA 6 (v6.425) with the Perugia 2011 tune [76], PHOJET (v1.12) [77], and EPOS LHC [60] event generators. The ratios between MC predictions and data are shown in the bottom panels. It is observed that PYTHIA 8 with the Monash 2013 tune and EPOS LHC are unable to reproduce the inclusive photon production in pp collisions at $\sqrt{s} = 5.02$ TeV at high multiplicities. On the other hand, the Perugia 2011 tune of PYTHIA 6 and PHOJET are found to fairly describe the data within uncertainties.

The multiplicity distributions in p–Pb collisions (Fig. 3(b)) are compared to predictions obtained from HIJING, DPMJET, and AMPT (v2.25) [78]. None of the models could reproduce the shape of the distributions in the full multiplicity range. HIJING is closer to the data points at intermediate multiplicities. The tail of the distribution at high multiplicity can be described by both HIJING and AMPT models within uncertainties. For Pb–p collisions (Fig. 3(c)), DPMJET strongly underestimates the data at high multiplicity. HIJING describes the data slightly better in the full range of multiplicity compared to AMPT.

These experimental results are expected to provide new constraints on inclusive photon production mechanisms implemented in theoretical models.

In Fig. 4, the measured multiplicity distributions are fitted by a weighted sum of two negative binomial distributions (NBDs) (Eq. 4) to extract the relative contributions from soft and semihard processes in the particle production mechanisms. In Eq. 4, the subscripts and superscripts in the parameters indicate the respective components. The $\alpha_{\text{soft}}$ parameter gives the fraction of soft events. The NBD distribution, $P_{\text{NBD}}$, has the following parameters: the average multiplicity $\langle n \rangle$, which is found to increase with increasing $\sqrt{s}$, and the shape parameter, $k$, which decreases with increasing $\sqrt{s}$ in hadronic collisions [16,22,23].

$$P(n) = \lambda[\alpha_{\text{soft}}P_{\text{NBD}}(n, \langle n^{\text{soft}} \rangle, k^{\text{soft}}) + (1 - \alpha_{\text{soft}})P_{\text{NBD}}(n, \langle n^{\text{semihard}} \rangle, k^{\text{semihard}})] \quad (4)$$

where

$$P_{\text{NBD}}(n, \langle n \rangle, k) = \frac{\Gamma(n + k)}{\Gamma(k)\Gamma(n + 1)} \frac{(\langle n \rangle/k)^n}{(1 + \langle n \rangle/k)^{n+k}} \quad (5)$$
<figure><img src="image_4.png" /><figcaption>**Fig. 3** Top panels: Inclusive photon multiplicity distributions measured within $2.3 < \eta_{\text{lab}} < 3.9$ in pp (a), p–Pb (b), and Pb–p (c) collisions at $\sqrt{s_{\text{NN}}}$ = 5.02 TeV. Results from various MC predictions are superimposed. Bottom panels: The ratios between MC results and data are shown. Shaded boxes represent the systematic uncertainties</figcaption></figure>

Double NBDs do not describe the value $P(0)$, therefore, the bin $N_{\gamma} = 0$ was excluded from the fit and a normalization factor $\lambda$ was introduced to account for this. The systematic uncertainties associated with (a) the change of photon-hadron discrimination thresholds, (b) the tuning of event generator for diffraction, (c) the effect of upstream material, and (d) the correction of event selection efficiency produce a correlated shift of the multiplicity distributions. These correlated uncertainties are therefore not considered in the fitting procedure. The other sources of systematic uncertainties (uncorrelated) are taken into account in the fit. The obtained parameters from these fits are printed in bold in Table 5. The other set of parameters reported in Table 5 is obtained from the fitting of the data by considering both correlated and uncorrelated uncertainties to provide an estimate of how much the

fit parameters change due to the presence of correlations in the systematic uncertainties.

The double NBD fit describes the data within uncertainties. The results show an increase in the mean multiplicity of inclusive photons with system size whereas the shape parameter decreases with system size. Moreover, it is observed that $\langle n^{\text{semihard}} \rangle \approx 3\langle n^{\text{soft}} \rangle$ which is consistent with the results of charged-particle production at central and forward rapidities in pp collisions [22,23].

## 9.2 Pseudorapidity distributions in NSD and INEL events

The pseudorapidity distributions of inclusive photons measured within $2.3 < \eta_{\text{lab}} < 3.9$ in pp, p–Pb, and Pb–p collisions are presented in Fig. 5. The gray bands represent
<figure><img src="image_5.png" /><figcaption>**Table 5** Double NBD fit parameters for the inclusive photon multiplicity distributions. The obtained parameters from fits to the data by considering only uncorrelated uncertainties are printed in bold. The other</figcaption></figure>

set of parameters is obtained when both correlated and uncorrelated uncertainties are considered in the fitting procedure

<table><thead><tr><th>Collision system</th><th>λ</th><th>α<sub>soft</sub></th><th>k<sup>soft</sup></th><th>⟨n<sup>soft</sup>⟩</th><th>k<sup>semihard</sup></th><th>⟨n<sup>semihard</sup>⟩</th></tr></thead><tbody><tr><td rowspan="2">pp</td><td><strong>0.79 ± 0.03</strong></td><td><strong>0.60 ± 0.19</strong></td><td><strong>2.97 ± 1.16</strong></td><td><strong>4.65 ± 1.01</strong></td><td><strong>4.56 ± 1.79</strong></td><td><strong>13.47 ± 2.60</strong></td></tr><tr><td>0.79 ± 0.04</td><td>0.61 ± 0.28</td><td>2.97 ± 1.74</td><td>4.67 ± 1.60</td><td>4.64 ± 2.76</td><td>13.55 ± 4.07</td></tr><tr><td rowspan="2">Pb–p</td><td><strong>0.99 ± 0.01</strong></td><td><strong>0.29 ± 0.10</strong></td><td><strong>1.65 ± 0.30</strong></td><td><strong>8.92 ± 2.10</strong></td><td><strong>4.09 ± 0.39</strong></td><td><strong>31.08 ± 1.31</strong></td></tr><tr><td>0.99 ± 0.01</td><td>0.31 ± 0.13</td><td>1.60 ± 0.33</td><td>9.70 ± 4.00</td><td>4.22 ± 0.65</td><td>31.56 ± 2.14</td></tr><tr><td rowspan="2">p–Pb</td><td><strong>0.99 ± 0.01</strong></td><td><strong>0.62 ± 0.10</strong></td><td><strong>1.10 ± 0.05</strong></td><td><strong>18.76 ± 2.50</strong></td><td><strong>3.89 ± 0.30</strong></td><td><strong>47.21 ± 1.28</strong></td></tr><tr><td>0.98 ± 0.01</td><td>0.60 ± 0.11</td><td>1.10 ± 0.10</td><td>18.80 ± 3.60</td><td>3.64 ± 0.40</td><td>46.18 ± 1.72</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Fig. 5** Top panels: pseudorapidity distribution of inclusive photons measured within $2.3 < \eta_{\text{lab}} < 3.9$ in pp (a), p–Pb (b), and Pb–p (c) collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Results from various MC predictions are superimposed. Bottom panels: the ratios between MC results and data are shown. Shaded boxes represent the systematic uncertainties</figcaption></figure>

the systematic uncertainty and the statistical uncertainty is smaller than the marker size. The measurements are compared to the predictions from MC event generators and the ratios between MC results and data are shown in the bottom panels.

In pp collisions (Fig. 5(a)), PHOJET and PYTHIA 6 with the Perugia 2011 tune models show better agreement with the data compared to EPOS LHC and PYTHIA 8 with the Monash 2013 tune. All model predictions are however found to be compatible within about 20% from the measured values, but the $\eta_{\text{lab}}$ dependence is flatter in the data than in the models.

In p–Pb collisions the value of pseudorapidity density of inclusive photons is found to be slightly larger than that for Pb–p collisions as expected due to the fact that the PMD measures particles produced in the Pb-going direction in the

former configuration and in the p-going direction in the latter, and also due to the different rapidity ranges covered by the PMD in the center-of-mass frame for p–Pb and Pb–p collisions as discussed in Sect. 9.1. The data are compared to several models with different descriptions of particle production, all shifted by $\Delta y = 0.465$ to take into account the shift to the laboratory system. It is observed that HIJING describes both p–Pb and Pb–p results within systematic uncertainties. AMPT overpredicts the data in the Pb-going side and reproduces the measurements within uncertainties in the p-going side. DPMJET underestimates the multiplicity in the p-going side. One can observe that the DPMJET curve lies in between AMPT and HIJING in the Pb-going side and goes below in the p-going side. This suggests that DPMJET predicts a slightly narrower distribution. It is also noted that all models lie within about 15% from the data points.
**Fig. 6 (a)** Comparison of $dN_γ/dη_{lab}$ measured at forward rapidity in INEL pp collisions for various center-of-mass energies, (b) average photon multiplicity within $2.3 < η_{lab} < 3.9$ as a function of collision energy in INEL pp collisions. The results at 0.9, 2.76, and 7 TeV are taken from Ref. [16]

<figure><img src="image_7.png" /><figcaption>**Table 6** The average photon multiplicity $\langle N_γ \rangle$ within $2.3 < η_{lab} < 3.9$ for various center-of-mass energies in INEL pp collisions. The quoted errors are systematic uncertainties. Statistical uncertainties are negligible. Data points at 0.9, 2.76, and 7 TeV are taken from Ref. [16]</figcaption></figure>

<table><thead><tr><th>$\sqrt{s}$ (TeV)</th><th>$\langle N_γ \rangle$</th></tr></thead><tbody><tr><td>0.9</td><td>4.09 ± 0.20</td></tr><tr><td>2.76</td><td>5.82 ± 0.26</td></tr><tr><td>5.02</td><td>6.44 ± 0.36</td></tr><tr><td>7</td><td>6.81 ± 0.35</td></tr></tbody></table>

The same models that do not reproduce the high-multiplicity part of the $P(N_γ)$ distribution in MB collisions, provide a fair description of the multiplicity of inclusive photons, $dN_γ/dη_{lab}$, at forward rapidities in pp, p-Pb, and Pb-p collisions.

The simultaneous comparison of the measured $P(N_γ)$ and $dN_γ/dη_{lab}$ distributions to the predictions of different event generators has therefore the potential of further constraining the models of inclusive photon production in high energy hadronic collisions.

The $dN_γ/dη_{lab}$ distributions measured in the rapidity interval $2.3 < η_{lab} < 3.9$ in INEL pp collisions at various center-of-mass energies are compared in Fig. 6(a). Two additional systematic uncertainties (diffraction shape and diffraction ratio) are considered in this analysis and they lead to a total systematic uncertainty that is slightly larger at $\sqrt{s} = 5.02$ TeV with respect to the other energies. However, the total systematic uncertainty at all energies mostly comes from the 10% uncertainty on the upstream material in front of the PMD. The material uncertainty shows essentially no $\sqrt{s}$ variation and is therefore fully correlated over all energies. A smooth increase of the inclusive photon multiplicity with increasing collision energy is observed. The average photon multiplicity $\langle N_γ \rangle$ within $2.3 < η_{lab} < 3.9$ as a function of $\sqrt{s}$ is shown in Fig. 6(b). The results of the measurements of $\langle N_γ \rangle$ are given in Table 6.

Since the dominant contribution to inclusive photon production comes from $\pi^0$ decays, the number of produced photons should be similar to the charged-particle multiplicity.

The comparison between the pseudorapidity distribution of inclusive photons in INEL pp and NSD p-Pb (combined with results of Pb-p mirrored with respect to $η_{lab} = 0$) events measured at forward rapidity and that of charged particles at midrapidity by ALICE [24,79] and CMS [41] are presented in Fig. 7. It is observed that the inclusive photon production at forward rapidity follows the trend of charged-particle production at midrapidity. The predictions from the different event generators are also displayed in Fig. 7 and they show similar values for charged-particle (dashed lines) and photon (solid lines) multiplicity at forward and backward pseudorapidity, while at midrapidity the photon and charged-particle pseudorapidity density differ. The origin of this difference is due to a mass effect in the transformation between $dN/dη$ and $dN/dη$ at $η \approx 0$. The measured $dN_{ch}/dη_{lab}$ in p-Pb collisions is well described by both HIJING and DPMJET event generators. HIJING reproduces better the $dN_γ/dη_{lab}$ compared to DPMJET. For pp collisions, both PYTHIA 8 with the Monash 2013 tune and EPOS LHC overpredict the photon and charged-particle multiplicity.

## 9.3 Pseudorapidity distributions in centrality classes of p-Pb collisions

The pseudorapidity distributions of inclusive photons as a function of pseudorapidity measured in p-Pb collisions are presented in Fig. 8 for several centrality classes estimated using two centrality estimators, CL1 (top) and ZNA (bottom). The analysis was not performed for the Pb-p configuration, hence only one rapidity interval is covered for the photon measurements, namely the one in the Pb-going direction, and the ZNA measures the energy of the neutrons emitted in the direction of the Pb beam. The multiplicity in the most central (0-5%) collisions when considering the CL1 (ZNA) estimator reaches ~3 (~2) times larger value compared to that measured in MB events. The lower values of $dN_γ/dη_{lab}$ in case of central events selected with the ZNA estimator is, most probably, due to the saturation of forward neutron emission [25]. The systematic uncertainty represented by open boxes mostly comes from the uncertainty on the upstream material in front of the PMD. To understand further the evo-
<figure><img src="image_8.png" /><figcaption>**Fig. 8** Pseudorapidity distributions of inclusive photons measured within $2.3 < \eta_{lab} < 3.9$ in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for several centrality classes and for two centrality estimators: CL1 (top) and ZNA (bottom)</figcaption></figure>

<caption>**Fig. 9** Ratios of $dN_{\gamma}/d\eta_{lab}$ distributions measured in different centrality classes to that in the most peripheral (80–100%) centrality class. The corresponding HIJING predictions are shown by dotted lines</caption>

lution of $dN_{\gamma}/d\eta_{lab}$ with centrality the distributions in each centrality interval are divided by the distribution in the most peripheral (80–100%) event class and presented in Fig. 9. With increasing pseudorapidity the ratios are observed to increase linearly with a slope whose magnitude increases for central events being highest for the most central (0–5%) collisions. HIJING shows a similar trend and explains the data within measured uncertainties except for the 0–5% and 40–60% event classes.

The results are compared in Fig. 10 to similar measurements of charged particles at midrapidity by ALICE [25] for selected centrality intervals (0–5%, 20–40% and 80–100%) determined with the CL1 (top panel) and ZNA (bottom panel) estimators. A clear asymmetric shape of $dN_{ch}/d\eta_{lab}$ is observed for most central collisions and the shape becomes symmetric (like in pp) in the most peripheral

event class. Figure 10 also reports the comparison of the measured $dN_{\gamma}/d\eta_{lab}$ with HIJING predictions for the three considered centrality classes. The $dN_{ch}/d\eta_{lab}$ distributions at midrapidity are better described by the model in the case of the ZNA estimator. In the most peripheral (80–100%) collisions, HIJING overpredicts (underpredicts) the $dN_{ch}/d\eta_{lab}$ for the CL1 (ZNA) estimator. The model describes the values of $dN_{\gamma}/d\eta_{lab}$ at forward rapidity to within 15% except for the 0–5% CL1 class and the 80–100% ZNA class where the model prediction is compatible with the data points within 5%.

The obtained values of average photon multiplicity $\langle N_{\gamma} \rangle$ at forward rapidity ($2.3 < \eta_{lab} < 3.9$) as a function of $\langle N_{part} \rangle$ for the two centrality estimators in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV are presented in Fig. 11(a) and listed in Table 7. The systematic uncertainties are represented by the color bands. The values of $\langle N_{part} \rangle$ for the different centrality classes for both the estimators are taken from [25]. The data point for pp collisions at $\sqrt{s}$ = 5.02 TeV is also included for reference. The average photon multiplicity $\langle N_{\gamma} \rangle$ divided by the average number of participants is presented as a function of $\langle N_{part} \rangle$ in Fig. 11(b). The $\langle N_{\gamma} \rangle / \langle N_{part} \rangle$ ratio has a steeper increase with $\langle N_{part} \rangle$ for the CL1 centrality estimator than for
<figure><img src="image_9.png" /><figcaption>**Fig. 10** Pseudorapidity distribution of inclusive photons measured within $2.3 < \eta_{\text{lab}} < 3.9$ in p–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV compared with the charged-particle measurements [25] and to HIJING predictions for three centrality classes and for two centrality estimators: CL1 (top), ZNA (bottom)</figcaption></figure>

the ZNA one. This could be attributed to the strong multiplicity bias as described in Ref. [25]. It is noted that the curve does not point towards the pp result. When the ZNA is used for the centrality estimation along with the hybrid method to determine the average number of participant nucleons ($N_{\text{part}}^{\text{high-}pT}$ or $N_{\text{part}}^{\text{Pb-side}}$), the inclusive photon multiplicity is found to scale with $\langle N_{\text{part}} \rangle$ within uncertainties and to point towards the pp data point at low $\langle N_{\text{part}} \rangle$. Moreover, the range in $\langle N_{\text{part}} \rangle$ for the ZNA centrality selection is more limited (that could be because of saturation of forward neutron emission) than what is obtained by particle-multiplicity-based centrality estimators. This effect is also emphasized in Fig. 11(c) where the same quantity $\langle N_{\gamma} \rangle/\langle N_{\text{part}} \rangle$ is presented as a function of $\langle N_{\gamma} \rangle$. Similar results were observed for charged-particle multiplicity at midrapidity reported by ALICE in Ref. [25].

## 10 Summary

In summary, we have presented results for inclusive photon multiplicity and pseudorapidity distributions at forward rapidities ($2.3 < \eta_{\text{lab}} < 3.9$) in pp collisions at $\sqrt{s} = 5.02$ TeV and p–Pb and Pb–p collisions at

$\sqrt{s_{NN}} = 5.02$ TeV. The centrality dependence of inclusive photon production is studied in p–Pb collisions for the CL1 and the ZNA centrality estimators. The evolution of the average photon multiplicity $\langle N_{\gamma} \rangle$ with $\sqrt{s}$ in pp collisions and with $\langle N_{\text{part}} \rangle$ in p–Pb collisions is presented. A comparison of inclusive photon production at forward rapidity with the production of charged particles at midrapidity is also discussed. The obtained results are compared to predictions from various event generators (PYTHIA 6 with the Perugia 2011 tune, PYTHIA 8 with the Monash 2013 tune, PHOJET, EPOS LHC for pp collisions and HIJING, DPMJET, AMPT for p–Pb and Pb–p collisions). The multiplicity distributions are well reproduced with double NBDs and results of the parameterization are provided.

The PYTHIA 8 generator with the Monash 2013 tune and the EPOS LHC generator do not reproduce the multiplicity distribution in pp collisions at high multiplicities whereas PYTHIA 6 with the Perugia 2011 tune and PHOJET better describe the data within uncertainties. For p–Pb collisions, none of the considered MC models could describe the multiplicity results in the full range. However, HIJING and AMPT are able to reproduce the tail of the distribution for both p–Pb and Pb–p collisions. Interestingly, compared to multiplicity distributions, the pseudorapidity distributions of inclusive photons are better described by all the considered MC models. The model predictions are found to describe the data within 20%. The pseudorapidity density of inclusive photons at forward rapidity follows the $\eta_{\text{lab}}$ dependence of charged-particle multiplicity at midrapidity in pp and for various centrality classes in p–Pb collisions. This can be attributed to the fact that the major contribution to inclusive photons comes from the decay of neutral pions, and therefore the number of produced photons should be similar to that of charged particles. The models give similar pseudorapidity densities for photons and charged particles only at forward rapidity where the difference between $\eta$ and $y$ is negligible.

An asymmetric $dN_{\text{ch}}/d\eta_{\text{lab}}$ distribution for the asymmetric p–Pb collision system is evident for most central collisions while the distributions become more symmetric for the most peripheral class. HIJING describes the centrality evolution of both photon and charged-particle production within about 10–12%. The average photon multiplicity $\langle N_{\gamma} \rangle$ divided by the average number of participants scales linearly with both $N_{\text{part}}^{\text{high-}pT}$ and $N_{\text{part}}^{\text{Pb-side}}$ estimated by the energy deposited in the ZDC and the trend with $\langle N_{\text{part}} \rangle$ is consistent with the data point from pp collisions at $\sqrt{s} = 5.02$ TeV.

These results of inclusive photon production in pp and p–Pb collisions provide important input for the tuning of theoretical models and MC event generators and help to establish the baseline measurements to interpret the Pb–Pb data.
**Fig. 11** (a) The average photon multiplicity $\langle N_Y \rangle$ within $2.3 < \eta_{lab} < 3.9$ as a function of the average number of participants for various centrality estimators in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Values of $\langle N_Y \rangle / \langle N_{part} \rangle$ are shown as a function of $\langle N_{part} \rangle$ (b) and $\langle N_Y \rangle$ (c). The lines between the data points are to guide the eye. The data point for pp collisions at $\sqrt{s} = 5.02$ TeV is also included for reference

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Table 7** The average photon multiplicity $\langle N_Y \rangle$ within $2.3 < \eta_{lab} < 3.9$ for various centrality classes defined using CL1 and ZNA centrality estimators in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The corresponding</figcaption></figure>

values of $\langle N_{part} \rangle$ are taken from Ref. [25]. The quoted errors are systematic uncertainties. Statistical uncertainties are negligible

<table><thead><tr><th rowspan="2">Centrality</th><th colspan="2">CL1</th><th colspan="3">ZNA</th></tr><tr><th>⟨N<sup>NBD-Glauber</sup><sub>part</sub>⟩</th><th>⟨N<sub>Y</sub>⟩</th><th>⟨N<sup>high-p-T</sup><sub>part</sub>⟩</th><th>⟨N<sup>Pb-side</sup><sub>part</sub>⟩</th><th>⟨N<sub>Y</sub>⟩</th></tr></thead><tbody><tr><td>0–5%</td><td>16.60 ± 1.66</td><td>84.44 ± 8.65</td><td>13.50 ± 1.08</td><td>14.30 ± 1.14</td><td>54.98 ± 5.59</td></tr><tr><td>5–10%</td><td>14.60 ± 1.46</td><td>65.82 ± 6.96</td><td>13.10 ± 1.05</td><td>13.30 ± 1.06</td><td>51.65 ± 5.09</td></tr><tr><td>10–20%</td><td>13.00 ± 1.30</td><td>54.23 ± 5.67</td><td>12.30 ± 0.98</td><td>12.40 ± 0.99</td><td>48.18 ± 4.82</td></tr><tr><td>20–40%</td><td>10.49 ± 0.94</td><td>39.95 ± 4.22</td><td>10.73 ± 0.86</td><td>10.60 ± 0.85</td><td>41.22 ± 4.22</td></tr><tr><td>40–60%</td><td>7.18 ± 0.52</td><td>26.28 ± 2.79</td><td>7.81 ± 0.62</td><td>7.74 ± 0.62</td><td>29.82 ± 3.20</td></tr><tr><td>60–80%</td><td>4.40 ± 0.88</td><td>15.21 ± 1.60</td><td>5.05 ± 0.40</td><td>5.00 ± 0.40</td><td>18.36 ± 1.93</td></tr><tr><td>80–100%</td><td>2.76 ± 0.63</td><td>6.34 ± 0.68</td><td>3.03 ± 0.24</td><td>3.06 ± 0.24</td><td>9.93 ± 1.05</td></tr></tbody></table>

**Acknowledgements** The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences, Austrian Science Fund (FWF): [M 2467-N36] and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de

Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) and Universidade Federal do Rio Grande do Sul (UFRGS), Brazil; Bulgarian Ministry of Education and Science, within the National Roadmap for Research Infrastructures 2020-2027 (object CERN), Bulgaria; Ministry of Education of China (MOEC), Ministry of Science & Technology of China (MSTC) and National Natural Science Foundation of China (NSFC), China; Ministry of Science and Education and Croatian Science Foundation, Croatia; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research | Natural Sciences, the VILLUM FONDEN and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l'Energie Atomique (CEA) and
Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung und Forschung (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Research and Religions, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE), Department of Science and Technology, Government of India (DST), University Grants Commission, Government of India (UGC) and Council of Scientific and Industrial Research (CSIR), India; National Research and Innovation Agency-BRIN, Indonesia; Istihtuto Nazionale di Fisica Nucleare (INFN), Italy; Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and Japan Society for the Promotion of Science (JSPS) KAKENHI, Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Education and Science, National Science Centre and WUT ID-UB, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics, Ministry of Research and Innovation and Institute of Atomic Physics and University Politehnica of Bucharest, Romania; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; Suranaree University of Technology (SUT), National Science and Technology Development Agency (NSTDA), Thailand Science Research and Innovation (TSRI) and National Science, Research and Innovation Fund (NSRF), Thailand; Turkish Energy, Nuclear and Mineral Research Agency (TENMAK),

Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America. In addition, individual groups or members have received support from: European Research Council, Strong 2020-Horizon 2020, Marie Skłodowska Curie (Grant Nos. 950692, 824093, 896850), European Union; Academy of Finland (Center of Excellence in Quark Matter) (Grant Nos. 346327, 346328), Finland; Programa de Apoyos para la Superación del Personal Académico, UNAM, Mexico.

**Data Availability Statement** This manuscript has no associated data or the data will not be deposited. [Authors' comment: Manuscript has associated data in the HEPData repository at https://www.hepdata.net/record/ins2637678.]

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Funded by SCOAP³. SCOAP³ supports the goals of the International Year of Basic Sciences for Sustainable Development.

## Appendix A

See Figs. 12 and 13.
<figure><img src="image_13.png" /><figcaption>**Fig. 13** The mean (a) and width (b) of the distribution of reconstructed photon multiplicity ($N_{\gamma-\text{like}}$) as a function of true photon multiplicity ($N_{\gamma-\text{true}}$) are presented for p–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The dotted lines correspond to $N_{\gamma-\text{true}} = 140$ up to which the results are reported in the paper</figcaption></figure>

## References

1. P. Braun-Munzinger, V. Koch, T. Schäfer, J. Stachel, Properties of hot and dense matter from relativistic heavy-ion collisions. Phys. Rep. **621**, 76–126 (2016). https://doi.org/10.1016/j.physrep.2015.12.003. arXiv:1510.00442 [nucl-th]

2. J. Schukraft, Heavy Ion physics with the ALICE experiment at the CERN LHC. Philos. Trans. R. Soc. Lond. **A 370**, 917–932 (2012). https://doi.org/10.1098/rsta.2011.0469. arXiv:1109.4291 [hep-ex]

3. **ALICE** Collaboration, The ALICE experiment—a journey through QCD. arXiv:2211.04384 [nucl-ex]

4. R. Xu, W.-T. Deng, X.-N. Wang, Nuclear modification of high-$p_T$ hadron spectra in p-A collisions at LHC. Phys. Rev. C **86**, 051901 (2012). https://doi.org/10.1103/PhysRevC.86.051901. arXiv:1204.1998 [nucl-th]

5. M. Arneodo, Nuclear effects in structure functions. Phys. Rep. **240**, 301–393 (1994). https://doi.org/10.1016/0370-1573(94)90048-5

6. CMS Collaboration, V. Khachatryan et al., Observation of long-range near-side angular correlations in proton-proton collisions at the LHC. JHEP **09**, 091 (2010). https://doi.org/10.1007/JHEP09(2010)091. arXiv:1009.4122 [hep-ex]

7. CMS Collaboration, V. Khachatryan et al., Measurement of long-range near-side two-particle angular correlations in pp collisions at $\sqrt{s}=13$ TeV. Phys. Rev. Lett. **116**, 172302 (2016). https://doi.org/10.1103/PhysRevLett.116.172302. arXiv:1510.03068 [nucl-ex]

8. **ATLAS** Collaboration, G. Aad et al., Observation of long-range elliptic azimuthal anisotropies in $\sqrt{s}=13$ and $2.76$ TeV pp collisions with the ATLAS detector. Phys. Rev. Lett. **116**, 172301 (2016). https://doi.org/10.1103/PhysRevLett.116.172301. arXiv:1509.04776 [hep-ex]

9. **ALICE** Collaboration, B. Abelev et al., Long-range angular correlations on the near and away side in p–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Phys. Lett. B **719**, 29–41 (2013). https://doi.org/10.1016/j.physletb.2013.01.012. arXiv:1212.2001 [nucl-ex]

10. **ATLAS** Collaboration, G. Aad et al., Observation of associated near-side and away-side long-range correlations in $\sqrt{s_{NN}} = 5.02$ TeV p–Pb Collisions with the ATLAS detector. Phys. Rev. Lett. **110**, 182302 (2013). https://doi.org/10.1103/PhysRevLett.110.182302. arXiv:1212.5198 [hep-ex]

11. CMS Collaboration, S. Chatrchyan et al., Observation of long-range near-side angular correlations in p-Pb collisions at the LHC. Phys. Lett. B **718**, 795–814 (2013). https://doi.org/10.1016/j.physletb.2012.11.025. arXiv:1210.5482 [nucl-ex]

12. **ATLAS** Collaboration, G. Aad et al., Measurement of long-range pseudorapidity correlations and azimuthal harmonics in $\sqrt{s_{NN}} = 5.02$ TeV p–Pb collisions with the ATLAS detector. Phys. Rev. C **90**, 044906 (2014). https://doi.org/10.1103/PhysRevC.90.044906. arXiv:1409.1792 [hep-ex]

13. **ALICE** Collaboration, J. Adam et al., Enhanced production of multi-strange hadrons in high-multiplicity proton–proton collisions. Nat. Phys. **13**, 535–539 (2017). https://doi.org/10.1038/nphys4111. arXiv:1606.07424 [nucl-ex]

14. D. Kharzeev, E. Levin, M. Nardi, Color glass condensate at the LHC: Hadron multiplicities in pp, p-A and A-A collisions. Nucl. Phys. A **747**, 609–629 (2005). https://doi.org/10.1016/j.nuclphysa.2004.10.018. arXiv:hep-ph/0408050

15. W.-T. Deng, X.-N. Wang, R. Xu, Hadron production in pp, p-Pb, and Pb-Pb collisions with the HIJING 2.0 model at energies available at the CERN Large Hadron Collider. Phys. Rev. C **83**, 014915 (2011). https://doi.org/10.1103/PhysRevC.83.014915. arXiv:1008.1841 [hep-ph]

16. **ALICE** Collaboration, B.B. Abelev et al., Inclusive photon production at forward rapidities in proton–proton collisions at $\sqrt{s} = 0.9, 2.76$ and $7$ TeV. Eur. Phys. J. C **75**, 146 (2015). https://doi.org/10.1140/epjc/s10052-015-3356-2. arXiv:1411.4981 [nucl-ex]

17. **STAR** Collaboration, J. Adams et al., Multiplicity and pseudorapidity distributions of photons in Au + Au collisions at $\sqrt{s_{NN}} = 62.4$ GeV. Phys. Rev. Lett. **95**, 062301 (2005). https://doi.org/10.1103/PhysRevLett.95.062301. arXiv:nucl-ex/0502008
18. **STAR** Collaboration, J. Adams et al., Multiplicity and pseudorapidity distributions of charged particles and photons at forward pseudorapidity in Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV. Phys. Rev. C **73**, 034906 (2006). https://doi.org/10.1103/PhysRevC.73.034906. [arXiv:nucl-ex/0511026]

19. **ALICE** Collaboration, K. Aamodt et al., First proton–proton collisions at the LHC as observed with the ALICE detector: measurement of the charged-particle pseudorapidity density at $\sqrt{s}$ = 900 GeV. Eur. Phys. J. C **65**, 111–125 (2010). https://doi.org/10.1140/epjc/s10052-009-1227-4. [arXiv:0911.5430 [hep-ex]]

20. **ALICE** Collaboration, K. Aamodt et al., Charged-particle multiplicity measurement in proton–proton collisions at $\sqrt{s}$ = 0.9 and 2.36 TeV with ALICE at LHC. Eur. Phys. J. C **68**, 89–108 (2010). https://doi.org/10.1140/epjc/s10052-010-1339-x. [arXiv:1004.3034 [hep-ex]]

21. **ALICE** Collaboration, K. Aamodt et al., Charged-particle multiplicity measurement in proton–proton collisions at $\sqrt{s}$ = 7 TeV with ALICE at LHC. Eur. Phys. J. C **68**, 345–354 (2010). https://doi.org/10.1140/epjc/s10052-010-1350-2. [arXiv:1004.3514 [hep-ex]]

22. **ALICE** Collaboration, J. Adam et al., Charged-particle multiplicities in proton–proton collisions at $\sqrt{s}$ = 0.9 to 8 TeV. Eur. Phys. J. C **77**, 33 (2017). https://doi.org/10.1140/epjc/s10052-016-4571-1. [arXiv:1509.07541 [nucl-ex]]

23. **ALICE** Collaboration, S. Acharya et al., Charged-particle multiplicity distributions over a wide pseudorapidity range in proton–proton collisions at $\sqrt{s}$ = 0.9, 7, and 8 TeV. Eur. Phys. J. C **77**, 852 (2017). https://doi.org/10.1140/epjc/s10052-017-5412-6. [arXiv:1708.01435 [hep-ex]]

24. **ALICE** Collaboration, B. Abelev et al., Pseudorapidity density of charged particles in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Phys. Rev. Lett. **110**, 032301 (2013). https://doi.org/10.1130/PhysRevLett.110.032301. [arXiv:1210.3615 [nucl-ex]]

25. **ALICE** Collaboration, J. Adam et al., Centrality dependence of particle production in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Phys. Rev. C **91**, 064905 (2015). https://doi.org/10.1130/PhysRevC.91.064905. [arXiv:1412.6828 [nucl-ex]]

26. **ALICE** Collaboration, S. Acharya et al., Charged-particle pseudorapidity density at midrapidity in p–Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV. Eur. Phys. J. C **79**, 307 (2019). https://doi.org/10.1140/epjc/s10052-019-6801-9. [arXiv:1812.01312 [nucl-ex]]

27. **ALICE** Collaboration, S. Acharya et al., Pseudorapidity distributions of charged particles as a function of mid- and forward rapidity multiplicities in pp collisions at $\sqrt{s}$ = 5.02, 7 and 13 TeV. Eur. Phys. J. C **81**, 630 (2021). https://doi.org/10.1140/epjc/s10052-021-09349-5. [arXiv:2009.09434 [nucl-ex]]

28. **ALICE** Collaboration, K. Aamodt et al., Charged-particle multiplicity density at midrapidity in central Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. Phys. Rev. Lett. **105**, 252301 (2010). https://doi.org/10.1130/PhysRevLett.105.252301. [arXiv:1011.3916 [nucl-ex]]

29. **ALICE** Collaboration, K. Aamodt et al., Centrality dependence of the charged-particle multiplicity density at midrapidity in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. Phys. Rev. Lett. **106**, 032301 (2011). https://doi.org/10.1130/PhysRevLett.106.032301. [arXiv:1012.1657 [nucl-ex]]

30. **ALICE** Collaboration, J. Adam et al., Centrality evolution of the charged-particle pseudorapidity density over a broad pseudorapidity range in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. Phys. Lett. B **754**, 373–385 (2016). https://doi.org/10.1016/j.physletb.2015.12.082. [arXiv:1509.07299 [nucl-ex]]

31. **ALICE** Collaboration, E. Abbas et al., Centrality dependence of the pseudorapidity density distribution for charged particles in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. Phys. Lett. B **726**, 610–622 (2013). https://doi.org/10.1016/j.physletb.2013.09.022. [arXiv:1304.0347 [nucl-ex]]

32. **ALICE** Collaboration, J. Adam et al., Centrality dependence of the charged-particle multiplicity density at midrapidity in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Phys. Rev. Lett. **116**, 222302 (2016). https://doi.org/10.1103/PhysRevLett.116.222302. [arXiv:1512.06104 [nucl-ex]]

33. **ALICE** Collaboration, J. Adam et al., Centrality dependence of the pseudorapidity density distribution for charged particles in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Phys. Lett. B **772**, 567–577 (2017). https://doi.org/10.1016/j.physletb.2017.07.017. [arXiv:1612.08966 [nucl-ex]]

34. **ATLAS** Collaboration, G. Aad et al., Charged-particle multiplicities in pp interactions at $\sqrt{s}$ = 900 GeV measured with the ATLAS detector at the LHC. Phys. Lett. B **688**, 21–42 (2010). https://doi.org/10.1016/j.physletb.2010.03.064. [arXiv:1003.3124 [hep-ex]]

35. **ATLAS** Collaboration, G. Aad et al., Charged-particle multiplicities in pp interactions measured with the ATLAS detector at the LHC. New J. Phys. **13**, 053033 (2011). https://doi.org/10.1088/1367-2630/13/5/053033. [arXiv:1012.5104 [hep-ex]]

36. **ATLAS** Collaboration, G. Aad et al., Charged-particle distributions in pp interactions at $\sqrt{s}$ = 8 TeV measured with the ATLAS detector. Eur. Phys. J. C **76**, 403 (2016). https://doi.org/10.1140/epjc/s10052-016-4203-9. [arXiv:1603.02439 [hep-ex]]

37. **ATLAS** Collaboration, G. Aad et al., Measurement of the centrality dependence of the charged-particle pseudorapidity distribution in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector. Eur. Phys. J. C **76**, 199 (2016). https://doi.org/10.1140/epjc/s10052-016-4002-3. [arXiv:1508.00848 [hep-ex]]

38. **ATLAS** Collaboration, G. Aad et al., Measurement of the centrality dependence of the charged-particle pseudorapidity distribution in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV with the ATLAS detector. Phys. Lett. B **710**, 363–382 (2012). https://doi.org/10.1016/j.physletb.2012.02.045. [arXiv:1108.6027 [hep-ex]]

39. **CMS** Collaboration, V. Khachatryan et al., Charged-particle multiplicities in pp interactions at $\sqrt{s}$ = 0.9, 2.36, and 7 TeV. JHEP **01**, 079 (2011). https://doi.org/10.1007/JHEP01(2011)079. [arXiv:1011.5531 [hep-ex]]

40. **CMS** Collaboration, V. Khachatryan et al., Pseudorapidity distribution of charged hadrons in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Phys. Lett. B **751**, 143–163 (2015). https://doi.org/10.1016/j.physletb.2015.10.004. [arXiv:1507.05915 [hep-ex]]

41. **CMS** Collaboration, A.M. Sirunyan et al., Pseudorapidity distributions of charged hadrons in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 and 8.16 TeV. JHEP **01**, 045 (2018). https://doi.org/10.1007/JHEP01(2018)045. [arXiv:1710.09355 [hep-ex]]

42. **CMS** Collaboration, S. Chatrchyan et al., Dependence on pseudorapidity and centrality of charged hadron production in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. JHEP **08**, 141 (2011). https://doi.org/10.1007/JHEP08(2011)141. [arXiv:1107.4800 [nucl-ex]]

43. **LHCb** Collaboration, R. Aaij et al., Measurement of charged-particle multiplicities and densities in pp collisions at $\sqrt{s}$ = 7 TeV in the forward region. Eur. Phys. J. C **74**, 2888 (2014). https://doi.org/10.1140/epjc/s10052-014-2888-1. [arXiv:1402.4430 [hep-ex]]

44. **ALICE** Collaboration, K. Aamodt et al., The ALICE experiment at the CERN LHC. JINST **3**, S08002 (2008). https://doi.org/10.1088/1748-0221/3/08/S08002

45. **ALICE** Collaboration, B.B. Abelev et al., Performance of the ALICE experiment at the CERN LHC. Int. J. Mod. Phys. A **29**, 1430044 (2014). https://doi.org/10.1142/S0217751X14300440. [arXiv:1402.4476 [nucl-ex]]

46. **ALICE** Collaboration, G. Dellacasa et al., ALICE technical design report: Photon multiplicity detector (PMD). CERN-LHCC-99-032. http://cds.cern.ch/record/451099

47. **ALICE** Collaboration, P. Cortese et al., ALICE: addendum to the technical design report of the photon multiplicity detector (PMD). CERN-LHCC-2003-038. https://cds.cern.ch/record/642177
48. **ALICE** Collaboration, K. Aamodt et al., Alignment of the ALICE Inner Tracking System with cosmic-ray tracks. JINST **5**, P03003 (2010). https://doi.org/10.1088/1748-0221/5/03/P03003. arXiv:1001.0502 [physics.ins-det]

49. ALICE Collaboration, P. Cortese et al., ALICE forward detectors: FMD, TO and VO: Technical Design Report. CERN-LHCC-2004-025. https://cds.cern.ch/record/781854

50. ALICE Collaboration, E. Abbas et al., Performance of the ALICE VZ ZERO system. JINST **8**, P10016 (2013). https://doi.org/10.1088/1748-0221/8/10/P10016. arXiv:1306.3130 [nucl-ex]

51. ALICE Collaboration, M. Gallio, W. Klempt, L. Leistam, J. De Groot, J. Schukraft, ALICE Zero Degree Calorimeter (ZDC): Technical Design Report. CERN-LHCC-99-005. https://cds.cern.ch/record/381433

52. ALICE Collaboration, L. Betev et al., Definition of the ALICE coordinate system and basic rules for sub-detector components numbering. https://edms.cern.ch/ui/#!master/navigator/document?D:1020137949:1020137949:subDocs

53. M.L. Miller, K. Reygers, S.J. Sanders, P. Steinberg, Glauber modeling in high energy nuclear collisions. Ann. Rev. Nucl. Part. Sci. **57**, 205-243 (2007). https://doi.org/10.1146/annurev.nucl.57.090506.123020. arXiv:nucl-ex/0701025

54. D. d'Enterria, C. Loizides, Progress in the Glauber model at collider energies. Ann. Rev. Nucl. Part. Sci. **71**, 315-344 (2021). https://doi.org/10.1146/annurev-nucl-102419-060007. arXiv:2011.14909 [hep-ph]

55. F. Sikler, Centrality control of hadron nucleus interactions by detection of slow nucleons. arXiv:hep-ph/0304065

56. B. Alessandro et al., ALICE: physics performance report, volume II. J. Phys. G Nucl. Part. Phys. **32**, 1295-2040 (2006). https://doi.org/10.1088/0954-3899/32/10/001

57. X.-N. Wang, M. Gyulassy, HIJING: a Monte Carlo model for multiple jet production in pp, p-A and A-A collisions. Phys. Rev. D **44**, 3501-3516 (1991). https://doi.org/10.1103/PhysRevD.44.3501

58. S. Roesler, R. Engel, J. Ranft, The Monte Carlo event generator DPMJET-III, in *International Conference on Advanced Monte Carlo for Radiation Physics, Particle Transport Simulation and Applications (MC 2000)*, vol. 12, pp. 1033-1038 (2000). https://doi.org/10.1007/978-3-642-18211-2_166. arXiv:hep-ph/0012252 [hep-ph]

59. P. Skands, S. Carrazza, J. Rojo, Tuning PYTHIA 8.1: the Monash 2013 Tune. Eur. Phys. J. C **74**, 3024 (2014). https://doi.org/10.1140/epjc/s10052-014-3024-y. arXiv:1404.5630 [hep-ph]

60. T. Pierog et al., EPOS LHC: test of collective hadronization with data measured at the CERN Large Hadron Collider. Phys. Rev. C **92**, 034906 (2015). https://doi.org/10.1103/PhysRevC.92.034906. arXiv:1306.0121 [hep-ph]

61. R. Brun et al., GEANT: Detector Description and Simulation Tool. CERN Program Library. CERN, Geneva (1993). https://doi.org/10.17181/CERN.MUHF.DMJ1http://cds.cern.ch/record/1082634

62. R. Brun et al., Computing in ALICE. Nucl. Instrum. Methods A **502**, 339-346 (2003). https://doi.org/10.1016/S0168-9002(03)00440-6

63. C. Bierlich et al., A comprehensive guide to the physics and usage of PYTHIA 8.3. arXiv:2203.11601 [hep-ph]

64. T. Sjöstrand, M. van Zijl, A multiple-interaction model for the event structure in hadron collisions. Phys. Rev. D **36**, 2019-2041 (1987). https://doi.org/10.1103/PhysRevD.36.2019

65. T. Sjöstrand et al., An introduction to PYTHIA 8.2. Comput. Phys. Commun. **191**, 159-177 (2015). https://doi.org/10.1016/j.cpc.2015.01.024. arXiv:1410.3012 [hep-ph]

66. J.R. Christiansen, P.Z. Skands, String formation beyond leading colour. JHEP **08**, 003 (2015). https://doi.org/10.1007/JHEP08(2015)003. arXiv:1505.01681 [hep-ph]

67. H. Pi, An event generator for interactions between hadrons and nuclei: FRITIOF version 7.0. Comput. Phys. Commun. **71**, 173-192 (1992). https://doi.org/10.1016/0010-4655(92)90082-A

68. A. Capella, U. Sukhatme, C. I. Tan, J. Tran Thanh Van, Dual parton model. Phys. Rep. **236**, 225-329 (1994) https://doi.org/10.1016/0370-1573(94)90064-7

69. P.D.B. Collins, An introduction to Regge theory and high-energy physics. Cambridge Monographs on Mathematical Physics. Cambridge Univ. Press, Cambridge, p5 (2009). https://doi.org/10.1017/CBO9780511897603

70. CDF Collaboration, F. Abe et al., Measurement of $p\bar{p}$ single diffractive dissociation at $\sqrt{s} = 546$ GeV and 1800 GeV. Phys. Rev. D **50**, 5535-5549 (1994). https://doi.org/10.1103/PhysRevD.50.5535

71. CDF Collaboration, F. Abe et al., Observation of rapidity gaps in $p\bar{p}$ collisions at 1.8 TeV. Phys. Rev. Lett. **74**, 855-859 (1995). https://doi.org/10.1103/PhysRevLett.74.855

72. ALICE Collaboration, B. Abelev et al., Measurement of inelastic, single- and double-diffraction cross sections in proton-proton collisions at the LHC with ALICE. Eur. Phys. J. C **73**, 2456 (2013). https://doi.org/10.1140/epjc/s10052-013-2456-0. arXiv:1208.4968 [hep-ex]

73. G. D'Agostini, A multidimensional unfolding method based on Bayes' theorem. Nucl. Instrum. Methods Phys. Res. Sect. A **362**, 487-498 (1995). https://doi.org/10.1016/0168-9002(95)00274-X

74. T. Adye, Unfolding algorithms and tests using RooUnfold. arXiv:1105.1160 [physics.data-an]

75. A. Hocker, V. Kartvelishvili, SVD approach to data unfolding. Nucl. Instrum. Methods A **372**, 469-481 (1996). https://doi.org/10.1016/0168-9002(95)01478-0. arXiv:hep-ph/9509307

76. P.Z. Skands, Tuning Monte Carlo generators: the Perugia tunes. Phys. Rev. D **82**, 074018 (2010). https://doi.org/10.1103/PhysRevD.82.074018. arXiv:1005.3457 [hep-ph]

77. F.W. Bopp, R. Engel, J. Ranft, Rapidity gaps and the PHOJET Monte Carlo, in *LAFEX International School on High-Energy Physics (LISHEP 98) Session A: Particle Physics for High School Teachers—Session B: Advanced School in HEP—Session C: Workshop on Diffractive Physics*, vol. 3, pp. 729-741 (1998). arXiv:hep-ph/9803437

78. Z.-W. Lin, C.M. Ko, B.-A. Li, B. Zhang, S. Pal, Multiphase transport model for relativistic heavy-ion collisions. Phys. Rev. C **72**, 064901 (2005). https://doi.org/10.1103/PhysRevC.72.064901. arXiv:nucl-th/0411110

79. ALICE Collaboration, Pseudorapidity densities of charged particles with transverse momentum thresholds in pp collisions at $\sqrt{s} = 5.02$ and 13 TeV. arXiv:2211.15364 [nucl-ex]
ALICE Collaboration*

S. Acharya¹²⁵, D. Adamová⁶⁶, A. Adler⁶⁹, G. Aglieri Rinella³², M. Agnello²⁹, N. Agrawal⁵⁰, Z. Ahammed¹³², S. Ahmad¹⁵, S. U. Ahn⁷⁰, I. Ahuja³⁷, A. Akindinov¹⁴⁰, M. Al-Turany⁹⁷, D. Aleksandrov¹⁴⁰, B. Alessandro⁵⁵, H. M. Alfanda⁶, R. Alfaro Molina⁶⁶, B. Ali¹⁵, A. Alici²⁵, N. Alizadehvandchali¹¹⁴, A. Alkin³², J. Alme²⁰, G. Alocco⁵¹, T. Alt⁶³, I. Altsybeev¹⁴⁰, M. N. Anaam⁶, C. Andrei⁴⁵, A. Andronic¹³⁵, V. Anguelov⁹⁴, F. Antinori⁵³, P. Antonioli⁵⁰, N. Apadula⁷⁴, L. Aphecetche¹⁰³, H. Appelshäuser⁶³, C. Arata⁷³, S. Arcelli²⁵, M. Aresti⁵¹, R. Arnaldi⁵⁵, J. G. M. C. A. Arneiro¹¹⁰, I. C. Arsene¹⁹, M. Arslandok¹³⁷, A. Augustinus³², R. Averbeck⁹⁷, M. D. Azmi¹⁵, A. Badalà⁵², J. Bae¹⁰⁴, Y. W. Baek⁴⁰, X. Bai¹¹⁸, R. Bailhache⁶³, Y. Bailung⁴⁷, A. Balbino²⁹, A. Baldisseri¹²⁸, B. Balis², D. Banerjee⁴, Z. Banoo⁹¹, R. Barbera²⁶, F. Barile³¹, L. Barioglio⁹⁵, M. Barlou⁷⁸, G. G. Barnaföldi¹³⁶, L. S. Barnby⁸⁵, V. Barret¹²⁵, L. Barreto¹¹⁰, C. Bartels¹¹⁷, K. Barth³², E. Bartsch⁶³, N. Bastid¹²⁵, S. Basu⁷⁵, G. Batigne¹⁰³, D. Battistini⁹⁵, B. Batyunya¹⁴¹, D. Bauri⁴⁶, J. L. Bazo Alba¹⁰¹, I. G. Bearden⁸³, C. Beattie¹³⁷, P. Becht⁹⁷, D. Behera⁴⁷, I. Belikov¹²⁷, A. D. C. Bell Hechavarria¹³⁵, F. Bellini²⁵, R. Bellwied¹¹⁴, S. Belokurova¹⁴⁰, V. Belyaev¹⁴⁰, G. Bencedi¹³⁶, S. Beole²⁴, A. Bercuci⁴⁵, Y. Berdnikov¹⁴⁰, A. Berdnikova⁹⁴, L. Bergmann⁹⁴, M. G. Besoiu⁶², L. Betev³², P. P. Bhaduri¹³², A. Bhasin⁹¹, M. A. Bhat⁴, B. Bhattacharjee⁴¹, L. Bianchi²⁴, N. Bianchi⁴⁸, J. Bielčík³⁵, J. Bielčíková⁶⁸, J. Biernat¹⁰⁷, A. P. Bigot¹²⁷, A. Bilandzic⁹⁵, G. Biro¹³⁶, S. Biswas⁴, N. Bize¹⁰³, J. T. Blair¹⁰⁸, D. Blau¹⁴⁰, M. B. Blidaru⁹⁷, N. Bluhme³⁸, C. Blume⁶³, G. Boca²¹,⁵⁴, F. Bock⁸⁷, T. Bodova²⁰, A. Bogdanov¹⁴⁰, S. Boi²², J. Bok⁵⁷, L. Boldizsár¹³⁶, M. Bombara³⁷, P. M. Bond³², G. Bonomi⁵⁴,¹³¹, H. Borel¹²⁸, A. Borissov¹⁴⁰, A. G. Borquez Carcamo⁹⁴, H. Bossi¹³⁷, E. Botta²⁴, Y. E. M. Bouziani⁶³, L. Bratrud⁶³, P. Braun-Munzinger⁹⁷, M. Bregant¹¹⁰, M. Broz³⁵, G. E. Bruno³¹,⁹⁶, M. D. Buckland²³, D. Budnikov¹⁴⁰, H. Buesching⁶³, S. Bufalino²⁹, P. Buhler¹⁰², Z. Buthelezi⁶⁷,¹²¹, A. Bylinkin²⁰, S. A. Bysiak¹⁰⁷, M. Cai⁶, H. Caines¹³⁷, A. Caliva⁹⁷, E. Calvo Villar¹⁰¹, J. M. M. Camacho¹⁰⁹, P. Camerini²³, F. D. M. Canedo¹¹⁰, M. Carabas¹²⁴, A. A. Carballo³², F. Carnesecchi³², R. Caron¹²⁶, L. A. D. Carvalho¹¹⁰, J. Castillo Castellanos¹²⁸, F. Catalano²⁴,³², C. Ceballos Sanchez¹⁴¹, I. Chakaberia⁷⁴, P. Chakraborty⁴⁶, S. Chandra¹³², S. Chapeland³², M. Chartier¹¹⁷, S. Chattopadhyay¹³², S. Chattopadhyay⁹⁹, T. G. Chavez⁴⁴, T. Cheng⁶,⁹⁷, C. Cheshkov¹²⁶, B. Cheynis¹²⁶, V. Chibante Barroso³², D. D. Chinellato¹¹¹, E. S. Chizzali⁹⁵,ª, J. Cho⁵⁷, S. Cho⁵⁷, P. Chochula³², P. Christakoglou⁸⁴, C. H. Christensen⁸³, P. Christiansen⁷⁵, T. Chujó¹²³, M. Ciacco²⁹, C. Cicalo⁵¹, F. Cindolo⁵⁰, M. R. Ciupke⁹⁷, G. Clai⁵⁰,ᵇ, F. Colamaria⁴⁹, J. S. Colburn¹⁰⁰, D. Colella³¹,⁹⁶, M. Colocci²⁵, G. Conesa Balbastre⁷³, Z. Conesa del Valle⁷², G. Contin²³, J. G. Contreras³⁵, M. L. Coquet¹²⁸, T. M. Cormier⁸⁷,f, P. Cortese⁵⁵,¹³⁰, M. R. Cosentino¹¹², F. Costa³², S. Costanza²¹,⁵⁴, C. Cot⁷², J. Crkovská⁹⁴, P. Crochet¹²⁵, R. Cruz-Torres⁷⁴, P. Cui⁶, A. Dainese⁵³, M. C. Danisch⁹⁴, A. Danu⁶², P. Das⁸⁰, P. Das⁴, S. Das⁴, A. R. Dash¹³⁵, S. Dash⁴⁶, R. M. H. David⁴⁴, S. De⁸⁰, A. De Caro²⁸, G. de Cataldo⁴⁹, J. de Cuveland³⁸, A. De Falco²², D. De Gruttola²⁸, N. De Marco⁵⁵, C. De Martin²³, S. De Pasquale²⁸, R. Deb¹³¹, S. Deb⁴⁷, R. J. Debski², K. R. Deja¹³³, R. Del Grande⁹⁵, L. Dello Stritto²⁸, W. Deng⁶, P. Dhankher¹⁸, D. Di Bari³¹, A. Di Mauro³², R. A. Diaz⁷,¹⁴¹, T. Dietel¹¹³, Y. Ding⁶, R. Divià³², D. U. Dixit¹⁸, Ø. Djuvsland²⁰, U. Dmitrieva¹⁴⁰, A. Dobrin⁶², B. Dönigus⁶³, J. M. Dubinski¹³³, A. Dubla⁹⁷, S. Dudi⁹⁰, P. Dupieux¹²⁵, M. Durkac¹⁰⁶, N. Dzalaiova¹², T. M. Eder¹³⁵, R. J. Ehlers⁷⁴, F. Eisenhut⁶³, D. Elia⁴⁹, B. Erazmus¹⁰³, F. Ercolessi²⁵, F. Erhardt⁸⁹, M. R. Ersdal²⁰, B. Espagnon⁷², G. Eulisse³², D. Evans¹⁰⁰, S. Evdokimov¹⁴⁰, L. Fabbietti⁹⁵, M. Faggin²⁷, J. Faivre⁷³, F. Fan⁶, W. Fan⁷⁴, A. Fantoni⁴⁸, M. Fasel⁸⁷, P. Fecchio²⁹, A. Feliciello⁵⁵, G. Feofilov¹⁴⁰, A. Fernández Téllez⁴⁴, L. Ferrandi¹¹⁰, M. B. Ferrer³², A. Ferrero¹²⁸, C. Ferrero⁵⁵, A. Ferretti²⁴, V. J. G. Feuillard⁹⁴, V. Filova³⁵, D. Finogeev¹⁴⁰, F. M. Fionda⁵¹, F. Flor¹¹⁴, A. N. Flores¹⁰⁸, S. Foertsch⁶⁷, I. Fokin⁹⁴, S. Fokin¹⁴⁰, E. Fragiacomo⁵⁶, E. Frajna¹³⁶, U. Fuchs³², N. Funicello²⁸, C. Furget⁷³, A. Furs¹⁴⁰, T. Fusayasu⁹⁸, J. J. Gaardhøje⁸³, M. Gagliardi²⁴, A. M. Gago¹⁰¹, C. D. Galvan¹⁰⁹, D. R. Gangadharan¹¹⁴, P. Ganoti⁷⁸, C. Garabatos⁹⁷, J. R. A. Garcia⁴⁴, E. Garcia-Solis⁹, C. Gargiulo³², K. Garner¹³⁵, P. Gasik⁹⁷, A. Gautam¹¹⁶, M. B. Gay Ducati⁶⁵, M. Germain¹⁰³, A. Ghimouz¹²³, C. Ghosh¹³², M. Giacalone²⁵,⁵⁰, P. Giubellino⁵⁵,⁹⁷, P. Giubilato²⁷, A. M. C. Glaenzer¹²⁸, P. Glässel⁹⁴, E. Glimos¹²⁰, D. J. Q. Goh⁷⁶, V. Gonzalez¹³⁴, M. Gorgon², S. Gotovac³³, V. Grabski⁶⁶, L. K. Graczykowski¹³³, E. Grecka⁸⁶, A. Grelli⁵⁸, C. Grigoras³², V. Grigoriev¹⁴⁰, S. Grigoryan¹,¹⁴¹, F. Grosa³², J. F. Grosse-Oetringhaus³², R. Grosso⁹⁷, D. Grund³⁵, G. G. Guardiano¹¹¹, R. Guernane⁷³, M. Guilbaud¹⁰³, K. Gulbrandsen⁸³, T. Gundem⁶³, T. Gunji¹²², W. Guo⁶, A. Gupta⁹¹, R. Gupta⁹¹, R. Gupta⁴⁷, S. P. Guzman⁴⁴, K. Gwizdziel¹³³, L. Gyulai¹³⁶, M. K. Habib⁹⁷, C. Hadjidakis⁷², F. U. Haider⁹¹, H. Hamagaki⁷⁶, A. Hamdi⁷⁴, M. Hamid⁶,
Y. Han¹³⁸, R. Hannigan¹⁰⁸, J. Hansen⁷⁵, M. R. Haque¹³³, J. W. Harris¹³⁷, A. Harton⁹, H. Hassan⁸⁷, D. Hatzifotiadou⁵⁰, P. Hauer⁴², L. B. Havener¹³⁷, S. T. Heckel⁹⁵, E. Hellbär⁹⁷, H. Helstrup³⁴, M. Hemmer⁶³, T. Herman³⁵, G. Herrera Corral⁸, F. Herrmann¹³⁵, S. Herrmann¹²⁶, K. F. Hetland³⁴, B. Heybeck⁶³, H. Hillemanns³², B. Hippolyte¹²⁷, F. W. Hoffmann⁶⁹, B. Hofman⁵⁸, B. Hohlweger⁸⁴, G. H. Hong¹³⁸, M. Horst⁹⁵, A. Horzyk², Y. Hou⁶, P. Hristov³², C. Hughes¹²⁰, P. Huhn⁶³, L. M. Huhta¹¹⁵, T. J. Humanic⁸⁸, A. Hutson¹¹⁴, D. Hutter³⁸, R. Ilkaev¹⁴⁰, H. Ilyas¹³, M. Inaba¹²³, G. M. Innocenti³², M. Ippolitov¹⁴⁰, A. Isakov⁸⁶, T. Isidori¹¹⁶, M. S. Islam⁹⁹, M. Ivanov¹², M. Ivanov⁹⁷, V. Ivanov¹⁴⁰, M. Jablonski², B. Jacak⁷⁴, N. Jacazio³², P. M. Jacobs⁷⁴, S. Jadlovska¹⁰⁶, J. Jadlovsky¹⁰⁶, S. Jaelani⁸², C. Jahnke¹¹¹, M. J. Jakubowska¹³³, M. A. Janik¹³³, T. Janson⁶⁹, M. Jercic⁸⁹, S. Jia¹⁰, A. A. P. Jimenez⁶⁴, F. Jonas⁸⁷, J. M. Jowett ³²,⁹⁷, J. Jung⁶³, M. Jung⁶³, A. Junique³², A. Jusko¹⁰⁰, M. J. Kabus³²,¹³³, J. Kaewjai¹⁰⁵, P. Kalinak⁵⁹, A. S. Kalteyer⁹⁷, A. Kalweit³², V. Kaplin¹⁴⁰, A. Karasu Uysal⁷¹, D. Karatovic⁸⁹, O. Karavichev¹⁴⁰, T. Karavicheva¹⁴⁰, P. Karczmarczyk¹³³, E. Karpechev¹⁴⁰, U. Kebschull⁶⁹, R. Keidel¹³⁹, D. L. D. Keijdener⁵⁸, M. Keil³², B. Ketzer⁴², S. S. Khade⁴⁷, A. M. Khan⁶,¹¹⁸, S. Khan¹⁵, A. Khanzadeev¹⁴⁰, Y. Kharlov¹⁴⁰, A. Khatun¹⁵,¹¹⁶, A. Khuntia¹⁰⁷, M. B. Kidson¹¹³, B. Kileng³⁴, B. Kim¹⁰⁴, C. Kim¹⁶, D. J. Kim¹¹⁵, E. J. Kim⁶⁸, J. Kim¹³⁸, J. S. Kim⁴⁰, J. Kim⁶⁸, M. Kim¹⁸, S. Kim¹⁷, T. Kim¹³⁸, K. Kimura⁹², S. Kirsch⁶³, I. Kisel³⁸, S. Kiselev¹⁴⁰, A. Kisiel¹³³, J. P. Kitowski², J. L. Klay⁵, J. Klein³², S. Klein⁷⁴, C. Klein-Bösing¹³⁵, M. Kleiner⁶³, T. Klemenz⁹⁵, A. Kluge³², A. G. Knospe¹¹⁴, C. Kobdaj¹⁰⁵, T. Kollegger⁹⁷, A. Kondratyev¹⁴¹, N. Kondratyeva¹⁴⁰, E. Kondratyuk¹⁴⁰, J. König⁶³, S. A. Königstorfer⁹⁵, P. J. Konopka³², G. Kornakov¹³³, S. D. Koryciak², A. Kotliarov⁸⁶, V. Kovalenko¹⁴⁰, M. Kowalski¹⁰⁷, V. Kozhuharov³⁶, I. Králik⁵⁹, A. Kravčáková³⁷, L. Krcal³²,³⁸, M. Krivda⁵⁹,¹⁰⁰, F. Krizek⁸⁶, K. Krizkova Gajdosova³², M. Kroesen⁹⁴, M. Krüger⁶³, D. M. Krupova³⁵, E. Kryshen¹⁴⁰, V. Kučera⁵⁷, C. Kuhn¹²⁷, P. G. Kuijer⁸⁴, T. Kumaoka¹²³, D. Kumar¹³², L. Kumar⁹⁰, N. Kumar⁹⁰, S. Kumar³¹, S. Kundu³², P. Kurashvili⁷⁹, A. Kurepin¹⁴⁰, A. B. Kurepin¹⁴⁰, A. Kuryakin¹⁴⁰, S. Kushpil⁸⁶, J. Kvapil¹⁰⁰, M. J. Kweon⁵⁷, J. Y. Kwon⁵⁷, Y. Kwon¹³⁸, S. L. La Pointe³⁸, P. La Rocca²⁶, A. Lakrathok¹⁰⁵, M. Lamanna³², R. Langoy¹¹⁹, P. Larionov³², E. Laudi³², L. Lautner³²,⁹⁵, R. Lavicka¹⁰², T. Lazareva¹⁴⁰, R. Lea⁵⁴,¹³¹, H. Lee¹⁰⁴, I. Legrand⁴⁵, G. Legras¹³⁵, J. Lehrbach³⁸, T. M. Lelek², R. C. Lemmon⁸⁵, I. León Monzón¹⁰⁹, M. M. Lesch⁹⁵, E. D. Lesser¹⁸, P. Lévai¹³⁶, X. Li¹⁰, X. L. Li⁶, J. Lien¹¹⁹, R. Lietava¹⁰⁰, I. Likmeta¹¹⁴, B. Lim²⁴, S. H. Lim¹⁶, V. Lindenstruth³⁸, A. Lindner⁴⁵, C. Lippmann⁹⁷, A. Liu¹⁸, D. H. Liu⁶, J. Liu¹¹⁷, I. M. Lofnes²⁰, C. Loizides⁸⁷, S. Lokos¹⁰⁷, J. Lomker⁵⁸, P. Loncar³³, J. A. Lopez⁹⁴, X. Lopez¹²⁵, E. López Torres⁷, P. Lu⁹⁷,¹¹⁸, J. R. Luhder¹³⁵, M. Lunardon²⁷, G. Luparello⁵⁶, Y. G. Ma³⁹, A. Maevskaya¹⁴⁰, M. Mager³², A. Maire¹²⁷, M. V. Makariev³⁶, M. Malaev¹⁴⁰, G. Malfattore²⁵, N. M. Malik⁹¹, Q. W. Malik¹⁹, S. K. Malik⁹¹, L. Malinina¹⁴¹, e, D. Mal'Kevich¹⁴⁰, D. Mallick⁸⁰, N. Mallick⁴⁷, G. Mandaglio³⁰,⁵², S. K. Mandal⁷⁹, V. Manko¹⁴⁰, F. Manso¹²⁵, V. Manzari⁴⁹, Y. Mao⁶, G. V. Margagliotti²³, A. Margotti⁵⁰, A. Marín⁹⁷, C. Markert¹⁰⁸, P. Martinengo³², J. L. Martinez¹¹⁴, M. I. Martínez⁴⁴, G. Martínez García¹⁰³, M. P. P. Martins¹¹⁰, S. Masciocchi⁹⁷, M. Masera²⁴, A. Masoni⁵¹, L. Massacrier⁷², A. Mastroserio⁴⁹,¹²⁹, O. Matonoha⁷⁵, P. F. T. Matuoka¹¹⁰, A. Matyja¹⁰⁷, C. Mayer¹⁰⁷, A. L. Mazuecos³², F. Mazzaschi²⁴, M. Mazzilli³², J. E. Mdhluli¹²¹, A. F. Mechler⁶³, Y. Melikyan⁴³,¹⁴⁰, A. Menchaca-Rocha⁶⁶, E. Meninno²⁸,¹⁰², A. S. Menon¹¹⁴, M. Meres¹², S. Mhlanga⁶⁷,¹¹³, Y. Miake¹²³, L. Micheletti⁵⁵, L. C. Migliorin¹²⁶, D. L. Mihaylov⁹⁵, K. Mikhaylov¹⁴⁰,¹⁴¹, A. N. Mishra¹³⁶, D. Miśkowiec⁹⁷, A. Modak⁴, A. P. Mohanty⁵⁸, B. Mohanty⁸⁰, M. Mohisin Khan¹⁵, e, M. A. Molander⁴³, Z. Moravcova⁸³, C. Mordasini⁹⁵, D. A. Moreira De Godoy¹³⁵, I. Morozov¹⁴⁰, A. Morsch³², T. Mrnjavac³², V. Muccifora⁴⁸, S. Muhuri¹³², J. D. Mulligan⁷⁴, A. Mulliri²², M. G. Munhoz¹¹⁰, R. H. Munzer⁶³, H. Murakami¹²², S. Murray¹¹³, L. Musa³², J. Musinsky⁵⁹, J. W. Myrcha¹³³, B. Naik¹²¹, A. I. Nambrath¹⁸, B. K. Nandi⁴⁶, R. Nania⁵⁰, E. Nappi⁴⁹, A. F. Nassirpour¹⁷,⁷⁵, A. Nath⁹⁴, C. Nattrass¹²⁰, T. K. Nayak⁸⁰, M. N. Naydenov³⁶, A. Neagu¹⁹, A. Negru¹²⁴, L. Nellen⁶⁴, S. V. Nesbo³⁴, G. Neskovic³⁸, D. Nesterov¹⁴⁰, B. S. Nielsen⁸³, E. G. Nielsen⁸³, S. Nikolaev¹⁴⁰, S. Nikulin¹⁴⁰, V. Nikulin¹⁴⁰, F. Noferini⁵⁰, S. Noh¹¹, P. Nomokonov¹⁴¹, J. Norman¹¹⁷, N. Novitzky¹²³, P. Nowakowski¹³³, A. Nyanin¹⁴⁰, J. Nystrand²⁰, M. Ogino⁷⁶, A. Ohlson⁷⁵, V. A. Okorokov¹⁴⁰, J. Oleniacz¹³³, A. C. Oliveira Da Silva¹²⁰, M. H. Oliver¹³⁷, A. Onnerstad¹¹⁵, C. Oppedisano⁵⁵, A. Ortiz Velasquez⁶⁴, J. Otwinowski¹⁰⁷, M. Oya⁹², K. Oyama⁷⁶, Y. Pachmayer⁹⁴, S. Padhan⁴⁶, D. Pagano⁵⁴,¹³¹, G. Paić⁶⁴, A. Palasciano⁴⁹, S. Panebianco¹²⁸, H. Park¹²³, H. Park¹⁰⁴, J. Park⁵⁷, J. E. Parkkila³², R. N. Patra⁹¹, B. Paul²², H. Pei⁶, T. Peitzmann⁵⁸, X. Peng⁶, M. Pennisi²⁴, L. G. Pereira⁶⁵, D. Peresunko¹⁴⁰, G. M. Perez⁷, S. Perrin¹²⁸, Y. Pestov¹⁴⁰, V. Petráček³⁵, V. Petrov¹⁴⁰, M. Petrovici⁴⁵, R. P. Pezzi⁶⁵,¹⁰³, S. Piano⁵⁶, M. Pikna¹², P. Pillot¹⁰³, O. Pinazza³²,⁵⁰, L. Pinsky¹¹⁴, C. Pinto⁹⁵, S. Pisano⁴⁸, M. Płoskoń⁷⁴,
M. Planinic⁸⁹, F. Pliquett⁶³, M. G. Poghosyan⁸⁷, B. Polichtchouk¹⁴⁰, S. Politano²⁹, N. Poljak⁸⁹, A. Pop⁴⁵, S. Porteboeuf-Houssais¹²⁵, V. Pozdnikov¹⁴¹, I. Y. Pozos⁴⁴, K. K. Pradhan⁴⁷, S. K. Prasad⁴, S. Prasad⁴⁷, R. Preghenella⁵⁰, F. Prino⁵⁵, C. A. Pruneau¹³⁴, I. Pshenichnov¹⁴⁰, M. Puccio³², S. Pucillo²⁴, Z. Pugelova¹⁰⁶, S. Qiu⁸⁴, L. Quaglia²⁴, R. E. Quishpe¹¹⁴, S. Ragoni¹⁴, A. Rakotozafindrabe¹²⁸, L. Ramello⁵⁵,¹³⁰, F. Rami¹²⁷, S. A. R. Ramirez⁴⁴, T. A. Rancien⁷³, M. Rasa²⁶, S. S. Räsänen⁴³, R. Rath⁵⁰, M. P. Rauch²⁰, I. Ravasenga⁸⁴, K. F. Read⁸⁷,¹²⁰, C. Reckziegel¹¹², A. R. Redelbach³⁸, K. Redlich⁷⁹,¹⁰, C. A. Reetz⁹⁷, A. Rehman²⁰, F. Reidt³², H. A. Reme-Ness³⁴, Z. Rescakova³⁷, K. Reygers⁹⁴, A. Riabov¹⁴⁰, V. Riabov¹⁴⁰, R. Ricci²⁸, M. Richter¹⁹, A. A. Riedel⁹⁵, W. Riegler³², C. Ristea⁶², M. Rodríguez Cahuantzi⁴⁴, K. Røed¹⁹, R. Rogalev¹⁴⁰, E. Rogochaya¹⁴¹, T. S. Rogoschinski⁶³, D. Rohr³², D. Röhrich²⁰, P. F. Rojas⁴⁴, S. Rojas Torres³⁵, P. S. Rokita¹³³, G. Romanenko¹⁴¹, F. Ronchetti⁴⁸, A. Rosano³⁰,⁵², E. D. Rosas⁶⁴, K. Roslon¹³³, A. Rossi⁵³, A. Roy⁴⁷, S. Roy⁴⁶, N. Rubini²⁵, O. V. Rueda¹¹⁴, D. Ruggiano¹³³, R. Rui²³, P. G. Russek², R. Russo⁸⁴, A. Rustamov⁸¹, E. Ryabinkin¹⁴⁰, Y. Ryabov¹⁴⁰, A. Rybicki¹⁰⁷, H. Rytkonen¹¹⁵, W. Rzesa¹³³, O. A. M. Saarimaki⁴³, R. Sadek¹⁰³, S. Sadhu³¹, S. Sadovsky¹⁴⁰, J. Saetre²⁰, K. Šafařík³⁵, P. Saha⁴¹, S. K. Saha⁴, S. Saha⁸⁰, B. Sahoo⁴⁶, B. Sahoo⁴⁷, R. Sahoo⁴⁷, S. Sahoo⁶⁰, D. Sahu⁴⁷, P. K. Sahu⁶⁰, J. Saini¹³², K. Sajdakova³⁷, S. Sakai¹²³, M. P. Salvan⁹⁷, S. Sambyal⁹¹, I. Sanna³²,⁹⁵, T. B. Saramela¹¹⁰, D. Sarkar¹³⁴, N. Sarkar¹³², P. Sarma⁴¹, V. Sarritzu²², V. M. Sarti⁹⁵, M. H. P. Sas¹³⁷, J. Schambach⁸⁷, H. S. Scheid⁶³, C. Schiaua⁴⁵, R. Schicker⁹⁴, A. Schmah⁹⁴, C. Schmidt⁹⁷, H. R. Schmidt⁹³, M. O. Schmidt³², M. Schmidt⁹³, N. V. Schmidt⁸⁷, A. R. Schmier¹²⁰, R. Schotter¹²⁷, A. Schröter³⁸, J. Schukraft³², K. Schwarz⁹⁷, K. Schweda⁹⁷, G. Scioli²⁵, E. Scomparin⁵⁵, J. E. Seger¹⁴, Y. Sekiguchi¹²², D. Sekihata¹²², I. Selyuzhenkov⁹⁷,¹⁴⁰, S. Senyukov¹²⁷, J. J. Seo⁵⁷, D. Serebryakov¹⁴⁰, L. Šerkšnytė⁹⁵, A. Sevcenco⁶², T. J. Shaba⁶⁷, A. Shabetai¹⁰³, R. Shahoyan³², A. Shangaraev¹⁴⁰, A. Sharma⁹⁰, B. Sharma⁹¹, D. Sharma⁴⁶, H. Sharma⁵³,¹⁰⁷, M. Sharma⁹¹, S. Sharma⁷⁶, S. Sharma⁹¹, U. Sharma⁹¹, A. Shatay⁷², O. Sheibani¹¹⁴, K. Shigaki⁹², M. Shimomura⁷⁷, J. Shin¹¹, S. Shirinkin¹⁴⁰, Q. Shou³⁹, Y. Sibiriak¹⁴⁰, S. Siddhanta⁵¹, T. Siemiarczuk⁷⁹, T. F. Silva¹¹⁰, D. Silvermyr⁷⁵, T. Simantathammakul¹⁰⁵, R. Simeonov³⁶, B. Singh⁹¹, B. Singh⁹⁵, R. Singh⁸⁰, R. Singh⁹¹, R. Singh⁴⁷, S. Singh¹⁵, V. K. Singh¹³², V. Singhal¹³², T. Sinha⁹⁹, B. Sitar¹², M. Sitta⁵⁵,¹³⁰, T. B. Skaali¹⁹, G. Skorodumov⁹⁴, M. Slupecki⁴³, N. Smirnov¹³⁷, R. J. M. Snellings⁵⁸, E. H. Solheim¹⁹, J. Song¹¹⁴, A. Songmoolnak¹⁰⁵, C. Sonnabend³²,⁹⁷, F. Soramel²⁷, A. B. Soto-hernandez⁸⁸, R. Spijkers⁸⁴, I. Sputowska¹⁰⁷, J. Staa⁷⁵, J. Stachel⁹⁴, I. Stan⁶², P. J. Steffanic¹²⁰, S. F. Stiefelmaier⁹⁴, D. Stocco¹⁰³, I. Storehaug¹⁹, P. Stratmann¹³⁵, S. Strazzi²⁵, C. P. Stylianidis⁸⁴, A. A. P. Suaide¹¹⁰, C. Suire⁷², M. Sukhanov¹⁴⁰, M. Suljic³², R. Sultanov¹⁴⁰, V. Sumberia⁹¹, S. Sumowidagdo⁸², S. Swain⁶⁰, I. Szarka¹², M. Szymkowski¹³³, S. F. Taghavi⁹⁵, G. Taillepied⁹⁷, J. Takahashi¹¹¹, G. J. Tambave⁸⁰, S. Tang⁶,¹²⁵, Z. Tang¹¹⁸, J. D. Tapia Takaki¹¹⁶, N. Tapus¹²⁴, L. A. Tarasovicova¹³⁵, M. G. Tarzila⁴⁵, G. F. Tassielli³¹, A. Tauro³², G. Tejeda Muñoz⁴⁴, A. Telesca³², L. Terlizzi²⁴, C. Terrevoli¹¹⁴, S. Thakur⁴, D. Thomas¹⁰⁸, A. Tikhonov¹⁴⁰, A. R. Timmins¹¹⁴, M. Tkacik¹⁰⁶, T. Tkacik¹⁰⁶, A. Toia⁶³, R. Tokumoto⁹², N. Topilskaya¹⁴⁰, M. Toppi⁴⁸, F. Torales-Acosta¹⁸, T. Tork⁷², A. G. Torres Ramos³¹, A. Trifiró³⁰,⁵², A. S. Triolo³⁰,³²,⁵², S. Tripathy⁵⁰, T. Tripathy⁴⁶, S. Trogolo³², V. Trubnikov³, W. H. Trzaska¹¹⁵, T. P. Trzcinski¹³³, A. Tumkin¹⁴⁰, R. Turrisi⁵³, T. S. Tveter¹⁹, K. Ullaland²⁰, B. Ulukutlu⁹⁵, A. Uras¹²⁶, M. Urioni⁵⁴,¹³¹, G. L. Usai²², M. Vala³⁷, N. Valle²¹, L. V. R. van Doremalen⁵⁸, M. van Leeuwen⁸⁴, C. A. van Veen⁹⁴, R. J. G. van Weelden⁸⁴, P. Vande Vyvre³², D. Varga¹³⁶, Z. Varga¹³⁶, M. Vasileiou⁷⁸, A. Vasiliev¹⁴⁰, O. Vázquez Doce⁴⁸, V. Vechernin¹⁴⁰, E. Vercellin²⁴, S. Vergara Limón⁴⁴, L. Vermunt⁹⁷, R. Vértesi¹³⁶, M. Verweij⁵⁸, L. Vickovic³³, Z. Vilakazi¹²¹, O. Villalobos Baillie¹⁰⁰, A. Villani²³, G. Vino⁴⁹, A. Vinogradov¹⁴⁰, T. Virgili²⁸, M. M. O. Virta¹¹⁵, V. Vislavicius⁷⁵, A. Vodopyanov¹⁴¹, B. Volkel³², M. A. Völkl⁹⁴, K. Voloshin¹⁴⁰, S. A. Voloshin¹³⁴, G. Volpe³¹, B. von Haller³², I. Vorobyev⁹⁵, N. Vozniuk¹⁴⁰, J. Vrláková³⁷, C. Wang³⁹, D. Wang³⁹, Y. Wang³⁹, A. Wegrzynek³², F. T. Weiglhofer³⁸, S. C. Wenzel³², J. P. Wessels¹³⁵, S. L. Weyhmiller¹³⁷, J. Wiechula⁶³, J. Wikne¹⁹, G. Wilk⁷⁹, J. Wilkinson⁹⁷, G. A. Willems¹³⁵, B. Windelband⁹⁴, M. Winn¹²⁸, J. R. Wright¹⁰⁸, W. Wu³⁹, Y. Wu¹¹⁸, R. Xu⁶, A. Yadav⁴², A. K. Yadav¹³², S. Yalcin⁷¹, Y. Yamaguchi⁹², S. Yang²⁰, S. Yano²⁰, Z. Yin⁶, I.-K. Yoo¹⁶, J. H. Yoon⁵⁷, S. Yuan²⁰, A. Yuncu⁹⁴, V. Zaccolo²³, C. Zampolli³², F. Zanone⁹⁴, N. Zardoshti³², A. Zarochentsev¹⁴⁰, P. Závada⁶¹, N. Zaviyalov¹⁴⁰, M. Zhalov¹⁴⁰, B. Zhang⁶, L. Zhang³⁹, S. Zhang³⁹, X. Zhang⁶, Y. Zhang¹¹⁸, Z. Zhang⁶, M. Zhao¹⁰, V. Zherebchevskii¹⁴⁰, Y. Zhi¹⁰, D. Zhou⁶, Y. Zhou⁸³, J. Zhu⁶,⁹⁷, Y. Zhu⁶, S. C. Zugravel⁵⁵, N. Zurlo⁵⁴,¹³¹

<footnote n="1">¹ A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, Yerevan, Armenia</footnote>
² AGH University of Science and Technology, Cracow, Poland

³ Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, Kiev, Ukraine

⁴ Department of Physics, Centre for Astroparticle Physics and Space Science (CAPSS), Bose Institute, Kolkata, India

⁵ California Polytechnic State University, San Luis Obispo, CA, USA

⁶ Central China Normal University, Wuhan, China

⁷ Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Havana, Cuba

⁸ Centro de Investigación y de Estudios Avanzados (CINVESTAV), Mexico City and Mérida, Mexico

⁹ Chicago State University, Chicago, IL, USA

¹⁰ China Institute of Atomic Energy, Beijing, China

¹¹ Chungbuk National University, Cheongju, Republic of Korea

¹² Faculty of Mathematics, Physics and Informatics, Comenius University Bratislava, Bratislava, Slovak Republic

¹³ COMSATS University Islamabad, Islamabad, Pakistan

¹⁴ Creighton University, Omaha, NE, USA

¹⁵ Department of Physics, Aligarh Muslim University, Aligarh, India

¹⁶ Department of Physics, Pusan National University, Pusan, Republic of Korea

¹⁷ Department of Physics, Sejong University, Seoul, Republic of Korea

¹⁸ Department of Physics, University of California, Berkeley, CA, USA

¹⁹ Department of Physics, University of Oslo, Oslo, Norway

²⁰ Department of Physics and Technology, University of Bergen, Bergen, Norway

²¹ Dipartimento di Fisica, Università di Pavia, Pavia, Italy

²² Dipartimento di Fisica dell'Università and Sezione INFN, Cagliari, Italy

²³ Dipartimento di Fisica dell'Università and Sezione INFN, Trieste, Italy

²⁴ Dipartimento di Fisica dell'Università and Sezione INFN, Turin, Italy

²⁵ Dipartimento di Fisica e Astronomia dell'Università and Sezione INFN, Bologna, Italy

²⁶ Dipartimento di Fisica e Astronomia dell'Università and Sezione INFN, Catania, Italy

²⁷ Dipartimento di Fisica e Astronomia dell'Università and Sezione INFN, Padua, Italy

²⁸ Dipartimento di Fisica 'E.R. Caianiello' dell'Università and Gruppo Collegato INFN, Salerno, Italy

²⁹ Dipartimento DISAT del Politecnico and Sezione INFN, Turin, Italy

³⁰ Dipartimento di Scienze MIFT, Università di Messina, Messina, Italy

³¹ Dipartimento Interateneo di Fisica 'M. Merlin' and Sezione INFN, Bari, Italy

³² European Organization for Nuclear Research (CERN), Geneva, Switzerland

³³ Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, University of Split, Split, Croatia

³⁴ Faculty of Engineering and Science, Western Norway University of Applied Sciences, Bergen, Norway

³⁵ Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Prague, Czech Republic

³⁶ Faculty of Physics, Sofia University, Sofia, Bulgaria

³⁷ Faculty of Science, P.J. Šafárik University, Kosice, Slovak Republic

³⁸ Frankfurt Institute for Advanced Studies, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt, Germany

³⁹ Fudan University, Shanghai, China

⁴⁰ Gangneung-Wonju National University, Gangneung, Republic of Korea

⁴¹ Department of Physics, Gauhati University, Guwahati, India

⁴² Helmholtz-Institut für Strahlen- und Kernphysik, Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn, Germany

⁴³ Helsinki Institute of Physics (HIP), Helsinki, Finland

⁴⁴ High Energy Physics Group, Universidad Autónoma de Puebla, Puebla, Mexico

⁴⁵ Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest, Romania

⁴⁶ Indian Institute of Technology Bombay (IIT), Mumbai, India

⁴⁷ Indian Institute of Technology Indore, Indore, India

⁴⁸ INFN, Laboratori Nazionali di Frascati, Frascati, Italy

⁴⁹ INFN, Sezione di Bari, Bari, Italy

⁵⁰ INFN, Sezione di Bologna, Bologna, Italy

⁵¹ INFN, Sezione di Cagliari, Cagliari, Italy

⁵² INFN, Sezione di Catania, Catania, Italy

⁵³ INFN, Sezione di Padova, Padua, Italy

⁵⁴ INFN, Sezione di Pavia, Pavia, Italy
55 INFN, Sezione di Torino, Turin, Italy

56 INFN, Sezione di Trieste, Trieste, Italy

57 Inha University, Incheon, Republic of Korea

58 Institute for Gravitational and Subatomic Physics (GRASP), Utrecht University/Nikhef, Utrecht, The Netherlands

59 Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak Republic

60 Institute of Physics, Homi Bhabha National Institute, Bhubaneswar, India

61 Institute of Physics of the Czech Academy of Sciences, Prague, Czech Republic

62 Institute of Space Science (ISS), Bucharest, Romania

63 Institut für Kernphysik, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt, Germany

64 Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Mexico City, Mexico

65 Instituto de Física, Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, Brazil

66 Instituto de Física, Universidad Nacional Autónoma de México, Mexico City, Mexico

67 iThemba LABS, National Research Foundation, Somerset West, South Africa

68 Jeonbuk National University, Jeonju, Republic of Korea

69 Johann-Wolfgang-Goethe Universität Frankfurt Institut für Informatik, Fachbereich Informatik und Mathematik, Frankfurt, Germany

70 Korea Institute of Science and Technology Information, Daejeon, Republic of Korea

71 KTO Karatay University, Konya, Turkey

72 Laboratoire de Physique des 2 Infinis, Irène Joliot-Curie, Orsay, France

73 Laboratoire de Physique Subatomique et de Cosmologie, CNRS-IN2P3, Université Grenoble-Alpes, Grenoble, France

74 Lawrence Berkeley National Laboratory, Berkeley, CA, USA

75 Division of Particle Physics, Department of Physics, Lund University, Lund, Sweden

76 Nagasaki Institute of Applied Science, Nagasaki, Japan

77 Nara Women's University (NWU), Nara, Japan

78 Department of Physics, School of Science, National and Kapodistrian University of Athens, Athens, Greece

79 National Centre for Nuclear Research, Warsaw, Poland

80 National Institute of Science Education and Research, Homi Bhabha National Institute, Jatni, India

81 National Nuclear Research Center, Baku, Azerbaijan

82 National Research and Innovation Agency-BRIN, Jakarta, Indonesia

83 Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark

84 Nikhef, National Institute for Subatomic Physics, Amsterdam, The Netherlands

85 Nuclear Physics Group, STFC Daresbury Laboratory, Daresbury, UK

86 Nuclear Physics Institute of the Czech Academy of Sciences, Husinec-Řež, Czech Republic

87 Oak Ridge National Laboratory, Oak Ridge, TN, USA

88 Ohio State University, Columbus, OH, USA

89 Physics Department, Faculty of science, University of Zagreb, Zagreb, Croatia

90 Physics Department, Panjab University, Chandigarh, India

91 Physics Department, University of Jammu, Jammu, India

92 Physics Program and International Institute for Sustainability with Knotted Chiral Meta Matter (SKCM2), Hiroshima University, Hiroshima, Japan

93 Physikalisches Institut, Eberhard-Karls-Universität Tübingen, Tübingen, Germany

94 Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany

95 Physik Department, Technische Universität München, Munich, Germany

96 Politecnico di Bari and Sezione INFN, Bari, Italy

97 Research Division and ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany

98 Saga University, Saga, Japan

99 Saha Institute of Nuclear Physics, Homi Bhabha National Institute, Kolkata, India

100 School of Physics and Astronomy, University of Birmingham, Birmingham, UK

101 Sección Física, Departamento de Ciencias, Pontificia Universidad Católica del Perú, Lima, Peru

102 Stefan Meyer Institut für Subatomare Physik (SMI), Vienna, Austria

103 SUBATECH, IMT Atlantique, CNRS-IN2P3, Nantes Université, Nantes, France

104 Sungkyunkwan University, Suwon City, Republic of Korea
105 Suranaree University of Technology, Nakhon Ratchasima, Thailand

106 Technical University of Košice, Kosice, Slovak Republic

107 The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Poland

108 The University of Texas at Austin, Austin, TX, USA

109 Universidad Autónoma de Sinaloa, Culiacán, Mexico

110 Universidade de São Paulo (USP), São Paulo, Brazil

111 Universidade Estadual de Campinas (UNICAMP), Campinas, Brazil

112 Universidade Federal do ABC, Santo Andre, Brazil

113 University of Cape Town, Cape Town, South Africa

114 University of Houston, Houston, TX, USA

115 University of Jyväskylä, Jyvaskyla, Finland

116 University of Kansas, Lawrence, KS, USA

117 University of Liverpool, Liverpool, UK

118 University of Science and Technology of China, Hefei, China

119 University of South-Eastern Norway, Kongsberg, Norway

120 University of Tennessee, Knoxville, TN, USA

121 University of the Witwatersrand, Johannesburg, South Africa

122 University of Tokyo, Tokyo, Japan

123 University of Tsukuba, Tsukuba, Japan

124 University Politehnică of Bucharest, Bucharest, Romania

125 CNRS/IN2P3, LPC, Université Clermont Auvergne, Clermont-Ferrand, France

126 CNRS/IN2P3, Institut de Physique des 2 Infinis de Lyon, Université de Lyon, Lyon, France

127 CNRS, IPHC UMR 7178, Université de Strasbourg, 67000 Strasbourg, France

128 Département de Physique Nucléaire (DPhN), IRFU, Université Paris-Saclay Centre d'Etudes de Saclay (CEA), Saclay, France

129 Università degli Studi di Foggia, Foggia, Italy

130 Università del Piemonte Orientale, Vercelli, Italy

131 Università di Brescia, Brescia, Italy

132 Variable Energy Cyclotron Centre, Homi Bhabha National Institute, Kolkata, India

133 Warsaw University of Technology, Warsaw, Poland

134 Wayne State University, Detroit, MI, USA

135 Institut für Kernphysik, Westfälische Wilhelms-Universität Münster, Münster, Germany

136 Wigner Research Centre for Physics, Budapest, Hungary

137 Yale University, New Haven, CT, USA

138 Yonsei University, Seoul, Republic of Korea

139 Zentrum für Technologie und Transfer (ZTT), Worms, Germany

140 Affiliated with an Institute Covered by a Cooperation Agreement with CERN, Geneva, Switzerland

141 Affiliated with an International Laboratory Covered by a Cooperation Agreement with CERN, Geneva, Switzerland

a Also at: Max-Planck-Institut für Physik, Munich, Germany

b Also at: Italian National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA), Bologna, Italy

c Also at: Department of Applied Physics, Aligarh Muslim University, Aligarh, India

d Also at: Institute of Theoretical Physics, University of Wroclaw, Wroclaw, Poland

e Also at: An Institution Covered by a Cooperation Agreement with CERN, Geneva, Switzerland

f Deceased