Changeset 5282


Ignore:
Timestamp:
May 18, 2022 2:24:07 PM (5 months ago)
Author:
vondreele
Message:

fixes to partial plots

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIdataGUI.py

    r5276 r5282  
    56545654                done = True
    56555655            if valList:
    5656                 phPartialFile = filename.replace('_part_N',
    5657                                                  '_part_'+str(hId)+'.csv')
     5656                phPartialFile = filename.replace('_part_N','_part_'+str(hId)+'.csv')
    56585657                fp1 = open(phPartialFile,'w')
    56595658                fp1.write(', '.join(lblList))
  • trunk/GSASIIplot.py

    r5276 r5282  
    299299    fil.write(line+'\n')
    300300
    301 
    302301class _tabPlotWin(wx.Panel):   
    303302    'Creates a basic tabbed plot window for GSAS-II graphics'
     
    309308        self.plotRequiresRedraw = True # delete plot if not updated
    310309        wx.Panel.__init__(self,parent,id=id,**kwargs)
    311        
    312310           
    313311class G2PlotMpl(_tabPlotWin):   
     
    321319        self.toolbar = GSASIItoolbar(self.canvas,publish=publish)
    322320        self.toolbar.Realize()
    323         self.plotStyle = {'qPlot':False,'dPlot':False,'sqrtPlot':False,'sqPlot':False,'logPlot':False,'exclude':False}
     321        self.plotStyle = {'qPlot':False,'dPlot':False,'sqrtPlot':False,'sqPlot':False,
     322            'logPlot':False,'exclude':False,'partials':True}
    324323       
    325324        sizer=wx.BoxSizer(wx.VERTICAL)
     
    368367        else:
    369368            self.canvas.SetToolTipString(text)
    370            
    371369       
    372370class G2Plot3D(_tabPlotWin):
     
    21792177                Page.plotStyle['qPlot'] = False
    21802178                Page.plotStyle['dPlot'] = False
     2179        elif (event.key == 'p' and 'PWDR' in plottype and G2frame.SinglePlot):
     2180            Page.plotStyle['partials'] = not Page.plotStyle['partials']
    21812181        elif (event.key == 'e' and 'PWDR' in plottype and G2frame.SinglePlot and ifLimits
    21822182                  and not G2frame.Contour):
     
    32393239                     'o: reset offset','f: select data',
    32403240                     '/: normalize']
     3241            else:
     3242                Page.Choice = Page.Choice+ ['p: toggle partials (if available)',]
    32413243        elif plottype in ['SASD','REFD']:
    32423244            Page.Choice = [' key press',
     
    38943896    # plot the partials. TODO: get partials to show up in publication plot
    38953897    plotOpt['lineList']  = ['obs','calc','bkg','zero','diff']
    3896     if 'PWDR' in plottype and G2frame.SinglePlot:
    3897         x, yb, ypList = G2frame.LoadPartial(data[0]['hId'])
    3898         if x is not None:
     3898    if 'PWDR' in plottype and G2frame.SinglePlot and Page.plotStyle['partials'] and 'hId' in data[0]:
     3899        x, yb, ypList = G2frame.LoadPartial(data[0]['hId'])           
     3900        if x is not None and len(ypList) > 1:
    38993901            if Page.plotStyle['qPlot']:
    39003902                x = 2.*np.pi/G2lat.Pos2dsp(Parms,x)
     
    39033905            olderr = np.seterr(invalid='ignore') #get around sqrt(-ve) error
    39043906            for ph in ypList:
    3905                 #ypList[ph] += yb    # looks better w/o adding background
     3907                if G2frame.SubBack:
     3908                    y = ypList[ph]
     3909                else:
     3910                    y = ypList[ph]+yb
    39063911                if Page.plotStyle['sqrtPlot']:
    3907                     y = np.where(ypList[ph]>=0.,np.sqrt(ypList[ph]),-np.sqrt(-ypList[ph]))
    3908                 else:
    3909                     y = ypList[ph]
    3910                 Plot.plot(x,y,Page.phaseColors.get(ph,'k'),picker=False,
    3911                           label=ph,linewidth=2.5)
     3912                    y = np.where(y>=0.,np.sqrt(y),-np.sqrt(-y))
     3913                Plot.plot(x,y,Page.phaseColors.get(ph,'k'),picker=False,label=ph,linewidth=1.5,dashes=(5,5))
    39123914                plotOpt['lineList'].append(ph)   # needed?
    39133915    if not newPlot:
Note: See TracChangeset for help on using the changeset viewer.