GFA Targets.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 245,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Using matplotlib backend: MacOSX\n",
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "&lt;Table length=59199&gt;\n",
       "<table id=\"table4539968032\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>RELEASE</th><th>BRICKID</th><th>BRICKNAME</th><th>OBJID</th><th>TYPE</th><th>RA</th><th>DEC</th><th>RA_IVAR</th><th>DEC_IVAR</th><th>DCHISQ [5]</th><th>EBV</th><th>FLUX_G</th><th>FLUX_R</th><th>FLUX_Z</th><th>FLUX_W1</th><th>FLUX_W2</th><th>FLUX_W3</th><th>FLUX_W4</th><th>FLUX_IVAR_G</th><th>FLUX_IVAR_R</th><th>FLUX_IVAR_Z</th><th>FLUX_IVAR_W1</th><th>FLUX_IVAR_W2</th><th>FLUX_IVAR_W3</th><th>FLUX_IVAR_W4</th><th>MW_TRANSMISSION_G</th><th>MW_TRANSMISSION_R</th><th>MW_TRANSMISSION_Z</th><th>MW_TRANSMISSION_W1</th><th>MW_TRANSMISSION_W2</th><th>MW_TRANSMISSION_W3</th><th>MW_TRANSMISSION_W4</th><th>NOBS_G</th><th>NOBS_R</th><th>NOBS_Z</th><th>NOBS_W1</th><th>NOBS_W2</th><th>NOBS_W3</th><th>NOBS_W4</th><th>RCHISQ_G</th><th>RCHISQ_R</th><th>RCHISQ_Z</th><th>RCHISQ_W1</th><th>RCHISQ_W2</th><th>RCHISQ_W3</th><th>RCHISQ_W4</th><th>FRACFLUX_G</th><th>FRACFLUX_R</th><th>FRACFLUX_Z</th><th>FRACFLUX_W1</th><th>FRACFLUX_W2</th><th>FRACFLUX_W3</th><th>FRACFLUX_W4</th><th>FRACMASKED_G</th><th>FRACMASKED_R</th><th>FRACMASKED_Z</th><th>FRACIN_G</th><th>FRACIN_R</th><th>FRACIN_Z</th><th>ANYMASK_G</th><th>ANYMASK_R</th><th>ANYMASK_Z</th><th>ALLMASK_G</th><th>ALLMASK_R</th><th>ALLMASK_Z</th><th>WISEMASK_W1</th><th>WISEMASK_W2</th><th>PSFSIZE_G</th><th>PSFSIZE_R</th><th>PSFSIZE_Z</th><th>PSFDEPTH_G</th><th>PSFDEPTH_R</th><th>PSFDEPTH_Z</th><th>GALDEPTH_G</th><th>GALDEPTH_R</th><th>GALDEPTH_Z</th><th>WISE_COADD_ID</th><th>FRACDEV</th><th>FRACDEV_IVAR</th><th>SHAPEDEV_R</th><th>SHAPEDEV_R_IVAR</th><th>SHAPEDEV_E1</th><th>SHAPEDEV_E1_IVAR</th><th>SHAPEDEV_E2</th><th>SHAPEDEV_E2_IVAR</th><th>SHAPEEXP_R</th><th>SHAPEEXP_R_IVAR</th><th>SHAPEEXP_E1</th><th>SHAPEEXP_E1_IVAR</th><th>SHAPEEXP_E2</th><th>SHAPEEXP_E2_IVAR</th></tr></thead>\n",
       "<thead><tr><th>int32</th><th>int32</th><th>str8</th><th>int32</th><th>str4</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>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>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>uint8</th><th>uint8</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>str8</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></tr></thead>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>0</td><td>PSF</td><td>66.1454199468</td><td>66.9084831807</td><td>4.72097e+09</td><td>4.43343e+10</td><td>1696.91 .. 0.0</td><td>0.5057</td><td>0.0</td><td>89.9539</td><td>0.0</td><td>-3.52</td><td>-6.55881</td><td>-68.2792</td><td>-208.378</td><td>0.0</td><td>0.196921</td><td>0.0</td><td>2.415</td><td>0.540222</td><td>0.00139922</td><td>3.02673e-05</td><td>0.223806</td><td>0.364807</td><td>0.568904</td><td>0.917868</td><td>0.948729</td><td>0.988838</td><td>0.995771</td><td>0</td><td>1</td><td>0</td><td>101</td><td>100</td><td>16</td><td>16</td><td>0.0</td><td>4.08859</td><td>0.0</td><td>16.8816</td><td>2.33363</td><td>0.549687</td><td>0.146833</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0798347</td><td>0.071718</td><td>0.684666</td><td>1.41097</td><td>0.0</td><td>0.889268</td><td>0.0</td><td>0.0</td><td>0.983889</td><td>0.0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>1</td><td>DEV</td><td>66.1253740885</td><td>66.9083976993</td><td>2.33287e+08</td><td>3.8971e+11</td><td>4220.77 .. 3351.96</td><td>0.509157</td><td>0.0</td><td>2341.74</td><td>0.0</td><td>-6.83844</td><td>-12.82</td><td>-1053.54</td><td>-3619.38</td><td>0.0</td><td>0.00426048</td><td>0.0</td><td>0.0385718</td><td>0.00983925</td><td>3.44313e-05</td><td>1.14029e-06</td><td>0.221527</td><td>0.362301</td><td>0.566714</td><td>0.917331</td><td>0.948388</td><td>0.988762</td><td>0.995742</td><td>0</td><td>1</td><td>0</td><td>100</td><td>100</td><td>15</td><td>16</td><td>0.0</td><td>29.6654</td><td>0.0</td><td>1.54533</td><td>0.260997</td><td>0.228251</td><td>0.0608615</td><td>0.0</td><td>0.159622</td><td>0.0</td><td>4.45708</td><td>1.9864</td><td>0.124554</td><td>0.329664</td><td>0.0</td><td>0.615911</td><td>0.0</td><td>0.0</td><td>0.196702</td><td>0.0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>1.0</td><td>1.0</td><td>146.407</td><td>146830.0</td><td>-0.998976</td><td>2.3774e+06</td><td>0.00477209</td><td>1.15825e+06</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>2</td><td>SIMP</td><td>66.119496885</td><td>66.9087245425</td><td>5.4273e+08</td><td>1.01007e+10</td><td>377.292 .. 0.0</td><td>0.510193</td><td>0.0</td><td>158.429</td><td>0.0</td><td>-6.94743</td><td>-3.93703</td><td>2.16558</td><td>229.356</td><td>0.0</td><td>0.0140095</td><td>0.0</td><td>2.38048</td><td>0.52552</td><td>0.00135402</td><td>2.88385e-05</td><td>0.220849</td><td>0.361554</td><td>0.56606</td><td>0.91717</td><td>0.948286</td><td>0.988739</td><td>0.995733</td><td>0</td><td>1</td><td>0</td><td>101</td><td>101</td><td>15</td><td>16</td><td>0.0</td><td>2.67833</td><td>0.0</td><td>3.21104</td><td>0.893549</td><td>0.805723</td><td>0.235319</td><td>0.0</td><td>0.136904</td><td>0.0</td><td>0.0558047</td><td>0.208612</td><td>42.9648</td><td>2.19395</td><td>0.0</td><td>0.947699</td><td>0.0</td><td>0.0</td><td>0.826869</td><td>0.0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.45</td><td>0.45</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>3</td><td>EXP</td><td>65.8779512588</td><td>66.9094045879</td><td>4.38489e+07</td><td>1.46808e+11</td><td>12519.8 .. 11703.1</td><td>0.557375</td><td>0.0</td><td>4122.78</td><td>0.0</td><td>24.2406</td><td>-2.26082</td><td>-377.466</td><td>-366.082</td><td>0.0</td><td>0.00378761</td><td>0.0</td><td>0.144967</td><td>0.0384864</td><td>0.000147072</td><td>4.89567e-06</td><td>0.192061</td><td>0.329089</td><td>0.537041</td><td>0.909866</td><td>0.943641</td><td>0.987704</td><td>0.995339</td><td>0</td><td>0</td><td>0</td><td>106</td><td>107</td><td>17</td><td>17</td><td>0.0</td><td>6.7321</td><td>0.0</td><td>16.982</td><td>1.75464</td><td>0.477412</td><td>0.101501</td><td>0.0</td><td>0.177873</td><td>0.0</td><td>1.39863</td><td>5.16148</td><td>0.250226</td><td>2.0485</td><td>0.0</td><td>0.743376</td><td>0.0</td><td>0.0</td><td>0.316783</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>38.1567</td><td>77625.0</td><td>-0.878348</td><td>2.12995e+06</td><td>0.00094475</td><td>44255.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>4</td><td>DEV</td><td>65.8843825997</td><td>66.9089025206</td><td>1.82127e+08</td><td>5.45169e+10</td><td>1194.18 .. 729.012</td><td>0.555777</td><td>0.0</td><td>433.258</td><td>0.0</td><td>-29.0117</td><td>-9.50801</td><td>-73.9846</td><td>-568.556</td><td>0.0</td><td>0.0175372</td><td>0.0</td><td>0.377651</td><td>0.0945586</td><td>0.00031626</td><td>9.04275e-06</td><td>0.192972</td><td>0.330139</td><td>0.537999</td><td>0.910112</td><td>0.943798</td><td>0.987739</td><td>0.995353</td><td>0</td><td>1</td><td>0</td><td>105</td><td>108</td><td>16</td><td>16</td><td>0.0</td><td>9.12232</td><td>0.0</td><td>13.7438</td><td>1.72827</td><td>0.629149</td><td>0.113117</td><td>0.0</td><td>1.23852</td><td>0.0</td><td>0.483345</td><td>0.222309</td><td>2.83223</td><td>0.524432</td><td>0.0</td><td>0.842701</td><td>0.0</td><td>0.0</td><td>0.579986</td><td>0.0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>1.0</td><td>1.0</td><td>20.9819</td><td>6423.28</td><td>-0.987801</td><td>26188.2</td><td>0.0685127</td><td>43082.1</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>5</td><td>DEV</td><td>65.886920988</td><td>66.9143779503</td><td>3.87727e+09</td><td>2.05513e+10</td><td>1513.68 .. 1534.54</td><td>0.557447</td><td>0.0</td><td>18.7607</td><td>0.0</td><td>25.5993</td><td>18.0026</td><td>-33.9599</td><td>-108.699</td><td>0.0</td><td>4.42238</td><td>0.0</td><td>2.36973</td><td>0.556949</td><td>0.0014191</td><td>2.92335e-05</td><td>0.19202</td><td>0.329041</td><td>0.536998</td><td>0.909854</td><td>0.943634</td><td>0.987703</td><td>0.995339</td><td>0</td><td>1</td><td>0</td><td>105</td><td>105</td><td>16</td><td>16</td><td>0.0</td><td>1.27761</td><td>0.0</td><td>61.5173</td><td>5.23349</td><td>0.649046</td><td>0.143124</td><td>0.0</td><td>2.37979e-05</td><td>0.0</td><td>0.0713721</td><td>0.132343</td><td>0.343749</td><td>1.32166</td><td>0.0</td><td>0.00927128</td><td>0.0</td><td>0.0</td><td>1.00524</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.43907</td><td>0.0</td><td>0.0</td><td>10.153</td><td>0.0</td><td>0.0</td><td>3.73677</td><td>0.0</td><td>0668p666</td><td>1.0</td><td>1.0</td><td>0.470398</td><td>728.064</td><td>-0.0778361</td><td>108.724</td><td>0.224055</td><td>181.639</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>6</td><td>DEV</td><td>65.7407595738</td><td>66.9113056262</td><td>6.51909e+08</td><td>5.73935e+10</td><td>4804.31 .. 7457.18</td><td>0.588289</td><td>0.0</td><td>209.79</td><td>0.0</td><td>75.1467</td><td>50.0695</td><td>-44.7712</td><td>55.9092</td><td>0.0</td><td>0.075517</td><td>0.0</td><td>0.912834</td><td>0.229799</td><td>0.000707816</td><td>1.85822e-05</td><td>0.175265</td><td>0.309415</td><td>0.518839</td><td>0.905111</td><td>0.940609</td><td>0.987027</td><td>0.995081</td><td>0</td><td>1</td><td>0</td><td>107</td><td>108</td><td>17</td><td>17</td><td>0.0</td><td>44.836</td><td>0.0</td><td>1035.23</td><td>62.7063</td><td>0.798494</td><td>0.139743</td><td>0.0</td><td>-3.5227e-11</td><td>0.0</td><td>2.46035</td><td>2.40414</td><td>1.36393</td><td>6.95074</td><td>0.0</td><td>0.757401</td><td>0.0</td><td>0.0</td><td>0.688342</td><td>0.0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>1.0</td><td>1.0</td><td>7.0184</td><td>101.636</td><td>-0.939724</td><td>6908.47</td><td>-0.122758</td><td>6419.66</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>7</td><td>DEV</td><td>65.718570237</td><td>66.9113793292</td><td>1.31202e+08</td><td>1.54705e+11</td><td>4598.37 .. 3980.31</td><td>0.593259</td><td>0.0</td><td>1284.52</td><td>0.0</td><td>-563.065</td><td>-337.725</td><td>-351.007</td><td>287.985</td><td>0.0</td><td>0.00407557</td><td>0.0</td><td>0.17652</td><td>0.0451231</td><td>0.000151452</td><td>4.61117e-06</td><td>0.172705</td><td>0.306364</td><td>0.515971</td><td>0.904349</td><td>0.940123</td><td>0.986918</td><td>0.99504</td><td>0</td><td>1</td><td>0</td><td>108</td><td>109</td><td>16</td><td>16</td><td>0.0</td><td>14.9423</td><td>0.0</td><td>3007.66</td><td>180.161</td><td>0.405844</td><td>0.100827</td><td>0.0</td><td>-4.0628e-09</td><td>0.0</td><td>3.48238</td><td>3.83555</td><td>1.25634</td><td>2.4523</td><td>0.0</td><td>0.845442</td><td>0.0</td><td>0.0</td><td>0.279163</td><td>0.0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>1.0</td><td>1.0</td><td>45.6157</td><td>31747.7</td><td>-0.983553</td><td>113342.0</td><td>-0.0318223</td><td>119569.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>8</td><td>SIMP</td><td>65.678707768</td><td>66.9113490784</td><td>8.56763e+08</td><td>1.02459e+10</td><td>1456.06 .. 1496.32</td><td>0.60102</td><td>0.0</td><td>36.2965</td><td>0.0</td><td>8.12262</td><td>4.55231</td><td>30.8309</td><td>43.976</td><td>0.0</td><td>0.355847</td><td>0.0</td><td>2.6601</td><td>0.602519</td><td>0.00145568</td><td>3.10341e-05</td><td>0.168783</td><td>0.301659</td><td>0.511524</td><td>0.90316</td><td>0.939364</td><td>0.986748</td><td>0.994975</td><td>0</td><td>1</td><td>0</td><td>112</td><td>113</td><td>17</td><td>17</td><td>0.0</td><td>40.5337</td><td>0.0</td><td>1323.86</td><td>14.8992</td><td>0.991515</td><td>0.206905</td><td>0.0</td><td>-1.02553e-09</td><td>0.0</td><td>0.0749197</td><td>0.0864372</td><td>0.101433</td><td>1.4362</td><td>0.0</td><td>0.784969</td><td>0.0</td><td>0.0</td><td>0.914711</td><td>0.0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.45</td><td>0.45</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635401</td><td>0658p670</td><td>9</td><td>PSF</td><td>65.9897961322</td><td>66.910932833</td><td>7.72493e+07</td><td>4.59567e+09</td><td>99.4883 .. 0.0</td><td>0.534237</td><td>0.0</td><td>8.59588</td><td>0.0</td><td>74.0008</td><td>42.4494</td><td>-69.8061</td><td>-16.6888</td><td>0.0</td><td>1.37427</td><td>0.0</td><td>2.47548</td><td>0.552212</td><td>0.0013634</td><td>2.83265e-05</td><td>0.205677</td><td>0.344628</td><td>0.551081</td><td>0.91344</td><td>0.945916</td><td>0.988212</td><td>0.995532</td><td>0</td><td>1</td><td>0</td><td>103</td><td>104</td><td>15</td><td>15</td><td>0.0</td><td>0.0303825</td><td>0.0</td><td>88.5364</td><td>7.71142</td><td>0.845486</td><td>0.304569</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0299704</td><td>0.0309096</td><td>0.0554335</td><td>4.04783</td><td>0.0</td><td>0.845932</td><td>0.0</td><td>0.0</td><td>0.993374</td><td>0.0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1284</td><td>PSF</td><td>65.3233694398</td><td>67.1239472917</td><td>3.06624e+10</td><td>1.25429e+11</td><td>6123.67 .. 0.0</td><td>0.70866</td><td>0.0</td><td>33.74</td><td>0.0</td><td>66.0719</td><td>33.2658</td><td>12.8436</td><td>317.969</td><td>0.0</td><td>5.37924</td><td>0.0</td><td>2.53352</td><td>0.565949</td><td>0.00155035</td><td>3.06422e-05</td><td>0.122729</td><td>0.243389</td><td>0.453654</td><td>0.886835</td><td>0.928899</td><td>0.984393</td><td>0.994078</td><td>0</td><td>1</td><td>0</td><td>104</td><td>105</td><td>17</td><td>17</td><td>0.0</td><td>0.789632</td><td>0.0</td><td>7.7163</td><td>1.15195</td><td>0.742796</td><td>0.243963</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0145637</td><td>0.0139571</td><td>0.103246</td><td>0.217924</td><td>0.0</td><td>0.00128149</td><td>0.0</td><td>0.0</td><td>1.00305</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1285</td><td>PSF</td><td>64.9916972979</td><td>67.124046637</td><td>1.75825e+11</td><td>6.79699e+11</td><td>31227.1 .. 31253.7</td><td>0.672425</td><td>0.0</td><td>107.834</td><td>0.0</td><td>70.9294</td><td>37.3654</td><td>-8.53234</td><td>294.975</td><td>0.0</td><td>2.68549</td><td>0.0</td><td>2.862</td><td>0.636959</td><td>0.0015424</td><td>3.16113e-05</td><td>0.136625</td><td>0.261625</td><td>0.472364</td><td>0.892297</td><td>0.932409</td><td>0.985185</td><td>0.99438</td><td>0</td><td>1</td><td>0</td><td>115</td><td>116</td><td>17</td><td>17</td><td>0.0</td><td>1.98048</td><td>0.0</td><td>8.861</td><td>1.79171</td><td>0.457314</td><td>0.141661</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.00665544</td><td>0.00979368</td><td>0.105022</td><td>0.376296</td><td>0.0</td><td>0.000217599</td><td>0.0</td><td>0.0</td><td>1.00407</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1286</td><td>PSF</td><td>65.3687565906</td><td>67.1238132712</td><td>3.42514e+08</td><td>1.38894e+09</td><td>82.0874 .. 0.0</td><td>0.702933</td><td>0.0</td><td>2.89285</td><td>0.0</td><td>-1.99261</td><td>6.13817</td><td>3.7048</td><td>98.6695</td><td>0.0</td><td>9.80942</td><td>0.0</td><td>2.56054</td><td>0.560574</td><td>0.00145441</td><td>2.95592e-05</td><td>0.124827</td><td>0.246184</td><td>0.456561</td><td>0.887696</td><td>0.929453</td><td>0.984518</td><td>0.994126</td><td>0</td><td>1</td><td>0</td><td>104</td><td>104</td><td>16</td><td>16</td><td>0.0</td><td>0.559269</td><td>0.0</td><td>796.636</td><td>35.9142</td><td>0.417957</td><td>0.125165</td><td>0.0</td><td>0.0</td><td>0.0</td><td>19.1612</td><td>2.57739</td><td>3.67555</td><td>0.842405</td><td>0.0</td><td>0.0154083</td><td>0.0</td><td>0.0</td><td>0.999168</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1287</td><td>PSF</td><td>65.2158544359</td><td>67.12445038</td><td>6.45784e+09</td><td>2.60891e+10</td><td>1425.73 .. 0.0</td><td>0.709012</td><td>0.0</td><td>13.8012</td><td>0.0</td><td>36.3678</td><td>20.0063</td><td>46.7205</td><td>285.244</td><td>0.0</td><td>7.4852</td><td>0.0</td><td>2.63095</td><td>0.597407</td><td>0.00146989</td><td>3.00888e-05</td><td>0.122601</td><td>0.243218</td><td>0.453476</td><td>0.886782</td><td>0.928865</td><td>0.984385</td><td>0.994075</td><td>0</td><td>1</td><td>0</td><td>109</td><td>110</td><td>16</td><td>16</td><td>0.0</td><td>0.67906</td><td>0.0</td><td>7.06216</td><td>0.780885</td><td>0.57506</td><td>0.249672</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0253799</td><td>0.0393728</td><td>0.0444715</td><td>0.159491</td><td>0.0</td><td>0.00225026</td><td>0.0</td><td>0.0</td><td>1.00242</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1288</td><td>PSF</td><td>65.1876739446</td><td>67.1240980751</td><td>5.30548e+08</td><td>2.03426e+09</td><td>126.79 .. 0.0</td><td>0.705701</td><td>0.0</td><td>3.62818</td><td>0.0</td><td>9.45561</td><td>3.54384</td><td>-91.3356</td><td>-179.871</td><td>0.0</td><td>9.63168</td><td>0.0</td><td>2.72391</td><td>0.619019</td><td>0.00147398</td><td>2.92692e-05</td><td>0.123809</td><td>0.244829</td><td>0.455154</td><td>0.88728</td><td>0.929185</td><td>0.984458</td><td>0.994103</td><td>0</td><td>1</td><td>0</td><td>111</td><td>111</td><td>16</td><td>16</td><td>0.0</td><td>0.55619</td><td>0.0</td><td>3.97758</td><td>0.787599</td><td>0.719726</td><td>0.190109</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.00275715</td><td>0.00187172</td><td>0.00517061</td><td>0.0575532</td><td>0.0</td><td>0.0137379</td><td>0.0</td><td>0.0</td><td>0.998547</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0668p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1289</td><td>DEV</td><td>65.4727570339</td><td>67.1241042311</td><td>1.99856e+08</td><td>1.05377e+09</td><td>112.276 .. 95.9568</td><td>0.697774</td><td>0.0</td><td>10.6853</td><td>0.0</td><td>127.884</td><td>93.099</td><td>51.1013</td><td>244.551</td><td>0.0</td><td>1.35321</td><td>0.0</td><td>1.47616</td><td>0.335514</td><td>0.000963151</td><td>2.38254e-05</td><td>0.126748</td><td>0.24873</td><td>0.459196</td><td>0.888472</td><td>0.929952</td><td>0.984631</td><td>0.994169</td><td>0</td><td>1</td><td>0</td><td>105</td><td>105</td><td>15</td><td>16</td><td>0.0</td><td>0.790545</td><td>0.0</td><td>85.5596</td><td>3.56815</td><td>0.882196</td><td>0.192182</td><td>0.0</td><td>4.79349e-10</td><td>0.0</td><td>0.0594682</td><td>0.0629412</td><td>0.216249</td><td>1.19884</td><td>0.0</td><td>0.11952</td><td>0.0</td><td>0.0</td><td>0.759118</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0668p666</td><td>1.0</td><td>1.0</td><td>1.98661</td><td>16.8262</td><td>-0.00192571</td><td>28.1853</td><td>0.0741796</td><td>84.5092</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1290</td><td>PSF</td><td>64.8907347785</td><td>67.1243826577</td><td>9.1484e+08</td><td>3.37207e+09</td><td>229.173 .. 0.0</td><td>0.647691</td><td>0.0</td><td>5.13661</td><td>0.0</td><td>13.2562</td><td>7.60401</td><td>15.2049</td><td>75.0356</td><td>0.0</td><td>8.68581</td><td>0.0</td><td>2.83573</td><td>0.653447</td><td>0.00156348</td><td>3.11861e-05</td><td>0.147004</td><td>0.274853</td><td>0.485578</td><td>0.896045</td><td>0.934812</td><td>0.985726</td><td>0.994586</td><td>0</td><td>1</td><td>0</td><td>116</td><td>118</td><td>17</td><td>17</td><td>0.0</td><td>0.75802</td><td>0.0</td><td>14.5024</td><td>1.89376</td><td>0.653331</td><td>0.131007</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.288851</td><td>0.282073</td><td>0.54828</td><td>1.73846</td><td>0.0</td><td>0.00942479</td><td>0.0</td><td>0.0</td><td>1.00332</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0630p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1291</td><td>PSF</td><td>64.8969161019</td><td>67.1245855192</td><td>2.42908e+09</td><td>9.10358e+09</td><td>592.75 .. 0.0</td><td>0.649154</td><td>0.0</td><td>8.55793</td><td>0.0</td><td>9.78951</td><td>0.847944</td><td>-5.97872</td><td>-15.7241</td><td>0.0</td><td>7.97132</td><td>0.0</td><td>2.83658</td><td>0.651973</td><td>0.00156741</td><td>3.09989e-05</td><td>0.146368</td><td>0.274052</td><td>0.484786</td><td>0.895823</td><td>0.93467</td><td>0.985694</td><td>0.994574</td><td>0</td><td>1</td><td>0</td><td>116</td><td>117</td><td>17</td><td>17</td><td>0.0</td><td>0.55794</td><td>0.0</td><td>13.5356</td><td>1.98146</td><td>0.689842</td><td>0.13457</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.403266</td><td>3.52936</td><td>1.49172</td><td>8.35876</td><td>0.0</td><td>0.00556062</td><td>0.0</td><td>0.0</td><td>1.00462</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0630p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1292</td><td>PSF</td><td>64.9455640848</td><td>67.1247950734</td><td>1.71846e+09</td><td>6.28892e+09</td><td>412.084 .. 0.0</td><td>0.660752</td><td>0.0</td><td>6.95395</td><td>0.0</td><td>14.5172</td><td>12.9472</td><td>31.3015</td><td>454.229</td><td>0.0</td><td>8.5403</td><td>0.0</td><td>2.79989</td><td>0.643768</td><td>0.00150858</td><td>3.18034e-05</td><td>0.141428</td><td>0.267786</td><td>0.478555</td><td>0.894064</td><td>0.933542</td><td>0.98544</td><td>0.994477</td><td>0</td><td>1</td><td>0</td><td>113</td><td>115</td><td>17</td><td>17</td><td>0.0</td><td>0.38592</td><td>0.0</td><td>23.2859</td><td>3.95242</td><td>0.726715</td><td>0.1616</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0112581</td><td>0.0075124</td><td>0.060389</td><td>0.126114</td><td>0.0</td><td>0.00584733</td><td>0.0</td><td>0.0</td><td>1.00378</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0630p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>4000</td><td>635400</td><td>0651p670</td><td>1293</td><td>PSF</td><td>64.8773228615</td><td>67.1245995638</td><td>1.7692e+08</td><td>6.32096e+08</td><td>45.0033 .. 0.0</td><td>0.644825</td><td>0.0</td><td>2.20656</td><td>0.0</td><td>28.3424</td><td>22.4107</td><td>21.0344</td><td>380.882</td><td>0.0</td><td>9.24301</td><td>0.0</td><td>2.88676</td><td>0.66026</td><td>0.00156719</td><td>3.16408e-05</td><td>0.148256</td><td>0.276428</td><td>0.487132</td><td>0.896481</td><td>0.935091</td><td>0.985789</td><td>0.99461</td><td>0</td><td>1</td><td>0</td><td>118</td><td>120</td><td>17</td><td>17</td><td>0.0</td><td>0.935654</td><td>0.0</td><td>6.97217</td><td>1.60143</td><td>0.694023</td><td>0.181933</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0319206</td><td>0.020021</td><td>0.349846</td><td>0.155053</td><td>0.0</td><td>0.0284189</td><td>0.0</td><td>0.0</td><td>0.995782</td><td>0.0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.0</td><td>1.46993</td><td>0.0</td><td>0.0</td><td>9.79258</td><td>0.0</td><td>0.0</td><td>3.34974</td><td>0.0</td><td>0630p666</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=59199>\n",
       "RELEASE BRICKID BRICKNAME ... SHAPEEXP_E1_IVAR SHAPEEXP_E2 SHAPEEXP_E2_IVAR\n",
       " int32   int32     str8   ...     float32        float32       float32     \n",
       "------- ------- --------- ... ---------------- ----------- ----------------\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...      2.12995e+06  0.00094475          44255.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "   4000  635401  0658p670 ...              0.0         0.0              0.0\n",
       "    ...     ...       ... ...              ...         ...              ...\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0\n",
       "   4000  635400  0651p670 ...              0.0         0.0              0.0"
      ]
     },
     "execution_count": 245,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Working on identifying which celestial objects land on the GFAs\n",
    "\n",
    "from astropy.io import fits\n",
    "from astropy.table import Table\n",
    "import desimodel.io\n",
    "import desimodel.footprint\n",
    "%pylab\n",
    "\n",
    "sweep = Table.read('sweep-050p030-060p035.fits')\n",
    "sweep1 = Table.read('sweep-060p065-070p070.fits')\n",
    "\n",
    "plot(sweep['RA'], sweep['DEC'], 'b.')\n",
    "tiles = desimodel.io.load_tiles()\n",
    "#print(tiles.dtype.names)\n",
    "#print(tiles['TILEID'])\n",
    "\n",
    "radec = desimodel.footprint.get_tile_radec(4)\n",
    "#print(radec)\n",
    "#print(tiles['RA'])\n",
    "#print(tiles['DEC'])\n",
    "for targets in sweep:\n",
    "    if(targets['FLUX_R'] > 0):\n",
    "        print(targets['FLUX_R'])\n",
    "#sweep\n",
    "sweep1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 237,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "\"\"\"newra, newdec = xy2radec(8.37, -10.65, -138.345, -333.179)\n",
    "8.927313423598427\n",
    "-9.324956250231294\"\"\"\n",
    "\n",
    "import desimodel.io\n",
    "gfatable = desimodel.io.load_gfa()\n",
    "gfatable\n",
    "import desimodel.focalplane\n",
    "tiles = desimodel.io.load_tiles()\n",
    "#print(indices)\n",
    "#indices = desimodel.footprint.find_tiles_over_point(tiles, 8.37, -10.65)\n",
    "#[8254, 15072, 198, 2213, 4225, 6242, 8255, 10265, 12279]\n",
    "#print(tiles[8254])\n",
    "#tile index 8254 has tileid 23658\n",
    "\n",
    "telra, teldec = desimodel.footprint.get_tile_radec(23658)\n",
    "#print(telra, teldec)\n",
    "ra, dec = desimodel.focalplane.xy2radec(7.11, -10.53, 116.279135121, -372.885546514)\n",
    "#print(ra, dec)\n",
    "#Certain gfa location is in ra, dec respectively 6.644525362152656 -9.055425745149217"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 238,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7.11 -10.53\n",
      "The target is on GFA 0\n",
      "The target is on GFA 8\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "&lt;Table length=40&gt;\n",
       "<table id=\"table4536110272\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>PETAL</th><th>CORNER</th><th>X</th><th>Y</th><th>Z</th><th>Q</th><th>RADIUS_DEG</th><th>RADIUS_MM</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th>mm</th><th>mm</th><th>mm</th><th>degrees</th><th></th><th></th></tr></thead>\n",
       "<thead><tr><th>int64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>0</td><td>0</td><td>116.279135121</td><td>-372.885546514</td><td>-17.053</td><td>287.31935812</td><td>1.54426530071</td><td>390.595017971</td></tr>\n",
       "<tr><td>0</td><td>1</td><td>121.005023896</td><td>-387.4297715</td><td>-18.487</td><td>287.345026987</td><td>1.5999816328</td><td>405.886737469</td></tr>\n",
       "<tr><td>0</td><td>2</td><td>91.784287159</td><td>-396.910065893</td><td>-18.631</td><td>283.02060759</td><td>1.60541232205</td><td>407.384285137</td></tr>\n",
       "<tr><td>0</td><td>3</td><td>87.0587074141</td><td>-382.364889847</td><td>-17.198</td><td>282.826719301</td><td>1.54995479057</td><td>392.150643917</td></tr>\n",
       "<tr><td>1</td><td>0</td><td>313.248421438</td><td>-233.323583313</td><td>-17.053</td><td>323.31935812</td><td>1.54426530071</td><td>390.595017971</td></tr>\n",
       "<tr><td>1</td><td>1</td><td>325.620626724</td><td>-242.312300771</td><td>-18.487</td><td>323.345026987</td><td>1.5999816328</td><td>405.886737469</td></tr>\n",
       "<tr><td>1</td><td>2</td><td>307.552931346</td><td>-267.157538162</td><td>-18.631</td><td>319.02060759</td><td>1.60541232205</td><td>407.384285137</td></tr>\n",
       "<tr><td>1</td><td>3</td><td>295.180417053</td><td>-258.167869637</td><td>-17.198</td><td>318.826719301</td><td>1.54995479057</td><td>392.150643917</td></tr>\n",
       "<tr><td>2</td><td>0</td><td>390.567457688</td><td>-4.6399416638</td><td>-17.053</td><td>359.31935812</td><td>1.54426530071</td><td>390.595017971</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>7</td><td>2</td><td>-405.846809099</td><td>35.3599112033</td><td>-18.631</td><td>175.02060759</td><td>1.60541232205</td><td>407.384285137</td></tr>\n",
       "<tr><td>7</td><td>3</td><td>-390.553240191</td><td>35.3594980285</td><td>-17.198</td><td>174.826719301</td><td>1.54995479057</td><td>392.150643917</td></tr>\n",
       "<tr><td>8</td><td>0</td><td>-318.703000001</td><td>-225.815999995</td><td>-17.053</td><td>215.31935812</td><td>1.54426530071</td><td>390.595017971</td></tr>\n",
       "<tr><td>8</td><td>1</td><td>-331.075000003</td><td>-234.805000001</td><td>-18.487</td><td>215.345026987</td><td>1.5999816328</td><td>405.886737469</td></tr>\n",
       "<tr><td>8</td><td>2</td><td>-349.121000002</td><td>-209.943999995</td><td>-18.631</td><td>211.02060759</td><td>1.60541232205</td><td>407.384285137</td></tr>\n",
       "<tr><td>8</td><td>3</td><td>-336.747999992</td><td>-200.955000001</td><td>-17.198</td><td>210.826719301</td><td>1.54995479057</td><td>392.150643917</td></tr>\n",
       "<tr><td>9</td><td>0</td><td>-125.10482863</td><td>-370.01790486</td><td>-17.053</td><td>251.31935812</td><td>1.54426530071</td><td>390.595017971</td></tr>\n",
       "<tr><td>9</td><td>1</td><td>-129.83038525</td><td>-384.56223777</td><td>-18.487</td><td>251.345026987</td><td>1.5999816328</td><td>405.886737469</td></tr>\n",
       "<tr><td>9</td><td>2</td><td>-159.04283509</td><td>-375.05643893</td><td>-18.631</td><td>247.02060759</td><td>1.60541232205</td><td>407.384285137</td></tr>\n",
       "<tr><td>9</td><td>3</td><td>-154.31646944</td><td>-360.51151824</td><td>-17.198</td><td>246.826719301</td><td>1.54995479057</td><td>392.150643917</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=40>\n",
       "PETAL CORNER       X        ...       Q         RADIUS_DEG    RADIUS_MM  \n",
       "                   mm       ...    degrees                               \n",
       "int64 int64     float64     ...    float64       float64       float64   \n",
       "----- ------ -------------- ... ------------- ------------- -------------\n",
       "    0      0  116.279135121 ...  287.31935812 1.54426530071 390.595017971\n",
       "    0      1  121.005023896 ... 287.345026987  1.5999816328 405.886737469\n",
       "    0      2   91.784287159 ...  283.02060759 1.60541232205 407.384285137\n",
       "    0      3  87.0587074141 ... 282.826719301 1.54995479057 392.150643917\n",
       "    1      0  313.248421438 ...  323.31935812 1.54426530071 390.595017971\n",
       "    1      1  325.620626724 ... 323.345026987  1.5999816328 405.886737469\n",
       "    1      2  307.552931346 ...  319.02060759 1.60541232205 407.384285137\n",
       "    1      3  295.180417053 ... 318.826719301 1.54995479057 392.150643917\n",
       "    2      0  390.567457688 ...  359.31935812 1.54426530071 390.595017971\n",
       "  ...    ...            ... ...           ...           ...           ...\n",
       "    7      2 -405.846809099 ...  175.02060759 1.60541232205 407.384285137\n",
       "    7      3 -390.553240191 ... 174.826719301 1.54995479057 392.150643917\n",
       "    8      0 -318.703000001 ...  215.31935812 1.54426530071 390.595017971\n",
       "    8      1 -331.075000003 ... 215.345026987  1.5999816328 405.886737469\n",
       "    8      2 -349.121000002 ...  211.02060759 1.60541232205 407.384285137\n",
       "    8      3 -336.747999992 ... 210.826719301 1.54995479057 392.150643917\n",
       "    9      0  -125.10482863 ...  251.31935812 1.54426530071 390.595017971\n",
       "    9      1  -129.83038525 ... 251.345026987  1.5999816328 405.886737469\n",
       "    9      2  -159.04283509 ...  247.02060759 1.60541232205 407.384285137\n",
       "    9      3  -154.31646944 ... 246.826719301 1.54995479057 392.150643917"
      ]
     },
     "execution_count": 238,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def is_on_gfa(tileid, ra, dec):\n",
    "    \"\"\"\n",
    "    Checks if a target is on any of the 10 GFAs given a tileid and an array of RA and Dec pointings. \n",
    "    Parameters\n",
    "    ----------\n",
    "    tileid: An integer representing the tile's unique ID which can retrieve its RA and Dec pointing\n",
    "    ra: An array of RA values for locations in the sky\n",
    "    dec: An array of declination values for locations in the sky\n",
    "    \"\"\"\n",
    "    import desimodel.footprint\n",
    "    import desimodel.io\n",
    "    # If any calculated area is under the threshold area, it is mathematically impossible\n",
    "    THRESHOLD_AREA = 469.7\n",
    "    TOLERANCE = 100\n",
    "    MIN_TOLERANCE = .001\n",
    "    # The area boundary's value is the area of the gfa plus some tolerance. \n",
    "    AREA_BOUNDARY = 469.8037221 + MIN_TOLERANCE + TOLERANCE\n",
    "    \n",
    "    targetindices = []\n",
    "    gfatable = desimodel.io.load_gfa()\n",
    "    telra, teldec = desimodel.footprint.get_tile_radec(tileid)\n",
    "    #print(telra, teldec)\n",
    "    targetx, targety = desimodel.focalplane.radec2xy(telra, teldec, ra, dec)\n",
    "    targetx = np.asarray(targetx)\n",
    "    targety = np.asarray(targety)\n",
    "    \n",
    "    # Method to check if point is inside the rectangle\n",
    "    for gfaid in range(0, 40, 4):\n",
    "        # a1 through a4 are edge lengths of the rectangle formed by corners of the GFAs\n",
    "        a1 = np.sqrt((gfatable['X'][gfaid] - gfatable['X'][gfaid + 1])**2\n",
    "                    + (gfatable['Y'][gfaid] - gfatable['Y'][gfaid + 1])**2)\n",
    "        a2 = np.sqrt((gfatable['X'][gfaid + 1] - gfatable['X'][gfaid + 2])**2\n",
    "                    + (gfatable['Y'][gfaid + 1] - gfatable['Y'][gfaid + 2])**2)\n",
    "        a3 = np.sqrt((gfatable['X'][gfaid + 2] - gfatable['X'][gfaid + 3])**2\n",
    "                    + (gfatable['Y'][gfaid + 2] - gfatable['Y'][gfaid + 3])**2)\n",
    "        a4 = np.sqrt((gfatable['X'][gfaid + 3] - gfatable['X'][gfaid])**2\n",
    "                    + (gfatable['Y'][gfaid + 3] - gfatable['Y'][gfaid])**2)\n",
    "        # b1 through b4 are the line segments from each corner to the target location\n",
    "        b1 = np.sqrt((gfatable['X'][gfaid] - targetx)**2\n",
    "                    + (gfatable['Y'][gfaid] - targety)**2)\n",
    "        b2 = np.sqrt((gfatable['X'][gfaid + 1] - targetx)**2\n",
    "                    + (gfatable['Y'][gfaid + 1] - targety)**2)\n",
    "        b3 = np.sqrt((gfatable['X'][gfaid + 2] - targetx)**2\n",
    "                    + (gfatable['Y'][gfaid + 2] - targety)**2)\n",
    "        b4 = np.sqrt((gfatable['X'][gfaid + 3] - targetx)**2\n",
    "                    + (gfatable['Y'][gfaid + 3] - targety)**2)\n",
    "        # Calculating areas of triangles using Heron's Formula\n",
    "        #gfaarea = a1 * a2\n",
    "        #gfaarea2 = a2 * a3\n",
    "        u1 = (a1 + b1 + b2) / 2.0\n",
    "        u2 = (a2 + b2 + b3) / 2.0\n",
    "        u3 = (a3 + b3 + b4) / 2.0\n",
    "        u4 = (a4 + b4 + b1) / 2.0\n",
    "        area1 = np.sqrt(u1 * (u1 - a1) * (u1 - b1) * (u1 - b2))\n",
    "        area2 = np.sqrt(u2 * (u2 - a2) * (u2 - b2) * (u2 - b3))\n",
    "        area3 = np.sqrt(u3 * (u3 - a3) * (u3 - b3) * (u3 - b4))\n",
    "        area4 = np.sqrt(u4 * (u4 - a4) * (u4 - b4) * (u4 - b1))\n",
    "        targetarea = area1 + area2 + area3 + area4\n",
    "        \n",
    "        assert np.all(targetarea > THRESHOLD_AREA)\n",
    "        #use np.where!!!!!!!!!\n",
    "        # This if statement may be redundant and thus unnecessary\n",
    "        if(any(targetarea < AREA_BOUNDARY) and any(targetarea > THRESHOLD_AREA)):\n",
    "            print(\"The target is on GFA\", int(gfaid / 4))  \n",
    "            newtargetindices = np.where(targetarea < AREA_BOUNDARY)\n",
    "            targetindices.append(newtargetindices[0])\n",
    "    \"\"\"targetindices is an array of targets with their respective indices in the \n",
    "    RA and Dec list passed in that fall on certain GFAs denoted by the index \n",
    "    in the targetindices array.\"\"\"\n",
    "    return targetindices\n",
    "        \n",
    "        \n",
    "#is_on_gfa(23658, 6.644525362152656, -9.055425745149217)\n",
    "import desimodel.focalplane\n",
    "#test point: (x,y) = (119, -380)\n",
    "#test tileid = 23658\n",
    "telra1, teldec1 = desimodel.footprint.get_tile_radec(23658)\n",
    "testra, testdec = desimodel.focalplane.xy2radec(telra1, teldec1, -335, -210)\n",
    "#is_on_gfa(23658, testra, testdec)\n",
    "\n",
    "#print(testra, testdec)\n",
    "testra2 = [testra, 6.644525362152656, 7]\n",
    "testdec2 = [testdec, -9.055425745149217, -10]\n",
    "is_on_gfa(23658, testra2, testdec2)\n",
    "\n",
    "gfatable"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 239,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7.11 -10.53\n",
      "The target is on GFA 0\n",
      "The target is on GFA 1\n",
      "The target is on GFA 2\n",
      "The target is on GFA 3\n",
      "The target is on GFA 4\n",
      "The target is on GFA 5\n",
      "The target is on GFA 6\n",
      "The target is on GFA 7\n",
      "The target is on GFA 8\n",
      "The target is on GFA 9\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "\"targetindices = is_on_gfa(23658, ra, dec)\\nprint(targetindices)\\n#print(tiles.dtype.names)\\nplt.plot(ra, dec, 'b.')\\nfor gfaid in targetindices:\\n    for indices in gfaid:\\n        plt.plot(ra[indices], dec[indices], 'g.')\\n        #ra.remove(ra[indices])\\n        #dec.remove(dec[indices])\""
      ]
     },
     "execution_count": 239,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "#telra and teldec = 7.11 -10.53\n",
    "\"\"\"mindec = -9.8\n",
    "maxdec = -9.6\n",
    "minra = 6.5\n",
    "maxra = 8.47\"\"\"\n",
    "mindec = -12.5\n",
    "maxdec = -7\n",
    "minra = 5.4\n",
    "maxra = 8.9\n",
    "ra = []\n",
    "dec = []\n",
    "import desimodel.focalplane\n",
    "import desimodel.footprint\n",
    "\n",
    "while mindec < maxdec:\n",
    "    startra = minra\n",
    "    while startra < maxra:\n",
    "        ra.append(startra)\n",
    "        dec.append(mindec)\n",
    "        startra += .05\n",
    "    mindec += .05\n",
    "\n",
    "import desimodel.io\n",
    "tiles = desimodel.io.load_tiles()\n",
    "ra = np.asarray(ra)\n",
    "dec = np.asarray(dec)\n",
    "\n",
    "\"\"\"index = 0\n",
    "for specific_tile in tiles:\n",
    "    if(tiles[index]['TILEID'] == 23658):\n",
    "        print(index)\n",
    "    index += 1\n",
    "#INDEX IS TILES[8254]\"\"\"\n",
    "points = desimodel.footprint.find_points_in_tiles(tiles[8254], ra, dec)\n",
    "\"\"\"for lists in points:\n",
    "    if lists:\n",
    "        print(lists)\"\"\"\n",
    "#print(points)\n",
    "new_ra = []\n",
    "new_dec = []\n",
    "for point_indices in points:\n",
    "    plt.plot(ra[point_indices], dec[point_indices], 'b.')\n",
    "    new_ra.append(ra[point_indices])\n",
    "    new_dec.append(dec[point_indices])\n",
    "targetindices = is_on_gfa(23658, new_ra, new_dec)\n",
    "for gfaid in targetindices:\n",
    "    for indices in gfaid:\n",
    "        plt.plot(new_ra[indices], new_dec[indices], 'g.')\n",
    "\"\"\"ra.append(8.452)\n",
    "dec.append(-9.6977)\n",
    "ra.append(6.644525362152656)\n",
    "dec.append(-9.055)\"\"\"\n",
    "\n",
    "\"\"\"targetindices = is_on_gfa(23658, ra, dec)\n",
    "print(targetindices)\n",
    "#print(tiles.dtype.names)\n",
    "plt.plot(ra, dec, 'b.')\n",
    "for gfaid in targetindices:\n",
    "    for indices in gfaid:\n",
    "        plt.plot(ra[indices], dec[indices], 'g.')\n",
    "        #ra.remove(ra[indices])\n",
    "        #dec.remove(dec[indices])\"\"\"\n",
    "#savefig('entire_gfa_tolerance_100.jpg')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.6.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}