K∗(892)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^{*}(892)^0$$\end{document} meson production in inelastic p+p interactions at 158 GeV/c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{ Ge }\text{ V }\!/\!c$$\end{document} beam momentum measured by NA61/SHINEat the CERN SPS

The measurement of K∗(892)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^{*}(892)^0$$\end{document} resonance production via its K+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^{+}\pi ^{-}$$\end{document} decay mode in inelastic p+p collisions at beam momentum 158 GeV/c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{ Ge }\text{ V }\!/\!c$$\end{document} (sNN=17.3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}}=17.3$$\end{document} GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{ Ge }\text{ V }$$\end{document}) is presented. The data were recorded by the NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The template method was used to extract the K∗(892)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^{*}(892)^0$$\end{document} signal and double-differential transverse momentum and rapidity spectra were obtained. The full phase-space mean multiplicity of K∗(892)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^{*}(892)^0$$\end{document} mesons was found to be (78.44±0.38(stat)±6.0(sys))·10-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(78.44 \pm 0.38 \mathrm {(stat)} \pm 6.0 \mathrm {(sys)) \cdot 10^{-3}}$$\end{document}. The NA61/SHINEresults are compared with the Epos1.99 and Hadron Resonance Gas models as well as with world data from p+p and nucleus–nucleus collisions.


Introduction and motivation
Strange hadron production is believed to be an important tool to study the dynamics of high-energy collisions. In collisions achieving high energy densities strangeness production was predicted to be enhanced [1] as a result of the decrease of the mass of strangeness carriers due to partial chiral symmetry restoration. The K * (892) 0 resonance state contains ans valence quark and is therefore sensitive to the level of strangeness production. Thus, the data on K * (892) 0 meson production provide a more complete understanding of hadron chemistry.
Measurements of the production of short-lived resonances are a unique tool to understand the less known aspects of high energy collisions, especially their time evolution. In heavy ion collisions the yields of resonances may help to distinguish between two possible freeze-out scenarios: the sudden and the gradual one [2]. Namely, the ratio of K * (892) 0 to charged kaon production may allow to estimate the time interval between chemical (end of inelastic collisions) and kinetic (end of elastic collisions) freeze-out. The lifetime of the K * (892) 0 resonance (≈ 4 fm/c) is comparable to the expected duration of the rescattering hadronic gas phase between the two freeze-out stages. Consequently, a certain fraction of K * (892) 0 resonances will decay inside the fireball. The momenta of their decay products are expected to be significantly modified by elastic scatterings, preventing the experimental reconstruction of the resonance via an invariant a e-mail: kperl@cern.ch mass analysis. In such a case a suppression of the observed K * (892) 0 yield is expected. Such an effect was indeed observed in nucleus-nucleus collisions at Super Proton Synchrotron (SPS) and Relativistic Heavy Ion Collider (RHIC) energies [3][4][5][6][7][8]. The ratio of K * /K production (K * stands for K * (892) 0 , K * (892) 0 or K * ± , and K denotes K + or K − ) showed a decrease with increasing system size as expected due to the increasing rescattering time between chemical and kinetic freeze-out. The same effect was recently reported also by the ALICE Collaboration at the Large Hadron Collider (LHC) [9][10][11][12].
When looking at the energy dependence of the K * 0 /K − ratio 1 in central Pb+Pb or Au+Au collisions, a bit larger suppression of K * 0 is observed for the 2.76 TeV 2 LHC energy [10] when compared to the top RHIC ( √ s N N = 200 GeV) energy [7], namely K * 0 /K − = 0.180± 0.027 (0.186 ± 0.027) for the 0-5% (5-10%) central Pb+Pb reactions at LHC and 0.20 ±0.04 for the 0-10% most central Au+Au interactions at RHIC. Those values can be compared with those for p+p interactions, which are 0.307 ± 0.043 at LHC [10] and 0.34 ± 0.05 at RHIC [7]. Thus, the K * 0 /K − ratio in central Pb+Pb collisions at LHC (2.76 TeV) drops to 59 (61)% of the value found for p+p interactions. For RHIC energies this drop is similar and equals 59%.
In the NA49 experiment at the CERN SPS K * (892) 0 and K * (892) 0 meson production was analyzed separately and the corresponding (almost 4π ) mean multiplicities obtained in the 23.5% most central Pb+Pb collisions at √ s N N = Thus, at SPS energy the resonance to non-resonance ratio in central Pb+Pb drops to about 43-44% of the value for p+p interactions. This effect is even stronger than the one observed at RHIC and LHC and might suggest that the lifetime (calculated in the K * 0 rest frame; see Eq. (21) in Sect. 5.4) of the hadron gas system created in central nucleusnucleus collisions at the SPS is longer than that at higher energies. Eventually, resonance regeneration processes start to play a role for higher energies counteracting the K * 0 suppression due to rescattering. It should also be pointed out that the whole picture assumes that the conditions at chemical freeze-out of p+p and Pb+Pb collisions are the same. More detailed calculations of the time between freeze-outs, both in the K * 0 rest frame and in collision center-of-mass reference system, are given in Sect. 5.4. The results for p+p collisions provide an important baseline for heavier nucleus-nucleus systems. So far the K * 0 /K − ratio for p+p interactions did not show large differences between the top RHIC and four LHC energies [10,12,15,16]. Most of the results at lower energies are less reliable due to large uncertainties, see the compilation in Ref. [15], and new points in Refs. [4,10,12,16]. This emphasizes the need to obtain high precision p+p data at energies lower than the top RHIC energy. Continuing considerations for p+p collisions, a very intriguing effect was reported in the most recent ALICE analysis of the multiplicity dependence in p+p collisions [17,18]. The K * 0 / K ± and K * 0 /K 0 S ratios decrease when going from low-multiplicity to high-multiplicity p+p interactions at the LHC energies. This may be an indication of a hadronic phase with significant non-zero lifetime even in p+p collisions.
The transverse mass spectra and yields of K * (892) 0 mesons are also important inputs for Blast-Wave models (determining kinetic freeze-out temperature and transverse flow velocity) and Hadron Resonance Gas models (determining chemical freeze-out temperature, baryochemical potential, strangeness under-saturation factor, system volume, etc.). Those models significantly contribute to our understanding of the phase diagram of strongly interacting matter. In principle, the precise determination of transverse flow velocity is attractive due to the fact, that recent LHC, RHIC and even SPS results suggest that dense and collectively behaving system may appear also in collisions of small nuclei, or even in elementary interactions. Finally, the study of resonances in elementary interactions contributes to the understanding of hadron production, due to the fact that products of resonance decays represent a large fraction of the final state particles. Resonance spectra and yields provide an important reference for tuning Monte Carlo string-hadronic models.
This paper reports measurements of K * (892) 0 resonance production via its K + π − decay mode in inelastic p+p collisions at beam momentum of 158 GeV/c ( √ s N N = 17.3 GeV) 3 [24]. The data were recorded by the NA61/SHINE hadron spectrometer [25] at the CERN SPS. Unlike in the previous NA49 analysis [4] at the same beam momentum, the template method was used to extract the K * (892) 0 signal. This method was found to allow a more precise background subtraction than the standard procedure. Moreover, the large statistics NA61/SHINE data (about 52.5M events recorded with the interaction trigger compared to 2.5M p+p events analysed in NA49 [26,27]) allowed to obtain high quality double-differential transverse momentum and rapidity spectra of K * (892) 0 mesons. The paper is organized as follows. Section 2 briefly describes the NA61/SHINE detector. Section 3 discusses the analysis procedures, including event and track cuts, method of signal extraction, corrections, and evaluation of uncertainties. The final results are presented in Sect. 4 and their comparison with world data and models in Sect. 5. A summary Sect. 6 closes the paper.

Experimental setup
The NA61/SHINE experiment [25] uses a large acceptance hadron spectrometer located in the CERN North Area. The schematic layout of the NA61/SHINE detector is shown in Fig. 1. The detailed description of the full detector can be found in Ref. [25]. Here only the detector components, which were used in this analysis, are described.
A set of scintillation and Cherenkov counters as well as beam position detectors (BPDs) upstream of the spectrometer provide timing reference, identification and position measurements of incoming beam particles. The trigger scintillator counter S4 placed downstream of the target is used to select events with collisions in the target area by the absence of a charged particle hit.
Secondary beams of positively charged hadrons at 158 GeV/c are produced from 400 GeV/c protons extracted from the SPS accelerator. Particles of the secondary hadron beam are identified by two Cherenkov counters, a CEDAR [30] (for 158 GeV/c beam CEDAR-N) and a threshold counter (THC). The CEDAR counter, using a coincidence of six Fig. 1 The schematic layout of the NA61/SHINE experiment at the CERN SPS (horizontal cut, not to scale). The beam and trigger detector configuration used for data taking in 2009 is shown in the inset (see Refs. [28,29] for detailed description). The chosen coordinate system is drawn on the lower left: its origin lies in the middle of the VTPC-2, on the beam axis   out of the eight photo-multipliers placed radially along the Cherenkov ring, provides positive identification of protons, while the THC, operated at pressure lower than the proton threshold, is used in anti-coincidence in the trigger logic. A selection based on signals from the Cherenkov counters allowed one to identify beam protons with a purity of about 99%. A consistent value for the purity was found by bending the beam into the TPCs with the full magnetic field and using identification based on its specific ionization energy loss dE/dx [31].
The main tracking devices of the spectrometer are four large volume Time Projection Chambers (TPCs). Two of them, the vertex TPCs (VTPC-1 and VTPC-2), are located in the magnetic fields of two super-conducting dipole magnets with a combined bending power of 9 Tm which corresponds to about 1.5 T and 1.1 T fields in the upstream and downstream magnets, respectively.
Two large main TPCs (MTPC-L and MTPC-R) are positioned downstream of the magnets symmetrically to the beam line. The fifth small TPC (GAP TPC) is placed between VTPC-1 and VTPC-2 directly on the beam line. It closes the gap between the beam axis and the sensitive volumes of the other TPCs. The TPCs are filled with Ar and CO 2 gas mixtures. Particle identification in the TPCs is based on measurements of the specific energy loss (dE/dx) in the chamber gas.
The p+p data sets, which are the topic of this paper, were recorded with the proton beam incident on a liquid hydrogen target (LHT), a 20 cm long cylinder positioned about 80 cm upstream of VTPC-1.

Data sets
The results for p+p interactions are based on high-statistics data runs (in years 2009, 2010, and 2011) which recorded about 56.65×10 6 collisions (52.53M×10 6 events selected by   Figure 2 (right) shows that for the 2009 production the ratio of the number of events in the target-removed sample to the number of events in the target-inserted sample (ratio calculated in the range − 590 < z < − 572 cm; histograms normalized in the range − 450 < z < − 300 cm) is on the level of 4.8%, and therefore no correction for non-target interactions was applied. An alternative method of analysis (see for example Ref. [28]) would be to measure and subtract the resonance yields in the target-removed data, but both the standard method and the template-fitting method used in this paper cannot be applied to data sets with small statistics such as the target removed data recorded by NA61/SHINE . In order to estimate the systematic biases related to the contamination by off-target interactions the window of z-position of the main vertex was varied (see Sect. 3.10). Table 1 presents the details of data sets collected in the three separate data taking periods. The number of events recorded with the interaction trigger, as well as the number of events selected for the analysis (see Sect. 3.3) are shown. One sees that only 44-56% of the events were used for the analysis. This drop is caused mainly by BPD reconstruction inefficiencies and off-target interactions accepted by the trigger. The number of tracks, given in the Table 1, refers to tracks registered in accepted events only. The agreement of the fractions of accepted tracks in the three analyzed data sets confirms the similarity of the data recorded in 2009, 2010 and 2011. For the analysis of K * (892) 0 production these three data sets were combined at the level of preparing invariant mass distributions (Sect. 3.6).

Analysis method
The details of NA61/SHINE calibration, track and vertex reconstruction procedures, as well as simulations used to correct the reconstructed data, are discussed in Refs. [28,29,32]. In the following section the analysis technique developed for the measurement of the K * (892) 0 spectra in p+p interactions is described. The procedure used for the data analysis consists of the following steps: (i) application of event and track selection criteria, (ii) selection of K + and π − candidates based on the measurement of their ionization energy loss (dE/dx) in the gas volume of the TPCs, (iii) creation of invariant mass distribution of K + π − pairs, (iv) creation of invariant mass distribution of K + π − pairs for mixed events and Monte Carlo templates, (v) extraction of K * (892) 0 signal, (vi) application of corrections (obtained from simulations) to the raw numbers of K * (892) 0 ; they include losses of inelastic p+p interactions due to the on-line and offline event selection as well as losses of K * (892) 0 due to track and pair selection cuts and the detector geometrical acceptance.
The details of the steps are described in the following subsections.

Event selection
Inelastic p+p interactions were selected by the following criteria: (i) an interaction was recognized by the trigger logic (see Refs. [28,29] for detailed description), (ii) no off-time beam particle was detected within ± 1 μs around the trigger (beam) particle, (iii) the trajectory of the beam particle was measured in at least one of BPD-1 or BPD-2 and in the BPD-3 detector and was well reconstructed, (iv) the primary interaction vertex fit converged, (v) the fit of the z-coordinate of the primary p+p interaction vertex (see Fig. 2) converged and the fitted z position was found between −590 cm and −572 cm, where the center of the LHT was at −580 cm. The range of this cut was selected to maximize the number of good events and minimize the contamination by off-target interactions, (vi) events with a single, well measured positively charged track with absolute momentum close to the beam momentum ( p > p beam − 1 GeV/c) were rejected.
The above event cuts select well measured inelastic p+p interactions. The background due to elastic interactions is removed (cuts (iv) and (vi)). The contribution of off-target interactions is reduced (cut (v)). The losses of inelastic interactions due to the event selection procedure were corrected using simulations (see below). The number of events after these cuts is 27.85 × 10 6 .

Track selection
After the event selection criteria a set of track quality cuts were applied to individual tracks. These cuts were used to ensure high reconstruction efficiency, proper identification of tracks and to reduce the contamination of tracks from secondary interactions, weak decays and off-time interactions. The individual tracks were selected by the following criteria: (i) the track fit including the interaction vertex converged, (ii) the total number of reconstructed points on the track should be greater than 30, (iii) the sum of the number of reconstructed points in VTPC-1 and VTPC-2 was greater than 15 or the number of reconstructed points in the GAP TPC was greater than 4, (iv) the distance between the track extrapolated to the interaction plane and the interaction point (impact parameter) should be smaller than 4 cm in the horizontal (bending) plane and 2 cm in the vertical (drift) plane, (v) the track momentum (in the laboratory reference system) is in the range 3 ≤ p lab ≤ 158 GeV/c, (vi) the track transverse momentum is required to be smaller than 1.5 GeV/c, (vii) dE/dx track cuts were applied to select K + and π − candidates (see Sect. 3.5).
The number of tracks left after these cuts is about 49.79 × 10 6 .

Selection of kaon and pion candidates
Charged particle identification in the NA61/SHINE experiment is based on the measurement of their ionization energy loss (dE/dx) in the gas of the TPCs and of the time of flight (tof) obtained from the ToF-L and ToF-R walls. For the region of the relativistic rise of the ionization at large momenta, the measurement of dE/dx alone allows identification. At lower momenta the dE/dx bands for different particle species overlap and the identification based only on measurements of dE/dx in the TPCs (this analysis) is not enough. For this reason the track cut (v) was applied. In Fig. 3 the dE/dx values as a function of total momentum ( p lab ), measured in the laboratory reference system, are shown for positively and negatively charged particles, separately. The K + and π − candidates were selected by requiring their dE/dx values to be within 1.5σ or 3.0σ around their nominal Bethe-Bloch values, respectively. Here σ represents the typical standard deviation of a Gaussian fitted to the dE/dx distribution of kaons and pions. Since only small variations of σ were observed for different total momentum and transverse momentum bins, fixed values σ = 0.044 were used for K + and σ = 0.052 for π − . The bands of selected K + and π − candidates are shown in the bottom panel of Fig. 3.
The raw numbers of K * (892) 0 are usually obtained by performing fits to the invariant mass spectra with the sum of a background and a signal function. The invariant mass is defined as: where E represents the total energy and p the momentum vector of daughter particles from K * (892) 0 decay. In the standard method (mixing method) the large combinatorial background is estimated by invariant mass spectra calculated for K + π − pairs originating from different events. Figures 4 and 5 (top, left) show combinatorial background histograms (red points) compared to the data histograms of m K + π − (blue points). Mixed events were normalized to the same number of pairs as in real data in the invariant mass range from 0.6 to 1.6 GeV. After subtracting the normalized mixed event background the blue points in Figs. 4, 5 (bottom, left) were obtained. The K * (892) 0 signal is prominently seen, but the histogram still shows a residual background, seen especially for low invariant mass values. This residual background probably comes from the products of other resonance decays, which are not properly accounted for by the event-mixing, and should be subtracted. The final fit (total fit 2) was performed with the function of Eq. (2) using an where d, e, f , and g are free parameters of the fit, and the Breit-Wigner (BW ) component is described by Eq. (3): where A is the normalization factor, and m K * and K * are also fitted. The initial values of the mass (m K * ) and width ( K * ) parameters of K * (892) 0 were taken from the Particle Data Group (PDG): m K * = m 0 = 0.89555 GeV and K * = 0 = 0.0473 GeV [35]. The red lines (polynomial background) in Figs. 4, 5 (bottom, left) show the fitted additional background component (Eq. (2) without BW ) and the brown lines (total fit 2) the total fit result (Eq. (2)).
In order to obtain a better background description compared to the mixing method, the template method was applied. Namely, the invariant mass spectra of the data (blue data points in Figs. 4, 5 (top, right)) were fitted with a function given by Eq. (4): is the background estimated based on the mixing method, which was discussed above. The T MC res template (MC stands for Monte Carlo) is the shape of background, which describes the contribution of K + π − pairs originating from: (i) combination of tracks that come from decays of resonances different than K * (892) 0 , for example one track from a ρ 0 meson and one from a K * + meson, (ii) combination of tracks where one comes from decay of a resonance and one comes from direct production in the primary interaction.
The T MC res templates were constructed by passing p+p interactions, generated with the Epos1.99 [33] hadronic interaction model using the CRMC 1.4 package [34], through the NA61/SHINE detector simulation chain and then through the same reconstruction routines as the data. The simulation keeps the history of particle production thus allowing to identify their identity and origin enabling the construction the proper templates. For the reconstructed MC events, the same event and track selection criteria, as for real data, were used. They also include the effects of the limited acceptance of the detector. Both the template and the data histograms were computed in bins of rapidity y (calculated in the center-ofmass reference system) and transverse momentum p T .
Finally, the signal (BW ) is described using the Breit-Wigner distribution Eq. For each m K + π − invariant mass bin in Fig. 4 or 5 (bottom, right), the bin content N bin (m K + π − ) was calculated as: where N raw (m K + π − ) is the raw production in a given m K + π − bin, and a, b, T MC res (m K + π − ) and T D AT A mi x (m K + π − ) are described in Eq. (4). The statistical uncertainty of N bin (m K + π − ) can be expressed as (the notation (m K + π − ) is omitted for simplifying the formula): where histograms the number of entries had to be properly normalized. Due to high statistics of data, Monte Carlo, and mixed events, the uncertainties of parameters a and b were neglected. Finally, for each bin of (y, p T ) in Fig. 6 the uncorrected number of K * (892) 0 , N K * (y, p T ), was calculated as the integral over the BW signal of total fit 2 in Figs. 4, 5 (bottom, right). The integral was obtained within the mass window m 0 ± 4 0 . The statistical uncertainty of the raw number of K * (892) 0 , N K * (y, p T ), was taken as the uncertainty of the integral calculated by the ROOT [36] package using covariance matrix of the fitted parameters.

Correction factors
In order to determine the number of K * (892) 0 mesons produced in inelastic p+p interactions, two corrections were applied to the extracted raw number of K * (892) 0 : (i) The loss of the K * (892) 0 due to the dE/dx requirement was corrected by a constant factor: where K + = 0.866 and π − = 0.997 are the probabilities (based on the cumulative Gaussian distribution) for K + or π − to lie within 1.5σ or 3σ around the nominal Bethe-Bloch value. (ii) A detailed Monte Carlo simulation was performed to correct for geometrical acceptance, reconstruction efficiency, losses due to the trigger bias, detector acceptance as well as the quality cuts applied in the analysis. The width of the K * (892) 0 resonance was simulated according to the known PDG value [37]. The correction factors are based on 227.9 × 10 6 inelastic p+p events produced by the Epos1.99 event generator [33]. The validity of these events for calculation of the corrections was verified in Refs. [28,38]. The particles in the generated events were tracked through the NA61/SHINE apparatus using the Geant3 package [39]. The TPC response was simulated by dedicated NA61/SHINE software packages which take into account all known detector effects. The simulated events were reconstructed with the same software as used for real events and the same selection cuts were applied (except the identification cuts: dE/dx and total momentum p lab ).
For each y and p T bin, the correction factor c MC (y, p T ) was calculated as: where: The reconstructed charged particles were matched to the simulated K + and π − based on cluster positions. Then the invariant mass was calculated for all K + π − pairs. The reconstructed number of K * (892) 0 was obtained by repeating the same steps (template method) as in raw data; they are described in Sect. 3.6, -N gen events is the number of generated inelastic p+p interactions (227.9 × 10 6 ), -N sel events is the number of accepted p+p events (140.1× 10 6 ).
The uncertainty of c MC (y, p T ) was calculated assuming that the denominator n sel (y, p T ) is a subset of the nominator n gen (y, p T ) and thus has a binomial distribution. The uncertainty of c MC (y, p T ) was calculated as follows: The values of correction factors c MC , together with statistical uncertainties, are presented in Fig. 7 for all analyzed (y, p T ) bins.
3.9 Corrected K * (892) 0 yields The double-differential yield of K * (892) 0 per inelastic event in a bin of (y, p T ) is calculated as follows: where: -B R = 2/3 is the branching ratio of K * (892) 0 decay into K + π − pairs (obtained [27] from the Clebsch-Gordan coefficients), -N K * (y, p T ) is the uncorrected number of K * (892) 0 , obtained by the signal extraction procedure described in Sect. 3.6, -N events is the number of events after cuts, c d E/dx , c MC (y, p T ) are the correction factors described above, y and p T are the bin widths.
The corrected double-differential yields of K * (892) 0 together with their uncertainties are presented in Sect. 4.

Statistical and systematic uncertainties
The statistical uncertainties of the corrected double-differential yields (see Eq. (10)) take into account the statistical uncertainties of c MC (y, p T ) (see Eq. (9)) and the statistical uncertainties N K * (y, p T ) (see Sect. 3.7) of the uncorrected number of K * (892) 0 . The correction c d E/dx has no statistical uncertainty. The final formula is expressed as follows: The systematic uncertainties were estimated taking into account two sources. The first group of effects is associated with the signal extraction procedure and the second with event and track quality cuts.
The considered sources of the systematic uncertainty and the corresponding modifications of the analysis method were the following: (I) The uncertainty due to the signal extraction procedure: (i) the lower limit of the invariant mass fitting range (see Figs. 4, 5 (top, right)) was changed from 0.66 GeV to 0.69 GeV, (ii) the initial value of the K * parameter of the signal function was changed by ± 8%, (iii) the initial value of the mass parameter of the Breit-Wigner distribution was changed by ± 0.3 MeV, (iv) the parameters a, b and c describing the contribution of the templates in the fitting function (see Eq. (4)) were changed by ± 10%, (v) the value of the K * parameter of the signal function was fixed at the PDG value 0 , (vi) the value of the m K * parameter of the signal function was fixed at the PDG value m 0 , (vii) in the final step of the background fit (see Figs. 4,5 (bottom, right)) the standard polynomial curve of the 2nd order was changed into a polynomial curve of the 3rd order, (viii) the invariant mass range over which the raw number of K * (892) 0 was integrated was changed from m 0 ± 4 0 to ± 3.5 0 and ± 4.5 0 , (ix) the raw number of K * (892) 0 was calculated as the sum of points (after 2nd order polynomial subtraction) instead of the BW signal integral.
(II) The effects of event and track quality cuts were checked by performing the analysis with the following cuts changed compared to the original values: (i) the window in which off-time beam particles are not allowed was increased from ± 1 µs to ± 1.5 µs around the trigger particle, (ii) the cut on the z-position of the interaction vertex was changed from [− 590; −572] cm to [− 591; − 571] cm and [− 589, − 573] cm, (iii) the standard dE/dx cuts (± 3σ for π − and ± 1.5σ for K + ) were modified to ± 2.5σ for π − , ± 1.0σ for K + (narrower cut) and ± 3.5σ for π − , ± 2.0σ for K + (wider cut), (iv) the minimum required total number of points in all TPCs for K * (892) 0 decay products was changed from 30 to 25 and 35, (v) the minimum required number of clusters in both VTPCs for K * (892) 0 decay products was changed from 15 to 12 and 18, (vi) the impact parameter cuts for the tracks were turned off.
For each of the possible sources described above the partial systematic uncertainty σ i was calculated as half of the difference between the lowest and the highest value obtained by varying the given parameter. Then, the final systematic uncertainty was taken as: σ sys = σ 2 i . The contributions of uncertainties σ i to the total uncertainty are negligible for I (ii), I (iii), and I (iv). The final systematic uncertainties are shown in the figures as light red shaded bands.

Mass and width of K * (892) 0
The values of mass and width of K * (892) 0 mesons were extracted from the fits to background subtracted invariant mass spectra (see Sect. 3.6). They are presented in Fig. 8 in different transverse momentum bins (numerical data are listed in Table 2). The results are shown for the rapidity range 0 < y < 0.5. Within uncertainties, the K * values are consistent with information provided by the PDG. However, one observes a slight increase of the m K * parameter with p T with an average close to the PDG value. The corresponding slope is significant since a large part of the shown systematic uncertainty is due to the magnetic field uncertainty (see below). The points (with their statistical uncertainties), presented in the left panel of Fig. 8, were fitted with a linear function resulting in the slope parameter value equal to 4.5 ± 1.2. The change of the m K * parameter with transverse momentum does not introduce a systematic variation of the K * (892) 0 yield since the parameter is fitted in each (y, p T ) bin, and the signal integration range (≈ 380 MeV) is much larger than the m K * change (≈ 6 MeV). The magnetic field strength was verified with a precision of better than 1% by studying the K 0 S and invariant mass distributions [40]. In order to check how the magnetic field calibration influences the results, the momentum components of K * (892) 0 decay products (kaons and pions) were varied by ± 1%. Such a change did not affect K * (892) 0 width and yield significantly. However, the resulting changes of the mass parameter are equal or larger than uncertainties described in Sect. 3.10, and they were taken into account in the calculation of the final uncertainty of the K * (892) 0 mass parameter shown in Fig. 8 (left) and Table 2.
The comparison of mass and width of K * (892) 0 mesons with other experiments is shown in Sect. 5. Table 2 Numerical values of mass and width of K * (892) 0 mesons fitted in 0 < y < 0.5 and presented in Fig. 8. The first uncertainty is statistical, while the second one is systematic 4.2 Double-differential K * (892) 0 spectra The double-differential yields of K * (892) 0 mesons in inelastic p+p interaction at 158 GeV/c in bins of (y, p T ) are presented in Fig. 9. The numerical values with statistical and systematic uncertainties are presented in Table 3. Figure 10 shows the double-differential yields of K * (892) 0 mesons as function of p T presented for separate rapidity bins. The corresponding numerical values are listed in Table 3.

Transverse momentum and transverse mass spectra
In order to measure the inverse slope parameter T of transverse momentum spectra and to estimate the yield of K * (892) 0 mesons in the unmeasured high p T region, the function: Fig. 9 Double-differential K * (892) 0 spectra in inelastic p+p interaction at 158 GeV/c in bins of (y, p T ) as obtained from Eq. (10). The numerical values are given in Table 3 was fitted to the measurements shown in Fig. 10. The inverse slope parameters obtained from the fits are cited in the figure legends.
The transverse mass (m T ≡ p 2 T + m 2 0 ) spectra 1 m T d 2 n dm T dy were calculated based on d 2 n dydp T spectra according to: The results are shown in Fig. 11 and the numerical values are presented in Table 4. For the mid-rapidity region (0 < y < 0.5) the inverse slope parameter of the transverse momentum spectrum was found to be equal to T = (173±3±9) MeV, where statistical uncertainty (the first one) is equal to the uncertainty of the fit parameter, and the systematic uncertainty was estimated in the way described in Sect. 3.10. The NA49 experiment measured the T parameter of the p T spectrum in the rapidity range 0.2 < y < 0.7 and reported a value T = (166 ± 11 ± 10) MeV [4].

p T -integrated and extrapolated rapidity distribution
The rapidity distribution dn dy was calculated by integrating and extrapolating (for the non-measured highp T region) the d 2 n dy dp T spectrum: dn dy = i d 2 n dy dp T · dp T + A p T I p T i d 2 n dy dp T · dp T , (14) Table 3 Numerical values of double-differential yields d 2 n dy dp T presented in Fig. 10, given in units of 10 −3 (GeV/c) −1 . The first uncertainty is statistical, while the second one is systematic   dy dp T for five bins of rapidity. The fitted function (solid line) is given by Eq. (12). The numerical values are listed in Table 3 and the fitted inverse slope parameters T for each bin are given in the legends where: T ⎞ ⎠ dp T , The parameters T were taken from the corresponding plots in Fig. 10. The statistical uncertainties of p T -integrated and     Table 4. The solid lines represent function given by Eqs. (12) and (13) with A and T parameters taken from Fig. 10 extrapolated dn dy values were calculated as follows: The p T -integrated and extrapolated dn dy spectrum of K * (892) 0 mesons is plotted in Fig. 12 and the numerical values are listed in Table 5.
A Gaussian function: dm T dy given in units of 10 −3 (GeV) −2 and presented in Fig. 11; the values of m T − m 0 specify the bin centers. The first uncertainty is statistical, while the second one is systematic  Fig. 12 The p T -integrated and extrapolated rapidity distribution. The fitted Gaussian function (solid line) is given by Eq. (17); the first point (with y < 0) was not included in the fit (see the text for details). The numerical data are listed in Table 5 was fitted to the data points to measure the width σ y of the K * (892) 0 rapidity distribution. The first point with y < 0 is plotted only to check the symmetry of the distribution and was not included in the fit. The fit was also used to determine the mean multiplicity K * (892) 0 (see Sect. 4.5 for details of the procedure). The statistical uncertainty of σ y was taken from the fit and the systematic uncertainty was estimated in the way described in Sect. 3.10. The numerical values of σ y and K * (892) 0 are listed in Table 5. Table 5 Numerical values of the p T -integrated and extrapolated dn dy distribution presented in Fig. 12. The first uncertainty is statistical, while the second one is systematic. Additionally, the width of the Gaussian fit to the dn dy distribution, as well as the mean multiplicity of K * (892) 0 mesons are shown (see the text for details) The mean multiplicity of K * (892) 0 mesons was calculated as the sum of measured points in Fig. 12 (the first point, with y < 0, was not included in the sum) and the integral of the fitted Gaussian function Eq. (17) in the unmeasured region assuming symmetry around y = 0: Fig. 13 The transverse momentum dependence of mass and width of K * (892) 0 (or K * 0 ) mesons obtained by NA61/SHINE , NA49 [4], ALICE [9] and STAR [5]. For ALICE and STAR the averaged (K * 0 ) measurements of K * (892) 0 and K * (892) 0 are shown. The horizontal lines represent PDG values [35] where: The statistical uncertainty of K * (892) 0 was obtained from the formula: and the systematic uncertainty was estimated in the way described in Sect. 3.10. The mean multiplicity of K * (892) 0 mesons, produced in inelastic p+p collisions at 158 GeV/c, is equal to (78.44±0.38 ± 6.0) ·10 −3 , where the first uncertainty is statistical and the second one is systematic.

Comparison with world data and model predictions
This section compares the NA61/SHINE measurements in inelastic p+p interactions at 158 GeV/c with publicly available world data as well as with predictions from microscopic and statistical models. 5.2 Comparison of rapidity spectra and yields with NA49 and Epos1.99 The NA61/SHINE measurements of the rapidity spectrum and mean multiplicity were also compared to those predicted by the model of hadron production Epos1.99 [33]. The results are presented in Fig. 14 and the numerical values of the multiplicity are listed in Table 6. One sees that the Epos1.99 model overestimates K * (892) 0 production in inelastic p+p collisions at 158 GeV/c. Table 6 also shows the comparison with the NA49 result [4] for the same collision system and beam momentum. Instead of analysing in separate p T bins, as in NA61/SHINE , the NA49 experiment used one wide p T bin (0 < p T < 1.5 GeV/c). The mean multiplicity of K * (892) 0 in NA49 was obtained as the integral under the Gaussian function in the range − 3 < y < 3 fitted to the dn dy distribution [27]. Within the uncertainties shown, the results of both experiments are consistent.

System size dependence of K * (892) 0 at 158A GeV/c and predictions of HGM
The statistical Hadron Resonance Gas Models (HGM) are commonly used to predict particle multiplicities in elementary and nucleus-nucleus collisions, using as adjustable parameters the chemical freeze-out temperature T chem , the baryochemical potential μ B , strangeness saturation parameter γ S , etc. In the following the measured K * (892) 0 multiplicities are compared with predictions of two HGM models described in Refs. [41,42]. In Ref. [41] the HGM results for K * (892) 0 multiplicities were calculated for two versions of the model fits to particle yields. The first one, called fit B, allowed for strangeness under-saturation so the usual parametrization with γ S was applied. For p+p interactions, the fit was carried out without including the multiplicities of and baryons. In the second fit, called A, the parameter γ S was replaced by the mean number of strange quark pairs ss . For p+p collisions fit A was performed without the φ meson. For both fits predicted multiplicities were calculated in the Canonical Ensemble (CE) [41]. The measured mean multiplicity of K * (892) 0 in 158 GeV/c inelastic p+p interactions was divided by HGM predictions based on fit A and B and compared with the value found by NA49 [4]. The results are shown in Fig. 15 for p+p interactions, as well as C+C, Si+Si, and Pb+Pb collisions measured by NA49 [4]. In Ref. [41] the S-Canonical Ensemble (SCE) with exact strangeness conservation and grandcanonical treatment of electric charge and baryon number was used for the heavier C+C and Si+Si systems, and the Grand Canonical Ensemble (GCE) was assumed for Pb+Pb collisions. For C+C and Si+Si interactions all available particles were used in the HGM fits, including φ meson and multi-strange baryons. For Pb+Pb data only the measured (1520) yield was removed from the fitted multiplicities. Note that the centrality of Pb+Pb collisions used in the HGM fits was 0-5% whereas the K * (892) 0 values in NA49 were obtained for the 0-23.5% most central interactions. Therefore, the HGM yields had to be scaled by a factor 262/362 Table 6 The mean multiplicities K * (892) 0 and the widths of the rapidity distributions σ y obtained from dn dy distributions (see the text for details). The first uncertainty is statistical and the second systematic  [42,43]. N W denotes the number of wounded nucleons taken from Ref. [4] corresponding to the respective number of wounded nucleons (see Table 8).
For heavier systems (including C+C and Si+Si), there is no significant difference between fit A and fit B, however the deviation between the HGM predictions and experimental data increases with increasing system size. The p+p measurements are very close to the HGM prediction but only in case of fit A, where the φ meson was excluded from the fit. In the most recent paper [42], where the HGM fits were done for the NA49 and the new NA61/SHINE measurement in p+p interactions, it is also stressed that at SPS energies the φ meson multiplicities in p+p collisions cannot be well fitted within the CE formulation of the HGM (the quality of CE fits becomes much worse when the φ meson yield is included). However, the mean multiplicity of K * (892) 0 mesons in inelastic p+p collisions at 158 GeV/c can also be compared to the HGM prediction based on the Grand Canonical Ensemble formulation [42]. The results for the NA49 and NA61/SHINE measurements are shown in Fig. 15 as closed cross and closed star symbols. Surprisingly, the GCE statistical model provides a good description of the K * (892) 0 yield in the small p+p system. The numerical values of the NA61/SHINE p+p measurement and the statistical models are presented in the Table 7. In Fig. 15 the total uncertainty of K * (892) 0 was taken as the square root of the sum of squares of statistical and systematic uncertainties. The uncer- tainty of the ratio shown on vertical axis was taken as the final uncertainty of K * (892) 0 divided by K * (892) 0 HGM .

K * over charged kaon ratios and time between freeze-outs
The K * to charged kaons ratios may allow to estimate the time interval between chemical and kinetic freeze-out in nucleus-nucleus collisions. The K * mesons have identical quark (anti-quark) content as K mesons, but different mass and relative orientation of quark spins. Thus, the K * (892) 0 / K − and K * (892) 0 / K + ratios are considered as the least model dependent ratios for studying the K * production properties as well as the freeze-out conditions. The system size dependence of the K * /K ratio at SPS, RHIC and LHC energies shows a strong decrease with increasing system size and/or multiplicity density (see Sect. 1 for a full list of references). The effect seems to be stronger at the SPS than at RHIC and LHC. Figure 16 presents this dependence at the SPS for the NA49 and NA61/SHINE results at 158A GeV/c. The numerical values are given in Table 8.
The NA61/SHINE K * (892) 0 / K + and K * (892) 0 / K − yield ratios for p+p interactions and the corresponding ratios in central Pb+Pb collisions from NA49 can be used to estimate the time interval between chemical and kinetic freeze-outs in Pb+Pb. Following Ref. [5]: where: -the ratio K * (892) 0 / K +/− in inelastic p+p interactions can be treated as the one at chemical freeze-out, -the ratio K * (892) 0 / K +/− for central Pb+Pb (NA49) interactions can be used as the one at kinetic freeze-out,  Fig. 16 The system size dependences of K * (892) 0 / K + and K * (892) 0 / K − yield ratios in p+p, C+C, Si+Si and Pb+Pb collisions at 158A GeV/c. N W denotes the number of wounded nucleons taken from Ref. [4]. The numerical values are listed in Table 8. For better visibility the NA61/SHINE points are shifted on the horizontal axis τ is the mean K * (892) 0 lifetime of approximately 4.17 fm/c [35], t is the time interval between chemical and kinetic freeze-outs calculated in the K * (892) 0 rest frame.
Assuming that the losses of K * (892) 0 before kinetic freeze-out are due to rescattering effects and that there are no regeneration processes, the time between chemical and kinetic freeze-outs (in the resonance rest frame) can be estimated as 3.7 ± 1.2 fm/c from the K * (892) 0 / K + ratio and 3.2 ± 1.2 fm/c from the K * (892) 0 / K − ratio. These numbers correspond to 23.5% of the most central Pb+Pb interactions but the time is similar when using 5% of the most central events.
Following Ref. [12], the above times may be expressed in the collision center-of-mass reference system using the multiplicative Lorentz factor: where p T can be used as an approximation for K * (892) 0 total momentum for the measurements at mid-rapidity. The NA49 experiment published the K * (892) 0 transverse momentum spectrum for 23.5% of the most central Pb+Pb interactions in the rapidity range 0.43 < y < 1.78 [4]. The p T can be obtained from the fitted exponential function in the range 0 < p T < 4 GeV/c. The average transverse momentum of K * (892) 0 mesons was found to be 0.908 GeV/c that results in γ ≈ 1.42. Finally, the Lorentz boosted time interval Table 8 The mean multiplicities of different particle species measured in nucleus-nucleus collisions at 158 A GeV/c by NA49 and NA61/SHINE . The total uncertainties of K * (892) 0 , K + and K − were taken as the square roots of the sums of squares of statistical and systematic uncertainties. For NA49 p+p data, the K + and K − results include statistical uncertainties only ( K + = 0.2267 ± 0.0006 and K − = 0.1303 ± 0.0004), whereas systematic uncertainties for total yields were not reported [45]. Therefore, NA61/SHINE K + and K − values were used in the K * (892) 0 / K + and K * (892) 0 / K − ratios. The numbers of K + and K − and their uncertainties for the 5% most central Pb+Pb collisions were multiplied by a factor 262/362 in order to estimate charged kaon multiplicities in the 23.5% most central Pb+Pb reactions 0.0784 ± 0.0060 0.234 ± 0.022 [14] 0.132 ± 0.014 [14] 0 --between chemical and kinetic freeze-outs can be estimated as 5.3 fm/c for the K * (892) 0 / K + ratio or 4.6 fm/c for the K * (892) 0 / K − ratio.
Similar calculations can be performed for the published RHIC ( √ s N N = 200 GeV) and LHC ( √ s N N = 2760 GeV) data. In the STAR experiment at RHIC, the K * 0 /K − ratio was found to be 0.20 ± 0.04 for the 10% most central Au+Au collisions, and 0.34 ± 0.05 for p+p interactions [7]. Thus, the time between freeze-outs (calculated in the K * 0 rest frame) is equal to 2.2 ± 1.0 fm/c. The average transverse momentum of K * 0 mesons in Au+Au collisions at mid-rapidity (|y| < 0.5) was found to be 1.09 GeV/c [7], which corresponds to γ ≈ 1.57. Therefore, the time interval between freeze-outs, determined in the collision center-of-mass reference system, can be estimated as 3.5 fm/c. This value is smaller than the ones obtained at SPS. In the ALICE experiment at LHC, the K * 0 /K − ratio was found to be 0.180 ± 0.027 for the 5% most central Pb+Pb collisions, and 0.307 ± 0.043 for p+p interactions [10]. Following Eq. (21), t can be evaluated as 2.2 ± 0.9 fm/c. The p T of K * 0 mesons in Pb+Pb collisions at mid-rapidity (|y| < 0.5) was found to be 1.310 GeV/c [10], which corresponds to γ ≈ 1.77. Thus, at LHC energy, the time interval between freeze-outs, determined in the collision center-of-mass reference system, can be estimated as 3.9 fm/c.
The above numbers may imply that, in central heavy ion collisions, the lifetime of the hadronic period of the fireball after chemical freeze-out is longer at SPS than at RHIC or even at LHC energies. One should, however, remember that such a conclusion is valid only under the assumption that there are no regeneration processes of K * 0 mesons before kinetic freeze-out. As the K * (892) 0 regeneration may happen at all energies, the obtained time interval values should be considered as lower limits of the time between chemical and kinetic freeze-outs.

Summary
In this paper the NA61/SHINE measurement of K * (892) 0 meson production via its K + π − decay mode in inelastic p+p collisions at beam momentum 158 GeV/c ( √ s N N = 17.3 GeV) was presented. The template method was used to extract raw K * (892) 0 signals. In this method the background is described as a sum of two contributions: background due to uncorrelated pairs modeled by event mixing and background of correlated pairs modeled by Epos1.99. For K * (892) 0 production the template method was found to provide a better background description than the standard one which relies on mixed events only. The mass and width of the K * (892) 0 were extracted from the fits to background subtracted invariant mass spectra. Their values, for different transverse momentum bins, are close to the PDG results, however, a slight increase of the K * (892) 0 mass with transverse momentum can be observed.
With the large statistics of NA61/SHINE data (52.53M events selected by the interaction trigger) it was possible to obtain double-differential transverse momentum and rapidity spectra of K * (892) 0 mesons. The full phase-space mean multiplicity of K * (892) 0 mesons, obtained from the p Tintegrated and extrapolated rapidity distribution, was found to be (78.44±0.38±6.0)·10 −3 , where the first uncertainty is statistical and the second one is systematic. The result agrees with the previous NA49 measurement for the same system and energy.
The NA61/SHINE result was compared with predictions of statistical Hadron Resonance Gas models in Canonical and Grand Canonical formulations. Surprisingly, the GCE model provides a good description of the NA61/SHINE measurement of the K * (892) 0 multiplicity in p+p collisions. The CE model also agrees provided that the φ meson is excluded from the fits.
Finally, based on the previous results of NA49 from central Pb+Pb collisions and the new measurements of NA61/SHINE on p+p interactions, an attempt was made to estimate the time between chemical and kinetic freeze-outs in central Pb+Pb reactions at 158 A GeV/c. This time was found to be larger than at RHIC, suggesting that either the system life-time between freeze-outs is indeed higher at SPS or the K * (892) 0 regeneration effects start to play a significant role at higher collision energies.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data (final numerical results) obtained during this study are included in this published article. There are no external data associated with the manuscript.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .