[fix] make sure sources smaller than the pix are masked

LiYunyang · LiYunyang · commit 66be5c101d17 · 2025-09-01T15:41:59.000-05:00
diff --git a/maps/src/maputils.cxx b/maps/src/maputils.cxx
@@ -564,6 +564,9 @@ MakePointSourceMask(const G3SkyMap &map, const std::vector<double> & ra,
 
 	for (size_t i = 0; i < ra.size(); i++) {
 		auto pixels = map.QueryDisc(ra[i], dec[i], radius[i]);
+		uint64_t _ipix = map.AngleToPixel(ra[i], dec[i]);
+		if (_ipix < (size_t)-1)
+			pixels.emplace_back(_ipix);
 		for (auto p: pixels)
 			(*mask)[p] = true;
 	}

Original file line number	Diff line number	Diff line change
`@@ -564,6 +564,9 @@ MakePointSourceMask(const G3SkyMap &map, const std::vector<double> & ra,`
`564`	`564`
`565`	`565`	`for (size_t i = 0; i < ra.size(); i++) {`
`566`	`566`	`auto pixels = map.QueryDisc(ra[i], dec[i], radius[i]);`
	`567`	`+ uint64_t _ipix = map.AngleToPixel(ra[i], dec[i]);`
	`568`	`+ if (_ipix < (size_t)-1)`
	`569`	`+ pixels.emplace_back(_ipix);`
`567`	`570`	`for (auto p: pixels)`
`568`	`571`	`(*mask)[p] = true;`
`569`	`572`	`}`