pi(0) and eta meson production in proton-proton collisions at root s=8 TeV

An invariant differential cross section measurement of inclusive π 0 and η meson production at mid-rapidity in pp collisions at √ s = 8 TeV was carried out by the ALICE experiment at the LHC. The spectra of π 0 and η mesons were measured in transverse momentum ranges of 0 . 3 < p T < 35 GeV/c and 0 . 5 < p T < 35 GeV/c, respectively. Next-to-leading order perturbative QCD calculations using fragmentation functions DSS14 for the π 0 and AESSS for the η overestimate the cross sections of both neutral mesons, although such calculations agree with the measured η/π 0 ratio within uncertainties. The results were also compared with PYTHIA 8.2 predictions for which the Monash 2013 tune yields the best agreement with the measured neutral meson spectra. The measurements conﬁrm a universal behavior of the η/π 0 ratio seen for NA27, PHENIX and ALICE data for pp collisions from √ s = 27 . 5 GeV to √ s = 8 TeV within experimental uncertainties. A relation between the π 0 and η production cross sections for pp collisions at √ s = 8 TeV is given by m T scaling for p T > 3 . 5 GeV/c. However, a deviation from this empirical scaling rule is observed for transverse momenta below p T < 3 . 5 GeV/c in the η/π 0 ratio with a signiﬁcance of 6 . 2 σ .


Introduction
Measuring identified particle production in proton-proton (pp) collisions over wide kinematic ranges is considered an informative probe of strong interactions at high energies. Quantum chromodynamics (QCD) is the fundamental theory of the strong interaction [1]. It succeeds in providing a qualitative description of a wide range of phenomena in hadronic collisions. At typical hadron collider energies its perturbative expansion (pQCD) permits a detailed quantitative comparison with experimental data. However, it remains a challenge to provide a consistent description of hadron spectra at all collision energies reached experimentally. In theoretical models, particle production is usually divided into e-mail: alice-publications@cern.ch two categories: the "soft" scattering regime describing particle production involving small momentum transfers and the "hard" scattering regime, responsible for producing particles with momenta of several GeV/c or more.
Only "hard" scattering processes with a sufficiently large transverse momentum transfer, Q 2 , can be calculated using methods based on pQCD. High-momentum particles originate from the fragmentation of partons produced in scattering processes with large Q 2 . The theoretical description of a "hard" scattering process can be factorized into parton distribution functions (PDFs), the QCD matrix element and fragmentation functions (FFs). PDFs describe the fraction of the proton's longitudinal momentum carried by a scattered parton, x, and FFs describe the ratio of the observed hadron momentum to the final-state parton momentum, z, respectively. Comprehensive parametrizations of PDFs and FFs are derived from global fits to the experimental data at various collision energies. The energies reached at the LHC [2] open up the domains in x and z not accessible at lower energy. In the past, experiments at the LHC consequently found discrepancies between the measured π 0 and η meson spectra [3][4][5] and pQCD calculations based on fragmentation functions, which include mostly data from experiments below the TeV scale [6]. Since the gluon contribution becomes more dominant with increased center of mass energy, √ s [7], π 0 and η meson spectra at LHC energies provide new constraints on the gluon to light-flavor hadron fragmentation functions. Recent progress in comprehensive global QCD analysis of parton-to-pion fragmentation functions at next-to-leading order (NLO) [8] derived from inclusive pion production in semi-inclusive electronpositron annihilation, deep-inelastic scattering and pp collisions over a wide energy range, including the LHC results [3], achieves a good and consistent description of pion spectra, including the latest measurements of π 0 and η spectra in pp collisions at √ s = 2.76 TeV [9] and 7 TeV [3]. One of the conclusions of that analysis was that meson production from gluon fragmentation is reduced, which turns out to be at tension with previously available data obtained at RHIC [10]. In the quark model, the π 0 consists of light-flavor quarkantiquark pairs, uū and dd, whereas the η additionally contains hidden strangeness, ss. Measurements of both neutral mesons are thus of particular interest due to their different quark content as they help to constrain the PDFs and FFs [11] of the s quark.
The majority of particles at low transverse momenta, p T , are produced in "soft" processes involving a small Q 2 . In this regime, the pQCD calculations are not applicable for description of the production mechanisms and phenomenological models are based on previous measurements of neutral meson production cross sections or other light mesons by other experiments at lower collision energies. Particle production measurements at transverse momenta down to a few hundred MeV/c, as reported here, are particularly important to further constrain such models.
The importance of precise identified particle production measurements is underlined by various empirical rules observed in relative particle yields which allow estimates of the hadronic background of rare probes such as direct photons, dileptons and heavy-quark production. Almost all lower-energy experiments from ISR to RHIC report the observation of such an empirical rule, so-called m T scaling, in particle production over wide p T ranges [12,13]. The practical use of m T scaling is the ability to derive the p Tdependent differential yields of most of particles from the well measured light-flavor mesons, like pions and kaons, by assuming that the meson spectra can be described as a function of transverse mass m T : Ed 3 σ/d p 3 where the function f (m T ) is universal for all hadron species, so that their spectra share the same shape up to a normalization factor C h [14]. In the context of rare probes, this empirical relation is hence widely used to estimate the various background sources, for which no measurements are available. However, phenomenological analyses of new data delivered by the LHC experiments show that m T scaling is violated at higher p T compared to lower collision energies [14,15]. Therefore, precise measurements of identified hadron spectra over wide transverse momentum ranges at different LHC energies are of particular importance for the quantitative description of particle production at the LHC.
In this paper, the differential invariant production cross sections, Ed 3 σ/d p 3 , of π 0 and η mesons and the particle production ratio η/π 0 are presented, measured over wide p T ranges at mid-rapidity in pp collisions at √ s = 8 TeV by ALICE. The new experimental results are compared with pQCD calculations using MSTW08 (PDF) [16] with DSS14 (FF) [8] for the π 0 and accordingly CTEQ6M5 (PDF) [17] with AESSS (FF) [11] for the η, as well as the PYTHIA8.210 Monte Carlo (MC) event generator [18] with the tunes Tune 4C [19] and Monash 2013 [20]. This paper is organized as follows: In Sect. 2, the ALICE experiment is briefly described with the focus on the detectors used in this analysis, namely the calorimeters and the central tracking systems. Section 3 describes the datasets, the event selection and also introduces the calorimeter triggers used in this analysis. In Sect. 4, the reconstruction principles for neutral mesons are introduced. Furthermore, the determination of correction factors, which are used to calculate the differential invariant cross sections from the measured raw yields, is described. Section 5 discusses the various contributions to the statistical and systematic uncertainties of the measurements. In Sect. 6, the p T differential invariant cross sections for π 0 and η meson production in pp collisions at √ s = 8 TeV are presented and compared with pQCD calculations. Subsequently, the measured ratio of η/π 0 is presented and compared to the same theoretical models. Section 7 concludes the paper with a summary of the obtained results.

Detector description
Neutral mesons, π 0 and η, decay into photons, which are reconstructed via two fundamentally different detection methods. The first method exploits the measurement of photons using electromagnetic calorimeters. Two such calorimeters are available in ALICE [21,22]: the electromagnetic calorimeter (EMCal) [23] and the photon spectrometer (PHOS) [24]. The second method of photon detection makes use of photons converted into e + e − pairs within the inner detector material located between the interaction point and a radius which corresponds to the midpoint between the inner and outer field cage of the time projection chamber (TPC) [25]. These electron-positron pairs, originating at secondary vertices (V 0 ), are reconstructed by the main tracking systems in ALICE centered at mid-rapidity and consisting of the inner tracking system (ITS) [26] and the TPC [25]. The aforementioned detectors are described below, noting the detector configurations during pp data taking at √ s = 8 TeV in 2012.
The EMCal detector [23] is a sampling electromagnetic calorimeter. Its active elements, called cells, are composed of 77 alternating layers of lead and plastic scintillator providing a radiation length of 20.1 X 0 . The scintillation light in each layer is collected by wavelength shifting fibers perpendicular to the face of each cell. The fibers are connected to 5 × 5 mm 2 active area Avalanche photo diodes (APDs) to detect the generated scintillation light. Each cell has a size of η × φ = 0.0143 × 0.0143 (≈ 6.0 × 6.0 cm 2 ), corresponding to approximately twice the Molière radius. Groups of 2×2 cells are combined into modules, which are further combined into arrays of 12×24 modules called supermodules. In total, there are ten active, full EMCal supermodules, covering φ = 100 • in azimuth and |η| < 0.7 in pseudorapidity with a total number of 11,520 cells. The EMCal is located at a radial distance of 4.28 m at the closest point from the nominal collision vertex. The intrinsic energy resolution of the EMCal is parametrized as σ E /E = 4.8%/E ⊕ 11.3%/ √ E ⊕ 1.7% with E in units of GeV [27]. The relative energy calibration of the detector is performed by measuring, in each cell, the reconstructed π 0 mass in the invariant mass distribution of photon pairs built with one photon in the given cell. The achieved calibration level is estimated to be 3% and adds up quadratically to the constant term of the energy resolution.
The PHOS [21,24] is a homogeneous electromagnetic calorimeter composed of lead tungstate, PbWO 4 . The size of its elementary active units, also called cells, is η × φ = 0.004 × 0.004 (≈ 2.2 × 2.2 cm 2 ). Thus, the lateral dimensions of the cells are slightly larger than the PbWO 4 Molière radius of 2 cm. APDs with an active area of 5 × 5 mm 2 detect the scintillation light generated within the detector cells. The spectrometer covers φ = 60 • in azimuth and |η| < 0.12 in pseudorapidity and is located at a distance of 4.6 m from the interaction point. It is operated at a temperature of − 25 • C, at which the light yield of PbWO 4 increases by about a factor of three compared to room temperature. The energy resolution of the PHOS is The fine granularity of the detector enables the measurement of π 0 candidates up to p T ≈ 50 GeV/c.
The ITS [26] consists of three sub-detectors each with two layers to measure the trajectories of charged particles and to reconstruct primary vertices. The two innermost layers are the silicon pixel detectors (SPD) positioned at radial distances of 3.9 and 7.6 cm. The middle two layers are silicon drift detectors (SDD) located at 15.0 and 23.9 cm relative to the beam line. The outer two layers are silicon strip detectors (SSD) located at radial distances of 38 and 43 cm. The two layers of SPD cover pseudorapidity ranges of |η| < 2 and |η| < 1.4, respectively. The SDD and SSD cover |η| < 0.9 and |η| < 1.0, accordingly.
The TPC [25] is a large (90 m 3 ) cylindrical drift detector filled with a gas mixture of Ne-CO 2 (90-10%). It covers a pseudorapidity range of |η| < 0.9 over full azimuth, providing up to 159 reconstructed space points per track. A magnetic field of B = 0.5 T is generated by a large solenoidal magnet surrounding the central barrel detectors. Charged tracks originating from the primary vertex can be reconstructed down to p T ≈ 100 MeV/c and charged secondaries down to p T ≈ 50 MeV/c [22]. The TPC provides particle identification via the measurement of energy loss, dE/dx, with a resolution of ≈ 5% [25]. Beyond the outer radius of the TPC, the transition radiation detector (TRD) and the time-of-flight detector (TOF) provide additional particle identification information, as well as allowing for improved momentum resolution and added triggering capability. The detectors represent most of the material between the TPC and the EMCal and hence dominate the material budget in front of the EMCal. These detectors are missing in front of PHOS in order to provide a minimal radiation length to profit from the high resolution of the spectrometer.
The V0 detector is made up of two scintillator arrays (V0A and V0C) [28] covering 2.8 < η < 5.1 and −3.7 < η < −1.7. It is used to provide a minimum bias (MB) trigger [29] and reduce background events [22]. It is also involved in the definition of calorimeter triggers [30,31] and is used for luminosity determination as described in the next section.
In addition, the T0 detector [32] was used for luminosity determination. It consists of two arrays of Cherenkov counters, T0A and T0C, which respectively cover 4.61 < η < 4.92 and −3.28 < η < −2.97. The T0 furthermore provides a precise timing signal to other detectors with a resolution of better than 50 ps, used as starting signal for the TOF detector for example.

Datasets and event selection
During the data taking period of pp collisions at √ s = 8 TeV in 2012, the LHC operated at high beam intensities of approximately 2 × 10 14 protons per beam. Collisions at the ALICE interaction point were realized using a so-called "mainsatellite" bunch scheme, which involved proton collisions between the high intensity main bunches and low intensity satellite bunches. The interaction probability per bunchsatellite crossing was about 0.01, corresponding to an average instantaneous luminosity of about 5 × 10 30 cm −2 s −1 . Background events caused by beam-gas interactions or detector noise are rejected in the analysis using the V0A and V0C timing information [22]. Pileup events, with more than one pp collision per bunch crossing, are rejected based on SPD pileup identification algorithms looking for multiple primary vertices in a single event [22]. Additionally, the SPD is used to reject background events by comparing the number of SPD clusters to the multiplicity of SPD track candidates found in the respective collision. Only events with a z-vertex position of |z| < 10 cm in the global ALICE coordinate system are accepted for the analysis.
Two different types of triggers were used during data taking to select the events to be recorded: the minimum bias (MB) trigger and the calorimeter triggers, which are provided by the EMCal and the PHOS, to enhance statistics at high p T by selectively recording events with high energy deposits in the calorimeters. The MB trigger is a hardware Level-0 (L0) trigger [29]. It requires at least one hit in each V0A and V0C [28]. Both calorimeters also provide L0 triggers: EMC-L0 [30] and PHOS-L0 [31]. These L0 calorimeter triggers are required to be in coincidence with the MB trigger and select events with a deposited energy exceeding a nominal threshold in 4 × 4 adjacent cells, which is set to E EMC-L0 ≈ 2 GeV and E PHOS-L0 ≈ 4 GeV, respectively. A software Level-1 (L1) trigger is also deployed for the EMCal which inspects events preselected by the EMC-L0 trigger [33]. The trigger algorithm is similar to the EMC-L0, but combines information from different trigger region units to enhance the trigger efficiency and overcome hardware boundary effects [33]. Additionally, a larger trigger threshold of E EMC-L1 ≈ 8.4 GeV is set to further obtain statistics at higher transverse momenta.
In order to correctly normalize each trigger, the trigger rejection factors (R F) are determined by constructing the ratio of cluster energy spectra from MB and calorimeter triggered events as a function of the cluster energy, E, which are shown in Fig. 1. The ratios are expected to follow a constant for high cluster energies, the so-called plateau region, assuming the triggers only enhance the rate of clusters but do not affect their reconstruction efficiency. To reduce the statistical uncertainties, the R Fs are always determined with respect to the next lower threshold trigger. The cluster energy ratios have a steep turn-on near the respective trigger threshold energies. Since the EMC-L0 trigger becomes fully efficient only above its triggering threshold of E EMC-L0 ≈ 2 GeV, there is a change of slope visible in the turn-on region of the EMC-L1 trigger. The turn-on curve of the PHOS-L0 trigger also changes its slope due to a non-uniformity of the channels hardware gains. However, only the R F plateau regions are mainly relevant for analysis, as they are needed to correctly normalize the triggered data, which are found to be: R F EMC-L0 = 67.0 ± 1.1, R F PHOS-L0 = (12.4 ± 1.5) × 10 3 and R F EMC-L1 = (14.9 ± 0.3) × 10 3 . The last factor is obtained by multiplying the two given rejection factors of the two EMCal triggers, see Fig. 1, as the R F for EMC-L1 to MB trigger is of interest.
The luminosity determination is based on the cross-section of the MB trigger condition, σ MB AND , measured in a van der Meer (vdM) scan [34,35]. The stability of the measured cross section throughout the whole data taking period is assessed by comparing the V0-based luminosity measurement with an independent luminosity signal, issued by the T0 detector. As discussed in Ref. [35], this comparison results in an overall normalization uncertainty of 2.6%, which includes contributions from both the vdM-based measurement and its stability over time. The integrated luminosity of each triggered sample is calculated with the number of analyzed events, N events , the respective rejection factors, R F, and the MB cross section, σ MB AND = 55.80 ± 1.45 (stat+sys) mb [35], given by: for which R F = 1 holds for the MB trigger. As the good run lists for each detection method do not coincide, integrated luminosities are individually quoted for all cases in Table 1.

Neutral meson reconstruction
Both π 0 and η mesons are reconstructed via their two-photon decay channels with branching ratios of 98.823 ± 0.034% and 39.31 ± 0.20% [36] by means of an invariant mass analysis. The neutral mesons are reconstructed using the two electromagnetic calorimeters, EMCal and PHOS, a photon conversion method (PCM) and a hybrid method, PCM-EMCal, which combines one photon candidate from the PCM and one from the EMCal, resulting in four (three) different methods for the reconstruction of π 0 (η) mesons. The reconstruction of η mesons is not accessible by PHOS due to the limited detector acceptance and, compared to the π 0 , the wider opening angle of the decay photons. The hybrid PCM-EMCal method benefits from the high momentum resolution of the PCM, a high reconstruction efficiency and, crucially, the triggering capabilities of the EMCal. Moreover, an extended p T coverage is achieved compared to the standalone EMCal measurement, as there is no limitation due to cluster merging effects, discussed later in this section. Photons and electrons/positrons generate electromagnetic showers when they enter an electromagnetic calorimeter. They usually spread their energy over multiple adjacent calorimeter cells. In order to reconstruct the full energy of impinging particles, those adjacent cells need to be grouped into clusters, which is realized by a clusterization algorithm. In the first step, the algorithm looks for the cell that recorded the highest energy in the event, exceeding the seed energy, E seed . After the identification of such a seed cell, adjacent cells with recorded energies above a minimum energy, E min , are added to the cluster. For the EMCal, the clusterization algorithm adds cells to the cluster as long as their recorded Table 1 The analyzed luminosities considering the individual statistics for the different reconstruction methods and triggers. The EMCal related measurements use the same list of good runs as indicated by the combined column. The uncertainties denoted with "sys" reflect the sys-tematical uncertainty of R F determination, whereas "norm" represents the uncertainties entering from the cross section determination of the MB trigger [35] L int (nb −1 ) MeV. Large clusters due to overlapping photon showers in the PHOS are separated into individual clusters by an unfolding method based on the knowledge of the lateral shape of the electromagnetic shower [37].
Cell energies are calibrated for both calorimeters to provide best estimates for the cluster energies. After the cell-bycell energy calibration of the EMCal [23,27], an improved correction for the relative energy scale as well as for the residual geometrical misalignment of the EMCal between data and MC simulations is derived by making use of the good momentum resolution of the PCM photon in the hybrid PCM-EMCal method. Using this method, the π 0 mass is evaluated as a function of EMCal cluster energy, E cluster , for data and MC. Therefrom, a cluster energy correction is deduced for the simulation, for which the reconstructed π 0 masses are adjusted to the measured mass positions in data. For E cluster ≈ 1 GeV, the correction is of the order of 2% and rises up to 4% for higher energies. Thus, a precise energy calibration scheme for the relevant energy regions is available which is found to be consistent for the EMCal and hybrid PCM-EMCal methods for π 0 as well as η mesons at the same time, hence demonstrating the validity of the procedure. After applying this calibration in the analysis, the π 0 and η mass values in data and MC are obtained for each p T bin and their ratio is computed. Then, the ratios are plotted versus p T and fitted with a constant, giving access to the residual miscalibration of the meson mass values between data and MC. Such residual offsets of 0.005 ± 0.043 and 0.14 ± 0.13% are found for π 0 and η mesons for the EMCal analysis, whereas 0.002 ± 0.042 and 0.02 ± 0.14% are obtained for PCM-EMCal, illustrating the performance of the calibration procedure. For the PHOS, the energy deposition in each cell is calibrated by adjusting the π 0 peak position in the invariant mass spectra of photon pairs to the true mass of the π 0 meson. The accuracy of this calibration procedure is estimated to be better than 1%. It is evaluated from a comparison of the π 0 peak width in calibrated data and MC simulations by introducing random, normal-distributed decalibration parameters to the MC simulation.
Photon identification criteria are applied to the sample of reconstructed clusters in order to primarily select clusters generated by true photon candidates. For the photon reconstruction with PHOS, relatively loose identification cuts are applied because the shower overlap is negligible and the combinatorial background is found to be small in pp collisions. A minimum cluster energy, E cluster > 0.3 GeV, as well as a minimum number of cells forming a cluster, N cell ≥ 3, are required in order to reject electronic noise and minimum ionizing particles which deposit about 270 MeV in the PHOS. For the EMCal, a minimum energy cut of E cluster > 0.7 GeV is applied and the minimum number of cells grouped in a cluster is set to N cell ≥ 2. Furthermore, the selection criteria of |η| < 0.67 and 1.40 < ϕ < 3.15 rad are imposed for EMCal clusters. Pileup from multiple events, which may occur within a readout interval of the front-end electronics, is rejected by applying a cluster timing cut relative to the collision time of −25 < t cluster < 25 ns for the PHOS and − 35 < t cluster < 25 ns for the EMCal. Thus, photon candidates from different bunch crossings are removed with high efficiency of > 99%. For the EMCal, all clusters matched with a primary charged track are rejected. This track matching procedure, referred to as general track matching, uses a track p T -dependent matching in η and ϕ, beginning from | η| < 0.04 and | ϕ| < 0.09 for very low track momenta of p T < 0.5 GeV/c and going down to | η| < 0.01 and | ϕ| < 0.015 for highest track momenta, using the p T -dependent matching conditions | η| < 0.01 + ( p T + 4.07) −2.5 and | ϕ| < 0.015 + ( p T + 3.65) −2 . Applying these conditions, a primary track to cluster matching efficiency of more than 95% is obtained over the full p T range, rising above 98% for the analyzed EMCal triggered datasets for p T beyond 10 GeV/c. To further enhance the photon purity and to reject neutral hadrons, a cluster shape cut of 0.1 ≤ σ 2 long ≤ 0.7 is applied for EMCal clusters, where σ 2 long stands for the smaller eigenvalue of the dispersion matrix of the shower shape ellipse defined by the responding cells and their energy contributions to the cluster [9,39]. The lower threshold of σ 2 long is chosen to remove contamination caused by neutrons hitting the APDs of the readout electronics.
Photons convert into lepton pairs within the detector material of ALICE with a probability of about 8.5%. The reconstruction of such photon conversion candidates using PCM may be divided into three major steps: (i) tracking of charged particles and secondary vertex (V 0 ) finding [37]; (ii) particle identification and (iii) photon candidate reconstruction and subsequent selection. The V 0 s used in this analysis are obtained during data reconstruction using all available tracking information, recalculating the momenta of the daughter tracks under the assumption that both daughters are created with parallel momentum vectors at the V 0 . The tracks associated with secondary vertices are required to have a minimum momentum of p track T > 50 MeV/c and at least 60% of clusters from the maximum possible number of clusters, that a particle track can create in the TPC along its path, need to be found. In order to reduce the contamination from Dalitz decays, conversion candidates are only considered with a vertex at a radial distance of at least R > 5 cm. In addition, a line-cut is applied to restrict the geometrical η distribution of the V 0 s in order to remove photon candidates that would otherwise appear outside the angular dimensions of the detector. The condition R conv > |Z conv |S ZR − 7 cm is applied with S ZR = tan (2 arctan(exp(−η cut ))) and η cut = 0.9, where the coordinates R conv and Z conv are determined with respect to the nominal center of the detector. Additional constraints are imposed on R conv < 180 cm and |Z conv | < 240 cm to ensure that the reconstruction of secondary tracks is performed inside the TPC. Electrons and positrons from photon conversions are identified via their energy deposit, dE/dx, in the TPC. The difference of the measured dE/dx value from the hypothesis of the electron/positron energy loss is used for particle identification. The dE/dx of measured charged tracks is required to be within − 3 < nσ e < 5 of the expected energy loss, which is a p T -dependent observable defined by nσ e = (dE/dx− dE/dx e )/σ e with the average energy loss of the electron/positron, dE/dx e , and the Gaussian width of the fit to the measured dE/dx distribution, σ e . To further reduce charged pion contamination as the pion dE/dxband begins to merge with the electron/positron dE/dx-band above p 4 GeV/c, a cut based on the separation from the hypothesis of charged pion energy loss is applied in nσ π , analog to the previous definition. Tracks with energy losses closer to the pion line than |nσ π | < 1 are removed, which is done up to 3.5 GeV/c. The non-photon V 0 can-didate contamination is further suppressed by a triangular two-dimensional cut, | pair | < pair,max (1 − χ 2 red /χ 2 red,max ), with χ 2 red,max = 30 and pair,max = 0.1. This cut is based on the reduced χ 2 of the Kalman-filter [40] hypothesis for the e + e − pair and on the angle pair between the plane perpendicular to the magnetic field of the ALICE magnet and the e + e − pair plane. Furthermore, a cut on the cosine of the pointing angle of cos(θ PA ) > 0.85 is performed, where the pointing angle, θ PA , is the angle between the reconstructed photon momentum vector and the vector joining the collision vertex. The remaining K 0 S , and contamination is removed by selecting q T < q T,max 1 − α 2 /α 2 max on the Armenteros-Podolanski plot [41] with q T,max = 0.05 GeV/c and α max = 0.95. Additionally, the PCM measurement requires an out-of-bunch pileup correction which estimates the contamination of photon candidates from multiple events overlapping in the TPC. The correction is based on a study of the distance of closest approach (DCA) of the conversion photon candidates which is the smallest distance in beam direction, z, between the primary vertex and the momentum vector of the photon candidate. Photon candidates from different events generate a broad underlying Gaussian-like DCA distribution, which is fitted in order to estimate the outof-bunch pileup contribution. The correction is found to be p T -dependent and ranges from 42% at low p T ≈ 0.35 GeV/c to 10% at high p T ≈ 11 GeV/c.
The hybrid PCM-EMCal method practically uses the same cuts on photon candidates as the respective standalone reconstruction methods. In context of the PCM, a wider cut of − 4 < nσ e < 5 concerning the electron/positron energy loss hypothesis is used for the hybrid method and the p T restriction of the charged pion dE/dx cut is loosened. Only the upper value of the cut on the short axis of the moment of the shower shape for the EMCal is changed and required to be σ 2 long ≤ 0.5 in order to further reject contamination of hadrons [9]. Due to the timing constraint of the EMCal restricting clusters to triggered bunch crossings, no DCA or additional out-of-bunch pileup rejection is needed for the hybrid method. In addition to the general matching of primary charged particles to EMCal clusters already described, a dedicated track matching procedure for the two charged V 0 daughters with respect to EMCal clusters is applied. This cluster-V 0 track matching is the most important ingredient for the hybrid analysis, as pairing one leg of the V 0 candidate with the EMCal cluster generated by one of these secondary charged tracks itself, leads to an auto-correlation and causes a broad peak between the masses of the π 0 and η mesons at around 300 MeV/c. The same parameters from the general track matching procedure are found to remove about 99% of such candidates.
Invariant mass distributions of photon pairs, shown in Figs. 2 and 3, include combinatorial background as well as  which have been subtracted from raw real events to obtain the signal displayed with red data points. The blue curves represent fits to the background-subtracted invariant mass spectra. Additional examples of invariant mass distributions for the different methods are given in Ref. [38] the neutral meson signal for photon candidate pairs from the same, real event. An opening angle cut of 17 mrad for the angle between the momentum vectors of the two paired photon candidates is applied for the EMCal measurement. Requiring a minimum separation between such pairs is needed to ensure a proper background description by event mixing, in which two clusters from different events might otherwise be separated by an arbitrarily small distance. In same events, such cluster configurations would overlap partially or even merge into single clusters, which has been explicitly considered for event mixing by not allowing the cells with largest deposited energies of respective clusters to be direct neighbors on the EMCal surface. For the PCM and hybrid PCM-EMCal methods, an opening angle cut of 5 mrad is further applied between the momentum vectors of the pair of conversion photon candidates and accordingly,  the pair of PCM and EMCal photon candidates. Furthermore, pairs are restricted to a rapidity of |y| < 0.12 for the PHOS and |y| < 0.8 for all other methods. The uncorrelated combinatorial background is estimated by using an event mixing technique, in which photon candidates from different events are paired in order to prevent correlations between the candidates. Different event pools are used for this purpose, binned by primary vertex position, photon candidate multiplicity and transverse momentum to ensure the mixing of similar events only. In contrast to sameevent combinations to extract the neutral meson signal, the mixed-event background is obtained with up to 80 different events, stored in each of the event pools, in order to minimize its statistical uncertainties. Therefore, the mixed-event background needs to be scaled to match the integral of the raw signal in the vicinity of the right side of the neutral meson peak, just outside the peak integration interval, after which it is subtracted from the raw distribution. The backgroundsubtracted signal is then fitted to determine the mass peak position and width of π 0 and η mesons for every p T bin. A function composed of a Gaussian modified by an exponential tail at the low mass side [42] is used for this purpose. The low mass tail accounts for late conversions of one or both photons for the EMCal method and for energy loss effects due to bremsstrahlung for the PCM and hybrid PCM-EMCal methods. To reflect the residual correlated background components which remain after the subtraction of the mixedevent background, the fitting is performed by including an additional first order polynomial function (deduced from MC simulations), which is also shown in Figs. 2 and 3 and which is further being subtracted from the invariant mass distribution. In contrast, a slightly different approach for the background description is followed for the PHOS as its limited acceptance results in a more complicated shape of the combinatorial background around the signal peak, especially at low p T . As both correlated and combinatorial backgrounds are influenced in the same manner, the ratio of the raw signal and mixed-event distributions is constructed and fitted with first or second order polynomial function outside the peak region. Then, the mixed-event distribution is scaled with the obtained polynomial function and subtracted from the raw signal, which can be followed in Fig. 2. A Crystal Ball function [43] is used as the main fit function for the PHOS method which also reproduces the tail at the low mass region to take into account the late conversion of photons in front of the calorimeter. The signal distribution is then obtained by subtracting the scaled mixed-event background from the raw invariant mass distribution. The resulting backgroundsubtracted signal distributions as well as raw signals from real events, the normalized mixed-event and residual background distributions are shown in Figs. 2, 3 and in Ref. [38] for the π 0 and η meson mass region, respectively, for given example p T bins, illustrating the meson reconstruction over the full reported p T range.
The neutral meson raw yields are extracted by integrating the background-subtracted invariant mass distributions. The integration windows are defined by the reconstructed mass position and width obtained by the respective fits of the signal distribution in a given p T bin. For the PHOS method, the integration range for π 0 is asymmetrically defined as [−5σ, +3σ ] around reconstructed peak position, where σ is the standard deviation of the Gaussian part of the Crystal Ball function to take the asymmetric shape into account. For the other methods, the integration windows for both neutral mesons are chosen to cover at least [−3σ, +3σ ] around the reconstructed peak position, where σ is the standard deviation of the Gaussian part of the fit function. For each reconstruction method, the peak position and width used for the signal extraction are shown in Fig. 4 as a function of reconstructed p T .
Corrections for geometrical acceptance and reconstruction efficiencies are evaluated using MC simulations. PYTHIA8 [18] and PHOJET [44] event generators with minimum bias processes are used for this purpose. The correction factors for both MC productions are found to be consistent and, hence, are combined. To generate enough statistics for high meson momenta to be able to correct the raw yields obtained with triggered data, a PYTHIA8 simulation is used that is enriched with jets, generated in bins of hard scatterings, p T,hard . Particles generated by the event generators are propagated through the ALICE detector using GEANT3 [45] which realistically reproduces interactions between the particles and the detector material. In the simulation, the same reconstruction algorithms and analysis cuts are applied as for real data. In Fig. 4, the reconstructed π 0 and η mass peak positions and widths are compared as a function of p T between data and MC to confirm a proper detector response in the simulation. The normalized correction factors, , for each method, containing the specific detector acceptances as well as the full reconstruction efficiencies, are shown in Fig. 5. For the EMCal analysis, the correction factor for the π 0 is observed to decrease for p T 10 GeV/c. This is due to the effect of cluster merging, as due to the Lorentz boost the opening angles of π 0 mesons become too small to resolve adjacent clusters given the finite segmentation of the calorimeter. While the dominant symmetric decays are first to merge, the asymmetric decay contributions become more relevant at higher momenta. Above a certain limiting momentum, it is no longer possible to separate the two decay photons of the π 0 , creating merged clusters that significantly reduce the reconstruction efficiency in the EMCal as seen in Fig. 5. Thus, the natural upper limit for the π 0 reconstruction with the EMCal is of the order of p π 0 T ≈ 20 GeV/c. In contrast, the PCM-EMCal hybrid approach overcomes the limitations of the EMCal cell segmentation and makes it possible to reconstruct π 0 mesons up to p T ≈ 35 GeV/c as reported in this paper. For the PHOS, such cluster merging effects are negligible for the reported p T range owing to the high granularity of the calorimeter. Since the opening angles of photons from η meson decays are much larger compared to the π 0 , merging effects are negligible for all approaches over the full reported p T range in this case.
The contributions of secondary π 0 from weak decays and hadronic interactions with the detector material are estimated and removed for the π 0 measurements. Weak decays of K 0 S represent the main source of secondaries. For all reconstruction methods, the spectra of the three main particles relevant for the secondary correction due to weak decays, K 0 S , K 0 L and , are obtained from Refs. [46][47][48] with extrapolation of spectra to √ s = 8 TeV, assuming a power law for each  The normalized correction factors, , for each reconstruction method for π 0 (left) and η mesons (right) plotted versus p T bins used in the analysis. The factors contain the detector acceptances and the respective reconstruction efficiencies, where acceptances are further normalized by the rapidity windows accessible with each method, y, and full azimuth coverage of 2π , in order to enable a direct comparison between the different methods p T bin as function of √ s. These spectra are used as weights in a PYTHIA6.4 simulation, where the respective particle decays are simulated on generator level, taking into account the full decay kinematics. Using this procedure, the invariant yields of secondary π 0 s from weakly decaying particles are obtained. From the full ALICE GEANT3 MC simulations, the acceptance and reconstruction efficiencies are cal-culated for these secondaries and multiplied with the respective invariant yields from the generator level MC simulation to arrive at the secondary π 0 raw yields from the different particles. On the other hand, the π 0 raw yield from interactions with the detector material is purely obtained from the full MC simulation, which is the only viable approach. All the estimated, secondary π 0 raw yields are subtracted from the reconstructed number of π 0 s, as indicated in Eq. 2. The corrections are of the order of 1-3% for K 0 S , < 0.5% for K 0 L , 0.02% for and 0.1-2% for material, varying within the given values for the different methods and triggers used.
As there are three different triggers available for the EMCal and hybrid PCM-EMCal methods, and two different ones for the PHOS measurement, each with its own statistical and systematic uncertainties, as well as correlations between the different systematical uncertainties, the results from each trigger class are properly combined in order to obtain the final result for each reconstruction method. Statistical uncertainties are ensured to be uncorrelated since different triggers use non-overlapping data samples. For the systematic uncertainties, the p T -dependent correlation coefficients are determined. Only a few systematic uncertainties are found to be uncorrelated, such as the uncertainty of signal extraction and partly "efficiency" and "trigger" related uncertainties, for which further details are contained in Sect. 5. The correlation coefficients are found to be generally above 0.8. The respective p T -dependent weights are calculated according to the BLUE algorithm [49][50][51][52][53], which are used to combine the spectra from each method.

Systematic uncertainties
Systematic uncertainties are respectively summarized in Tables 2, 3 and 4 for the neutral mesons π 0 , η and their ratio η/π 0 . The values are given in percent and refer to relative systematic uncertainties of the measured values. Three different example p T bins are listed for each reconstruction method in order to illustrate their relative strengths. An additional, more detailed description of the systematic sources and the determination of uncertainties for all methods except PHOS may be found in Ref. [9], which is fully applicable to this paper.
For the π 0 measurement by PHOS, the systematic uncertainty related to the signal extraction is evaluated by varying the fitting range and the assumptions about the mass peak and background shapes. The systematic uncertainty related to the material budget is taken from Ref. [3], which is estimated by comparing the results of the analysis with and without magnetic field in the ALICE solenoid. Photons, which converted to e + e − pairs within the detector material, are most likely being reconstructed as two clusters in the presence of a magnetic field. Without a field, the secondary tracks from photon conversions are less separated and can be dominantly detected as single clusters, building the correct invariant masses for π 0 s in a di-cluster analysis. Therefore, comparing the π 0 spectra from data and MC with nominal and zero magnetic fields is a straightforward method to evaluate the uncertainty of the material budget description in simulations. Systematic uncertainties due to the cluster energy calibration are decomposed into the uncertainty of the energy scale of clusters and non-linearity effects. The energy scale Table 2 Summary of relative systematic uncertainties in percent for selected p T bins for the reconstruction of π 0 mesons. The statistical uncertainties are given in addition to the total systematic uncertainties for each bin. Moreover, the combined statistical and systematic uncertainties are also listed, obtained by applying the BLUE method [49][50][51][52][53] for all reconstruction methods available in the given p T bin, considering the uncertainty correlations for the different methods (see Sect. 6 for further details). The uncertainty from σ MBAND determination of 2.6%, see Ref. [35], is independent from the reported measurements and is separately indicated in the following plots below   Table 3 Summary of relative systematic uncertainties in percent for selected p T bins for the reconstruction of η mesons, see  Table 4 Summary of relative systematic uncertainties in percent for selected p T bins for the determination of the η/π 0 ratio. The statistical uncertainties are given in addition to the total systematic uncertainties for each bin. Moreover, the combined statistical and systematic uncertainties are listed as well, see also explanations in caption of uncertainty of 0.1% is estimated from a comparison of the π 0 mass peak position for the two-photon invariant mass spectra in data and MC. This energy uncertainty is translated to an uncertainty of the π 0 yield by convolution with the shape of the p T spectrum. The systematic uncertainty due to the non-linearity correction is evaluated by introducing different non-linearity correction schemes and calibration parameters for the MC simulation, whereas the p T dependence of the π 0 peak position and width is always kept consistent with data. The efficiency uncertainty consists of acceptance variations and differences between MC event generators. The acceptance uncertainty is estimated by changing the good cluster selection criteria, and the MC generator-dependent uncertainty is evaluated by comparing efficiencies of MB MC generators and single particle MC simulation which generates events containing single neutral mesons with realistic transverse momentum and rapidity distributions. Moreover, it includes the trigger efficiency uncertainty for the high energy photon trigger analysis, which is estimated by comparing the trigger turn-on curve from data with MC simulations. "Trigger normalization & pileup" summarizes systematic uncertainties due to the trigger normalization factor and pileup effects. The uncertainty related to the trigger normalization factor is estimated by changing the range of the fit to determine the rejection factor (R F). Furthermore, the out-of-bunch pileup contribution is evaluated by varying the timing cut to accept clusters. For the PCM measurement, the main source of systematic uncertainty is the material budget, for which the same value is used as previously calculated in Ref. [3]. The signal extraction uncertainty is estimated by changing the integration window around the invariant mass peak, the normalization range of the mixed-event background and by using different order polynomials as well as other fit functions to evaluate the remaining background contribution. "Track reconstruction" summarizes the systematic uncertainties found by requiring different numbers of TPC clusters and by applying different minimum transverse momentum cuts on tracks. The systematic uncertainties due to the electron identification ("electron PID" and "PCM photon PID") are determined by varying the PID cuts, which are elaborated in Sect. 4, and by comparing the respective results. For PCM, the "trigger normalization & pileup" uncertainty is dominated by the uncertainty of the DCAz background description for the out-of-bunch pileup estimation. Furthermore, it contains the systematic uncertainty due to the pileup rejection by the SPD due to its finite efficiency to remove pileup events.
For the EMCal, one main systematic uncertainty arises from the knowledge of the outer material budget, which is defined by all detector components from the radial center of the TPC up to the EMCal. The uncertainty is assessed by running the analysis only with/without TRD modules in front of the EMCal, since part of the data taking in 2012 occurred with the EMCal only partially obscured by the TRD. Since TRD and TOF have similar material budgets, the same uncertainty is assigned to the TOF as well, which covered the full polar angle so that a similar assessment as for the TRD is not feasible. Both uncertainties are quadratically combined to arrive at the given uncertainties which are listed in the tables. The signal extraction uncertainty contains the systematic uncertainties obtained from variations of the background normalization region, the choice of the background fit function and integration intervals, analog to the PCM method, as well as from variations of the minimum opening angle cut on the meson level. The systematic uncertainty denoted as "cluster description" reflects the mismatch of the description of the clusterization process between data and MC simulations, giving rise to modified reconstruction efficiencies, which includes the following cluster related quantities: minimum energy, shower shape, number of cells and timing as well as variations of the energy thresholds used for the clusterization process.
Moreover, cell timing cut variations are also included in this category. "Cluster energy calibration" considers the systematic uncertainties due to non-linearity effects and the energy scale of clusters. Different non-linearity schemes are used in this analysis from which this uncertainty is obtained. Moreover, the energy scale uncertainty is determined by obtaining the residual differences of reconstructed meson mass values from data and MC simulations. The systematic uncertainty induced by the charged particle veto on cluster level, introduced as "general track matching" in Sect. 4, is determined by variations of the matching residuals. The "efficiency" uncertainty reflects differences between MB MC generators for the calculation of reconstruction efficiencies. Moreover, it contains the uncertainty of the actual trigger turn-on, obtained by comparing the turn-on curves in data and MC. The uncertainties from the determination of trigger rejection factors (R F) as well as from the pileup rejection by the SPD, which has a finite efficiency for pileup removal, are summarized with "trigger normalization & pileup".
For the hybrid method PCM-EMCal, the same cut variations are performed as for the standalone methods. However, given the fact that only one photon candidate of each system is used, most systematic uncertainties are found to be of different size or behavior, e.g. the minimum opening angle cut variations. The "track matching to cluster" uncertainty reflects the V 0 -track to cluster matching, which is assessed by varying the matching residuals.
As indicated in Table 4, many uncertainties cancel for the η/π 0 ratio, such as the material-related systematics. For the remaining categories, the respective uncertainties of the π 0 and η measurements are added quadratically and canceled partially beforehand, if applicable.

Results
The invariant differential cross sections of π 0 and η production are obtained by where N π 0 (η) is the number of reconstructed π 0 (η) meson in a given p T bin, N π 0 sec represents the estimated number of secondary π 0 mesons, L int is the integrated luminosity, A · ε rec is the product of the geometrical acceptance and reconstruction efficiency, also referred to as in Fig. 5, Br π 0 (η)→γ γ is the branching ratio for the two-gamma decay channel and y p T is the bin width in rapidity and transverse momentum. For the measurement of π 0 mesons by PCM, the out-  Fig. 6 Ratios of the fully corrected π 0 (left) and η (right) spectra for each reconstruction method to the TCM fit of the combined spectrum of-bunch pileup correction has to be noted for completeness and to be applied as well.
The invariant differential cross sections are independently calculated for each method. The final spectra are obtained by combining the results in the overlap regions using the BLUE method [49][50][51][52][53], properly taking into account the correlations of the systematic uncertainties of the different methods. Possible statistical correlations between the measurements, for instance due to the conversions at small distances relative to the beam axis, are negligible due to the small conversion probability and the small likelihood of reconstructing the respective electron in the calorimeters leading to a meson candidate which finally ends up in the respective integration window. As there are no common uncertainties present for PCM, EMCal and PHOS, all systematic uncertainties are considered to be completely uncorrelated in those cases. On the other hand, the correlations introduced by including the hybrid PCM-EMCal measurement have to be taken into account. By construction, there are different numbers of conversion photons entering the two methods. Thus, all systematic uncertainty sources from PCM are found to be partially correlated in the PCM-EMCal method. Half of the size of the material budget uncertainty, for example, is assumed to be uncorrelated. Furthermore, the uncorrelated systematic uncertainties from PCM-EMCal with respect to PCM are, with full size, all the calorimeter related uncertainties as well as trigger and efficiency uncertainties.
Due to finite bin widths of the measured production cross sections, the neutral meson spectra are shifted along the horizontal axis [54]. All bin width corrections are of the order of 1% and below. In contrast, the reported η/π 0 ratios are shifted along the vertical axis, as otherwise the ratio could not be computed and the different measurements could not be combined. The correction is below 1% for p T > 2 GeV/c, but becomes significant for smaller momenta and rises to 8% for the lowest bin. The combined invariant cross sections of inclusive π 0 and η meson production cover transverse momentum ranges of 0.3 < p T < 35 and 0.5 < p T < 35 GeV/c, respectively. The total uncertainties of the measurements, obtained by quadratically adding the combined statistical and systematic uncertainties, are of the order of 5% for the π 0 and 10% for the η meson for most of the p T bins covered, increasing for lowest and highest momenta due to statistical limitations as well as systematic effects. Both combined neutral meson spectra are fitted with a two-component model (TCM), proposed in Ref. [55], by using the total uncertainties for each p T bin. The functional form of the TCM is a combination of a Boltzmann component and a power-law part, which, in general, should be the dominant components at low and high p T , respectively. The fit function is able to reproduce the spectra over the full p T range and is described as: where E T,kin = p 2 T + m 2 − m is the transverse kinematic energy with the meson rest mass m and A e , A, T e , T as well as n are the free parameters. To compare the different methods, the ratios of spectra measured by each reconstruction method to the TCM fit of the combined spectrum are shown in Fig. 6. The vertical error bars represent the statistical uncertainties, whereas the boxes quantify the bin widths in horizontal direction and the systematic uncertainties in vertical direction. All measurements agree within uncertainties over the full p T range.
The π 0 and η meson cross sections are also fitted with a Tsallis function [56], which has been used in previous measurements of π 0 and η meson production in pp collisions reported by ALICE [3,4]: Parameters of the fits to the π 0 and η invariant differential cross sections using the TCM fit [55]   where C, n and T are free parameters of the fit with m and m T being the rest as well as the transverse mass of the meson. The fit parameters extracted from both the TCM and Tsallis fits are summarized in Table 5. The TCM is chosen as the standard fit function, since it better describes the spectra at low and high p T than the Tsallis counterpart [38]. This is also reflected in the smaller values obtained for the reduced χ 2 red of the respective fits, which are also recorded in Table 5. These values are calculated without assuming any correlation of systematic uncertainties and are found to be rather small for both fits, as the total uncertainties of meson spectra are used for their calculation. A direct comparison of TCM and Tsallis fits can be found in Fig. 7, where both fits are plotted, in addition to the measured spectra and theory calculations.
The measured invariant differential cross sections are compared with NLO pQCD calculations [8,11] using MSTW08 (PDF) [16] with DSS14 (FF) [8] for π 0 and CTEQ6M5 (PDF) [17] with AESSS (FF) [11] for the η meson. The same factorization scale value, μ, (0.5 p T < μ < 2 p T ) is chosen for the factorization, renormalisation and fragmentation scales used in the NLO pQCD calculations. For the π 0 the NLO PDF, pQCD and FF combination describes the RHIC data rather well [57], whereas for √ s = 2.76 TeV pQCD overpredicts ALICE data by 30% at moderate p T and agrees at higher p T [9]. The ratios of data and NLO pQCD predictions to the TCM fits of neutral meson spectra are shown in Fig. 7. The largest uncertainty of the NLO pQCD calculation is due to the choice of μ. For all μ values, these calculations overestimate the measured data for both π 0 and η mesons. This is also observed for meson measurements at √ s = 2.76 TeV by ALICE [9], although better description of data is achieved for μ = 2 p T , for which calculations are above data by 10−40% depending on p T . It has to be noted that FF uncertainties of NLO pQCD calculations have been considerably reduced after including the published π 0 measurement at √ s = 7 TeV [3] for DSS14. Including precise new data for η meson production measured at √ s = 2.76, 7 and 8 TeV [3,9] will also help to considerably reduce NLO pQCD uncertainty bands in that case. In addition, the reported neutral meson measurements at √ s = 8 TeV are compared to PYTHIA8.210 [18] references; Tune 4C [19] and Monash 2013 tune [20]. To enable a proper comparison of the PYTHIA tunes with the measured neutral meson spectra, π 0 from decays of long-living strange particles (K 0 S , , and ) are excluded. The Tune 4C calculation is about 30% above the π 0 measurement for p T > 1.5 GeV/c. In contrast, the Monash 2013 tune reproduces the π 0 spectrum within 10% for almost the complete transverse momentum range, although both tunes are not able to describe the shape of the measured spectrum indicated by the bump at approximately 3 GeV/c. Concerning the η meson, both tunes reproduce the measured spectrum for p T > 1.5 GeV/c within uncertainties. At lower momenta p T < 1.5 GeV/c, both tunes follow the same trend and deviate significantly in magnitude and shape from data. The tuning parameters of the soft QCD part of PYTHIA apparently fail to describe the measured η meson spectrum below p T < 1.5 GeV/c, whereas there is further tension up to p T ≈ 3.5 GeV/c. On the other hand, both PYTHIA tunes are consistent within uncertainties with the measured π 0 spectrum for the low transverse momentum interval 0.3 < p T < 1.5 GeV/c.
The mean transverse momenta, p T , are determined for the neutral meson spectra shown in Fig. 7. Three different fit functions are used in this context: a TCM, Eq. 3, a Tsallis, Eq. 4, and a modified Hagedorn [58] fit that is used as the default fit function, since it yields the best agreement with data at lowest p T measured [38]. The obtained values for π 0 and η mesons are listed in Table 6, where statistical and systematic uncertainties are quoted. The additional uncertainty term denoted with "fit sys" reflects the choice of the fitting function. Moreover, the introduced fit functions are also used to calculate the integrated yields, dN /dy| y ≈ 0 , for both neutral mesons in inelastic events. The cross section for inelastic pp collisions at √ s = 8 TeV, σ INEL = 74.7 ± 1.7 mb [59], is used for this purpose. The obtained yields are given in Table  6, which are based on extrapolation fractions, F extpol , of about 45% for the π 0 and about 34% for the η meson. Additionally, the integrated η/π 0 ratio is estimated and can be found in Table 6 as well. For the recent paper by ALICE on neutral meson production in pp collisions at √ s = 2.76 TeV [9], the mean p T as well as the integrated yields are also calculated for the reported spectra, which are furthermore added to Table  6. The inelastic pp cross section at √ s = 2.76 TeV, quoted in Ref. [9] as well, is used to calculate the integrated yields which include extrapolation fractions of about 59% for the π 0 and about 52% for the η meson. The obtained values for p T and dN /dy| y ≈ 0 for both neutral mesons are compared with measurements of average transverse momenta of charged particles [60] and with results concerning charged-particle multiplicity [61]. Due to a large extrapolation fraction of the π 0 and η meson spectra with respect to charged particles and the given systematics for the lowest transverse momenta, the uncertainties of p T and dN /dy| y ≈ 0 are found to be larger. Hence, the integrated η/π 0 ratios are also affected. Nevertheless, all values quoted in this paragraph are consistent within experimental uncertainties with the results from charged particle measurements [47,62]. Within their substantial uncertainties, the η/π 0 ratios at both pp energies are found to be consistent as well.
Both meson spectra, which are shown in Fig. 7, exhibit a similar power-law behavior, Ed 3 σ/d p 3 ∝ p −n T , with n π 0 = 5.936 ± 0.012(stat) ± 0.023(sys) and n η = 5.930 ± 0.029(stat) ± 0.044(sys) for high momenta of p T > 3.5 GeV/c. This is also reflected in the η/π 0 ratio which is shown in Fig. 8. The ratio is flat for p T > 3.5 GeV/c with a constant value of C η/π 0 = 0.455 ± 0.006(stat) ± 0.014(sys). Despite of the inability of NLO calculations to describe individual π 0 and η meson spectra, the η/π 0 ratio is reproduced fairly well, as it can be followed from left part of Fig. 8. It has to be noted that a different FF for the π 0 is used to compile the theory curve, namely DSS07, since there is no recent calculation for the η meson available which could be compared to the recent DSS14 π 0 prediction. The agreement of pQCD calculations with the data can be viewed as an indication that the η/π 0 ratio is driven by the π 0 and η meson FFs in the fac- Table 6 The mean transverse momenta, p T , and integrated yields, dN /dy| y ≈ 0 , for ALICE measurements of π 0 and η mesons at √ s = 2.76 and 8 TeV are summarized [38]. It has to be noted that the uncertainties from the measurements of the inelastic cross sections are not included for the given numbers, which are +3.9% −6.4% (model) ± 2.0(lumi)% for √ s = 2.76 TeV [29] and ± 2.3% for 8 TeV [59]. Moreover, the integrated η/π 0 ratios are quoted for the different energies  Fig. 8 Left: η/π 0 ratio compared to NLO pQCD predictions using PDF CTEQ6M5 and FFs DSS07 for the π 0 and AESSS for the η, and PYTHIA8.210 calculations using Tune 4C and Monash 2013 tune. The total uncertainties of the measured η/π 0 ratio are of the order of 10% for most of the p T bins covered, increasing for lower and higher momenta due to limited statistics as well as systematic effects. Right: Comparison of the η/π 0 ratio to related, previous ALICE measurements as well as other experiments at lower collision energies, for which total uncertainties are drawn. Furthermore, a comparison to the η/π 0 ratio obtained with m T scaling is added torized picture of pQCD. A comparison of the reported η/π 0 ratio to the different PYTHIA tunes indicates an agreement within uncertainties down to p T ≈ 1.5 GeV/c, although the shape, as well as the ratio, cannot be fully reproduced below p T < 1.5 GeV/c due to the already mentioned deviations of PYTHIA tunes from data in this region.
The validity of m T scaling is tested by means of the η/π 0 ratio. For this purpose, the TCM parameterization of the measured π 0 spectrum, given in Table 5, is used to obtain the η spectrum via the application of m T scaling by replacing the π 0 mass with the η mass and using the normalization ratio C η /C π 0 = 0.455. From these two spectra, the η/π 0 ratio is constructed, plotted as blue curve in the right part of Fig. 8. The measured η/π 0 ratio is consistent with the m T scaling prediction (blue curve) above p T > 3.5 GeV/c. However, for smaller transverse momenta of p T < 3.5 GeV/c the ratio of the measured η/π 0 ratio over the η/π 0 ratio obtained with m T scaling, which can be found in Ref. [38], constantly decreases and reaches about 45% at around 1 GeV/c. For the region below 3.5 GeV/c, m T scaling is observed to be clearly broken with a significance of 6.2σ . Given this observation, the measured η/π 0 ratios in pp collisions at √ s = 2.76 TeV and 7 TeV, previously reported by ALICE [3,9], are re-evaluated. Whereas there is indication for a m T scaling violation with 2.1σ for 2.76 TeV, we also observe a significant disagreement of 5.7σ for 7 TeV. Hence, both ratios are found to be consistent with our observation at 8 TeV. Furthermore, both η/π 0 ratios are fitted with a constant for p T > 3.5 GeV/c, yielding values of C η/π 0 = 0.474 ± 0.015(stat) ± 0.024(sys) for 2.76 TeV and C η/π 0 = 0.476 ± 0.020(stat) ± 0.020(sys) for 7 TeV. They are consistent within uncertainties with the measured η/π 0 ratio at 8 TeV for the given p T range. Therefore, all three ALICE measurements are simultaneously fitted with a constant for p T > 3.5 GeV/c in order to obtain a combined value of C η/π 0 = 0.459 ± 0.006(stat) ± 0.011(sys). For the region p T < 3.5 GeV/c, all collision energies covered by ALICE also agree within experimental uncertainties, so that η/π 0 ratios may be claimed to be consistent within accuracy for ALICE measurements in pp collisions at √ s = 2.76, 7 and 8 TeV.
Before the LHC era, the precision of η/π 0 measurements was not sufficient to probe m T scaling over broad ranges of p T with high statistics. PHENIX and NA27 provide the η/π 0 ratio with highest accuracy at high and low p T and therefore are compared to the reported measurement. PHENIX measurements for pp collisions at √ s = 200 GeV are available only for p T region > 2.25 GeV/c [63], where π 0 and η spectra are already described by m T scaling. However, PHENIX notably does not apply any secondary π 0 correction concerning weak decays, which further has to be taken into account when comparing with data points from ALICE. Measurements of π 0 and η spectra in pp collisions at √ s = 27.5 GeV from NA27 [64] are used to obtain the η/π 0 ratio in the p T range of 0.4 < p T < 1.6 GeV/c. The paper does not mention a secondary correction of π 0 spectrum; however, it cannot significantly change the conclusions to be drawn from the measurement. Although the NA27 measurement provides the world's most precise published data points for the η/π 0 ratio at low p T < 2.0 GeV/c in the pre-LHC era for pp collisions, it is not conclusive concerning m T scaling violation. The first NA27 points at p T < 1 GeV/c are consistent with both the m T scaling curve and the new data from pp collisions at √ s = 2.76, 7 and 8 TeV within uncertainties, whereas for p T > 1 GeV/c the results of NA27 show a tendency to be above the m T scaling prediction, although uncertainties become significant. A clearer confirmation of the m T scaling at low p T is observed for other particle species, such as kaons, φ and J/ψ in pp collisions at √ s = 200 GeV [13]. Whether the magnitude of m T scaling violation depends on the collision energy can be clarified in future by ongoing analysis of hadron spectra measurements in pp collisions at √ s = 13 TeV delivered by the LHC.

Conclusion
The invariant differential cross sections of π 0 and η mesons in pp collisions at √ s = 8 TeV have been measured at mid-rapidity over a wide p T range by the ALICE experiment, using four different reconstruction methods for the π 0 , and three for the η meson. NLO pQCD calculations with MSTW08 (PDF) with DSS14 (FF) for the π 0 and CTEQ6M5 (PDF) with AESSS (FF) for the η meson overestimate the measured spectra of both neutral mesons. This discrepancy is also reported for pp collisions at √ s = 2.76 TeV by ALICE. However, the ratio of η/π 0 is reproduced by NLO pQCD calculations within uncertainties, which is a test for the FFs in the factorized picture of pQCD. The prediction from PYTHIA8.2 Tune 4C describes the η spectrum within uncertainties for p T > 1.5 GeV/c, but it is about 30% larger than the measured π 0 production cross section. On the other hand, the Monash 2013 tune agrees with the reported neutral meson measurements within 10% for p T > 1.5 GeV/c. Both PYTHIA8.2 tunes are able to reproduce the π 0 spectrum below p T < 1.5 GeV/c within uncertainties, but fail to describe the η spectrum in this region. The η/π 0 ratio is described by m T scaling for p T > 3.5 GeV/c, whereas a deviation from this empirical scaling law is found for p T < 3.5 GeV/c with a significance of 6.2σ . Within experimental uncertainties, the η/π 0 ratios measured by NA27, PHENIX and ALICE are in agreement for the covered transverse momentum intervals of each measurement, representing pp collisions starting at center of mass energies of √ s = 27.5 GeV up to √ s = 8 TeV.