def plot_tnos(ax, name=True, ossos=True):
# plot the star of the show first
plt.errorbar(740, 50, xerr=70,
fmt='*', ms=12,
mfc='r', mec=palette[3], ecolor='r', elinewidth=1, capsize=0, zorder=0,
label='_nolegend_')
plt.text(630, 50.5, 'L91', fontsize='small')
tnos = Table.read('/Users/bannisterm/Dropbox/Papers_in_progress/OSSOS/high_q_uo3l91/data/all_tnos_no_centaurs.csv',
format='csv', data_start=1)
plt.errorbar(tnos['a'], tnos['q'],
xerr=tnos['sigma_a'],
fmt='o', ms=6,
mfc=palette[3], mec=palette[3], ecolor=palette[3], capsize=0, elinewidth=0.5, zorder=0,
label='_nolegend_')
if name:
for id in ['GP136', 'Sedna', 'VP113', 'GB174', 'SR349', 'VN112', 'CR105', 'FT28', 'FE72']:
for tno in tnos:
if tno['full_name'].__contains__(id):
xoffset = 15
yoffset = 0.5
if id in ['SR349']: # pesky overlapping labels
xoffset = -55
if id in ['VN112']:
yoffset = -0.2
plt.text(tno['a']+xoffset, tno['q']+yoffset, id, fontsize='xx-small')
if ossos:
ossos_tnos = parsers.ossos_release_parser(table=True)
ossos_tnos['q'] = ossos_tnos['a'] * (1. - ossos_tnos['e'])
ossos_tnos['e_q'] = (ossos_tnos['e_a'] / ossos_tnos['a']) + (ossos_tnos['e_e'] / ossos_tnos['e'])
plt.errorbar(ossos_tnos['a'], ossos_tnos['q'],
xerr=ossos_tnos['e_a'],
yerr=ossos_tnos['e_q'],
fmt='o', ms=6,
capsize=0, elinewidth=0.5, mfc=palette[2],
mec=palette[3], ecolor=palette[2], zorder=1,
label='_nolegend_'
)
for tno in ossos_tnos:
if tno['object'] in ['o5t52', 'o5m52', 'o5p060', 'o3l83', 'o5m72']:
plt.text(tno['a'] + 5, tno['q'] + 0.3, tno['object'], fontsize='xx-small')
return
def plot_tnos(ax, name=True, ossos=True):
# plot the star of the show first
plt.errorbar(740, 50, xerr=70,
fmt='*', ms=12,
mfc='r', mec=palette[3], ecolor='r', elinewidth=1, capsize=0, zorder=0,
label='_nolegend_')
plt.text(630, 50.5, 'L91', fontsize='small')
tnos = Table.read('/Users/bannisterm/Dropbox/Papers_in_progress/OSSOS/high_q_uo3l91/data/all_tnos_no_centaurs.csv',
format='csv', data_start=1)
plt.errorbar(tnos['a'], tnos['q'],
xerr=tnos['sigma_a'],
fmt='o', ms=6,
mfc=palette[3], mec=palette[3], ecolor=palette[3], capsize=0, elinewidth=0.5, zorder=0,
label='_nolegend_')
if name:
for id in ['GP136', 'Sedna', 'VP113', 'GB174', 'SR349', 'VN112', 'CR105', 'FT28', 'FE72']:
for tno in tnos:
if tno['full_name'].__contains__(id):
xoffset = 15
yoffset = 0.5
if id in ['SR349']: # pesky overlapping labels
xoffset = -55
if id in ['VN112']:
yoffset = -0.2
plt.text(tno['a']+xoffset, tno['q']+yoffset, id, fontsize='xx-small')
if ossos:
ossos_tnos = parsers.ossos_release_parser(table=True)
ossos_tnos['q'] = ossos_tnos['a'] * (1. - ossos_tnos['e'])
ossos_tnos['e_q'] = (ossos_tnos['e_a'] / ossos_tnos['a']) + (ossos_tnos['e_e'] / ossos_tnos['e'])
plt.errorbar(ossos_tnos['a'], ossos_tnos['q'],
xerr=ossos_tnos['e_a'],
yerr=ossos_tnos['e_q'],
fmt='o', ms=6,
capsize=0, elinewidth=0.5, mfc=palette[2],
mec=palette[3], ecolor=palette[2], zorder=1,
label='_nolegend_'
)
for tno in ossos_tnos:
if tno['object'] in ['o5t52', 'o5m52', 'o5p060', 'o3l83', 'o5m72']:
plt.text(tno['a'] + 5, tno['q'] + 0.3, tno['object'], fontsize='xx-small')
return
"specify blocks to be plotted, e.g. 13AE. Without specifying, will do all eight blocks,"
"making eight separate plots.")
args = parser.parse_args()
print(args)
file_id = ""
if args.blocks: # Coordinates at discovery-opposition of currently-determined OSSOS blocks (2013: 4; 2015: 8
# total).
blocks = {}
for b in args.blocks:
assert b in parameters.BLOCKS
blocks[b] = parameters.BLOCKS[b]
else: # do them all!
blocks = parameters.BLOCKS
discoveries = parsers.ossos_release_parser(table=True)
# extent = plot_extents['13A'] # hardwired for doing both at once
for blockname, block in list(blocks.items()):
extent = plot_extents[blockname]
if args.opposition:
# 13AE is a month different in opposition from 13AO, but 13B are both in October
date = parameters.OPPOSITION_DATES[blockname]
file_id = 'opposition-'
if args.discovery:
# E block discovery triplets are April 4,9,few 19; O block are May 7,8.
# BL is 09-29 AND 10-31. Discoveries in 10-31 are placed at predicted locations, to accommodate split.
date = parameters.DISCOVERY_DATES[blockname]
ax[2].set_ylabel('eccentricity')
ax[2].grid(True, alpha=grid_alpha)
# ax[2].set_xlabel('semimajor axis (AU)')
plt.xlabel('semimajor axis (AU)')
# ax[3].grid(True, alpha=grid_alpha)
# ax[3].set_xlabel('perihelion (AU)')
# ax[3].set_ylabel('arg. peri. (degrees)')
plt.draw()
# ofile = 'OSSOS+CFEPS+NGVS_aei_argperi.pdf'
outfile = 'OSSOS_aei_{}.pdf'.format(parameters.RELEASE_VERSION)
plt.savefig(outfile, transparent=True, bbox_inches='tight')
if __name__ == '__main__':
tnos = parsers.ossos_release_parser(
table=True) # return as an astropy.Table.Table
mag = tnos[numpy.where(tnos['mag'] < 23.5)]
ossos = []
for t in mag:
if t['object'].startswith('o'):
ossos.append(t)
print(len(ossos))
sys.exit()
# argperi_a(parameters.REAL_KBO_AST_DIR)
# fov = ['medium', 'classical', 'thousand_au']
# full_aei(tnos, fov[0])
# classicals_qi(tnos)
def main():
# parsers.output_discoveries_for_animation()
discoveries = parsers.ossos_release_parser(table=True)
ossos_sequence(discoveries)
ax[2].set_ylabel('eccentricity')
ax[2].grid(True, alpha=grid_alpha)
# ax[2].set_xlabel('semimajor axis (AU)')
plt.xlabel('semimajor axis (AU)')
# ax[3].grid(True, alpha=grid_alpha)
# ax[3].set_xlabel('perihelion (AU)')
# ax[3].set_ylabel('arg. peri. (degrees)')
plt.draw()
# ofile = 'OSSOS+CFEPS+NGVS_aei_argperi.pdf'
outfile = 'OSSOS_aei_{}.pdf'.format(parameters.RELEASE_VERSION)
plt.savefig(outfile, transparent=True, bbox_inches='tight')
if __name__ == '__main__':
tnos = parsers.ossos_release_parser(table=True) # return as an astropy.Table.Table
mag = tnos[numpy.where(tnos['mag'] < 23.5)]
ossos = []
for t in mag:
if t['object'].startswith('o'):
ossos.append(t)
print len(ossos)
sys.exit()
# argperi_a(parameters.REAL_KBO_AST_DIR)
# fov = ['medium', 'classical', 'thousand_au']
# full_aei(tnos, fov[0])
# classicals_qi(tnos)
