docs: def query(self, coords, **kwargs):
          """
          Returns E(B-V) (or a different Planck dust inference, depending on how
          the class was intialized) at the specified location(s) on the sky.
  
          Args:
              coords (:obj:`astropy.coordinates.SkyCoord`): The coordinates to query.
  
          Returns:
              A float array of the selected Planck component, at the given
              coordinates. The shape of the output is the same as the shape of the
              coordinates stored by ``coords``. If extragalactic E(B-V), tau_353
              or radiance was chosen, then the output has units of magnitudes of
              E(B-V). If the selected Planck component is temperature (or
              temperature error), then an :obj:`astropy.Quantity` is returned, with
              units of Kelvin. If beta (or beta error) was chosen, then the output
              is unitless.
          """
          return self._scale * super(PlanckQuery, self).query(coords, **kwargs)