def b_cov_TEE(lmax=200,frequency=353, psky=70, apo=2): """this is theory-based""" ls, (cltt, clee, clbb, clte) = fg.get_theory_fg(f=frequency, lmax=lmax, psky=psky, apo=apo) cov = calc_b_cov_TEB(cltt, clee, clee) hs = sg.get_hs(lmax=lmax) iszero = np.isclose(hs,0.) cov[iszero] = 0. return cov
#!/usr/bin/env python import sys import numpy as np from process_fullsky import simulate_alms, get_theory_cmb, PLANCK_DATA_PATH, FGS_SIM_PATH from foregrounds import get_theory_fg lmax = 2000 nside = 2048 filebase = sys.argv[1] filename_cmb = FGS_SIM_PATH + 'alms/' + '{}_cmb_alm_lmax{}'.format( filebase, lmax) filename_fg = FGS_SIM_PATH + 'alms/' + '{}_fg_alm_lmax{}'.format( filebase, lmax) ls, cls_theory_cmb = get_theory_cmb(lmax=lmax) ls, cls_theory_fg = get_theory_fg(lmax=lmax) Tlm_cmb, Elm_cmb, Blm_cmb = simulate_alms(cls_theory_cmb, nside=nside, lmax=lmax) Tlm_fg, Elm_fg, Blm_fg = simulate_alms(cls_theory_fg, nside=nside, lmax=lmax) np.savez(filename_cmb, Tlm=Tlm_cmb, Elm=Elm_cmb, Blm=Blm_cmb) np.savez(filename_fg, Tlm=Tlm_fg, Elm=Elm_fg, Blm=Blm_fg)
#!/usr/bin/env python import sys import numpy as np from process_fullsky import simulate_alms,get_theory_cmb,PLANCK_DATA_PATH,FGS_SIM_PATH from foregrounds import get_theory_fg lmax = 2000 nside = 2048 filebase = sys.argv[1] filename_cmb = FGS_SIM_PATH + 'alms/' + '{}_cmb_alm_lmax{}'.format(filebase,lmax) filename_fg = FGS_SIM_PATH + 'alms/' + '{}_fg_alm_lmax{}'.format(filebase,lmax) ls, cls_theory_cmb = get_theory_cmb(lmax=lmax) ls, cls_theory_fg = get_theory_fg(lmax=lmax) Tlm_cmb, Elm_cmb, Blm_cmb = simulate_alms(cls_theory_cmb, nside=nside, lmax=lmax) Tlm_fg, Elm_fg, Blm_fg = simulate_alms(cls_theory_fg, nside=nside, lmax=lmax) np.savez(filename_cmb, Tlm=Tlm_cmb, Elm=Elm_cmb, Blm=Blm_cmb) np.savez(filename_fg, Tlm=Tlm_fg, Elm=Elm_fg, Blm=Blm_fg)