Production of Sigma (1385)(+/-) and Xi (1530)(0) in proton-proton collisions at root s=7

The production of the strange and double-strange baryon resonances ( (cid:3)( 1385 ) ± , (cid:4)( 1530 ) 0 ) has been measured at mid-rapidity ( | y | < 0 . 5) in proton–proton collisions at √ s = 7 TeV with the ALICE detector at the LHC. Trans-verse momentum spectra for inelastic collisions are compared to QCD-inspired models, which in general underpre-dict the data. A search for the φ( 1860 ) pentaquark, decaying in the (cid:4)π channel, has been carried out but no evidence is seen.


Introduction
The study of strange baryon resonances in proton-proton (pp) collisions contributes to the understanding of hadron production mechanisms and provides a reference for tuning QCD-inspired event generators. The strange-quark content makes these baryons a valuable tool in understanding production mechanisms, since the initial state colliding projectiles contain no strange valence quarks and therefore all strange particles are created in the collision.
In addition, a measurement of resonance production in the pp system serves as a reference for understanding resonance production in heavy-ion collisions, where resonances, due to their lifetime of a few fm/c being comparable to the lifetime of the hadronic phase, are sensitive probes of the dynamical evolution of the fireball. Previous measurements at a collision energy of √ s = 0.2 TeV with the STAR detector at the RHIC have shown that the yields of (1385) in Au-Au in comparison to pp collisions indicate the presence of rescattering and regeneration in the time span between chemical and kinetic freezeout [1]. Forthcoming analysis of strange baryon resonances in Pb-Pb collisions by the ALICE collaboration will further explore those effects at higher energy and density of the colliding system. The results for the (1385) ± and (1530) 0 baryons in pp collisions will therefore serve as benchmark. e-mail: alice-publications@cern.ch Measurements of differential (d 2 N /(dyd p T )) and integrated (dN /dy) yields of the (1385) ± and (1530) 0 baryons are presented at mid-rapidity (|y|< 0.5) in inelastic (INEL) pp collisions at √ s = 7 TeV, collected with the ALICE detector [2] at the LHC. The differential spectra are compared to Monte Carlo (MC) event generators. The mean transverse momentum p T is compared to those of other particles measured in pp collisions with the ALICE detector at both √ s = 7 TeV and √ s = 0.9 TeV, and with the STAR detector at √ s = 0.2 TeV. The (1530) reconstruction channel π is additionally analysed to investigate evidence of the φ(1860) pentaquark, previously reported by the NA49 experiment [3]. No such signal was observed by other experiments at different energies and with different beams and reactions [4][5][6][7][8][9][10][11][12][13][14].
This article is organized as follows. Section 2 gives a brief description of the main detectors used for this analysis and the experimental conditions. Section 2.1 describes track and topological selections. Signal extraction methods are presented in Sect. 2.2, and the efficiency corrections in Sect. 2.3. The evaluation of systematic uncertainties is discussed in Sect. 2.4. In Sect. 3, the p T spectra and the integrated yields of the studied particle species are given and compared to model predictions. In Sect. 4 the search for the φ(1860) pentaquark is discussed. Conclusions are presented in Sect. 5.

Experiment and data analysis
The ALICE detector [2] is designed to study a variety of colliding systems, including pp and lead-lead (Pb-Pb) collisions, at TeV-scale energies. The sub-detectors used in this analysis are described in the following. A six-layer silicon inner tracking system (ITS) [15] and a large-volume time projection chamber (TPC) [16] enable charged particle reconstruction with excellent momentum and spatial resolution in full azimuth down to a p T of 100 MeV/c in the pseudorapidity range |η| < 0.9. The primary interaction vertex is determined with the TPC and ITS detectors with a resolution of 200 μm Table 1 Particles involved in this analysis and their PDG parameters [17]. Antiparticles are not listed for brevity. From [17], (1530) 0 −→ + π has a branching ratio of ∼ 100 %, then (1530) 0 −→ − +π + has a branching ratio of ∼ 66. 7  for events with few tracks (N ch 3) and below 100 μm for events with higher multiplicity (N ch 25). In addition, both detectors are able to provide particle identification (PID) via energy-loss measurements. The data analysis is carried out using a sample of ∼ 250 million minimum-bias pp collisions at √ s = 7 TeV collected during 2010.
During the data-taking period, the luminosity at the interaction point was kept in the range 0.6−1.2 × 10 29 cm −2 s −1 . Runs with a mean pile-up probability per event larger than 2.9 % are excluded from the analysis. The vertex of each collision is required to be within ±10 cm of the detector's centre along the beam direction. The event vertex range is selected to optimize the reconstruction efficiency of particle tracks within the ITS and TPC acceptance.

Particle selections
The resonances are reconstructed via their hadronic decay channel, shown in Table 1 together with the branching ratio (BR).
(1530) 0 is measured together with its antiparticle ( (1530) 0 ) due to limited statistics. Therefore in this paper, unless otherwise specified, Note that, for brevity, antiparticles are not listed and the selection criteria, described in the following, are discussed for particles; equivalent criteria hold for antiparticles.
Several quality criteria, summarized in Table 2, are used for track selection.
Charged pions from the strong decay of both (1385) and (1530) 0 are not distinguishable from primary particles and therefore primary track selections are used. They are requested to have a distance of closest approach (DCA) to the primary interaction vertex of less than 2 cm along the beam direction and a DCA in the transverse plane smaller than 7 σ DCA ( p T ), where σ DCA ( p T ) = (0.0026 + 0.0050 GeV/c × p T −1 ) cm is the parametrization which accounts for the p Tdependent resolution of the DCA in the transverse plane [18]. Primary tracks are also required to have at least one hit in one of the two innermost layers of the ITS (silicon pixel detector, SPD) and at least 70 reconstructed clusters in the TPC out of the maximum 159 available, which keeps the contamination from secondary and fake tracks small, while ensuring a high efficiency and good dE/dx resolution. Tracks close to the TPC edge or with transverse momentum p T < 0.15 GeV/c are rejected because the resolution of track reconstruction deteriorates.
In the (1385) analysis, PID is implemented for π ± and p from .
Particles are identified based on a comparison of the energy deposited in the TPC drift gas and an expected value computed using a Bethe-Bloch parametrization [19]. The filter is set to 3 σ TPC , where σ is the resolution estimated by averaging over reconstructed tracks. An averaged value of σ TPC = 6.5 % is found over all reconstructed tracks [20].
PID selection criteria are not applied in the (1530) analysis as the combinatorial background is sufficiently removed through topological selection.
produced in the decay of (1385) decays weakly into π − p with cτ = 7.89 cm [17]. These pions and protons do not originate from the primary collision vertex, and thus they  [19]. Selection criteria for used in the (1385) analysis are summarized in Table 3. − produced in the decay of the (1530) 0 decays weakly into π − with cτ = 4.91 cm [17]. Pions are selected from tracks with a DCA to the interaction point greater than 0.05 cm. Pions and protons from are required to have a DCA to the interaction point greater than 0.04 cm. All pions and protons are requested to have at least 70 reconstructed clusters in the TPC. Decay topologies for − and are used as described in [19]. Selection criteria are summarized in Table 4.
All these criteria are optimized to obtain maximum signal significance. Values for the significance are presented in Sect. 2.2.3.

Combinatorial background and event-mixing
Due to their very short lifetime of a few fm/c, resonance decay products originate from a position that is indistinguishable from the primary vertex. Thus, the computation of invariant mass distributions for potential resonance decay candidates has significant combinatorial background that has to be subtracted to ensure reliable yield determination. This is shown in the left panels of Figs Fig. 2 the peak from − −→ + π − is visible. The combinatorial background distributions are obtained and subtracted from the invariant mass distribution by means of a mixed-event technique, in which a reference background distribution is built with uncorrelated candidates from different events. To avoid mismatch due to different acceptances and to ensure a similar event structure, only tracks from events with similar vertex positions z ( z < 1 cm) and track multiplicities n ( n < 10) are mixed. In order to reduce statistical uncertainties, each event is mixed with several other events (5 in the (1385) analysis and > 20 in the (1530) 0 analysis), so that the total number of entries in the mixed-event invariant mass distribution is higher than the total number of entries in the distribution from the same event. Thus the mixed-event distribution needs to be scaled before it can be used to describe the background in the sameevent distribution. For (1385), the regions for the normalization of the mixed-event distribution are selected in the rightmost part of the invariant mass window, where the residual background is absent (see Sect. 2.2.2 for a description of the residual background). These regions are different for the different p T bins, ranging from 1.48 < M < 2.0 GeV/c 2 , for the lowest p T bin, to 1.95 < M < 2.0 GeV/c 2 , for the highest p T bin (M being the invariant mass of (1385) and 2.0 GeV/c 2 being the upper extreme of the invariant mass window). The reason for this p T -dependent choice is due to the reach of the residual background, which is higher in invariant mass for higher p T . Fixed regions, 1.6 < M < 1.8 GeV/c 2 and 1.8 < M < 2.0 GeV/c 2 , have also been tried, giving a systematic uncertainty of ∼ 1 %. For (1530) 0 a fixed region 1.49 < M < 1.51 GeV/c 2 , just at the left of the signal, is selected. A fixed region can be selected because for all p T intervals the background shape is similar and the invariant mass resolution on the reconstructed peak is the same. The uncertainty in the normalization (∼ 1 %), which is included in the quoted systematic uncertainty for signal extraction, is estimated by using another normalization region, 1.56 < M < 1.58 GeV/c 2 , just at the right of the signal. The open squares in the left panels of Figs. 1, 2 and 3 correspond to the properly scaled mixed-event invariant mass distribution.
The right panels show the signals for each resonance after the mixed-event combinatorial background is subtracted.

Same-event pairs
Mixed-event background

Residual correlated background
The mixed-event technique removes only uncorrelated background pairs in the invariant mass spectrum. The consequence is that residual correlations near the signal mass range are not subtracted by the mixed-event spectrum and correlated background pairs remain [21]. This is especially dominant for (1385) (see Figs. 1, 2, right), for which the correlated residual background takes contributions from two dominant sources: Table 5 Potential sources of contamination in the reconstruction of (1385). Checkmarks show which species is potentially affected. Checkboxes further indicate whether the source gives a significant contamination (see text). A similar scheme, not shown for sake of brevity, is valid for the antiparticles • Type A: correlated π pairs coming from the decays of other particles which have and π among the decay products. • Type B: correlated π pairs which come from the dynamics of the collision and are not removed from the subtraction of the mixed-event background.
All these contributions are present in the MC, albeit with potentially incorrect proportions. Thus, simulations are used to determine the shapes of such contributions in invariant mass space and then these contributions are renormalized using data, as described later.
All the sources of contamination of Type A, which can potentially produce correlated π pairs, are listed in Table 5. A similar scheme, not discussed for sake of brevity, is valid for the antiparticles (e.g. the + −→ π + decay channel affects the reconstruction of (1385) + ). Only sources A1, A5 and A6 in Table 5 give a significant contribution to the correlated residual background of Type A. This is discussed in the following. Source A1 in Table 5 is due to the primary − which decays weakly to π − , affecting the reconstruction of (1385) − . Since the − hyperon is metastable, it shows up in the π − invariant mass spectrum as a very narrow peak at around the − mass, M − = 1321.71 MeV/c 2 [17], just on the left tail of the (1385) − signal. The − peak is clearly seen in Fig. 2. This contribution, which is expected to be important since the yield of − is comparable to the yield of (1385) − , is in fact suppressed, by an order of magnitude, because the filter on the DCA to the primary vertex of both and π filters out most of the π pairs from − . Indeed, the filter on the DCA to the primary vertex is optimized for the (1385) decay products, which are not distinguishable from primary particles (see Sect. 2.1), whereas and π from − come from a secondary vertex, centimetres away from the primary vertex. Only a small percentage of the − yield survives the filter on the DCA. Source A1 is taken into account by adding a Gaussian function, with the mean value fixed to the − mass and the width and normalization left free, to the combined fit of the invariant mass spectrum in the reconstruction of (1385) − . The contamination from − reaches about 5-10 % of the raw (1385) − signal and varies little with p T .
Sources A2, A3 and A4 give a negligible contribution. Sources A2 and A3 are due to the hadronic decay channels of (1530) − , with BR = 33.3 % and BR = 66.7 %, respectively 1 , and, like A1, affect only the (1385) − reconstruction. Source A4 is due to (1530) 0 and potentially affects the reconstruction of both (1385) + and (1385) − , since it involves two opposite-sign pions. The same topological considerations hold for A2 as they do for A1, since it involves a − . Indeed, this − comes from the strong decay of (1530) − , therefore it is practically not distinguishable from the (primary) − in A1. Unlike contribution A1, a further suppression, of about an order of magnitude with respect to A1, comes from both the smaller yield of (1530) − with respect to the primary − , and the BR of the (1530) − → − π 0 channel. This further suppression makes contribution A2 practically negligible. Similar conclusions hold for contributions A3 and A4.
Source A5 in Table 5 is related to the second (1385) decay channel, (1385) ± → 0 π ± (BR = 5.8 % 2 ), with 0 → γ (BR 100 % [17]). from 0 is paired with π ± from (1385) ± . This gives a Gaussian-like peak at around 1.306 GeV/c 2 , with a width of ∼ 0.059 GeV/c 2 (FWHM). This peak is used in the combined fit to the signal (see below) with a relative normalization with respect to the signal which accounts for the ratio (=0.067) between the BR (=5.8 %) for the (1385) ± → 0 π ± channel and the BR (=87 %) for the (1385) ± → π ± channel. Source A6 in Table 5 is due to the (1520) → π ± π ∓ channel (BR = 5 % 3 ). The positive (negative) pion, paired with , produces a Gaussian-like peak, which contaminates the invariant mass distribution of (1385) + ( (1385) − ). This peak is centred at ∼ 1.315 GeV/c 2 and has a width of ∼ 0.076 GeV/c 2 (FWHM). The peak is used in the combined fit to the signal. The normalization of the peak is kept free in the fit since the (1520) yield is not measured. The contamination from (1520) decreases with increasing p T , ranging from about 75 % of the raw (1385) − signal in the first p T interval, down to 0 for p T > 4 GeV/c.
A third-degree polynomial is used to fit the residual background of Type B in the MC. The fit to MC data is performed in the region from 1.26 GeV/c 2 (just left of the signal region) to the lower edge of the event-mixing normalization region. The fitting function is then normalized to the residual background in real data; the normalization is done in the region from 1.46 GeV/c 2 (just right of the signal region) to the lower edge of the event-mixing normalization region, where other sources of contamination are absent. The lower point of the normalization region is the same for all p T intervals since the mean, the width and the invariant mass resolution on the reconstructed peak stay the same over all the p T range considered. Comparable results are obtained from using different event generators (PYTHIA 6.4, tune Perugia 0 [22], and PHOJET [23]) and other degrees for the polynomial (second and fourth). The differences of about 2 % are included in the systematic uncertainties.
The invariant mass distribution is fitted with a combined fit function: a (non-relativistic) Breit-Wigner peak plus the functions that make up the residual background (Figs. 1, 2, right). The Breit-Wigner width is kept fixed to the PDG value to improve the stability of the fit.
For (1530) 0 , the residual background after the mixedevent background subtraction is fitted with a first-degree polynomial. The fitting procedure is done in three stages. First, the background is fitted alone from 1.48 to 1.59 GeV/c 2 while excluding the (1530) 0 mass region from 1.51 to 1.56 GeV/c 2 . Second, a combined fit for signal and background is performed over the full range with the background polynomial fixed to the results from the first fit stage; a Voigtian function-a convolution of Breit-Wigner and Gaussian functions-is used for the signal. The Gaussian part accounts for detector resolution. Third, a fit is redone over the full range again with all parameters free but set initially to the values from the second stage.

Counting signal and signal characteristics
The above procedure is applied for 10 (8) p T bins for (1385) ( (1530) 0 ), from 0.7 to 6.0 (0.8 to 5.6) GeV/c. For (1385), the fit is repeated leaving the Breit-Wigner width free to move, and, for each p T interval, the difference in the yield is included in the systematic uncertainties (∼ 4 % maximum contribution). The widths of both (1385) and (1530) 0 are consistent with the PDG values for all p T intervals. In the (1385) − analysis, a Gaussian function, centred at 1.321 GeV/c 2 and with a starting value for the width of 2 MeV/c 2 , is used to help the combined fit around the (1321) − peak (Fig. 2). The value of 2 MeV/c 2 is obtained from the analysis of (1321) − [19] and is related to the mass resolution. Since the (1385) mass binning of 8 MeV/c 2 , which is optimised for the χ 2 of the combined fit, is larger than the mass resolution, only a rough description of the (1321) − peak is possible. For (1530) 0 , the standard deviation of the Gaussian component of the Voigtian peak is found to be ∼ 2 MeV/c 2 , which is consistent with the detector resolution, as obtained from the MC simulation. At low p T , the fitted mass values for (1385) are found to be slightly lower (by ∼ 5 MeV/c 2 ) than the PDG value, which is attributed to imperfections in the corrections for energy loss in the detector material. For (1530) 0 , the reconstructed masses are found to be in agreement with the PDG value within the statistical uncertainties.
The raw yields N RAW are obtained by integrating the Breit-Wigner function. As an alternative, N RAW is calculated by integrating the invariant mass histogram after the subtraction of the event-mixing background and subtracting the integral of the residual background (bin-counting method). The difference between the two methods of integration is lower than 2 % on average.
Significance values (defined as S/ √ S + B, where S is the signal and B the background) for (1385) + ( (1530) 0 ) are found to be 16.6 (16.5) in the lowest p T interval, and 20.9 (22.8) in the highest p T interval, and reached 24.2 (52.4) in the intermediate p T interval. Significance values comparable to those of (1385) + are obtained for the other (1385) species.

Correction and normalization
In order to extract the baryon yields, N RAW are corrected for BR, the geometrical acceptance ( A), the detector efficiency ( ) and the correction factor which accounts for the GEANT3 overestimation of thep cross sections ( GEANT3/FLUKA ) [24] The product of acceptance and efficiency (A × ) is determined from MC simulations with the PYTHIA 6.4 event generator (tune Perugia 0 [22]) and a GEANT3-based simulation of the ALICE detector response [25]. The GEANT3/FLUKA correction factor is equal to 0.99 for the protons from   Fig. 4 The product of acceptance, efficiency and branching ratio of (1385) + and (1530) 0 , obtained with PYTHIA 6.4 [22] and GEANT3 [25], as function of p T in |y|<0.5. Only statistical uncertainties are reported. The dashed-and the dash-dotted lines indicate the overall branching ratio for the two reconstruction channels those generated in |y|< 0.5, differentially as a function of transverse momentum, as shown in Fig. 4.
The drop in efficiency at low p T is due to the loss of slow pions involved in the decay chain. As a cross-check, the efficiency × acceptance has also been assessed with PHO-JET [23] as event generator. The relative difference of the resulting A × , averaged over the various p T intervals, is below 1 %.
Finally, corrections for the trigger inefficiency ( trigger ) and the loss of candidates outside of the z-vertex range ( vert ) are applied via where N cor and N MB are the number of reconstructed (1385) or (1530) and the total number of minimum bias triggers, respectively. y and p T are the rapidity window width and the p T bin width, respectively. The trigger selection efficiency for inelastic collisions trigger is equal to 0.852 +0.062 −0.030 [26]. The loss of resonances due to the trigger selection, estimated by MC simulations, is negligible, less than 0.2 %. The vert correction factor accounts for resonance losses (∼ 7 %) due to the requirement to have a primary vertex z position in the range ±10 cm.

Systematic uncertainties of p T spectra
Two types of systematic uncertainties in the particle spectra are considered: p T -dependent systematic uncertainties, which are due to the selection efficiency and signal extraction at a given p T , and p T -independent uncertainties due to the normalization to inelastic collisions and other corrections. The minimum and maximum values of the major contributions to the point-to-point systematic uncertainties are listed in Table 6.
The uncertainties introduced by tracking, topology selection and PID are obtained by varying the selection criteria for the decay products. To this purpose, the selection criteria listed in Tables 2, 3 and 4 are changed by a certain amount which varies the raw yield in real data by ±10 %. The maximum difference between the default yield and the alternate value obtained by varying the selection, is taken as systematic uncertainty. The uncertainties introduced by the signal extraction come from several sources: normalization of the event-mixing background, fitting function and range of the residual background, signal fitting and integration. For (1385), the contamination from the (1520) introduced the largest contribution (∼ 8 %). All the sources are combined by summing in quadrature the uncertainties for each p T .
Among the p T -independent uncertainties, the INEL normalization leads to a +7.3 % and −3.5 % uncertainty [26], the determination of the material thickness traversed by the particles (material budget) introduces a 4 % uncertainty and the use of FLUKA [27,28] to correct the antiproton absorption cross section in GEANT3 leads to a further 2 % uncertainty [24]. For (1385), a further 1.5 % comes from the uncertainty in the branching ratio. A summary of the p Tindependent uncertainties is presented in Table 6.
The vertical error bars in Fig. 5 represent the sum in quadrature of the statistical and systematic uncertainties, excluding the p T -independent uncertainties, which affect only the normalization. All spectra are fitted with a Lévy-Tsallis function [29], which is used for most of the identified particle spectra in pp collisions [19,20,[30][31][32], where m T = m 2 0 + p 2 T and m 0 denotes the PDG particle mass. This function, quantified by the inverse slope parameter C and the exponent parameter n, describes both the exponential shape of the spectrum at low p T and the power law distribution at large p T . The parameter dN /dy represents the particle yield per unit rapidity per INEL event. dN /dy, C and n are the free parameters considered for this function. Table 7 presents the parameter outcome of the Lévy-Tsallis fit, together with the mean transverse momentum, p T , and the reduced χ 2 .
The values of dN /dy in Table 7 are obtained by adding the integral of the experimental spectrum in the measured range and the extrapolations with the fitted Lévy-Tsallis function to both p T = 0 and high p T . The contribution of the low-p T extrapolation to the total dN /dy is ∼ 30 % for both (1385) and (1530) 0 . The contribution of the highp T extrapolation is negligible.
For each species considered here, such a composite dN /dy differs very little (< 1 %) from the value of dN /dy as the first free parameter returned by the fit, i.e. from the integration of the fit function from 0 to infinity.
In order to obtain the systematic uncertainty on the parameters of the Lévy-Tsallis fit (dN /dy, C and n) and on the mean transverse momentum ( p T ), the Lévy-Tsallis fit is repeated for each p T spectrum obtained by varying separately the selection criteria in each source of systematic uncertainties. Only statistical uncertainties on the points of the p T spectrum are used for the fit. The values for dN /dy, C, n and p T , obtained for each source, are compared to those from the fit to the reference p T spectrum, obtained with default selection criteria. The fit to the reference p T spectrum is also done with statistical uncertainties only. The statistically significant differences are summed in quadrature to contribute to the overall systematic uncertainties on dN /dy, C, n and p T .
Although the Lévy-Tsallis function describes the spectra both at low and at large p T , other functions (e.g. m T exponential or p T power law) are likely to reproduce the lowp T behaviour and are suitable for the low-p T extrapolation. These functions are fitted to the low-p T part of the spectrum below 3 GeV/c and used to evaluate the low-p T contribution Table 7 Parameters extracted from the Lévy-Tsallis (LT) fits (Eq. 3) to the transverse momentum spectra. The values of dN /dy are calculated using the spectra in the measured range and the extrapolation of the fitted Lévy-Tsallis function outside the measured range. Systematic uncertainties quoted here are the ones derived from Lévy-Tsallis fit only (see text) where A is the normalization factor and C is the inverse slope parameter, gives values for dN /dy which are ∼ 5-6 % lower and values for p T which are ∼ 3 % higher than those obtained with the Lévy-Tsallis function. A p T power law functional form gives values for dN /dy which are ∼ 10-15 % higher and values for p T which are ∼ 9-11 % lower than those obtained with the Lévy-Tsallis function. Arithmetic averages of the values obtained with the three functions (Lévy-Tsallis, m T exponential, p T power law) are taken for dN /dy and p T and the unbiased estimators of standard deviation are considered as systematic uncertainties associated to the low-p T extrapolation. These systematic uncertainties are summed in quadrature to contribute to the overall systematic uncertainties on dN /dy and p T . Table 8 summaries the results. The anti-baryon to baryon ratios, (1385) − / (1385) + and (1385) + / (1385) − , are compatible with unity, although the large uncertainties leave very little predictive power on the mechanisms of baryon-number transport [33].
The latest release of PYTHIA 6 (6.427) is used. One of its latest tunes (Perugia 2011, tune 350 [22]) is compared with The transverse momentum spectrum of (1385) + is compared to standard tunes of PYTHIA 6 [34] and PYTHIA 8 [35], the latest release of HERWIG (6.521) [36], and SHERPA release 1.4.6 [37]. The MC data are binned according to the data. Spectra points are represented at the centre of the p T interval. The lower panel shows the ratio data/MC. p T -independent uncertainties are not shown takes into account some of the early LHC minimum-bias and underlying-event data at 0.9 and 7 TeV (see [22] and references therein) and describes the 7 TeV pp charged particle spectra reasonably well [30]. The multi-strange baryon yields are also better described by the Perugia 2011 tune, even if it still underpredicts the data [31]. Similar conclusions hold for the strange meson resonances φ and K * [20]. For both (1385) and (1530) 0 , the Perugia 2011 tune underestimates the data, though it gives a better description with respect to Perugia 0. Also the Perugia 2012 tune of PYTHIA 6 (tune 370 [38]) has been tested with no significant improvement in the predictions for both (1385)  The latest release of PYTHIA 8 (8.176) is used. The standard 4C tune (CTEQ6L1 [35]) gives a worse description with respect to the Perugia 2011 tune of PYTHIA 6. The 4C tune has color reconnection (CR) enabled by default: switching CR off gives a worse description, as expected [39]. ATLAS tunes A2-MSTW2008LO and AU2-CTEQ6L1 have been considered as alternatives to the standard 4C tune (CTEQ6L1). The A2-MSTW2008LO utilizes a different parton distribution function and the AU2-CTEQ6L1 is better tuned for underlying events. None of them performs better than the 4C tune; therefore, they are not reported in Figs. 6 and 7.
Also shown in Figs. 6 and 7 are the results from HER-WIG (release 6.521) [36] and SHERPA (release 1.4.6) [37]. HERWIG predicts a much softer production than for both the other models and the data, for both (1385) and (1530) 0 . For (1385), HERWIG is likely to describe the data at low p T , but it underpredicts the data by a factor ∼ 2−4 in the intermediate-p T region 2 < p T < 3 GeV/c, and more than one order of magnitude at higher p T . For (1530) 0 , HER-WIG fails both at low p T , where the predictions are overestimated by a factor ∼ 2−4, and at high p T , where the predictions are underestimated by more than one order of magnitude. SHERPA gives a better description of the spectral shape for both (1385) and (1530) 0 , but the overall production cross sections are largely underestimated.
The integrated yields dN /dy are also compared to thermal model calculations by Becattini et al. [40], tuned on the yields measured by the ALICE experiment at √ s = 7 TeV for π + , K * 0 , φ, ± and ± [20,30,31], giving a temperature of T = 160 MeV. The other parameters, as obtained from the fit to the ALICE data, are the strangeness suppression factor, γ S = 0.72, the normalization parameter, A = 0.0355, and V T 3 = 231.2, where V is the volume.
The comparison is done for the ratios (1385) + / and (1530) 0 / − , which are sensitive to the temperature T . The experimental yields of and − are from [31,41]

Mean transverse momentum p T
The mean transverse momentum p T serves as a single variable to characterize the soft part of the measured particle spectra. Figure 8 shows the p T as a function of the particle mass, covering a wide range of hadron mass up to the − .
The plot includes (1385) + and (1530) 0 from this analysis, and other particles measured in pp collisions at √ s = 0.9 TeV and √ s = 7 TeV with the ALICE experiment [19,20,[30][31][32]. The STAR pp data at √ s = 0.2 TeV [1,[42][43][44][45] are added for comparison. The dashed line in Fig. 8 is the ISR parametrization, an empirical curve proposed originally [46] to describe the ISR [47] and FNAL [48] data for π , K and p only, at √ s = 0.025 TeV. For STAR data, the ISR parametrization still works relatively well for lower-mass particles up to ∼ 1 GeV/c 2 [44], despite the jump in the collision energy by nearly an order of magnitude with respect to previous experiments, but it fails to describe the dependence of p T for higher-mass particles. At the RHIC energies, this was attributed to an increasing contribution to the transverse momentum spectra from minijet production [49]. In particular, it was noted that strange baryon resonances ( (1385) and (1520)) follow a steeper increase, similar to the trend of heavier mass particles [1].
For ALICE data, the ISR parametrization fails to fit the lower-mass particles already at the collision energy of √ s = 0.9 TeV and the dependence of p T with the mass is even steeper at √ s = 7 TeV. Unlike STAR, strange baryon resonances follow the same trend as the lower-mass particles.
At the LHC energies, flow-like effects in pp collisions are investigated [39,50] which might explain the harder behaviour of transverse momentum spectra, specially for higher mass particles.
The ALICE points at √ s = 7 TeV are fitted with a function similar to the ISR parametrization,  Fig. 8 is the prediction from SHERPA, which is also softer than data.

Search for the φ(1860) pentaquark
In order to explore the existence of the φ(1860) pentaquark, reported by the NA49 experiment [3], the − π + invariant mass spectrum in Fig. 3 was extended up to above 2 GeV/c 2 , as shown in Fig. 9.
The arrow and the shaded area give the region where the φ(1860) pentaquark is expected and where the search was performed. From MC studies with reconstructed particles, the detector mass resolution of the (1530) 0 is ∼ 2 MeV/c 2 and no significant worsening is expected at masses around 1860 MeV/c 2 . The expected theoretical width of the φ(1860) is quite narrow ( 10 MeV/c 2 [3]) so that, eventually, the detector resolution should not affect the measurement. Also in Fig. 9 the like-sign, − π − , invariant mass distribution is presented. Both channels could potentially  Fig. 9 clearly demonstrate the lack of significant evidence for the φ(1860) pentaquark.
No signal of the pentaquark was found by many other experiments [4][5][6][7][8][9][10][11][12][13][14]. A measure of the maximum likely yield of the φ(1860) has been made according to the procedure used by the COMPASS experiment [7]. The background is first estimated by fitting the like-sign distribution with a 4 thorder polynomial from 1.6 to 2.2 GeV/c 2 while excluding the supposed pentaquark range from 1.825 to 1.895 GeV/c 2 . The signal is counted by integrating the entries in the likesign distribution in a 28 MeV/c 2 interval centred around 1.860 GeV/c 2 . The maximum likely signal estimated at the 3 σ (99 %) confidence level is where the counted signal and background are given by s and b, respectively. The ratio of the integrated (1530) 0 yield to the pentaquark yield, S φ(1860) , is to be compared to other experiments. This is shown in Table 9 for the φ(1860) − . The acceptance effects largely cancel in the ratio. The pentaquark search was also performed moving the centre of the search interval by 10 MeV/c 2 to the left and to the right; the same result is obtained. Similar results for S φ(1860) are obtained for φ(1860) 0 .

Conclusions
The transverse momentum spectra of the baryon resonances (1385) and (1530) 0 have been measured by the ALICE collaboration in pp collisions at an energy in the centre of mass of √ s = 7 TeV. A Lévy-Tsallis function describes the spectra well. The mean transverse momentum p T of both (1385) and (1530) 0 , when plotted as a function of the particle mass, follows the trend of other particles measured with the ALICE experiment in pp collisions at √ s = 7 TeV. The differential spectra have been compared to several MC event generators, e.g. standard tunes of PYTHIA 6 and PYTHIA 8, HERWIG and SHERPA. PYTHIA 6 Perugia 2011 (tune 350) performs better than any other tested generator, still underpredicting the data by a factor ∼ 2-3 in the intermediate-p T region 2 < p T < 3 GeV/c.
The search for the φ(1860) 0 and φ(1860) − pentaquark states in the π charged channels has shown no evidence for the existence of such exotic particles.
Acknowledgments 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 acknowledges the following funding agencies for their support in building and running the ALICE detector: State Committee of Science, World Federation of Scientists (WFS) and Swiss Fonds Kidagan, Armenia, Conselho Nacional de Desenvolvimento Científico e Tecnológico Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3