  The pseudorapidity distributions of charged particles measured in p+p($\rm
\overline{p}$) collisions for energies ranging from $\sqrt{s_{NN}}=23.6$ GeV to
13 TeV, d+Au collisions at $\sqrt{s_{NN}}=200$ GeV, p+Pb collisions at
$\sqrt{s_{NN}}= 5.02$ TeV and A+A collisions at RHIC and LHC are investigated
in the fireball model with Tsallis thermodynamics. We assume that the rapidity
axis is populated with fireballs following q-Gaussian distribution and the
charged particles follow the Tsallis distribution in the fireball. The
theoretical results are in good agreement with the experimental data for all
the collision systems and centralities investigated. The collision energy and
centrality dependence of the central position $y_0$ and its width $\sigma$ of
the fireball distribution are also investigated. A possible application of the
model to predict the charged particle pseudorapidity distributions for the
system size scan program proposed recently for the STAR experiment at RHIC is
proposed.

Nuclear Theory