Changeset 4653 for trunk/GSASIIplot.py
- Timestamp:
- Nov 8, 2020 12:47:28 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIplot.py
r4652 r4653 1059 1059 if Type in ['|DFsq|/sig','|DFsq|>sig','|DFsq|>3sig']: 1060 1060 if A > 0.0: 1061 Plot.add_artist(Circle(xy,radius=A,ec='g',fc='w',pick er=1.,gid=hid))1061 Plot.add_artist(Circle(xy,radius=A,ec='g',fc='w',pickradius=1.,gid=hid)) 1062 1062 else: 1063 Plot.add_artist(Circle(xy,radius=-A,ec='r',fc='w',pick er=1.,gid=hid))1063 Plot.add_artist(Circle(xy,radius=-A,ec='r',fc='w',pickradius=1.,gid=hid)) 1064 1064 else: 1065 1065 if A > 0.0 and A > B: 1066 1066 Plot.add_artist(Circle(xy,radius=A,ec='g',fc='w')) 1067 1067 if B: 1068 Plot.add_artist(Circle(xy,radius=B,ec='b',fc='w',pick er=1.,gid=hid))1068 Plot.add_artist(Circle(xy,radius=B,ec='b',fc='w',pickradius=1.,gid=hid)) 1069 1069 if A < B: 1070 1070 Plot.add_artist(Circle(xy,radius=A,ec='g',fc='w')) … … 2372 2372 backDict['FixedPoints'].append(xy) 2373 2373 Plot = Page.figure.gca() 2374 Plot.plot(event.xdata,event.ydata,'rD',clip_on=Clip_on,pick er=3.)2374 Plot.plot(event.xdata,event.ydata,'rD',clip_on=Clip_on,pickradius=3.) 2375 2375 Page.canvas.draw() 2376 2376 return … … 3062 3062 if x > xmax: 3063 3063 continue 3064 magMarkers.append(Plot.axvline(x,color='0.5',dashes=(1,1),pick er=2.,label='_magline'))3064 magMarkers.append(Plot.axvline(x,color='0.5',dashes=(1,1),pickradius=2.,label='_magline')) 3065 3065 lbl = Plot.annotate("x{}".format(ml), xy=(x, tpos), xycoords=("data", "axes fraction"), 3066 3066 verticalalignment='bottom',horizontalalignment=halign,label='_maglbl') … … 3103 3103 elif Page.plotStyle['dPlot'] and 'PWDR' in plottype and not ifLimits: 3104 3104 lims = G2lat.Pos2dsp(Parms,lims) 3105 Lines.append(Plot.axvline(lims[0],color='g',dashes=(5,5),pick er=3.))3106 Lines.append(Plot.axvline(lims[1],color='r',dashes=(5,5),pick er=3.))3105 Lines.append(Plot.axvline(lims[0],color='g',dashes=(5,5),pickradius=3.)) 3106 Lines.append(Plot.axvline(lims[1],color='r',dashes=(5,5),pickradius=3.)) 3107 3107 for i,item in enumerate(limits[2:]): 3108 Lines.append(Plot.axvline(item[0],color='m',dashes=(5,5),pick er=3.))3109 Lines.append(Plot.axvline(item[1],color='m',dashes=(5,5),pick er=3.))3108 Lines.append(Plot.axvline(item[0],color='m',dashes=(5,5),pickradius=3.)) 3109 Lines.append(Plot.axvline(item[1],color='m',dashes=(5,5),pickradius=3.)) 3110 3110 exclLines += [2*i+2,2*i+3] 3111 3111 if G2frame.Contour: … … 3199 3199 else: 3200 3200 DZ = (xye[1]-xye[3])*np.sqrt(wtFactor*xye[2]) 3201 DifLine = Plot1.plot(X,DZ,colors[3],pick er=1.,label=incCptn('diff')) #(Io-Ic)/sig(Io)3201 DifLine = Plot1.plot(X,DZ,colors[3],pickradius=1.,label=incCptn('diff')) #(Io-Ic)/sig(Io) 3202 3202 Plot1.axhline(0.,color='k') 3203 3203 … … 3205 3205 if 'PWDR' in plottype: 3206 3206 Plot.set_yscale("log",nonposy='mask') 3207 Plot.plot(X,Y,marker=pP,color=colors[0],pick er=3.,clip_on=Clip_on,label=incCptn('obs'))3207 Plot.plot(X,Y,marker=pP,color=colors[0],pickradius=3.,clip_on=Clip_on,label=incCptn('obs')) 3208 3208 if G2frame.SinglePlot or G2frame.plusPlot: 3209 3209 Plot.plot(X,Z,colors[1],picker=False,label=incCptn('calc')) … … 3216 3216 if Page.plotStyle['sqPlot']: 3217 3217 Plot.errorbar(X,YB,yerr=X**4*Sample['Scale'][0]*np.sqrt(1./(Pattern[0]['wtFactor']*xye[2])), 3218 ecolor=colors[0],pick er=3.,clip_on=Clip_on)3218 ecolor=colors[0],pickradius=3.,clip_on=Clip_on) 3219 3219 else: 3220 3220 Plot.errorbar(X,YB,yerr=Sample['Scale'][0]*np.sqrt(1./(Pattern[0]['wtFactor']*xye[2])), 3221 ecolor=colors[0],pick er=3.,clip_on=Clip_on,label=incCptn('obs'))3221 ecolor=colors[0],pickradius=3.,clip_on=Clip_on,label=incCptn('obs')) 3222 3222 else: 3223 Plot.plot(X,YB,marker=pP,color=colors[0],pick er=3.,clip_on=Clip_on,label=incCptn('obs'))3223 Plot.plot(X,YB,marker=pP,color=colors[0],pickradius=3.,clip_on=Clip_on,label=incCptn('obs')) 3224 3224 Plot.plot(X,W,colors[1],picker=False,label=incCptn('bkg')) #const. background 3225 3225 Plot.plot(X,ZB,colors[2],picker=False,label=incCptn('calc')) … … 3231 3231 CalcLine = Plot.plot(X,Z-W,colors[1],picker=False,label=incCptn('calc-bkg')) #Ic-Ib 3232 3232 else: 3233 Plot.plot(X,YB,color=colors[0],marker=pP,pick er=3.,clip_on=Clip_on,label=incCptn('obs'))3233 Plot.plot(X,YB,color=colors[0],marker=pP,pickradius=3.,clip_on=Clip_on,label=incCptn('obs')) 3234 3234 Plot.plot(X,ZB,colors[2],picker=False,label=incCptn('calc')) 3235 3235 else: 3236 3236 if 'PWDR' in plottype: 3237 ObsLine = Plot.plot(Xum,Y,color=colors[0],marker=pP,pick er=3.,clip_on=Clip_on,label=incCptn('obs')) #Io3237 ObsLine = Plot.plot(Xum,Y,color=colors[0],marker=pP,pickradius=3.,clip_on=Clip_on,label=incCptn('obs')) #Io 3238 3238 CalcLine = Plot.plot(X,Z,colors[1],picker=False,label=incCptn('calc')) #Ic 3239 3239 else: 3240 Plot.plot(X,YB,color=colors[0],marler=pP,pick er=3.,clip_on=Clip_on,label=incCptn('obs'))3240 Plot.plot(X,YB,color=colors[0],marler=pP,pickradius=3.,clip_on=Clip_on,label=incCptn('obs')) 3241 3241 Plot.plot(X,ZB,colors[2],picker=False,label=incCptn('calc')) 3242 3242 if 'PWDR' in plottype and (G2frame.SinglePlot and G2frame.plusPlot): 3243 3243 BackLine = Plot.plot(X,W,colors[2],picker=False,label=incCptn('bkg')) #Ib 3244 3244 if not G2frame.Weight and np.any(Z): 3245 DifLine = Plot.plot(X,D,colors[3],pick er=1.,label=incCptn('diff')) #Io-Ic3245 DifLine = Plot.plot(X,D,colors[3],pickradius=1.,label=incCptn('diff')) #Io-Ic 3246 3246 Plot.axhline(0.,color='k',label='_zero') 3247 3247 Page.SetToolTipString('') … … 3261 3261 Ni = N 3262 3262 if Page.plotStyle['qPlot']: 3263 Lines.append(Plot.axvline(2.*np.pi/G2lat.Pos2dsp(Parms,item[0]),color='b',pick er=2.))3263 Lines.append(Plot.axvline(2.*np.pi/G2lat.Pos2dsp(Parms,item[0]),color='b',pickradius=2.)) 3264 3264 elif Page.plotStyle['dPlot']: 3265 Lines.append(Plot.axvline(G2lat.Pos2dsp(Parms,item[0]),color='b',pick er=2.))3265 Lines.append(Plot.axvline(G2lat.Pos2dsp(Parms,item[0]),color='b',pickradius=2.)) 3266 3266 else: 3267 Lines.append(Plot.axvline(item[0],color='b',pick er=2.))3267 Lines.append(Plot.axvline(item[0],color='b',pickradius=2.)) 3268 3268 if Ni == N+1: 3269 3269 Lines[-1].set_lw(Lines[-1].get_lw()+1) … … 3343 3343 plsym = Page.phaseColors[phase]+'|' # yellow should never happen! 3344 3344 if Page.plotStyle['qPlot']: 3345 Page.tickDict[phase],j = Plot.plot(2*np.pi/peak.T[0],pos,plsym,mew=w,ms=l,pick er=3.,label=phase)3345 Page.tickDict[phase],j = Plot.plot(2*np.pi/peak.T[0],pos,plsym,mew=w,ms=l,pickradius=3.,label=phase) 3346 3346 elif Page.plotStyle['dPlot']: 3347 Page.tickDict[phase],j = Plot.plot(peak.T[0],pos,plsym,mew=w,ms=l,pick er=3.,label=phase)3347 Page.tickDict[phase],j = Plot.plot(peak.T[0],pos,plsym,mew=w,ms=l,pickradius=3.,label=phase) 3348 3348 else: 3349 Page.tickDict[phase],j = Plot.plot(peak.T[1],pos,plsym,mew=w,ms=l,pick er=3.,label=phase)3349 Page.tickDict[phase],j = Plot.plot(peak.T[1],pos,plsym,mew=w,ms=l,pickradius=3.,label=phase) 3350 3350 handles,legends = Plot.get_legend_handles_labels() #got double entries in the phase legends for some reason 3351 3351 if handles: … … 3401 3401 else: 3402 3402 break 3403 Plot.plot(x,y,'rD',clip_on=Clip_on,pick er=3.)3403 Plot.plot(x,y,'rD',clip_on=Clip_on,pickradius=3.) 3404 3404 if not newPlot: 3405 3405 # this restores previous plot limits (but I'm not sure why there are two .push_current calls) … … 4782 4782 X = XYlist[0].T[0] 4783 4783 Y = XYlist[0].T[1] 4784 Plot.plot(X,Y,color='b',pick er=3)4784 Plot.plot(X,Y,color='b',pickradius=3) 4785 4785 if 'calc' in Peaks and len(Peaks['calc']): 4786 4786 XC,YC= Peaks['calc'] 4787 4787 Plot.plot(XC,YC,color='g') 4788 G2frame.Lines.append(Plot.axvline(peaks['Limits'][0],color='g',dashes=(5,5),pick er=2.))4789 G2frame.Lines.append(Plot.axvline(peaks['Limits'][1],color='r',dashes=(5,5),pick er=2.))4788 G2frame.Lines.append(Plot.axvline(peaks['Limits'][0],color='g',dashes=(5,5),pickradius=2.)) 4789 G2frame.Lines.append(Plot.axvline(peaks['Limits'][1],color='r',dashes=(5,5),pickradius=2.)) 4790 4790 for peak in Peaks['Peaks']: 4791 G2frame.Lines.append(Plot.axvline(peak[0],color='r',pick er=2.))4791 G2frame.Lines.append(Plot.axvline(peak[0],color='r',pickradius=2.)) 4792 4792 Xb = [0.,peaks['Limits'][1]] 4793 4793 Yb = [0.,Xb[1]*peaks['Background'][1][1]] … … 4872 4872 Plot.errorbar(X,Y,ecolor='k',yerr=E) 4873 4873 if ixy: 4874 Plot.plot(X,Y,'kx',pick er=3)4874 Plot.plot(X,Y,'kx',pickradius=3) 4875 4875 else: 4876 4876 Plot.plot(X,Y,'kx',label='peak') … … 5289 5289 colors = list(np.where(np.array(Probs1)>thresh[0][1],'r','b')) 5290 5290 resNums = np.arange(len(resNames)) 5291 Plot1.bar(resNums,Probs1,color=colors,linewidth=0,pick er=1)5291 Plot1.bar(resNums,Probs1,color=colors,linewidth=0,pickradius=1) 5292 5292 if thresh is not None: 5293 5293 for item in thresh[0]: … … 5297 5297 Plot2.set_xlabel(r'Residue',fontsize=14) 5298 5298 colors = list(np.where(np.array(Probs2)>thresh[1][1],'r','b')) 5299 Plot2.bar(resNums,Probs2,color=colors,linewidth=0,pick er=1)5299 Plot2.bar(resNums,Probs2,color=colors,linewidth=0,pickradius=1) 5300 5300 if thresh is not None: 5301 5301 for item in thresh[1]: … … 6043 6043 acolor = mpl.cm.get_cmap(G2frame.ContourColor) 6044 6044 Img = Plot.imshow(Z.T,aspect='equal',cmap=acolor,extent=[-1,1,-1,1]) 6045 Plot.plot(y,x,'+',pick er=3)6045 Plot.plot(y,x,'+',pickradius=3) 6046 6046 Page.figure.colorbar(Img) 6047 6047 Plot.axis('off') … … 6238 6238 Ydet = list(npsind(Pdet)*Rdet) 6239 6239 for i,[x,y] in enumerate(zip(Xdet,Ydet)): 6240 Plot.plot(x,-y,'k+',pick er=5,gid=textureData['det Angles'][i][0])6240 Plot.plot(x,-y,'k+',pickradius=5,gid=textureData['det Angles'][i][0]) 6241 6241 h,k,l = SHData['PFhkl'] 6242 6242 Plot.axis('off') … … 6520 6520 if len(Angles): 6521 6521 Eval = np.array([-G2mth.calcTorsionEnergy(x,Coeff)[1] for x in Angles]) 6522 Plot.plot(Angles,Eval,'ro',pick er=5)6522 Plot.plot(Angles,Eval,'ro',pickradius=5) 6523 6523 Plot.set_xlim((0.,360.)) 6524 6524 Plot.set_title('Torsion angles for '+TorName+' in '+phaseName) … … 6601 6601 PhiPsi = np.where(PhiPsi>180.,PhiPsi-360.,PhiPsi) 6602 6602 Phi,Psi = PhiPsi.T 6603 Plot.plot(Phi,Psi,'ro',pick er=5)6603 Plot.plot(Phi,Psi,'ro',pickradius=5) 6604 6604 Plot.set_xlim((-180.,180.)) 6605 6605 Plot.set_ylim((-180.,180.)) … … 6613 6613 if len(PhiPsi): 6614 6614 Phi,Psi = PhiPsi.T 6615 Plot.plot(Phi,Psi,'ro',pick er=5)6615 Plot.plot(Phi,Psi,'ro',pickradius=5) 6616 6616 Plot.set_xlim((0.,360.)) 6617 6617 Plot.set_ylim((0.,360.)) … … 7028 7028 spot = [event.xdata,event.ydata,G2frame.spotSize] 7029 7029 Masks['Points'].append(spot) 7030 artist = Circle(spot[:2],radius=spot[2]/2,fc='none',ec='r',pick er=3)7030 artist = Circle(spot[:2],radius=spot[2]/2,fc='none',ec='r',pickradius=3) 7031 7031 Page.figure.gca().add_artist(artist) 7032 7032 artist.itemNumber = len(Masks['Points'])-1 … … 7460 7460 spot = [Xpos,Ypos,sig] 7461 7461 Masks['Points'].append(spot) 7462 artist = Circle((Xpos,Ypos),radius=spot[2]/2,fc='none',ec='r',pick er=3)7462 artist = Circle((Xpos,Ypos),radius=spot[2]/2,fc='none',ec='r',pickradius=3) 7463 7463 Page.figure.gca().add_artist(artist) 7464 7464 artist.itemNumber = len(Masks['Points'])-1 … … 7757 7757 xyI = np.array(xyI) 7758 7758 arcxI,arcyI = xyI.T 7759 Plot.plot(arcxI,arcyI,pick er=3,label='Itth')7759 Plot.plot(arcxI,arcyI,pickradius=3,label='Itth') 7760 7760 if ellO: 7761 7761 xyO = [] … … 7768 7768 xyO = np.array(xyO) 7769 7769 arcxO,arcyO = xyO.T 7770 Plot.plot(arcxO,arcyO,pick er=3,label='Otth')7770 Plot.plot(arcxO,arcyO,pickradius=3,label='Otth') 7771 7771 if ellO and ellI and len(arcxO): 7772 Plot.plot([arcxI[0],arcxO[0]],[arcyI[0],arcyO[0]],pick er=3,label='Lazm')7773 Plot.plot([arcxI[-1],arcxO[-1]],[arcyI[-1],arcyO[-1]],pick er=3,label='Uazm')7772 Plot.plot([arcxI[0],arcxO[0]],[arcyI[0],arcyO[0]],pickradius=3,label='Lazm') 7773 Plot.plot([arcxI[-1],arcxO[-1]],[arcyI[-1],arcyO[-1]],pickradius=3,label='Uazm') 7774 7774 for i in range(Nazm): 7775 7775 cake = LRAzim[0]+i*delAzm-AzmthOff … … 7787 7787 dspO = wave/(2.0*sind(IOtth[1]/2.0)) 7788 7788 xyO = G2img.GetDetectorXY(dspO,azm,Data) 7789 Plot.plot([xyI[0],xyO[0]],[xyI[1],xyO[1]],pick er=3,label='linescan')7789 Plot.plot([xyI[0],xyO[0]],[xyI[1],xyO[1]],pickradius=3,label='linescan') 7790 7790 7791 7791 if G2frame.PickId and G2frame.GPXtree.GetItemText(G2frame.PickId) in ['Image Controls',]: 7792 7792 for xring,yring in Data['ring']: 7793 Plot.plot(xring,yring,'r+',pick er=3)7793 Plot.plot(xring,yring,'r+',pickradius=3) 7794 7794 if Data['setRings']: 7795 7795 N = 0 … … 7815 7815 if len(spot): 7816 7816 x,y,d = spot 7817 artist = Circle((x,y),radius=d/2,fc='none',ec='r',pick er=3)7817 artist = Circle((x,y),radius=d/2,fc='none',ec='r',pickradius=3) 7818 7818 Plot.add_artist(artist) 7819 7819 artist.itemNumber = i … … 7826 7826 wave = Data['wavelength'] 7827 7827 (x1,y1),(x2,y2) = ComputeArc(tth-thick/2.,tth+thick/2.,wave) 7828 artistO, = Plot.plot(x1,y1,'r',pick er=3)7828 artistO, = Plot.plot(x1,y1,'r',pickradius=3) 7829 7829 artistO.itemNumber = iring 7830 7830 artistO.itemType = 'RingOuter' 7831 artistI, = Plot.plot(x2,y2,'r',pick er=3)7831 artistI, = Plot.plot(x2,y2,'r',pickradius=3) 7832 7832 artistI.itemNumber = iring 7833 7833 artistI.itemType = 'RingInner' … … 7842 7842 (x1,y1),(x2,y2) = ComputeArc(tth-thick/2.,tth+thick/2.,wave,azm[0],azm[1]) 7843 7843 arcList = [] 7844 arcList.append(Plot.plot(x2,y2,'r',pick er=3)[0]) # 'inner'7844 arcList.append(Plot.plot(x2,y2,'r',pickradius=3)[0]) # 'inner' 7845 7845 arcList[-1].itemNumber = iarc 7846 7846 arcList[-1].itemType = 'ArcInner' 7847 arcList.append(Plot.plot(x1,y1,'r',pick er=3)[0]) # 'outer'7847 arcList.append(Plot.plot(x1,y1,'r',pickradius=3)[0]) # 'outer' 7848 7848 arcList[-1].itemNumber = iarc 7849 7849 arcList[-1].itemType = 'ArcOuter' 7850 arcList.append(Plot.plot([x1[0],x2[0]],[y1[0],y2[0]],'r',pick er=3)[0]) # 'lower'7850 arcList.append(Plot.plot([x1[0],x2[0]],[y1[0],y2[0]],'r',pickradius=3)[0]) # 'lower' 7851 7851 arcList[-1].itemNumber = iarc 7852 7852 arcList[-1].itemType = 'ArcLower' 7853 arcList.append(Plot.plot([x1[-1],x2[-1]],[y1[-1],y2[-1]],'r',pick er=3)[0]) # 'upper'7853 arcList.append(Plot.plot([x1[-1],x2[-1]],[y1[-1],y2[-1]],'r',pickradius=3)[0]) # 'upper' 7854 7854 arcList[-1].itemNumber = iarc 7855 7855 arcList[-1].itemType = 'ArcUpper' … … 7865 7865 G2frame.polyList.append(Plot.plot(xl,yl,'r')[0]) # line 7866 7866 for i,(x,y) in enumerate(zip(xl[:-1],yl[:-1])): 7867 artist = Plot.plot(x,y,'r+',pick er=10)[0] # point (plus sign)7867 artist = Plot.plot(x,y,'r+',pickradius=10)[0] # point (plus sign) 7868 7868 artist.itemNumber = ipoly 7869 7869 artist.itemType = 'Polygon' … … 7879 7879 G2frame.frameArtist = Plot.plot(xl,yl,'g')[0] 7880 7880 for i,(x,y) in enumerate(zip(xl[:-1],yl[:-1])): 7881 artist = Plot.plot(x,y,'g+',pick er=10)[0] # point (plus sign)7881 artist = Plot.plot(x,y,'g+',pickradius=10)[0] # point (plus sign) 7882 7882 artist.itemType = 'Frame' 7883 7883 artist.pointNumber = i
Note: See TracChangeset
for help on using the changeset viewer.