code: def add_catalogue_cluster(self, catalogue, vcl, flagvector,
                                cluster_id=None, overlay=True):
          
          # Create simple magnitude scaled point basemap
          self.add_size_scaled_points(catalogue.data['longitude'],
                                      catalogue.data['latitude'],
                                      catalogue.data['magnitude'],
                                      shape="o",
                                      alpha=0.8,
                                      colour=(0.5, 0.5, 0.5),
                                      smin=1.0,
                                      sscale=1.5,
                                      overlay=True)
          # If cluster ID is not specified just show mainshocks
          if cluster_id is None:
              idx = flagvector == 0
              self.add_size_scaled_points(catalogue.data['longitude'][idx],
                                          catalogue.data['latitude'][idx],
                                          catalogue.data['magnitude'][idx],
                                          shape="o",
                                          colour="r",
                                          smin=1.0,
                                          sscale=1.5,
                                          overlay=overlay)
              return
          if not isinstance(cluster_id, collections.Iterable):
              cluster_id = [cluster_id]
          for iloc, clid in enumerate(cluster_id):
              if iloc == (len(cluster_id) - 1):
                  # On last iteration set overlay to function overlay
                  temp_overlay = overlay
              else:
                  temp_overlay = True
              idx = vcl == clid
              self.add_size_scaled_points(
                  catalogue.data["longitude"][idx],
                  catalogue.data["latitude"][idx],
                  catalogue.data["magnitude"][idx],
                  shape="o",
                  colour=DISSIMILAR_COLOURLIST[(iloc + 1) % NCOLS],
                  smin=1.0,
                  sscale=1.5,
                  overlay=temp_overlay)
docs: Creates a plot of a catalogue showing where particular clusters exist