code: def set_maxrad(self,maxrad, distribution_skip=True):
          
          self.maxrad = maxrad
          self.apply_constraint(UpperLimit(self.Rsky,maxrad,
                                           name='Max Rsky'),
                                overwrite=True,
                                distribution_skip=distribution_skip)
docs: Adds a constraint that rejects everything with Rsky > maxrad
  
          Requires ``Rsky`` attribute, which should always have units.
  
          :param maxrad:
              The maximum angular value of Rsky.
          :type maxrad:
              :class:`astropy.units.Quantity`
  
          :param distribution_skip:
              This is by default ``True``.  *To be honest, I'm not
              exactly sure why.  Might be important, might not
              (don't remember).*