source: sasview/sansguiframe/src/sans/guiframe/local_perspectives/plotting/Plotter2D.py @ cce182ae

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since cce182ae was 2a09236, checked in by Jae Cho <jhjcho@…>, 13 years ago

tick label works w/ 2D too.

  • Property mode set to 100644
File size: 28.1 KB
RevLine 
[1bf33c1]1
[d955bf19]2################################################################################
3#This software was developed by the University of Tennessee as part of the
4#Distributed Data Analysis of Neutron Scattering Experiments (DANSE)
5#project funded by the US National Science Foundation.
6#
7#See the license text in license.txt
8#
9#copyright 2008, University of Tennessee
10################################################################################
[1bf33c1]11
12
13import wx
[32c0841]14import sys
[ed8ad21]15import os
[32c0841]16import math
[003fa4e]17import numpy
[0d9dae8]18import pylab
[1bf33c1]19import danse.common.plottools
20from danse.common.plottools.PlotPanel import PlotPanel
[4ac8556]21from danse.common.plottools.plottables import Graph
[49a470d]22from danse.common.plottools.TextDialog import TextDialog
[55a0dc1]23from sans.guiframe.events import EVT_NEW_PLOT
24from sans.guiframe.events import EVT_SLICER_PARS
25from sans.guiframe.events import StatusEvent
26from sans.guiframe.events import NewPlotEvent
[a07e72f]27from sans.guiframe.events import PanelOnFocusEvent
[55a0dc1]28from sans.guiframe.events import SlicerEvent
[0d9dae8]29from sans.guiframe.utils import PanelMenu
[1bf33c1]30from binder import BindArtist
31from Plotter1D import ModelPanel1D
[32c0841]32from danse.common.plottools.toolbar import NavigationToolBar
[4ac8556]33from sans.guiframe.dataFitting import Data1D
[49a470d]34from matplotlib.font_manager import FontProperties
[32c0841]35(InternalEvent, EVT_INTERNAL) = wx.lib.newevent.NewEvent()
[0d9dae8]36
37DEFAULT_QMAX = 0.05
[1bf33c1]38DEFAULT_QSTEP = 0.001
39DEFAULT_BEAM = 0.005
[ef0c170]40BIN_WIDTH = 1.0
[32c0841]41
[d955bf19]42
[ac8671e]43class NavigationToolBar2D(NavigationToolBar):
[d955bf19]44    """
45    """
[ac8671e]46    def __init__(self, canvas, parent=None):
47        NavigationToolBar.__init__(self, canvas=canvas, parent=parent)
48       
49    def delete_option(self):
50        """
[d955bf19]51        remove default toolbar item
[ac8671e]52        """
53        #delete reset button
54        self.DeleteToolByPos(0) 
55        #delete dragging
[88ca6db]56        self.DeleteToolByPos(2) 
[ac8671e]57        #delete unwanted button that configures subplot parameters
58        self.DeleteToolByPos(4)
59       
60    def add_option(self):
61        """
[d955bf19]62        add item to the toolbar
[ac8671e]63        """
[dc51a7f]64        #add button
65        id_context = wx.NewId()
66        context_tip = 'Graph Menu'
67        context =  wx.ArtProvider.GetBitmap(wx.ART_LIST_VIEW, wx.ART_TOOLBAR)
68        self.InsertSimpleTool(0, id_context, context, 
69                                   context_tip, context_tip)
70        wx.EVT_TOOL(self, id_context,  self.parent.onToolContextMenu)
71        self.InsertSeparator(1)
[ac8671e]72        #add print button
73        id_print = wx.NewId()
74        print_bmp =  wx.ArtProvider.GetBitmap(wx.ART_PRINT, wx.ART_TOOLBAR)
75        self.AddSimpleTool(id_print, print_bmp,
76                           'Print', 'Activate printing')
77        wx.EVT_TOOL(self, id_print, self.on_print)
78       
79       
80class ModelPanel2D(ModelPanel1D):
[1bf33c1]81    """
[d955bf19]82    Plot panel for use with the GUI manager
[1bf33c1]83    """
84   
85    ## Internal name for the AUI manager
86    window_name = "plotpanel"
87    ## Title to appear on top of the window
88    window_caption = "Plot Panel"
89    ## Flag to tell the GUI manager that this panel is not
90    #  tied to any perspective
91    ALWAYS_ON = True
92    ## Group ID
93    group_id = None
94   
95   
[32c0841]96    def __init__(self, parent, id=-1, data2d=None, color = None,
97                 dpi=None, style=wx.NO_FULL_REPAINT_ON_RESIZE, **kwargs):
[1bf33c1]98        """
[d955bf19]99        Initialize the panel
[1bf33c1]100        """
[32c0841]101        ModelPanel1D.__init__(self, parent, id=id, style=style, **kwargs)
[1bf33c1]102       
103        ## Reference to the parent window
104        self.parent = parent
[6c0568b]105        ## Dictionary containing Plottables
[1bf33c1]106        self.plots = {}
[6c0568b]107        ## Save reference of the current plotted
108        self.data2D = data2d
[1bf33c1]109        ## Unique ID (from gui_manager)
110        self.uid = None
111        ## Action IDs for internal call-backs
112        self.action_ids = {}
[6c0568b]113        ## Create Artist and bind it
[1bf33c1]114        self.connect = BindArtist(self.subplot.figure)
[6c0568b]115        ## Beam stop
[1bf33c1]116        self.beamstop_radius = DEFAULT_BEAM
[6c0568b]117        ## to set the order of lines drawn first.
[f15ed33]118        self.slicer_z = 5
[6c0568b]119        ## Reference to the current slicer
[1bf33c1]120        self.slicer = None
[6c0568b]121        ## event to send slicer info
[d468daa]122        self.Bind(EVT_INTERNAL, self._onEVT_INTERNAL)
[1bf33c1]123       
[6c0568b]124        self.axes_frozen = False
125        ## panel that contains result from slicer motion (ex: Boxsum info)
[32c0841]126        self.panel_slicer = None
[49a470d]127        self.title_label = None
128        self.title_font = None
129        self.title_color = 'black'
[1bf33c1]130        ## Graph       
131        self.graph = Graph()
132        self.graph.xaxis("\\rm{Q}", 'A^{-1}')
[32c0841]133        self.graph.yaxis("\\rm{Intensity} ", "cm^{-1}")
[1bf33c1]134        self.graph.render(self)
[8dfdd20]135        ## store default value of zmin and zmax
136        self.default_zmin_ctl = self.zmin_2D
137        self.default_zmax_ctl = self.zmax_2D
[ac8671e]138       
[a07e72f]139    def onLeftDown(self, event): 
140        """
141        left button down and ready to drag
142       
143        """
144        # Check that the LEFT button was pressed
[53cf669]145        PlotPanel.onLeftDown(self, event)   
[e85e2bc]146        ax = event.inaxes
147        if ax != None:
148            # data coordinate position
149            pos_x = "%8.3g"% event.xdata
150            pos_y = "%8.3g"% event.ydata
151            position = "x: %s    y: %s" % (pos_x, pos_y)
152            wx.PostEvent(self.parent, StatusEvent(status=position))
[a07e72f]153        self.plottable_selected(self.data2D.id)
154        self._manager.set_panel_on_focus(self)
155        wx.PostEvent(self.parent, PanelOnFocusEvent(panel=self))
156       
[ac8671e]157    def add_toolbar(self):
158        """
[d955bf19]159        add toolbar
[ac8671e]160        """
161        self.enable_toolbar = True
[32c0841]162        self.toolbar = NavigationToolBar2D(parent=self, canvas=self.canvas)
[ac8671e]163        self.toolbar.Realize()
164        # On Windows platform, default window size is incorrect, so set
165        # toolbar width to figure width.
166        tw, th = self.toolbar.GetSizeTuple()
167        fw, fh = self.canvas.GetSizeTuple()
168        self.toolbar.SetSize(wx.Size(fw, th))
[32c0841]169        self.sizer.Add(self.toolbar, 0, wx.LEFT|wx.EXPAND)
[ac8671e]170        # update the axes menu on the toolbar
171        self.toolbar.update()
172         
[a07e72f]173    def plot_data(self, data):
[1bf33c1]174        """
[d955bf19]175        Data is ready to be displayed
176       
177        :TODO: this name should be changed to something more appropriate
178             Don't forget that changing this name will mean changing code
179             in plotting.py
180         
181        :param event: data event
[1bf33c1]182        """
[e4a703a]183        xlo = None
184        xhi = None
185        ylo = None 
186        yhi = None
[b837e1f]187        if data.__class__.__name__ == 'Data1D':
188            return
[6c0568b]189        ## Update self.data2d with the current plot
[857d00f]190        self.data2D = data
[a07e72f]191        if data.id in self.plots.keys():
192            #replace
[e4a703a]193            xlo, xhi = self.subplot.get_xlim()
194            ylo, yhi = self.subplot.get_ylim()
[a07e72f]195            self.graph.replace(data)
196            self.plots[data.id] = data
[ab8f936]197        else:
[a07e72f]198            self.plots[data.id] = data
199            self.graph.add(self.plots[data.id]) 
[ac9a5f6]200            # update qmax with the new xmax of data plotted
[a07e72f]201            self.qmax = data.xmax
[4b91fd1]202           
[32c0841]203        self.slicer = None
[1bf33c1]204        # Check axis labels
205        #TODO: Should re-factor this
[6c0568b]206        ## render the graph with its new content
[7fff5cd]207               
208        #data2D: put 'Pixel (Number)' for axis title and unit in case of having no detector info and none in _units
[49a470d]209        if len(data.detector) < 1: 
[a07e72f]210            if len(data._xunit) < 1 and len(data._yunit) < 1:
211                data._xaxis = '\\rm{x}'
212                data._yaxis = '\\rm{y}'
213                data._xunit = 'pixel'
214                data._yunit = 'pixel'
[857d00f]215       
216        # graph properties
[a07e72f]217        self.graph.xaxis(data._xaxis, data._xunit)
218        self.graph.yaxis(data._yaxis, data._yunit)
[36288ca]219        if self._is_changed_legend_label:   
220                data.label = self.title_label
[49a470d]221       
[fe48fcc]222        if data.label == None:
223            data.label = data.name
[49a470d]224           
225        if not self.title_font:
226            self.graph.title(data.label)
227            self.graph.render(self)
[857d00f]228            # Set the axis labels on subplot
229            self._set_axis_labels()
[49a470d]230            self.draw_plot()
231        else:
232            self.graph.render(self)
233            self.draw_plot()
234            self.subplot.set_title(label=data.label,
235                                   fontproperties=self.title_font,
236                                   color=self.title_color)
237            self.subplot.figure.canvas.draw_idle() 
238       
[857d00f]239       
[8dfdd20]240        ## store default value of zmin and zmax
241        self.default_zmin_ctl = self.zmin_2D
242        self.default_zmax_ctl = self.zmax_2D
[7022fdc]243        # Check if zoomed
244        toolbar_zoomed = self.toolbar.GetToolEnabled(self.toolbar._NTB2_BACK)
245        if not self.is_zoomed and not toolbar_zoomed:
246            return
247       
[e4a703a]248        # Recover the x,y limits
249        if (xlo and xhi and ylo and yhi) != None:
250            if (xlo > data.xmin and xhi < data.xmax and\
251                        ylo > data.ymin and yhi < data.ymax):
252                self.subplot.set_xlim((xlo, xhi))     
[fe48fcc]253                self.subplot.set_ylim((ylo, yhi)) 
254            else: 
255                self.toolbar.update()
[7022fdc]256                self._is_zoomed = False
[53cf669]257               
258        # Update Graph menu and help string       
259        pos = self.parent._window_menu.FindItem(self.window_caption)
260        helpString = 'Show/Hide Graph: '
[d74fbff]261        helpString += (' ' + str(data.label) +';')
[53cf669]262        self.parent._window_menu.SetHelpString(pos, helpString)
[a4964b8e]263
[857d00f]264    def _set_axis_labels(self):
265        """
266        Set axis labels
267        """
268        data = self.data2D
269        # control axis labels from the panel itself
270        yname, yunits = data.get_yaxis()
271        if self.yaxis_label != None:
272            yname = self.yaxis_label
273            yunits = self.yaxis_unit
274        else:
275            self.yaxis_label = yname
276            self.yaxis_unit = yunits
277        xname, xunits = data.get_xaxis()
278        if self.xaxis_label != None:
279            xname = self.xaxis_label
280            xunits = self.xaxis_unit
281        else:
282            self.xaxis_label = xname
283            self.xaxis_unit = xunits
[2a09236]284        self.xaxis(xname, xunits, self.xaxis_font, 
285                   self.xaxis_color, self.xaxis_tick)
286        self.yaxis(yname, yunits, self.yaxis_font, 
287                   self.yaxis_color, self.yaxis_tick)
[857d00f]288       
[1bf33c1]289    def onContextMenu(self, event):
290        """
[d955bf19]291        2D plot context menu
292       
293        :param event: wx context event
[15550f4]294       
[d955bf19]295        """
[1bf33c1]296        slicerpop = PanelMenu()
297        slicerpop.set_plots(self.plots)
298        slicerpop.set_graph(self.graph)
[9a585d0]299             
300        id = wx.NewId()
[6d727ae]301        slicerpop.Append(id, '&Save Image')
[9a585d0]302        wx.EVT_MENU(self, id, self.onSaveImage)
303       
304        id = wx.NewId()
[6d727ae]305        slicerpop.Append(id,'&Print Image', 'Print image')
[9a585d0]306        wx.EVT_MENU(self, id, self.onPrint)
307       
[1ce365f8]308        id = wx.NewId()
[6d727ae]309        slicerpop.Append(id,'&Print Preview', 'Print preview')
[1ce365f8]310        wx.EVT_MENU(self, id, self.onPrinterPreview)
[003fa4e]311
312        id = wx.NewId()
313        slicerpop.Append(id, '&Copy to Clipboard', 'Copy to the clipboard')
314        wx.EVT_MENU(self, id, self.OnCopyFigureMenu)
[c5a769e]315        slicerpop.AppendSeparator()
[ed8ad21]316        # saving data
317        plot = self.data2D
318        id = wx.NewId()
319        name = plot.name
[c17760d]320        slicerpop.Append(id, "&Save as a File (DAT)" )
[ed8ad21]321        self.action_ids[str(id)] = plot
322        wx.EVT_MENU(self, id, self._onSave)
323
[9a585d0]324        slicerpop.AppendSeparator()
[7fff5cd]325        if len(self.data2D.detector) == 1:       
[0e13148]326           
[a07e72f]327            item_list = self.parent.get_context_menu(self)
[6d727ae]328            if (not item_list == None) and (not len(item_list) == 0) and\
329                self.data2D.name.split(" ")[0] != 'Residuals': 
330                # The line above; Not for trunk
[32c0841]331                for item in item_list:
332                    try:
333                        id = wx.NewId()
334                        slicerpop.Append(id, item[0], item[1])
335                        wx.EVT_MENU(self, id, item[2])
336                    except:
337                        msg = "ModelPanel1D.onContextMenu: "
338                        msg += "bad menu item  %s"%sys.exc_value
339                        wx.PostEvent(self.parent, StatusEvent(status=msg))
340                        pass
341                slicerpop.AppendSeparator()
[0e13148]342           
[15550f4]343            id = wx.NewId()
[c17760d]344            slicerpop.Append(id, '&Perform Circular Average')
[6d727ae]345            wx.EVT_MENU(self, id, self.onCircular) \
346            # For Masked Data
347            if not plot.mask.all():
348                id = wx.NewId()
[c17760d]349                slicerpop.Append(id, '&Masked Circular Average')
[6d727ae]350                wx.EVT_MENU(self, id, self.onMaskedCircular) 
[15550f4]351            id = wx.NewId()
[c17760d]352            slicerpop.Append(id, '&Sector [Q View]')
[15550f4]353            wx.EVT_MENU(self, id, self.onSectorQ) 
354            id = wx.NewId()
[c17760d]355            slicerpop.Append(id, '&Annulus [Phi View ]')
[15550f4]356            wx.EVT_MENU(self, id, self.onSectorPhi) 
[92c2345]357            id = wx.NewId()
[15550f4]358            slicerpop.Append(id, '&Box Sum')
359            wx.EVT_MENU(self, id, self.onBoxSum) 
360            id = wx.NewId()
[c17760d]361            slicerpop.Append(id, '&Box Averaging in Qx')
[15550f4]362            wx.EVT_MENU(self, id, self.onBoxavgX) 
363            id = wx.NewId()
[c17760d]364            slicerpop.Append(id, '&Box Averaging in Qy')
[15550f4]365            wx.EVT_MENU(self, id, self.onBoxavgY) 
[32c0841]366            if self.slicer != None:
[eba08f1a]367                id = wx.NewId()
[c17760d]368                slicerpop.Append(id, '&Clear Slicer')
[15550f4]369                wx.EVT_MENU(self, id,  self.onClearSlicer) 
[32c0841]370                if self.slicer.__class__.__name__  != "BoxSum":
[15550f4]371                    id = wx.NewId()
372                    slicerpop.Append(id, '&Edit Slicer Parameters')
373                    wx.EVT_MENU(self, id, self._onEditSlicer) 
374            slicerpop.AppendSeparator() 
[fe48fcc]375           
376        id = wx.NewId()
[dc51a7f]377        slicerpop.Append(id, '&Edit Graph Label', 'Edit Graph Label')
[fe48fcc]378        wx.EVT_MENU(self, id, self.onEditLabels)
379        slicerpop.AppendSeparator()
380       
[0e13148]381        id = wx.NewId()
[5637362]382        slicerpop.Append(id, '&Edit Y Axis Label')
[857d00f]383        wx.EVT_MENU(self, id, self._on_yaxis_label)     
384        id = wx.NewId()
[5637362]385        slicerpop.Append(id, '&Edit X Axis Label')
[857d00f]386        wx.EVT_MENU(self, id, self._on_xaxis_label)
[5637362]387        id = wx.NewId()
388        slicerpop.Append(id, '&Toggle Grid On/Off', 'Toggle Grid On/Off')
389        wx.EVT_MENU(self, id, self.onGridOnOff)
[857d00f]390        slicerpop.AppendSeparator()
391       
392        id = wx.NewId()
[52f3c98]393        slicerpop.Append(id, '&2D Color Map')
[32c0841]394        wx.EVT_MENU(self, id, self._onEditDetector)
[857d00f]395        slicerpop.AppendSeparator()
396       
[1bf33c1]397        id = wx.NewId()
[c17760d]398        slicerpop.Append(id, '&Toggle Linear/Log Scale')
[1bf33c1]399        wx.EVT_MENU(self, id, self._onToggleScale) 
[53cf669]400       
401       
402        slicerpop.AppendSeparator()
403        id = wx.NewId()
404        slicerpop.Append(id, '&Window Title')
405        wx.EVT_MENU(self, id, self.onChangeCaption)
406       
[dc51a7f]407        try:
408            pos_evt = event.GetPosition()
409            pos = self.ScreenToClient(pos_evt)
410        except:
411            pos_x, pos_y = self.toolbar.GetPositionTuple()
412            pos = (pos_x, pos_y + 5)
[1bf33c1]413        self.PopupMenu(slicerpop, pos)
[37c36d9]414       
[fe48fcc]415           
416    def onEditLabels(self, event):
417        """
418        Edit legend label
419        """
420        selected_plot = self.plots[self.graph.selected_plottable]
421        label = selected_plot.label
[49a470d]422        dial = TextDialog(None, -1, 'Change Label', label)
[fe48fcc]423        if dial.ShowModal() == wx.ID_OK:
[49a470d]424            try:
425                FONT = FontProperties()
426                newlabel = dial.getText()
427                font = FONT.copy()
428                font.set_size(dial.getSize())
429                font.set_family(dial.getFamily())
430                font.set_style(dial.getStyle())
431                font.set_weight(dial.getWeight())
432                colour = dial.getColor()
433                if len(newlabel) > 0:
[53cf669]434                    # update Label
[49a470d]435                    selected_plot.label = newlabel
[53cf669]436                    self.graph.title(newlabel)
[49a470d]437                    self.title_label = selected_plot.label
438                    self.title_font = font
439                    self.title_color = colour
440                    ## render the graph
441                    self.subplot.set_title(label=self.title_label,
442                                           fontproperties=self.title_font,
443                                           color=self.title_color)
[36288ca]444                    self._is_changed_legend_label = True
[49a470d]445                    self.subplot.figure.canvas.draw_idle() 
446            except:
447                if self.parent != None:
448                    from sans.guiframe.events import StatusEvent
449                    msg= "Add Text: Error. Check your property values..."
450                    wx.PostEvent(self.parent, StatusEvent(status = msg ))
451                else:
452                    raise
[fe48fcc]453        dial.Destroy()
[53cf669]454       
455        # Update Graph menu and help string
456        if self.title_label != None:     
457            pos = self.parent._window_menu.FindItem(self.window_caption)
458            helpString = 'Show/Hide Graph: '
[e2160ab]459            helpString += (' ' + str(self.title_label) +';')
[53cf669]460            self.parent._window_menu.SetHelpString(pos, helpString)
[49a470d]461
[fe48fcc]462       
[ea290ee]463    def _onEditDetector(self, event):
[6d920cd]464        """
[d955bf19]465        Allow to view and edits  detector parameters
466       
467        :param event: wx.menu event
468       
[6d920cd]469        """
[ea290ee]470        import detector_dialog
[8dfdd20]471        dialog = detector_dialog.DetectorDialog(self, -1,base=self.parent,
472                       reset_zmin_ctl =self.default_zmin_ctl,
473                       reset_zmax_ctl = self.default_zmax_ctl,cmap=self.cmap)
[6c0568b]474        ## info of current detector and data2D
[ea290ee]475        xnpts = len(self.data2D.x_bins)
476        ynpts = len(self.data2D.y_bins)
477        xmax = max(self.data2D.xmin, self.data2D.xmax)
478        ymax = max(self.data2D.ymin, self.data2D.ymax)
[32c0841]479        qmax = math.sqrt(math.pow(xmax, 2) + math.pow(ymax, 2))
[ea290ee]480        beam = self.data2D.xmin
[6c0568b]481        ## set dialog window content
[ea290ee]482        dialog.setContent(xnpts=xnpts,ynpts=ynpts,qmax=qmax,
483                           beam=self.data2D.xmin,
484                           zmin = self.zmin_2D,
485                          zmax = self.zmax_2D)
486        if dialog.ShowModal() == wx.ID_OK:
487            evt = dialog.getContent()
488            self.zmin_2D = evt.zmin
489            self.zmax_2D = evt.zmax
[32c0841]490            self.cmap = evt.cmap
[ea290ee]491        dialog.Destroy()
[6c0568b]492        ## Redraw the current image
[32c0841]493        self.image(data=self.data2D.data,
[20b6760]494                   qx_data=self.data2D.qx_data,
495                   qy_data=self.data2D.qy_data,
[ea290ee]496                   xmin= self.data2D.xmin,
497                   xmax= self.data2D.xmax,
498                   ymin= self.data2D.ymin,
499                   ymax= self.data2D.ymax,
500                   zmin= self.zmin_2D,
501                   zmax= self.zmax_2D,
[8dfdd20]502                   cmap= self.cmap,
[32c0841]503                   color=0, symbol=0, label=self.data2D.name)
[ea290ee]504        self.subplot.figure.canvas.draw_idle()
[1bf33c1]505       
506    def freeze_axes(self):
[d955bf19]507        """
508        """
[1bf33c1]509        self.axes_frozen = True
510       
511    def thaw_axes(self):
[d955bf19]512        """
513        """
[1bf33c1]514        self.axes_frozen = False
515       
516    def onMouseMotion(self,event):
[d955bf19]517        """
518        """
[1bf33c1]519        pass
[d955bf19]520   
[1bf33c1]521    def onWheel(self, event):
[d955bf19]522        """
523        """
[6c0568b]524        pass 
525     
[1bf33c1]526    def update(self, draw=True):
527        """
[d955bf19]528        Respond to changes in the model by recalculating the
529        profiles and resetting the widgets.
[1bf33c1]530        """
[6d727ae]531        self.draw_plot()
[1bf33c1]532       
533    def _getEmptySlicerEvent(self):
[6c0568b]534        """
[d955bf19]535        create an empty slicervent
[6c0568b]536        """
[32c0841]537        return SlicerEvent(type=None, params=None, obj_class=None)
[d955bf19]538       
[1bf33c1]539    def _onEVT_INTERNAL(self, event):
540        """
[d955bf19]541        Draw the slicer
542       
543        :param event: wx.lib.newevent (SlicerEvent) containing slicer
[6c0568b]544            parameter
[d955bf19]545           
[1bf33c1]546        """
547        self._setSlicer(event.slicer)
548           
549    def _setSlicer(self, slicer):
[6c0568b]550        """
[d955bf19]551        Clear the previous slicer and create a new one.Post an internal
552        event.
553       
554        :param slicer: slicer class to create
555       
[6c0568b]556        """
557        ## Clear current slicer
[1bf33c1]558        if not self.slicer == None: 
559            self.slicer.clear()           
[6c0568b]560        ## Create a new slicer   
[1bf33c1]561        self.slicer_z += 1
562        self.slicer = slicer(self, self.subplot, zorder=self.slicer_z)
[240c805]563        self.subplot.set_ylim(self.data2D.ymin, self.data2D.ymax)
564        self.subplot.set_xlim(self.data2D.xmin, self.data2D.xmax)
[6c0568b]565        ## Draw slicer
[1bf33c1]566        self.update()
567        self.slicer.update()
[32c0841]568        msg = "Plotter2D._setSlicer  %s"%self.slicer.__class__.__name__
569        wx.PostEvent(self.parent, StatusEvent(status=msg))
[1bf33c1]570        # Post slicer event
571        event = self._getEmptySlicerEvent()
572        event.type = self.slicer.__class__.__name__
573        event.obj_class = self.slicer.__class__
574        event.params = self.slicer.get_params()
[d468daa]575        wx.PostEvent(self, event)
[6d727ae]576       
577    def onMaskedCircular(self, event):
578        """
579        perform circular averaging on Data2D with mask if it exists
580       
581        :param event: wx.menu event
582       
583        """
584        self.onCircular(event, True)
585       
586    def onCircular(self, event, ismask=False):
[1bf33c1]587        """
[d955bf19]588        perform circular averaging on Data2D
589       
590        :param event: wx.menu event
591       
[1bf33c1]592        """
[003fa4e]593        # Find the best number of bins
594        npt = math.sqrt(len(self.data2D.data[numpy.isfinite(self.data2D.data)]))
595        npt = math.floor(npt)
[fe857e2]596        from sans.dataloader.manipulations import CircularAverage
[6c0568b]597        ## compute the maximum radius of data2D
[32c0841]598        self.qmax = max(math.fabs(self.data2D.xmax), 
599                        math.fabs(self.data2D.xmin))
600        self.ymax = max(math.fabs(self.data2D.ymax),
601                        math.fabs(self.data2D.ymin))
602        self.radius = math.sqrt(math.pow(self.qmax, 2)+ math.pow(self.ymax, 2)) 
[6c0568b]603        ##Compute beam width
[003fa4e]604        bin_width = (self.qmax + self.qmax)/npt
[6c0568b]605        ## Create data1D circular average of data2D
[32c0841]606        Circle = CircularAverage(r_min=0, r_max=self.radius, 
607                                 bin_width=bin_width)
[6d727ae]608        circ = Circle(self.data2D, ismask=ismask)
[1bf33c1]609        from sans.guiframe.dataFitting import Data1D
[32c0841]610        if hasattr(circ, "dxl"):
611            dxl = circ.dxl
[1bf33c1]612        else:
[32c0841]613            dxl = None
614        if hasattr(circ, "dxw"):
615            dxw = circ.dxw
[1bf33c1]616        else:
[32c0841]617            dxw = None
[003fa4e]618
619        new_plot = Data1D(x=circ.x, y=circ.y, dy=circ.dy, dx=circ.dx)
[32c0841]620        new_plot.dxl = dxl
621        new_plot.dxw = dxw
622        new_plot.name = "Circ avg " + self.data2D.name
623        new_plot.source = self.data2D.source
624        #new_plot.info = self.data2D.info
[1bf33c1]625        new_plot.interactive = True
[32c0841]626        new_plot.detector = self.data2D.detector
[6d727ae]627       
[6c0568b]628        ## If the data file does not tell us what the axes are, just assume...
[32c0841]629        new_plot.xaxis("\\rm{Q}", "A^{-1}")
[2ca51f44]630        if hasattr(self.data2D, "scale") and \
631                    self.data2D.scale == 'linear':
632            new_plot.ytransform = 'y'
633            new_plot.yaxis("\\rm{Residuals} ", "normalized")
[6d727ae]634        else:
635            new_plot.yaxis("\\rm{Intensity} ", "cm^{-1}")
636
[32c0841]637        new_plot.group_id = "Circ avg " + self.data2D.name
638        new_plot.id = "Circ avg " + self.data2D.name
639        new_plot.is_data = True
[13382fc7]640        self.parent.update_theory(data_id=self.data2D, \
641                                       theory=new_plot)
[32c0841]642        wx.PostEvent(self.parent, 
643                     NewPlotEvent(plot=new_plot, title=new_plot.name))
[ef0c170]644       
[1bf33c1]645    def _onEditSlicer(self, event):
[6c0568b]646        """
[d955bf19]647        Is available only when a slicer is drawn.Create a dialog
648        window where the user can enter value to reset slicer
649        parameters.
650       
651        :param event: wx.menu event
652       
[6c0568b]653        """
[32c0841]654        if self.slicer != None:
[1bf33c1]655            from SlicerParameters import SlicerParameterPanel
[4f8a00c]656            dialog = SlicerParameterPanel(self, -1, "Slicer Parameters")
[1bf33c1]657            dialog.set_slicer(self.slicer.__class__.__name__,
658                            self.slicer.get_params())
659            if dialog.ShowModal() == wx.ID_OK:
660                dialog.Destroy() 
661       
662    def onSectorQ(self, event):
663        """
[d955bf19]664        Perform sector averaging on Q and draw sector slicer
[1bf33c1]665        """
[ef0c170]666        from SectorSlicer import SectorInteractor
[1bf33c1]667        self.onClearSlicer(event)
[32c0841]668        wx.PostEvent(self, InternalEvent(slicer=SectorInteractor))
[1bf33c1]669       
670    def onSectorPhi(self, event):
671        """
[d955bf19]672        Perform sector averaging on Phi and draw annulus slicer
[1bf33c1]673        """
[ef0c170]674        from AnnulusSlicer import AnnulusInteractor
[1bf33c1]675        self.onClearSlicer(event)
[32c0841]676        wx.PostEvent(self, InternalEvent(slicer=AnnulusInteractor))
[1bf33c1]677       
[d955bf19]678    def onBoxSum(self, event):
679        """
680        """
[7ab9241]681        from boxSum import BoxSum
682        self.onClearSlicer(event)
[54cc36a]683        self.slicer_z += 1
684        self.slicer =  BoxSum(self, self.subplot, zorder=self.slicer_z)
685        self.subplot.set_ylim(self.data2D.ymin, self.data2D.ymax)
686        self.subplot.set_xlim(self.data2D.xmin, self.data2D.xmax)
687        self.update()
688        self.slicer.update()
[6c0568b]689        ## Value used to initially set the slicer panel
690        type = self.slicer.__class__.__name__
691        params = self.slicer.get_params()
692        ## Create a new panel to display results of summation of Data2D
[8a7a21b]693        from slicerpanel import SlicerPanel
[32c0841]694        new_panel = SlicerPanel(parent=self.parent, id=-1,
695                                    base=self, type=type,
696                                    params=params, style=wx.RAISED_BORDER)
697       
698        new_panel.window_caption = self.slicer.__class__.__name__ + " " + \
699                                    str(self.data2D.name)
700        new_panel.window_name = self.slicer.__class__.__name__+ " " + \
701                                    str(self.data2D.name)
[6c0568b]702        ## Store a reference of the new created panel
[32c0841]703        self.panel_slicer = new_panel
[6c0568b]704        ## save the window_caption of the new panel in the current slicer
[32c0841]705        self.slicer.set_panel_name(name=new_panel.window_caption)
[6c0568b]706        ## post slicer panel to guiframe to display it
[55a0dc1]707        from sans.guiframe.events import SlicerPanelEvent
[32c0841]708        wx.PostEvent(self.parent, SlicerPanelEvent(panel=self.panel_slicer,
709                                                    main_panel=self))
[6c0568b]710
[8a7a21b]711    def onBoxavgX(self,event):
[6c0568b]712        """
[d955bf19]713        Perform 2D data averaging on Qx
714        Create a new slicer .
715       
716        :param event: wx.menu event
[6c0568b]717        """
[8a7a21b]718        from boxSlicer import BoxInteractorX
[38224f10]719        self.onClearSlicer(event)
[32c0841]720        wx.PostEvent(self, InternalEvent(slicer=BoxInteractorX))
[d468daa]721       
[8a7a21b]722    def onBoxavgY(self,event):
[6c0568b]723        """
[d955bf19]724        Perform 2D data averaging on Qy
725        Create a new slicer .
726       
727        :param event: wx.menu event
728       
[6c0568b]729        """
[8a7a21b]730        from boxSlicer import BoxInteractorY
731        self.onClearSlicer(event)
[32c0841]732        wx.PostEvent(self, InternalEvent(slicer=BoxInteractorY))
[6c0568b]733       
[1bf33c1]734    def onClearSlicer(self, event):
735        """
[d955bf19]736        Clear the slicer on the plot
[1bf33c1]737        """
[32c0841]738        if not self.slicer == None:
[1bf33c1]739            self.slicer.clear()
740            self.subplot.figure.canvas.draw()
741            self.slicer = None
742            # Post slicer None event
743            event = self._getEmptySlicerEvent()
[d468daa]744            wx.PostEvent(self, event)
[ed8ad21]745           
746    def _onSave(self, evt):
747        """
748        Save a data set to a dat(text) file
749       
750        :param evt: Menu event
751       
752        """
753        id = str(evt.GetId())
[c553b18]754        if self.parent != None:
755            self._default_save_location = self.parent._default_save_location
[3858e0f]756        default_name = self.plots[self.graph.selected_plottable].label
757        if default_name.count('.') > 0:
758            default_name = default_name.split('.')[0]
759        default_name += "_out"
[ed8ad21]760        if id in self.action_ids:         
761            path = None
[176fbf1]762            self.parent.save_data2d(self.data2D, default_name)
763       
Note: See TracBrowser for help on using the repository browser.