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