code: def run_analysis(self, argv):
          
          args = self._parser.parse_args(argv)
  
          exttype = splitext(args.infile)[-1]
          if exttype in ['.fits', '.npy']:
              castro_data = CastroData.create_from_sedfile(args.infile)
          elif exttype in ['.yaml']:
              castro_data = CastroData.create_from_yamlfile(args.infile)
          else:
              raise ValueError("Can not read file type %s for SED" % extype)
  
          ylims = [1e-8, 1e-5]
  
          plot = plotCastro(castro_data, ylims)
          if args.outfile:
              plot[0].savefig(args.outfile)
docs: Run this analysis