This work provides a description of the critical threshold phenomenon in multidimensional restricted Euler-Poisson (REP) equations, introduced in [H. Liu and E. Tadmor, Comm. Math. Phys. 228 (2002), 435-466]. For three-dimensional REP equations, we identified both upper-thresholds for finite time blow up of solutions and sub-thresholds for global existence of solutions, with thresholds depending on the relative size of the eigenvalues of the initial velocity gradient matrix and the initial density. For attractive forcing case, these one-sided threshold conditions of initial configurations are optimal, and the corresponding results also hold for arbitrary n-dimensions (n ≥ 3).