In the continuation of our previous work, the transverse momentum (pT ) spectra and nuclear modification factor (RAA) are derived using the relaxation time approximation of Boltzmann Transport Equation (BTE). The initial pT -distribution used to describe p + p collisions has been studied with the pQCD inspired power-law distribution, Hagedorn's empirical formula and with the Tsallis nonextensive statistical distribution. The non-extensive Tsallis distribution is observed to describe the complete range of the transverse momentum spectra. The Boltzmann-Gibbs Blast Wave (BGBW) distribution is used as the equilibrium distribution in the present formalism, to describe the pTdistribution and nuclear modification factor in nucleus-nucleus collisions. The experimental data for Pb+Pb collisions at √ sNN = 2.76 TeV at the Large Hadron Collider at CERN have been analyzed for pions, kaons, protons, K * 0 and φ. It is observed that the present formalism while explaining the transverse momentum spectra up to 5 GeV/c, explains the nuclear modification factor very well up to 8 GeV/c in pT for all these particles except for protons. RAA is found to be independent of the degree of non-extensivity, qpp after pT ∼ 8 GeV/c.