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)