Measurements of inclusive spectra and mean multiplicities of $\pi^\pm$, K$^\pm$, p and $\bar{\textrm{p}}$ produced in inelastic p+p interactions at incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7, 8.8, 12.3 and 17.3 GeV, respectively) were performed at the CERN Super Proton Synchrotron using the large acceptance NA61/SHINE hadron spectrometer. Spectra are presented as function of rapidity and transverse momentum and are compared to predictions of current models. The measurements serve as the baseline in the NA61/SHINE study of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter.

We present measurements of $\rho^0$, $\omega$ and K$^{*0}$ spectra in $\pi^{-} + $C production interactions at 158 GeV/c and $\rho^0$ spectra at 350 GeV/c using the NA61/SHINE spectrometer at the CERN SPS. Spectra are presented as a function of the Feynman's variable $x_\text{F}$ in the range $0 < x_\text{F} < 1$ and $0 < x_\text{F} < 0.5$ for 158 GeV/c and 350 GeV/c respectively. Furthermore, we show comparisons with previous measurements and predictions of several hadronic interaction models. These measurements are essential for a better understanding of hadronic shower development and for improving the modeling of cosmic ray air showers.

Results on two-particle $\Delta\eta\Delta\phi$ correlations in inelastic p+p interactions at 20, 31, 40, 80, and 158~GeV/c are presented. The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The data show structures which can be attributed mainly to effects of resonance decays, momentum conservation, and quantum statistics. The results are compared with the EPOS and UrQMD models.

Measurements of particle emission from a replica of the T2K 90 cm-long carbon target were performed in the NA61/SHINE experiment at CERN SPS, using data collected during a high-statistics run in 2009. An efficient use of the long-target measurements for neutrino flux predictions in T2K requires dedicated reconstruction and analysis techniques. Fully-corrected differential yields of $\pi^\pm$-mesons from the surface of the T2K replica target for incoming 31 GeV/c protons are presented. A possible strategy to implement these results into the T2K neutrino beam predictions is discussed and the propagation of the uncertainties of these results to the final neutrino flux is performed.

Measurements of multiplicity and transverse momentum fluctuations of charged particles were performed in inelastic p+p interactions at 20, 31, 40, 80 and 158 GeV/c beam momentum. Results for the scaled variance of the multiplicity distribution and for three strongly intensive measures of multiplicity and transverse momentum fluctuations \$\Delta[P_{T},N]\$, \$\Sigma[P_{T},N]\$ and \$\Phi_{p_T}\$ are presented. For the first time the results on fluctuations are fully corrected for experimental biases. The results on multiplicity and transverse momentum fluctuations significantly deviate from expectations for the independent particle production. They also depend on charges of selected hadrons. The string-resonance Monte Carlo models EPOS and UrQMD do not describe the data. The scaled variance of multiplicity fluctuations is significantly higher in inelastic p+p interactions than in central Pb+Pb collisions measured by NA49 at the same energy per nucleon. This is in qualitative disagreement with the predictions of the Wounded Nucleon Model. Within the statistical framework the enhanced multiplicity fluctuations in inelastic p+p interactions can be interpreted as due to event-by-event fluctuations of the fireball energy and/or volume.

Production of $d$, $t$, and $^3$He nuclei in central Pb+Pb interactions was studied at five collision energies ($\sqrt{s_{NN}}=$ 6.3, 7.6, 8.8, 12.3, and 17.3 GeV) with the NA49 detector at the CERN SPS. Transverse momentum spectra, rapidity distributions, and particle ratios were measured. Yields are compared to predictions of statistical models. Phase-space distributions of light nuclei are discussed and compared to those of protons in the context of a coalescence approach. The coalescence parameters $B_2$ and $B_3$, as well as coalescence radii for $d$ and $^3$He were determined as a function of transverse mass at all energies.

Inclusive production of $\Lambda$-hyperons was measured with the large acceptance NA61/SHINE spectrometer at the CERN SPS in inelastic p+p interactions at beam momentum of 158~\GeVc. Spectra of transverse momentum and transverse mass as well as distributions of rapidity and x$_{_F}$ are presented. The mean multiplicity was estimated to be $0.120\,\pm0.006\;(stat.)\,\pm 0.010\;(sys.)$. The results are compared with previous measurements and predictions of the EPOS, UrQMD and FRITIOF models.

Measurements of hadron production in p+C interactions at 31 GeV/c are performed using the NA61/ SHINE spectrometer at the CERN SPS. The analysis is based on the full set of data collected in 2009 using a graphite target with a thickness of 4% of a nuclear interaction length. Inelastic and production cross sections as well as spectra of $\pi^\pm$, $K^\pm$, p, $K^0_S$ and $\Lambda$ are measured with high precision. These measurements are essential for improved calculations of the initial neutrino fluxes in the T2K long-baseline neutrino oscillation experiment in Japan. A comparison of the NA61/SHINE measurements with predictions of several hadroproduction models is presented.

We look for fluctuations expected for the QCD critical point using an intermittency analysis in the transverse momentum phase space of protons produced around midrapidity in the 12.5\% most central C+C, Si+Si and Pb+Pb collisions at the maximum SPS energy of 158$A$~GeV. We find evidence of power-law fluctuations for the Si+Si data. The fitted power-law exponent $\phi_{2} = 0.96^{+0.38}_{-0.25}\text{ (stat.)}$ $\pm 0.16\text{ (syst.)}$ is consistent with the va\-lue expected for critical fluctuations. Power-law fluctuations had previously also been observed in low-mass $\pi^+ \pi^-$ pairs in the same Si+Si collisions.

Results from the NA49 experiment at the CERN SPS are presented on event-by-event transverse momentum and multiplicity fluctuations of charged particles, produced at forward rapidities in central Pb+Pb interactions at beam momenta 20$A$, 30$A$, 40$A$, 80$A$, and 158$A$ GeV/c, as well as in systems of different size ($p+p$, C+C, Si+Si, and Pb+Pb) at 158$A$ GeV/c. This publication extends the previous NA49 measurements of the strongly intensive measure $\Phi_{p_T}$ by a study of the recently proposed strongly intensive measures of fluctuations $\Delta[P_T, N]$ and $\Sigma[P_T, N]$. In the explored kinematic region transverse momentum and multiplicity fluctuations show no significant energy dependence in the SPS energy range. However, a remarkable system size dependence is observed for both $\Delta[P_T, N]$ and $\Sigma[P_T, N]$, with the largest values measured in peripheral Pb+Pb interactions. The results are compared with NA61/SHINE measurements in $p+p$ collisions, as well as with predictions of the UrQMD and EPOS models.

We present experimental results on inclusive spectra and mean multiplicities of negatively charged pions produced in inelastic p+p interactions at incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7, 8.8, 12.3 and 17.3 GeV, respectively). The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. Two-dimensional spectra are determined in terms of rapidity and transverse momentum. Their properties such as the width of rapidity distributions and the inverse slope parameter of transverse mass spectra are extracted and their collision energy dependences are presented. The results on inelastic p+p interactions are compared with the corresponding data on central Pb+Pb collisions measured by the NA49 experiment at the CERN SPS. The results presented in this paper are part of the NA61/SHINE ion program devoted to the study of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter. They are required for interpretation of results on nucleus-nucleus and proton-nucleus collisions.

A novel approach, the identity method, was used for particle identification and the study of fluctuations of particle yield ratios in Pb+Pb collisions at the CERN Super Proton Synchrotron (SPS). This procedure allows to unfold the moments of the unknown multiplicity distributions of protons (p), kaons (K), pions ($\pi$) and electrons (e). Using these moments the excitation function of the fluctuation measure $\nu_{\text{\text{dyn}}}$[A,B] was measured, with A and B denoting different particle types. The obtained energy dependence of $\nu_{\text{dyn}}$ agrees with previously published NA49 results on the related measure $\sigma_{\text{dyn}}$. Moreover, $\nu_{\text{dyn}}$ was found to depend on the phase space coverage for [K,p] and [K,$\pi$] pairs. This feature most likely explains the reported differences between measurements of NA49 and those of STAR in central Au+Au collisions.

The T2K long-baseline neutrino oscillation experiment in Japan needs precise predictions of the initial neutrino flux. The highest precision can be reached based on detailed measurements of hadron emission from the same target as used by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. The corresponding data were recorded in 2007-2010 by the NA61/SHINE experiment at the CERN SPS using a replica of the T2K graphite target. In this paper details of the experiment, data taking, data analysis method and results from the 2007 pilot run are presented. Furthermore, the application of the NA61/SHINE measurements to the predictions of the T2K initial neutrino flux is described and discussed.

Measurements of charged pion and kaon production are presented in centrality selected Pb+Pb collisions at 40A GeV and 158A GeV beam energy as well as in semi-central C+C and Si+Si interactions at 40A GeV. Transverse mass spectra, rapidity spectra and total yields are determined as a function of centrality. The system-size and centrality dependence of relative strangeness production in nucleus-nucleus collisions at 40A GeV and 158A GeV beam energy are derived from the data presented here and published data for C+C and Si+Si collisions at 158A GeV beam energy. At both energies a steep increase with centrality is observed for small systems followed by a weak rise or even saturation for higher centralities. This behavior is compared to calculations using transport models (UrQMD and HSD), a percolation model and the core-corona approach.

Spectra of positively charged kaons in p+C interactions at 31 GeV/c were measured with the NA61/SHINE spectrometer at the CERN SPS. The analysis is based on the full set of data collected in 2007 with a graphite target with a thickness of 4% of a nuclear interaction length. Interaction cross sections and charged pion spectra were already measured using the same set of data. These new measurements in combination with the published ones are required to improve predictions of the neutrino flux for the T2K long baseline neutrino oscillation experiment in Japan. In particular, the knowledge of kaon production is crucial for precisely predicting the intrinsic electron neutrino component and the high energy tail of the T2K beam. The results are presented as a function of laboratory momentum in 2 intervals of the laboratory polar angle covering the range from 20 up to 240 mrad. The kaon spectra are compared with predictions of several hadron production models. Using the published pion results and the new kaon data, the K+/\pi+ ratios are computed.

Kaons and protons carry large parts of two conserved quantities, strangeness and baryon number. It is argued that their correlation and thus also fluctuations are sensitive to conditions prevailing at the anticipated parton-hadron phase boundary. Fluctuations of the $(\mathrm{K}^+ + \mathrm{K}^-)/(\mathrm{p}+\bar{\mathrm{p}})$ and $\mathrm{K}^+/\mathrm{p}$ ratios have been measured for the first time by NA49 in central Pb+Pb collisions at 5 SPS energies between $\sqrt{s_{NN}}$= 6.3 GeV and 17.3 GeV. Both ratios exhibit a change of sign in $\sigma_{\mathrm{dyn}}$, a measure of non-statistical fluctuations, around $\sqrt{s_{NN}}$ = 8 GeV. Below this energy, $\sigma_{\mathrm{dyn}}$ is positive, indicating higher fluctuation compared to a mixed event background sample, while for higher energies, $\sigma_{\mathrm{dyn}}$ is negative, indicating correlated emission of kaons and protons. The results are compared to UrQMD calculations which which give a good description at the higher SPS energies, but fail to reproduce the transition to positive values.

The momentum correlation between protons and lambda particles emitted from central Pb+Pb collisions at sqrt(s_{NN}) = 17.3 GeV was studied by the NA49 experiment at the CERN SPS. A clear enhancement is observed for small relative momenta (q_{inv} < 0.2 GeV). By fitting a theoretical model, which uses the strong interaction between the proton and the lambda in a given pair, to the measured data a value for the effective source size is deduced. Assuming a static Gaussian source distribution we derive an effective radius parameter of R_G = 3.02 \pm 0.20$(stat.)^{+0.44}_{-0.16}(syst.) fm.

The yields of (anti-)protons were measured by the NA49 Collaboration in centrality selected Pb+Pb collisions at 40A GeV and 158A GeV. Particle identification was obtained in the laboratory momentum range from 5 to 63 GeV/c by the measurement of the energy loss dE/dx in the TPC detector gas. The corresponding rapidity coverage extends 1.6 units from mid-rapidity into the forward hemisphere. Transverse mass spectra, the rapidity dependences of the average transverse mass, and rapidity density distributions were studied as a function of collision centrality. The values of the average transverse mass as well as the midrapidity yields of protons when normalized to the number of wounded nucleons show only modest centrality dependences. In contrast, the shape of the rapidity distribution changes significantly with collision centrality, especially at 40A GeV. The experimental results are compared to calculations of the HSD and UrQMD transport models.