-
Notifications
You must be signed in to change notification settings - Fork 21
Extend TOAST/sotodlib pointing comparison #1460
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
Conversation
tskisner
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks for improving these tests. We can't merge this in until the upstream toast PR is merged and a new tag is made.
|
Turns out that the 1/2 pixel offset was only needed before a recent @Bai-Chiang noticed that the exact weather parameters used in |
| box = np.radians([[dec_start, ra_start], [dec_stop, ra_stop]]) | ||
| reso = np.radians(1.0 / 60) | ||
| shape, wcs = enmap.geometry(box, res=reso, proj="car") | ||
| wcs.wcs.crpix[1] -= 0.5 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is this half pixel offset removed or not?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Here we create a geometry with non-integer center pixel. It is separate from TOAST applying an offset in translating between pixel numbers and pointing angles.
The half-pixel offset turned out to be not needed in TOAST.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Also it's better change the default weather to "toco" in mlmapmaker
sotodlib/sotodlib/toast/ops/mlmapmaker.py
Line 205 in 9fd7beb
| weather = Unicode("vacuum", help="Weather to assume when making maps") |
and add so3g_compat_mode to workflow script
sotodlib/sotodlib/toast/workflows/sim_observe.py
Lines 103 to 108 in 9fd7beb
| toast.ops.SimGround( | |
| name="sim_ground", | |
| weather="atacama", | |
| detset_key="pixel", | |
| session_split_key="wafer_slot", | |
| enabled=False, |
|
Excellent points! Implemented both. We won't be able to merge this PR until there is a tagged version of TOAST with |
0918323 to
3c0ab01
Compare


Check that TOAST and sotodlib assign same samples to same pixels given the same WCS configuration.
Unit tests will not pass until hpc4cmb/toast#894 is merged.