n Mathematica habe ich eine Liste von Punktkoordinat

size = 50;
points = Table[{RandomInteger[{0, size}], RandomInteger[{0, size}]}, {i, 1, n}];

und eine Liste von Cluster-Indizes, zu denen diese Punkte gehören

clusterIndices = {1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1};

Was ist der einfachste Weg, um die Punkte basierend auf den clusterIndices-Werten in zwei separate Listen aufzuteilen?

EDIT: Die Lösung, die ich mir ausgedacht habe:

pointIndices =
  Map[#[[2]] &,
    GatherBy[MapIndexed[{#1, #2[[1]]} &, clusterIndices], First],
    {2}];
pointsByCluster = Map[Part[points, #] &, pointIndices];

Es gibt einen besseren Weg, dies zu tun?