import numpy as np

def get_background_full(norm_adata, threshold, n_cells_low, n_cells_high):

    thr_per_gene = np.sum(norm_adata.X > threshold, axis=0).flat
    genes_filter = np.logical_and(thr_per_gene >= n_cells_low, thr_per_gene <= n_cells_high)
    print("Background genes: " + str(np.sum(genes_filter)))

    norm_adata.var.drop(columns="CIARA_background", inplace=True, errors='ignore')
    norm_adata.var.insert(0, "CIARA_background", genes_filter.tolist())

    return
