In [2]:
import numpy as np
from astropy.table import Table, join
import os
from matplotlib import pyplot as plt
In [4]:
old = Table.read("../lofar_pw_pdf.fits")
In [5]:
new = Table.read("../lofar_pw_corSigma.fits")
In [6]:
len(old)
Out[6]:
318769
In [7]:
len(new)
Out[7]:
325694
In [8]:
old.colnames
Out[8]:
['Source_Name',
'RA',
'E_RA',
'DEC',
'E_DEC',
'Peak_flux',
'E_Peak_flux',
'Total_flux',
'E_Total_flux',
'Maj',
'E_Maj',
'Min',
'E_Min',
'PA',
'E_PA',
'Isl_rms',
'S_Code',
'Mosaic_ID',
'AllWISE',
'objID',
'ra',
'dec',
'raErr',
'decErr',
'W1mag',
'W1magErr',
'i',
'iErr',
'colour',
'category',
'lr',
'lr_dist']
In [9]:
pwl = join(old, new,
join_type='left',
keys='Source_Name',
uniq_col_name='{col_name}{table_name}',
table_names=['_old', '_new'])
In [10]:
A = pwl["lr_new"]-pwl["lr_old"]
plt.hist(A[~np.isnan(A)], bins=50)
Out[10]:
(array([1.00000e+00, 0.00000e+00, 0.00000e+00, 1.00000e+00, 0.00000e+00,
1.00000e+00, 2.00000e+00, 5.00000e+00, 2.00000e+00, 1.30000e+01,
6.00000e+00, 7.00000e+00, 1.80000e+01, 3.70000e+01, 5.70000e+01,
8.90000e+01, 1.35000e+02, 2.16000e+02, 3.89000e+02, 7.63000e+02,
1.32600e+03, 2.74000e+03, 6.54800e+03, 1.66990e+04, 5.02190e+04,
1.87843e+05, 1.12860e+04, 5.62400e+03, 3.65700e+03, 2.55500e+03,
1.78100e+03, 1.24700e+03, 9.87000e+02, 6.21000e+02, 4.85000e+02,
3.42000e+02, 3.08000e+02, 2.37000e+02, 2.07000e+02, 1.57000e+02,
9.60000e+01, 9.30000e+01, 7.60000e+01, 8.00000e+01, 6.50000e+01,
3.90000e+01, 2.20000e+01, 1.90000e+01, 1.20000e+01, 6.00000e+00]),
array([-3736.78453907, -3590.04282495, -3443.30111083, -3296.5593967 ,
-3149.81768258, -3003.07596845, -2856.33425433, -2709.5925402 ,
-2562.85082608, -2416.10911196, -2269.36739783, -2122.62568371,
-1975.88396958, -1829.14225546, -1682.40054134, -1535.65882721,
-1388.91711309, -1242.17539896, -1095.43368484, -948.69197071,
-801.95025659, -655.20854247, -508.46682834, -361.72511422,
-214.98340009, -68.24168597, 78.50002815, 225.24174228,
371.9834564 , 518.72517053, 665.46688465, 812.20859878,
958.9503129 , 1105.69202702, 1252.43374115, 1399.17545527,
1545.9171694 , 1692.65888352, 1839.40059764, 1986.14231177,
2132.88402589, 2279.62574002, 2426.36745414, 2573.10916827,
2719.85088239, 2866.59259651, 3013.33431064, 3160.07602476,
3306.81773889, 3453.55945301, 3600.30116713]),
<a list of 50 Patch objects>)
In [11]:
cond = (pwl['ra_old'] == pwl['ra_new']) & (pwl['dec_old'] == pwl['dec_new'])
In [12]:
# Galaxies with the same id (or lack of id)
np.sum(cond)
Out[12]:
304332
In [13]:
len(pwl) - np.sum(cond)
Out[13]:
14439
In [14]:
pwl[cond & np.isnan(pwl["lr_new"])]
Out[14]:
<Table masked=True length=12679>
Source_Name RA_old E_RA_old DEC_old E_DEC_old Peak_flux_old E_Peak_flux_old Total_flux_old E_Total_flux_old Maj_old E_Maj_old Min_old E_Min_old PA_old E_PA_old Isl_rms_old S_Code_old Mosaic_ID_old AllWISE_old objID_old ra_old dec_old raErr_old decErr_old W1mag_old W1magErr_old i_old iErr_old colour_old category_old lr_old lr_dist_old RA_new E_RA_new DEC_new E_DEC_new Peak_flux_new E_Peak_flux_new Total_flux_new E_Total_flux_new Maj_new E_Maj_new Min_new E_Min_new PA_new E_PA_new Isl_rms_new S_Code_new Mosaic_ID_new AllWISE_new objID_new ra_new dec_new raErr_new decErr_new W1mag_new W1magErr_new i_new iErr_new colour_new category_new lr_new lr_dist_new deg arcsec deg arcsec mJy / beam mJy / beam mJy mJy arcsec arcsec arcsec arcsec deg deg mJy / beam str22 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 str1 str11 str19 int64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 str1 str8 str19 int64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 float64 ILTJ104358.55+513354.6 160.99396170310192 0.3804751732950043 51.56517118870972 0.35519395523884045 0.7504358543374862 0.11003057760648959 0.7238821644431407 0.1929865180123449 6.052143496251171 0.9112500467753947 5.741264073240394 0.8197450714237257 112.23586180232195 114.35902138229292 0.11110970808658749 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 160.99396170310192 0.3804751732950043 51.56517118870972 0.35519395523884045 0.7504358543374862 0.11003057760648959 0.7238821644431407 0.1929865180123449 6.052143496251171 0.9112500467753947 5.741264073240394 0.8197450714237257 112.23586180232195 114.35902138229292 0.11110970808658749 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104406.52+513113.1 161.02716295073938 0.4011857454685031 51.52030892960141 0.2956757648741586 0.7270173923847574 0.10214151036916275 0.6074119548975849 0.1909636448440653 6.268112922677453 1.0178900114002445 4.801343377156147 0.5835977430401842 60.57505752252567 21.443186110238656 0.10890992416534573 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.02716295073938 0.4011857454685031 51.52030892960141 0.2956757648741586 0.7270173923847574 0.10214151036916275 0.6074119548975849 0.1909636448440653 6.268112922677453 1.0178900114002445 4.801343377156147 0.5835977430401842 60.57505752252567 21.443186110238656 0.10890992416534573 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104418.69+522635.8 161.07788017007587 0.1392443079737803 52.44328024264552 0.15365161017254791 1.8803603269742886 0.10945499310990824 1.8244981225922174 0.1916173331372365 6.083887081701809 0.36475149468082985 5.744506125661684 0.3246577502745283 13.250364279525344 39.6429350698305 0.11037105286959559 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.07788017007587 0.1392443079737803 52.44328024264552 0.15365161017254791 1.8803603269742886 0.10945499310990824 1.8244981225922174 0.1916173331372365 6.083887081701809 0.36475149468082985 5.744506125661684 0.3246577502745283 13.250364279525344 39.6429350698305 0.11037105286959559 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104435.58+523504.2 161.1482631171491 0.08568213610331692 52.58450341354888 0.0830429609908379 3.208150878287284 0.10862472207765825 3.065275865605974 0.19105289251250496 5.914684271250655 0.2018966289669094 5.818499590408209 0.19541653134798362 79.7336792381118 84.64457548122874 0.10991691669914871 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.1482631171491 0.08568213610331692 52.58450341354888 0.0830429609908379 3.208150878287284 0.10862472207765825 3.065275865605974 0.19105289251250496 5.914684271250655 0.2018966289669094 5.818499590408209 0.19541653134798362 79.7336792381118 84.64457548122874 0.10991691669914871 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104436.25+524735.2 161.15103385405467 1.2199617729336236 52.79312678166211 0.8644133872023749 0.3540097948942407 0.11011454494707756 0.5677568296586294 0.16598952916402196 8.661663562496447 2.9797587427210432 6.669417101265727 1.8757150060463206 107.53881858862735 53.67163232115898 0.10323095921194181 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.15103385405467 1.2199617729336236 52.79312678166211 0.8644133872023749 0.3540097948942407 0.11011454494707756 0.5677568296586294 0.16598952916402196 8.661663562496447 2.9797587427210432 6.669417101265727 1.8757150060463206 107.53881858862735 53.67163232115898 0.10323095921194181 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104444.07+511936.1 161.18360938477204 0.24009202735296814 51.326707066776294 0.27964659773273015 1.2490819284670485 0.11463794713995668 1.4963173845205207 0.18778394744352095 7.047989575236918 0.6893780970359981 6.122431132096382 0.5272378282351275 150.46446886540957 26.18714403351791 0.11065466242143884 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.18360938477204 0.24009202735296814 51.326707066776294 0.27964659773273015 1.2490819284670485 0.11463794713995668 1.4963173845205207 0.18778394744352095 7.047989575236918 0.6893780970359981 6.122431132096382 0.5272378282351275 150.46446886540957 26.18714403351791 0.11065466242143884 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104459.26+522548.7 161.24691576951935 0.15319821110047088 52.43019508431113 0.12370772586806064 2.078805373091584 0.11003620345428822 2.1631742131679084 0.1884511595995359 6.468146652899722 0.36115852258770487 5.7943461124981095 0.2907568683985586 92.4561145877511 18.4087689539321 0.10922510409727693 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.24691576951935 0.15319821110047088 52.43019508431113 0.12370772586806064 2.078805373091584 0.11003620345428822 2.1631742131679084 0.1884511595995359 6.468146652899722 0.36115852258770487 5.7943461124981095 0.2907568683985586 92.4561145877511 18.4087689539321 0.10922510409727693 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104505.08+530140.3 161.27118146589584 1.601209138617367 53.0278699329833 2.4242716705861094 0.3674352418591535 0.06448618761641643 7.117435452768374 0.06773581332541503 nan nan nan nan nan nan 0.10877009481191635 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.27118146589584 1.601209138617367 53.0278699329833 2.4242716705861094 0.3674352418591535 0.06448618761641643 7.117435452768374 0.06773581332541503 31.952561397523215 5.717130991109009 21.83663454888869 3.759109988309131 0.752301043472845 23.92850660622191 0.10877009481191635 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ104509.34+514531.1 161.28890899746335 0.3215004244592514 51.75865000558583 0.3087822159947033 0.7099355545937909 0.09271305269034723 0.5981837795488666 0.17132508811718583 6.001053684442445 0.8599360027680633 5.057016646140941 0.6018657789493743 46.16780908101893 31.05828598601022 0.09783437417354435 S P4Hetdex16 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 161.28890899746335 0.3215004244592514 51.75865000558583 0.3087822159947033 0.7099355545937909 0.09271305269034723 0.5981837795488666 0.17132508811718583 6.001053684442445 0.8599360027680633 5.057016646140941 0.6018657789493743 46.16780908101893 31.05828598601022 0.09783437417354435 S P4Hetdex N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ILTJ152029.96+495917.8 230.12484821403467 0.03722151202149774 49.98829606972382 0.06506072273737616 8.682316929307623 0.12054750550305471 15.149193800545861 0.28763633993271837 8.749106824614495 0.15795677661758914 5.965037482290683 0.07876764057836 76.0001172927061 2.001076903456291 0.12054750550305471 M P227+50 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.12484821403467 0.03722151202149774 49.98829606972382 0.06506072273737616 8.682316929307623 0.12054750550305471 15.149193800545861 0.28763633993271837 8.749106824614495 0.15795677661758914 5.965037482290683 0.07876764057836 76.0001172927061 2.001076903456291 0.12054750550305471 M P227+50 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152053.26+534624.4 230.22191225318585 0.6145211744323705 53.773453742727135 0.4788175314375589 1.3602823602543328 0.2397597629766541 1.8583078292366353 0.37867980106313837 7.91180361112263 1.5440828431314768 6.219759032679513 0.9900859453991127 64.9979317875435 35.76486982396624 0.22824716870673 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.22191225318585 0.6145211744323705 53.773453742727135 0.4788175314375589 1.3602823602543328 0.2397597629766541 1.8583078292366353 0.37867980106313837 7.91180361112263 1.5440828431314768 6.219759032679513 0.9900859453991127 64.9979317875435 35.76486982396624 0.22824716870673 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152129.05+515755.4 230.37105039987154 0.025647087724637 51.96539795833628 0.04288014871828056 17.335821453611853 0.21187863603699952 21.18741786497023 0.48045551988396806 7.0827176947700385 0.10099830969432191 5.448543696816709 0.060355343450877265 90.58148999838481 2.2013324017455536 0.21187863603699952 M P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.37105039987154 0.025647087724637 51.96539795833628 0.04288014871828056 17.335821453611853 0.21187863603699952 21.18741786497023 0.48045551988396806 7.0827176947700385 0.10099830969432191 5.448543696816709 0.060355343450877265 90.58148999838481 2.2013324017455536 0.21187863603699952 M P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152138.26+520831.4 230.40942772072216 0.2876561755724256 52.14207737034642 0.19003551227481258 1.7000159560797052 0.16478288589794218 1.4825424052304481 0.3022180511716143 6.2922282081368905 0.6899083013735393 4.991657902092878 0.42773632024198793 106.10673820373907 20.89622609414001 0.17283174383919686 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.40942772072216 0.2876561755724256 52.14207737034642 0.19003551227481258 1.7000159560797052 0.16478288589794218 1.4825424052304481 0.3022180511716143 6.2922282081368905 0.6899083013735393 4.991657902092878 0.42773632024198793 106.10673820373907 20.89622609414001 0.17283174383919686 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152143.21+520106.4 230.43004053429638 0.5368960156911345 52.01846426571062 0.4613945714089374 5.143750660809247 0.9123487645832231 5.654568489739902 1.5452412617900468 7.107430204151674 1.4159681920356553 5.5707355958286975 0.8796980166224263 127.23608538297238 35.3340175742053 0.9024523315019906 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.43004053429638 0.5368960156911345 52.01846426571062 0.4613945714089374 5.143750660809247 0.9123487645832231 5.654568489739902 1.5452412617900468 7.107430204151674 1.4159681920356553 5.5707355958286975 0.8796980166224263 127.23608538297238 35.3340175742053 0.9024523315019906 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152154.98+511925.1 230.4790780474188 0.19593362179459206 51.32366165729939 0.23063508724146933 2.349538681872924 0.18072055602370704 2.7621324178625897 0.2973111788475373 6.790805945467852 0.5434531804998639 6.236394104978554 0.46104423540206063 3.4206139792392594 40.348336172848 0.1746469351928681 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.4790780474188 0.19593362179459206 51.32366165729939 0.23063508724146933 2.349538681872924 0.18072055602370704 2.7621324178625897 0.2973111788475373 6.790805945467852 0.5434531804998639 6.236394104978554 0.46104423540206063 3.4206139792392594 40.348336172848 0.1746469351928681 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152157.74+524715.2 230.49058845820875 0.643207077538589 52.787570266293606 0.8588062676925968 0.7124078161691768 0.1815158055864309 0.9330066833732312 0.2893476067704043 7.422524898829871 2.0238277522084624 6.3548237272972194 1.5130496950761525 1.291509683295203 74.34904271920153 0.17277093138545752 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.49058845820875 0.643207077538589 52.787570266293606 0.8588062676925968 0.7124078161691768 0.1815158055864309 0.9330066833732312 0.2893476067704043 7.422524898829871 2.0238277522084624 6.3548237272972194 1.5130496950761525 1.291509683295203 74.34904271920153 0.17277093138545752 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152240.73+525728.8 230.66970948634108 1.7006681337892429 52.958001046417834 0.9610494654018784 0.7324545221121065 0.18003331337788836 2.289819815390047 0.23308638262876863 14.86681935851982 4.320545653000943 7.5741541986764735 1.5745556939920569 116.13082467588823 25.12824128707611 0.1816209842218086 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.66970948634108 1.7006681337892429 52.958001046417834 0.9610494654018784 0.7324545221121065 0.18003331337788836 2.289819815390047 0.23308638262876863 14.86681935851982 4.320545653000943 7.5741541986764735 1.5745556939920569 116.13082467588823 25.12824128707611 0.1816209842218086 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152337.32+523922.7 230.90550269444554 0.17370906033346753 52.656307196951765 0.15638177004415502 3.222109364760001 0.20299577806930055 3.3891084586057296 0.3459729485049387 6.346237865424089 0.4116454173654634 5.9701601936769695 0.3653192568938889 106.81620972635247 46.03068876041295 0.20068121375516057 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 230.90550269444554 0.17370906033346753 52.656307196951765 0.15638177004415502 3.222109364760001 0.20299577806930055 3.3891084586057296 0.3459729485049387 6.346237865424089 0.4116454173654634 5.9701601936769695 0.3653192568938889 106.81620972635247 46.03068876041295 0.20068121375516057 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan ILTJ152419.22+523536.1 231.0800664187888 0.14067545366080925 52.59336774314981 0.11938343216287275 4.461245306752376 0.21571539424938707 4.916950454229268 0.3623642833263276 6.604212351624555 0.333671267857945 6.011786561646778 0.27819918126441257 79.93510531121898 24.68290977990122 0.2112516958732158 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan 231.0800664187888 0.14067545366080925 52.59336774314981 0.11938343216287275 4.461245306752376 0.21571539424938707 4.916950454229268 0.3623642833263276 6.604212351624555 0.333671267857945 6.011786561646778 0.27819918126441257 79.93510531121898 24.68290977990122 0.2112516958732158 S P227+53 N/A 999999 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 1e+20 nan nan
In [ ]:
Content source: nudomarinero/mltier1
Similar notebooks: