Skip to content

Instantly share code, notes, and snippets.

@adrn
Last active June 15, 2026 17:23
Show Gist options
  • Select an option

  • Save adrn/a27a23850a943532c4e9355817384356 to your computer and use it in GitHub Desktop.

Select an option

Save adrn/a27a23850a943532c4e9355817384356 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": null,
"id": "98970e21",
"metadata": {},
"outputs": [],
"source": [
"import astropy.coordinates as coord\n",
"import astropy.table as at\n",
"import astropy.units as u\n",
"import matplotlib as mpl\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "d2a28730",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING: hdu= was not specified but multiple tables are present, reading in first available table (hdu=1) [astropy.io.fits.connect]\n"
]
}
],
"source": [
"apogee = at.Table.read(\"~/data/APOGEE/DR17/allStar-dr17-synspec_rev1.fits\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "a73d9aa9",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"['FILE',\n",
" 'APOGEE_ID',\n",
" 'TARGET_ID',\n",
" 'APSTAR_ID',\n",
" 'ASPCAP_ID',\n",
" 'TELESCOPE',\n",
" 'LOCATION_ID',\n",
" 'FIELD',\n",
" 'ALT_ID',\n",
" 'RA',\n",
" 'DEC',\n",
" 'GLON',\n",
" 'GLAT',\n",
" 'J',\n",
" 'J_ERR',\n",
" 'H',\n",
" 'H_ERR',\n",
" 'K',\n",
" 'K_ERR',\n",
" 'SRC_H',\n",
" 'WASH_M',\n",
" 'WASH_M_ERR',\n",
" 'WASH_T2',\n",
" 'WASH_T2_ERR',\n",
" 'DDO51',\n",
" 'DDO51_ERR',\n",
" 'IRAC_3_6',\n",
" 'IRAC_3_6_ERR',\n",
" 'IRAC_4_5',\n",
" 'IRAC_4_5_ERR',\n",
" 'IRAC_5_8',\n",
" 'IRAC_5_8_ERR',\n",
" 'IRAC_8_0',\n",
" 'IRAC_8_0_ERR',\n",
" 'WISE_4_5',\n",
" 'WISE_4_5_ERR',\n",
" 'TARG_4_5',\n",
" 'TARG_4_5_ERR',\n",
" 'WASH_DDO51_GIANT_FLAG',\n",
" 'WASH_DDO51_STAR_FLAG',\n",
" 'TARG_PMRA',\n",
" 'TARG_PMDEC',\n",
" 'TARG_PM_SRC',\n",
" 'AK_TARG',\n",
" 'AK_TARG_METHOD',\n",
" 'AK_WISE',\n",
" 'SFD_EBV',\n",
" 'APOGEE_TARGET1',\n",
" 'APOGEE_TARGET2',\n",
" 'APOGEE2_TARGET1',\n",
" 'APOGEE2_TARGET2',\n",
" 'APOGEE2_TARGET3',\n",
" 'APOGEE2_TARGET4',\n",
" 'TARGFLAGS',\n",
" 'SURVEY',\n",
" 'PROGRAMNAME',\n",
" 'NVISITS',\n",
" 'SNR',\n",
" 'SNREV',\n",
" 'STARFLAG',\n",
" 'STARFLAGS',\n",
" 'ANDFLAG',\n",
" 'ANDFLAGS',\n",
" 'VHELIO_AVG',\n",
" 'VSCATTER',\n",
" 'VERR',\n",
" 'RV_TEFF',\n",
" 'RV_LOGG',\n",
" 'RV_FEH',\n",
" 'RV_ALPHA',\n",
" 'RV_CARB',\n",
" 'RV_CHI2',\n",
" 'RV_CCFWHM',\n",
" 'RV_AUTOFWHM',\n",
" 'RV_FLAG',\n",
" 'N_COMPONENTS',\n",
" 'MEANFIB',\n",
" 'SIGFIB',\n",
" 'MIN_H',\n",
" 'MAX_H',\n",
" 'MIN_JK',\n",
" 'MAX_JK',\n",
" 'GAIAEDR3_SOURCE_ID',\n",
" 'GAIAEDR3_PARALLAX',\n",
" 'GAIAEDR3_PARALLAX_ERROR',\n",
" 'GAIAEDR3_PMRA',\n",
" 'GAIAEDR3_PMRA_ERROR',\n",
" 'GAIAEDR3_PMDEC',\n",
" 'GAIAEDR3_PMDEC_ERROR',\n",
" 'GAIAEDR3_PHOT_G_MEAN_MAG',\n",
" 'GAIAEDR3_PHOT_BP_MEAN_MAG',\n",
" 'GAIAEDR3_PHOT_RP_MEAN_MAG',\n",
" 'GAIAEDR3_DR2_RADIAL_VELOCITY',\n",
" 'GAIAEDR3_DR2_RADIAL_VELOCITY_ERROR',\n",
" 'GAIAEDR3_R_MED_GEO',\n",
" 'GAIAEDR3_R_LO_GEO',\n",
" 'GAIAEDR3_R_HI_GEO',\n",
" 'GAIAEDR3_R_MED_PHOTOGEO',\n",
" 'GAIAEDR3_R_LO_PHOTOGEO',\n",
" 'GAIAEDR3_R_HI_PHOTOGEO',\n",
" 'ASPCAP_GRID',\n",
" 'FPARAM_GRID',\n",
" 'CHI2_GRID',\n",
" 'FPARAM',\n",
" 'FPARAM_COV',\n",
" 'ASPCAP_CHI2',\n",
" 'PARAM',\n",
" 'PARAM_COV',\n",
" 'PARAMFLAG',\n",
" 'ASPCAPFLAG',\n",
" 'ASPCAPFLAGS',\n",
" 'FRAC_BADPIX',\n",
" 'FRAC_LOWSNR',\n",
" 'FRAC_SIGSKY',\n",
" 'FELEM',\n",
" 'FELEM_ERR',\n",
" 'X_H',\n",
" 'X_H_ERR',\n",
" 'X_M',\n",
" 'X_M_ERR',\n",
" 'ELEM_CHI2',\n",
" 'ELEMFRAC',\n",
" 'ELEMFLAG',\n",
" 'EXTRATARG',\n",
" 'MEMBERFLAG',\n",
" 'MEMBER',\n",
" 'X_H_SPEC',\n",
" 'X_M_SPEC',\n",
" 'TEFF',\n",
" 'TEFF_ERR',\n",
" 'LOGG',\n",
" 'LOGG_ERR',\n",
" 'M_H',\n",
" 'M_H_ERR',\n",
" 'ALPHA_M',\n",
" 'ALPHA_M_ERR',\n",
" 'VMICRO',\n",
" 'VMACRO',\n",
" 'VSINI',\n",
" 'TEFF_SPEC',\n",
" 'LOGG_SPEC',\n",
" 'C_FE',\n",
" 'C_FE_SPEC',\n",
" 'C_FE_ERR',\n",
" 'C_FE_FLAG',\n",
" 'CI_FE',\n",
" 'CI_FE_SPEC',\n",
" 'CI_FE_ERR',\n",
" 'CI_FE_FLAG',\n",
" 'N_FE',\n",
" 'N_FE_SPEC',\n",
" 'N_FE_ERR',\n",
" 'N_FE_FLAG',\n",
" 'O_FE',\n",
" 'O_FE_SPEC',\n",
" 'O_FE_ERR',\n",
" 'O_FE_FLAG',\n",
" 'NA_FE',\n",
" 'NA_FE_SPEC',\n",
" 'NA_FE_ERR',\n",
" 'NA_FE_FLAG',\n",
" 'MG_FE',\n",
" 'MG_FE_SPEC',\n",
" 'MG_FE_ERR',\n",
" 'MG_FE_FLAG',\n",
" 'AL_FE',\n",
" 'AL_FE_SPEC',\n",
" 'AL_FE_ERR',\n",
" 'AL_FE_FLAG',\n",
" 'SI_FE',\n",
" 'SI_FE_SPEC',\n",
" 'SI_FE_ERR',\n",
" 'SI_FE_FLAG',\n",
" 'P_FE',\n",
" 'P_FE_SPEC',\n",
" 'P_FE_ERR',\n",
" 'P_FE_FLAG',\n",
" 'S_FE',\n",
" 'S_FE_SPEC',\n",
" 'S_FE_ERR',\n",
" 'S_FE_FLAG',\n",
" 'K_FE',\n",
" 'K_FE_SPEC',\n",
" 'K_FE_ERR',\n",
" 'K_FE_FLAG',\n",
" 'CA_FE',\n",
" 'CA_FE_SPEC',\n",
" 'CA_FE_ERR',\n",
" 'CA_FE_FLAG',\n",
" 'TI_FE',\n",
" 'TI_FE_SPEC',\n",
" 'TI_FE_ERR',\n",
" 'TI_FE_FLAG',\n",
" 'TIII_FE',\n",
" 'TIII_FE_SPEC',\n",
" 'TIII_FE_ERR',\n",
" 'TIII_FE_FLAG',\n",
" 'V_FE',\n",
" 'V_FE_SPEC',\n",
" 'V_FE_ERR',\n",
" 'V_FE_FLAG',\n",
" 'CR_FE',\n",
" 'CR_FE_SPEC',\n",
" 'CR_FE_ERR',\n",
" 'CR_FE_FLAG',\n",
" 'MN_FE',\n",
" 'MN_FE_SPEC',\n",
" 'MN_FE_ERR',\n",
" 'MN_FE_FLAG',\n",
" 'FE_H',\n",
" 'FE_H_SPEC',\n",
" 'FE_H_ERR',\n",
" 'FE_H_FLAG',\n",
" 'CO_FE',\n",
" 'CO_FE_SPEC',\n",
" 'CO_FE_ERR',\n",
" 'CO_FE_FLAG',\n",
" 'NI_FE',\n",
" 'NI_FE_SPEC',\n",
" 'NI_FE_ERR',\n",
" 'NI_FE_FLAG',\n",
" 'CU_FE',\n",
" 'CU_FE_SPEC',\n",
" 'CU_FE_ERR',\n",
" 'CU_FE_FLAG',\n",
" 'CE_FE',\n",
" 'CE_FE_SPEC',\n",
" 'CE_FE_ERR',\n",
" 'CE_FE_FLAG',\n",
" 'YB_FE',\n",
" 'YB_FE_SPEC',\n",
" 'YB_FE_ERR',\n",
" 'YB_FE_FLAG',\n",
" 'VISIT_PK']"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"apogee.colnames # what column names are in the table"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "997a0c4c",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<Column name='APOGEE_ID' dtype='bytes30' length=4>\n",
"<table>\n",
"<tr><td>2M00000222+5625359</td></tr>\n",
"<tr><td>2M00000233+1452324</td></tr>\n",
"<tr><td>2M00000317+5821383</td></tr>\n",
"<tr><td>2M00000446+5854329</td></tr>\n",
"</table>"
],
"text/plain": [
"<Column name='APOGEE_ID' dtype='bytes30' length=4>\n",
"2M00000222+5625359\n",
"2M00000233+1452324\n",
"2M00000317+5821383\n",
"2M00000446+5854329"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"apogee[\"APOGEE_ID\"][10:14] # unique identifier per source in APOGEE"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "de98c73b",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"'VESTA'"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"apogee[\"APOGEE_ID\"][0] # an asteroid! for a spectrum of the sun. We can ignore"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "6faec5b5",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" RA \n",
"--------\n",
"0.000103\n",
"0.000832\n",
"0.001335\n",
" DEC \n",
"----------\n",
" 74.285408\n",
"-19.413851\n",
" 57.61953\n"
]
}
],
"source": [
"# Sky coordinates - Equatorial or ICRS\n",
"print(apogee[\"RA\"][1:4]) # right ascension in degrees\n",
"print(apogee[\"DEC\"][1:4]) # declination in degrees"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "547fd538",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" GLON \n",
"------------------\n",
"119.40180722431822\n",
" 63.39412165971994\n",
"116.06537094825998\n",
" GLAT \n",
"------------------\n",
"11.767413541486242\n",
"-75.90639677002453\n",
"-4.564767514078821\n"
]
}
],
"source": [
"# Sky coordinates - Galactic\n",
"print(apogee[\"GLON\"][1:4]) # Galactic Longitude\n",
"print(apogee[\"GLAT\"][1:4]) # Galactic Latitude"
]
},
{
"cell_type": "markdown",
"id": "a29fcf03",
"metadata": {},
"source": [
"Stellar parameters:\n",
"- `TEFF` = Effective temperature (in Kelvin)\n",
"- `LOGG` = log(surface gravity) in log(g cm/s^2)\n",
"- `M_H` = bulk metallicity $[M/H]$\n",
"- `ALPHA_M` = bulk alpha $[\\alpha/M]$\n",
"\n",
"All (or most) columns with measurements have an associated `_ERR` column that contains the uncertainty on that measurement."
]
},
{
"cell_type": "markdown",
"id": "82fa532b",
"metadata": {},
"source": [
"### Homework 1: Make a 2D histogram of all APOGEE stars - TEFF vs. LOGG\n",
"\n",
"This is related to the HR diagram, but is a spectroscopic version of it. It is sometimes called the Kiel diagram."
]
},
{
"cell_type": "markdown",
"id": "d4ae0bd7",
"metadata": {},
"source": [
"### Homework 2: Make a 2D histogram of all APOGEE stars - `M_H` vs. `ALPHA_M`"
]
},
{
"cell_type": "markdown",
"id": "9d6b7b07",
"metadata": {},
"source": [
"---"
]
},
{
"cell_type": "markdown",
"id": "68dd0ba0",
"metadata": {},
"source": [
"# 2026-06-15\n",
"\n",
"### Cross-match to El-Badry & Rix wide binary catalog"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "0577990b",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING: UnitsWarning: 'log(cm.s**-2)' did not parse as fits unit: 'log' is not a recognized function If this is meant to be a custom unit, define it with 'u.def_unit'. To have it recognized inside a file reader or other code, enable it with 'u.add_enabled_units'. For details, see https://docs.astropy.org/en/latest/units/combining_and_defining.html [astropy.units.core]\n",
"WARNING: UnitsWarning: 'log(cm.s**-2)' did not parse as fits unit: 'log' is not a recognized function If this is meant to be a custom unit, define it with 'u.def_unit'. To have it recognized inside a file reader or other code, enable it with 'u.add_enabled_units'. For details, see https://docs.astropy.org/en/latest/units/combining_and_defining.html [astropy.units.core]\n"
]
},
{
"data": {
"text/html": [
"<div><i>Table length=4</i>\n",
"<table id=\"table4389851104\" class=\"table-striped table-bordered table-condensed\">\n",
"<thead><tr><th>solution_id1</th><th>solution_id2</th><th>source_id1</th><th>source_id2</th><th>random_index1</th><th>random_index2</th><th>ref_epoch1</th><th>ref_epoch2</th><th>ra1</th><th>ra2</th><th>ra_error1</th><th>ra_error2</th><th>dec1</th><th>dec2</th><th>dec_error1</th><th>dec_error2</th><th>parallax1</th><th>parallax2</th><th>parallax_error1</th><th>parallax_error2</th><th>parallax_over_error1</th><th>parallax_over_error2</th><th>pm1</th><th>pm2</th><th>pmra1</th><th>pmra2</th><th>pmra_error1</th><th>pmra_error2</th><th>pmdec1</th><th>pmdec2</th><th>pmdec_error1</th><th>pmdec_error2</th><th>ra_dec_corr1</th><th>ra_dec_corr2</th><th>ra_parallax_corr1</th><th>ra_parallax_corr2</th><th>ra_pmra_corr1</th><th>ra_pmra_corr2</th><th>ra_pmdec_corr1</th><th>ra_pmdec_corr2</th><th>dec_parallax_corr1</th><th>dec_parallax_corr2</th><th>dec_pmra_corr1</th><th>dec_pmra_corr2</th><th>dec_pmdec_corr1</th><th>dec_pmdec_corr2</th><th>parallax_pmra_corr1</th><th>parallax_pmra_corr2</th><th>parallax_pmdec_corr1</th><th>parallax_pmdec_corr2</th><th>pmra_pmdec_corr1</th><th>pmra_pmdec_corr2</th><th>astrometric_n_obs_al1</th><th>astrometric_n_obs_al2</th><th>astrometric_n_obs_ac1</th><th>astrometric_n_obs_ac2</th><th>astrometric_n_good_obs_al1</th><th>astrometric_n_good_obs_al2</th><th>astrometric_n_bad_obs_al1</th><th>astrometric_n_bad_obs_al2</th><th>astrometric_gof_al1</th><th>astrometric_gof_al2</th><th>astrometric_chi2_al1</th><th>astrometric_chi2_al2</th><th>astrometric_excess_noise1</th><th>astrometric_excess_noise2</th><th>astrometric_excess_noise_sig1</th><th>astrometric_excess_noise_sig2</th><th>astrometric_params_solved1</th><th>astrometric_params_solved2</th><th>astrometric_primary_flag1</th><th>astrometric_primary_flag2</th><th>nu_eff_used_in_astrometry1</th><th>nu_eff_used_in_astrometry2</th><th>pseudocolour1</th><th>pseudocolour2</th><th>pseudocolour_error1</th><th>pseudocolour_error2</th><th>ra_pseudocolour_corr1</th><th>ra_pseudocolour_corr2</th><th>dec_pseudocolour_corr1</th><th>dec_pseudocolour_corr2</th><th>parallax_pseudocolour_corr1</th><th>parallax_pseudocolour_corr2</th><th>pmra_pseudocolour_corr1</th><th>pmra_pseudocolour_corr2</th><th>pmdec_pseudocolour_corr1</th><th>pmdec_pseudocolour_corr2</th><th>astrometric_matched_transits1</th><th>astrometric_matched_transits2</th><th>visibility_periods_used1</th><th>visibility_periods_used2</th><th>astrometric_sigma5d_max1</th><th>astrometric_sigma5d_max2</th><th>matched_transits1</th><th>matched_transits2</th><th>new_matched_transits1</th><th>new_matched_transits2</th><th>matched_transits_removed1</th><th>matched_transits_removed2</th><th>ipd_gof_harmonic_amplitude1</th><th>ipd_gof_harmonic_amplitude2</th><th>ipd_gof_harmonic_phase1</th><th>ipd_gof_harmonic_phase2</th><th>ipd_frac_multi_peak1</th><th>ipd_frac_multi_peak2</th><th>ipd_frac_odd_win1</th><th>ipd_frac_odd_win2</th><th>ruwe1</th><th>ruwe2</th><th>scan_direction_strength_k11</th><th>scan_direction_strength_k12</th><th>scan_direction_strength_k21</th><th>scan_direction_strength_k22</th><th>scan_direction_strength_k31</th><th>scan_direction_strength_k32</th><th>scan_direction_strength_k41</th><th>scan_direction_strength_k42</th><th>scan_direction_mean_k11</th><th>scan_direction_mean_k12</th><th>scan_direction_mean_k21</th><th>scan_direction_mean_k22</th><th>scan_direction_mean_k31</th><th>scan_direction_mean_k32</th><th>scan_direction_mean_k41</th><th>scan_direction_mean_k42</th><th>duplicated_source1</th><th>duplicated_source2</th><th>phot_g_n_obs1</th><th>phot_g_n_obs2</th><th>phot_g_mean_flux1</th><th>phot_g_mean_flux2</th><th>phot_g_mean_flux_error1</th><th>phot_g_mean_flux_error2</th><th>phot_g_mean_flux_over_error1</th><th>phot_g_mean_flux_over_error2</th><th>phot_g_mean_mag1</th><th>phot_g_mean_mag2</th><th>phot_bp_n_obs1</th><th>phot_bp_n_obs2</th><th>phot_bp_mean_flux1</th><th>phot_bp_mean_flux2</th><th>phot_bp_mean_flux_error1</th><th>phot_bp_mean_flux_error2</th><th>phot_bp_mean_flux_over_error1</th><th>phot_bp_mean_flux_over_error2</th><th>phot_bp_mean_mag1</th><th>phot_bp_mean_mag2</th><th>phot_rp_n_obs1</th><th>phot_rp_n_obs2</th><th>phot_rp_mean_flux1</th><th>phot_rp_mean_flux2</th><th>phot_rp_mean_flux_error1</th><th>phot_rp_mean_flux_error2</th><th>phot_rp_mean_flux_over_error1</th><th>phot_rp_mean_flux_over_error2</th><th>phot_rp_mean_mag1</th><th>phot_rp_mean_mag2</th><th>phot_bp_n_contaminated_transits1</th><th>phot_bp_n_contaminated_transits2</th><th>phot_bp_n_blended_transits1</th><th>phot_bp_n_blended_transits2</th><th>phot_rp_n_contaminated_transits1</th><th>phot_rp_n_contaminated_transits2</th><th>phot_rp_n_blended_transits1</th><th>phot_rp_n_blended_transits2</th><th>phot_proc_mode1</th><th>phot_proc_mode2</th><th>phot_bp_rp_excess_factor1</th><th>phot_bp_rp_excess_factor2</th><th>bp_rp1</th><th>bp_rp2</th><th>bp_g1</th><th>bp_g2</th><th>g_rp1</th><th>g_rp2</th><th>dr2_radial_velocity1</th><th>dr2_radial_velocity2</th><th>dr2_radial_velocity_error1</th><th>dr2_radial_velocity_error2</th><th>dr2_rv_nb_transits1</th><th>dr2_rv_nb_transits2</th><th>dr2_rv_template_teff1</th><th>dr2_rv_template_teff2</th><th>dr2_rv_template_logg1</th><th>dr2_rv_template_logg2</th><th>dr2_rv_template_fe_h1</th><th>dr2_rv_template_fe_h2</th><th>l1</th><th>l2</th><th>b1</th><th>b2</th><th>ecl_lon1</th><th>ecl_lon2</th><th>ecl_lat1</th><th>ecl_lat2</th><th>pairdistance</th><th>sep_AU</th><th>binary_type</th><th>Sigma18</th><th>R_chance_align</th><th>dr2_source_id1</th><th>dr2_source_id2</th><th>dr2_parallax1</th><th>dr2_parallax2</th><th>dr2_parallax_error1</th><th>dr2_parallax_error2</th><th>dr2_pmra1</th><th>dr2_pmra2</th><th>dr2_pmdec1</th><th>dr2_pmdec2</th><th>dr2_pmra_error1</th><th>dr2_pmra_error2</th><th>dr2_pmdec_error1</th><th>dr2_pmdec_error2</th><th>dr2_ruwe1</th><th>dr2_ruwe2</th><th>ID</th></tr></thead>\n",
"<thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th>yr</th><th>yr</th><th>deg</th><th>deg</th><th>mas</th><th>mas</th><th>deg</th><th>deg</th><th>mas</th><th>mas</th><th>mas</th><th>mas</th><th>mas</th><th>mas</th><th></th><th></th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th>mas / yr</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mas</th><th>mas</th><th></th><th></th><th></th><th></th><th></th><th></th><th>1 / um</th><th>1 / um</th><th>1 / um</th><th>1 / um</th><th>1 / um</th><th>1 / um</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mas</th><th>mas</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>deg</th><th>deg</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mag</th><th>mag</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mag</th><th>mag</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mag</th><th>mag</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mag</th><th>mag</th><th>mag</th><th>mag</th><th>mag</th><th>mag</th><th>km / s</th><th>km / s</th><th>km / s</th><th>km / s</th><th></th><th></th><th>K</th><th>K</th><th>log(cm.s**-2)</th><th>log(cm.s**-2)</th><th></th><th></th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th>deg</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead>\n",
"<thead><tr><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>bool</th><th>bool</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>bool</th><th>bool</th><th>int16</th><th>int16</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>int16</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>int16</th><th>int16</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>bytes4</th><th>float64</th><th>float64</th><th>int64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>int64</th></tr></thead>\n",
"<tr><td>1636042515805110273</td><td>1636042515805110273</td><td>4282339100022417152</td><td>4282339095720590720</td><td>1360492561</td><td>1271845445</td><td>2016.0</td><td>2016.0</td><td>282.80943113005753</td><td>282.8088579357593</td><td>0.012759874</td><td>0.28606945</td><td>5.692464376105104</td><td>5.688739805905678</td><td>0.012760569</td><td>0.24214806</td><td>2.466308813127245</td><td>3.129725645547363</td><td>0.015508463</td><td>0.28737023</td><td>159.02988</td><td>10.890918</td><td>23.447868</td><td>23.732042</td><td>-18.500443947672455</td><td>-18.279317121880627</td><td>0.015738279</td><td>0.3337244</td><td>-14.406112556212221</td><td>-15.135270582024575</td><td>0.01592136</td><td>0.30464804</td><td>0.0135885</td><td>0.31395996</td><td>-0.07543358</td><td>-0.39887094</td><td>0.0022587231</td><td>-0.02041166</td><td>-0.037667114</td><td>-0.11068297</td><td>-0.22045594</td><td>-0.012907138</td><td>0.042067286</td><td>-0.14564936</td><td>-0.01791944</td><td>0.17442913</td><td>-0.24544272</td><td>-0.1367347</td><td>0.18651143</td><td>-0.22729401</td><td>0.19416681</td><td>0.47724915</td><td>238</td><td>201</td><td>0</td><td>0</td><td>238</td><td>201</td><td>0</td><td>0</td><td>2.6347468</td><td>-0.69246155</td><td>294.6977</td><td>215.23299</td><td>0.023164036</td><td>0.0</td><td>0.32207066</td><td>0.0</td><td>31</td><td>95</td><td>False</td><td>False</td><td>1.4765114</td><td>1e+20</td><td>1e+20</td><td>1.1112912</td><td>1e+20</td><td>0.07932555</td><td>1e+20</td><td>0.11098927</td><td>1e+20</td><td>-0.18573934</td><td>1e+20</td><td>0.23411377</td><td>1e+20</td><td>-0.24445736</td><td>1e+20</td><td>-0.17332205</td><td>27</td><td>23</td><td>16</td><td>13</td><td>0.02391526</td><td>0.5400338</td><td>27</td><td>23</td><td>7</td><td>6</td><td>0</td><td>0</td><td>0.061279934</td><td>0.04265072</td><td>83.45006</td><td>168.73517</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1.123015</td><td>0.96345276</td><td>0.27786088</td><td>0.26865715</td><td>0.12818311</td><td>0.14001864</td><td>0.20724767</td><td>0.22512177</td><td>0.7795465</td><td>0.7824155</td><td>-98.74266</td><td>-111.90673</td><td>-38.26392</td><td>-39.064552</td><td>-35.498913</td><td>-27.067425</td><td>40.009457</td><td>41.823524</td><td>False</td><td>False</td><td>233</td><td>198</td><td>67666.58939786452</td><td>348.18941563607075</td><td>12.298133</td><td>1.4432361</td><td>5502.184</td><td>241.25603</td><td>13.611431</td><td>19.332829</td><td>25</td><td>21</td><td>30478.356684246726</td><td>55.88765065725574</td><td>27.55253</td><td>6.282352</td><td>1106.1908</td><td>8.895976</td><td>14.128563</td><td>20.970253</td><td>25</td><td>21</td><td>53484.867591955146</td><td>453.3957013079742</td><td>27.917492</td><td>5.7025685</td><td>1915.8192</td><td>79.50728</td><td>12.927319</td><td>18.106703</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1.2408372</td><td>1.4626617</td><td>1.2012444</td><td>2.8635502</td><td>0.5171318</td><td>1.6374245</td><td>0.68411255</td><td>1.2261257</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>0</td><td>0</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>37.97896667366811</td><td>37.975385523516714</td><td>2.636656992696798</td><td>2.6354736772757086</td><td>284.5371627955313</td><td>284.53609195487655</td><td>28.48760913726913</td><td>28.483960591336587</td><td>0.0037679894569219168</td><td>5500.025776463481</td><td>MSMS</td><td>1599.825487959732</td><td>0.15309618110360945</td><td>4282339100022417152</td><td>4282339095720590720</td><td>2.466708691330006</td><td>3.0406992428300974</td><td>0.02155052097944453</td><td>0.3581828570273455</td><td>-18.392622882035198</td><td>-18.703935951842688</td><td>-14.236095644480766</td><td>-14.871651131577826</td><td>0.03657619222459426</td><td>0.7373521857596478</td><td>0.03823012027556794</td><td>0.7127700515631703</td><td>1.0173349</td><td>1.0269454</td><td>0</td></tr>\n",
"<tr><td>1636042515805110273</td><td>1636042515805110273</td><td>4089436931798712576</td><td>4089437344115586432</td><td>761856756</td><td>699795040</td><td>2016.0</td><td>2016.0</td><td>276.31365516286564</td><td>276.32563147342296</td><td>0.0833635</td><td>0.21712646</td><td>-23.55099339732246</td><td>-23.53250430327995</td><td>0.07513903</td><td>0.19658527</td><td>2.4230962227012838</td><td>2.92015611211368</td><td>0.089352295</td><td>0.23494595</td><td>27.118456</td><td>12.429055</td><td>16.392246</td><td>16.321041</td><td>-4.344694168896968</td><td>-3.844595829572494</td><td>0.10976257</td><td>0.2791153</td><td>-15.805990832000269</td><td>-15.86176044805748</td><td>0.0855414</td><td>0.22076848</td><td>0.09113432</td><td>0.1280763</td><td>-0.058927197</td><td>-0.062295474</td><td>-0.30838892</td><td>-0.27625978</td><td>-0.025096338</td><td>-0.08854862</td><td>-0.33874184</td><td>-0.29305327</td><td>0.008507379</td><td>-0.03977013</td><td>-0.4753635</td><td>-0.44567153</td><td>-0.03957939</td><td>-0.09652548</td><td>0.1924429</td><td>0.22701822</td><td>0.037752945</td><td>0.036217965</td><td>221</td><td>205</td><td>0</td><td>0</td><td>221</td><td>187</td><td>0</td><td>18</td><td>9.096802</td><td>10.193055</td><td>485.0724</td><td>541.3772</td><td>0.51799184</td><td>1.5900779</td><td>10.650072</td><td>35.430706</td><td>31</td><td>95</td><td>False</td><td>False</td><td>1.3405141</td><td>1e+20</td><td>1e+20</td><td>1.2622169</td><td>1e+20</td><td>0.045261778</td><td>1e+20</td><td>-0.022542885</td><td>1e+20</td><td>0.08699399</td><td>1e+20</td><td>-0.20559515</td><td>1e+20</td><td>0.01587107</td><td>1e+20</td><td>-0.042962488</td><td>25</td><td>22</td><td>16</td><td>14</td><td>0.15480213</td><td>0.39342427</td><td>26</td><td>26</td><td>11</td><td>11</td><td>0</td><td>0</td><td>0.03635774</td><td>0.10744859</td><td>1.9436731</td><td>162.18494</td><td>0</td><td>3</td><td>0</td><td>3</td><td>1.4664384</td><td>1.5789028</td><td>0.21004309</td><td>0.19065337</td><td>0.22895774</td><td>0.18879564</td><td>0.3447228</td><td>0.24035944</td><td>0.8034685</td><td>0.8496756</td><td>-118.56969</td><td>-137.35849</td><td>-8.152375</td><td>-7.749687</td><td>-31.397596</td><td>-33.579506</td><td>42.427635</td><td>41.564384</td><td>False</td><td>False</td><td>228</td><td>213</td><td>4495.049815890725</td><td>1791.380810424936</td><td>6.4739532</td><td>7.7872605</td><td>694.3284</td><td>230.03992</td><td>16.55553</td><td>17.554398</td><td>21</td><td>21</td><td>1363.7568860350852</td><td>546.6146365890555</td><td>8.314398</td><td>8.711877</td><td>164.02353</td><td>62.743614</td><td>17.5017</td><td>18.494339</td><td>21</td><td>22</td><td>5290.561792175467</td><td>2639.088121765906</td><td>13.429629</td><td>24.940964</td><td>393.94696</td><td>105.8134</td><td>15.439141</td><td>16.194262</td><td>0</td><td>0</td><td>21</td><td>21</td><td>0</td><td>0</td><td>21</td><td>22</td><td>0</td><td>0</td><td>1.480366</td><td>1.7783504</td><td>2.0625582</td><td>2.3000774</td><td>0.9461689</td><td>0.9399414</td><td>1.1163893</td><td>1.360136</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>0</td><td>0</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>9.048744643468892</td><td>9.070298449174286</td><td>-5.116972437022136</td><td>-5.118207820054266</td><td>275.78594499148744</td><td>275.7977236291275</td><td>-0.23843192962063267</td><td>-0.22044131436892278</td><td>0.021503402427766518</td><td>31947.65771771943</td><td>MSMS</td><td>3247.7157687332165</td><td>1.1065654664144589</td><td>4089436931798712576</td><td>4089437344115586432</td><td>2.7156720716700447</td><td>2.973328646283168</td><td>0.1357684499083049</td><td>0.18723981922559285</td><td>-4.461936422952587</td><td>-2.6124374196703797</td><td>-15.582909772985671</td><td>-17.414057477749957</td><td>0.2319180727363056</td><td>0.3366228285711499</td><td>0.19365308742074325</td><td>0.288403784877605</td><td>1.5857931</td><td>1.2596724</td><td>1</td></tr>\n",
"<tr><td>1636042515805110273</td><td>1636042515805110273</td><td>4282722760856124416</td><td>4282722790916299648</td><td>291010844</td><td>652726121</td><td>2016.0</td><td>2016.0</td><td>283.2152167318123</td><td>283.19616572286793</td><td>0.07057537</td><td>0.18889558</td><td>5.949548004855735</td><td>5.946408290921299</td><td>0.067139946</td><td>0.18544991</td><td>1.043858388313062</td><td>1.129170131567232</td><td>0.074837394</td><td>0.20538272</td><td>13.948353</td><td>5.497883</td><td>5.834742</td><td>5.547782</td><td>-2.044045573754162</td><td>-2.223293512377718</td><td>0.095851004</td><td>0.24849989</td><td>-5.464987837879764</td><td>-5.082799378502095</td><td>0.0833044</td><td>0.2140873</td><td>-0.016794901</td><td>0.005652371</td><td>-0.06883118</td><td>0.034949817</td><td>-0.3903171</td><td>-0.34274393</td><td>0.09097812</td><td>0.028798917</td><td>-0.19159089</td><td>-0.25179547</td><td>0.080729</td><td>0.044025928</td><td>-0.37003598</td><td>-0.37776563</td><td>0.028635109</td><td>-0.059663367</td><td>-0.1395068</td><td>-0.059960846</td><td>-0.13877364</td><td>-0.08365903</td><td>234</td><td>214</td><td>0</td><td>0</td><td>234</td><td>214</td><td>0</td><td>0</td><td>0.4069359</td><td>0.3564547</td><td>248.60141</td><td>235.54224</td><td>0.08711514</td><td>0.0</td><td>0.19021663</td><td>0.0</td><td>31</td><td>95</td><td>False</td><td>False</td><td>1.3148338</td><td>1e+20</td><td>1e+20</td><td>1.2797575</td><td>1e+20</td><td>0.0564599</td><td>1e+20</td><td>-0.21444717</td><td>1e+20</td><td>0.04149634</td><td>1e+20</td><td>0.031542715</td><td>1e+20</td><td>-0.09258738</td><td>1e+20</td><td>-0.016667336</td><td>27</td><td>24</td><td>15</td><td>15</td><td>0.139839</td><td>0.3553256</td><td>27</td><td>24</td><td>11</td><td>10</td><td>0</td><td>0</td><td>0.017543955</td><td>0.012991707</td><td>96.186745</td><td>162.54947</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1.0176105</td><td>1.015916</td><td>0.26805305</td><td>0.26782247</td><td>0.12179362</td><td>0.16312487</td><td>0.14784658</td><td>0.11431878</td><td>0.74657273</td><td>0.7420139</td><td>-130.41455</td><td>-125.53641</td><td>16.231266</td><td>27.461136</td><td>-28.650194</td><td>-30.53759</td><td>39.35053</td><td>38.98983</td><td>False</td><td>False</td><td>230</td><td>213</td><td>2447.0742595051647</td><td>548.7967944372285</td><td>2.3995965</td><td>1.4743094</td><td>1019.78577</td><td>372.2399</td><td>17.21575</td><td>18.838839</td><td>23</td><td>21</td><td>612.1657759089306</td><td>100.80162749969595</td><td>8.513242</td><td>6.786092</td><td>71.90748</td><td>14.85415</td><td>18.37137</td><td>20.329874</td><td>22</td><td>23</td><td>2705.6832777904297</td><td>674.9737645698359</td><td>14.810265</td><td>8.481548</td><td>182.68973</td><td>79.581436</td><td>16.167204</td><td>17.674679</td><td>0</td><td>0</td><td>2</td><td>1</td><td>1</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>1.3558432</td><td>1.4135932</td><td>2.2041664</td><td>2.6551952</td><td>1.1556206</td><td>1.4910355</td><td>1.0485458</td><td>1.1641598</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>0</td><td>0</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>38.391788615675836</td><td>38.38037008158468</td><td>2.393826957497478</td><td>2.4092784109100647</td><td>285.0244060434098</td><td>285.00254873799764</td><td>28.70217240185595</td><td>28.701012648899912</td><td>0.019206806115865043</td><td>66239.35084610071</td><td>MSMS</td><td>1632.6114062366623</td><td>0.9797521845707992</td><td>4282722760856124416</td><td>4282722790916299648</td><td>1.1612547920623557</td><td>0.9349971137078245</td><td>0.107813843419762</td><td>0.27852800764198726</td><td>-1.99265369953549</td><td>-3.6482952573971734</td><td>-5.398722276484872</td><td>-3.9345624494829297</td><td>0.231762645997847</td><td>0.6310041696938047</td><td>0.20443003304744498</td><td>0.5958823499872595</td><td>1.0708025</td><td>0.9836256</td><td>2</td></tr>\n",
"<tr><td>1636042515805110273</td><td>1636042515805110273</td><td>4128188325747199104</td><td>4128188325747198080</td><td>281631074</td><td>1344115143</td><td>2016.0</td><td>2016.0</td><td>256.87034730707893</td><td>256.8701664814668</td><td>0.015697097</td><td>0.025262747</td><td>-20.23579945116703</td><td>-20.23719657986973</td><td>0.010571578</td><td>0.016789839</td><td>1.6715541895959316</td><td>1.6474244682595314</td><td>0.016625766</td><td>0.026787495</td><td>100.53998</td><td>61.499756</td><td>1.507807</td><td>1.3825816</td><td>-0.5260189837622264</td><td>-0.6637311157992818</td><td>0.021589695</td><td>0.034472495</td><td>-1.4130767700283122</td><td>-1.2128449819335219</td><td>0.0147116585</td><td>0.023518544</td><td>0.07724272</td><td>0.08398727</td><td>0.021930734</td><td>0.04200924</td><td>-0.13659847</td><td>-0.11972336</td><td>-0.25750425</td><td>-0.2896075</td><td>-0.18570586</td><td>-0.21287854</td><td>-0.26072368</td><td>-0.26702547</td><td>-0.2816742</td><td>-0.23805536</td><td>0.020924134</td><td>-0.043921072</td><td>0.058815546</td><td>0.045249768</td><td>0.14092436</td><td>0.111938976</td><td>273</td><td>255</td><td>273</td><td>0</td><td>270</td><td>254</td><td>3</td><td>1</td><td>-2.3143704</td><td>-1.1270227</td><td>526.63513</td><td>228.53296</td><td>0.082437895</td><td>0.0</td><td>7.8279963</td><td>0.0</td><td>31</td><td>31</td><td>False</td><td>False</td><td>1.6391888</td><td>1.4466388</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>31</td><td>30</td><td>16</td><td>16</td><td>0.030412057</td><td>0.048369117</td><td>31</td><td>30</td><td>13</td><td>12</td><td>0</td><td>0</td><td>0.026347542</td><td>0.018985057</td><td>158.05692</td><td>125.00534</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.89995915</td><td>0.94860876</td><td>1e+20</td><td>0.111483134</td><td>1e+20</td><td>0.42406237</td><td>1e+20</td><td>0.2005713</td><td>1e+20</td><td>0.43129733</td><td>1e+20</td><td>-92.99507</td><td>1e+20</td><td>-0.064303614</td><td>1e+20</td><td>-30.90789</td><td>1e+20</td><td>-41.977512</td><td>False</td><td>False</td><td>264</td><td>250</td><td>980659.1912538076</td><td>23215.96307283315</td><td>294.92294</td><td>13.060419</td><td>3325.137</td><td>1777.5818</td><td>10.708571</td><td>14.772901</td><td>30</td><td>28</td><td>610902.9329543022</td><td>9989.856380440262</td><td>263.2021</td><td>76.26344</td><td>2321.0413</td><td>130.99141</td><td>10.873611</td><td>15.339644</td><td>28</td><td>24</td><td>542188.2571678016</td><td>20083.66415302557</td><td>255.86897</td><td>92.109055</td><td>2119.0073</td><td>218.04224</td><td>10.41252</td><td>13.990788</td><td>0</td><td>0</td><td>0</td><td>17</td><td>0</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>1.1758327</td><td>1.2953812</td><td>0.46109104</td><td>1.348856</td><td>0.16504002</td><td>0.56674385</td><td>0.29605103</td><td>0.7821121</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>0</td><td>0</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>1e+20</td><td>2.5069057662082868</td><td>2.5056376331611387</td><td>12.017129623246857</td><td>12.016464565884577</td><td>257.6805819951459</td><td>257.6805393572747</td><td>2.643737119052615</td><td>2.642330370928953</td><td>0.0014073927687894647</td><td>3031.0796976715646</td><td>MSMS</td><td>1733.833950043108</td><td>0.00018070210517276975</td><td>4128188325747199104</td><td>4128188325747198080</td><td>1.6694597921324998</td><td>1.6354206782265281</td><td>0.04536595779583257</td><td>0.03732093869415257</td><td>-0.6805072008255597</td><td>-0.6873758177346445</td><td>-1.268069533191858</td><td>-1.1844411191277335</td><td>0.09526324583714986</td><td>0.0793840482073535</td><td>0.05703277274203961</td><td>0.046327022794392145</td><td>0.927876</td><td>1.0243161</td><td>3</td></tr>\n",
"</table></div>"
],
"text/plain": [
"<Table length=4>\n",
" solution_id1 solution_id2 source_id1 ... dr2_ruwe2 ID \n",
" ... \n",
" int64 int64 int64 ... float64 int64\n",
"------------------- ------------------- ------------------- ... --------- -----\n",
"1636042515805110273 1636042515805110273 4282339100022417152 ... 1.0269454 0\n",
"1636042515805110273 1636042515805110273 4089436931798712576 ... 1.2596724 1\n",
"1636042515805110273 1636042515805110273 4282722760856124416 ... 0.9836256 2\n",
"1636042515805110273 1636042515805110273 4128188325747199104 ... 1.0243161 3"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"binaries = at.Table.read(\"~/data/Misc/ElBadry2021/all_columns_catalog.fits\")\n",
"binaries[\"ID\"] = np.arange(len(binaries), dtype=int)\n",
"binaries[:4]"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "219331af",
"metadata": {},
"outputs": [],
"source": [
"star1 = at.unique(\n",
" at.join(\n",
" apogee,\n",
" binaries,\n",
" keys_left=\"GAIAEDR3_SOURCE_ID\",\n",
" keys_right=\"source_id1\",\n",
" join_type=\"inner\",\n",
" ),\n",
" keys=\"APOGEE_ID\",\n",
")\n",
"star2 = at.unique(\n",
" at.join(\n",
" apogee,\n",
" binaries,\n",
" keys_left=\"GAIAEDR3_SOURCE_ID\",\n",
" keys_right=\"source_id2\",\n",
" join_type=\"inner\",\n",
" ),\n",
" keys=\"APOGEE_ID\",\n",
")"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "a1a8b736",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(391, 391)"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mask1 = np.isin(star1[\"ID\"], star2[\"ID\"])\n",
"star1 = star1[mask1]\n",
"star1.sort(\"ID\")\n",
"\n",
"mask2 = np.isin(star2[\"ID\"], star1[\"ID\"])\n",
"star2 = star2[mask2]\n",
"star2.sort(\"ID\")\n",
"\n",
"len(star1), len(star2)"
]
},
{
"cell_type": "code",
"execution_count": 21,
"id": "5aa1c754",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"&lt;MaskedColumn name=&apos;MG_FE&apos; dtype=&apos;float32&apos; length=391&gt;\n",
"<table>\n",
"<tr><td>-0.08541149</td></tr>\n",
"<tr><td>0.2592385</td></tr>\n",
"<tr><td>0.03301151</td></tr>\n",
"<tr><td>0.120337486</td></tr>\n",
"<tr><td>-0.16841048</td></tr>\n",
"<tr><td>0.041224018</td></tr>\n",
"<tr><td>0.0339365</td></tr>\n",
"<tr><td>0.108945504</td></tr>\n",
"<tr><td>-0.02710449</td></tr>\n",
"<tr><td>-0.0375815</td></tr>\n",
"<tr><td>0.0827294</td></tr>\n",
"<tr><td>-0.050447498</td></tr>\n",
"<tr><td>...</td></tr>\n",
"<tr><td>-0.08567151</td></tr>\n",
"<tr><td>0.06441358</td></tr>\n",
"<tr><td>-0.114105485</td></tr>\n",
"<tr><td>-0.2237685</td></tr>\n",
"<tr><td>0.027359223</td></tr>\n",
"<tr><td>0.0927175</td></tr>\n",
"<tr><td>--</td></tr>\n",
"<tr><td>-0.03652849</td></tr>\n",
"<tr><td>0.035951495</td></tr>\n",
"<tr><td>--</td></tr>\n",
"<tr><td>-0.031656496</td></tr>\n",
"</table>"
],
"text/plain": [
"<MaskedColumn name='MG_FE' dtype='float32' length=391>\n",
" -0.08541149\n",
" 0.2592385\n",
" 0.03301151\n",
" 0.120337486\n",
" -0.16841048\n",
" 0.041224018\n",
" 0.0339365\n",
" 0.108945504\n",
" -0.02710449\n",
" -0.0375815\n",
" 0.0827294\n",
"-0.050447498\n",
" ...\n",
" -0.08567151\n",
" 0.06441358\n",
"-0.114105485\n",
" -0.2237685\n",
" 0.027359223\n",
" 0.0927175\n",
" --\n",
" -0.03652849\n",
" 0.035951495\n",
" --\n",
"-0.031656496"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"star1[\"MG_FE\"]"
]
},
{
"cell_type": "markdown",
"id": "bb5d87b7",
"metadata": {},
"source": [
"### TODO:\n",
"- Plot the stellar parameters (TEFF and LOGG) for both stars. I like to make a scatter plot where star1 is colored red and star2 is blue, and they are connected by lines. Maybe this is too many to plot on one figure, so maybe we want to split the sample into a few panels?\n",
"- Pick a few pairs of stars where both stars have similar stellar parameters (TEFF and LOGG) - let's say TEFF within 200 K and LOGG within 0.3. Make a plot like the one in Oh et al - x axis are elements, and y axis are the differences in those elements between the two stars."
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "ff4750f4",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "demo-notebooks (3.12.10)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.10"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment