Changeset 1339


Ignore:
Timestamp:
May 9, 2014 2:56:31 PM (8 years ago)
Author:
vondreele
Message:

make single pattern plot the default
begin slit smearing implementation - doesn't work just yet
some tweaking of pattern plotting

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1334 r1339  
    19031903        self.Contour = False
    19041904        self.Legend = False
    1905         self.SinglePlot = False
     1905        self.SinglePlot = True
    19061906        self.SubBack = False
    19071907        self.seqReverse = False
     
    32193219                        Id = item
    32203220                    item, cookie = self.PatternTree.GetNextChild(self.root, cookie)               
     3221                if Id:
     3222                    self.PatternTree.SelectItem(Id)
    32213223                if parentName:
    32223224                    parentId = G2gd.GetPatternTreeItemId(self, self.root, parentName)
     
    32263228                        itemId = G2gd.GetPatternTreeItemId(self, self.root, oldName)
    32273229                    self.PatternTree.SelectItem(itemId)
    3228                 elif Id:
    3229                     self.PatternTree.SelectItem(Id)
    32303230        finally:
    32313231            dlg.Destroy()
  • trunk/GSASIIplot.py

    r1334 r1339  
    444444        newPlot = False
    445445        if event.key == 'w' and 'PWDR' in plottype:
    446             G2frame.ErrorBars = not G2frame.ErrorBars
     446            G2frame.Weight = not G2frame.Weight
    447447            if not G2frame.Weight:
    448448                G2frame.SinglePlot = True
     
    928928                X += G2frame.Offset[1]*.005*N
    929929            Xum = ma.getdata(X)
     930           
    930931            if ifpicked:
    931932                Z = xye[3]+offset*N
    932                 if 'PWDR' in plottype:  #powder background
     933                if 'PWDR' in plottype:
    933934                    W = xye[4]+offset*N
    934                     D = xye[5]-Ymax*G2frame.delOffset
     935                    D = xye[5]-Ymax*G2frame.delOffset  #powder background
    935936                elif 'SASD' in plottype:
    936937                    if G2frame.sqPlot:
     
    947948                        ZB = Z+B
    948949                    Plot.set_yscale("log",nonposy='mask')
    949                     Plot.set_ylim(bottom=np.min(np.trim_zeros(Y))/2.,top=np.max(Y)*2.)
     950                    Plot.set_ylim(bottom=np.min(np.trim_zeros(YB))/2.,top=np.max(Y)*2.)
    950951                if G2frame.logPlot:
    951952                    if 'PWDR' in plottype:
     
    10081009                    Page.canvas.SetToolTipString(tip)
    10091010                    data = G2frame.LimitsTable.GetData()
    1010             else:
     1011                   
     1012            else:   #not picked
    10111013                if G2frame.logPlot:
    10121014                    if 'PWDR' in plottype:
     
    10201022                        Plot.loglog(X,Y,colors[N%6],picker=False,nonposy='mask')
    10211023                        Plot.set_ylim(bottom=np.min(np.trim_zeros(Y))/2.,top=np.max(Y)*2.)
    1022             if G2frame.logPlot:
     1024                       
     1025            if G2frame.logPlot and 'PWDR' in plottype:
    10231026                Plot.set_ylim(bottom=np.min(np.trim_zeros(Y))/2.,top=np.max(Y)*2.)
    10241027    if PickId and not G2frame.Contour:
  • trunk/GSASIIpwdGUI.py

    r1337 r1339  
    8484        'Thick':1.0,'Contrast':[0.0,0.0],       #contrast & anomalous contrast
    8585        'Trans':1.0,                            #measured transmission
    86         'SlitLen':0.0,                          #Slit length - units?
     86        'SlitLen':0.0,                          #Slit length - in Q(A-1)
    8787        }
    8888def SetupSampleLabels(histName,dataType):
     
    108108        parms.append(['Thick','Sample thickness (mm)',[10,3]])
    109109        parms.append(['Trans','Transmission (meas)',[10,3]])
    110         parms.append(['SlitLen','Slit length',[10,3]])
     110        parms.append(['SlitLen',u'Slit length (Q,\xc5'+Pwrm1+')',[10,3]])
    111111    parms.append(['Omega','Goniometer omega:',[10,3]])
    112112    parms.append(['Chi','Goniometer chi:',[10,3]])
  • trunk/GSASIIsasd.py

    r1337 r1339  
    10011001    return (G*P/Rg**P)*np.exp(scsp.gammaln(P/2))
    10021002   
    1003 def SmearData(Ic,Q,slitLength):
    1004     return Ic   #for now
     1003def SmearData(Ic,Q,slitLen):
     1004    Qtemp = np.concatenate([Q,20*Q])
     1005    Ictemp = np.concatenate([Ic,np.zeros_like(Ic)])
     1006    print Ictemp
     1007    Icsm = np.zeros_like(Qtemp)
     1008    Np = Q.shape[0]
     1009    Qsm = 2*slitLen*(np.interp(np.arange(2*Np)/2.,np.arange(Np),Q)-Q[0])/(Q[-1]-Q[0])
     1010    Sp = np.searchsorted(Qsm,slitLen)
     1011    for i in range(Np):
     1012        Ism = np.interp(np.sqrt(Q[i]**2+Qsm**2),Qtemp,Ictemp)
     1013        print Ism
     1014        raise Exception
     1015        Icsm[i] = np.sum(Ism[:Sp])
     1016    Icsm /= slitLen
     1017    return Icsm   #for now
    10051018   
    10061019###############################################################################
     
    12101223                    parmDict[cid+'PkSig'],parmDict[cid+'PkGam'],Q)
    12111224        Ic += parmDict['Back']  #/parmDict['Scale']
     1225        slitLen = Sample['SlitLen']
     1226        if slitLen:
     1227            Ic = SmearData(Ic,Q,slitLen)
    12121228        return Ic
    12131229       
     
    13591375            Dist.append([])
    13601376    Ic[Ibeg:Ifin] += Back[0]
     1377    slitLen = Sample['SlitLen']
     1378    if slitLen:
     1379        Ic[Ibeg:Ifin] = SmearData(Ic,Q,slitLen)[Ibeg:Ifin]
    13611380    sasdData['Size Calc'] = [Rbins,Dist]
    13621381   
Note: See TracChangeset for help on using the changeset viewer.