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

Problem with Log Pol Tens formalism intensity transformation ? #109

Closed giuspugl closed 2 years ago

giuspugl commented 2 years ago

Recently, an issue has been raised on the transformation of the intensity component within the Log Pol tens formalism . I address it here

Screen Shot 2022-02-15 at 6 21 06 PM below the maps how they look like image

Finally I asserted that :

i would conclude that this issue shouldn't be worrisome for time being. what do you think ?

delabrou commented 2 years ago

I agree with all of this, the only issue is if the generation of the small scales depends locally on the total i (or i'). Indeed, typically, when one injects small scales, those are weighted by the local total dust emission level, and an offset in the map changes this. Unless you weight the generation of small scales with some different weighting, or you do some re-scaling later on, there should be a difference. We don't see anything in full-sky maps. Just to check, do you think you could show the original map, the small scales as generated (using the same seed), and the sum in two cases: when your original I and P maps are in uK_RJ, and when they are in K_RJ, specifically at the North Galactic pole, for a gnomonic projection?

giuspugl commented 2 years ago

Below the plot as @delabrou prescribed , don't see any difference apart from the fact that maps have different units... Maps are 9x9 deg2 stack in the Gal. north pole

image