Closed segasai closed 5 months ago
@segasai: Thanks! I'll take this code and run with it so I can make it fit some of the desitarget
conventions. I'll also add the cross-match between Gaia and DR9 of the Legacy Surveys. I'll also think about how to handle the dependency on gaiadr3_zeropoint
.
I'm traveling over the U.S. long weekend. But, at some point next week I'll put together a PR based on this code and send it to you for review.
Great ! Thank you @geordie666 !
Just a small remark, ideally it'd be good if functions rotate_pm, sphere_rotate, correct_pm were not changed (as those are essentially external functions copied here, so if their behaviour was different from what we use elsewhere, it'll be source of pain/confusion). Also if we ever add more streams, likely the basic structure of the selection code will be the same and it'll still rely on those utility functions.
Added in #814.
Hi,
This is the starting point to iterate over the integration of the selection of GD1 targets into desitarget following the discussion yesterday. (@geordie666 , @sazabi4 )
The code below is based on my notebook and is self-contained, in the sense that I integrated several dependencies there (i.e. various rotation, proper motion transformations from github.com/segasai/astrolibpy, specifically ( i.e.https://github.com/segasai/astrolibpy/blob/master/my_utils/rotate_pm.py https://github.com/segasai/astrolibpy/blob/master/my_utils/sphere_rotate.py ).
The main selection code is in the main section. The reading of the data is in read_data() and that's what need to be replaced by the DECals X Gaia DR3 match from desitarget.
The only non-standard dependency that I cannot easily get away from is the gaia_zeropoint https://gitlab.com/icc-ub/public/gaiadr3_zeropoint which is necessary to fix the offset in Gaia parallaxes.
I can probably create a PR based on this code, provided there are no major objections to how things are structured now.