Changeset 1439


Ignore:
Timestamp:
Jul 24, 2014 9:39:05 AM (7 years ago)
Author:
vondreele
Message:

add difB as a term in d3 for TOF - trial
a patch for old files for missing mapDataFlip? flag
in Unit cells list the Show hkl pos now updates the Index peak list; this can be useful for trying out unit cells by hand, e.g. for contaminating phases, etc.

Location:
trunk
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1420 r1439  
    984984                return [G2IO.makeInstDict(names,data,codes),{}]
    985985            elif 'T' in DataType:
    986                 names = ['Type','2-theta','difC','difA','Zero','alpha','beta-0','beta-1',
     986                names = ['Type','2-theta','difC','difA','difB','Zero','alpha','beta-0','beta-1',
    987987                    'beta-q','sig-0','sig-1','sig-q','X','Y','Azimuth']
    988                 codes = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
     988                codes = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
    989989                azm = 0.
    990990                if 'INS  1DETAZM' in Iparm:
     
    993993                data.extend([G2IO.sfloat(s[1]),])               #2-theta for bank
    994994                s = Iparm['INS  1 ICONS'].split()
    995                 data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),G2IO.sfloat(s[2])])    #difC, difA, Zero
     995                data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),0.0,G2IO.sfloat(s[2])])    #difC,difA,difB,Zero
    996996                if 'INS  1PRCF  ' in Iparm:
    997997                    s = Iparm['INS  1PRCF  '].split()
  • trunk/GSASIIlattice.py

    r1402 r1439  
    417417def Pos2dsp(Inst,pos):
    418418    ''' convert powder pattern position (2-theta or TOF, musec) to d-spacing
    419     ignores secondary effects (e.g. difA in TOF)
     419    ignores secondary effects (e.g. difA,difB in TOF)
    420420    '''
    421421    if 'C' in Inst['Type'][0]:
    422422        wave = G2mth.getWave(Inst)
    423423        dsp = wave/(2.0*sind((pos-Inst['Zero'][1])/2.0))
    424     else:   #'T'OF - ignore difA
     424    else:   #'T'OF - ignore difA, difB
    425425        dsp = (pos-Inst['Zero'][1])/Inst['difC'][1]
    426426    return dsp
     
    433433        wave = G2mth.getWave(Inst)
    434434        pos = 2.0*asind(wave/(2.*dsp))+Inst['Zero'][1]             
    435     else:   #'T'OF - ignore difA
    436         pos = Inst['difC'][1]*dsp+Inst['Zero'][1]
     435    else:   #'T'OF
     436        pos = Inst['difC'][1]*dsp+Inst['Zero'][1]+Inst['difA'][1]*dsp**2+Inst.get('difB',[0,0,False])[1]*dsp**3
    437437    return pos             
    438438   
  • trunk/GSASIIplot.py

    r1420 r1439  
    42054205    Page.SetFocus()
    42064206    Page.Choice = None
    4207     if mapData['Flip']:
     4207    if mapData.get('Flip',False):
    42084208        choice = [' save as/key:','jpeg','tiff','bmp','c: center on 1/2,1/2,1/2',
    42094209            'u: roll up','d: roll down','l: roll left','r: roll right']
  • trunk/GSASIIpwd.py

    r1415 r1439  
    12351235                        parmDict[sigName] = G2mth.getCWsig(parmDict,pos)
    12361236                    else:
    1237                         dsp = G2mth.Pos2dsp(Inst,pos)
     1237                        dsp = G2lat.Pos2dsp(Inst,pos)
    12381238                        parmDict[sigName] = G2mth.getTOFsig(parmDict,dsp)
    12391239                gamName = 'gam'+str(iPeak)
     
    12421242                        parmDict[gamName] = G2mth.getCWgam(parmDict,pos)
    12431243                    else:
    1244                         dsp = G2mth.Pos2dsp(Inst,pos)
     1244                        dsp = G2lat.Pos2dsp(Inst,pos)
    12451245                        parmDict[gamName] = G2mth.getTOFgamma(parmDict,dsp)
    12461246                iPeak += 1
  • trunk/GSASIIpwdGUI.py

    r1428 r1439  
    909909    G2frame.dataDisplay.SetMargins(0,0)
    910910    G2frame.dataDisplay.AutoSizeColumns(False)
    911     G2frame.dataFrame.setSizePosLeft([230,260])
    912                                
     911    G2frame.dataFrame.setSizePosLeft([230,260])                               
    913912   
    914913################################################################################
     
    925924            if key in ['Type','U','V','W','X','Y','SH/L','I(L2)/I(L1)','alpha',
    926925                'beta-0','beta-1','beta-q','sig-0','sig-1','sig-q','Polariz.',
    927                 'Lam','Azimuth','2-theta','difC','difA','Zero','Lam1','Lam2']:
     926                'Lam','Azimuth','2-theta','difC','difA','difB','Zero','Lam1','Lam2']:
    928927                good.append(key)
    929928        return good
     
    12421241                    subSizer.Add(wx.StaticText(G2frame.dataDisplay,-1,'  alpha, beta: fixed by table'),0,WACV)
    12431242                else:
    1244                     Items = ['difC','difA','Zero','alpha','beta-0','beta-1','beta-q','sig-0','sig-1','sig-q','X','Y']
     1243                    Items = ['difC','difA','difB','Zero','alpha','beta-0','beta-1','beta-q','sig-0','sig-1','sig-q','X','Y']
    12451244                mainSizer.Add(subSizer)
    12461245                for item in Items:
     
    12591258                    elemKeysLst.append([item,1])
    12601259                    dspLst.append(nDig)
    1261                     if not ifHisto and item in ['difC','difA','Zero',]:
     1260                    if not ifHisto and item in ['difC','difA','difB','Zero',]:
    12621261                        refFlgElem.append(None)
    12631262                        instSizer.Add((5,5),0)
     
    13491348                insDef['Azimuth'] = 0.0
    13501349                insRef['Azimuth'] = False
     1350        if 'T' in insVal['Type']:
     1351            if 'difB' not in insVal:
     1352                insVal['difB'] = 0.0
     1353                insDef['difB'] = 0.0
     1354                insRef['difB'] = False
    13511355    #end of patch
    13521356    if 'P' in insVal['Type']:                   #powder data menu commands
     
    18611865                G2frame.dataDisplay.SetReadOnly(r,c,isReadOnly=True)
    18621866            if 'PNT' in Inst['Type'][0] and data[r][3]:
    1863                 X = data[r][8]*Inst['difC'][1]
     1867                X = G2lat.Dsp2pos(Inst,data[r][8])
    18641868                Y = data[r][0]
    1865                 XY.append([X,Y])
     1869                XY.append([X,Y-X])
    18661870    G2frame.dataDisplay.Bind(wg.EVT_GRID_CELL_LEFT_CLICK, RefreshIndexPeaksGrid)
    18671871    G2frame.dataDisplay.Bind(wx.EVT_KEY_DOWN, KeyEditPickGrid)                 
     
    19962000        PatternId = G2frame.PatternId
    19972001        PickId = G2frame.PickId   
     2002        peaks = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId, 'Index Peak List'))
    19982003        limits = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId, 'Limits'))[1]
    19992004        controls,bravais,cells,dmin = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId, 'Unit Cells List'))
     
    20052010        if 'C' in Inst['Type'][0]:
    20062011            dmin = G2lat.Pos2dsp(Inst,limits[1])
    2007         else:
     2012        else:   #TOF - use other limit!
    20082013            dmin = G2lat.Pos2dsp(Inst,limits[0])
    20092014        G2frame.HKL = G2pwd.getHKLpeak(dmin,SGData,A)
     2015        G2indx.IndexPeaks(peaks,G2frame.HKL)
    20102016        for hkl in G2frame.HKL:
    20112017            hkl.append(G2lat.Dsp2pos(Inst,hkl[3])++controls[1])
Note: See TracChangeset for help on using the changeset viewer.