wstr = '$W (1/\mu K \\times 10^7)$' else: wstr = '$W$ (dimensionless)' #pl = io.Plotter(xyscale='loglin',xlabel='$\\ell$',ylabel=wstr,ftsize=16) pl = io.Plotter(xyscale='linlin', xlabel='$\\ell$', ylabel=wstr, ftsize=16) # !!! for i in range(len(qids)): col = cols[i] qid = qids[i] lmin, lmax, hybrid, radial, friend, cfreq, fgroup, wrfit = aspecs(qid) if tutils.is_lfi(qid): ls = "-." lab = "LFI %d GHz" % cfreq elif tutils.is_hfi(qid): ls = "--" lab = "HFI %d GHz" % cfreq else: ls = "-" aind = qid.split("_")[1] lab = actmap[qid] #"ACT_%s %d GHz" % (aind,cfreq ) mul = 1e7 if comp == 'comptony' else 1 cents, w1d = np.loadtxt("weights_%s_%s_%s.txt" % (comp, version, lab), unpack=True) pl.add(cents, w1d * mul, label=lab if comp == 'comptony' else None, ls=ls, color=col)
pl.add(ncents,n1d/fbeam1(ncents)/fbeam2(ncents)) pl._ax.set_ylim(1,1e8) pl.done("%ss1d_%s_%s.png" % (opath,qid1,qid2),verbose=False) #res = s1d/fbeam1(ccents)/fbeam2(ccents) - stheory.get_theory_cls(f1,f2,a_cmb=1,a_gal=0,a_cibp=0,a_cibc=0,a_radps=0,a_ksz=0,a_tsz=1) res = s1d - stheory.get_theory_cls(f1,f2,a_cmb=1,a_gal=0,a_cibp=0,a_cibc=0,a_radps=0,a_ksz=0,a_tsz=1) if tutils.is_planck(qid1) and tutils.is_planck(qid2): flmin = 20 else: flmin = 1000 if not(tutils.is_planck(qid1)) and not(tutils.is_planck(qid2)): flmax = 5000 elif tutils.is_hfi(qid1) and tutils.is_hfi(qid2): flmax = 2000 elif (not(tutils.is_planck(qid1)) and tutils.is_hfi(qid2)) or (not(tutils.is_planck(qid2)) and tutils.is_hfi(qid1)): flmax = 3000 elif (qid1=='p03' or qid2=='p03') and (qid1 not in ['p01','p02']) and (qid2 not in ['p01','p02']): flmax = 1000 else: flmax = 300 if flmax<=flmin: #!!!!!! io.save_cols("%sfgcov_%s_%s.txt" % (opath,qid1,qid2),(fells,fells*0)) else: print("Rank %d doing task %d for array %s x %s with lmin %d and lmax %d ..." % (rank,task,qids[0],qids[1],flmin,flmax))