Welcome! This is a small repository containing the final selection function for the APOGEE survey, along with some simple tutorials for how to apply it for your science case.
If you have any questions, please reach out to me at [email protected]!
apogeeCombinedSF.dat.zip
- the pre-computed raw (2D) selection function for APOGEE DR17 computed on December 16, 2022.TUTORIAL1.ipynb
- a jupyter notebook tutorial for how to work with the raw selection function.TUTORIAL2.ipynb
- a jupyter notebook tutorial for how to calculate the effective (3D) selection function from the raw selection function.JHK_isochrones.txt
- sample PARSEC isochrones used inTUTORIAL2
Additional Files Coming Soon:
effective_selfunc.dat
(coming soon!) - the pre-computed effective selection function used in Imig et al. 2023 (in prep)
- This selection function is only valid for APOGEE main survey targets, stars that were targeted randomly from 2MASS. You can easily select main survey targets in the APOGEE
allStar
file using the criterionallStar['EXTRATARG']==0
. Do not use this selection function for non-main survey targets, as it will be an underestimate of the true selection fraction for those stars. - Some of the halo stars were not targeted on 2MASS photometry, but on additional Washington+DDO51 photometry (more detail found here). Therefore, the 2MASS-based selection function here may not be accurate for all halo fields. If you plan to use this selection function to study the Galactic halo, be cautious of this! In most cases, this results in an overestimate of the selection fraction for halo fields, in a few cases even exceeding 100%. Relevant stars are flagged with the
APOGEE_TARGET1
targeting bitsAPOGEE2_WASH_GIANT
orAPOGEE2_WASH_DWARF
(documented more here) to easily identify the affected fields. - There is a minor known issue for 14 (out of 1937) fields in the current selection function where the color limit returns NaN. A fix coming soon!!
You'll need an installation of the apogee
and mwdust
packages for this, along with some standard python packages (numpy
, astropy
, matplotlib
, etc.)
If you use this selection function for published work, please cite Imig et al. 2023 (in prep) for this selection function, and Bovy et al. 2016a or Bovy et al. 2016b for the original apogee
code used to calculate the selection function.
The raw selection function, which I like to think of as the 2D selection function, reflects the targeting strategies of APOGEE. Targets in APOGEE were selected randomly from the 2MASS catalog1, so the raw selection function is the ratio of the number of observed stars in APOGEE divided by the total number of possible targets in 2MASS, for each location on the 2D sky (APOGEE field).
1: This description is a vast oversimplification of the targeting strategies! For more detail, see the APOGEE targeting webpage, or the citations Zasowski et al. 2013, Zasowski et al. 2017, Beaton et al. 2021, and Santana et al. 2021. The targeting strategies, and therefore the raw selection function, is split into different cohorts (defined by apparent magnitudes) and different color bins, which both vary between fields. The final product returns the selection fraction given the APOGEE field, apparent H magnitude, and (J-K) color of a given star.
Figure 1: The raw selection function for APOGEE, showing the selection fraction of observations for each APOGEE field (location on the sky). Near the plane of the disk, the selection fraction is generally low, simply because there are more stars in those fields. Near the galactic poles in the halo, the selection fraction is generally high.
The effective (or 3D) selection function is the ratio of the observed number of stars in APOGEE divided by the intrinsic number of stars in the Milky Way as a function of 3D Heliocentric position. This includes the effects of the raw selection function, with some additional ingredients of distances, isochrones (to assume some intrinsic distribution on the HR diagram) and a 3D Milky Way dust map. This is a more complex application of the selection function and will largely depend on your planned science case.
For example, you may want to split up the effective selection for different metallicities (as done in Mackereth & Bovy 2020), or by both metallicity and stellar ages (as in Mackereth et al. 2017). You may want smaller or larger sized bins than adopted by previous work, or a different range of distances sampled. Depending on how you define your data sample in APOGEE (cuts in surface gravity, kinematics, etc.), this should also be reflected in the effective selection function for best results.
Figure 2: The effective selection function for APOGEE, showing the selection fraction along different lines of sight. Generally, the selection fraction is larger close to the Sun, and decreases with distance.
A collection of additional resources related to the APOGEE selection function.
- The Stellar Population Structure of the Galactic Disk (Bovy et al. 2016a)
- On Galactic Density Modeling in the Presence of Dust Extinction (Bovy et al. 2016b)
- Effects of the selection function on metallicity trends in spectroscopic surveys of the Milky Way (Nandakumar et al. 2017)
- The age-metallicity structure of the Milky Way disc using APOGEE (Mackereth et al. 2017)
- Weighing the stellar constituents of the galactic halo with APOGEE red giant stars (Mackereth & Bovy 2020)
- The Milky Way tomography with APOGEE: intrinsic density distribution and structure of mono-abundance populations (Lian et al. 2022)
- https://github.com/jobovy/apogee: original
apogee
code for computing the raw selection function - https://github.com/jobovy/mwdust: repository containing several 3D dust maps for computing the effective selection function
- https://github.com/jmackereth/halo-mass: code for reproduction of Mackereth & Bovy 2020