def p7anal(event=None, filedir='.', topdir=None, idl=False, islak=False, isinteractive=True): import p7anal as p7 #reload(p7) #global numevents, nummodels, isinteractive if event == None: event = p6Restore(filedir, topdir, idl) cwd = os.getcwd().split("/") if cwd[-1] == event[0].modeldir: os.chdir('..') printout = printoutput.init(event[0].params.printout, event) numevents = len(event) nummodels = np.zeros(numevents, dtype=int) for j in range(numevents): nummodels[j] = len(event[j].params.model) for j in range(numevents): print("\n" + event[j].eventname, file=printout) for i in range(nummodels.min()): print("\nCurrent model = " + str(event[j].params.model[i]), file=printout) p7.stdanal(event[j], event[j].fit[i], j * nummodels.min() + i, printout, islak) p7Save(event[j]) printoutput.close(printout) if isinteractive == False: plt.close('all') else: plt.show() return
def p7anal(events=None, directory='./', idl=False, islak=False, isinteractive=isinteractive): if events == None: events = p6Restore(directory, idl) nummodels = np.array([len(event.params.model) for event in events], dtype=int) printout = po.init(events[0].params.printout, events) for j, event in enumerate(events): print("\n" + event.eventname, file=printout) minnum = nummodels.min() for i in range(minnum): print("\nCurrent model = " + str(event.params.model[i]), file=printout) p7.stdanal(event, event.fit[i], j*minnum+i, printout, islak) p7Save(event, directory) po.close(printout) if isinteractive == False: plt.close('all') else: plt.show() return