def score_programs(adata_here,program_name='bulk.guide.program',
                  pref=None,copy=False):
    
    """Get program score for each cell, by averaging program genes

    """

    if copy: adata_here = adata_here.copy()
    
    if pref is None:
        pref=program_name
    programs=list(set(adata_here.var[program_name]))
    for pro in programs:
        print('scoring',pro)
        pro_genes=adata_here.var_names[adata_here.var[program_name]==pro]
        adata_here.obs[pref+str(pro)]=adata_here[:,pro_genes].X.mean(axis=1)
        
    if copy:
        return(adata_here)