  The pseudorapidity distribution of charged hadron over a wide $\eta$ range
gives us crucial information about the dynamics of particle production.
Constraint on the detector acceptance, particularly at forward rapidities,
demands a proper distribution function to extrapolate the pseudorapidity
distribution to large $\eta$. In this work, we have proposed a phenomenological
model based on Pearson statistical framework to study the pseudorapidity
distribution. We have analyzed and fit data of charged hadrons produced in
$Pb-Pb$ collision at $2.76$ TeV and $Xe-Xe$ collision at $5.44$ TeV using the
proposed model.
