docs: def inFootprint(footprint,ra,dec):
      """
      Check if set of ra,dec combinations are in footprint.
      Careful, input files must be in celestial coordinates.
      
      filename : Either healpix map or mangle polygon file
      ra,dec   : Celestial coordinates
  
      Returns:
      inside   : boolean array of coordinates in footprint
      """
      if footprint is None:
          return np.ones(len(ra),dtype=bool)
      
      try:
          if isinstance(footprint,str) and os.path.exists(footprint):
              filename = footprint
              #footprint = hp.read_map(filename,verbose=False)
              #footprint = fitsio.read(filename)['I'].ravel()
              footprint = read_map(filename)
          nside = hp.npix2nside(len(footprint))
          pix = ang2pix(nside,ra,dec)
          inside = (footprint[pix] > 0)
      except IOError:
          logger.warning("Failed to load healpix footprint; trying to use mangle...")
          inside = inMangle(filename,ra,dec)
      return inside