ObservationGroupAxis('ZEN_PNT', zenith, fmt='edges'), ObservationGroupAxis('N_TELS', ntels, fmt='values'), ]) obs_table zenith = Angle([0, 30, 40, 50], 'deg') ntels = [3, 4] obs_groups = ObservationGroups([ ObservationGroupAxis('ZEN_PNT', zenith, fmt='edges'), ObservationGroupAxis('N_TELS', ntels, fmt='values'), ]) obs1=obs_groups.apply(datastore.obs_table) obs1 print(obs1) print(datastore.obs_table) obs_gr=obs_groups.get_group_of_observations(obs1,2) obs_gr datastore datastore.obs_table datastore.obs_table.Column datastore.obs_table.field crab=SkyCoord.from_name("Crab") crab crab_table=datastore.obs_table(SkyCoord.separation(crab)) crab_table=datastore.obs_table(SkyCoord.separation(crab,2.0)) get_ipython().magic(u'pinfo SkyCoord.separation') SkyCoord.separation(datastore.obs_table,crab) SkyCoord.separation(crab,datastore.obs_table) geminga=SkyCoord.from_name("geminga") SkyCoord.separation(geminga,crab) get_ipython().magic(u'pinfo datastore.obs_table.group_by')