
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 stringresonance 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 eventbyevent fluctuations of the fireball energy
and/or volume.

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.
Twodimensional 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 nucleusnucleus and protonnucleus collisions.

Spectra of K0S mesons and Lambda hyperons were measured in p+C interactions
at 31 GeV/c with the large acceptance NA61/SHINE spectrometer at the CERN SPS.
The data were collected with an isotropic graphite target with a thickness of
4% of a nuclear interaction length. Interaction cross sections, charged pion
spectra, and charged kaon spectra were previously measured using the same data
set. Results on K0S and Lambda production in p+C interactions serve as
reference for the understanding of the enhancement of strangeness production in
nucleusnucleus collisions. Moreover, they provide important input for the
improvement of neutrino flux predictions for the T2K long baseline neutrino
oscillation experiment in Japan. Inclusive production cross sections for K0S
and Lambda are presented as a function of laboratory momentum in intervals of
the laboratory polar angle covering the range from 0 up to 240 mrad. The
results are compared with predictions of several hadron production models. The
K0S mean multiplicity in production processes <n_K0S> and the inclusive cross
section for K0S production were measured and amount to 0.127 + 0.005 (stat) +
0.022 (sys) and 29.0 + 1.6 (stat) + 5.0 (sys) mb, respectively.

The T2K longbaseline 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 20072010 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.

By direct MonteCarlo simulations it is shown that the anisotropic flows can
be successfully described in the colour string picture with fusion and
percolation provided anisotropy of particle emission from the fused string is
taken into account. Quenching of produced particles in the strong colour field
of the string is the basic mechanism for this anisotropy. The concrete
realization of this mechanism is borrowed from the QED. Due to dependence of
this mechanism on the external field strength the found flows grow with energy,
with values for $v_2$ at LHC energies greater by ~15% than at RHIC energies.

In the framework of the classical Glauber approach, the analytical
expressions for the variance of the number of wounded nucleons and binary
collisions in AA interactions at a given centrality are presented. Along with
the optical approximation term, they contain additional contact terms arising
only in the case of nucleusnucleus collisions. The magnitude of the additional
contributions, e.g., for PbPb collisions at SPS energies, is larger than the
contribution of the optical approximation at some values of the impact
parameter. The sum of the additional contributions is in good agreement with
the results of independent Monte Carlo simulations of this process. Due to
these additional terms, the variance of the total number of participants for
peripheral PbPb collisions and the variance of the number of collisions at all
values of the impact parameter exceed several multiples of the Poisson
variances. The correlator between the numbers of participants in colliding
nuclei at fixed centrality is also analytically calculated.

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.

Interaction cross sections and charged pion spectra in p+C interactions at 31
GeV/c were measured with the large acceptance NA61/SHINE spectrometer at the
CERN SPS. These data are required to improve predictions of the neutrino flux
for the T2K long baseline neutrino oscillation experiment in Japan. A set of
data collected during the first NA61/SHINE run in 2007 with an isotropic
graphite target with a thickness of 4% of a nuclear interaction length was used
for the analysis. The measured p+C inelastic and production cross sections are
257.2 + 1.9 + 8.9 mb and 229.3 + 1.9 + 9.0 mb, respectively. Inclusive
production cross sections for negatively and positively charged pions are
presented as a function of laboratory momentum in 10 intervals of the
laboratory polar angle covering the range from 0 up to 420 mrad. The spectra
are compared with predictions of several hadron production models.

The correlation between multiplicities in two separated rapidity windows, the
socalled longrange correlation (LRC), is studied in the framework of the
model with independent identical emitters. It's shown that the LRC coefficient,
defined for the scaled (relative) variables, nevertheless depends on the
absolute width of the forward rapidity window and does not depend on the width
of the backward one. The dependence of the LRC coefficient on the forward
rapidity acceptance is explicitly found with only one theoretical parameter.
The preliminary comparison with ALICE 7TeV pp collisions data shows that the
multiplicity LRC in the data can be described in the framework of the suggested
approach.

It is shown that naive two stage scenario of the soft multiparticle
production in hadronic and nuclear collisions at high energy, when at first
stage the colour strings are formed and at the second stage these strings, or
some other (higher colour) strings formed due to fusion of primary strings, are
decaying, emitting observed particles, encounters some difficulties at the
attempt to analyse the spacetime picture of the process. Simple analysis shows
the dominant is the process when the formation and the decay of a string occur
in parallel  a string breaks into two parts already at rather small length
(about 12 fm in its c.m. system), then the process repeats in the pieces and
so on. Nevertheless it is proved to be possible to agree the string fusion idea
with the spacetime picture of a string decay. In the framework of the
ArtruMennessier model of a string fragmentation the simple interpretation of
the homogeneity of the rapidity distribution for hadrons produced from the
decay of a single string at high energy is presented and the analytical
estimate for the density of this rapidity distribution is obtained.

In the framework of the classical Glauber approach the exact analytical
expression for the variance of the number of participants (wounded nucleons)
for given centrality AA interactions is presented. It's shown, that in the case
of nucleusnucleus collisions along with the optical approximation term the
additional "contact" term appears. The numerical calculations for PbPb
collisions at SPS energies show that at intermediate values of the impact
parameter the "optical" and "contact" terms contributions to the variance of
the participants number are of the same order and their sum is in a good
agreement with the results of independent MC simulations of this process. The
correlation between the numbers of participants in colliding nuclei is taken
into account. In particular it's demonstrated that in PbPb collisions at SPS
energies the variance of the total number of participants approximately three
times exceeds the Poisson one in the impact parameter region 1012 Fm. The
fluctuations of the number of collisions are also discussed.

The theoretical description of the correlations between observables in two
separated rapidity intervals in relativistic nuclear collisions is presented.
It is shown, that the eventbyevent ptpt correlation defined as the
correlation between event mean values of transverse momenta of all particles
emitted in two different rapidity intervals does not decrease to zero with the
increase of the number of strings in contrast with two particle ptpt
correlation  the correlation between the transverse momenta of single
particles produced in these two rapidity windows.
In the idealized case with the homogeneous string distribution in the
transverse plane in the framework of the cellular analog of string fusion model
(SFM) the asymptotic of ptpt correlation coefficient is analytically
calculated and compared with the results of the MonteCarlo (MC) calculations
fulfilled both in the framework of the original SFM and in the framework of its
cellular analog, which enables to control the MC algorithms.
In the case with the realistic nucleon distribution density of colliding
nuclei the results of the MC calculations of the ptpt correlation function for
minimum bias nuclear collisions at SPS, RHIC and LHC energies are presented and
analysed.

The theoretical description of the correlations between observables in two
separated rapidity intervals for AAinteractions at high energies is presented.
In the case with the realistic nucleon distribution density of colliding nuclei
the MC calculations of the longrange correlation functions at different values
of impact parameter are done.
For nn and ptn correlations it is shown that the impact parameter
fluctuations at a level of a few fermi, unavoidable in any experiment,
significantly change the magnitude of correlation coefficients. The rise of
ptn and especially ptpt correlation coefficients is found when one passes
from SPS to RHIC and LHC energies.

The longrange $p_t$ and multiplicity($n$) correlations in highenergy
nuclear collisions are studied in the framework of a simple cellular analog of
the string fusion model.
Two cases with local and global string fusion is considered. The $p_t$$n$
and $n$$n$ correlation functions and correlation coefficients are calculated
analytically in some asymptotic cases using suggested Gauss approximation.
It's shown that at large string density the $p_t$$n$ and $n$$n$
correlation coefficients are connected and the scaling takes place. The
behavior of the correlations at small string density is also studied.
The asymptotic results are compared with results of the numerical
calculations in the framework of proposed cellular approach.

A simple cellular model for the description of the longrange multiplicity
and $p_t$ correlations in highenergy nuclear collisions originating from the
string fusion model is proposed. Three versions of the model: without fusion,
with local and with global string fusion are formulated.
A Gauss approximation which enables explicit analytical calculations of the
correlation functions in some asymptotic cases in the framework of the model is
developed. The assumptions of the model and the validity of a Gauss
approximation are checked up in the simplest (no fusion) case when the explicit
solution of the model can be found.
The role of a size of cells is analysed. The modification of the results in
the case of nonPoissonian distributions is also discussed.