Inclusive photon production at forward rapidities in pp and p–Pb collisions at √ s NN = 5 . 02 TeV

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 √ s NN = 5 . 02 TeV using the ALICE detector in the forward pseudorapidity region 2.3 < η lab < 3.9 is presented. Measurements in p–Pb collisions are reported for two beam conﬁgurations 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 whicharedeﬁnedbasedondifferenteventactivityestimators, i.e., the charged-particle multiplicity measured at midrapid-ity 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 classesinp–Pbcollisions.Theresultsarecomparedtopredic-tions 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.


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][2][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][7][8][9][10][11][12][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 (η) density of produced particles (dN/dη) 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 (π 0 ) [16][17][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 at central and forward rapidities.Charged-particle multiplicity measurements at the LHC were previously performed by ALICE, ATLAS, CMS, and LHCb experiments in pp, p-Pb, and Pb-Pb collisions .Inclusive photon production at forward rapidity was studied in pp collisions at √ 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 √ s [20][21][22] whereas the mean multiplicity of inclusive photons at forward rapidity can be described by both a logarithmic and a power law dependence with √ 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 √ 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 √ 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 Sec. 2. The data samples used in this analysis and the selection of events are discussed in Sec. 3. Section 4 discusses the details of the used centrality estimators.A discussion on Monte Carlo (MC) event generators and the simulation framework is given in Sec. 5.The reconstruction of inclusive photons is presented in Sec. 6.The correction for instrumental effects using the unfolding method is described in Sec. 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 Sec.10.

Experimental setup
A comprehensive description of the ALICE detectors and their performances can be found in Refs.[44,45].The Photon 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 reference 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 < η lab < 3.9 with full azimuth coverage.The PMD consists of two fine granular planes: the Charged Particle Veto plane and the preshower (PRE) plane.A lead converter with a thickness of three radiation lengths (3X 0 ) is placed between these two planes.The PMD consists of 184,320 hexagonal cells of size 0.22 cm 2 and depth of 0.5 cm with a copper honey-comb shaped cathode extended towards a 20 µ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 presence of the lead converter, the incident photons produce electromagnetic 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 Sec. 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] surrounding the beam pipe.The two layers cover the pseudorapidity ranges |η lab | < 2 (inner layer) and |η lab | < 1.4 (outer layer).For this analysis, the SPD is mainly used in the determination of the position of the interaction vertex, the selection of the minimum bias (MB) trigger (combined with V0), and the centrality estimation [25][26][27].The V0 detector [49,50] is made of two scintillator arrays placed on either side of the IP at z = 330 cm (V0A, covering the pseudorapidity interval 2.8 < η lab < 5.1), and z = −90 cm (V0C, covering −3.7 < η 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 ± 112.5 m from the IP, on both sides.The ZPA and ZPC calorimeters cover the pseudorapidity range 6.5 < |η lab | < 7.4, while the ZNA and ZNC calorimeters have a geometric coverage |η 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].

Data sample and event selection
The p-Pb and pp data samples used in this analysis were collected 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 proton 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 √ 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 ∆y = 0.465 with respect to symmetric pp collisions in the direction of the proton beam.In the following, the variable η 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 ± 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.

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 heavyion 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 (SNM) 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: : The yield of high-p T charged particles at midrapidity is proportional to the number of inelastic nucleon-nucleon collisions (N coll ); N Pb-side part : The charged-particle multiplicity in the Pb-going direction is proportional to the number of wounded target nucleons.

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 → 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, photonhadron, 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.

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 γ-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 γ-like clusters.The N γ-like clusters provide a photonrich 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 γ-like clusters are therefore required to be corrected for these effects as discussed in the next section.

Correction procedure
The measured distributions of N γ-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 γ 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 Pomerons.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 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 γ-true ) from the event generator and the reconstructed photon multiplicity (N γ-like ) in the simulated events.The mean ( N γ-like ) and width (σ N γ-like ) of the distribution of N γ-like are presented in  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 η 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 One can therefore recover the true distribution (T ) for given M as However, the matrix R mt may be singular and one can calculate inverse matrix R −1 mt 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 RooUnfold software package [74] is used to correct the measured N γ-like distributions following the description outlined in Refs.[16,[20][21][22][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 γ-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 γ-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 fluctuations 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 = U SV 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 (ε Trig ) and vertex reconstruction (ε Vtx ) efficiencies.A similar procedure as discussed in [16] is adopted to estimate ε Trig and ε Vtx using MC simulations.These efficiencies can be defined as 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 ε Trig and ε 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].

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.

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 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 η 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 √ s = 0.9, 2.76, and 7 TeV [16].

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.

Unfolding procedure
As discussed in Sec. 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.

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.

Diffraction mass distributions
As described in Sec. 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.

Results and discussions
The multiplicity and pseudorapidity distributions of inclusive photons in the forward pseudorapidity interval 2.3 < η lab < 3.9 measured using the PMD for pp collisions at √ s = 5.02 TeV and p-Pb and Pb-p collisions at √ 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.

Multiplicity distributions
The measured multiplicity distributions of inclusive photons are presented as probability distributions (P(N γ )) as a function of N γ 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 < η lab < 3.9.The average photon multiplicity for pp collisions is 6.44 ± 0.36 (sys), for p-Pb collisions is 29.38 ± 1.88 (sys), and for Pb-p collisions is 24.59 ± 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 ∆y shift of ± 0.465 units between p-Pb and Pb-p collisions.
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 √ 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 α soft parameter gives the fraction of soft events.The NBD dis- tribution, P NBD , has the following parameters: the average multiplicity n , which is found to increase with increasing √ s, and the shape parameter, k, which decreases with increasing √ s in hadronic colli- sions [16,22,23]. where (5) Double NBDs do not describe the value P(0), therefore, the bin N γ = 0 was excluded from the fit and a normalization factor λ 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 n semihard ≈ 3 n soft which is consistent with the results of chargedparticle production at central and forward rapidities in pp collisions [22,23].

Pseudorapidity distributions in NSD and INEL events
The pseudorapidity distributions of inclusive photons measured within 2.3 < η lab < 3.9 in pp, p-Pb, and Pb-p collisions are presented in Fig. 5.The gray bands represent 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 η 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 Sec.9.1.The data are compared to several models with different descriptions of particle production, all shifted by ∆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.
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 √ 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 √ 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 N γ within 2.3 < η lab < 3.9 as a function of √ s is shown in Fig. 6(b).The results of the measurements of N γ are given in Table 6. 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].
Table 6: The average photon multiplicity N γ 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]. Figure 7: The pseudorapidity distributions of inclusive photons are compared with charged-particle measurements at midrapidity by ALICE [79] in pp collisions (a) and both ALICE [24] and CMS [41] in p-Pb collisions at √ s NN = 5.02 TeV (b).Mirrored data points of pp and Pb-p results and predictions from various MC models are superimposed.
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 pseudorapidities, 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/dy and dN/dη at η ≈ 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.

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 evolution of dN γ /dη 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η 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 γ /dη lab with HIJING predictions for the three considered centrality classes.The dN ch /dη 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η lab for the CL1 (ZNA) estimator.The model describes the values of dN γ /dη 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 N γ at forward rapidity (2.3 < η lab < 3.9) as a function of N part for the two centrality estimators in p-Pb collisions at √ 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 N part for the different centrality classes for both the estimators are taken from [25].The data point for pp collisions at √ s = 5.02 TeV is also included for reference.The average photon multiplicity N γ divided by the average number of participants is presented as a function of N part in Fig. 11(b).The N γ / N part ratio has a steeper increase with N part for the CL1 centrality estimator than for 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 high-p T part or N Pb-side part ), the inclusive photon multiplicity is found to scale with N part within uncertainties and to point towards the pp data point at low N part .Moreover, the range in N part for the ZNA centrality selection is more limited (that could be because of saturation of forward neutron emission) than what is obtained by particlemultiplicity-based centrality estimators.This effect is also emphasized in Fig. 11(c) where the same quantity N γ / N part is presented as a function of N γ .Similar results were observed for chargedparticle multiplicity at midrapidity reported by ALICE in Ref. [25].

Summary
In summary, we have presented results for inclusive photon multiplicity and pseudorapidity distributions at forward rapidities (2.3 < η lab < 3.9) in pp collisions at √ s = 5.02 TeV and p-Pb and Pb-p collisions at √ 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 N γ with √ s in pp collisions and with N part 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 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 η 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 η and y is negligible.
An asymmetric dN ch /dη 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 N γ divided by the average number of participants scales linearly with both N high-p T part and N Pb-side part estimated by the energy deposited in the ZDC and the trend with N part is consistent with the data point from pp collisions at √ 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.
Figure1shows 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 γ-true ) from the event generator and the reconstructed photon multiplicity (N γ-like ) in the simulated events.The mean ( N γ-like ) and width (σ N γ-like ) of the distribution of N γ-like are presented in Fig. A.1 and Fig. A.2 as a function of N γ-true for pp and p-Pb collisions at √ s NN = 5.02 TeV.The resolution of photon multiplicity reconstruction in pp (p-Pb) collisions at N γ-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 η 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.(1)

Figure 1 :
Figure 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 γ-like and N γ-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.

Figure 2 :
Figure 2: MC closure test of the unfolding method for pp (a), p-Pb (b), and Pb-p (c) collisions at √ 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.

Figure 3 :
Figure 3: Top panels: Inclusive photon multiplicity distributions measured within 2.3 < η lab < 3.9 in pp (a), p-Pb (b), and Pb-p (c) collisions at √ 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.

Figure 4 :
Figure 4: Top panels: Multiplicity distribution of inclusive photons fitted to double NBDs for pp (a), p-Pb (b), and Pb-p (c) collisions at √ s NN = 5.02 TeV.Bottom panels: The ratio of the data to the fit is presented.Shaded boxes represent the systematic uncertainties.

Figure 5 :
Figure 5: Top panels: Pseudorapidity distribution of inclusive photons measured within 2.3 < η lab < 3.9 in pp (a), p-Pb (b), and Pb-p (c) collisions at √ 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.

Figure 6 :
Figure6: (a) Comparison of dN γ /dη lab measured at forward rapidity in INEL pp collisions for various center-ofmass 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].
the dominant contribution to inclusive photon production comes from π 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

Figure 8 :
Figure 8: Pseudorapidity distributions of inclusive photons measured within 2.3 < η lab < 3.9 in p-Pb collisions at √ s NN = 5.02 TeV for several centrality classes and for two centrality estimators: CL1 (top) and ZNA (bottom).

Figure 9 :
Figure 9: Ratios of dN γ /dη 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.

Figure 10 :
Figure10: Pseudorapidity distribution of inclusive photons measured within 2.3 < η lab < 3.9 in p-Pb collisions at √ 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).

Figure 11 :Figure A. 1 :Figure A. 2 :
Figure 11: (a) The average photon multiplicity N γ within 2.3 < η lab < 3.9 as a function of the average number of participants for various centrality estimators in p-Pb collisions at √ s NN = 5.02 TeV.Values of N γ / N part are shown as a function of N part (b) and N γ (c).The lines between the data points are to guide the eye.The data point for pp collisions at √ s = 5.02 TeV is also included for reference.

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.

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, numbers are given at multiplicity values of 0, the mean m and when P(N γ ) = 10 -3 .

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.

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.

Table 7 :
[25]average photon multiplicity N γ within 2.3 < η lab < 3.9 for various centrality classes defined using CL1 and ZNA centrality estimators in p-Pb collisions at √ s NN = 5.02 TeV.The corresponding values of N part are taken from Ref.[25].The quoted errors are systematic uncertainties.Statistical uncertainties are negligible.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.