# got some sdss sources from Vizier, a 12 minute radius with rmag<21
# about 773 sources found
sources.run(bench, open('/envoy5/raid2/Vim/data/benches/44892761645094959059548390744299/sources.votable').read())

# went and found 2mass in the registry and upload
textPositions.run(bench, open('/envoy5/raid2/Vim/data/benches/44892761645094959059548390744299/newresources.5041').read())

# run 2mass prox search at 2 arc seconds
proxsearch.run(bench, survey='2MASS-PSC', sr=0.000555555555556, tablename='_2mass', verb='1', maxthreads='10')

# merge with the sdss
join.run(bench, tablename='_2mass', match='best')

# select rows with non-null j magnitude (i.e. thos in 2mass)
# we get abouyt 148 remaining
selectrows.run(bench, tablename='sources', criterion='! NULL__2mass_j_m')

# compute 2mass-J - SDSS-r
addcol.run(bench, tablename='sources', colname='JminusR', expression='_2mass_j_m - rmag')

# compute 2mass-H - SDSS-z
addcol.run(bench, tablename='sources', colname='Hminusz', expression='_2mass_h_m - zmag')

# a couple of wierd outliers
selectrows.run(bench, tablename='sources', criterion='Hminusz > -5')
selectrows.run(bench, tablename='sources', criterion='Hminusz < 1')

# then download sources.votable from Utility/ViewAllFiles, and plot with topcat