Plotting the Scattering Covariance coefficients

[giuspugl]

Hello,
first of thanks a lot for making this package publicly available and documented! Nice work!
I'd like to visualize the 4 SC S1, S2, S3, S4 coefficients and compare with different fields in the context of post component separated CMB B-mode maps .
I went through the notebook you provided
https://github.com/astro-informatics/s2scat/blob/main/notebooks/auto_compression.ipynb
and am trying to figure out the right way to plot em all, perhaps I 've already figured out how to plot S1 and S2 (which are the easiest), starting from a healpix map with a power law power spectrum, than I estimate the coefficients as in this notebook
https://colab.research.google.com/drive/1rCrOjVCIYsqF7RHUp2gh2vckQaXUD-wY?usp=sharing
and got this plots (for a given \gamma direction)

Do you think am finding anything meaningful ?
any of your insight and/or suggestions for S3 and S4 are more than welcome too...