From 595f5a2be0e0257c1ef0168ecf7d31025c61c608 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Tue, 13 Dec 2016 18:46:05 +0800 Subject: galactic/snr.py: Use HEALPix pixel area if object smaller than it --- fg21sim/galactic/snr.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/fg21sim/galactic/snr.py b/fg21sim/galactic/snr.py index bc39c7e..d331a88 100644 --- a/fg21sim/galactic/snr.py +++ b/fg21sim/galactic/snr.py @@ -88,8 +88,9 @@ class SuperNovaRemnants: self.checksum = self.configs.getn("output/checksum") self.clobber = self.configs.getn("output/clobber") self.nside = self.configs.getn("common/nside") + self.pixsize = hp.nside2resol(self.nside, arcmin=True) / 60.0 + self.pixarea = self.pixsize ** 2 # [ deg^2 ] self.freq_unit = au.Unit(self.configs.getn("frequency/unit")) - # logger.info("Loaded and set up configurations") def _load_catalog(self): @@ -224,6 +225,9 @@ class SuperNovaRemnants: freq_ref = self.catalog_flux_freq # [ MHz ] Fnu = flux * (frequency / freq_ref) ** (-specindex) # [ Jy ] omega = size[0] * size[1] # [ deg^2 ] + if omega < self.pixarea: + # The object is smaller than a pixel, so round up to a pixel area + omega = self.pixarea Tb = Fnu_to_Tb_fast(Fnu, omega, frequency) return Tb -- cgit v1.2.2