We present the doubly differential invariant yield of primary charged particles π + , π − , K + , K − and p , p ̄ in Xe–Xe collisions at s N N = 5.44 TeV, as a function of p T in kinematic region of 0 < p T < 1 5 GeV / c and for the pseudorapidity interval of | η | < 0. 8 by using PYTHIA8 tuned to different Parton Showers like Simple Showers (Time-Like, Space-Like), DIRE and VINCIA Showers. The predictions of hard and soft QCD processes of Parton Showers are then compared with the ALICE data at different centrality classes of 0–5%, 50–60% and 70–80%. The observed deviations are connected with the kinematics involved and are mostly due to the color coherence and soft gluon emission from parton pairs. As PYTHIA8 does not take into account the Cronin effect, thus the deviation at low p T is due to Cronin effect. The nuclear modification factor is constructed and it is observed that the signal is less than unity or suppression is seen which is coming from Parton Showers. [ABSTRACT FROM AUTHOR]