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)