def slice_clusters_idx(rm_pd, z_bins):
    """
    Returns all clusters from the DataFrame of clusters between
    a range of specified redshift limits.

    Parameters
    -----------
    rm_pd : `pd.DataFrame`
        DataFrame of cluster centers with ['RA','DEC','Z_LAMBDA'] information.

    z_bins : `numpy.ndarray`
        Array of the left and right bin edges of the different redshift bins.

    Returns
    -----------
    cluster_idx_arr : `list`
        List of indices of the cluster centers for each of the different
        redshift bins.
    """
    # Array of cluster indices
    cluster_idx_arr = [[] for x in range(len(z_bins))]
    # Populating array with cluser indices
    for kk, (zz_low, zz_high) in enumerate(z_bins):
        # Cluster indices
        cluster_idx_kk = rm_pd.loc[ (rm_pd['Z_LAMBDA'] >= zz_low) &
                                    (rm_pd['Z_LAMBDA'] < zz_high)].index.values
        # Saving to list
        cluster_idx_arr[kk] = cluster_idx_kk

    return cluster_idx_arr