def plotDensHistPlots(self,
getVoxelStats = False,
perAtmDensHist = False):
# histogram & kde plots of number of voxels per atom
self.startTimer()
self.printStepNumber()
self.lgwrite(ln = 'Plotting histogram plots of voxels per atom...')
self.lgwrite(ln = 'Plots written to "{}plots"'.format(self.filesOut))
stats = plotVxlsPerAtm(pdbName = self.pdbName,
where = self.filesOut,
vxlsPerAtom = self.vxlsPerAtom,
plotType = 'both',
returnStats = getVoxelStats)
if stats != '':
print 'mean: {}\nstd: {}\nmax: {}\nmin: {}'.format(*stats)
if perAtmDensHist:
plotDensForAtm(pdbName = self.pdbName,
where = self.filesOut,
vxlsPerAtom = self.vxlsPerAtom,
plotType = 'both',
PDBarray = self.PDBarray)
self.stopTimer()
def plotDensHistPlots(self):
# histogram & kde plots of number of voxels per atom
for plotType in ('histogram','kde'):
plotVxlsPerAtm(self.pdbname,self.filesOut,self.vxlsPerAtom,plotType)
