Dividir una lista basada en otros valores de lista en Mathematica

En Mathematica tengo una lista de coordenadas de puntos

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

y una lista de índices de clúster a los que pertenecen estos puntos

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

¿Cuál es la forma más fácil de dividir los puntos en dos listas separadas según los valores de clusterIndices?

EDITAR: La solución que se me ocurrió:

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

¿Hay una mejor manera de hacer esto?

Respuestas a la pregunta(6)

Su respuesta a la pregunta