def analyze_prior_results(self, analysis_source): """ Runs analysis of previous grid search / integration results, used in an analyze-only mode """ from prime.iota.iota_analysis import Analyzer from libtbx import easy_pickle as ep if os.path.isdir(analysis_source): int_folder = os.path.abspath(analysis_source) else: try: int_folder = os.path.abspath(os.path.join(os.curdir, 'integration/{}/image_objects'.format(analysis_source))) except ValueError: print 'Run #{} not found'.format(analysis_source) if os.path.isdir(int_folder): int_list = [os.path.join(int_folder, i) for i in os.listdir(int_folder)] img_objects = [ep.load(i) for i in int_list] analysis = Analyzer(img_objects, None, self.iver, self.now) analysis.print_results() analysis.unit_cell_analysis(self.params.analysis.cluster_threshold, int_folder, False) analysis.print_summary(None) analysis.show_heatmap() else: print 'No results found in {}'.format(int_folder)
misc.iota_exit(iota_version, True) else: if os.path.isfile('iota.ini'): init = ep.load('iota.ini') else: print 'This run not initialized. Use "--mpi init" option' misc.iota_exit(iota_version, True) if 'ana' in args.mpi: img_objects = [ep.load(os.path.join(init.obj_base, i)) for i in os.listdir(init.obj_base)] int_objects = [i for i in img_objects if i.fail == None] if len(int_objects) != 0: # Analysis of integration results analysis = Analyzer(int_objects, init.logfile, iota_version, init.now) analysis.print_results() analysis.unit_cell_analysis(init.params.analysis.cluster_threshold, init.int_base) analysis.print_summary(init.int_base) analysis.make_prime_input(init.int_base) else: print "No images integrated" else: if 'imp' in args.mpi: inp_list = init.input_list msg = "Importing {} images".format(len(inp_list)) elif 'gri' in args.mpi: obj_list = [ep.load(os.path.join(init.obj_base, i)) for i in os.listdir(init.obj_base)]