In [11]:

    

# Identify the list of GFAs with the 100 brightest objects on them
def findBrightestObjects(sweepname):
    from astropy.io import fits
    from astropy.table import Table
    sweep = Table.read(sweepname)
    import desimodel.io
    import desimodel.footprint
    tiles = desimodel.io.load_tiles()
    points = desimodel.footprint.find_points_in_tiles(tiles, sweep['RA'], sweep['DEC'], 1.651)
    counter = 0
    for lists in points:
        if lists:
            #test if on gfa
        counter += 1
sweepname = 'sweep-280p050-290p055.fits'
findBrightestObjects(sweepname)


    

    




[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]
[1]