$\Sigma(1385)^{\pm}$ resonance production in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV

Hadronic resonances are used to probe the hadron gas produced in the late stage of heavy-ion collisions since they decay on the same timescale, of the order of 1 to 10 fm/$c$, as the decoupling time of the system. In the hadron gas, (pseudo)elastic scatterings among the products of resonances that decayed before the kinetic freeze-out and regeneration processes counteract each other, the net effect depending on the resonance lifetime, the duration of the hadronic phase, and the hadronic cross sections at play. In this context, the $\Sigma(1385)^{\pm}$ particle is of particular interest as models predict that regeneration dominates over rescattering despite its relatively short lifetime of about 5.5 fm/$c$. The first measurement of the $\Sigma(1385)^{\pm}$ resonance production at midrapidity in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}= 5.02$ TeV with the ALICE detector is presented in this Letter. The resonances are reconstructed via their hadronic decay channel, $\Lambda\pi$, as a function of the transverse momentum ($p_{\rm T}$) and the collision centrality. The results are discussed in comparison with the measured yield of pions and with expectations from the statistical hadronization model as well as commonly employed event generators, including PYTHIA8/Angantyr and EPOS3 coupled to the UrQMD hadronic cascade afterburner. None of the models can describe the data. For $\Sigma(1385)^{\pm}$, a similar behaviour as ${\rm K}^{*} (892)^{0}$ is observed in data unlike the predictions of EPOS3 with afterburner.


Introduction
Relativistic heavy-ion (A-A) collisions provide an excellent tool to study nuclear matter under extreme conditions of temperature and density and the phase transition between hadronic matter and a deconfined state of quarks and gluons, the quark-gluon plasma (QGP), predicted by lattice quantum chromodynamics (QCD) calculations [1].As the QGP produced in a heavy-ion collision expands, it cools down until a phase transition occurs that confines quarks and gluons inside hadrons around a temperature of 155-158 MeV [2].Soon after hadronisation, the resulting dense and hot gas of stable hadrons and resonances reaches the chemical freeze-out.Afterwards, hadrons keep interacting (pseudo-)elastically, and thus exchanging momentum, until the final decoupling at kinetic freeze-out where their momentum distribution is determined.In this stage, the hadron gas temperature decreases from around 150 MeV to approximately 100 MeV.The late hadronic stage of heavy-ion collisions represents a unique environment for the study of a hadronic system in such conditions of high temperature and density.The understanding of the processes occurring in the excited hadron gas is therefore of primary relevance for the interpretation of observables that are employed to characterize the chemical and kinetic freeze-outs, the hadronic phase and more broadly speaking, the time evolution of a heavy-ion collision.In this respect, short-lived resonances are sensitive probes to study the hadronic phase, which is formed after the chemical freeze-out.This is due to their different lifetimes that are comparable to the expected time duration of the hadronic phase (≈ 1-10 fm/c) [3].
The resonances with the shortest lifetime such as the ρ 0 (lifetime cτ ≈ 1 fm [4]) decay during the hadronic phase and their decay products are subject to rescattering effects in the dense hadronic matter, losing memory of the resonance from which they have decayed.Due to the lack of correlation among the decay products, the result is that the resonance cannot be reconstructed via the usual invariant-mass analysis and the measured yield is suppressed with respect to expectations [5].The resonances with the longest lifetime such as the φ (cτ ≈ 45 fm [4]) are likely to survive the hadronic phase and decay in vacuum after the kinetic freeze-out.The scenario is further complicated by regeneration effects for which two particles in the hadronic phase interact via resonance formation, producing an enhancement with respect to the primordial production of that resonance.The Σ(1385) ± (cτ ≈ 5.01 fm for Σ(1385) − and cτ ≈ 5.48 fm for Σ(1385) + [4]) provides a crucial test case in this context.The typical modeling of the hadronic phase is based on transport codes in which hadronic cross sections are implemented to describe the known hadronic interactions.The widely used hadronic cascade simulator UrQMD [6] is therefore coupled to event generators such as EPOS [7,8] in order to predict the production of hadronic resonances.So far, this picture has been tested mostly with the K * (892) 0 resonance which has a lifetime of cτ ≈ 4.2 fm and its production has been measured to be strongly suppressed [9] due to the dominance of rescattering effects.Results from ρ 0 and Λ(1520) further support this picture of rescattering [5,10].Despite its relatively short lifetime, the Σ(1385) ± is expected to be essentially unaffected by the hadronic phase, likely due to larger cross sections for the regeneration processes [11,12].In this Letter, the first measurement of Σ(1385) ± in heavy-ion collisions at the LHC is reported, based on a data sample of Pb-Pb collisions at a centre-of-mass energy per nucleon pair of √ s NN = 5.02 TeV.The measurement presented in this Letter provides an essential test of our current understanding of the evolution of hadronisation induced by heavy-ion collisions.This Letter is organised as follows.After a brief description of the ALICE detector and the data analysis in Section 2, the measured transverse-momentum (p T ) spectra and p T -integrated yields and ratios are presented in Section 3 together with a discussion of our findings.Our conclusions are then summarized in Section 4.
2 Experiment and data analysis

Experimental setup and event selection
A detailed description of the ALICE detector and its performance can be found in Refs.[13,14].This apparatus is optimized for providing particle identification (PID) in a wide momentum range (0.1-20 GeV/c) and high track-density environment by using different techniques.For this analysis, the Inner Tracking System (ITS) and the Time Projection Chamber (TPC) are used for vertex determination and tracking, while the TPC and the Time-Of-Flight (TOF) systems are employed for PID.These detectors cover the full azimuth over a pseudorapidity region |η| < 0.9 and are located inside a large solenoidal magnet providing a field of 0.5 T. The ITS [15], located at a radial distance 3.9 < r < 43 cm from the beam axis, consists of two layers of silicon pixels (SPD), two layers of silicon drift chambers and two layers of silicon strips.The SPD, in particular, is used to reconstruct the track segments that serve to determine the primary vertex of the collision.The TPC [16] is a large cylindrical drift chamber covering a radial distance 85 < r < 247 cm and the main tracking device in the central barrel.The TOF [17] consists of a cylindrical array of MRPCs located at a radial distance of about 380 cm from the beam axis, with an intrinsic resolution of 50 ps.Charged particles can be identified via their specific energy loss, dE/dx, measured in the TPC with a resolution of 5%, and via their time-of-flight measured by TOF.
The data samples analysed were recorded in 2018 using Pb-Pb collisions at √ s NN = 5.02 TeV.The minimum bias (MB) interaction trigger during the data collection requires a coincidence of signals in the V0A and V0C.The central and semi-central triggers require in addition to the MB trigger, an online selection on the total signal amplitude in the V0 detectors, corresponding to collision centralities of 0-10% and 30-50%, respectively.Background events are rejected using the timing information from the ZDC and the V0 detectors [20].For the pile-up removal, the correlations between the number of reconstructed space points (clusters) in the SPD, the number of reconstructed SPD track segments and the total signal in the V0A and V0C detectors are considered.
Events are selected only if they contain at least two tracks and a reconstructed primary vertex located within ±10 cm with respect to the nominal interaction point along the beam axis (z axis).A total of 26 × 10 7 collision events have been accepted for the analysis.Events are further classified into centrality classes, expressed in terms of percentiles of the total hadronic Pb-Pb cross section, using the signal amplitudes in the V0 arrays [20][21][22].The analysis is performed in three centrality classes: 0-10%, 30-50% and 50-90%.The 0-10% and 50-90% centrality classes correspond to the most central and the most peripheral Pb-Pb collisions with small and large impact parameters, respectively.

Resonance reconstruction and yield extraction
The resonances Σ(1385) + and Σ(1385) − are measured separately via their hadronic decay channel to Λ and π ± with a branching ratio of 87% [4], by calculating their invariant mass.The Λ decays weakly into a pπ pair with a branching ratio of 63.9% and a characteristic V 0 topology, which is exploited for its reconstruction.In the following, when referring to Σ(1385) + or Σ(1385) − , the sum of the particle and antiparticle (cc, charge conjugate) is considered.
Primary charged tracks are required to fulfill the standard criteria for good reconstruction quality described in Ref. [23].These selections affect the π ± from the Σ(1385) ± decays, for which it is additionally required to have p T > 120 MeV/c.Pion identification is based on the requirement that the Σ(1385) ± in Pb-Pb ALICE Collaboration dE/dx and the time-of-flight measured in the TPC and in the TOF, respectively, are compatible with their expected values within a fixed number of standard deviations (σ TPC , σ TOF ) [24].In particular, the signals are requested to be within ±5σ TPC , ±3σ TPC and ±2σ TPC for tracks with p T < 0.35 GeV/c, 0.35 < p T < 0.5 GeV/c, and 0.5 < p T < 20 GeV/c, respectively.If the measurement of the time-of-flight is available, this information is used in combination with the TPC PID by applying a selection based on a ±3σ TOF range if p T < 1.5 GeV/c and on a ±2.5σTOF range otherwise.
For the decay products (daughters) of the Λ, a subset of the standard track quality criteria are applied with the addition of p T > 150 MeV/c.Furthermore, the candidate Λs fulfill the criteria of the V 0 decay topology listed in Table 1.The distance of closest approach (DCA) of the Λ daughters is measured in standard deviations of this distribution.The pointing angle is the angle between the direction of the momentum of the Λ and the line connecting the secondary to the primary vertex.The invariant mass window is defined so as to take into account the resolution of the reconstructed Λ and R Λ is the allowed radial distance within which the Λ decay is considered.The Λ daughters are identified as pions and protons based on a 4σ TPC criterium for all p T .|y Λ | < 0.5 DCA of the proton daughter to primary vertex > 0.12 cm DCA of the pion daughter to primary vertex > 0.3 cm DCA between Λ decay products < 0.8 standard deviations Λ cosine of pointing angle > 0.98 The Λπ pairs are reconstructed within the region |η| < 0.8 and |y| < 0.5 by combining candidates from the same event.The uncorrelated combinatorial background is estimated via the mixed-event technique, by combining Λ and pions belonging to different events of the same centrality class with similar primary vertex position along the z-axis and charged-particle multiplicity.Specifically, the difference among the z-coordinate of the vertices of the events being mixed must be less than 1 cm and the difference of the charged-particle multiplicities less than five.To reduce the statistical fluctuations in the mixed-event background distribution, each event is mixed with nine other events.
As the Σ(1385) ± decay products have a large mass difference, a selection on their momenta is applied which helps to reduce the combinatorial background under the peak.This is achieved by requiring the momentum asymmetry of the decay products (defined as (p Λ − p π )/(p Λ + p π )) to be between 0.3 and 0.95.The asymmetry selection results in the partial removal of the correlated background remaining after the event-mixing background subtraction.
The criteria used to select the Λπ pairs are applied to both same-event and mixed-event invariant mass distributions, which are then normalized in the region 1.8−2.0GeV/c 2 .The normalised mixed-event background distributions are subtracted from the same-event ones.The resulting invariant mass (M Λπ ) distributions exhibit the signal peak on top of a residual background consisting of correlated Λπ pairs from other particles or misidentified ones [25].The residual background has a smooth shape that is described by the following function: where m π and m Λ are the π and Λ mass, respectively, taken from Ref. [4]   The mixed-event subtracted invariant mass distribution is fitted with a Breit-Wigner function for the signal and Eq. 1 for the residual background.The width of the resonances is kept fixed at the PDG value of 36 MeV/c 2 for Σ(1385) + and 39 MeV/c 2 for Σ(1385) − [4].This procedure is repeated in nine intervals in the pair p T , from 1 to 9 GeV/c, and for each centrality class.The fit range for each p T interval varies in order to achieve a better χ 2 per degree of freedom.The lower values vary between 1.26 and 1.30 GeV/c 2 and the upper boundaries from 1.55 to 1.70 GeV/c 2 .The yields of Σ(1385) ± are extracted in each p T interval and centrality class by integrating the Breit-Wigner function in the range [M p − 5Γ, M p + 5Γ], where M p is the peak position resulting from the fit and Γ the width of the resonance.In Fig. 1, the mixed-event subtracted invariant mass distribution is presented for Σ(1385) + (left) in the 0-10% centrality class for 3.5 < p T < 4.5 GeV/c and for Σ(1385) − (right) in the 30-50% centrality class for 2.5 < p T < 3.5 GeV/c.In the latter, the Ξ peak is visible at ≈ 1.321 GeV/c 2 , requiring for Σ(1385) − an additional gaussian function to be used in the fit to take into account the Ξ particle.

Corrections
The extracted raw yields of Σ(1385) ± are normalized to the number of events of the corresponding centrality class and corrected for the detector acceptance (A), the reconstruction efficiency (ε rec ) and the branching ratio (BR) [4].The detector acceptance and reconstruction efficiency (called efficiency) are determined from a Monte Carlo simulation based on the HIJING event generator [26].The Σ(1385) ± signals are injected with a flat p T distribution in the range 0−10 GeV/c 2 into HIJING events.The generated particles and their decay products are propagated through the detector material using GEANT3 as transport code [27].The same criteria are applied to select the resonance decay products and the pairs as for the data analysis.The product A × ε rec ×BR is calculated as the fraction of the generated Σ(1385) ± at midrapidity, reconstructed and identified after the application of all selection criteria.The calculation is performed in each centrality class and in Fig. 2, the A × ε rec ×BR (denoted as efficiency×BR) is shown as a function of p T for the three centrality classes.It depends on the centrality due to the centrality dependence of the efficiency of both π and Λ particles.Due to the flat input p T spectrum of the injected Σ(1385) ± resonances, an iterative procedure is applied for re-weighting it, to remove input spectrum shape effects.The efficiencies obtained from the above described procedure are compared with those obtained from standard HIJING simulations and their ratio is compatible with unity.

Systematic uncertainties
The sources of systematic uncertainties in the measurement of the Σ(1385) ± resonances are the yield extraction, the tracking efficiency of the pion daughter, the selection criteria of Λ particles, the pion PID, the knowledge of the detector material budget, the asymmetry selection of the resonance daughter particles and the branching ratio to the decay channel used in the analysis.The uncertainties, except those of the asymmetry and the branching ratio, are p T and centrality dependent.The uncertainty on the raw yield extraction comprises variations of the fit range, of the background fit function and of the mixed-event background normalization range, the use of bin counting instead of the integration of the signal function, as well as letting the width of the Breit-Wigner free in the fit.From all variations, the maximum deviation from the nominal value is assigned as systematic uncertainty.The difference between the global tracking efficiency in the data and the Monte Carlo simulation contributes to the total uncertainty by affecting the pion daughter and it varies from 5% at p T = 1 GeV/c to 3% at p T = 6 GeV/c.The uncertainty due to the Λ topological selection is calculated by varying the DCA of the Λ daughters to the primary vertex, the DCA of the Λ daughters and the cosine of pointing angle.This uncertainty has a value ranging from 2-3% at the lowest p T to 4-5% at the highest p T .The uncertainty associated with the pion identification is quantified by varying the selections in the TPC and the TOF detectors, and reaches a maximum value of 8% for Σ(1385) − in central collisions.The uncertainty on the yields occurring from the implementation of the material budget in the detector simulation was evaluated by increasing and decreasing the material amount within its systematic uncertainty [14], resulting in a variation of ±4.5%.Reconstruction in data and Monte Carlo was done in the two extreme cases and the systematic uncertainty has been estimated on the final observables, which are the p T distributions of the resonances.The asymmetry uncertainty is 5%, regardless of p T and centrality, and is evaluated by varying the accepted asymmetry range.Finally, the branching ratio uncertainty is 1.1% [28].In Table 2, the uncertainties entering the measurement are given for two p T intervals in 0-10% central collisions.The total systematic uncertainty of Σ(1385) − is slightly higher than the uncertainty of Σ(1385) + mainly due to the fact that in the yield extraction of the former an additional gaussian function has to be taken into account in the fit to parametrise the Ξ peak.In the other centrality classes the uncertainties are lower.

Results and Discussion
The transverse-momentum distributions of Σ(1385) + and Σ(1385) − in Pb-Pb collisions are reported in Fig. 3 and compared to the measurements in inelastic pp collisions at √ s = 7 TeV [25] and in non-single diffractive (NSD) p-Pb collisions at √ s NN = 5.02 TeV [28].Each of the Pb-Pb distributions in this figure is fitted individually with a Blast-Wave function [29].The p T -integrated yields and the mean transverse momentum, p T , are calculated based on the data where the spectra are measured and on the extrapolation of the Blast-Wave function at low p T .The fraction of the integrated yields in the high-p T extrapolation region is negligible, while this fraction for the low-p T one ranges from 28% to 37%, depending on the collision centrality.The statistical and systematic uncertainties of the p T -integrated yields and the mean p T are evaluated by repeating the fit after moving the spectra within their statistical and systematic uncertainties, respectively.The procedure of the fitting and the calculation of the p T -integrated yields and the mean p T described above is repeated with Boltzmann, Fermi-Dirac, m T -exponential and Lévy-Tsallis [30] functions to calculate an additional systematic uncertainty defined as the maximum deviation of these quantities from the ones obtained with the Blast-Wave fit.The dN/dy and p T are reported in Tables 3 and 4, respectively, for the three centrality classes of this analysis.As a part of the total systematic uncertainty is correlated across multiplicity, the uncorrelated part of the uncertainties has been calculated by following a similar procedure as the calculation of the total uncertainty.The multiplicity-uncorrelated uncertainties are represented by the third number in Tables 3 and 4. The mean charged-particle multiplicity density for each centrality class, dN ch /dη , is also reported in Table 3.Both dN/dy and p T increase with dN ch /dη in Pb-Pb collisions.
From pp to central heavy-ion collisions, the spectra become harder with increasing charged-particle multiplicity, following a similar trend as observed for the other measured light-flavour particles and resonances [9,24].The p T increases from about 1.15 GeV/c in pp collisions to about 1.37 GeV/c in p-Pb and further increases with centrality in Pb-Pb as reported in Table 4.

Class
p T (GeV/c) Σ(1385) + + cc Σ(1385) − + cc 0-10% (17.185 ± 0.007 ± 0.812 ± 0.748) × 10 −1 (16.974 ± 0.007 ± 1.067 ± 0.974) × 10 −1 30-50% 1.556±0.001±0.068±0.0571.587±0.001±0.082±0.07550-90% 1.465±0.001±0.086±0.0791.551±0.002±0.103±0.070 Figure 4 presents the ratios of the measured spectra of the summed Σ(1385) ± resonance states to the distributions obtained from the EPOS3 and PYTHIA8/Angantyr Monte Carlo event generators in the three Pb-Pb centrality classes.The EPOS3 model [7,8,31] describes the evolution of a heavy-ion collision with the reaction volume being divided into a core and a corona part.For high string densities, the model does not allow the strings to decay independently, instead, if the energy density from string segments is high enough, these fuse into the so-called "core" region, which evolves hydrodynamically.The low energy density region forms the "corona", which hadronizes according to the unmodified string fragmentation.After hadronization, hadrons are fed into the UrQMD hadron cascade afterburner, which describes hadronic interactions in a microscopic approach [6].Previous ALICE measurements of ρ(770) 0 , K * (892) 0 and Λ(1520) in Pb-Pb collisions at √ s NN = 2.76 TeV were compared with predictions from the EPOS 3.107 event generator [5,10,11], indicating that the UrQMD afterburner is necessary for the model to describe the spectral shape of these resonances in central collisions, especially at low p T .Remarkably, the model calculation with UrQMD in [11] predicted that regeneration effects could balance rescattering in the hadronic phase for the Σ(1385) ± decay products, resulting in no suppression for this resonance as a function of centrality.The new Σ(1385) ± data reported in this Letter are compared with the latest version of EPOS3 (EPOS 3.4) both with and without coupling it to the UrQMD afterburner in Fig. 4. No significant difference is observed between the calculation with the UrQMD afterburner and without it in semicentral and peripheral Pb-Pb collisions at √ s NN = 5.02 TeV.A difference in the two predictions is observed for most central collisions for 1 < p T < 5 GeV/c.The model largely overestimates the production of Σ(1385) ± resonances for p T < 5 GeV/c in 0-10% central Pb-Pb collisions, ≈ 60% in the lowest p T interval, whereas it describes the data within 20-30% in semicentral and peripheral collisions.PYTHIA8/Angantyr [32] is an extension of the PYTHIA 8 [33] event generator to the case of heavyion collisions.PYTHIA describes nucleon-nucleon interaction at the parton level, based on multiple partonic interactions and Lund-string hadronization.In Angantyr, PYTHIA is extended to model nucleon-nucleus and nucleus-nucleus collisions in a three steps procedure.First, the number of participating nucleons in a collision is calculated from a Glauber-like model.Second, PYTHIA is used to simulate the collision at the partonic level according to the results of the Glauber calculation as a sum of incoherent inelastic, diffractive and/or elastic collisions.Third, the partonic state is allowed to have final state interactions before it is hadronized according to the Lund-string model in PYTHIA.Hadronic rescattering and regeneration processes like those implemented in UrQMD are not considered in the version of PYTHIA8/Angantyr employed here.The data-to-model ratios reported in Fig. 4 show that PYTHIA8/Angantyr underpredicts the Σ(1385) ± production at low momenta by a factor of 3 to 4 and exhibits a softer spectrum.It only tends to describe the data reasonably well at momenta above p T 7 GeV/c.Σ(1385) ± in Pb-Pb

ALICE Collaboration
The p T -integrated yield ratios of Σ(1385) ± to pions are shown in Fig. 5 for different collision systems as measured by the ALICE [25,28] and the STAR collaborations [34,35].In general, no particular trend with multiplicity is observed given the uncertainties.A fit with a zeroth order polynomial of all data points reported in Fig. 5 yields a χ 2 /NDF of 20.2/10, whereas the exclusion of the Pb-Pb most central point from the fit leads to a χ 2 /NDF of 14.6/9.In addition, for the ratio of the most central Pb-Pb point to the pp one, we obtain 0.86 ± 0.16 by taking into account both the statistical and systematic uncertainties, with the latter having the largest contribution to the error of this ratio.[25,28] together with the STAR measurements [34,35] in various collision systems and energies are reported as a function of the charged-particle multiplicity density at midrapidity, dN ch /dη lab |η lab <0.5| .Data are compared with model calculations for LHC energies from the GSI-Heidelberg grand canonical statistical hadronization model [2], PYTHIA 8.2 [33], PYTHIA8/Angantyr [32] and EPOS3 with and without UrQMD [31].
In the same figure, models applicable to different charged-particle multiplicities are also presented.The comparison to the models, concerning especially the higher multiplicities, can give insight to the dominating mechanism in the hadronic phase that results in the observed behaviour.The statistical hadronization model [2] (indicated as GSI-Heidelberg in Fig. 5) describes the process of hadron formation at the scale where perturbative QCD is no longer applicable.It is assumed that near hadronization the fireball created in heavy-ion collisions is close to thermal equilibrium and hadron yields can be characterized by a grand canonical partition function.The prediction compared to data in Fig. 5 is obtained for a chemical freeze-out temperature of T=156 MeV, which results from a fit to the light-flavour hadron and nucleus yields measured by ALICE at the LHC [2].For the ρ(770) 0 [5], K * (892) 0 [9] and Λ(1520) [10] shortlived resonances, a suppression with respect to the grand canonical statistical hadronization model expectation is observed for central Pb-Pb collisions, while the data in peripheral collisions are well described.This behaviour, which is also observed for Σ(1385) ± in Fig. 5, is typically attributed to rescattering in the hadronic phase.In 0-10% central Pb-Pb collisions, the Σ(1385) ± to pion ratio is 3.6σ lower than Σ(1385) ± in Pb-Pb ALICE Collaboration the statistical hadronization model prediction, whereas the difference of the zeroth order polynomial fit (without taking into account the most central event class) from the same model is a 3.3σ effect.The calculations from EPOS3 with UrQMD, which is able to describe qualitatively the suppression of the K * (892) 0 /K and Λ(1520)/Λ yield ratios from central to peripheral Pb-Pb collisions [5,9,10], reproduce qualitatively also the weak centrality dependence of the Σ(1385) ± /π yield ratio.Its magnitude, however, is overpredicted for all centralities.
It is to be noted that the p T -integrated Σ(1385) ± /π yield ratio (see Fig. 5) in the model exhibits a discrepancy with respect to data that is consistent with the one observed for the Σ(1385) ± spectrum alone.In both cases, EPOS3 with UrQMD overestimates the production.
The centrality dependence of the Σ(1385) ± /π ratio is qualitatively captured also by PYTHIA8/Angantyr, within the current uncertainties.As PYTHIA8/Angantyr underpredicts the production at low transverse momenta of both Σ(1385) ± and pions, it consequently underpredicts the p T -integrated production.Interestingly, this discrepancy is already present in pp collisions [25], indicating that it is not related to the heavy-ion modeling part.The semicentral Pb-Pb point is closer to PYTHIA8/Angantyr prediction because the pion yield is underestimated to a greater degree than the sigma yield.

Conclusions
We have presented the first measurement of Σ(1385) ± production in heavy-ion collisions at LHC energies.An evidence of suppression with respect to the grand canonical thermal model is observed in central collisions as for the K * (892) 0 meson, which has a similar lifetime.While the EPOS3 model coupled to the UrQMD afterburner describes the centrality dependence of the K * (892) 0 data, it clearly overestimates the production of Σ(1385) ± .This may either be caused by a missing element in the model that is common to all centralities or due to the fact that the rescattering and (the seemingly dominant) regeneration effects that EPOS3 withUrQMD predicts for Σ(1385) ± do not manifest themselves in the data.The current implementation of PYTHIA8/Angantyr does not reproduce the Σ(1385) ± yield either, while capturing the centrality dependence of Σ(1385) ± /π within the uncertainties.It remains to be seen if ongoing developments of the model to include hadron rescattering will provide a better agreement with the data.
On the experimental side, future higher precision measurements will clarify if a suppression with respect to pp or peripheral Pb-Pb collisions is present and thus will allow for a model-independent investigation of the rescattering and regeneration picture.In addition, detailed comparisons with future Λ measurements will elucidate if the observed data-model discrepancies are driven by the strangeness content of the hadron under study or the modeling of the hadronic phase.

Figure 1 :
Figure 1: Left:Invariant mass distribution of Λπ pairs for Σ(1385) + after subtraction of the mixed-event background in 0-10% central Pb-Pb collisions and 3.5 < p T < 4.5 GeV/c.Right:Invariant mass distribution of Λπ pairs for Σ(1385) − after subtraction of the mixed-event background in 30-50% centrality class and 2.5 < p T < 3.5 GeV/c.The black curves represent the combined fit using a signal (continuous blue line) plus residual background (red dashed line) model, as described in the text.

Figure 2 :
Figure 2: The product of efficiency×BR for Σ(1385) + (left) and Σ(1385) − (right) in the three centrality classes used for the analysis.The branching ratio, BR, is included in the correction.The error bars represent the statistical uncertainties.

Figure 3 :
Figure 3: p T spectra of Σ(1385) + (left) and Σ(1385) − (right) resonances in Pb-Pb collisions at √ s NN = 5.02 TeV in three centrality classes.Empty and shaded boxes depict the total and uncorrelated uncertainties, respectively.The Blast-Wave fit functions are plotted up to 7 GeV/c.The measurements in pp collisions at √ s = 7 TeV [25] and in p-Pb collisions at √ s = 5.02 TeV (NSD) [28] are quoted for comparisons.

Figure 4 :
Figure 4: Ratio of the measured p T distributions of summed Σ(1385) ± to model predictions from PYTHIA8/Angantyr (green) and EPOS with (blue) and without (red) the UrQMD afterburner.Each panel corresponds to one of the three centrality classes in Pb-Pb collisions at √ s NN = 5.02 TeV.Shaded bands represent the model statistical uncertainty, while bars correspond to the data total uncertainties.

Table 4 :
p T of Σ(1385) ± per centrality class.The first, second, and third uncertainty indicate the statistical, the total systematic uncertainty and the multiplicity-uncorrelated systematic uncertainty, respectively.