  In this article, we have studied the nuclear modification factor measured in
Pb-Pb collisions ($R_{PbPb}$) for $\pi^{\pm}$, $K^{\pm}$, $p+\bar{p}$, $K^{*0}
+ \bar{K^{*0}}$, $\phi$ and in p-Pb collisions ($R_{pPb}$) for $\pi^{\pm}$,
$K^{\pm}$, $p+\bar{p}$ at Large hadron collider (LHC) energy of $\sqrt{s_{NN}}$
= 5.02 TeV for the most central and peripheral collisions. We have also
analysed the experimental data of transverse momentum ($p_T$) spectra for these
identified hadrons at LHC for Pb-Pb as well as for p-Pb collisions. We have
used Boltzmann transport equation (BTE) in relaxation time approximation (RTA)
for this analysis. The Tsallis statistics is used as an initial distribution
function and The Tsallis blast wave (TBW) model is employed as an equilibrium
distribution in BTE. The present model fits the measured transverse momentum
spectra, $R_{PbPb}$, and $R_{pPb}$ successfully upto $p_T$ = 8 GeV with a
reasonable $\chi^2/ndf$ for all the considered hadrons at various centralities.
The experimental data for $R_{pPb}$ are generated using the particle yields at
pPb and pp collisions where number of binary collisions are taken from Glauber
model calculations. We find that the average transverse flow velocity
($<\beta_r>$) follows the mass and centrality ordering and decreases with the
mass as well as when one move from the central collisions to peripheral
collisions. These findings are inline with the results of the hydrodynamical
calculations.
