Unified mapmaking for an anisotropic stochastic gravitational wave background

A stochastic gravitational wave background (SGWB), created by the superposition of signals from unresolved astrophysical sources, may be detected in the next few years. Estimating the variation of intensity across the sky can, therefore, play a key role in improving our understanding of astrophysical models. Skymaps have been produced for all the data-taking runs of the advanced ground-based interferometric detectors. While these maps are being produced in pixel and SpH bases, to probe, respectively, localized and diffuse astrophysical and cosmological sources, with algorithms that employ cross-correlation as the common strategy, the underlying algebra and numerical implementation remain different. As a consequence, there was a need for producing skymaps in both bases in those analyses. We show that these manifestly redundant methods could indeed be unified to a single analysis that can probe very different scales and demonstrate it by applying them on real data. We first develop the algebra to show that the results in two different bases are easily transformable. We then incorporate both the schemes in the now-standard analysis pipeline for anisotropic SGWB, PyStoch. This will enable SGWB anisotropy searches in SpH basis also to take full advantage of integrated \hpx tools and makes it computationally feasible to perform the search in every frequency bin. We, however, follow a different approach for direct estimation of the SpH moments. We show that the results obtained from these different methods match very well; the differences are less than $0.1$\% for the SpH moments and less than $0.01$\% for the Fisher information matrices. Thus we conclude that a single skymap will be sufficient to describe the anisotropies in a stochastic background. The multiple capabilities of PyStoch will be useful for estimating and constraining various measures that characterize an anisotropic background.