Abstract Measurements of elliptic (v 2) and triangular (v 3) flow coefficients of π ±, K±, p+ p ¯ $$ \overline{\mathrm{p}} $$ , K S 0 $$ {\mathrm{K}}_{\mathrm{S}}^0 $$ , and Λ+ Λ ¯ $$ \overline{\Lambda} $$ obtained with the scalar product method in Xe-Xe collisions at s NN $$ \sqrt{s_{\mathrm{NN}}} $$ = 5.44 TeV are presented. The results are obtained in the rapidity range |y| < 0.5 and reported as a function of transverse momentum, p T, for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for p T < 3 GeV/c, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3 < p T < 8 GeV/c). The magnitude of the baryon v 2 is larger than that of mesons up to p T = 6 GeV/c. The centrality dependence of the shape evolution of the p T-differential v 2 is studied for the various hadron species. The v 2 coefficients of π ±, K±, and p+ p ¯ $$ \overline{\mathrm{p}} $$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for p T < 1 GeV/c. A comparison with v n measurements in the corresponding centrality intervals in Pb-Pb collisions at s NN $$ \sqrt{s_{\mathrm{NN}}} $$ = 5.02 TeV yields an enhanced v 2 in central collisions and diminished value in semicentral collisions.