Changeset 1402


Ignore:
Timestamp:
Jul 1, 2014 3:29:14 PM (8 years ago)
Author:
vondreele
Message:

make pseudovariables work for strain seq. results
fix Pos2dsp for TOF
allow plot of TOF vs d-spacing for indexed peaks - might be useful for getting TOF to d-space coefficients - uses PlotXY

Location:
trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIgrid.py

    r1400 r1402  
    44924492            parmDict = data[name].get('parmDict')
    44934493            if parmDict:
    4494                 constraintInfo = data[name].get('constraintInfo')
     4494                constraintInfo = data[name].get('constraintInfo',[[],[],{},[],seqnum])
    44954495                groups,parmlist,constrDict,fixedList,ihst = constraintInfo
    44964496                varyList = data[name]['varyList']
  • trunk/GSASIIimgGUI.py

    r1399 r1402  
    15611561                        variables += item['Emat']
    15621562                        sig += item['Esig']
    1563                         varyList += ['%d%s%s'%(i,':',Name) for Name in varyNames]
     1563                        varylist = ['%d%s%s'%(i,':',Name) for Name in varyNames]
     1564                        varyList += varylist
     1565                        parmDict.update(dict(zip(varylist,item['Emat'])))
    15641566                        parmDict['%d:Dcalc'%(i)] = item['Dcalc']
    15651567                    SeqResult[name] = {'variables':variables,'varyList':varyList,'sig':sig,'Rvals':[],
    1566                         'covMatrix':[],'title':name,'parmDict':parmDict}
     1568                        'covMatrix':np.eye(len(variables)),'title':name,'parmDict':parmDict}
    15671569                else:
    15681570                    SeqResult['histNames'] = goodnames
  • trunk/GSASIIlattice.py

    r1372 r1402  
    423423        dsp = wave/(2.0*sind((pos-Inst['Zero'][1])/2.0))
    424424    else:   #'T'OF - ignore difA
    425         dsp = (pos-Inst['Zero'][1])/difC
     425        dsp = (pos-Inst['Zero'][1])/Inst['difC'][1]
    426426    return dsp
    427427   
  • trunk/GSASIIobj.py

    r1390 r1402  
    13161316        'PkGam': 'Bragg peak gamma',
    13171317        'e([12][12])' : 'strain tensor e\1',   # strain vars e11, e22, e12
     1318        'Dcalc': 'Calc. d-spacing',
    13181319        }.items():
    13191320        VarDesc[key] = value
  • trunk/GSASIIplot.py

    r1401 r1402  
    14471447################################################################################
    14481448           
    1449 def PlotXY(G2frame,XY,XY2=None,labelX=None,labelY=None,newPlot=False,type=''):
     1449def PlotXY(G2frame,XY,XY2=None,labelX=None,labelY=None,newPlot=False,Title=''):
    14501450    '''simple plot of xy data, used for diagnostic purposes
    14511451    '''
     
    14561456            Page.canvas.SetCursor(wx.CROSS_CURSOR)
    14571457            try:
    1458                 G2frame.G2plotNB.status.SetStatusText('X =%9.3f %s =%9.3f'%(xpos,type,ypos),1)                   
     1458                G2frame.G2plotNB.status.SetStatusText('X =%9.3f %s =%9.3f'%(xpos,Title,ypos),1)                   
    14591459            except TypeError:
    1460                 G2frame.G2plotNB.status.SetStatusText('Select '+type+' pattern first',1)
     1460                G2frame.G2plotNB.status.SetStatusText('Select '+Title+' pattern first',1)
    14611461
    14621462    try:
    1463         plotNum = G2frame.G2plotNB.plotList.index(type)
     1463        plotNum = G2frame.G2plotNB.plotList.index(Title)
    14641464        Page = G2frame.G2plotNB.nb.GetPage(plotNum)
    14651465        if not newPlot:
     
    14701470    except ValueError:
    14711471        newPlot = True
    1472         Plot = G2frame.G2plotNB.addMpl(type).gca()
    1473         plotNum = G2frame.G2plotNB.plotList.index(type)
     1472        Plot = G2frame.G2plotNB.addMpl(Title).gca()
     1473        plotNum = G2frame.G2plotNB.plotList.index(Title)
    14741474        Page = G2frame.G2plotNB.nb.GetPage(plotNum)
    14751475        Page.canvas.mpl_connect('motion_notify_event', OnMotion)
     
    14781478    Page.SetFocus()
    14791479    G2frame.G2plotNB.status.DestroyChildren()
    1480     Plot.set_title(type)
    1481     if xLabel:
    1482         Plot.set_xlabel(r''+xLabel,fontsize=14)
     1480    Plot.set_title(Title)
     1481    if labelX:
     1482        Plot.set_xlabel(r''+labelX,fontsize=14)
    14831483    else:
    14841484        Plot.set_xlabel(r'X',fontsize=14)
    1485     if yLabel:
    1486         Plot.set_ylabel(r''+yLabel,fontsize=14)
     1485    if labelY:
     1486        Plot.set_ylabel(r''+labelY,fontsize=14)
    14871487    else:
    14881488        Plot.set_ylabel(r'Y',fontsize=14)
     
    14901490    for ixy,xy in enumerate(XY):
    14911491        X,Y = xy
    1492         Plot.plot(X,Y,color(ixy%6)+'+',picker=False)
     1492        Plot.plot(X,Y,colors[ixy%6]+'+',picker=False)
    14931493    if len(XY2):
    14941494        for ixy,xy in enumerate(XY2):
    14951495            X,Y = xy
    1496             Plot.plot(X,Y,color(ixy%6),picker=False)
     1496            Plot.plot(X,Y,colors[ixy%6],picker=False)
    14971497    if not newPlot:
    14981498        Page.toolbar.push_current()
  • trunk/GSASIIpwdGUI.py

    r1391 r1402  
    17901790    G2frame.dataDisplay = G2gd.GSGrid(parent=G2frame.dataFrame)               
    17911791    G2frame.dataDisplay.SetTable(G2frame.IndexPeaksTable, True)
     1792    XY = []
    17921793    for r in range(G2frame.dataDisplay.GetNumberRows()):
    17931794        for c in range(G2frame.dataDisplay.GetNumberCols()):
     
    17961797            else:
    17971798                G2frame.dataDisplay.SetReadOnly(r,c,isReadOnly=True)
     1799            if 'PNT' in Inst['Type'][0] and data[r][3]:
     1800                X = data[r][8]*Inst['difC'][1]
     1801                Y = data[r][0]
     1802                XY.append([X,Y])
    17981803    G2frame.dataDisplay.Bind(wg.EVT_GRID_CELL_LEFT_CLICK, RefreshIndexPeaksGrid)
    17991804    G2frame.dataDisplay.Bind(wx.EVT_KEY_DOWN, KeyEditPickGrid)                 
     
    18011806    G2frame.dataDisplay.AutoSizeColumns(False)
    18021807    G2frame.dataFrame.setSizePosLeft([490,300])
    1803  
     1808    if len(XY): #NB: only for TOF
     1809        G2plt.PlotXY(G2frame,np.array(XY),XY2=[],labelX='T-calc',labelY='Tobs-Tcalc',
     1810            newPlot=True,Title='Diffractometer const')
     1811     
    18041812################################################################################
    18051813#####  Unit cells
Note: See TracChangeset for help on using the changeset viewer.