galsci / pysm

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

Simulate CO line emission map #113

Closed mousset closed 2 years ago

mousset commented 2 years ago

Hi, I am trying to simulate the CO emission line at 230.5 GHz, in µK_CMB. Here is the code I use:

$ sky = pysm3.Sky(nside=128, preset_strings=["co1"], output_unit=u.uK_CMB) $ co_emission = sky.get_emission(freq=[230.4, 230.6]*u.GHz)

I think I am not doing the right thing because the intensity map I get is very small compared to the dust map I made using the d1 model. Could you please tell me how to generate such a map? Thanks in advance.

zonca commented 2 years ago

Which version of PySM are you using?

mousset commented 2 years ago

I am using version 3.4.0b1.

zonca commented 2 years ago

thanks @mousset, I confirm there was a bug I introduced recently in unit conversion in CO (and messed up the tests as well so they agreed!)

It should be fixed in https://github.com/galsci/pysm/pull/114

zonca commented 2 years ago

@mousset merged into main, could you please confirm if it is fixed?

mousset commented 2 years ago

Yes I confirm, it works perfectly now. Thanks!