Production of \pi^{0} and \eta mesons up to high transverse momentum in pp collisions at 2.76 TeV

ALICE Collaborationthanks: See Appendix A for the list of collaboration members

July 19, 2019

The invariant differential cross sections for inclusive and mesons at midrapidity were measured in pp collisions at  TeV for transverse momenta  GeV/ and  GeV/, respectively, using the ALICE detector. This large range in was achieved by combining various analysis techniques and different triggers involving the electromagnetic calorimeter (EMCal). In particular, a new single-cluster, shower-shape based method was developed for the identification of high- neutral pions, which exploits that the showers originating from their decay photons overlap in the EMCal. Above  GeV/, the measured cross sections are found to exhibit a similar power-law behavior with an exponent of about . Next-to-leading-order perturbative QCD calculations differ from the measured cross sections by about % for the , and between % for the meson, while generator-level simulations with PYTHIA 8.2 describe the data to better than %, except at  GeV/. The new data can therefore be used to further improve the theoretical description of and meson production.

1 Introduction

Measurements of identified hadron spectra in proton-proton (pp) collisions are well suited to constrain predictions from Quantum Chromodynamics (QCD) [Gross:1973ju]. Such predictions are typically calculated in the pertubative approximation of QCD (pQCD) based on the factorization of the elementary short-range scattering processes (such as quark–quark, quark–gluon and gluon–gluon scatterings) involving large momentum transfer () and long-range universal properties of QCD that need to be experimentally constrained. The universal properties are typically modeled by parton distribution functions (PDFs), which describe the kinematic distributions of quarks and gluons within the proton in the collinear approximation, and fragmentation functions (FFs), which describe the probability for a quark or gluon to fragment into hadrons of a certain type. The cross section for the production of a given hadron of type H can be written as a sum over parton types

(1)

where denotes the proton PDF of parton carrying a fraction of the proton’s longitudinal momentum, the FF of parton into hadron H carrying a fraction of the parton’s momentum, and the inclusive short-distance scattering cross section of partons and into  (see e.g. [Brock:1994er]).

Measurements of hadron production provide constraints on the PDFs and FFs, which are crucial for pQCD predictions, and at LHC energies probe rather low values of and . The neutral pion ( is of special interest because as the lightest hadron it is abundantly produced, and at LHC collision energies below a transverse momentum () of GeV/ dominantly originates from gluon fragmentation. While the collision energy () dependence of cross sections has been useful for guiding the parametrization of the FFs [deFlorian:2014xna], experimental data for neutral pions [Abelev:2012cn, Abelev:2014ypa] at the LHC are not available above  GeV/, where quark fragmentation starts to play a role. The new data presented in this paper extend our previous measurement [Abelev:2014ypa] in pp collisions at  TeV to values of  GeV/ allowing one to investigate the dependence of the cross section at high transverse momentum. In addition, we present the cross section of the meson, which due to its strange quark content provides access to the study of possible differences of fragmentation functions with and without strange quarks [Aidala:2010bn]. Furthermore, the meson constitutes the second most important source of decay photons and electrons after the . Hence, and meson spectra over a large range are needed for a precise characterization of the decay photon (electron) background for direct photon (semileptonic open charm and beauty) measurements.

The new measurement of the cross section is a result of five analyses using data from various ALICE detector systems and different identification techniques. The decay photons are either measured directly in the Electromagnetic Calorimeter (EMCal), the Photon Spectrometer (PHOS) or via the photon conversion method (PCM). In the PCM measurement, the photons are reconstructed via their conversions into pairs within the detector material, where the pairs are reconstructed with the charged-particle tracking systems. The is reconstructed statistically using the invariant mass technique. At high , where the decay photons are too close together to be resolved individually, the can still be measured via the characteristic shape of their energy deposition in the EMCal. We combine statistically independent analyses where (i) both photons are individually resolved in the EMCal (EMC), (ii) one photon is identified in the EMCal and one is reconstructed via its conversion to (PCM-EMC), and (iii) the photon pair’s energy is merged in the EMCal (mEMC). Finally, the previously published measurements based on methods where both photons are reconstructed with (iv) PHOS or (v) PCM are included as well [Abelev:2014ypa]. The addition of the EMCal based measurements extends the reach from  GeV/ to  GeV/, the highest for identified hadrons achieved so far. The meson cross section that was previously not available at TeV is measured in the range from to  GeV/ using the PCM, PCM-EMC and EMC methods. Consequently, the ratio is measured in the same range.

The article is organized as follows: Section 2 briefly describes the experimental setup. Section 3 describes the data samples and event selection. Section 4 describes the neutral meson reconstruction techniques and corresponding corrections for the cross section measurements. Section 5 discusses the systematic uncertainties of the various measurements. Section 6 presents the data and comparison with calculations and Section 7 provides a summary.

2 ALICE detector

A detailed description of the ALICE detector systems and their performance can be found in Refs. [Aamodt:2008zz, Abelev:2014ffa]. The new measurements primarily use the Electromagnetic Calorimeter (EMCal), the Inner Tracking System (ITS), and the Time Projection Chamber (TPC) at mid-rapidity, which are positioned within a  T solenoidal magnetic field. Two forward scintillator arrays (V0A and V0C) subtending a pseudorapidity () range of and , respectively, provided the minimum bias trigger, which will be further discussed in the next section.

The ITS [Aamodt:2008zz] consists of two layers of Silicon Pixel Detectors (SPD) positioned at a radial distance of 3.9  and 7.6 , two layers of Silicon Drift Detectors (SDD) at 15.0  and 23.9 , and two layers of Silicon Strip Detectors (SSD) at 38.0  and 43.0  from the beamline. The two SPD layers cover a pseudorapidity range of and , respectively. The SDD and the SSD subtend and , respectively. The primary vertex can be reconstructed with a precision of , where  (m, for the longitudinal () and transverse () directions, respectively,  m and .

The TPC [Alme:2010ke] is a large (90 m) cylindrical drift detector filled with a Ne/CO gas mixture. It covers a pseudorapidity range of over the full azimuthal angle for the maximum track length of 159 reconstructed space points. The ITS and the TPC were aligned with respect to each other to a precision better than using tracks from cosmic rays and proton-proton collisions [Aamodt:2010aa]. The combined information of the ITS and TPC allows one to determine the momenta of charged particles in the range of to GeV/ with a resolution between % at low and % at high . In addition, the TPC provides particle identification via the measurement of the specific energy loss (d/d) with a resolution of 5%. The tracking detectors are complemented by the Transition Radiation Detector (TRD) and a large time-of-flight (TOF) detector. These detectors were used to estimate the systematic uncertainty resulting from the non-perfect knowledge of the material in front of the EMCal.

The EMCal [Cortese:2008zza] is a layered lead-scintillator sampling calorimeter with wavelength shifting fibers for light collection. The overall EMCal covers in azimuth and in pseudorapidity. The detector consists of cells (also called towers) with a size of corresponding to about twice the effective Molière radius; the cells are read out individually. With a depth of  cm, or radiation lengths, cells comprise a physical module. The modules are arranged in full-sized and one-third-sized supermodules, consisting of and modules, respectively, of which only the full-sized modules, corresponding to an azimuthal coverage of , were readout for the data recorded in 111The detector was installed in its complete configuration by early 2012, while and full-sized supermodules were present in 2010 and 2011, respectively. The modules are installed with a radial distance to the nominal collision vertex of  m at the closest point, and assembled to be approximately projective in . The scintillation light from each cell is collected with wavelength shifting fibers that are connected to a  mm active-area avalanche photodiode. The relative energy and position resolutions improve with rising incident energy of the particle [Abeysekara:2010ze]. The energy resolution can be described by a constant and two energy dependent terms parametrized as % with , , and in GeV. The position resolution is linear as a function of and parametrized as with in GeV. Starting with the highest cell  GeV, the energy depositions from directly adjacent EMCal cells with  GeV are combined to form clusters representing the total energy and physical position of incident particles [Abelev:2014ffa]. The clustering algorithm allows only one local energy maximum in a cluster; if a second is found a new cluster is initiated. Each cell is restricted to only be part of one cluster. Individual cells were calibrated using the mass peak position evaluated cell-by-cell, achieving a relative variation of below %.

3 Data samples and event selection

The data presented in this paper were recorded during the 2011 and 2013 periods with pp collisions at . Various EMCal triggers were employed and, while the majority of the minimum bias data were recorded in 2011, the 2013 running period took advantage of higher threshold EMCal triggers to collect a notable high- data sample. For the pp data collected in 2011, the minimum bias trigger (MB) required a hit in either V0 detector or a hit in the SPD, while it required hits in both V0 detectors for the data collected in 2013 (MB). The respective cross sections were determined based on van-der-Meer scans, and found to be  mb with and  [Abelev:2012sea]. For the normalisation of the 2013 data, for which there was no vdM scan, the uncertainty was conservatively increased to %, to account for possible variations of the MB trigger efficiency between 2011 and 2013. The resulting uncertainty due to the luminosity determination is % for both datasets together.

The EMCal issues triggers at two different levels, Level 0 (L0) and Level 1 (L1). The events accepted at L0 are further processed at L1. The L0 decision, issued latest  s after the collision, is based on the analog charge sum of adjacent cells evaluated with a sliding window algorithm within each physical Trigger Region Unit (TRU) spanning cells in coincidence with a minimum bias trigger. The L1 trigger decision, which must be taken within  s after the collision, can incorporate additional information from different TRUs, as well as other triggers or detectors. The data presented in this paper used the photon (EG) trigger at L1, which extends the sliding window search across neighboring TRUs, resulting in a % larger trigger area than the L0 trigger.

Fig. 1: Energy dependence of ratios between cluster spectra for EMC1/INT1, EMC7/INT7, EG2/EMC7 and EG1/EG2. The trigger names INT1 and INT7 denote the minimum bias triggers MB and MB respectively. The trigger names EMC1, EMC7, EG2 and EG1 denote the EMCal triggers at L0 in 2011 and 2013, and the EMCal triggers at L1 in 2013 with increasing threshold respectively. The individual trigger rejection factors and their respective fit ranges in the plateau region are indicated as well. The final rejection factors with respect to the minimum bias trigger are given in Tab. 1.

In 2011, only the L0 trigger was used with one threshold (EMC1), while in 2013, one L0 (EMC7) and two L1 triggers (EG1, EG2) with different thresholds were used, as summarized in Tab. 1. The lower L1 trigger threshold in 2013 was set to approximately match the L0 threshold in 2011 for consistency. In case an event was associated with several triggers, the trigger with the lowest threshold was retained.

However, the thresholds are configured in the hardware via analog values, not actual units of energy. Their transformation into energy values directly depends on the energy calibration of the detector. For a reliable normalization of each trigger, the Trigger Rejection Factor () is used. The takes into account a combination of the efficiency, acceptance and the downscaling of the respective triggers. It can be obtained from the ratio  of the number of clusters reconstructed in EMCal triggered events to those in minimum bias events at high cluster energy where should be approximately constant (plateau region), assuming the trigger does not affect the cluster reconstruction efficiency, but only the overall rate of clusters. To reduce the statistical uncertainties on the normalization for the higher threshold triggers, was always estimated with respect to the trigger with the next lower threshold in the EMCal or the respective minimum bias trigger if no lower EMCal trigger was available. By consecutively multiplying the individual rejection factors up to the minimum bias trigger, the final was obtained with respect to the minimum bias trigger. The energy dependence of the ratios between cluster spectra of the relevant trigger combinations (EMC1/INT1, EMC7/INT7, EG2/EMC7 and EG1/EG2) are shown in Fig. 1. At low , there is a minimum at roughly the threshold of the lower-level trigger for EG2/EMC7 and EG1/EG2, while at high there is a pronounced plateau for every trigger combination. The averages above the threshold in the plateau region, which represent for the respective trigger combinations, are indicated by a line whose width represents the respective statistical uncertainty. The corresponding systematic uncertainties were obtained by varying the range for the fit of the plateau region. Finally, the values for the average trigger rejection factors above the threshold with respect to the corresponding minimum bias triggers are given in Tab. 1. For the PCM-EMC and EMC analyses, all available triggers were used, while for mEMC only the EMC1, EG2 and EG1 triggers were included. The collected integrated luminosities for minimum bias and EMCal triggers

(2)

where refers to for 2011 and for 2013, are summarized in Tab. 1. The statistical uncertainties on are treated as systematic uncertainties on the integrated luminosity.

Year Trigger Trigger Approx. Trigger rejection
name threshold factor ()
2011 MB INT1
EMCal L0 EMC1  GeV
2013 MB INT7
EMCal L0 EMC7  GeV
EMCal L1 (G2) EG2  GeV
EMCal L1 (G1) EG1  GeV
Table 1: Approximate trigger threshold and corresponding trigger rejection factor for EMCal triggers, as well as integrated luminosity for minimum bias and various EMCal triggers.

Monte Carlo (MC) samples were generated using PYTHIA8 [Sjostrand:2007gs] and PHOJET [Engel:1995sb]. The correction factors obtained independently from the two MC samples were found to be consistent, and hence combined. For mesons with  GeV/, as in the triggered or merged cluster analyses, PYTHIA6 [Sjostrand:2006za] simulations enriched with jets generated in bins of the hard scattering () were used. All MC simulations were obtained for a full ALICE detector description using the GEANT3 [Brun:1987ma] framework and reconstructed with the same algorithms as for the data processing.

The different triggers of the EMCal affect the properties of the reconstructible mesons, like the energy asymmetry () of the decay photons, and hence significantly alter the reconstruction efficiency above the trigger threshold in the trigger turn-on region. The efficiency biases induced by the triggers were simulated using the approximate thresholds and their spread for different TRUs. The bias was defined as the ratio of the or reconstruction efficiency in triggered events over that in minimum bias events. Figure 2 shows the dependence of for different triggers and reconstruction methods for the and meson. While is unity for the mEMC analysis in the considered kinematic range, it is significantly below one for the PCM-EMC and EMC neutral meson reconstruction, and reaches only at about twice the trigger threshold. The corresponding correction factors are found to be larger for the PCM-EMC compared to the EMC method, and larger for the than the meson. This is a consequence of the much lower energy threshold imposed on the photons reconstructed with PCM, which leads to wider opening angle and asymmetry distributions of the reconstructible mesons. At low , also exhibits the effect of the trigger on subleading particles, for which the efficiency in triggered events is strongly reduced. However, the various triggers are only used if the meson momentum is at least 1.5 times the trigger threshold, thus the effect on the subleading particles is neglible.

In the offline analysis, only events with a reconstructed vertex with with respect to the nominal interaction vertex position along the beam direction were used. The finite primary vertex reconstruction efficiency for the MB(MB) trigger of about () is taken into account in the normalization of the respective minimum bias triggers. Furthermore, only events with exactly one reconstructed vertex were accepted to remove pileup from in- and out-of-bunch collisions. While the in-bunch pileup is negligible after the vertex selection, the out-of-bunch pileup accumulating in the TPC due to its readout time of , needs to be subtracted statistically for the mesons measured with PCM, as described in LABEL:Abelev:2014ypa. For the  () mesons reconstructed with PCM the out-of-bunch pileup correction ranges from  () at low to about above  GeV/. Analyses involving the EMCal are not affected because contributions of clusters from different bunch crossings are suppressed by a suitable selection of clusters within a certain time window around the main bunch crossing.

Fig. 2: Efficiency bias induced by different triggers (EMC1, EMC7 and EG1) for neutral pions (left panel) and mesons (right panel) for PCM-EMC (open symbols) and EMC (closed symbols).

4 Neutral meson reconstruction

Neutral mesons decaying into two photons fulfill

(3)

where is the reconstructed mass of the meson, and are the measured energies of two photons, and is the opening angle between the photons measured in the laboratory frame. Photon candidates are measured either by a calorimeter or by PCM. Neutral meson candidates are then obtained by correlating photon candidates measured either by EMC, PHOS or PCM exclusively, or by a combination of them (PCM-EMC). The corresponding and meson measurements are described in Sect. 4.1. The typical opening angle decreases with increasing of the meson due to the larger Lorentz boost. For mesons with above GeV/, the decay photons become close enough so that their electromagnetic showers overlap in neighboring calorimeter cells of the EMCal. At above GeV/, the clustering algorithm can no longer efficiently distinguish the individual showers in the EMCal, and mesons can be measured by inspecting the shower shape of single clusters, referred to as “merged” clusters and explained in Sect. 4.2.

To be able to directly compare the reconstruction performances of the various measurement techniques and triggers, the invariant differential neutral meson cross sections were expressed as

(4)

with the inverse of the normalized efficiency

(5)

and integrated luminosity (see Eq. 2). The measured cross sections were obtaind by correcting the reconstructed meson yield  for reconstruction efficiency , purity  and acceptance , efficiency bias , integrated luminosity , as well as for the and interval ranges, and , respectively, and the decay branching ratio . For invariant mass methods, the effect of reconstructed photon impurities on the meson purity are significantly reduced due to the subtraction of the combinatorial background, and hence the resulting meson impurities were neglected. For the mEMC method, the purity correction was obtained from MC simulations tuned to data. In the case of neutral pions, the contribution from secondary was subtracted from before applying the corrections. The contribution from weak decays was estimated for the different methods by simulating the decays of the K and using their measured spectra [Abelev:2014laa], taking into account the reconstruction efficiencies, as well as resolution and acceptance effects for the respective daughter particles The contribution from neutral pions produced by hadronic interactions in the detector material was estimated based on the full detector simulations using GEANT3. Finally, the results were not reported at the center of the intervals used for the measurements, but following the prescription in LABEL:Lafferty:1994cj at slightly lower values, in order to take into account the effect of the finite bin width . The correction was found to be less than % in every interval for the , and between % for the meson.

4.1 Invariant mass analyses

Applying Eq. 3, the invariant mass distribution is obtained by correlating all pairs of photon candidates per event. The neutral meson yield is then statistically extracted using the distinct mass line shape for identification of the signal and a model of the background. In the following, only the new measurements are described. Details of the PCM and PHOS measurements can be found in Refs. [Abelev:2012cn, Abelev:2014ypa].

Track selection
Track quality selection GeV/
Electron selection
Pion rejection for GeV/,
for GeV/ (PCM)
for GeV/ (PCM-EMC)
Photon criteria
Conversion point
 cm  cm
 cm
Photon quality ,
with and
Armenteros-Podolanski ,
with GeV/ and
Table 2: Criteria for photon candidate selection for PCM.

For the reconstruction of photons with PCM, only tracks from secondary vertices without kinks with a minimum momentum of  GeV/ were taken into account. The tracks had to be reconstructed within the fiducial acceptance of the TPC and ITS and with at least % of the reconstructible track points in the TPC. The photon momentum resolution is better than % at low , resulting from the precise determination of the track momenta by the TPC. Furthermore, the associated energy loss measured in the TPC was required to be within of the electron expectation, where with and the average energy loss and resolution for particle , respectively. The contamination from charged pions was suppressed by excluding all track candidates within of the pion expectation. The charged pion rejection was applied for track momenta between  GeV/ for PCM and  GeV/ for PCM-EMC, while for PCM it was released to above  GeV/. Only conversions which were pointing to the primary vertex and could be reconstructed with a conversion point with  cm within the acceptance of the ITS and TPC were considered. Compared to previous PCM standalone measurements [Abelev:2014ypa], the photon candidate selection criteria were optimized in order to reduce the combinatorial background. In particular, a two dimensional selection on the reduced of the photon conversion fit and the angle between the plane defined by the conversion pair and the magnetic field was introduced to suppress random pairs. Furthermore, the selection in the Armenteros-Podolanski variables [podolanski:1954] was tightened to reduce the contamination from K and decays. A summary of the conversion photon selection criteria is given in Tab. 2.

Clusters in the EMCal were reconstructed by aggregating cells with  GeV to a leading cell energy with at least  GeV, and were required to have only one local maximum. Photon candidates were obtained from reconstructed clusters by requiring a cluster energy of  GeV to ensure acceptable timing and energy resolution and to remove contamination from minimum-ionizing ( MeV) and low-energy hadrons. Furthermore, a cluster had to contain at least two cells to ensure a minimum cluster size and to remove single cell electronic noise fluctuations. Clusters which could be matched to a track propagated to the average shower depth in the EMCal (at  cm) within and criteria that depend on track as given in Tab. 3, were rejected to further reduce contamination by charged particles. The track-to-cluster matching efficiency amounts to about % for primary charged hadrons at cluster energies of GeV, decreasing slowly to % for clusters of GeV. The removal of matched tracks is particularly important for the PCM-EMC method as otherwise a severe auto-correlation between the clusters originating from one of the conversion electrons and the conversion photon would be introduced. Such auto-correlated pairs strongly distort the shape of the invariant mass distribution between the and mass peak region. The standard track matching applied to each conversion leg allowed for the removal of these auto-correlation pairs with an efficiency of more than since the corresponding track was already found. An additional distinction between clusters from mainly photons, electrons and neutrons is based on their shower shape. The shower shape can be characterized by the larger eigenvalue squared of the cluster’s energy decomposition in the EMCal plane. It is expressed as

(6)

where and are weighted over all cells associated with the cluster in the or direction. The weights logarithmically depend on the ratio of the energy of a given cell to the cluster energy, as , and  [AWES1992130]. Nuclear interactions, in particular for neutrons, create an abnormal signal when hitting the corresponding avalanche photodiodes for the readout of the scintillation light. Such a signal is mainly localized in one high-energy cell with a few surrounding low-energy cells, and can be removed by requiring . While the showers from electrons and photons tend to be similar, they can be distinguished based on their elongation, as most of the low- electrons will hit the EMCal surface at an angle due to the bending in the magnetic field. Most of the pure photons are reconstructed with a ; only late conversions elongate the showers beyond this. Thus, rejecting clusters with   for EMC (PCM-EMC) rejects the contamination from late conversion electrons significantly. At very high transverse momenta ( GeV/), it also rejects part of the contamination from neutral pions for which both photons have been reconstructed in a single cluster. Contributions of clusters from different bunch crossings were suppressed by a suitable selection of clusters within a certain time window around the main bunch crossing. A summary of the selection criteria for EMCal photon candidates is given in Tab. 3.

Cluster reconstruction
Minimum cell energy GeV
Minimum leading cell energy GeV
Cluster selection
Selection in ,
Minimum cluster energy GeV
Minimum number of cells
Cluster-shape parameter (PCM-EMC)
(EMC)
(mEMC)
Cluster time  ns (2011)
 ns  ns (2013)
Cluster–track matching
Table 3: Criteria for photon candidate selection for EMCal-based methods.

The good momentum resolution for the PCM photon was exploited to derive an improved correction for the relative energy scale, as well as for the residual misalignment of the EMCal between data and simulation. The neutral pion mass was evaluated for the PCM-EMC method as a function of the EMCal photon energy for data and simulation. A correction for the cluster energy was deduced which for a given simulation adjusts the neutral pion mass peak position to the measured position in the data as a function of the cluster energy. Above  GeV, the corrections for the various MC datasets are typically about 3%.

Fig. 3: Invariant mass distributions in the peak region for INT1 (left panels) and EG1 (right panels) triggers and EMC (top panels) and PCM-EMC (bottom panels) methods.
Fig. 4: Invariant mass distributions in the peak region for INT1 (left panels) and EG1 (right panels) triggers and EMC (top panels) and PCM-EMC (bottom panels) methods.

Example invariant mass distributions obtained by correlating photons reconstructed with EMCal or by one photon from PCM and one from EMCal are shown in Fig. 3 for neutral pions and Fig. 4 for mesons. The combinatorial background was calculated using the mixed event technique [Kopylov:1974th] using event pools binned by primary vertex position, multiplicity and transverse momentum. The mixed-event background has been normalized to the right side of the peak. Additionally, a residual correlated background estimated using a linear fit was subtracted. Only pairs with a minimum opening angle of  () mrad for EMC (PCM and PCM-EMC) methods were considered for signal and background construction. Finally, pairs are restricted to rapidity of .

Fig. 5: Neutral pion (left panels) and meson (right panels) mass position (bottom panels) and width (top panels) for the PCM, PCM-EMC and EMC methods. The performance of PHOS for is taken from LABEL:Abelev:2014ypa. Data are displayed as closed symbols, simulations as open symbols.

A Gaussian with an exponential tail on the left side was fitted to the subtracted invariant mass distributions, in order to determine the mass position and width of the peak. The results of the fits for the mass position and widths of neutral pions and mesons are shown in Fig. 5. The performance of PHOS from LABEL:Abelev:2014ypa in the case of is added for completeness. For all systems, the data for both and are reproduced by the MC simulations to a precision on average better than 0.3% for the mass position. For EMC, the -dependence of the mass position is especially pronounced, due to non-linearity effects for low clusters, shower merging and shower overlaps, and decay asymmetry enhanced by the employed triggers at high . The widths of the meson peaks are similarly well described, with the expected ordering for the various methods. In particular, the peak widths of the PCM-EMC fits are between the standalone measurements of PCM and EMC and are comparable to the PHOS measurement above 7 GeV/. This illustrates that the inclusion of one photon from PCM significantly improves the resolution of the neutral meson measurements.

The neutral meson raw yield was extracted by integrating the background-subtracted invariant mass distributions around the measured peak mass. The integration windows for the different reconstruction techniques were adjusted based on the average width of the meson peaks and their signal shape: (, ), (, ) for PCM, (, ), (, ) for PCM-EMC, and (, ), (, ) for EMC. For both mesons, an asymmetric range around the measured mass position was used to account for the low mass tail originating not only from the bremsstrahlung energy loss of conversion electrons and positrons, but also from additional missing energy in the EMCal due to the partial reconstruction of the photon.

Fig. 6: Normalized efficiency for different methods of neutral pion (left panel) and meson (right panel) reconstruction methods. The values for PHOS are taken from [Abelev:2014ypa].

The corrections for the geometric acceptance and reconstruction efficiency for the different mesons were calculated using MC simulations as mentioned in Sect. 3. The acceptance for the EMCal reconstruction techniques was calculated as the fraction of  (), whose decay photons point to the EMCal surface (), compared to the  () generated with . In the case of PCM-EMC, only one photon was required to point to the EMCal surface, while the other was required to be within the acceptance of the TPC (). The output from the full event MC simulations was reconstructed and analyzed in the same way as the data. The reconstruction efficiency was calculated as the fraction of reconstructed mesons compared to the mesons whose decay photons passed the acceptance criteria. The normalized efficiency  (see Eq. 5) as a function of meson is shown in Fig. 6 for the various methods. For EMC, rises at low and reaches its maximum at about at  GeV/. Subsequently, drops due to the merging of the two clusters, and is already a factor of smaller at about  GeV/. In the case of the , the efficiency at GeV/ is not yet affected by the cluster merging due to its higher mass. The efficiency for PCM-EMC is approximately a factor smaller than for EMC for both mesons due to the conversion probability of about in the respective pseudorapidity window. For the , it is similar to that of PHOS. The small decrease at higher for the PCM-EMC results from shower overlaps of the EMC photon with one of the conversion legs, and thus a stronger rejection of the EMCal photons due to track matching. Relative to PCM-EMC, for PCM is suppressed by the conversion probability affecting both decay photons.

The correction for secondaries from hadronic interactions depends on for the EMC-related methods. It ranges from % at the lowest to % (%) above  GeV/ for the PCM-EMC (EMC) method. For PCM, the correction amounts to less than % independent of . However, the contribution of the neutral pions from K is strongly dependent due to the tight selection criteria forcing the photons to point to the primary vertex. The correction drops quickly from about % to less than % at  GeV/. For the PCM-EMC and EMC, the corresponding correction amounts to % and %, respectively, independent of in the measured range. Contributions from other weak decays are below 0.1% and thus neglected for all reconstruction techniques.

4.2 Single cluster analysis

At high the showers induced by the two decay photons from a neutral pion merge into a single EMCal cluster, and therefore are unidentifiable in an invariant mass analysis. Hence, for above 15 GeV/ we use a different approach, namely to reconstruct and identify based only on single clusters, exploiting that clusters at high mostly originate from merged decay photons.

Merged clusters from decays tend to be more elongated than clusters from photons and electrons, and their deformation is reflected by the shower shape , defined in Eq. 6. The shower shape distributions are shown for data and MC in Fig. 7 for candidates, i.e. clusters fulfilling the selection criteria listed in Tab. 3 except . The distribution is found to be fairly well described by the MC, in particular for . For , the dominant contribution to candidates is from merged showers, while for clusters dominate where only the energy of one decay photon contributed. The most significant background is from decay photons of the meson and direct photons, located mainly at . Hence, for the mEMC measurement, candidates are simply required to have in order to discriminate from decay and direct photons. Only candidates with a rapidity of are considered.

The corrections for the geometric acceptance, reconstruction efficiency, and purity were calculated using MC simulations as described in Sect. 3. The resulting efficiency is shown in Fig. 6 compared to the other neutral pion reconstruction techniques. At high , mEMC clearly has an advantage due to its larger coverage compared to PHOS, and the exploitation of merging of the decay photons in the EMCal.

Fig. 7: Shower shape () distributions for candidates with  GeV/ compared in data and MC (left panel), and corresponding signal and background contributions in MC (right panel).

The reconstruction efficiency was calculated by comparing the reconstructed with generator-level distributions within a rapidity of . By comparing measured and generated of the neutral pion, the resolution correction is included in the inefficiency correction. The resolution is significantly different for candidate clusters containing all or only parts of the decay products, i.e. single photons or conversions. If all decay products contribute to the cluster, the mean momentum difference between reconstructed and generated is smaller than with an RMS of % above  GeV/. Otherwise, the mean momentum difference can reach up to % depending on the fraction of decay particles which could be reconstructed and whether they converted in the detector material.

The purity represents the fraction of reconstructed clusters that pass all the selections and are from a decay. For GeV/, it is almost constant at around % with variations of %. As can be seen in Fig. 7, the largest contamination in the considered window originates from the meson decay ( after fine-tuning the ratio to the measured value), closely followed by the hadronic background consisting mainly of charged pions () and K (). The contamination from mesons rises by about % towards higher momenta, while the contamination from the other two sources decrease by about %. Fragmentation photons contribute to the background about %. Their contribution was additionally scaled up by up to a factor , given by the ratio of fragmentation photons to direct photons according to NLO pQCD calculations [Gordon:1993qc, Vogelsang:1997cq], to account for direct photons which are not included in generator. Lastly, prompt electrons contribute to the contamination about %.

The correction for secondary pions from K decays amounts to approximately %, as their reconstruction efficiency is very similar to that of primary , albeit with worse resolution. In addition, corrections for from weak decays from K and  (together only about 0.3%) and from secondary hadronic interactions (%) were applied.

5 Systematic uncertainties

The sources of systematic uncertainties associated with the various measurement techniques and their magnitude in different ranges, chosen to reflect the strengths of the various methods, are given in Tab. 4 for the meson, in Tab. 5 for meson and in Tab. 6 for the ratio. Since the measurements obtained with PCM-EMC, EMC and mEMC are a combination of multiple triggers, the systematic uncertainties associated with each method reflect the contribution of different triggered data samples weighted by their statistical uncertainties. The uncertainties for the were evaluated directly on the ratio in order to cancel correlated uncertainties between the and measurements. In the following, we first describe the uncertainties on photon candidates reconstructed with EMC and PCM, then those on the meson level, and finally those related to the overall normalization, in the same order as given in the tables.

EMCal clustering: The uncertainty on clustering quantifies the mismatch in the description of the clusterization process between data and simulation. It incorporates the uncertainties arising from the variation of the minimum energy and time on cluster and cell level, the minimum number of cells per cluster as well as the variation of the selection on the clusters. For mEMC, varying the selection on is especially important since it quantifies the uncertainty of how well the distributions of the background are described in the simulation, and was varied from to and . The corresponding uncertainties range between % and % depending on and method.

EMCal cluster energy calibration: To estimate the uncertainty of the cluster energy calibration, the remaining relative difference between data and simulation in the mass position of the neutral pion was used. On average, the difference is %, which leads to an uncertainty on the spectra of about % taking into account that they approximately fall with . In addition, the correction of the simulations for relative energy scale and residual misalignment, described in Section 4.1, was varied by changing the underlying parametrization of the mass position correction with . We chose only correction factors where the measured neutral pion mass position could be reproduced by the simulation to better than % over all . The overall resulting uncertainties range between % and % depending on and method. For the meson ( ratio), the uncertainties are approximately a factor  () larger at similar due to lower photon energies entering at the same meson .

Track matching to cluster: The uncertainty introduced by the imperfection of the cluster-track matching procedure was studied by repeating the measurements with different track-matching parameters. The criteria were varied from tight selections, which removed only centrally matched clusters, to rather loose selections allowing a distance of cells depending on and . At low the uncertainties on the measurement are below %, while with increasing higher track densities due to the jettier environment become more important and lead to uncertainties of about %. In the case of the , the uncertainties are generally larger, between % and %, due to the worse signal-to-background ratio. For the ratio, the uncertainty of the alone is used, since part of the uncertainty is expected to cancel.

interval (GeV/) 1.4–1.6 3.0–3.5 16–20 30–35
Method PCM P-E EMC PCM P-E EMC P-E EMC mEMC mEMC
EMCal clustering - 2.4% 4.9% - 2.1% 2.3% 6.2% 4.4% 4.6% 5.9%
EMCal energy calib. - 2.0% 4.9% - 2.1% 2.5% 5.4% 5.5% 4.2% 4.8%
Track matching - 0.9% 1.8% - 1.4% 1.7% 6.9% 6.7% 5.4% 6.1%
Secondary track reco. 1.6% 1.1% - 0.9% 0.8% - 5.7% - - -
Electron PID 1.3% 0.7% - 1.5% 0.6% - 12.7% - - -
PCM photon PID 1.7% 1.4% - 2.3% 1.1% - 13.4% - - -
Signal extraction 1.9% 1.5% 2.4% 4.0% 1.9% 1.5% 3.4% 14.1% - -
Efficiency - 2.0% 2.0% - 3.6% 2.5% 2.1% 2.1% 8.4% 7.1%
Secondary correction - - - - - - - - 1.8% 1.8%
Inner material 9.0% 4.5% - 9.0% 4.5% - 4.5% - - -
Outer material - 4.2% 4.2% - 4.2% 4.2% 4.2% 4.2% 4.2% 4.2%
Trigger norm.+pileup 0.8% - - 0.4% 1.1% 0.5% 7.5% 5.5% 8.0% 8.8%
Tot. sys. uncertainty 9.6% 7.6% 8.9% 10.3% 8.3% 6.5% 24.5% 18.6% 14.9% 15.6%
Stat. uncertainty 2.8% 2.0% 6.5% 5.1% 3.3% 2.8% 14.8% 15.6% 5.7% 11.3%
Table 4: Systematic uncertainty for various sources and methods assigned to the measurement at different intervals. For comparison, the total systematic and the statistical uncertainties are also given. P-E stands for PCM-EMC.

Secondary track reconstruction: The uncertainty on the secondary track reconstruction quantifies the uncertainty related to secondary track finding used in PCM. It is estimated by variation of the TPC found-over-findable cluster selection and the minimum cut as well as reducing the acceptance for the conversion photons in  requiring them to approximately point towards the EMCal direction. The uncertainty depends on the precision of the relative alignment and track matching efficiency between TPC and ITS in different sectors of the TPC, and hence can vary for different data taking periods and trigger conditions. For the EMCal triggers, for instance, the conversion photons are mainly sampled in the region directly in front of the EMCal, where the ITS had larger inefficiencies than in other areas. The uncertainties range from % to %.

Electron PID: Systematic uncertainty on the electron identification for the PCM photon reconstruction was estimated by varying the TPC -based electron inclusion as well as the pion rejection selections. The corresponding uncertainties are small at low  (%), where there is good separation between electrons and pions, but reach up to % at high , where electrons and pions can not be efficiently separated any longer.

PCM photon PID: The uncertainty assigned to the PCM photon reconstruction combines the contributions from varying the criteria for the photon quality and Armenteros-Podolanski selections. The uncertainties are slightly larger than those on the electron PID, with similar dependence, since both the electron and the photon PID selections attempt to reduce the contamination which increases with increasing . For the  ratio, it is one of the dominant uncertainties, in particular at high , as only a small fraction cancels in the ratio due to the different decay kinematics of the two mesons.

interval (GeV/) 1–1.5 3–4 10–12
Method PCM PCM-EMC PCM PCM-EMC EMC PCM-EMC EMC
EMCal clustering - 3.1% - 3.1% 2.7% 3.6% 3.1%
EMCal energy calib. - 3.0% - 3.2% 4.5% 5.0% 6.8%
Track matching - 8.9% - 4.9% 5.7% 6.6% 8.8%
Secondary track reco. 3.7% 3.3% 1.6% 3.3% - 4.1% -
Electron PID 2.1% 2.5% 2.4% 2.2% - 5.2% -
PCM photon PID 3.9% 7.7% 3.9% 7.3% - 11.2% -
Signal extraction 6.0% 16.4% 6.0% 8.1% 9.3% 11.8% 3.5%
Efficiency - 5.0% - 5.0% 5.7% 5.8% 5.3%
Inner material 9.0% 4.5% 9.0% 4.5% - 4.5% -
Outer material - 4.2% - 4.2% 4.2% 4.2% 4.2%
Trigger norm.+pileup 1.8% - 1.9% - 2.8% 7.0% 7.2%
Tot. sys. uncertainty 12.3% 22.5% 11.9% 15.5% 14.3% 22.6% 15.5%
Stat. uncertainty 20.4% 43.4% 17.2% 16.7% 10.8% 21.3% 8.9%
Table 5: Systematic uncertainty for various sources and methods assigned to the measurement at different intervals. For comparison, the total systematic and the statistical uncertainties are also given.
interval (GeV/) 1–1.5 3–4 10–12
Method PCM PCM-EMC PCM PCM-EMC EMC PCM-EMC EMC
EMCal clustering - 4.1% - 4.2% 2.4% 6.0% 2.8%
EMCal energy calib. - 4.1% - 4.3% 4.6% 6.6% 7.6%
Track matching - 8.9% - 4.9% 5.7% 6.6% 9.0%
Secondary track reco. 3.7% 4.5% 1.6% 4.2% - 8.1% -
Electron PID 2.1% 3.3% 2.4% 3.2% - 7.0% -
PCM photon PID 3.9% 7.7% 4.0% 6.5% - 12.7% -
Signal extraction 6.1% 16.6% 7.0% 9.1% 9.3% 10.5% 8.5%
Efficiency - 5.4% - 5.4% 3.8% 7.0% 4.3%
Tot. sys. uncertainty 8.4% 22.5% 8.5% 15.6% 12.6% 23.8% 15.4%
Stat. uncertainty 20.4% 44.1% 17.7% 17.9% 10.9% 22.1% 8.8%
Table 6: Systematic uncertainty for various sources and methods assigned to the measurement at different intervals. For comparison, the total systematic and the statistical uncertainties are also given.

Signal extraction: The uncertainties arising from the signal extraction for the invariant mass analyses were estimated by varying the integration window, the background normalization region as well as the minimum opening angle, and requiring a mild asymmetry of the decay photons. For the neutral pion, the signal extraction uncertainty for PCM ranges from % at low to % at higher , due to the good momentum resolution of the tracks. For PCM-EMC, the equivalent uncertainty ranges from % to % at low and high , respectively, while for EMC it ranges from % at low to % at intermediate and % at high . Above  GeV/ the signal extraction uncertainty for the EMC arises from the merging of the two photon clusters, and the exact dependence of the corresponding description in the simulation. For the meson the signal extraction uncertainty generally is larger since the signal-to-background ratio is smaller, particularly at low . For PCM the uncertainty is %, for PCM-EMC it ranges from % to % to % and for EMC from % to % GeV/ at low, intermediate and high , respectively. Unlike in the case of the , the uncertainty for EMC decreases with increasing since the merging of the clusters for the meson only sets in at much higher  (around  GeV/). For the  ratio, the signal extraction uncertainties of the and mesons contribute independently.

Efficiency: The uncertainties on the efficiency were estimated using different MC generators to vary the input spectrum for the efficiency calculation, to quantify effects affecting the resolution. Also, the uncertainties on the modeling of the efficiency bias in the simulation were included. For PCM-EMC and EMC the uncertainties range from % to % depending on for the , while they are between % and % for the meson. For the  measurement, the uncertainties were added quadratically, without including the trigger-related uncertainties, which largely cancel. In the case of mEMC, the uncertainty on the resolution is particularly important, since it strongly depends on whether the neutral pion could be reconstructed with all decay particles contributing to the single cluster or just some of them. To estimate the uncertainty due to a possible imperfection of the MC simulation in the contribution of the various possibilities, the fractions of the respective reconstruction possibilities were varied by 20% each, leading to an uncertainty on the efficiency of % at mid ( GeV/) and % at high  ( GeV/).

Secondary correction: The correction for secondary was estimated applying the efficiency and acceptance from the full ALICE GEANT3 simulation to a fast MC simulation of the decay kinematics based on the parametrized K (K) and spectra [Abelev:2014laa]. The corresponding uncertainty was obtained by varying the kaon and yield within their measured uncertainties. Since the correction due to the secondaries is only %, for all but the mEMC reconstruction technique, even a variation of % on the input yields leads to a negligible contribution compared to other uncertainties. For mEMC, where the correction is about %, an uncertainty of was obtained. In addition, % were added to the uncertainty to account for the limited precision in the shape and size of the correction factors of the full simulations for the pions from K, K and , which was estimated by varying the parametrization underlying the efficiencies for secondary .

Inner material: The uncertainty related to the knowledge of the inner (radius ) material budget reflects the uncertainty of the conversion probability of photons, and hence dominantly affects the PCM measurements. It was estimated to be independent of based on detailed comparison between simulation and data for pp collisions at  [Abelev:2012cn]. Thus, it affects the PCM meson measurements with %, while it only contributes 4.5% to PCM-EMC. In , the uncertainty cancels as both mesons are affected in the same way.

Outer material: For the reconstructed photons in the EMCal, a possible mismatch between the material present in reality and assumed in the simulation in front of the EMCal may cause an error in the absorption rate or the production of secondary pions. In most cases, however, the photon simply converts and at least one of its daughter electrons can be reconstructed in the EMCal so that the likely will be reconstructed as well, although with degraded resolution. The probability to still reconstruct the neutral meson increases with increasing conversion radius, i.e. the closer the conversion happens to the surface of the EMCal. Most of the material is located at most  m away from the EMCal, namely the TPC outer wall, the Transition Radiation Detector (TRD) and the Time-Of-Flight (TOF) detector plus their support structures. The TRD was only fully installed in the LHC shutdown period after 2013. For the 2011 and 2013 data there were regions in without TRD modules in front of the EMCal. Hence, the net-effect of the material in front of the EMCal could be studied by comparing fully corrected yields for different regions with or without the TRD in front of the EMCal. From the observed difference measured using the EMC and PCM-EMC measurements, an uncertainty on the neutral meson yields of % independent of was derived, and assigned to all measurements involving the EMCal. For  the uncertainty is assumed to cancel as both mesons should be affected in a similar way.

Trigger normalization and pileup: The uncertainties for the trigger normalization were calculated by varying the range for the fit of the plateau region (see Fig. 1) for the different trigger combinations, leading to the respective rejection factors with their uncertainties given in Tab. 1. Since the final spectra for each measurement technique using the EMCal are composed of several triggers, the contributions of the respective trigger rejection uncertainties enter the final measurement with different magnitudes depending on . The uncertainties range between % and %. For  the uncertainties cancel as the ratio was measured per trigger and reconstruction method and combined afterwards. For PCM only minimum bias triggers were used, and hence no uncertainty due to the trigger rejection was assigned. However, an uncertainty of % to % was taken into account for the out-of-bunch pileup subtraction described in [Abelev:2014ypa]. The pileup uncertainty is about % for the meson. It largely cancels in the  ratio, however, and the remaining error can be neglected compared to other error sources.

Meson (pb GeV) (GeV) (pb GeV) (GeV)
Table 7: Parameters of the two-component model, Eq. 7 [Bylinkin:2014fta, Bylinkin:2014qea], which are used to parametrize the neutral pion and meson spectra, respectively, for the comparisons to models and among the different methods.

6 Results

Since the meson measurements with PHOS, PCM, EMC, PCM-EMC and mEMC have partly uncorrelated systematic uncertainties, their combination will increase the precision of the respective cross section measurements. The BLUE method [Lyons:1988rp, Valassi:2003mu, Valassi:2013bga] was used to calculate the combined spectra of the and mesons as well as the  ratio. For the combination of the spectra, the full correlation matrix was taken into account by estimating the correlated and uncorrelated part of the systematics for all pairs of measurements versus . Correlations are most apparent between the three EMC related measurements (EMC, PCM-EMCand mEMC), as well as for the PCM-EMC and PCM results. At high , for instance, the uncertainties are dominated by the uncertainty on which is largely common between the EMCal triggered analyses. Uncertainties between PHOS, PCM, and EMC (mEMC) are uncorrelated. The combined spectra were fitted with a two-component model (TCM)

(7)

introduced by Bylinkin and Rostovtsev [Bylinkin:2014fta, Bylinkin:2014qea], which serves as convenient parametrization of the data without aiming for a physics interpretation. The parameters for the and fits are given in Tab. 7 for values of better than taking statistical and systematic uncertainties in quadrature. Unlike for Tsallis [Tsallis:1987eu] and power-law distributions, which at high and low , respectively, systematically deviate from the data, the TCM parameterization describes the data over the full measured range to better than 10%.

Figure 8 shows a comparison of the individual measurements in their respective measured ranges summarized in Tab. 8 to the two-component model fits for the and mesons. As already mentioned above, the spectrum in pp collisions at TeV has been measured by ALICE using the PHOS and PCM [Abelev:2014ypa]. The new results obtained with the different EMC measurements and with the hybrid PCM-EMC method are consistent with these earlier results, and the combination with the former measurements improves the precision of the data. The figure also demonstrates an approximately fourfold extension of the reach of the measurement by using the EMCal. The measurement, which is the first such measurement at TeV, spans from  GeV/ to  GeV/. There is good agreement within the statistical uncertainties among the different detection techniques. Above  GeV/, the result is dominated by the EMCal measurements.

Fig. 8: Comparison of the individual measurements in their respective measured transverse momentum ranges relative to the two-component model fits [Bylinkin:2014fta, Bylinkin:2014qea] of the final spectra. The final spectra are obtained by combining the individual measurements in the overlapping  regions with the highest granularity using the full correlation matrix as defined in the BLUE-algorithm [Lyons:1988rp, Valassi:2003mu, Valassi:2013bga].
Method
PCM 0.4–8.0 0.5–6.0 0.5–6.0
PHOS 0.8–12.0 n/a n/a
EMC 1.4–20.0 2.0–20.0 2.0–20.0
PCM-EMC 0.8–20.0 1.0–16.0 1.0–16.0
mEMC 16.0–40.0 n/a n/a
Table 8: Summary of the  reach (in GeV/) of the various reconstruction methods for , and .
Fig. 9: Invariant differential cross section of the (left, top panel) and meson (right, top panel) for pp collisions at TeV. The data are compared to PYTHIA 8.2 [Sjostrand:2014zea] generator-level simulations using the Monash 2013 tune as well as recent NLO pQCD calculations [deFlorian:2014xna, Aidala:2010bn]. The ratios of the data and the calculations to the respective two-component model fits [Bylinkin:2014fta, Bylinkin:2014qea] to the data are shown in the lower panels. The horizontal error bars denote statistical, the boxes systematic uncertainties.

Figure 9 shows the combined and cross sections in pp collisions at  TeV, and Fig. 10 the corresponding  ratio. As mentioned earlier, the data were parameterized with a two-component model of Bylinkin and Rostovtsev [Bylinkin:2014qea] (see Tab. 7) and compared to recent NLO pQCD calculations [deFlorian:2014xna, Aidala:2010bn], and PYTHIA 8.2 [Sjostrand:2014zea] generator-level simulations using the widely-used Monash 2013 tune [Skands:2014pea]. A large fraction of hadrons at low is produced in pp collisions via soft parton interactions and from resonance decays, which cannot be well described within the framework of pQCD, but are taken into account in the event-generator approach. For the , the pQCD calculation [deFlorian:2014xna], which uses the DSS14 fragmentation functions seems to have a different shape than the data. It overpredicts the data by about % at intermediate ( GeV/ GeV/), while it agrees with the data at higher . The PYTHIA 8.2 calculation describes the data well, except below  GeV/, where it overpredicts the data by up to %. For above  GeV/ PYTHIA has a tendency to underpredict the data by about %; however this slight difference is covered by the uncertainties of the measurement. For the meson, the data and the NLO pQCD calculation [Aidala:2010bn], which uses the AESSS fragmentation functions, agree within the uncertainties for for factorization and fragmentation scale, while for the calculation overpredicts the data by up to a factor of , leaving room for future improvements in the understanding of the strange versus non-strange quark fragmentation functions. The PYTHIA 8.2 simulation with the Monash 2013 tune performs slightly worse for the than for the , in particular for  GeV/ where it underpredicts the data by about %. In the  ratio, parts of the systematic uncertainties cancel not only for the data but also for the NLO pQCD calculation. Thus, even the predictions using older fragmentation functions for the  [deFlorian:2007aj] and the  [Aidala:2010bn], which can not reproduce the individual spectra [Abelev:2014ypa], are in good agreement for the  measurement. PYTHIA 8.2 using the Monash 2013 tune can reproduce the dependence of the ratio; however it underpredicts the ratio by about % above  GeV/, albeit still in agreement with the data to within . The measured  ratio is found to agree with previous measurements in pp collisions at TeV [Adare:2010cy] and TeV [Abelev:2012cn] suggesting that is collision-energy independent. Above  GeV/, both mesons exhibit a similar power-law behavior with and with of below . This is also reflected in the  ratio, which above  GeV/ reaches a value of .

Fig. 10: Measured ratio in pp collisions at TeV compared to NLO pQCD calculations [deFlorian:2007aj, Aidala:2010bn] and PYTHIA 8.2  [Sjostrand:2007gs] generator-level simulations using the Monash 2013 tune. The horizontal error bars denote statistical, the boxes systematic uncertainties. The data at  TeV [Adare:2010cy] and  TeV [Abelev:2012cn] are shown with statistical and systematic uncertainties added in quadrature.

7 Summary

The invariant differential cross sections for inclusive and production at midrapidity in pp collisions at  TeV were measured over a large range in transverse momentum of  GeV/ and  GeV/, respectively. To achieve these measurements, for the  () five (three) different reconstruction techniques and multiple higher-level triggers involving the EMCal in ALICE were exploited. In particular, a new single-cluster, shower-shape based method was developed to identify high- neutral pions whose decay photons overlap in the EMCal. Above  GeV/, both the and cross sections are found to exhibit a similar power-law behavior with an exponent of about . The data were compared to state-of-the-art NLO pQCD calculations which are found to reproduce the neutral pion cross section within , while the deviations for the meson are significantly larger. Calculations using PYTHIA 8.2 at generator-level with the Monash 2013 tune turn out to be consistent with the measurement, except below  GeV/, where the calculation overpredicts the data by up to %. For the , the agreement is slightly worse than for the , in particular for  GeV/ where the calculation underpredicts the data by about %. The  ratio, which was found to be described by the calculations to within , is above  GeV/, consistent with previous measurements. The new data provide significant constraints for future calculations of hadron spectra over a large range in .

Acknowledgments

We thank Werner Vogelsang and Marco Stratmann for providing us with their calculations.

The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Universidade Federal do Rio Grande do Sul (UFRGS), Financiadora de Estudos e Projetos (Finep) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil; Ministry of Science & Technology of China (MSTC), National Natural Science Foundation of China (NSFC) and Ministry of Education of China (MOEC) , China; Ministry of Science, Education and Sport and Croatian Science Foundation, Croatia; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research | Natural Sciences, the Carlsberg Foundation and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Energie Atomique (CEA) and Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; Ministry of Education, Research and Religious Affairs, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE) and Council of Scientific and Industrial Research (CSIR), New Delhi, India; Indonesian Institute of Science, Indonesia; Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionale di Fisica Nucleare (INFN), Italy; Institute for Innovative Science and Technology , Nagasaki Institute of Applied Science (IIST), Japan Society for the Promotion of Science (JSPS) KAKENHI and Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nationaal instituut voor subatomaire fysica (Nikhef), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education and National Science Centre, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics and Romanian National Agency for Science, Technology and Innovation, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education and Science of the Russian Federation and National Research Centre Kurchatov Institute, Russia; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba, Ministerio de Ciencia e Innovacion and Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Spain; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; National Science and Technology Development Agency (NSDTA), Suranaree University of Technology (SUT) and Office of the Higher Education Commission under NRU project of Thailand, Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America.

References

A The ALICE Collaboration

S. Acharya, D. Adamová, M.M. Aggarwal, G. Aglieri Rinella, M. Agnello, N. Agrawal, Z. Ahammed, N. Ahmad, S.U. Ahn, S. Aiola, A. Akindinov, S.N. Alam, D.S.D. Albuquerque, D. Aleksandrov, B. Alessandro, D. Alexandre, R. Alfaro Molina, A. Alici,,, A. Alkin, J. Alme, T. Alt, I. Altsybeev, C. Alves Garcia Prado, M. An, C. Andrei, H.A. Andrews, A. Andronic, V. Anguelov, C. Anson, T. Antičić, F. Antinori, P. Antonioli, R. Anwar, L. Aphecetche, H. Appelshäuser, S. Arcelli, R. Arnaldi, O.W. Arnold,, I.C. Arsene, M. Arslandok, B. Audurier, A. Augustinus, R. Averbeck, T. Awes, M.D. Azmi, A. Badalà, Y.W. Baek, S. Bagnasco, R. Bailhache, R. Bala, A. Baldisseri, M. Ball, R.C. Baral, A.M. Barbano, R. Barbera, F. Barile,, L. Barioglio, G.G. Barnaföldi, L.S. Barnby,, V. Barret, P. Bartalini, K. Barth, J. Bartke,i, E. Bartsch, M. Basile, N. Bastid, S. Basu, B. Bathen, G. Batigne, A. Batista Camejo, B. Batyunya, P.C. Batzing, I.G. Bearden, H. Beck, C. Bedda, N.K. Behera, I. Belikov, F. Bellini, H. Bello Martinez, R. Bellwied, L.G.E. Beltran, V. Belyaev, G. Bencedi, S. Beole, A. Bercuci, Y. Berdnikov, D. Berenyi, R.A. Bertens,, D. Berzano, L. Betev, A. Bhasin, I.R. Bhat, A.K. Bhati, B. Bhattacharjee, J. Bhom, L. Bianchi, N. Bianchi, C. Bianchin, J. Bielčík, J. Bielčíková, A. Bilandzic,, G. Biro, R. Biswas, S. Biswas, J.T. Blair, D. Blau, C. Blume, G. Boca, F. Bock,, A. Bogdanov, L. Boldizsár, M. Bombara, G. Bonomi, M. Bonora, J. Book, H. Borel, A. Borissov, M. Borri, E. Botta, C. Bourjau, P. Braun-Munzinger, M. Bregant, T.A. Broker, T.A. Browning, M. Broz, E.J. Brucken, E. Bruna, G.E. Bruno, D. Budnikov, H. Buesching, S. Bufalino, P. Buhler, S.A.I. Buitron, P. Buncic, O. Busch, Z. Buthelezi, J.B. Butt, J.T. Buxton, J. Cabala, D. Caffarri, H. Caines, A. Caliva, E. Calvo Villar, P. Camerini, A.A. Capon, F. Carena, W. Carena, F. Carnesecchi,, J. Castillo Castellanos, A.J. Castro, E.A.R. Casula,, C. Ceballos Sanchez, P. Cerello, B. Chang, S. Chapeland, M. Chartier, J.L. Charvet, S. Chattopadhyay, S. Chattopadhyay, A. Chauvin,, M. Cherney, C. Cheshkov, B. Cheynis, V. Chibante Barroso, D.D. Chinellato