def run(self, maxiter = 4, verbose = False):
        """
        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.
        """
            
        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