Changeset 3045 for trunk/GSASIIplot.py


Ignore:
Timestamp:
Sep 4, 2017 12:46:53 PM (4 years ago)
Author:
vondreele
Message:

plot both protein vlidators (#2 still in erroe)
remove config value for choice - not needed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r3042 r3045  
    31513151##### PlotHist
    31523152################################################################################
    3153 def PlotAAProb(G2frame,resNames,Probs,Title='',thresh=None,newPlot=True):
     3153def PlotAAProb(G2frame,resNames,Probs1,Probs2,Title='',thresh=None,newPlot=True):
    31543154
    31553155
     
    31703170        Plot.clear()
    31713171        Plot.set_title(Title)
    3172         Plot.set_xlabel(r'Error score',fontsize=14)
    3173         Plot.set_ylabel(r'Residue',fontsize=14)
    3174         try:
    3175             Plot.barh(np.arange(len(resNames)),Probs,tick_label=resNames,)
    3176         except AttributeError:  #tick_label not in mpl 1.4
    3177             Plot.barh(np.arange(len(resNames)),Probs)
     3172        Plot.set_axis_off()
     3173        Plot1 = Page.figure.add_subplot(211)
     3174        Plot1.set_xlabel(r'Error score 1',fontsize=14)
     3175        Plot1.set_ylabel(r'Residue',fontsize=14)
     3176        colors = list(np.where(np.array(Probs1)>thresh[0][1],'r','b'))
     3177        Plot1.barh(np.arange(len(resNames)),Probs1,color=colors,linewidth=0)
    31783178        if thresh is not None:
    3179             for item in thresh:
    3180                 Plot.axvline(item,dashes=(5,5),picker=False)
    3181         if not newPlot:
    3182             Page.toolbar.push_current()
    3183             Plot.set_xlim(xylim[0])
    3184             Plot.set_ylim(xylim[1])
    3185             xylim = []
    3186             Page.toolbar.push_current()
    3187             Page.toolbar.draw()
    3188             Page.canvas.draw()
    3189         else:
    3190             Page.canvas.draw()
     3179            for item in thresh[0]:
     3180                Plot1.axvline(item,dashes=(5,5),picker=False)
     3181        Plot2 = Page.figure.add_subplot(212)
     3182        Plot2.set_xlabel(r'Error score 2',fontsize=14)
     3183        Plot2.set_ylabel(r'Residue',fontsize=14)       
     3184        colors = list(np.where(np.array(Probs2)>thresh[1][1],'r','b'))
     3185        Plot2.barh(np.arange(len(resNames)),Probs2,color=colors,linewidth=0)
     3186        if thresh is not None:
     3187            for item in thresh[1]:
     3188                Plot2.axvline(item,dashes=(5,5),picker=False)
     3189        Page.canvas.draw()
    31913190   
    31923191    new,plotNum,Page,Plot,lim = G2frame.G2plotNB.FindPlotTab(Title,'mpl')
Note: See TracChangeset for help on using the changeset viewer.