code: def match_default_fit(imglist, reference_catalog):
      
      log.info("-------------------- STEP 5b: (match_default_fit) Cross matching and fitting ---------------------------")
      # Specify matching algorithm to use
      match = tweakwcs.TPMatch(searchrad=250, separation=0.1,
                               tolerance=100, use2dhist=False)
      # Align images and correct WCS
      tweakwcs.align_wcs(imglist, reference_catalog, match=match, expand_refcat=False) #TODO: turn on 'expand_refcat' option in future development
  
      # Interpret RMS values from tweakwcs
      interpret_fit_rms(imglist, reference_catalog)
  
      return imglist
docs: Perform cross-matching and final fit using 2dHistogram matching
  
      Parameters
      ----------
      imglist : list
          List of input image `~tweakwcs.tpwcs.FITSWCS` objects with metadata and source catalogs
  
      reference_catalog : Table
          Astropy Table of reference sources for this field
  
      Returns
      --------
      imglist : list
          List of input image `~tweakwcs.tpwcs.FITSWCS` objects with metadata and source catalogs