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))