def make_covTD(d): """ Usually coverage map is provided in a separate file. But if not the case, this method can compute a coverage map Parameters: d: Qubic dictionary Return: cov: coverage map in a.QubicMultibandAcquisition shape (nfreq, npix). """ pointing = qubic.get_pointing(d) q= qubic.QubicMultibandInstrument(d) s= qubic.QubicScene(d) nf_sub_rec = d['nf_recon'] _, nus_edge, nus, _, _, _ = qubic.compute_freq(d['filter_nu'] / 1e9, nf_sub_rec, d['filter_relative_bandwidth']) arec = qubic.QubicMultibandAcquisition(q, pointing, s, d, nus_edge) cov = arec.get_coverage() return cov
# about closed horns.
from __future__ import division, print_function
from matplotlib.pyplot import *
import healpy as hp
from astropy.io import fits
import qubic
basedir = '/home/louisemousset/QUBIC/MyGitQUBIC'
dictfilename = basedir + '/qubic/qubic/scripts/global_source.dict'
d = qubic.qubicdict.qubicDict()
d.read_from_file(dictfilename)
# ========== Sweeping pointing in a fix frame ===========
p = qubic.get_pointing(d)
nptg = len(p.index) # Pointing number
el_angles = np.unique(p.elevation) # all different elevation angles
az_angles = np.unique(p.azimuth) # all different azimuth angles
print('el_angles = ' + str(el_angles))
print('az_angles = ' + str(az_angles))
figure('pointing')
subplot(3, 1, 1)
plot(p.index, p.azimuth, '.')
xlabel('index pointing')
ylabel('azimuth angle')
subplot(3, 1, 2)
plot(p.index, p.elevation, '.')
xlabel('index pointing')
