{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Band-aid for Gaia refcat errors for Jointcal on HSC RC2\n", "\n", "\n", "January 4-6 2021\n", "\n", "After switching to the Gaia refcat on RC2 to we noticed that the solution walks away from the reference catalog. Hypothesis: The refcat errors are too large and thus the refcat is being underwieghted. \n", "\n", "See: DM-28173\n", "\n", "Symptoms are the stellar locus width increased and the measured positions - ref cat positions increased. This didn't appear in the relative astrometry metric because the shifts are coherent.\n", "\n", "\n", "## Intro to metrics computed\n", "\n", "This is basically the same notebook as https://github.com/yalsayyad/dm_notebooks/blob/master/calibration/run21i_validate_drp.ipynb\n", "but with HSC RC2 with the Gaia refcat." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We care about both relative astrometry and absolute astrometry.\n", "\n", "### Relative astrometry\n", "\n", "Validate_drp computes astrometric and photometric repeatabliity metrics per the SRD. See ls.st/LPM-17 for more info.\n", "The astrometric repeatability metrics `AM1`, `AM2`, and `AM3` is the median RMS of distances between star pairs. \n", "\n", "The ELI5 version courtesy of John Parejko is: \"Take a pair of stars ~5' apart on the sky. Any time that pair is seen in a single visit, record their separation. Compute the RMS of all those separations...that RMS is one datapoint. Do that for all such pairs separated by ~5': that is the histogram shown\" below. The median of all those RMS's is the AM1. \n", "\n", "The 1, 2, 3 correspond to pairs 5', 20' and 200' apart to represent scales of an individual sensor, a raft, and\n", "the camera. \n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The other validate_drp astrometry plot prints the naive distribution of the RMS variation in the position of stars. This is the metric that one thinks of first when they hear \"repeatabiliy\" and is the analog of PA1. \n", "\n", "Code that computes AMx can be found in:\n", "* `https://github.com/lsst/validate_drp/blob/tickets/DM-18638/python/lsst/validate/drp/calcsrd/amx.py`\n", "\n", "\n", "### Absolute astrometry\n", "\n", "We don't have a squash metric for this yet, but `visitAnalysis` in the pipe_analysis package computes this.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Data (links to QA plots and configs)\n", "\n", "Ran `visitAnalysis` and `validateDRP` on the following 2 configurations.\n", "\n", "\n", "---\n", "| link to QA plots | title | Description | Why | \n", "|-------------------------------------------------------------------------|-------|----------------------------|------|\n", "| https://lsst.ncsa.illinois.edu/~emorgan2/w_2020_50_qaplots | w50 | w50 (Gaia) | | \n", "| https://lsst.ncsa.illinois.edu/~yusra/RC_QA/w_2020_50_ps1 | w50_ps1 | w50 w/ jointcal and obs_subaru on w49 (PS1) | Control. To compare with just the refcat change. Effect of reverting to PS1 | \n", "| https://lsst.ncsa.illinois.edu/~yusra/RC_QA/w_2020_50_refErr10 | w50_refErr10 | w50 `astrometryReferenceErr=10.0` (Gaia) | bandaid | \n", "| https://lsst.ncsa.illinois.edu/~yusra/RC_QA/w_2020_50_fixErr | w50_fixErr | w01 + DM-28233 (Gaia) | Proposed fix | \n", "\n", "\n", "Expectation: One change that came with the PS1->Gaia is that Gaia comes with per-object position errors. We know from https://github.com/yalsayyad/dm_notebooks/blob/master/calibration/run21i_validate_drp.ipynb that jointcal the reference errors affect how strictly it tries to adhere to the reference catalog. Too roomy and the fit walks away from the reference catalog. \n", "\n", "My hypothesis is that this is what happened here and I want to find a way to dial it back. \n", "\n", "Update: Eli found refcat errors are in the wrong units and posted fix on DM-28233. Tested as: `w50_fixErr` " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import sys\n", "import os\n", "import numpy as np\n", "import copy\n", "import pickle\n", "import json\n", "import matplotlib.pyplot as plt\n", "import scipy.integrate as integrate\n", "import seaborn as sns\n", "import pandas as pd\n", "%matplotlib inline \n", "sns.set_palette('Dark2')\n", "sns.set(context='talk',\n", " style='whitegrid',\n", " palette='deep',\n", " font='sans-serif',\n", " font_scale=1,\n", " color_codes=True,\n", " rc={'text.usetex': False})" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9813/HSC-G/matchedVisit_HSC-G.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9813/HSC-R/matchedVisit_HSC-R.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9813/HSC-I/matchedVisit_HSC-I.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9813/HSC-Z/matchedVisit_HSC-Z.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9813/HSC-Y/matchedVisit_HSC-Y.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9697/HSC-G/matchedVisit_HSC-G.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9697/HSC-R/matchedVisit_HSC-R.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9697/HSC-I/matchedVisit_HSC-I.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9697/HSC-Z/matchedVisit_HSC-Z.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9697/HSC-Y/matchedVisit_HSC-Y.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9615/HSC-G/matchedVisit_HSC-G.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9615/HSC-R/matchedVisit_HSC-R.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9615/HSC-I/matchedVisit_HSC-I.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9615/HSC-Z/matchedVisit_HSC-Z.json\n", "loaded\n", "/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics/9615/HSC-Y/matchedVisit_HSC-Y.json\n", "loaded\n" ] } ], "source": [ "DIR_DICT_DC2 = {\n", " #\"w50_refErr10\": \"/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_refErr10/plots/validateDrp/matchedVisitMetrics\",\n", " #\"w50_ps1\": \"/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_ps1/plots/validateDrp/matchedVisitMetrics\",\n", " #\"w50\": \"/datasets/hsc/repo/rerun/RC/w_2020_50/DM-28140/validateDrp/matchedVisitMetrics\",\n", " \"w50_fixErr\": \"/datasets/hsc/repo/rerun/private/yusra/RC2/w_2020_50_fixErr/plots/validateDrp/matchedVisitMetrics\",\n", "}\n", "\n", "rowList = []\n", "\n", "for TRACT in (9813, 9697, 9615):\n", " for FILTER in ['HSC-%s'%(band) for band in 'GRIZY']:\n", " for dirName, DIR in DIR_DICT_DC2.items():\n", " dictDF = {}\n", " fname = \"matchedVisit_%s\" % (FILTER)\n", " path = os.path.join(DIR, str(TRACT), FILTER, fname + '.json')\n", " print(path)\n", " try:\n", " with open(path) as json_file:\n", " data = json.load(json_file)\n", " except:\n", " continue\n", " print(\"loaded\")\n", " values = [d['value'] for d in data['measurements']]\n", " metrics = [d['metric'].replace('validate_drp.', '') for d in data['measurements']]\n", " for metric, value in zip(metrics, values):\n", " dictDF[metric] = value\n", " dictDF['tract'] = TRACT\n", " dictDF['filter'] = FILTER\n", " dictDF['config'] = dirName\n", " rowList.append(dictDF)\n", "\n", "dffixErr = pd.DataFrame(rowList)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "unitsDict = {}\n", "for name, unit, disc in [(d['name'], d['unit'], d['description']) for d in data['metrics']]:\n", " unitsDict[name.replace('validate_drp.', '')] = (unit, disc)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'PA1': ('mmag',\n", " 'The maximum rms of the unresolved source magnitude distribution around the mean value.'),\n", " 'PF1_minimum_gri': ('%',\n", " 'The maximum fraction of magnitudes deviating by more than PA2_minimum_gri from mean.'),\n", " 'PF1_design_gri': ('%',\n", " 'The maximum fraction of magnitudes deviating by more than PA2_design_gri from mean.'),\n", " 'PF1_stretch_gri': ('%',\n", " 'The maximum fraction of magnitudes deviating by more than PA2_stretch_gri from mean.'),\n", " 'PF1_minimum_uzy': ('%',\n", " 'The maximum fraction of magnitudes deviating by more than PA2_minimum_uzy from mean.'),\n", " 'PF1_design_uzy': ('%',\n", " 'The maximum fraction of magnitudes deviating by more than PA2_design_uzy from mean.'),\n", " 'PF1_stretch_uzy': ('%',\n", " 'The maximum fraction of magnitudes deviating by more than PA2_stretch_uzy from mean.'),\n", " 'PA2_minimum_gri': ('mmag',\n", " 'At most PF1_minimum_gri % of magnitudes may deviate by more than this magnitude from the mean.'),\n", " 'PA2_design_gri': ('mmag',\n", " 'At most PF1_design_gri % of magnitudes may deviate by more than this magnitude from the mean.'),\n", " 'PA2_stretch_gri': ('mmag',\n", " 'At most PF1_stretch_gri % of magnitudes may deviate by more than this magnitude from the mean.'),\n", " 'PA2_minimum_uzy': ('mmag',\n", " 'At most PF1_minimum_uzy % of magnitudes may deviate by more than this magnitude from the mean.'),\n", " 'PA2_design_uzy': ('mmag',\n", " 'At most PF1_design_uzy % of magnitudes may deviate by more than this magnitude from the mean.'),\n", " 'PA2_stretch_uzy': ('mmag',\n", " 'At most PF1_stretch_uzy % of magnitudes may deviate by more than this magnitude from the mean.'),\n", " 'AM1': ('marcsec',\n", " 'The maximum rms of the astrometric distance distribution for stellar pairs with separations of D=5 arcmin (repeatability).'),\n", " 'AM2': ('marcsec',\n", " 'The maximum rms of the astrometric distance distribution for stellar pairs with separations of D=20 arcmin (repeatability).'),\n", " 'AM3': ('marcsec',\n", " 'The maximum rms of the astrometric distance distribution for stellar pairs with separations of D=200 arcmin (repeatability).'),\n", " 'AF1_minimum': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD1 milliarcsec (see AM1).'),\n", " 'AF1_design': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD1 milliarcsec (see AM1).'),\n", " 'AF1_stretch': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD1 milliarcsec (see AM1).'),\n", " 'AF2_minimum': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD2 milliarcsec (see AM2).'),\n", " 'AF2_design': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD2 milliarcsec (see AM2).'),\n", " 'AF2_stretch': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD2 milliarcsec (see AM2).'),\n", " 'AF3_minimum': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD3 milliarcsec (see AM3).'),\n", " 'AF3_design': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD3 milliarcsec (see AM3).'),\n", " 'AF3_stretch': ('%',\n", " 'The maximum fraction of astrometric distances which deviate by more than AD3 milliarcsec (see AM3).'),\n", " 'AD1_minimum': ('marcsec',\n", " 'No more than AF1 of the astrometric distances will deviate by more than this from the median (see AM1, AF1).'),\n", " 'AD1_design': ('marcsec',\n", " 'No more than AF1 of the astrometric distances will deviate by more than this from the median (see AM1, AF1).'),\n", " 'AD1_stretch': ('marcsec',\n", " 'No more than AF1 of the astrometric distances will deviate by more than this from the median (see AM1, AF1).'),\n", " 'AD2_minimum': ('marcsec',\n", " 'No more than AF2 of the astrometric distances will deviate by more than this from the median (see AM2, AF2).'),\n", " 'AD2_design': ('marcsec',\n", " 'No more than AF2 of the astrometric distances will deviate by more than this from the median (see AM2, AF2).'),\n", " 'AD2_stretch': ('marcsec',\n", " 'No more than AF2 of the astrometric distances will deviate by more than this from the median (see AM2, AF2).'),\n", " 'AD3_minimum': ('marcsec',\n", " 'No more than AF3 of the astrometric distances will deviate by more than this from the median (see AM3, AF3).'),\n", " 'AD3_design': ('marcsec',\n", " 'No more than AF3 of the astrometric distances will deviate by more than this from the median (see AM3, AF3).'),\n", " 'AD3_stretch': ('marcsec',\n", " 'No more than AF3 of the astrometric distances will deviate by more than this from the median (see AM3, AF3).'),\n", " 'AB1': ('marcsec',\n", " 'The rms difference of separations measured between the r-band and any other filter will be no larger than this.'),\n", " 'AB1_minimum': ('marcsec',\n", " 'The rms difference of separations measured between the r-band and any other filter will be no larger than this.'),\n", " 'AB1_design': ('marcsec',\n", " 'The rms difference of separations measured between the r-band and any other filter will be no larger than this.'),\n", " 'AB1_stretch': ('marcsec',\n", " 'The rms difference of separations measured between the r-band and any other filter will be no larger than this.'),\n", " 'AB2': ('marcsec',\n", " 'No more than ABF1% of the separation measurements between the r-band and any other filter will deviate by more than this amount from the mean.'),\n", " 'AB2_minimum': ('marcsec',\n", " 'No more than ABF1% of the separation measurements between the r-band and any other filter will deviate by more than this amount from the mean.'),\n", " 'AB2_design': ('marcsec',\n", " 'No more than ABF1% of the separation measurements between the r-band and any other filter will deviate by more than this amount from the mean.'),\n", " 'AB2_stretch': ('marcsec',\n", " 'No more than ABF1% of the separation measurements between the r-band and any other filter will deviate by more than this amount from the mean.'),\n", " 'ABF1': ('%',\n", " 'No more than this percentage of separations between r-band and any other filter will deviate by more than AB2 from the mean.'),\n", " 'ABF1_minimum': ('%',\n", " 'No more than this percentage of separations between r-band and any other filter will deviate by more than AB2 from the mean.'),\n", " 'ABF1_design': ('%',\n", " 'No more than this percentage of separations between r-band and any other filter will deviate by more than AB2 from the mean.'),\n", " 'ABF1_stretch': ('%',\n", " 'No more than this percentage of separations between r-band and any other filter will deviate by more than AB2 from the mean.'),\n", " 'AA1': ('marcsec',\n", " 'The median error in absolute astrometric positions in each coordinate will be no larger than this.'),\n", " 'AA1_minimum': ('marcsec',\n", " 'The median error in absolute astrometric positions in each coordinate will be no larger than this.'),\n", " 'AA1_design': ('marcsec',\n", " 'The median error in absolute astrometric positions in each coordinate will be no larger than this.'),\n", " 'AA1_stretch': ('marcsec',\n", " 'The median error in absolute astrometric positions in each coordinate will be no larger than this.'),\n", " 'TE1': ('',\n", " 'The averaged E1, E2, and Ex residual PSF ellipticity correlations must have median absolute value less than this for separations < 1 arcmin.'),\n", " 'TE2': ('',\n", " 'The averaged E1, E2, and Ex residual PSF ellipticity correlations must have median absolute value less than this for separations > 5 arcmin.'),\n", " 'modelPhotRepGal1': ('mmag',\n", " 'The maximum rms of the resolved source model magnitude distribution around the mean value in the 1st faintest bin.'),\n", " 'modelPhotRepGal2': ('mmag',\n", " 'The maximum rms of the resolved source model magnitude distribution around the mean value in the 2nd faintest bin.'),\n", " 'modelPhotRepGal3': ('mmag',\n", " 'The maximum rms of the resolved source model magnitude distribution around the mean value in the 3rd faintest bin.'),\n", " 'modelPhotRepGal4': ('mmag',\n", " 'The maximum rms of the resolved source model magnitude distribution around the mean value in the 4th faintest bin.'),\n", " 'modelPhotRepStar1': ('mmag',\n", " 'The maximum rms of the unresolved source model magnitude distribution around the mean value in the 1st faintest bin.'),\n", " 'modelPhotRepStar2': ('mmag',\n", " 'The maximum rms of the unresolved source model magnitude distribution around the mean value in the 2nd faintest bin.'),\n", " 'modelPhotRepStar3': ('mmag',\n", " 'The maximum rms of the unresolved source model magnitude distribution around the mean value in the 3rd faintest bin.'),\n", " 'modelPhotRepStar4': ('mmag',\n", " 'The maximum rms of the unresolved source model magnitude distribution around the mean value in the 4th faintest bin.'),\n", " 'psfPhotRepStar1': ('mmag',\n", " 'The maximum rms of the unresolved source psf magnitude distribution around the mean value in the 1st faintest bin.'),\n", " 'psfPhotRepStar2': ('mmag',\n", " 'The maximum rms of the unresolved source psf magnitude distribution around the mean value in the 2nd faintest bin.'),\n", " 'psfPhotRepStar3': ('mmag',\n", " 'The maximum rms of the unresolved source psf magnitude distribution around the mean value in the 3rd faintest bin.'),\n", " 'psfPhotRepStar4': ('mmag',\n", " 'The maximum rms of the unresolved source psf magnitude distribution around the mean value in the 4th faintest bin.')}" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "unitsDict" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [], "source": [ "df2 = pd.concat([df, dffixErr])" ] }, { "cell_type": "code", "execution_count": 40, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | AM1 | \n", "AM2 | \n", "AM3 | \n", "config | \n", "
---|---|---|---|---|
0 | \n", "4.599104 | \n", "4.674071 | \n", "NaN | \n", "w50_refErr10 | \n", "
1 | \n", "4.634609 | \n", "4.683638 | \n", "NaN | \n", "w50_ps1 | \n", "
2 | \n", "4.987492 | \n", "5.083376 | \n", "NaN | \n", "w50 | \n", "
3 | \n", "5.209379 | \n", "4.649905 | \n", "NaN | \n", "w50_refErr10 | \n", "
4 | \n", "5.136466 | \n", "4.668592 | \n", "NaN | \n", "w50_ps1 | \n", "
5 | \n", "4.992864 | \n", "4.971158 | \n", "NaN | \n", "w50 | \n", "
6 | \n", "4.704458 | \n", "5.144493 | \n", "NaN | \n", "w50_refErr10 | \n", "
7 | \n", "4.729630 | \n", "5.074105 | \n", "NaN | \n", "w50_ps1 | \n", "
8 | \n", "4.791790 | \n", "4.860627 | \n", "NaN | \n", "w50 | \n", "
9 | \n", "5.633133 | \n", "5.770187 | \n", "NaN | \n", "w50_refErr10 | \n", "
10 | \n", "5.574288 | \n", "5.751263 | \n", "NaN | \n", "w50_ps1 | \n", "
11 | \n", "5.261944 | \n", "5.521375 | \n", "NaN | \n", "w50 | \n", "
12 | \n", "6.176154 | \n", "5.506783 | \n", "NaN | \n", "w50_refErr10 | \n", "
13 | \n", "6.275905 | \n", "6.050604 | \n", "NaN | \n", "w50_ps1 | \n", "
14 | \n", "6.029512 | \n", "6.047009 | \n", "NaN | \n", "w50 | \n", "
15 | \n", "5.889564 | \n", "6.134081 | \n", "5.552543 | \n", "w50_refErr10 | \n", "
16 | \n", "5.887495 | \n", "6.153822 | \n", "5.752080 | \n", "w50_ps1 | \n", "
17 | \n", "5.930031 | \n", "6.172185 | \n", "4.564986 | \n", "w50 | \n", "
18 | \n", "6.069636 | \n", "5.749155 | \n", "6.361017 | \n", "w50_refErr10 | \n", "
19 | \n", "6.248940 | \n", "5.867251 | \n", "6.300364 | \n", "w50_ps1 | \n", "
20 | \n", "7.319564 | \n", "6.585689 | \n", "7.296884 | \n", "w50 | \n", "
21 | \n", "3.569945 | \n", "3.726240 | \n", "3.384210 | \n", "w50_refErr10 | \n", "
22 | \n", "3.584043 | \n", "3.740764 | \n", "3.394108 | \n", "w50_ps1 | \n", "
23 | \n", "3.568332 | \n", "3.709551 | \n", "3.171652 | \n", "w50 | \n", "
24 | \n", "4.624966 | \n", "4.753859 | \n", "4.330667 | \n", "w50_refErr10 | \n", "
25 | \n", "4.671559 | \n", "4.768336 | \n", "4.365674 | \n", "w50_ps1 | \n", "
26 | \n", "4.496229 | \n", "4.557803 | \n", "4.531334 | \n", "w50 | \n", "
27 | \n", "5.094351 | \n", "5.120109 | \n", "5.152024 | \n", "w50_refErr10 | \n", "
28 | \n", "5.097568 | \n", "5.123462 | \n", "5.172630 | \n", "w50_ps1 | \n", "
29 | \n", "5.885524 | \n", "5.977577 | \n", "5.036576 | \n", "w50 | \n", "
30 | \n", "3.544733 | \n", "3.569312 | \n", "3.061433 | \n", "w50_refErr10 | \n", "
31 | \n", "3.509842 | \n", "3.565447 | \n", "3.121952 | \n", "w50_ps1 | \n", "
32 | \n", "3.419372 | \n", "3.335032 | \n", "2.598327 | \n", "w50 | \n", "
33 | \n", "3.408693 | \n", "3.385798 | \n", "3.168869 | \n", "w50_refErr10 | \n", "
34 | \n", "3.376569 | \n", "3.369218 | \n", "3.163556 | \n", "w50_ps1 | \n", "
35 | \n", "3.140067 | \n", "3.072959 | \n", "2.609368 | \n", "w50 | \n", "
36 | \n", "4.077353 | \n", "3.966167 | \n", "3.619436 | \n", "w50_refErr10 | \n", "
37 | \n", "4.069645 | \n", "3.976692 | \n", "3.650448 | \n", "w50_ps1 | \n", "
38 | \n", "4.712038 | \n", "4.463905 | \n", "3.326634 | \n", "w50 | \n", "
39 | \n", "5.028232 | \n", "5.282562 | \n", "4.679492 | \n", "w50_refErr10 | \n", "
40 | \n", "5.023942 | \n", "5.283519 | \n", "4.713056 | \n", "w50_ps1 | \n", "
41 | \n", "4.560869 | \n", "4.680561 | \n", "4.432526 | \n", "w50 | \n", "
42 | \n", "7.873271 | \n", "7.691543 | \n", "7.018579 | \n", "w50_refErr10 | \n", "
43 | \n", "7.933537 | \n", "7.729503 | \n", "7.050278 | \n", "w50_ps1 | \n", "
44 | \n", "8.122135 | \n", "8.041493 | \n", "7.024956 | \n", "w50 | \n", "
0 | \n", "4.836269 | \n", "4.876988 | \n", "NaN | \n", "w50_fixErr | \n", "
1 | \n", "5.237869 | \n", "4.547787 | \n", "NaN | \n", "w50_fixErr | \n", "
2 | \n", "4.723724 | \n", "5.031397 | \n", "NaN | \n", "w50_fixErr | \n", "
3 | \n", "5.645630 | \n", "5.793335 | \n", "NaN | \n", "w50_fixErr | \n", "
4 | \n", "6.104800 | \n", "5.646678 | \n", "NaN | \n", "w50_fixErr | \n", "
5 | \n", "6.232383 | \n", "6.341734 | \n", "5.705591 | \n", "w50_fixErr | \n", "
6 | \n", "6.364804 | \n", "5.911020 | \n", "6.486985 | \n", "w50_fixErr | \n", "
7 | \n", "3.588911 | \n", "3.730127 | \n", "3.387058 | \n", "w50_fixErr | \n", "
8 | \n", "4.700888 | \n", "4.798699 | \n", "4.386471 | \n", "w50_fixErr | \n", "
9 | \n", "5.122123 | \n", "5.146675 | \n", "5.213211 | \n", "w50_fixErr | \n", "
10 | \n", "3.626137 | \n", "3.706853 | \n", "3.272397 | \n", "w50_fixErr | \n", "
11 | \n", "3.443340 | \n", "3.436378 | \n", "3.391356 | \n", "w50_fixErr | \n", "
12 | \n", "4.104107 | \n", "4.008136 | \n", "3.651038 | \n", "w50_fixErr | \n", "
13 | \n", "5.099280 | \n", "5.334973 | \n", "4.768760 | \n", "w50_fixErr | \n", "
14 | \n", "7.975516 | \n", "7.775370 | \n", "7.186471 | \n", "w50_fixErr | \n", "
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |
rerun | \n", "VS. Gaia | \n", "VS. PS1 | \n", "
---|---|---|
w50 (Gaia) | \n",
" \n", " | \n", " |
w50_ps1 (PS1) | \n",
" \n", " | \n", " |
w50_refErr10 (Gaia) | \n",
" \n", " | \n", " |
w50_fixErr (Gaia) | \n",
" \n", " | \n", " |