def MDtest_get_glm(Nreal=1,minreal=0,Ndesbins=[2,3],nvss=True):
cldat=MDtest_get_Cl(justread=True,Ndesbins=Ndesbins,nvss=nvss)
print 'cldat.bintaglist:',cldat.bintaglist
#cldat seems to be associated with correct crosspairs
# also, manipulating cl data has expected effects on rho hist
rlzns=np.arange(minreal,minreal+Nreal)
#leaving reclist empty means just gen gal and isw maps
au.getmaps_fromCl(cldat,rlzns=rlzns)
def MDtest_get_glm_and_rec(Nreal=1,minreal=0,justgetrho=0,dorho=1,Ndesbins=[2,3],lmin=3,lmax=80,rhofiletag='',nvss=True):
#adding this block of text fixed bug when passing two sets of analysis. why?
maptypes=MDtest_get_maptypelist(Ndesbins=Ndesbins,nvss=nvss)
mapsfor=[mt.tag for mt in maptypes] #tags for maps we want to make
mapsfor.append('isw')
allcldat=MDtest_get_Cl(justread=True,Ndesbins=Ndesbins,nvss=nvss)
#cldat=get_reduced_cldata(allcldat,dothesemaps=mapsfor)#basically just reorders
cldat=allcldat
makeplots = Nreal==1
rlzns=np.arange(minreal,minreal+Nreal)
reclist=MDtest_get_reclist(Ndesbins=Ndesbins,lmin=lmin,lmax=lmax,nvss=nvss)
au.getmaps_fromCl(cldat,rlzns=rlzns,reclist=reclist,justgetrho=justgetrho,dorho=dorho,dos=False,dochisq=False,rhofiletag=rhofiletag)
