In order to characterize the higher order moments of the particle multiplicity, we implement the linear-sigma model with Polyakov-loop correction. We first studied the critical phenomena and estimated some thermodynamic quantities. Then, we compared all these results with the first--principle lattice QCD calculations. Then, the extensive study of non-normalized four moments is followed by investigating their thermal and density dependence. We repeat this for moments normalized to temperature and chemical potential. The fluctuations of the second order moment is used to estimate the chiral phase--transition. Then, we implement all these in mapping out the chiral phase transition, which shall be compared with the freeze-out parameters estimated from the lattice QCD simulations and the thermal models are compared with the chiral phase--diagram.