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)