def calcResolution(self, partlist, stackfile, apix):
### group particles by refnum
reflistsdict = {}
for partdict in partlist:
refnum = partdict['template']
partnum = partdict['num']
if not refnum in reflistsdict:
reflistsdict[refnum] = []
reflistsdict[refnum].append(partnum)
### get resolution
self.resdict = {}
boxsizetuple = apFile.getBoxSize(stackfile)
boxsize = boxsizetuple[0]
for refnum in reflistsdict.keys():
partlist = reflistsdict[refnum]
esttime = 3e-6 * len(partlist) * boxsize**2
apDisplay.printMsg("Ref num %d; %d parts; est time %s"
%(refnum, len(partlist), apDisplay.timeString(esttime)))
frcdata = apFourier.spectralSNRStack(stackfile, apix, partlist, msg=False)
frcfile = "frcplot-%03d.dat"%(refnum)
apFourier.writeFrcPlot(frcfile, frcdata, apix, boxsize)
res = apFourier.getResolution(frcdata, apix, boxsize)
self.resdict[refnum] = res
return
def calcResolution(self, partlist, stackfile, apix):
### group particles by refnum
reflistsdict = {}
for partdict in partlist:
refnum = partdict['template']
partnum = partdict['num']
if not refnum in reflistsdict:
reflistsdict[refnum] = []
reflistsdict[refnum].append(partnum)
### get resolution
self.resdict = {}
boxsizetuple = apFile.getBoxSize(stackfile)
boxsize = boxsizetuple[0]
for refnum in reflistsdict.keys():
partlist = reflistsdict[refnum]
esttime = 3e-6 * len(partlist) * boxsize**2
apDisplay.printMsg("Ref num %d; %d parts; est time %s"
%(refnum, len(partlist), apDisplay.timeString(esttime)))
frcdata = apFourier.spectralSNRStack(stackfile, apix, partlist, msg=False)
frcfile = "frcplot-%03d.dat"%(refnum)
apFourier.writeFrcPlot(frcfile, frcdata, apix, boxsize)
res = apFourier.getResolution(frcdata, apix, boxsize)
self.resdict[refnum] = res
return