
Using a simple quasiparticle model of QCD matter, presented some time ago in
the literature, in which interactions are modelled by some effective fugacities
z, we investigate the interplay between the dynamical content of fugacities z
and effects induced by nonextensivity in situations when this model is used in
a nonextensive environment characterized by some nonextensive parameter q
$\neq$ 1 (for the usual extensive case q=1). This allows for a better
understanding of the role of nonextensivity in the more complicated
descriptions of dense hadronic and QCD matter recently presented (in which
dynamics is defined by a lagrangian, the form of which is specific to a given
model).

Recently the CERN ALICE experiment, in its dedicated cosmic ray run, observed
muon bundles of very high multiplicities, thereby confirming similar findings
from the LEP era at CERN (in the CosmoLEP project). Significant evidence for
anisotropy of arrival directions of the observed high multiplicity muonic
bundles is found. Estimated directionality suggests their possible
extragalactic provenance. We argue that muonic bundles of highest multiplicity
are produced by strangelets, hypothetical stable lumps of strange quark matter
infiltrating our Universe.

Recently the CERN ALICE experiment, in its dedicated cosmic ray run, observed
muon bundles of very high multiplicities, thereby confirming similar findings
from the LEP era at CERN (in the CosmoLEP project). Originally it was argued
that they apparently stem from the primary cosmic rays with a heavy masses. We
propose an alternative possibility arguing that muonic bundles of highest
multiplicity are produced by strangelets, hypothetical stable lumps of strange
quark matter infiltrating our Universe. We also address the possibility of
additionally deducing their directionality which could be of astrophysical
interest. Significant evidence for anisotropy of arrival directions of the
observed high multiplicity muonic bundles is found. Estimated directionality
suggests their possible extragalactic provenance.

We discuss a Tsallis distribution with complex nonextensivity parameter $q$.
In this case the usual distribution is decorated with a logperiodic
oscillating factor (apparently, such oscillations can bee seen in recently
measured transverse momentum distributions in collisions at very high
energies). Complex $q$ also means complex heat capacity which shall also be
briefly discussed.

It has been recently noticed that transverse momenta ($p_T$) distributions
observed in high energy production processes exhibit remarkably universal
scaling behavior. This is seen when they are in some suitable variable,
replacing the usual $p_T$. On the other hand, it is also known that transverse
momenta distributions in general follow a powerlike Tsallis distribution,
rather than an exponential BoltzmannGibbs, with a (generally energy dependent)
nonextensivity parameter $q$. We now show that it is possible to choose a
suitable variable such that $p_T$ distributions of particles produced in
protonproton interactions in a wide energy range can be fitted by the same
Tsallis distribution (with the same, energy independent, value of the
$q$parameter). Similar scaling behavior in nucleusnucleus collisions is also
observed. The possible dynamical origin of the q parameter used in these fits
will be discussed.

The ALICE Collaboration has measured inclusive J/psi production in pp
collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results
presented in this Letter refer to the rapidity ranges y<0.9 and 2.5<y<4 and
have been obtained by measuring the electron and muon pair decay channels,
respectively. The integrated luminosities for the two channels are L^e_int=1.1
nb^1 and L^mu_int=19.9 nb^1, and the corresponding signal statistics are
N_J/psi^e+e=59 +/ 14 and N_J/psi^mu+mu=1364 +/ 53. We present
dsigma_J/psi/dy for the two rapidity regions under study and, for the forwardy
range, d^2sigma_J/psi/dydp_t in the transverse momentum domain 0<p_t<8 GeV/c.
The results are compared with previously published results at sqrt(s)=7 TeV and
with theoretical calculations.

Energy spectrum of cosmic rays (CR) exhibits powerlike behavior with very
characteristic "knee" structure. We consider a generalized statistical model
for the production process of cosmic rays which accounts for such behavior in a
natural way either by assuming existence of temperature fluctuations in the
source of CR, or by assuming specific temperature distribution of CR sources.
Both possibilities yield the so called Tsallis statistics and lead to the
powerlike distribution. We argue that the "knee" structure arises as result of
abrupt change of fluctuations in the source of CR. Its possible origin is
briefly discussed.

The nonextensive onedimensional version of a hydrodynamical model for
multiparticle production processes is proposed and discussed. It is based on
nonextensive statistics assumed in the form proposed by Tsallis and
characterized by a nonextensivity parameter $q$. In this formulation the
parameter $q$ characterizes some specific form of local equilibrium which is
characteristic for the nonextensive thermodynamics and which replaces the usual
local thermal equilibrium assumption of the usual hydrodynamical models. We
argue that there is correspondence between the perfect nonextensive
hydrodynamics and the usual dissipative hydrodynamics. It leads to simple
expression for dissipative entropy current and allows for predictions for the
ratio of bulk and shear viscosities to entropy density, $\zeta/s$ and $\eta/s$,
to be made.

Examples of joint probability distributions are studied in terms of Tsallis'
nonextensive statistics both for correlated and uncorrelated variables, in
particular it is explicitely shown how correlations in the system can make
Tsallis entropy additive and that the effective nonextensivity parameter $q_N$
decreases towards unity when the number of variables $N$ increases. We
demonstrate that Tsallis distribution of energies of particles in a system
leads in natural way to the Negative Binomial multiplicity distribution in this
system.

We have modelled parton distribution in nuclei using a suitably modified
nuclear Fermi motion. The modifications concern the nucleon rest energy which
changes the Bj\"orken x in nuclear medium. We also introduce final state
interactions between the scattered nucleon and the rest of the nucleus. The
energymomentum sum rule is saturated. Good agreement with experimental data of
the EMC effect for $x> 0.15$ and nuclear lepton pair production data has been
obtained.

Stochastic networks represent very important subject of research because they
have been found in almost all branches of modern science, including also
sociology and economy. We provide a information theory point of view, mostly
based on its nonextensive version, on their most characteristic properties
illustrating it with some examples.

The Interacting Gluon Model (IGM) is a tool designed to study energy flow,
especially stopping and leading particle spectra, in high energy hadronic
collisions. In this model, valence quarks fly through and the gluon clouds of
the hadrons interact strongly both in the soft and in the semihard regime.
Developing this picture we arrive at a simple description of energy loss, given
in terms of few parameters, which accounts for a wide variety of experimental
data. This text is a survey of our main results and predictions.

Recently new data from the CosmoLEP project appeared, this time from DELPHI
detector. They essentially confirm the findings reported some time ago by
ALEPH, namely the appearence of bundles of muons with unexpectedly high
multiplicities, which so far cannot be accounted by present day models. We
argue, using arguments presented by us some time ago, that this phenomenon
could be regarded as one more candidate for the presence in the flux of cosmic
rays entering the Earth's atmosphere from outer space nuggets od Strange Quark
Matter (SQM) in form of so called strangelets.

Assuming that cosmic rays entering the Earth's atmosphere contain a small
admixture of nuggets of strange quark matter in form of strangelets one can
explain a number of apparently "strange" effects observed in different cosmic
rays experiments. We shall demonstrate here that the mass spectrum of such
strangelets filles the "nuclear desert" gap existing between the heaviest
elements observed in Universe and the next "nuclearlike objects" represented
by neutron and strange stars.

We demonstrate that hadron production viewed as formation of specific
stochastic network can explain in natural way the powerlaw distributions of
transverse mass spectra of pions found recently, which seem to substitute the
expected Boltzmann statistical factor.

We have calculated the mass distribution ($d\sigma/dM_X$) as observed by UA8
Collaboration in the inclusive reaction $p \bar p \to p X \bar p$ at $\sqrt s =
630 GeV$, using the Interacting Gluon Model (IGM) with Double Pomeron Exchange
included. The only new parameter is the $IP IP$ cross section, which we can
extract from fitting experimental data. We compare our results with the values
obtained in the UA8 study. Assuming a constant PomeronPomeron total cross
section ($\sigma_{IP IP} = 1 mb$), we make predictions for $d\sigma/dM_X$ at
Tevatron and LHC energies.

We propose simple Monte Carlo method for calculating parton distribution in
nuclei. The EMC effect is automatically included by means of two parameters,
which characterize the change of nuclear pion field. Good agreement with
experimental data in the broad range of variable x is obtained.

We have proposed a new method for the Coulomb wave function correction which
includes the momentum resolution for charged hadron pairs and applied it to the
precise data on pi^+pi^ correlations obtained in p + Ta reaction at 70 GeV/c.
It is found that interaction regions of this reaction (assuming Gaussian source
function) are 9.8+5.8 fm and 7.7+4.8 fm for the thicknesses of the target 8
micro m and 1.4 micro m, respectively. We have also analyzed the data by the
formula including strong interaction effect. The physical picture of the source
size obtained in this way is discussed.

Using a QCD inspired model developed by our group for particle production,
the Interacting Gluon Model (IGM), we have made a systematic analysis of all
available data on leading particle spectra. These data include diffractive
collisions and photoproduction at HERA. With a small number of parameters
(essentially only the nonperturbative gluongluon cross section and the
fraction of diffractive events) good agreement with data is found. We show that
the difference between pion and proton leading spectra is due to their
different gluon distributions. We predict a universality in the diffractive
leading particle spectra in the large momentum region, which turns out to be
independent of the incident energy and of the projectile type.

Recent DESYHERA data on $J/\Psi$ elasticity distribution show that it
emerges mostly as a fast particle. Interpreting photoproduction as a collision
between a preformed charmed hadron and the proton, the outcoming $J/\Psi$ is a
leading particle of the collision. We analyse these data using a model
formulated to describe energy flow in hadronhadron reactions. The measured
$J/\Psi$ spectrum can be successfully described in terms of this model. We
conclude that the observed transparency of the charmed hadronproton collisions
arises because of the particularly small gluonic content of the initial $c 
\bar c$ state.

We have applied the Interacting Gluon Model (IGM) to calculate diffractive
mass spectra measured recently in ep collisions at HERA. We show that it is
possible to treat them in terms of gluongluon collisions in the same way as
was done before for hadronic collisions. Analysis of available data is
performed.

We have extended the Interacting Gluon Model (IGM) to calculate diffractive
mass spectra generated in hadronic collisions. We show that it is possible to
treat both diffractive and nondiffractive events on the same footing, in terms
of gluongluon collisions. A systematic analysis of available data is
performed. The energy dependence of diffractive mass spectra is addressed. They
show a moderate narrowing at increasing energies. Predictions for LHC energies
are presented.

The requirements for observing Centaurolike phenomena and their role as
possible probes of deeply penetrating component in Cosmic Rays are discussed.

We discuss leading charm production in connection with energy deposition in
the central rapidity region. Special attention is given to the correlation
between production in central and fragmentation regions. If the fraction of the
reaction energy released in the central region increases the asymmetry in the
$x_F$ distributions of charmed mesons will become smaller. We illustrate this
quantitatively with simple calculations performed using the Interacting Gluon
Model. Leading beauty production is also considered.

We present improved Coulomb correction formulae for BoseEinstein
correlations including also exchange term and use them to calculate appropriate
correction factors for several source functions. It is found that Coulomb
correction to the exchange function in the BoseEinstein correlations cannot be
neglected.