Changeset 1368


Ignore:
Timestamp:
May 29, 2014 11:19:07 AM (8 years ago)
Author:
vondreele
Message:

implement sqrt(intensity) powder plots

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1365 r1368  
    19401940        self.qPlot = False
    19411941        self.sqPlot = False
     1942        self.SqrtPlot = False
    19421943        self.ErrorBars = False
    19431944        self.Contour = False
  • trunk/GSASIIplot.py

    r1366 r1368  
    462462                    G2frame.Offset[0] = 0
    463463                newPlot = True
     464        elif event.key == 's' and 'PWDR' in plottype:
     465            if G2frame.Contour:
     466                choice = [m for m in mpl.cm.datad.keys() if not m.endswith("_r")]
     467                choice.sort()
     468                dlg = wx.SingleChoiceDialog(G2frame,'Select','Color scheme',choice)
     469                if dlg.ShowModal() == wx.ID_OK:
     470                    sel = dlg.GetSelection()
     471                    G2frame.ContourColor = choice[sel]
     472                else:
     473                    G2frame.ContourColor = 'Paired'
     474                dlg.Destroy()
     475            elif G2frame.SinglePlot:
     476                G2frame.SqrtPlot = not G2frame.SqrtPlot
     477                if G2frame.SqrtPlot:
     478                    G2frame.delOffset = .002
     479                    G2frame.refOffset = -10.0
     480                    G2frame.refDelt = .001
     481                else:
     482                    G2frame.delOffset = .02
     483                    G2frame.refOffset = -100.0
     484                    G2frame.refDelt = .01
     485            newPlot = True
    464486        elif event.key == 'u' and (G2frame.Contour or not G2frame.SinglePlot):
    465487            if G2frame.Contour:
     
    495517                G2frame.sqPlot = not G2frame.sqPlot       
    496518        elif event.key == 'm':
    497             if G2frame.Contour:
    498                 choice = [m for m in mpl.cm.datad.keys() if not m.endswith("_r")]
    499                 choice.sort()
    500                 dlg = wx.SingleChoiceDialog(G2frame,'Select','Color scheme',choice)
    501                 if dlg.ShowModal() == wx.ID_OK:
    502                     sel = dlg.GetSelection()
    503                     G2frame.ContourColor = choice[sel]
    504                 else:
    505                     G2frame.ContourColor = 'Paired'
    506                 dlg.Destroy()
    507             else:
    508                 G2frame.SinglePlot = not G2frame.SinglePlot               
     519            G2frame.SqrtPlot = False
     520            G2frame.SinglePlot = not G2frame.SinglePlot               
    509521            newPlot = True
    510522        elif event.key == '+':
     
    790802                if G2frame.SinglePlot:
    791803                    Page.Choice = (' key press',
    792                         'b: toggle subtract background','n: log(I) on','c: contour on',
     804                        'b: toggle subtract background','n: log(I) on','s: toggle sqrt plot','c: contour on',
    793805                        'q: toggle q plot','m: toggle multidata plot','w: toggle divide by sig','+: no selection')
    794806                else:
     
    879891        if 'C' in ParmList[0]['Type'][0]:
    880892            if 'PWDR' in plottype:
    881                 Plot.set_ylabel('$Intensity$',fontsize=16)
     893                if G2frame.SqrtPlot:
     894                    Plot.set_ylabel(r'$\sqrt{Intensity}$',fontsize=16)
     895                else:
     896                    Plot.set_ylabel('$Intensity$',fontsize=16)
    882897            elif 'SASD' in plottype:
    883898                if G2frame.sqPlot:
     
    921936            lenX = len(X)
    922937        if 'PWDR' in plottype:
    923             Y = xye[1]+offset*N
     938            if G2frame.SqrtPlot:
     939                Y = np.where(xye[1]>=0.,np.sqrt(xye[1]),-np.sqrt(-xye[1]))
     940            else:
     941                Y = xye[1]+offset*N
    924942        elif 'SASD' in plottype:
    925943            B = xye[5]
     
    957975            DifLine = ['']
    958976            if ifpicked:
    959                 Z = xye[3]+offset*N
     977                if G2frame.SqrtPlot:
     978                    Z = np.where(xye[3]>=0.,np.sqrt(xye[3]),-np.sqrt(-xye[3]))
     979                else:
     980                    Z = xye[3]+offset*N
    960981                if 'PWDR' in plottype:
    961                     W = xye[4]+offset*N
    962                     D = xye[5]-Ymax*G2frame.delOffset  #powder background
     982                    if G2frame.SqrtPlot:
     983                        W = np.where(xye[4]>=0.,np.sqrt(xye[4]),-np.sqrt(-xye[4]))
     984                        D = np.where(xye[5],(Y-Z),0.)-Ymax*G2frame.delOffset
     985                    else:
     986                        W = xye[4]+offset*N
     987                        D = xye[5]-Ymax*G2frame.delOffset  #powder background
    963988                elif 'SASD' in plottype:
    964989                    if G2frame.sqPlot:
Note: See TracChangeset for help on using the changeset viewer.