-
Notifications
You must be signed in to change notification settings - Fork 22
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 d12 and color corrected HFI353 #128
Comments
thanks @erussier, can you please repeat the test at Nside 2048 and update the notebook? |
also to test the implementation I am comparing with this reference notebook with a simple analysis ran at Nside 8: can you check there if you find any errors? |
@erussier can you please adapt the notebook to run at NSIDE 2048 instead of 512? |
Hi @zonca, I ran the notebook at NSIDE 2048 except for the MKD map since the native NSIDE of the MKD map I have is 512. So if you want to run this notebook at NSIDE 2048, please change the NSIDE in the section "Plot map of Planck MKD model" and the path to the MKD map. Here it is: https://gist.github.com/erussier/89b656bba5c65331b5764d39754b7d38 |
the MKD is the most important map to have at 2048, @delabrou ran the MKD model at Nside 2048, those are the maps that I am using as inputs to PySM. |
Hi @zonca, I found the MKD map at NSIDE = 2048 and indeed we get the same results between the MKD model at 353 GHz and d12. I guess we can close this issue. |
Good news! Thanks |
In fact, the mkd map I was using at NSIDE 2048 and used in the PySM pipeline, was generated without the color correction implemented. We need to multiply the mkd map by 0.92 to take into account the color correction. Probably the best thing to do is to change the input maps. @delabrou will generate new maps to use as inputs for PySM. |
Sure, I can then replace the inputs |
@erussier any news on the new maps? |
@erussier @delabru what if I just apply the |
Sorry for the late reply, we think that is okay to apply the 0.911 factor inside PySM, thank you. |
Yes, I confirm. You multiply all IQU templates by the colour correction factor (in the PSM I use 0.92). |
Do not do anything to the temperature and spectral index maps. |
Hello, I plotted the PySM d12 dust power spectra and compared to the HFI353GHz. There seems to be an offset at large scales that we do not find when comparing the MKD map (without color correction) from here (https://apc.u-paris.fr/Downloads/Planck/PSM-Data/MKD-MODEL-PySM-2048/skyinbands/HFI/detector_F353/) and the HFI353GHz without color correction.
![download](https://user-images.githubusercontent.com/54247542/197905131-9f899138-c2eb-43af-8e40-aeb19e11b0dc.png)
![download](https://user-images.githubusercontent.com/54247542/197905121-0f6eb97a-7600-4869-8a93-26462b88d75c.png)
This is what I get for Emodes and Bmodes.
In principle, the PSM MKD and HFI353 without color correction should be compatible and d12 and color corrected HFI353 should also be compatible at least on the largest scales.
Here is the jupyter notebook where I did this comparison: comp_dust_PS_d12_MKD_map.ipynb
The text was updated successfully, but these errors were encountered: