Use the latest catalogue

nudomarinero · Nov 15, 2018 · affb4ac · affb4ac
1 parent 1d49e01
commit affb4ac
Showing 1 changed file with 32 additions and 38 deletions.
diff --git a/Match_LOFAR_combined_final.ipynb b/Match_LOFAR_combined_final.ipynb
@@ -130,7 +130,7 @@
    },
    "outputs": [],
    "source": [
-    "idp = \"idata/final_pdf_v0.8\""
+    "idp = \"idata/final_pdf_v0.9\""
    ]
   },
   {
@@ -245,8 +245,8 @@
    },
    "outputs": [],
    "source": [
-    "#lofar_all = Table.read(\"data/LOFAR_HBA_T1_DR1_catalog_v0.9.srl.fits\")\n",
-    "lofar_all = Table.read(\"data/LOFAR_HBA_T1_DR1_merge_ID_optical_v0.8.fits\")"
+    "lofar_all = Table.read(\"data/LOFAR_HBA_T1_DR1_catalog_v0.9.srl.fixed.fits\")\n",
+    "#lofar_all = Table.read(\"data/LOFAR_HBA_T1_DR1_merge_ID_optical_v0.8.fits\")"
    ]
   },
   {
@@ -260,8 +260,7 @@
      "data": {
       "text/plain": [
        "array(['AllWISE', 'objID', 'ra', 'dec', 'raErr', 'decErr', 'W1mag',\n",
-       "       'W1magErr', 'i', 'iErr'],\n",
-       "      dtype='<U8')"
+       "       'W1magErr', 'i', 'iErr'], dtype='<U8')"
       ]
      },
      "execution_count": 17,
@@ -283,26 +282,11 @@
     {
      "data": {
       "text/plain": [
-       "array(['Source_Name', 'RA', 'E_RA', 'DEC', 'E_DEC', 'Peak_flux',\n",
-       "       'E_Peak_flux', 'Total_flux', 'E_Total_flux', 'Maj', 'E_Maj', 'Min',\n",
-       "       'E_Min', 'DC_Maj', 'E_DC_Maj', 'DC_Min', 'E_DC_Min', 'PA', 'E_PA',\n",
-       "       'DC_PA', 'E_DC_PA', 'Isl_rms', 'S_Code', 'Mosaic_ID',\n",
-       "       'Masked_Fraction', 'ID_flag', 'ID_name', 'ID_ra', 'ID_dec', 'ML_LR',\n",
-       "       'LGZ_Size', 'LGZ_Width', 'LGZ_PA', 'LGZ_Assoc', 'LGZ_Assoc_Qual',\n",
-       "       'LGZ_ID_Qual', 'AllWISE', 'objID', 'gFApFlux', 'gFApFluxErr',\n",
-       "       'gFApMag', 'gFApMagErr', 'rFApFlux', 'rFApFluxErr', 'rFApMag',\n",
-       "       'rFApMagErr', 'iFApFlux', 'iFApFluxErr', 'iFApMag', 'iFApMagErr',\n",
-       "       'zFApFlux', 'zFApFluxErr', 'zFApMag', 'zFApMagErr', 'yFApFlux',\n",
-       "       'yFApFluxErr', 'yFApMag', 'yFApMagErr', 'w1Flux', 'w1FluxErr',\n",
-       "       'w1Mag', 'w1MagErr', 'w2Flux', 'w2FluxErr', 'w2Mag', 'w2MagErr',\n",
-       "       'w3Flux', 'w3FluxErr', 'w3Mag', 'w3MagErr', 'w4Flux', 'w4FluxErr',\n",
-       "       'w4Mag', 'w4MagErr', 'XrayClass', '2RXS_ID', 'XMMSL2_ID', 'IRClass',\n",
-       "       'z_spec', 'z_source', 'specAGN', 'mqcAGN', 'EBV', 'objName',\n",
-       "       'z_best', 'z_best_source', 'z1_median', 'z1_min', 'z1_max',\n",
-       "       'z1_area', 'z2_median', 'z2_min', 'z2_max', 'z2_area', 'chi_r_eazy',\n",
-       "       'chi_r_atlas', 'chi_r_cosmos', 'chi_r_stellar', 'stellar_type',\n",
-       "       'z_gpz', 'z_gpz_err'],\n",
-       "      dtype='<U15')"
+       "array(['Source_Name', 'RA', 'E_RA', 'E_RA_tot', 'DEC', 'E_DEC',\n",
+       "       'E_DEC_tot', 'Peak_flux', 'E_Peak_flux', 'E_Peak_flux_tot',\n",
+       "       'Total_flux', 'E_Total_flux', 'E_Total_flux_tot', 'Maj', 'E_Maj',\n",
+       "       'Min', 'E_Min', 'PA', 'E_PA', 'Isl_rms', 'S_Code', 'Mosaic_ID',\n",
+       "       'Isl_id'], dtype='<U16')"
       ]
      },
      "execution_count": 18,
@@ -531,10 +515,10 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "/disk2/jsm/prog/anaconda/envs/py36/lib/python3.6/site-packages/astropy/table/column.py:928: RuntimeWarning: invalid value encountered in less\n",
-      "  return getattr(self.data, oper)(other)\n",
-      "/disk2/jsm/prog/anaconda/envs/py36/lib/python3.6/site-packages/astropy/table/column.py:928: RuntimeWarning: invalid value encountered in greater_equal\n",
-      "  return getattr(self.data, oper)(other)\n"
+      "/home/jsm/programs/soft/anaconda3/lib/python3.6/site-packages/astropy/table/column.py:954: RuntimeWarning: invalid value encountered in less\n",
+      "  return getattr(self.data, op)(other)\n",
+      "/home/jsm/programs/soft/anaconda3/lib/python3.6/site-packages/astropy/table/column.py:954: RuntimeWarning: invalid value encountered in greater_equal\n",
+      "  return getattr(self.data, op)(other)\n"
      ]
     }
    ],
@@ -636,9 +620,7 @@
   },
   {
    "cell_type": "markdown",
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "source": [
     "## Maximum Likelihood"
    ]
@@ -685,7 +667,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 39,
+   "execution_count": null,
    "metadata": {
     "collapsed": true
    },
@@ -697,7 +679,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 40,
+   "execution_count": null,
    "metadata": {
     "collapsed": true
    },
@@ -743,11 +725,23 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 41,
+   "execution_count": 1,
    "metadata": {
     "collapsed": true
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "ename": "NameError",
+     "evalue": "name 'search_around_sky' is not defined",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
+      "\u001b[0;32m<ipython-input-1-240a63d46585>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m idx_lofar, idx_i, d2d, d3d = search_around_sky(\n\u001b[0m\u001b[1;32m      2\u001b[0m     coords_lofar, coords_combined[selection], radius*u.arcsec)\n",
+      "\u001b[0;31mNameError\u001b[0m: name 'search_around_sky' is not defined"
+     ]
+    }
+   ],
    "source": [
     "idx_lofar, idx_i, d2d, d3d = search_around_sky(\n",
     "    coords_lofar, coords_combined[selection], radius*u.arcsec)"
@@ -1223,9 +1217,9 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.6.1"
+   "version": "3.6.4"
   }
  },
  "nbformat": 4,
- "nbformat_minor": 1
+ "nbformat_minor": 2
 }