def run_check(config, verbose=True):
if config.noise_only:
if not config.do_pencil_beam:
if verbose:
prompter.prompt_warning("No beam convolution required for noise only simulation.")
config.do_pencil_beam = True
if config.do_pencil_beam:
if config.oversampling_rate != 1:
if verbose:
prompter.prompt_warning("No oversampling required for pencil beam")
config.oversampling_rate = 1
def get_sky_map(self):
if self.config.sim_pol_type == "noise_only":
self.sky_map = None
elif self.config.sim_pol_type == "T":
self.sky_map = hp.read_map(self.config.input_map, verbose=False)
elif self.config.sim_pol_type == "QU":
self.sky_map = hp.read_map(self.config.input_map, field=(0,1), verbose=False)
elif self.config.sim_pol_type == "_QU":
self.sky_map = hp.read_map(self.config.input_map, field=(1,2), verbose=False)
else:
self.sky_map = hp.read_map(self.config.input_map, field=(0,1,2), verbose=False)
if not self.config.sim_pol_type == "noise_only":
map_nside = hp.get_nside(self.sky_map)
if map_nside != self.config.nside_in:
prompter.prompt_warning("NSIDE of config does not match NSIDE of map")
sys.exit()
self.npix = hp.nside2npix(self.config.nside_in)
