Changeset 837 for trunk/GSASIIplot.py


Ignore:
Timestamp:
Jan 22, 2013 7:18:17 PM (9 years ago)
Author:
vondreele
Message:

more on rigid bodies - imput of text file of coordinates
work on HKL plot
add import of F2 from hkl file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r836 r837  
    243243            self.G2plotNB.status.SetFields(['','HKL, Fosq, sig, Fcsq, delFsq/sig = '+HKLtxt])
    244244                                 
    245     def OnSCKeyPress(event):
    246         print event.key
    247 
    248245    try:
    249246        plotNum = self.G2plotNB.plotList.index('Structure Factors')
     
    259256        plotNum = self.G2plotNB.plotList.index('Structure Factors')
    260257        Page = self.G2plotNB.nb.GetPage(plotNum)
    261 #        Page.canvas.mpl_connect('key_press_event', OnSCKeyPress)
    262258        Page.canvas.mpl_connect('pick_event', OnSCPick)
    263259        Page.canvas.mpl_connect('motion_notify_event', OnSCMotion)
     
    275271    FosqMax = Data['FoMax']
    276272    FoMax = math.sqrt(FosqMax)
    277     ifFc = Data['ifFc']
    278273    xlabel = ['k, h=','h, k=','h, l=']
    279274    ylabel = ['l','l','k']
     
    286281    for refl in HKLref:
    287282        H = np.array(refl[:3])
    288         sig,Fosq,Fcsq = refl[7:10]
     283        Fosq,sig,Fcsq = refl[5:8]
    289284        HKL.append(H)
    290285        HKLF.append([Fosq,sig,Fcsq])
    291286        if H[izone] == Data['Layer']:
     287            A = 0
    292288            B = 0
    293289            if Type == 'Fosq':
     
    300296                C = abs(A-B)
    301297            elif Type == '|DFsq|/sig':
    302                 A = abs(Fosq-Fcsq)/(scale*sig)
     298                if sig > 0.:
     299                    A = abs(Fosq-Fcsq)/(scale*sig)
     300                B = 0
    303301            elif Type == '|DFsq|>sig':
    304                 A = abs(Fosq-Fcsq)/(scale*sig)
    305                 if A < 1.0: A = 0                   
     302                if sig > 0.:
     303                    A = abs(Fosq-Fcsq)/sig
     304                if A < 1.0: A = 0
     305                A /= scale
     306                B = 0                   
    306307            elif Type == '|DFsq|>3sig':
    307                 A = abs(Fosq-Fcsq)/(scale*sig)
    308                 if A < 3.0: A = 0                   
     308                if sig > 0.:
     309                    A = abs(Fosq-Fcsq)/sig
     310                if A < 3.0: A = 0
     311                A /= scale
     312                B = 0                   
    309313            xy = (H[pzone[izone][0]],H[pzone[izone][1]])
    310314            if A > 0.0:
     
    324328    Plot.set_xlim((HKLmin[pzone[izone][0]],HKLmax[pzone[izone][0]]))
    325329    Plot.set_ylim((HKLmin[pzone[izone][1]],HKLmax[pzone[izone][1]]))
    326     if not newPlot:
    327         Page.toolbar.push_current()
    328         Plot.set_xlim(xylim[0])
    329         Plot.set_ylim(xylim[1])
    330         xylim = []
    331         Page.toolbar.push_current()
    332         Page.toolbar.draw()
    333     else:
    334         Page.canvas.draw()
     330    Page.canvas.draw()
     331#    if not newPlot:
     332#        Page.toolbar.push_current()
     333#        Plot.set_xlim(xylim[0])
     334#        Plot.set_ylim(xylim[1])
     335#        xylim = []
     336#        Page.toolbar.push_current()
     337#        Page.toolbar.draw()
     338#    else:
     339#        Page.canvas.draw()
    335340       
    336341################################################################################
     
    17591764    G2frame.G2plotNB.status.SetFields(['',''])   
    17601765    acolor = mpl.cm.get_cmap(G2frame.VcovColor)
    1761     Img = Plot.imshow(covArray,aspect='equal',cmap=acolor,interpolation='nearest',origin='lower')
     1766    Img = Plot.imshow(covArray,aspect='equal',cmap=acolor,interpolation='nearest',origin='lower',
     1767        vmin=-1.,vmax=1.)
    17621768    imgAx = Img.get_axes()
    17631769    ytics = imgAx.get_yticks()
     
    26322638        mode = cb.GetValue()
    26332639        if mode in ['jpeg','bmp','tiff',]:
    2634             Fname = os.path.joint(Mydir,generalData['Name']+'.'+mode)
     2640            Fname = os.path.join(Mydir,generalData['Name']+'.'+mode)
    26352641            size = Page.canvas.GetSize()
    26362642            glPixelStorei(GL_UNPACK_ALIGNMENT, 1)
     
    36023608       
    36033609    try:
    3604         plotNum = G2frame.G2plotNB.plotList.index(rbData['RBname'])
     3610        plotNum = G2frame.G2plotNB.plotList.index('Rigid body')
    36053611        Page = G2frame.G2plotNB.nb.GetPage(plotNum)       
    36063612    except ValueError:
    3607         Plot = G2frame.G2plotNB.addOgl(rbData['RBname'])
    3608         plotNum = G2frame.G2plotNB.plotList.index(rbData['RBname'])
     3613        Plot = G2frame.G2plotNB.addOgl('Rigid body')
     3614        plotNum = G2frame.G2plotNB.plotList.index('Rigid body')
    36093615        Page = G2frame.G2plotNB.nb.GetPage(plotNum)
    36103616        Page.views = False
Note: See TracChangeset for help on using the changeset viewer.