Pointed scans

[mlincett]

I have been chatting with @G-Sommani about how to implement "pointed" scans. This would allow to investigate, for example, different sky regions compared to the one to which the scan iterations converge.

The relevant logic is in start_scan.py, specifically the PixelsToReco class.

The simplest way of doing this is to pass a pair of coordinates to start_scan and select only the first-NSIDE pixel pointed by those coordinates (instead of scanning/refining all the pixels). So the set of changes could be fairly minimal.

SkyDriver will have to support the additional input.

A logic of the type "scan a given radius around given coordinates" to determine the corresponding appropriate NSIDE sequence can be implemented in SkyMist.

[mlincett]

The logic is actually already implemented here:

skymap_scanner/skymap_scanner/server/pixels.py

Lines 168 to 174 in 11081a4

	def choose_pixels_to_reconstruct(
	nsides_dict: NSidesDict,
	nside_progression: NSideProgression,
	ang_dist: float = 2.,
	mc_ra_dec: Optional[Tuple[float, float]] = None,
	) -> Set[Tuple[icetray.I3Int, icetray.I3Int]]:
	"""Get more pixels to reconstruct/refine.

It scans at max NSIDE around given coordinates and with a given radius. It is a bit more rigid than what I had in mind but also does not require any additional logic.

All there is to do is back-propagate the arguments up to start_scan.py.

[ric-evans]

All there is to do is back-propagate the arguments up to start_scan.py.

Yup, I left this in. If you add the needed CI arguments, I can hook this up to SkyDriver.