Changeset 1806


Ignore:
Timestamp:
Apr 24, 2015 10:03:39 AM (7 years ago)
Author:
vondreele
Message:

add 'L' option to plot/not plot lines between points in seqPlots

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1805 r1806  
    21482148        self.SubBack = False
    21492149        self.seqReverse = False
     2150        self.seqLines = True #draw lines between points
    21502151        self.plotView = 0
    21512152        self.Image = 0
  • trunk/GSASIImath.py

    r1800 r1806  
    16561656        for hist in Gangls.keys():
    16571657            Refs = refData[hist]
     1658            wt = 1.0    #   np.sqrt(np.abs(Refs[:,5]))?
    16581659            sumObs += np.sum(Refs[:,5])
    16591660            Refs[:,6] = 1.
     
    16681669                    Lnorm = G2lat.Lnorm(L)
    16691670                    Refs[:,6] += parmDict[item]*Lnorm*Kcl*Ksl
    1670             mat = Refs[:,5]-Refs[:,6]
     1671            mat = wt*(Refs[:,5]-Refs[:,6])
    16711672            Mat = np.concatenate((Mat,mat))
    16721673        sumD = np.sum(np.abs(Mat))
     
    16831684            Refs = refData[hist]
    16841685            H = Refs[:,:3]
     1686            wt = 1.0    #   np.sqrt(np.abs(Refs[:,5]))?
    16851687            phi,beta = G2lat.CrsAng(H,cell,SGData)
    16861688            psi,gam,dPdA,dGdA = G2lat.SamAng(Refs[:,3]/2.,Gangls[hist],Sangls,False) #assume not Bragg-Brentano!
     
    16911693                    Ksl,dKdp,dKdg = G2lat.GetKsl(L,M,shModel,psi,gam)
    16921694                    Lnorm = G2lat.Lnorm(L)
    1693                     mat[j] = -Lnorm*Kcl*Ksl
     1695                    mat[j] = -wt*Lnorm*Kcl*Ksl
    16941696                    for k,itema in enumerate(['Sample omega','Sample chi','Sample phi']):
    16951697                        try:
    16961698                            l = varyList.index(itema)
    1697                             mat[l] -= parmDict[item]*Lnorm*Kcl*(dKdp*dPdA[k]+dKdg*dGdA[k])
     1699                            mat[l] -= parmDict[item]*wt*Lnorm*Kcl*(dKdp*dPdA[k]+dKdg*dGdA[k])
    16981700                        except ValueError:
    16991701                            pass
  • trunk/GSASIIphsGUI.py

    r1804 r1806  
    35603560            wx.CallAfter(G2plt.PlotTexture,G2frame,data)
    35613561           
     3562        def OnPlaneSel(event):
     3563            Obj = event.GetEventObject()
     3564            textureData['PFproj'] = Obj.GetValue()
     3565            wx.CallAfter(G2plt.PlotTexture,G2frame,data)
     3566           
     3567        def OnRevPF(event):
     3568            Obj = event.GetEventObject()
     3569            textureData['PFrev'] = Obj.GetValue()
     3570            wx.CallAfter(G2plt.PlotTexture,G2frame,data)
     3571           
    35623572        def OnAngRef(event):
    35633573            Obj = event.GetEventObject()
     
    36953705                style=wx.CB_READONLY|wx.CB_DROPDOWN)
    36963706            colorSel.Bind(wx.EVT_COMBOBOX,OnColorSel)
    3697             PTSizer.Add(colorSel,0,WACV)       
     3707            PTSizer.Add(colorSel,0,WACV)
    36983708        mainSizer.Add(PTSizer,0,WACV)
    36993709        mainSizer.Add((0,5),0)
  • trunk/GSASIIplot.py

    r1804 r1806  
    28232823                Plot.set_zlabel(r'Z, MRD')
    28242824        else:
     2825            PFproj = textureData.get('PFproj','XY')
     2826            PRrev = textureData.get('PFrev',False)
    28252827            X,Y = np.meshgrid(np.linspace(1.,-1.,npts),np.linspace(-1.,1.,npts))
    28262828            R,P = np.sqrt(X**2+Y**2).flatten(),npatan2d(X,Y).flatten()
     
    32513253            dlg.Destroy()
    32523254            Draw()
     3255        elif event.key == 'l':
     3256            G2frame.seqLines = not G2frame.seqLines
     3257            wx.CallAfter(Draw)
    32533258           
    32543259    def Draw():
    32553260        global Title,xLabel,yLabel
    32563261        Page.SetFocus()
    3257         G2frame.G2plotNB.status.SetStatusText('press s to select X axis, t to change titles',1)
     3262        G2frame.G2plotNB.status.SetStatusText(  \
     3263            'press L to toggle lines, S to select X axis, T to change titles (reselect column to show?)',1)
    32583264        Plot.clear()
    32593265        if G2frame.seqXaxis is not None:   
     
    32773283                    Ynew = Ynew[::-1]
    32783284                    Ysnew = Ysnew[::-1]
    3279                 Plot.errorbar(Xnew,Ynew,yerr=Ysnew,label=name)
     3285                if G2frame.seqLines:
     3286                    Plot.errorbar(Xnew,Ynew,yerr=Ysnew,label=name)
     3287                else:
     3288                    Plot.errorbar(Xnew,Ynew,yerr=Ysnew,label=name,linestyle='None',marker='x')
    32803289            else:
    32813290                if G2frame.seqReverse and not G2frame.seqXaxis:
     
    32893298           
    32903299        Plot.legend(loc='best')
    3291         print Title,xLabel,yLabel
    32923300        if Title:
    32933301            Plot.set_title(Title)
     
    33023310        else:
    33033311            Plot.set_ylabel('Parameter values')
    3304         Page.canvas.draw()           
     3312        Page.canvas.draw()
    33053313           
    33063314    G2frame.seqXselect = SelectX
     
    33273335        Page.canvas.mpl_connect('key_press_event', OnKeyPress)
    33283336        Page.canvas.mpl_connect('motion_notify_event', OnMotion)
    3329     Page.Choice = ['s - select x-axis','t - change titles',]
     3337    Page.Choice = ['l - toggle lines','s - select x-axis','t - change titles',]
    33303338    Page.keyPress = OnKeyPress
    33313339    Page.seqYaxisList = ColumnList
Note: See TracChangeset for help on using the changeset viewer.