code: def run(self, maxiter = 4, verbose = False):
          
              
          if self.satlevel > 0 and self.satstars == None:
              self.findsatstars(verbose=True)
              
          print "Starting %i L.A.Cosmic iterations ..." % maxiter
          for i in range(1, maxiter+1):
              print "Iteration %i" % i
              
              iterres = self.lacosmiciteration(verbose=verbose)
              print "%i cosmic pixels (%i new)" % (iterres["niter"], iterres["nnew"])
              
              #self.clean(mask = iterres["mask"]) # No, we want clean to operate on really clean pixels only !
              # Thus we always apply it on the full mask, as lacosmic does :
              self.clean(verbose=verbose)
              # But note that for huge cosmics, one might want to revise this.
              # Thats why I added a feature to skip saturated stars !
              
              if iterres["niter"] == 0:
                  break
docs: Full artillery :-)
              - Find saturated stars
              - Run maxiter L.A.Cosmic iterations (stops if no more cosmics are found)
      
          Stops if no cosmics are found or if maxiter is reached.