for sim_idx in subtasks:
log.info("processing %d" %sim_idx)
if plot_only: continue
tmap, qmap, umap = get_sim(sim_idx)
tmap -= np.mean(tmap)
qmap -= np.mean(qmap)
umap -= np.mean(umap)
qmap *= -1.
tmap *= taper
qmap *= taper
umap *= taper
emap, bmap = qu2eb(qmap, umap)
if sim_idx == 0:
io.high_res_plot_img(tmap, output_path('tmap_unlen_%d.png'%sim_idx), down=3)
io.high_res_plot_img(qmap, output_path('qmap_unlen_%d.png'%sim_idx), down=3)
io.high_res_plot_img(umap, output_path('umap_unlen_%d.png'%sim_idx), down=3)
io.high_res_plot_img(emap, output_path('emap_unlen_%d.png'%sim_idx), down=3)
io.high_res_plot_img(bmap, output_path('bmap_unlen_%d.png'%sim_idx), down=3)
cmb_dict = {'t': tmap, 'e': emap, 'b': bmap}
# take take tt and te spectra
for polcomb in polcombs:
if polcomb in ['ee', 'eb', 'bb', 'pp']: continue
if st.has_data('frac%s_deconv'%polcomb, sim_idx): continue
emap1 = cmb_dict[polcomb[0]]
emap2 = cmb_dict[polcomb[1]]
log.info('[add_spectra]: ' + polcomb)
l, cl = cusps_fc.get_power(emap1, emap2=emap2, polcomb=polcomb.upper())
for sim_idx in subtasks:
log.info("processing %d" % sim_idx)
if st.has_data(sim_idx, 'fracbb_deconv'): continue
tmap, qmap, umap = get_sim(sim_idx)
tmap -= np.mean(tmap)
qmap -= np.mean(qmap)
umap -= np.mean(umap)
tmap *= taper
qmap *= taper
umap *= taper
tmap, emap, bmap = tqu2teb(tmap, qmap, umap)
if sim_idx == 0:
io.high_res_plot_img(tmap,
output_path('tmap_unlen_%d_%s.png' %
(sim_idx, coords_str)),
down=3)
io.high_res_plot_img(qmap,
output_path('qmap_unlen_%d_%s.png' %
(sim_idx, coords_str)),
down=3)
io.high_res_plot_img(umap,
output_path('umap_unlen_%d_%s.png' %
(sim_idx, coords_str)),
down=3)
io.high_res_plot_img(emap,
output_path('emap_unlen_%d_%s.png' %
(sim_idx, coords_str)),
down=3)
io.high_res_plot_img(bmap,
output_path('bmap_unlen_%d_%s.png' %
def make_plots(prefix, emaps):
for idx, cmb_type in enumerate(['I', 'Q', 'U']):
file_name = "{}_{}.png".format(prefix, cmb_type)
io.high_res_plot_img(emaps[idx], file_name, down=3)