galsci / pysm

PySM 3: Sky emission simulations for Cosmic Microwave Background experiments
https://pysm3.readthedocs.io/
BSD 3-Clause "New" or "Revised" License
33 stars 23 forks source link

Synchrotron beta with local-spectra modulation #151

Closed zonca closed 1 year ago

zonca commented 1 year ago

This is related to #148.

I am applying the same modulation generated for the Synchrotron templates also for the synchrotron spectral index.

image

If I disable the modulation I get: image

The same plot for the latest beta release, from https://pysm3.readthedocs.io/en/latest/preprocess-templates/synchrotron_beta.html:

image

The modulation map itself is between a factor of 2 and a factor of 3, shouldn't it have a mean of 1?:

image

@giuspugl @brandonshensley @seclark should I just normalize it? However if I normalize it for beta should I also normalize it for the templates?

zonca commented 1 year ago

The synchrotron beta notebook: https://gist.github.com/7d60503160fa512987474707f7006853 The synchrotron templates notebook: https://github.com/galsci/pysm/blob/9db4fcb4dcb04f411b6f552501f45a0cb7af78ef/docs/preprocess-templates/synchrotron_template_logpoltens.ipynb

zonca commented 1 year ago

The GNILC temperature modulation for comparison: image

zonca commented 1 year ago

also @NicolettaK @delabrou

giuspugl commented 1 year ago

Hi @zonca , i was expecting this to happen! I wouldn't touch the modulation map of Beta_s as this uses large scale templates that are derived from the old s1 + small scales fitted from SPASS data. If you agree i'd use the old modulation template of the previous beta release. That is derived from the large scale of the intensity template . What do you and the others think ?

zonca commented 1 year ago

ok @giuspugl, implemented in https://github.com/galsci/pysm/pull/152