Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Discrepancy between synchrotron models s4, s5 and s7 and Planck LFI 30 GHz data #130

Open
1cosmologist opened this issue Oct 29, 2022 · 8 comments

Comments

@1cosmologist
Copy link

This is an effort to validate the PySM3 synchrotron models by comparing with Planck LFI 30 GHz (NPIPE Full, A and B splits) and WMAP 23 GHz (9year) maps. For this I generate synchrotron maps for models s4, s5, and s7 at 30 GHz at NSIDE=512. A color correction factor of 0.969 is multiplied to the LFI 30GHz, and a factor of 0.427 is multiplied to the WMAP K band. I also use an inverse noise variance weighted computation of the synchrotron map at 30 GHz by Jacques Delabrouille. All maps are reconvolved with a 2 degree beam. The power spectrum is computed for three cases: A. 70% galactic synchrotron mask, B. 33% mask of the northern part of the galactic mask and C. 37% mask of the southern part of galactic mask.
lfi30R4_70pc_filledmask

All The EE and BB power spectra computed from models and maps for 70% mask:
DlEE_PySMmodel_vs_crossinv
DlBB_PySMmodel_vs_crossinv

For north sky:
DlEE_PySMmodel_vs_crossinv_northsky
DlBB_PySMmodel_vs_crossinv_northsky

For south sky:
DlEE_PySMmodel_vs_crossinv_southsky
DlBB_PySMmodel_vs_cross_southsky

Clearly, the models show discrepancies when compared with the data. This discrepency is very much present in the multipole range important for primordial B-mode science.

@brandonshensley
Copy link
Contributor

Thanks @1cosmologist, a lot to dig into here!

First question: where does the factor of 0.969 come from? I.e., can you point me to a paper?

Second question: overall, the agreement looks pretty good to me. Are you most concerned with the BB spectra in the northern sky? Have you done any map-based comparisons?

@seclark
Copy link

seclark commented Oct 29, 2022

It would also be helpful to see the masks used for each of these. And a caution for our interpretation that ~1/2 the points in these plots are below the beam scale.

@seclark
Copy link

seclark commented Oct 29, 2022

And let me also add my thanks and I'm eager to dig into this together at Monday's hack session!

@1cosmologist
Copy link
Author

Thanks @1cosmologist, a lot to dig into here!

First question: where does the factor of 0.969 come from? I.e., can you point me to a paper?

Second question: overall, the agreement looks pretty good to me. Are you most concerned with the BB spectra in the northern sky? Have you done any map-based comparisons?

Hey @brandonshensley, the color correction factor is from the Planck 2013 V: LFI calibration (https://arxiv.org/abs/1505.08022), table 1. Note this bit different from the factor Jacques calculated using the LFI 30 GHz bandpass. We are sticking to the LFI team's value for average color correction factor here.

I have indeed done map based comparisons. The figure below is the relative difference of the polarised intensity s5-LFI30/LFI30.
s5-LFI30_relpolintensity

We also look at scatter of the log(pol intensity) of s5 vs LFI30 at NSIDE=16. We would expect a slope of 1 but there is a clear tilt. The black line is slope 1, red is the actual trend line.
s5_v_30_scatter

I think the agreement needs to be looked at bit closely, in particular in the patches of interest for CMB experiments. While the general agreement is reasonable, except in the north sky BB, are we OK with the agreement for 20 < ell < 100?

@1cosmologist
Copy link
Author

It would also be helpful to see the masks used for each of these. And a caution for our interpretation that ~1/2 the points in these plots are below the beam scale.

Hi @seclark, here are the masks:
lfi30_34pc_northskymask
lfi30_39pc_southskymask

I agree that the beam is perhaps too large. It might be helpful to a power spectrum comparison at native resolution. It might be better to get a bandpass integrated synchrotron map for the models. I will add it to my list.

@1cosmologist
Copy link
Author

And let me also add my thanks and I'm eager to dig into this together at Monday's hack session!

How long is the hack session planned for? It is a bit early for me, but I will try to join in 30 min or an hour later.

@seclark
Copy link

seclark commented Oct 29, 2022

How long is the hack session planned for? It is a bit early for me, but I will try to join in 30 min or an hour later.

It's scheduled for an hour, 7-8am PT. I will have to jump at 8 but discussions could certainly continue longer if people are keen. If you're able to have some overlap in that timeframe that would be great, but I understand it's early!

@1cosmologist
Copy link
Author

the color correction factor is from the Planck 2013 V: LFI calibration (https://arxiv.org/abs/1505.08022), table 1.

CORRECTION: https://arxiv.org/abs/1303.5066

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants