Measurements of two- and multiparticle azimuthal correlations provide valuable information on the properties of the quark--gluon plasma created in ultrarelativistic heavy-ion collisions. Further input can be achieved using the event shape engineering (ESE) technique which allows to select events corresponding to a specific initial geometry in a given centrality class. The elliptic ($v_2$) and triangular ($v_3$) flow coefficients of inclusive charged hadrons, $\pi^{\pm}$, $\rm{K}^{\pm}$, p+$\rm\overline{p}$, ${\rm K}^{0}_{\rm{S}}$, $\Lambda$+$\overline{\Lambda}$, and $\Xi^{-}$+$\overline{\Xi}^{+}$ are measured with ESE in Pb--Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$~TeV recorded by ALICE. In addition, the $v_2$ of $\pi^{\pm}$, $\rm{K}^{\pm}$, p+$\rm\overline{p}$, ${\rm K}^{0}_{\rm{S}}$, $\Lambda$+$\overline{\Lambda}$ from Xe--Xe collisions at $\sqrt{s_{\rm NN}} = 5.44$~TeV is presented. The results obtained using the scalar product method are reported as a function of transverse momentum, $p_{\rm T}$, within the pseudorapidity (rapidity) range |$\eta$| < 0.8 (|$y$| < 0.5) for various centrality classes. The flow coefficients exhibit a particle mass dependence for $p_{\rm T} < 3$~GeV/$c$ and a grouping according to particle type (i.e. mesons and baryons) for $3<p_{\rm T}<10$~GeV/$c$. The $v_2$ of $\pi^{\pm}$, $\rm{K}^{\pm}$, p+$\rm\overline{p}$ is described by MUSIC hydrodynamic calculations for $p_{\rm T}<1$~GeV/$c$. The ratios between $v_{\rm n}$ for shape selected and unbiased events show a weak dependence on transverse momentum for $p_{\rm T}<3$~GeV/$c$ and are almost flat at higher $p_{\rm T}$, independent on particle species. These measurements provide new and independent constraints for the initial conditions and system properties of nuclear matter created in heavy-ion collisions.