docs: def bands(self):
          """
          Bandpasses for which StarPopulation has magnitude data
          """
          bands = []
          for c in self.stars.columns:
              if re.search('_mag',c):
                  bands.append(c)
          return bands