code: def run_analysis(self, argv):
          
          args = self._parser.parse_args(argv)
  
          if not HAVE_ST:
              raise RuntimeError(
                  "Trying to run fermipy analysis, but don't have ST")
  
          gta = GTAnalysis(args.config, logging={'verbosity': 3},
                           fileio={'workdir_regex': '\.xml$|\.npy$'})
  
          gta.setup(overwrite=False)
  
          baseline_roi_fit(gta, make_plots=args.make_plots,
                           minmax_npred=[1e3, np.inf])
  
          localize_sources(gta, nstep=5, dtheta_max=0.5, update=True,
                           prefix='base', make_plots=args.make_plots)
  
          gta.find_sources(sqrt_ts_threshold=5.0, search_skydir=gta.roi.skydir,
                           search_minmax_radius=[1.0, np.nan])
          gta.optimize()
          gta.print_roi()
          gta.print_params()
  
          gta.free_sources(skydir=gta.roi.skydir, distance=1.0, pars='norm')
          gta.fit(covar=True)
          gta.print_roi()
          gta.print_params()
  
          gta.write_roi(args.roi_baseline, make_plots=args.make_plots)
docs: Run this analysis