From 595f5a2be0e0257c1ef0168ecf7d31025c61c608 Mon Sep 17 00:00:00 2001
From: Aaron LI <aaronly.me@outlook.com>
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