source: sasview/sansview/perspectives/fitting/basepage.py @ 2a62d5c

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 2a62d5c was 1258aa3, checked in by Jae Cho <jhjcho@…>, 13 years ago

disabled bookmark when there is no data set up

  • Property mode set to 100644
File size: 107.2 KB
Line 
1
2import sys
3import os
4import wx
5import numpy
6import time
7import copy 
8import math
9import string
10from sans.guiframe.panel_base import PanelBase
11from wx.lib.scrolledpanel import ScrolledPanel
12from sans.guiframe.utils import format_number,check_float
13from sans.guiframe.events import PanelOnFocusEvent
14from sans.guiframe.events import StatusEvent
15from sans.guiframe.events import AppendBookmarkEvent
16import pagestate
17from pagestate import PageState
18(PageInfoEvent, EVT_PAGE_INFO)   = wx.lib.newevent.NewEvent()
19(PreviousStateEvent, EVT_PREVIOUS_STATE)   = wx.lib.newevent.NewEvent()
20(NextStateEvent, EVT_NEXT_STATE)   = wx.lib.newevent.NewEvent()
21
22_BOX_WIDTH = 76
23_QMIN_DEFAULT = 0.0005
24_QMAX_DEFAULT = 0.5
25_NPTS_DEFAULT = 200
26#Control panel width
27if sys.platform.count("darwin")==0:
28    PANEL_WIDTH = 450
29    FONT_VARIANT = 0
30    ON_MAC = False
31else:
32    PANEL_WIDTH = 500
33    FONT_VARIANT = 1
34    ON_MAC = True
35
36
37
38class BasicPage(ScrolledPanel, PanelBase):
39    """
40    This class provide general structure of  fitpanel page
41    """
42     ## Internal name for the AUI manager
43    window_name = "Fit Page"
44    ## Title to appear on top of the window
45    window_caption = "Fit Page "
46   
47    def __init__(self, parent,color='blue', **kwargs):
48        """
49        """
50        ScrolledPanel.__init__(self, parent, **kwargs)
51        PanelBase.__init__(self, parent)
52        self.SetupScrolling()
53        #Set window's font size
54        self.SetWindowVariant(variant=FONT_VARIANT)
55     
56        self.SetBackgroundColour(color)
57        ## parent of the page
58        self.parent = parent
59        ## manager is the fitting plugin
60        ## owner of the page (fitting plugin)
61        self.event_owner = None
62         ## current model
63        self.model = None
64        ## data
65        self.data = None
66        self.mask = None
67        self.uid = None
68        ## Q range
69        self.qmin = None
70        self.qmax = None
71        self.qmax_x = _QMAX_DEFAULT
72        self.qmin_x = _QMIN_DEFAULT
73        self.npts_x = _NPTS_DEFAULT
74        ## total number of point: float
75        self.npts = None
76        ## default fitengine type
77        self.engine_type = 'scipy'
78        ## smear default
79        self.smearer = None
80        self.current_smearer = None
81        ## 2D smear accuracy default
82        self.smear2d_accuracy = 'Low'
83        ## slit smear:
84        self.dxl = None
85        self.dxw = None
86        ## pinhole smear
87        self.dx_min = None
88        self.dx_max = None
89       
90        self.disp_cb_dict = {}
91   
92        self.state = PageState(parent=parent)
93        ## dictionary containing list of models
94        self.model_list_box = {}
95       
96        ## Data member to store the dispersion object created
97        self._disp_obj_dict = {}
98        ## selected parameters to apply dispersion
99        self.disp_cb_dict ={}
100
101        ## smearer object
102        self.smearer = None
103        self.enable2D = False
104       
105        ##list of model parameters. each item must have same length
106        ## each item related to a given parameters
107        ##[cb state, name, value, "+/-", error of fit, min, max , units]
108        self.parameters = []
109        # non-fittable parameter whose value is astring
110        self.str_parameters = []
111        ## list of parameters to fit , must be like self.parameters
112        self.param_toFit = []
113        ## list of looking like parameters but with non fittable parameters info
114        self.fixed_param = []
115        ## list of looking like parameters but with  fittable parameters info
116        self.fittable_param = []
117        ##list of dispersion parameters
118        self.disp_list = []
119        self.disp_name = ""
120       
121        ## list of orientation parameters
122        self.orientation_params = []
123        self.orientation_params_disp = []
124        if self.model != None:
125            self.disp_list = self.model.getDispParamList()
126       
127        ##enable model 2D draw
128        self.enable2D = False
129        ## check that the fit range is correct to plot the model again
130        self.fitrange = True
131        ## Create memento to save the current state
132        self.state = PageState(parent=self.parent,
133                               model=self.model, data=self.data)
134        ## flag to determine if state has change
135        self.state_change = False
136        ## save customized array
137        self.values = []
138        self.weights = []
139        ## retrieve saved state
140        self.number_saved_state = 0
141        ## dictionary of saved state
142        self.saved_states = {} 
143        ## Create context menu for page
144        self.popUpMenu = wx.Menu()
145   
146        id = wx.NewId()
147        self._keep = wx.MenuItem(self.popUpMenu,id,"BookMark",
148                                 " Keep the panel status to recall it later")
149        self.popUpMenu.AppendItem(self._keep)
150        self._keep.Enable(True)
151        self._set_bookmark_flag(False)
152        self._set_save_flag(False)
153        wx.EVT_MENU(self, id, self.on_bookmark)
154        self.popUpMenu.AppendSeparator()
155   
156        ## Default locations
157        self._default_save_location = os.getcwd()     
158        ## save initial state on context menu
159        #self.onSave(event=None)
160        self.Bind(wx.EVT_CONTEXT_MENU, self.onContextMenu)
161       
162        # bind key event
163        self.Bind(wx.EVT_LEFT_DOWN, self.on_left_down)
164
165        ## create the basic structure of the panel with empty sizer
166        self.define_page_structure()
167        ## drawing Initial dispersion parameters sizer
168        self.set_dispers_sizer()
169       
170        ## layout
171        self.set_layout()
172       
173   
174       
175    def on_set_focus(self, event):
176        """
177        """
178        if self._manager is not None:
179            wx.PostEvent(self._manager.parent, PanelOnFocusEvent(panel=self))
180       
181    class ModelTextCtrl(wx.TextCtrl):
182        """
183        Text control for model and fit parameters.
184        Binds the appropriate events for user interactions.
185        Default callback methods can be overwritten on initialization
186       
187        :param kill_focus_callback: callback method for EVT_KILL_FOCUS event
188        :param set_focus_callback:  callback method for EVT_SET_FOCUS event
189        :param mouse_up_callback:   callback method for EVT_LEFT_UP event
190        :param text_enter_callback: callback method for EVT_TEXT_ENTER event
191       
192        """
193        ## Set to True when the mouse is clicked while the whole string is selected
194        full_selection = False
195        ## Call back for EVT_SET_FOCUS events
196        _on_set_focus_callback = None
197       
198        def __init__(self, parent, id=-1, 
199                     value=wx.EmptyString, 
200                     pos=wx.DefaultPosition, 
201                     size=wx.DefaultSize,
202                     style=0, 
203                     validator=wx.DefaultValidator,
204                     name=wx.TextCtrlNameStr,
205                     kill_focus_callback=None,
206                     set_focus_callback=None,
207                     mouse_up_callback=None,
208                     text_enter_callback = None):
209             
210            wx.TextCtrl.__init__(self, parent, id, value, pos,
211                                  size, style, validator, name)
212           
213            # Bind appropriate events
214            self._on_set_focus_callback = parent.onSetFocus \
215                      if set_focus_callback is None else set_focus_callback
216            self.Bind(wx.EVT_SET_FOCUS, self._on_set_focus)
217            self.Bind(wx.EVT_KILL_FOCUS, self._silent_kill_focus \
218                      if kill_focus_callback is None else kill_focus_callback)               
219            self.Bind(wx.EVT_TEXT_ENTER, parent._onparamEnter \
220                      if text_enter_callback is None else text_enter_callback)
221            if not ON_MAC :
222                self.Bind(wx.EVT_LEFT_UP,    self._highlight_text \
223                          if mouse_up_callback is None else mouse_up_callback)
224           
225        def _on_set_focus(self, event):
226            """
227            Catch when the text control is set in focus to highlight the whole
228            text if necessary
229           
230            :param event: mouse event
231           
232            """
233            event.Skip()
234            self.full_selection = True
235            return self._on_set_focus_callback(event)
236       
237 
238           
239        def _highlight_text(self, event):
240            """
241            Highlight text of a TextCtrl only of no text has be selected
242           
243            :param event: mouse event
244           
245            """
246            # Make sure the mouse event is available to other listeners
247            event.Skip()
248            control  = event.GetEventObject()
249            if self.full_selection:
250                self.full_selection = False
251                # Check that we have a TextCtrl
252                if issubclass(control.__class__, wx.TextCtrl):
253                    # Check whether text has been selected,
254                    # if not, select the whole string
255                    (start, end) = control.GetSelection()
256                    if start==end:
257                        control.SetSelection(-1,-1)
258                       
259        def _silent_kill_focus(self,event):
260            """
261            Save the state of the page
262            """
263           
264            event.Skip()
265            pass
266   
267    def set_page_info(self, page_info):
268        """
269        set some page important information at once
270        """
271       ##window_name
272        self.window_name = page_info.window_name
273        ##window_caption
274        self.window_caption = page_info.window_caption
275        ## manager is the fitting plugin
276        self._manager= page_info.manager
277        ## owner of the page (fitting plugin)
278        self.event_owner= page_info.event_owner
279         ## current model
280        self.model = page_info.model
281        ## data
282        self.data = page_info.data
283        ## dictionary containing list of models
284        self.model_list_box = page_info.model_list_box
285        ## Data member to store the dispersion object created
286        self.populate_box(dict=self.model_list_box)
287       
288    def onContextMenu(self, event): 
289        """
290        Retrieve the state selected state
291        """
292        # Skipping the save state functionality for release 0.9.0
293        #return
294   
295        pos = event.GetPosition()
296        pos = self.ScreenToClient(pos)
297       
298        self.PopupMenu(self.popUpMenu, pos) 
299     
300       
301    def onUndo(self, event):
302        """
303        Cancel the previous action
304        """
305        event = PreviousStateEvent(page = self)
306        wx.PostEvent(self.parent, event)
307       
308    def onRedo(self, event):
309        """
310        Restore the previous action cancelled
311        """
312        event = NextStateEvent(page= self)
313        wx.PostEvent(self.parent, event)
314   
315    def define_page_structure(self):
316        """
317        Create empty sizer for a panel
318        """
319        self.vbox  = wx.BoxSizer(wx.VERTICAL)
320        self.sizer0 = wx.BoxSizer(wx.VERTICAL)
321        self.sizer1 = wx.BoxSizer(wx.VERTICAL)
322        self.sizer2 = wx.BoxSizer(wx.VERTICAL)
323        self.sizer3 = wx.BoxSizer(wx.VERTICAL)
324        self.sizer4 = wx.BoxSizer(wx.VERTICAL)
325        self.sizer5 = wx.BoxSizer(wx.VERTICAL)
326        self.sizer6 = wx.BoxSizer(wx.VERTICAL)
327       
328        self.sizer0.SetMinSize((PANEL_WIDTH,-1))
329        self.sizer1.SetMinSize((PANEL_WIDTH,-1))
330        self.sizer2.SetMinSize((PANEL_WIDTH,-1))
331        self.sizer3.SetMinSize((PANEL_WIDTH,-1))
332        self.sizer4.SetMinSize((PANEL_WIDTH,-1))
333        self.sizer5.SetMinSize((PANEL_WIDTH,-1))
334        self.sizer6.SetMinSize((PANEL_WIDTH,-1))
335       
336        self.vbox.Add(self.sizer0)
337        self.vbox.Add(self.sizer1)
338        self.vbox.Add(self.sizer2)
339        self.vbox.Add(self.sizer3)
340        self.vbox.Add(self.sizer4)
341        self.vbox.Add(self.sizer5)
342        self.vbox.Add(self.sizer6)
343       
344    def set_layout(self):
345        """
346        layout
347        """
348        self.vbox.Layout()
349        self.vbox.Fit(self) 
350        self.SetSizer(self.vbox)
351        self.Centre()
352 
353    def set_owner(self,owner):
354        """
355        set owner of fitpage
356       
357        :param owner: the class responsible of plotting
358       
359        """
360        self.event_owner = owner   
361        self.state.event_owner = owner
362       
363    def get_state(self):
364        """
365        """
366        return self.state
367    def get_data(self):
368        """
369        return the current data
370        """
371        return self.data 
372   
373    def set_manager(self, manager):
374        """
375        set panel manager
376       
377        :param manager: instance of plugin fitting
378       
379        """
380        self._manager = manager 
381        self.state.manager = manager
382       
383    def populate_box(self, dict):
384        """
385        Store list of model
386       
387        :param dict: dictionary containing list of models
388       
389        """
390        self.model_list_box = dict
391        self.state.model_list_box = self.model_list_box
392        self.initialize_combox()
393       
394    def initialize_combox(self): 
395        """
396        put default value in the combobox
397        """ 
398        ## fill combox box
399        if self.model_list_box is None:
400            return
401        if len(self.model_list_box) > 0:
402            self._populate_box(self.formfactorbox,
403                               self.model_list_box["Shapes"])
404       
405        if len(self.model_list_box) > 0:
406            self._populate_box(self.structurebox,
407                                self.model_list_box["Structure Factors"])
408            self.structurebox.Insert("None", 0, None)
409            self.structurebox.SetSelection(0)
410            self.structurebox.Hide()
411            self.text2.Hide()
412            self.structurebox.Disable()
413            self.text2.Disable()
414             
415            if self.model.__class__ in self.model_list_box["P(Q)*S(Q)"]:
416                self.structurebox.Show()
417                self.text2.Show()
418                self.structurebox.Enable()
419                self.text2.Enable()           
420               
421    def set_dispers_sizer(self):
422        """
423        fill sizer containing dispersity info
424        """
425        self.sizer4.Clear(True)
426        name="Polydispersity and Orientational Distribution"
427        box_description= wx.StaticBox(self, -1,name)
428        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
429        #----------------------------------------------------
430        self.disable_disp = wx.RadioButton(self, -1, 'Off', (10, 10),
431                                            style=wx.RB_GROUP)
432        self.enable_disp = wx.RadioButton(self, -1, 'On', (10, 30))
433
434        self.disp_help_bt = wx.Button(self,wx.NewId(),'?', size=(20,15))
435        self.disp_help_bt.Bind(wx.EVT_BUTTON, 
436                        self.on_pd_help_clicked,id= self.disp_help_bt.GetId())
437        self.disp_help_bt.SetToolTipString("Helps for Polydispersion.")       
438       
439        self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param,
440                     id=self.disable_disp.GetId())
441        self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param,
442                   id=self.enable_disp.GetId())
443        #MAC needs SetValue
444        self.disable_disp.SetValue(True)
445        sizer_dispersion = wx.BoxSizer(wx.HORIZONTAL)
446        sizer_dispersion.Add((20,20))
447        name=""#Polydispersity and \nOrientational Distribution "
448        sizer_dispersion.Add(wx.StaticText(self,-1,name))
449        sizer_dispersion.Add(self.enable_disp )
450        sizer_dispersion.Add((20,20))
451        sizer_dispersion.Add(self.disable_disp )
452        sizer_dispersion.Add((25,20))
453        sizer_dispersion.Add(self.disp_help_bt)
454       
455        ## fill a sizer for dispersion         
456        boxsizer1.Add( sizer_dispersion,0,
457                wx.TOP|wx.BOTTOM|wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE,border=5)
458        #boxsizer1.Add( (10,10) )
459        #boxsizer1.Add( sizer_select_dispers )
460        self.sizer4_4 = wx.GridBagSizer(6,5)
461
462        boxsizer1.Add( self.sizer4_4  )
463        #-----------------------------------------------------
464        self.sizer4.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
465        self.sizer4_4.Layout()
466        self.sizer4.Layout()
467        self.Layout()
468     
469        self.Refresh()
470        ## saving the state of enable dispersity button
471        self.state.enable_disp= self.enable_disp.GetValue()
472        self.state.disable_disp= self.disable_disp.GetValue()
473        self.SetupScrolling()
474
475   
476    def onResetModel(self, event):
477        """
478        Reset model state
479        """
480        menu = event.GetEventObject()
481        ## post help message for the selected model
482        msg = menu.GetHelpString(event.GetId())
483        msg +=" reloaded"
484        wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
485       
486        name = menu.GetLabel(event.GetId())
487        self._on_select_model_helper()
488       
489        if name in self.saved_states.keys():
490            previous_state = self.saved_states[name]
491            ## reset state of checkbox,textcrtl  and  regular parameters value
492            self.reset_page(previous_state)   
493               
494    def on_preview(self, event):
495        """
496        Report the current fit results
497        """   
498        # Get plot image from plotpanel
499        images, canvases = self.get_images()
500        # get the report dialog
501        self.state.report(images, canvases)
502       
503         
504    def on_save(self, event):   
505        """
506        Save the current state into file
507        """ 
508        self.save_current_state()
509        new_state = self.state.clone()
510        # Ask the user the location of the file to write to.
511        path = None
512        dlg = wx.FileDialog(self, "Choose a file", self._default_save_location,
513                                                 "", "*.fitv", wx.SAVE)
514        if dlg.ShowModal() == wx.ID_OK:
515            path = dlg.GetPath()
516            self._default_save_location = os.path.dirname(path)
517        else:
518            return None
519        #the manager write the state into file
520        self._manager.save_fit_state(filepath=path, fitstate=new_state)
521        return new_state 
522   
523    def _get_time_stamp(self):
524        """
525        return time and date stings
526        """
527        # date and time
528        year, month, day,hour,minute,second,tda,ty,tm_isdst= time.localtime()
529        current_time= str(hour)+":"+str(minute)+":"+str(second)
530        current_date= str( month)+"/"+str(day)+"/"+str(year)
531        return current_time, current_date
532     
533    def on_bookmark(self, event):
534        """
535        save history of the data and model
536        """
537        if self.model==None:
538            msg="Can not bookmark; Please select Data and Model first..."
539            wx.MessageBox(msg, 'Info')
540            return 
541        self.save_current_state()
542        new_state = self.state.clone()
543        ##Add model state on context menu
544        self.number_saved_state += 1
545        current_time, current_date = self._get_time_stamp()
546        #name= self.model.name+"[%g]"%self.number_saved_state
547        name = "Fitting: %g]" % self.number_saved_state
548        name += self.model.__class__.__name__
549        name += "bookmarked at %s on %s" % (current_time, current_date)
550        self.saved_states[name]= new_state
551       
552        ## Add item in the context menu
553        msg =  "Model saved at %s on %s"%(current_time, current_date)
554         ## post help message for the selected model
555        msg +=" Saved! right click on this page to retrieve this model"
556        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
557       
558        id = wx.NewId()
559        self.popUpMenu.Append(id,name,str(msg))
560        wx.EVT_MENU(self, id, self.onResetModel)
561        wx.PostEvent(self.parent.parent, 
562                     AppendBookmarkEvent(title=name, 
563                                         hint=str(msg), handler=self.onResetModel))
564       
565    def old_on_bookmark(self, event):
566        """
567        save history of the data and model
568        """
569        if self.model==None:
570            msg="Can not bookmark; Please select Data and Model first..."
571            wx.MessageBox(msg, 'Info')
572            return 
573        if hasattr(self,"enable_disp"):
574            self.state.enable_disp = copy.deepcopy(self.enable_disp.GetValue())
575        if hasattr(self, "disp_box"):
576            self.state.disp_box = copy.deepcopy(self.disp_box.GetSelection())
577
578        self.state.model.name= self.model.name
579       
580        #Remember fit engine_type for fit panel
581        if self.engine_type == None: 
582            self.engine_type = "scipy"
583        if self._manager !=None:
584            self._manager._on_change_engine(engine=self.engine_type)
585       
586            self.state.engine_type = self.engine_type
587
588        new_state = self.state.clone()
589        new_state.model.name = self.state.model.name
590       
591        new_state.enable2D = copy.deepcopy(self.enable2D)
592        ##Add model state on context menu
593        self.number_saved_state += 1
594        #name= self.model.name+"[%g]"%self.number_saved_state
595        name= self.model.__class__.__name__+"[%g]"%self.number_saved_state
596        self.saved_states[name]= new_state
597       
598        ## Add item in the context menu
599       
600        year, month, day,hour,minute,second,tda,ty,tm_isdst= time.localtime()
601        my_time= str(hour)+" : "+str(minute)+" : "+str(second)+" "
602        date= str( month)+"|"+str(day)+"|"+str(year)
603        msg=  "Model saved at %s on %s"%(my_time, date)
604         ## post help message for the selected model
605        msg +=" Saved! right click on this page to retrieve this model"
606        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
607       
608        id = wx.NewId()
609        self.popUpMenu.Append(id,name,str(msg))
610        wx.EVT_MENU(self, id, self.onResetModel)
611       
612    def onSetFocus(self, evt):
613        """
614        highlight the current textcrtl and hide the error text control shown
615        after fitting
616        """
617        return
618   
619    def read_file(self, path):
620        """
621        Read two columns file
622       
623        :param path: the path to the file to read
624       
625        """
626        try:
627            if path==None:
628                wx.PostEvent(self.parent.parent, StatusEvent(status=\
629                            " Selected Distribution was not loaded: %s"%path))
630                return None, None
631            input_f = open(path, 'r')
632            buff = input_f.read()
633            lines = buff.split('\n')
634           
635            angles = []
636            weights=[]
637            for line in lines:
638                toks = line.split()
639                try:
640                    angle = float(toks[0])
641                    weight = float(toks[1])
642                except:
643                    # Skip non-data lines
644                    pass
645                angles.append(angle)
646                weights.append(weight)
647            return numpy.array(angles), numpy.array(weights)
648        except:
649            raise 
650
651    def createMemento(self):
652        """
653        return the current state of the page
654        """
655        return self.state.clone()
656   
657   
658    def save_current_state(self):
659        """
660        Store current state
661        """
662        self.state.engine_type = copy.deepcopy(self.engine_type)
663        ## save model option
664        if self.model!= None:
665            self.disp_list= self.model.getDispParamList()
666            self.state.disp_list= copy.deepcopy(self.disp_list)
667            self.state.model = self.model.clone()
668        #save radiobutton state for model selection
669        self.state.shape_rbutton = self.shape_rbutton.GetValue()
670        self.state.shape_indep_rbutton = self.shape_indep_rbutton.GetValue()
671        self.state.struct_rbutton = self.struct_rbutton.GetValue()
672        self.state.plugin_rbutton = self.plugin_rbutton.GetValue()
673        #model combobox
674        self.state.structurebox = self.structurebox.GetSelection()
675        self.state.formfactorbox = self.formfactorbox.GetSelection()
676       
677        self.state.enable2D = copy.deepcopy(self.enable2D)
678        self.state.values= copy.deepcopy(self.values)
679        self.state.weights = copy.deepcopy( self.weights)
680        ## save data   
681        self.state.data= copy.deepcopy(self.data)
682        self.state.qmax_x = self.qmax_x
683        self.state.qmin_x = self.qmin_x
684     
685        if hasattr(self,"enable_disp"):
686            self.state.enable_disp= self.enable_disp.GetValue()
687            self.state.disable_disp = self.disable_disp.GetValue()
688           
689        self.state.smearer = copy.deepcopy(self.smearer)
690        if hasattr(self,"enable_smearer"):
691            self.state.enable_smearer = \
692                                copy.deepcopy(self.enable_smearer.GetValue())
693            self.state.disable_smearer = \
694                                copy.deepcopy(self.disable_smearer.GetValue())
695
696        self.state.pinhole_smearer = \
697                                copy.deepcopy(self.pinhole_smearer.GetValue())
698        self.state.slit_smearer = copy.deepcopy(self.slit_smearer.GetValue()) 
699                 
700        if len(self._disp_obj_dict)>0:
701            for k , v in self._disp_obj_dict.iteritems():
702                self.state._disp_obj_dict[k]= v
703                       
704           
705            self.state.values = copy.deepcopy(self.values)
706            self.state.weights = copy.deepcopy(self.weights)
707        ## save plotting range
708        self._save_plotting_range()
709       
710        self.state.orientation_params = []
711        self.state.orientation_params_disp = []
712        self.state.parameters = []
713        self.state.fittable_param = []
714        self.state.fixed_param = []
715        self.state.str_parameters = []
716
717       
718        ## save checkbutton state and txtcrtl values
719        self._copy_parameters_state(self.str_parameters, 
720                                    self.state.str_parameters)
721        self._copy_parameters_state(self.orientation_params,
722                                     self.state.orientation_params)
723        self._copy_parameters_state(self.orientation_params_disp,
724                                     self.state.orientation_params_disp)
725       
726        self._copy_parameters_state(self.parameters, self.state.parameters)
727        self._copy_parameters_state(self.fittable_param,
728                                     self.state.fittable_param)
729        self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
730        #save chisqr
731        self.state.tcChi = self.tcChi.GetValue()
732       
733    def save_current_state_fit(self):
734        """
735        Store current state for fit_page
736        """
737        ## save model option
738        if self.model!= None:
739            self.disp_list= self.model.getDispParamList()
740            self.state.disp_list= copy.deepcopy(self.disp_list)
741            self.state.model = self.model.clone()
742        if hasattr(self, "engine_type"):
743            self.state.engine_type = copy.deepcopy(self.engine_type)
744           
745        self.state.enable2D = copy.deepcopy(self.enable2D)
746        self.state.values= copy.deepcopy(self.values)
747        self.state.weights = copy.deepcopy( self.weights)
748        ## save data   
749        self.state.data= copy.deepcopy(self.data)
750       
751        if hasattr(self,"enable_disp"):
752            self.state.enable_disp= self.enable_disp.GetValue()
753            self.state.disable_disp = self.disable_disp.GetValue()
754           
755        self.state.smearer = copy.deepcopy(self.smearer)
756        if hasattr(self,"enable_smearer"):
757            self.state.enable_smearer = \
758                                copy.deepcopy(self.enable_smearer.GetValue())
759            self.state.disable_smearer = \
760                                copy.deepcopy(self.disable_smearer.GetValue())
761           
762        self.state.pinhole_smearer = \
763                                copy.deepcopy(self.pinhole_smearer.GetValue())
764        self.state.slit_smearer = copy.deepcopy(self.slit_smearer.GetValue()) 
765           
766        if hasattr(self,"disp_box"):
767            self.state.disp_box = self.disp_box.GetCurrentSelection()
768
769            if len(self.disp_cb_dict) > 0:
770                for k, v in self.disp_cb_dict.iteritems():
771         
772                    if v == None :
773                        self.state.disp_cb_dict[k] = v
774                    else:
775                        try:
776                            self.state.disp_cb_dict[k] = v.GetValue()
777                        except:
778                            self.state.disp_cb_dict[k] = None
779           
780            if len(self._disp_obj_dict) > 0:
781                for k , v in self._disp_obj_dict.iteritems():
782     
783                    self.state._disp_obj_dict[k] = v
784                       
785           
786            self.state.values = copy.deepcopy(self.values)
787            self.state.weights = copy.deepcopy(self.weights)
788           
789        ## save plotting range
790        self._save_plotting_range()
791       
792        ## save checkbutton state and txtcrtl values
793        self._copy_parameters_state(self.orientation_params,
794                                     self.state.orientation_params)
795        self._copy_parameters_state(self.orientation_params_disp,
796                                     self.state.orientation_params_disp)
797        self._copy_parameters_state(self.parameters, self.state.parameters)
798        self._copy_parameters_state(self.fittable_param,
799                                             self.state.fittable_param)
800        self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
801   
802         
803    def check_invalid_panel(self): 
804        """
805        check if the user can already perform some action with this panel
806        """ 
807        flag = False
808        if self.data is None:
809            self.disable_smearer.SetValue(True)
810            self.disable_disp.SetValue(True)
811            msg = "Please load Data and select Model to start..."
812            wx.MessageBox(msg, 'Info')
813            return  True
814       
815    def set_model_state(self, state):
816        """
817        reset page given a model state
818        """
819        self.disp_cb_dict = state.disp_cb_dict
820        self.disp_list = state.disp_list
821     
822        ## set the state of the radio box
823        self.shape_rbutton.SetValue(state.shape_rbutton )
824        self.shape_indep_rbutton.SetValue(state.shape_indep_rbutton)
825        self.struct_rbutton.SetValue(state.struct_rbutton)
826        self.plugin_rbutton.SetValue(state.plugin_rbutton)
827       
828        ## fill model combobox
829        self._show_combox_helper()
830        #select the current model
831        self.formfactorbox.Select(int(state.formfactorcombobox))
832        self.structurebox.SetSelection(state.structurecombobox )
833        if state.multi_factor != None:
834            self.multifactorbox.SetSelection(state.multi_factor)
835           
836         ## reset state of checkbox,textcrtl  and  regular parameters value
837        self._reset_parameters_state(self.orientation_params_disp,
838                                     state.orientation_params_disp)
839        self._reset_parameters_state(self.orientation_params,
840                                     state.orientation_params)
841        self._reset_parameters_state(self.str_parameters,
842                                     state.str_parameters)
843        self._reset_parameters_state(self.parameters,state.parameters)
844         ## display dispersion info layer       
845        self.enable_disp.SetValue(state.enable_disp)
846        self.disable_disp.SetValue(state.disable_disp)
847       
848        if hasattr(self, "disp_box"):
849           
850            self.disp_box.SetSelection(state.disp_box) 
851            n= self.disp_box.GetCurrentSelection()
852            dispersity= self.disp_box.GetClientData(n)
853            name = dispersity.__name__     
854
855            self._set_dipers_Param(event=None)
856       
857            if name == "ArrayDispersion":
858               
859                for item in self.disp_cb_dict.keys():
860                   
861                    if hasattr(self.disp_cb_dict[item], "SetValue") :
862                        self.disp_cb_dict[item].SetValue(\
863                                                    state.disp_cb_dict[item])
864                        # Create the dispersion objects
865                        from sans.models.dispersion_models import ArrayDispersion
866                        disp_model = ArrayDispersion()
867                        if hasattr(state,"values")and\
868                                 self.disp_cb_dict[item].GetValue() == True:
869                            if len(state.values)>0:
870                                self.values=state.values
871                                self.weights=state.weights
872                                disp_model.set_weights(self.values,
873                                                        state.weights)
874                            else:
875                                self._reset_dispersity()
876                       
877                        self._disp_obj_dict[item] = disp_model
878                        # Set the new model as the dispersion object
879                        #for the selected parameter
880                        self.model.set_dispersion(item, disp_model)
881                   
882                        self.model._persistency_dict[item] = \
883                                                [state.values, state.weights]
884                   
885            else:
886                keys = self.model.getParamList()
887                for item in keys:
888                    if item in self.disp_list and \
889                        not self.model.details.has_key(item):
890                        self.model.details[item] = ["", None, None]
891                for k,v in self.state.disp_cb_dict.iteritems():
892                    self.disp_cb_dict = copy.deepcopy(state.disp_cb_dict) 
893                    self.state.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
894         ## smearing info  restore
895        if hasattr(self, "enable_smearer"):
896            ## set smearing value whether or not the data
897            #contain the smearing info
898            self.enable_smearer.SetValue(state.enable_smearer)
899            self.disable_smearer.SetValue(state.disable_smearer)
900            self.onSmear(event=None)           
901        self.pinhole_smearer.SetValue(state.pinhole_smearer)
902        self.slit_smearer.SetValue(state.slit_smearer)
903        ## we have two more options for smearing
904        if self.pinhole_smearer.GetValue(): self.onPinholeSmear(event=None)
905        elif self.slit_smearer.GetValue(): self.onSlitSmear(event=None)
906       
907        ## reset state of checkbox,textcrtl  and dispersity parameters value
908        self._reset_parameters_state(self.fittable_param,state.fittable_param)
909        self._reset_parameters_state(self.fixed_param,state.fixed_param)
910       
911        ## draw the model with previous parameters value
912        self._onparamEnter_helper()
913        self.select_param(event=None) 
914        #Save state_fit
915        self.save_current_state_fit()
916        self._lay_out()
917        self.Refresh()
918       
919    def reset_page_helper(self, state):
920        """
921        Use page_state and change the state of existing page
922       
923        :precondition: the page is already drawn or created
924       
925        :postcondition: the state of the underlying data change as well as the
926            state of the graphic interface
927        """
928        if state == None:
929            #self._undo.Enable(False)
930            return 
931        # set data, etc. from the state
932        self.set_data(state.data)
933        self.enable2D= state.enable2D
934        self.engine_type = state.engine_type
935
936        self.disp_cb_dict = state.disp_cb_dict
937        self.disp_list = state.disp_list
938     
939        ## set the state of the radio box
940        self.shape_rbutton.SetValue(state.shape_rbutton )
941        self.shape_indep_rbutton.SetValue(state.shape_indep_rbutton)
942        self.struct_rbutton.SetValue(state.struct_rbutton)
943        self.plugin_rbutton.SetValue(state.plugin_rbutton)
944       
945        ## fill model combobox
946        self._show_combox_helper()
947        #select the current model
948        self.formfactorbox.Select(int(state.formfactorcombobox))
949        self.structurebox.SetSelection(state.structurecombobox )
950        if state.multi_factor != None:
951            self.multifactorbox.SetSelection(state.multi_factor)
952
953        #reset the fitting engine type
954        self.engine_type = state.engine_type
955        #draw the pnael according to the new model parameter
956        self._on_select_model(event=None)
957       
958        if self._manager !=None:
959            self._manager._on_change_engine(engine=self.engine_type)
960        ## set the select all check box to the a given state
961        self.cb1.SetValue(state.cb1)
962     
963        ## reset state of checkbox,textcrtl  and  regular parameters value
964        self._reset_parameters_state(self.orientation_params_disp,
965                                     state.orientation_params_disp)
966        self._reset_parameters_state(self.orientation_params,
967                                     state.orientation_params)
968        self._reset_parameters_state(self.str_parameters,
969                                     state.str_parameters)
970        self._reset_parameters_state(self.parameters,state.parameters)   
971         ## display dispersion info layer       
972        self.enable_disp.SetValue(state.enable_disp)
973        self.disable_disp.SetValue(state.disable_disp)
974        # If the polydispersion is ON
975        if state.enable_disp:
976            # reset dispersion according the state
977            self._set_dipers_Param(event=None)
978            self._reset_page_disp_helper(state)
979        ##plotting range restore   
980        self._reset_plotting_range(state)
981        ## smearing info  restore
982        if hasattr(self, "enable_smearer"):
983            ## set smearing value whether or not the data
984            #contain the smearing info
985            self.enable_smearer.SetValue(state.enable_smearer)
986            self.disable_smearer.SetValue(state.disable_smearer)
987            self.onSmear(event=None)           
988        self.pinhole_smearer.SetValue(state.pinhole_smearer)
989        self.slit_smearer.SetValue(state.slit_smearer)
990        ## we have two more options for smearing
991        if self.pinhole_smearer.GetValue(): self.onPinholeSmear(event=None)
992        elif self.slit_smearer.GetValue(): self.onSlitSmear(event=None)
993       
994        ## reset state of checkbox,textcrtl  and dispersity parameters value
995        self._reset_parameters_state(self.fittable_param,state.fittable_param)
996        self._reset_parameters_state(self.fixed_param,state.fixed_param)
997       
998        ## draw the model with previous parameters value
999        self._onparamEnter_helper()
1000        #reset the value of chisqr when not consistent with the value computed
1001        self.tcChi.SetValue(str(self.state.tcChi))
1002        ## reset context menu items
1003        self._reset_context_menu()
1004       
1005        ## set the value of the current state to the state given as parameter
1006        self.state = state.clone() 
1007   
1008    def _reset_page_disp_helper(self, state):
1009        """
1010        Help to rest page for dispersions
1011        """
1012        keys = self.model.getParamList()
1013        for item in keys:
1014            if item in self.disp_list and \
1015                not self.model.details.has_key(item):
1016                self.model.details[item] = ["", None, None]
1017        #for k,v in self.state.disp_cb_dict.iteritems():
1018        self.disp_cb_dict = copy.deepcopy(state.disp_cb_dict) 
1019        self.state.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1020        self.values = copy.deepcopy(state.values)
1021        self.weights = copy.deepcopy(state.weights)
1022       
1023        for key, disp in state._disp_obj_dict.iteritems():
1024            # From saved file, disp_model can not be sent in model obj.
1025            # it will be sent as a string here, then converted to model object.
1026            if disp.__class__.__name__ == 'str':
1027                com_str  = "from sans.models.dispersion_models "
1028                com_str += "import %s as disp_func"
1029                exec com_str % disp
1030                disp_model = disp_func()
1031            else:
1032                disp_model = disp
1033
1034            self._disp_obj_dict[key] = disp_model
1035            param_name = key.split('.')[0]
1036            # Try to set dispersion only when available
1037            # for eg., pass the orient. angles for 1D Cal
1038            try:
1039                self.model.set_dispersion(param_name, disp_model)
1040                self.model._persistency_dict[key] = \
1041                                 [state.values, state.weights]
1042            except:
1043                pass
1044            selection = self._find_polyfunc_selection(disp_model)
1045            for list in self.fittable_param:
1046                if list[1] == key and list[7] != None:
1047                    list[7].SetSelection(selection)
1048                    # For the array disp_model, set the values and weights
1049                    if selection == 1:
1050                        disp_model.set_weights(self.values[key], 
1051                                              self.weights[key])
1052                        try:
1053                            # Diables all fittable params for array
1054                            list[0].SetValue(False)
1055                            list[0].Disable()
1056                            list[2].Disable()
1057                            list[5].Disable()
1058                            list[6].Disable()
1059                        except:
1060                            pass
1061            # For array, disable all fixed params
1062            if selection == 1:
1063                for item in self.fixed_param:
1064                    if item[1].split(".")[0] == key.split(".")[0]:
1065                        # try it and pass it for the orientation for 1D
1066                        try:
1067                            item[2].Disable()
1068                        except:
1069                            pass
1070   
1071        # Make sure the check box updated when all checked
1072        if self.cb1.GetValue():
1073            self.select_all_param(None)       
1074     
1075    def _selectDlg(self):
1076        """
1077        open a dialog file to selected the customized dispersity
1078        """
1079        import os
1080        dlg = wx.FileDialog(self, "Choose a weight file",
1081                                self._default_save_location , "", 
1082                                "*.*", wx.OPEN)
1083        path = None
1084        if dlg.ShowModal() == wx.ID_OK:
1085            path = dlg.GetPath()
1086        dlg.Destroy()
1087        return path
1088
1089    def _reset_context_menu(self):
1090        """
1091        reset the context menu
1092        """
1093        for name, state in self.state.saved_states.iteritems():
1094            self.number_saved_state += 1
1095            ## Add item in the context menu
1096            id = wx.NewId()
1097            msg = 'Save model and state %g' % self.number_saved_state
1098            self.popUpMenu.Append(id, name, msg)
1099            wx.EVT_MENU(self, id, self.onResetModel)
1100   
1101    def _reset_plotting_range(self, state):
1102        """
1103        Reset the plotting range to a given state
1104        """
1105        # if self.check_invalid_panel():
1106        #    return
1107        self.qmin.SetValue(str(state.qmin))
1108        self.qmax.SetValue(str(state.qmax)) 
1109
1110    def _save_typeOfmodel(self):
1111        """
1112        save radiobutton containing the type model that can be selected
1113        """
1114        self.state.shape_rbutton = self.shape_rbutton.GetValue()
1115        self.state.shape_indep_rbutton = self.shape_indep_rbutton.GetValue()
1116        self.state.struct_rbutton = self.struct_rbutton.GetValue()
1117        self.state.plugin_rbutton = self.plugin_rbutton.GetValue()
1118        self.state.structurebox= self.structurebox.GetCurrentSelection()
1119        self.state.formfactorbox = self.formfactorbox.GetCurrentSelection()
1120       
1121        #self._undo.Enable(True)
1122        ## post state to fit panel
1123        event = PageInfoEvent(page = self)
1124        wx.PostEvent(self.parent, event)
1125       
1126    def _save_plotting_range(self ):
1127        """
1128        save the state of plotting range
1129        """
1130        self.state.qmin = self.qmin_x
1131        self.state.qmax = self.qmax_x
1132        self.state.npts = self.npts_x
1133           
1134    def _onparamEnter_helper(self):
1135        """
1136        check if values entered by the user are changed and valid to replot
1137        model
1138        """
1139        # Flag to register when a parameter has changed.   
1140        is_modified = False
1141        self.fitrange = True
1142        is_2Ddata = False
1143        #self._undo.Enable(True)
1144        # check if 2d data
1145        if self.data.__class__.__name__ == "Data2D":
1146            is_2Ddata = True
1147        if self.model !=None:
1148            try:
1149                is_modified = self._check_value_enter(self.fittable_param,
1150                                                     is_modified)
1151                is_modified = self._check_value_enter(self.fixed_param,
1152                                                      is_modified)
1153                is_modified = self._check_value_enter(self.parameters,
1154                                                      is_modified) 
1155            except:
1156                pass
1157            #if is_modified:
1158
1159            # Here we should check whether the boundaries have been modified.
1160            # If qmin and qmax have been modified, update qmin and qmax and
1161            # set the is_modified flag to True
1162            if self._validate_qrange(self.qmin, self.qmax):
1163                tempmin = float(self.qmin.GetValue())
1164                if tempmin != self.qmin_x:
1165                    self.qmin_x = tempmin
1166                    is_modified = True
1167                tempmax = float(self.qmax.GetValue())
1168                if tempmax != self.qmax_x:
1169                    self.qmax_x = tempmax
1170                    is_modified = True
1171           
1172                if is_2Ddata:
1173                    # set mask   
1174                    is_modified = self._validate_Npts()
1175                   
1176            else:
1177                self.fitrange = False   
1178           
1179            ## if any value is modify draw model with new value
1180            if not self.fitrange:
1181                #self.btFit.Disable()
1182                if is_2Ddata: self.btEditMask.Disable()
1183            else:
1184                #self.btFit.Enable(True)
1185                if is_2Ddata: self.btEditMask.Enable(True)
1186            if is_modified and self.fitrange:
1187                if self.data == None:
1188                    # Theory case: need to get npts value to draw
1189                    self.npts_x = float(self.Npts_total.GetValue())
1190                self.state_change= True
1191                self._draw_model() 
1192                self.Refresh()
1193        return is_modified
1194   
1195    def _update_paramv_on_fit(self):
1196        """
1197        make sure that update param values just before the fitting
1198        """
1199        #flag for qmin qmax check values
1200        flag = True
1201        self.fitrange = True
1202        is_modified = False
1203
1204        #wx.PostEvent(self._manager.parent, StatusEvent(status=" \
1205        #updating ... ",type="update"))
1206
1207        ##So make sure that update param values on_Fit.
1208        #self._undo.Enable(True)
1209        if self.model !=None:           
1210            ##Check the values
1211            self._check_value_enter( self.fittable_param ,is_modified)
1212            self._check_value_enter( self.fixed_param ,is_modified)
1213            self._check_value_enter( self.parameters ,is_modified)
1214
1215            # If qmin and qmax have been modified, update qmin and qmax and
1216             # Here we should check whether the boundaries have been modified.
1217            # If qmin and qmax have been modified, update qmin and qmax and
1218            # set the is_modified flag to True
1219            self.fitrange = self._validate_qrange(self.qmin, self.qmax)
1220            if self.fitrange:
1221                tempmin = float(self.qmin.GetValue())
1222                if tempmin != self.qmin_x:
1223                    self.qmin_x = tempmin
1224                tempmax = float(self.qmax.GetValue())
1225                if tempmax != self.qmax_x:
1226                    self.qmax_x = tempmax
1227                if tempmax == tempmin:
1228                    flag = False   
1229                temp_smearer = None
1230                if not self.disable_smearer.GetValue():
1231                    temp_smearer= self.current_smearer
1232                    if self.slit_smearer.GetValue():
1233                        flag = self.update_slit_smear()
1234                    elif self.pinhole_smearer.GetValue():
1235                        flag = self.update_pinhole_smear()
1236                    else:
1237                        self._manager.set_smearer(smearer=temp_smearer,
1238                                                  uid=self.uid,
1239                                                     qmin=float(self.qmin_x),
1240                                                      qmax=float(self.qmax_x),
1241                                                      draw=False)
1242                elif not self._is_2D():
1243                    self._manager.set_smearer(smearer=temp_smearer,
1244                                              qmin=float(self.qmin_x),
1245                                              uid=self.uid, 
1246                                                 qmax= float(self.qmax_x))
1247                    if self.data != None:
1248                        index_data = ((self.qmin_x <= self.data.x)&\
1249                                      (self.data.x <= self.qmax_x))
1250                        val = str(len(self.data.x[index_data==True]))
1251                        self.Npts_fit.SetValue(val)
1252                    else:
1253                        # No data in the panel
1254                        try:
1255                            self.npts_x = float(self.Npts_total.GetValue())
1256                        except:
1257                            flag = False
1258                            return flag
1259                    flag = True
1260                if self._is_2D():
1261                    # only 2D case set mask 
1262                    flag = self._validate_Npts()
1263                    if not flag:
1264                        return flag
1265            else: flag = False
1266        else: 
1267            flag = False
1268
1269        #For invalid q range, disable the mask editor and fit button, vs.   
1270        if not self.fitrange:
1271            #self.btFit.Disable()
1272            if self._is_2D():self.btEditMask.Disable()
1273        else:
1274            #self.btFit.Enable(True)
1275            if self._is_2D():self.btEditMask.Enable(True)
1276
1277        if not flag:
1278            msg = "Cannot Plot or Fit :Must select a "
1279            msg += " model or Fitting range is not valid!!!  "
1280            wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
1281       
1282        self.save_current_state()
1283   
1284        return flag                           
1285               
1286    def _is_modified(self, is_modified):
1287        """
1288        return to self._is_modified
1289        """
1290        return is_modified
1291                       
1292    def _reset_parameters_state(self, listtorestore, statelist):
1293        """
1294        Reset the parameters at the given state
1295        """
1296        if len(statelist) == 0 or len(listtorestore) == 0:
1297            return
1298        if len(statelist) !=  len(listtorestore):
1299            return
1300
1301        for j in range(len(listtorestore)):
1302            item_page = listtorestore[j]
1303            item_page_info = statelist[j]
1304            ##change the state of the check box for simple parameters
1305            if item_page[0]!=None:   
1306                item_page[0].SetValue(item_page_info[0])
1307            if item_page[2]!=None:
1308                item_page[2].SetValue(item_page_info[2])
1309                if item_page[2].__class__.__name__ == "ComboBox":
1310                   if self.model.fun_list.has_key(item_page_info[2]):
1311                       fun_val = self.model.fun_list[item_page_info[2]]
1312                       self.model.setParam(item_page_info[1],fun_val)
1313            if item_page[3]!=None:
1314                ## show or hide text +/-
1315                if item_page_info[2]:
1316                    item_page[3].Show(True)
1317                else:
1318                    item_page[3].Hide()
1319            if item_page[4]!=None:
1320                ## show of hide the text crtl for fitting error
1321                if item_page_info[4][0]:
1322                    item_page[4].Show(True)
1323                    item_page[4].SetValue(item_page_info[4][1])
1324                else:
1325                    item_page[3].Hide()
1326            if item_page[5]!=None:
1327                ## show of hide the text crtl for fitting error
1328                item_page[5].Show(item_page_info[5][0])
1329                item_page[5].SetValue(item_page_info[5][1])
1330               
1331            if item_page[6]!=None:
1332                ## show of hide the text crtl for fitting error
1333                item_page[6].Show(item_page_info[6][0])
1334                item_page[6].SetValue(item_page_info[6][1])
1335
1336                   
1337    def _reset_strparam_state(self, listtorestore, statelist):
1338        """
1339        Reset the string parameters at the given state
1340        """
1341        if len(statelist) == 0:
1342            return
1343
1344        listtorestore = copy.deepcopy(statelist)
1345       
1346        for j in range(len(listtorestore)):
1347            item_page = listtorestore[j]
1348            item_page_info = statelist[j]
1349            ##change the state of the check box for simple parameters
1350           
1351            if item_page[0] != None:   
1352                item_page[0].SetValue(format_number(item_page_info[0], True))
1353
1354            if item_page[2] != None:
1355                param_name = item_page_info[1]
1356                value = item_page_info[2]
1357                selection = value
1358                if self.model.fun_list.has_key(value):
1359                    selection = self.model.fun_list[value]
1360                item_page[2].SetValue(selection)
1361                self.model.setParam(param_name, selection)
1362                                     
1363    def _copy_parameters_state(self, listtocopy, statelist):
1364        """
1365        copy the state of button
1366       
1367        :param listtocopy: the list of check button to copy
1368        :param statelist: list of state object to store the current state
1369       
1370        """
1371        if len(listtocopy)==0:
1372            return
1373       
1374        for item in listtocopy:
1375 
1376            checkbox_state = None
1377            if item[0]!= None:
1378                checkbox_state= item[0].GetValue()
1379            parameter_name = item[1]
1380            parameter_value = None
1381            if item[2]!=None:
1382                parameter_value = item[2].GetValue()
1383            static_text = None
1384            if item[3]!=None:
1385                static_text = item[3].IsShown()
1386            error_value = None
1387            error_state = None
1388            if item[4]!= None:
1389                error_value = item[4].GetValue()
1390                error_state = item[4].IsShown()
1391               
1392            min_value = None
1393            min_state = None
1394            if item[5]!= None:
1395                min_value = item[5].GetValue()
1396                min_state = item[5].IsShown()
1397               
1398            max_value = None
1399            max_state = None
1400            if item[6]!= None:
1401                max_value = item[6].GetValue()
1402                max_state = item[6].IsShown()
1403            unit=None
1404            if item[7]!=None:
1405                unit = item[7].GetLabel()
1406               
1407            statelist.append([checkbox_state, parameter_name, parameter_value,
1408                              static_text ,[error_state, error_value],
1409                                [min_state, min_value],
1410                                [max_state, max_value], unit])
1411           
1412    def _set_model_sizer_selection(self, model):
1413        """
1414        Display the sizer according to the type of the current model
1415        """
1416        if model == None:
1417            return
1418        if hasattr(model ,"s_model"):
1419           
1420            class_name = model.s_model.__class__
1421            name = model.s_model.name
1422            flag = (name != "NoStructure")
1423            if flag and \
1424                (class_name in self.model_list_box["Structure Factors"]):
1425                self.structurebox.Show()
1426                self.text2.Show()               
1427                self.structurebox.Enable()
1428                self.text2.Enable()
1429                items = self.structurebox.GetItems()
1430                self.sizer1.Layout()
1431               
1432                for i in range(len(items)):
1433                    if items[i]== str(name):
1434                        self.structurebox.SetSelection(i)
1435                        break
1436                   
1437        if hasattr(model ,"p_model"):
1438            class_name = model.p_model.__class__
1439            name = model.p_model.name
1440            self.formfactorbox.Clear()
1441           
1442            for k, list in self.model_list_box.iteritems():
1443                if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
1444                    self.shape_rbutton.SetValue(True)
1445                    ## fill the form factor list with new model
1446                    self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
1447                    items = self.formfactorbox.GetItems()
1448                    ## set comboxbox to the selected item
1449                    for i in range(len(items)):
1450                        if items[i]== str(name):
1451                            self.formfactorbox.SetSelection(i)
1452                            break
1453                    return
1454                elif k == "Shape-Independent":
1455                    self.shape_indep_rbutton.SetValue(True)
1456                elif k == "Structure Factors":
1457                     self.struct_rbutton.SetValue(True)
1458                elif  k == "Multi-Functions":
1459                    continue
1460                else:
1461                    self.plugin_rbutton.SetValue(True)
1462               
1463                if class_name in list:
1464                    ## fill the form factor list with new model
1465                    self._populate_box(self.formfactorbox, list)
1466                    items = self.formfactorbox.GetItems()
1467                    ## set comboxbox to the selected item
1468                    for i in range(len(items)):
1469                        if items[i]== str(name):
1470                            self.formfactorbox.SetSelection(i)
1471                            break
1472                    break
1473        else:
1474
1475            ## Select the model from the menu
1476            class_name = model.__class__
1477            name = model.name
1478            self.formfactorbox.Clear()
1479            items = self.formfactorbox.GetItems()
1480   
1481            for k, list in self.model_list_box.iteritems():         
1482                if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
1483                    if class_name in self.model_list_box["P(Q)*S(Q)"]:
1484                        self.structurebox.Show()
1485                        self.text2.Show()
1486                        self.structurebox.Enable()
1487                        self.structurebox.SetSelection(0)
1488                        self.text2.Enable()
1489                    else:
1490                        self.structurebox.Hide()
1491                        self.text2.Hide()
1492                        self.structurebox.Disable()
1493                        self.structurebox.SetSelection(0)
1494                        self.text2.Disable()
1495                       
1496                    self.shape_rbutton.SetValue(True)
1497                    ## fill the form factor list with new model
1498                    self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
1499                    items = self.formfactorbox.GetItems()
1500                    ## set comboxbox to the selected item
1501                    for i in range(len(items)):
1502                        if items[i]== str(name):
1503                            self.formfactorbox.SetSelection(i)
1504                            break
1505                    return
1506                elif k == "Shape-Independent":
1507                    self.shape_indep_rbutton.SetValue(True)
1508                elif k == "Structure Factors":
1509                    self.struct_rbutton.SetValue(True)
1510                elif  k == "Multi-Functions":
1511                    continue
1512                else:
1513                    self.plugin_rbutton.SetValue(True)
1514                if class_name in list:
1515                    self.structurebox.SetSelection(0)
1516                    self.structurebox.Disable()
1517                    self.text2.Disable()                   
1518                    ## fill the form factor list with new model
1519                    self._populate_box(self.formfactorbox, list)
1520                    items = self.formfactorbox.GetItems()
1521                    ## set comboxbox to the selected item
1522                    for i in range(len(items)):
1523                        if items[i]== str(name):
1524                            self.formfactorbox.SetSelection(i)
1525                            break
1526                    break
1527   
1528    def _draw_model(self, update_chisqr=True):
1529        """
1530        Method to draw or refresh a plotted model.
1531        The method will use the data member from the model page
1532        to build a call to the fitting perspective manager.
1533       
1534        :param chisqr: update chisqr value [bool]
1535        """
1536        #if self.check_invalid_panel():
1537        #    return
1538        if self.model !=None:
1539            temp_smear=None
1540            if hasattr(self, "enable_smearer"):
1541                if not self.disable_smearer.GetValue():
1542                    temp_smear= self.current_smearer
1543            toggle_mode_on = self.model_view.IsEnabled()
1544            self._manager.draw_model(self.model, 
1545                                    data=self.data,
1546                                    smearer= temp_smear,
1547                                    qmin=float(self.qmin_x), 
1548                                    qmax=float(self.qmax_x),
1549                                    qstep= float(self.npts_x),
1550                                    page_id=self.uid,
1551                                    toggle_mode_on=toggle_mode_on, 
1552                                    state = self.state,
1553                                    enable2D=self.enable2D,
1554                                    update_chisqr=update_chisqr)
1555       
1556       
1557    def _on_show_sld(self, event=None):
1558        """
1559        Plot SLD profile
1560        """
1561        # get profile data
1562        x,y=self.model.getProfile()
1563
1564        from danse.common.plottools import Data1D
1565        #from sans.perspectives.theory.profile_dialog import SLDPanel
1566        from sans.guiframe.local_perspectives.plotting.profile_dialog \
1567        import SLDPanel
1568        sld_data = Data1D(x,y)
1569        sld_data.name = 'SLD'
1570        sld_data.axes = self.sld_axes
1571        self.panel = SLDPanel(self, data=sld_data,axes =self.sld_axes,id =-1 )
1572        self.panel.ShowModal()   
1573       
1574    def _set_multfactor_combobox(self, multiplicity=10):   
1575        """
1576        Set comboBox for muitfactor of CoreMultiShellModel
1577        :param multiplicit: no. of multi-functionality
1578        """
1579        # build content of the combobox
1580        for idx in range(0,multiplicity):
1581            self.multifactorbox.Append(str(idx),int(idx))
1582            #self.multifactorbox.SetSelection(1)
1583        self._hide_multfactor_combobox()
1584       
1585    def _show_multfactor_combobox(self):   
1586        """
1587        Show the comboBox of muitfactor of CoreMultiShellModel
1588        """ 
1589        if not self.mutifactor_text.IsShown():
1590            self.mutifactor_text.Show(True)
1591            self.mutifactor_text1.Show(True)
1592        if not self.multifactorbox.IsShown():
1593            self.multifactorbox.Show(True) 
1594             
1595    def _hide_multfactor_combobox(self):   
1596        """
1597        Hide the comboBox of muitfactor of CoreMultiShellModel
1598        """ 
1599        if self.mutifactor_text.IsShown():
1600            self.mutifactor_text.Hide()
1601            self.mutifactor_text1.Hide()
1602        if self.multifactorbox.IsShown():
1603            self.multifactorbox.Hide()   
1604
1605       
1606    def _show_combox_helper(self):
1607        """
1608        Fill panel's combo box according to the type of model selected
1609        """
1610        if self.shape_rbutton.GetValue():
1611            ##fill the combobox with form factor list
1612            self.structurebox.SetSelection(0)
1613            self.structurebox.Disable()
1614            self.formfactorbox.Clear()
1615            self._populate_box( self.formfactorbox,self.model_list_box["Shapes"])
1616        if self.shape_indep_rbutton.GetValue():
1617            ##fill the combobox with shape independent  factor list
1618            self.structurebox.SetSelection(0)
1619            self.structurebox.Disable()
1620            self.formfactorbox.Clear()
1621            self._populate_box( self.formfactorbox,
1622                                self.model_list_box["Shape-Independent"])
1623        if self.struct_rbutton.GetValue():
1624            ##fill the combobox with structure factor list
1625            self.structurebox.SetSelection(0)
1626            self.structurebox.Disable()
1627            self.formfactorbox.Clear()
1628            self._populate_box( self.formfactorbox,
1629                                self.model_list_box["Structure Factors"])
1630        if self.plugin_rbutton.GetValue():
1631            ##fill the combobox with form factor list
1632            self.structurebox.Disable()
1633            self.formfactorbox.Clear()
1634            self._populate_box( self.formfactorbox,
1635                                self.model_list_box["Customized Models"])
1636       
1637    def _show_combox(self, event=None):
1638        """
1639        Show combox box associate with type of model selected
1640        """
1641        #if self.check_invalid_panel():
1642        #    self.shape_rbutton.SetValue(True)
1643        #    return
1644
1645        self._show_combox_helper()
1646        self._on_select_model(event=None)
1647        self._save_typeOfmodel()
1648        self.sizer4_4.Layout()
1649        self.sizer4.Layout()
1650        self.Layout()
1651        self.Refresh()
1652 
1653    def _populate_box(self, combobox, list):
1654        """
1655        fill combox box with dict item
1656       
1657        :param list: contains item to fill the combox
1658            item must model class
1659        """
1660        for models in list:
1661            model= models()
1662            name = model.__class__.__name__
1663            if models.__name__!="NoStructure":
1664                if hasattr(model, "name"):
1665                    name = model.name
1666                combobox.Append(name,models)
1667        return 0
1668   
1669    def _onQrangeEnter(self, event):
1670        """
1671        Check validity of value enter in the Q range field
1672       
1673        """
1674        tcrtl = event.GetEventObject()
1675        #Clear msg if previously shown.
1676        msg = ""
1677        wx.PostEvent(self.parent, StatusEvent(status=msg))
1678        # Flag to register when a parameter has changed.
1679        is_modified = False
1680        if tcrtl.GetValue().lstrip().rstrip() != "":
1681            try:
1682                value = float(tcrtl.GetValue())
1683                tcrtl.SetBackgroundColour(wx.WHITE)
1684                # If qmin and qmax have been modified, update qmin and qmax
1685                if self._validate_qrange(self.qmin, self.qmax):
1686                    tempmin = float(self.qmin.GetValue())
1687                    if tempmin != self.qmin_x:
1688                        self.qmin_x = tempmin
1689                    tempmax = float(self.qmax.GetValue())
1690                    if tempmax != self.qmax_x:
1691                        self.qmax_x = tempmax
1692                else:
1693                    tcrtl.SetBackgroundColour("pink")
1694                    msg = "Model Error:wrong value entered : %s" % sys.exc_value
1695                    wx.PostEvent(self.parent, StatusEvent(status=msg))
1696                    return 
1697            except:
1698                tcrtl.SetBackgroundColour("pink")
1699                msg = "Model Error:wrong value entered : %s" % sys.exc_value
1700                wx.PostEvent(self.parent, StatusEvent(status=msg))
1701                return 
1702            #Check if # of points for theory model are valid(>0).
1703            if self.npts != None:
1704                if check_float(self.npts):
1705                    temp_npts = float(self.npts.GetValue())
1706                    if temp_npts !=  self.num_points:
1707                        self.num_points = temp_npts
1708                        is_modified = True
1709                else:
1710                    msg = "Cannot Plot :No npts in that Qrange!!!  "
1711                    wx.PostEvent(self.parent, StatusEvent(status=msg))
1712        else:
1713           tcrtl.SetBackgroundColour("pink")
1714           msg = "Model Error:wrong value entered!!!"
1715           wx.PostEvent(self.parent, StatusEvent(status=msg))
1716        #self._undo.Enable(True)
1717        self.save_current_state()
1718        event = PageInfoEvent(page=self)
1719        wx.PostEvent(self.parent, event)
1720        self.state_change = False
1721        #Draw the model for a different range
1722        self._draw_model()
1723                   
1724    def _theory_qrange_enter(self, event):
1725        """
1726        Check validity of value enter in the Q range field
1727        """
1728       
1729        tcrtl= event.GetEventObject()
1730        #Clear msg if previously shown.
1731        msg= ""
1732        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1733        # Flag to register when a parameter has changed.
1734        is_modified = False
1735        if tcrtl.GetValue().lstrip().rstrip()!="":
1736            try:
1737                value = float(tcrtl.GetValue())
1738                tcrtl.SetBackgroundColour(wx.WHITE)
1739
1740                # If qmin and qmax have been modified, update qmin and qmax
1741                if self._validate_qrange(self.theory_qmin, self.theory_qmax):
1742                    tempmin = float(self.theory_qmin.GetValue())
1743                    if tempmin != self.theory_qmin_x:
1744                        self.theory_qmin_x = tempmin
1745                    tempmax = float(self.theory_qmax.GetValue())
1746                    if tempmax != self.qmax_x:
1747                        self.theory_qmax_x = tempmax
1748                else:
1749                    tcrtl.SetBackgroundColour("pink")
1750                    msg= "Model Error:wrong value entered : %s"% sys.exc_value
1751                    wx.PostEvent(self._manager.parent, StatusEvent(status = msg ))
1752                    return 
1753            except:
1754                tcrtl.SetBackgroundColour("pink")
1755                msg= "Model Error:wrong value entered : %s"% sys.exc_value
1756                wx.PostEvent(self._manager.parent, StatusEvent(status = msg ))
1757                return 
1758            #Check if # of points for theory model are valid(>0).
1759            if self.Npts_total.IsEditable() :
1760                if check_float(self.Npts_total):
1761                    temp_npts = float(self.Npts_total.GetValue())
1762                    if temp_npts !=  self.num_points:
1763                        self.num_points = temp_npts
1764                        is_modified = True
1765                else:
1766                    msg= "Cannot Plot :No npts in that Qrange!!!  "
1767                    wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1768        else:
1769           tcrtl.SetBackgroundColour("pink")
1770           msg = "Model Error:wrong value entered!!!"
1771           wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
1772        #self._undo.Enable(True)
1773        self.save_current_state()
1774        event = PageInfoEvent(page = self)
1775        wx.PostEvent(self.parent, event)
1776        self.state_change= False
1777        #Draw the model for a different range
1778        self._draw_model()
1779                   
1780    def _on_select_model_helper(self): 
1781        """
1782        call back for model selection
1783        """
1784        ## reset dictionary containing reference to dispersion
1785        self._disp_obj_dict = {}
1786        self.disp_cb_dict ={}
1787       
1788        f_id = self.formfactorbox.GetCurrentSelection()
1789        #For MAC
1790        form_factor = None
1791        if f_id >= 0:
1792            form_factor = self.formfactorbox.GetClientData(f_id)
1793
1794        if not form_factor in  self.model_list_box["multiplication"]:
1795            self.structurebox.Hide()
1796            self.text2.Hide()           
1797            self.structurebox.Disable()
1798            self.structurebox.SetSelection(0)
1799            self.text2.Disable()
1800        else:
1801            self.structurebox.Show()
1802            self.text2.Show()
1803            self.structurebox.Enable()
1804            self.text2.Enable()
1805           
1806        if form_factor != None:   
1807            # set multifactor for Mutifunctional models   
1808            if form_factor().__class__ in self.model_list_box["Multi-Functions"]:
1809                m_id = self.multifactorbox.GetCurrentSelection()
1810                multiplicity = form_factor().multiplicity_info[0]
1811                self.multifactorbox.Clear()
1812                #self.mutifactor_text.SetLabel(form_factor().details[])
1813                self._set_multfactor_combobox(multiplicity)
1814                self._show_multfactor_combobox()
1815                #ToDo:  this info should be called directly from the model
1816                text = form_factor().multiplicity_info[1]#'No. of Shells: '
1817
1818                #self.mutifactor_text.Clear()
1819                self.mutifactor_text.SetLabel(text)
1820                if m_id > multiplicity -1:
1821                    # default value
1822                    m_id = 1
1823                   
1824                self.multi_factor = self.multifactorbox.GetClientData(m_id)
1825                if self.multi_factor == None: self.multi_factor =0
1826                form_factor = form_factor(int(self.multi_factor))
1827                self.multifactorbox.SetSelection(m_id)
1828                # Check len of the text1 and max_multiplicity
1829                text = ''
1830                if form_factor.multiplicity_info[0] == len(form_factor.multiplicity_info[2]):
1831                    text = form_factor.multiplicity_info[2][self.multi_factor]
1832                self.mutifactor_text1.SetLabel(text)
1833                # Check if model has  get sld profile.
1834                if len(form_factor.multiplicity_info[3]) > 0:
1835                    self.sld_axes = form_factor.multiplicity_info[3]
1836                    self.show_sld_button.Show(True)
1837                else:
1838                    self.sld_axes = ""
1839
1840            else:
1841                self._hide_multfactor_combobox()
1842                self.show_sld_button.Hide()
1843                form_factor = form_factor()
1844                self.multi_factor = None
1845        else:
1846            self._hide_multfactor_combobox()
1847            self.show_sld_button.Hide()
1848            self.multi_factor = None 
1849             
1850        s_id = self.structurebox.GetCurrentSelection()
1851        struct_factor = self.structurebox.GetClientData( s_id )
1852       
1853        if  struct_factor !=None:
1854            from sans.models.MultiplicationModel import MultiplicationModel
1855            self.model= MultiplicationModel(form_factor,struct_factor())
1856           
1857        else:
1858            if form_factor != None:
1859                self.model= form_factor
1860            else:
1861                self.model = None
1862                return self.model
1863           
1864
1865        ## post state to fit panel
1866        self.state.parameters =[]
1867        self.state.model =self.model
1868        self.state.qmin = self.qmin_x
1869        self.state.multi_factor = self.multi_factor
1870        self.disp_list =self.model.getDispParamList()
1871        self.state.disp_list = self.disp_list
1872        self.Layout()     
1873       
1874    def _validate_qrange(self, qmin_ctrl, qmax_ctrl):
1875        """
1876        Verify that the Q range controls have valid values
1877        and that Qmin < Qmax.
1878       
1879        :param qmin_ctrl: text control for Qmin
1880        :param qmax_ctrl: text control for Qmax
1881       
1882        :return: True is the Q range is value, False otherwise
1883       
1884        """
1885        qmin_validity = check_float(qmin_ctrl)
1886        qmax_validity = check_float(qmax_ctrl)
1887        if not (qmin_validity and qmax_validity):
1888            return False
1889        else:
1890            qmin = float(qmin_ctrl.GetValue())
1891            qmax = float(qmax_ctrl.GetValue())
1892            if qmin <= qmax:
1893                #Make sure to set both colours white. 
1894                qmin_ctrl.SetBackgroundColour(wx.WHITE)
1895                qmin_ctrl.Refresh()
1896                qmax_ctrl.SetBackgroundColour(wx.WHITE)
1897                qmax_ctrl.Refresh()
1898            else:
1899                qmin_ctrl.SetBackgroundColour("pink")
1900                qmin_ctrl.Refresh()
1901                qmax_ctrl.SetBackgroundColour("pink")
1902                qmax_ctrl.Refresh()
1903                msg= "Invalid Q range: Q min must be smaller than Q max"
1904                wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
1905                return False
1906        return True
1907   
1908    def _validate_Npts(self): 
1909        """
1910        Validate the number of points for fitting is more than 10 points.
1911        If valid, setvalues Npts_fit otherwise post msg.
1912        """
1913        #default flag
1914        flag = True
1915        # Theory
1916        if self.data == None and self.enable2D:
1917            return flag
1918
1919        # q value from qx and qy
1920        radius= numpy.sqrt( self.data.qx_data * self.data.qx_data + 
1921                            self.data.qy_data * self.data.qy_data )
1922        #get unmasked index
1923        index_data = (float(self.qmin.GetValue()) <= radius) & \
1924                        (radius <= float(self.qmax.GetValue()))
1925        index_data = (index_data) & (self.data.mask) 
1926        index_data = (index_data) & (numpy.isfinite(self.data.data))
1927
1928        if len(index_data[index_data]) < 10:
1929            # change the color pink.
1930            self.qmin.SetBackgroundColour("pink")
1931            self.qmin.Refresh()
1932            self.qmax.SetBackgroundColour("pink")
1933            self.qmax.Refresh()
1934            msg= "Cannot Plot :No or too little npts in that data range!!!  "
1935            wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1936            self.fitrange = False
1937            flag = False
1938        else:
1939            self.Npts_fit.SetValue(str(len(self.data.mask[index_data==True])))
1940            self.fitrange = True
1941           
1942        return flag
1943   
1944    def _check_value_enter(self, list, modified):
1945        """
1946        :param list: model parameter and panel info
1947        :Note: each item of the list should be as follow:
1948            item=[check button state, parameter's name,
1949                paramater's value, string="+/-",
1950                parameter's error of fit,
1951                parameter's minimum value,
1952                parrameter's maximum value ,
1953                parameter's units]
1954        """ 
1955        is_modified =  modified
1956        if len(list)==0:
1957            return is_modified
1958        for item in list:
1959            #skip angle parameters for 1D
1960            if not self.enable2D:#self.data.__class__.__name__ !="Data2D":
1961                if item in self.orientation_params:
1962                    continue
1963            #try:
1964            name = str(item[1])
1965           
1966            if string.find(name,".npts") ==-1 and \
1967                                        string.find(name,".nsigmas")==-1:     
1968                ## check model parameters range             
1969                param_min= None
1970                param_max= None
1971               
1972                ## check minimun value
1973                if item[5]!= None and item[5]!= "":
1974                    if item[5].GetValue().lstrip().rstrip()!="":
1975                        try:
1976                           
1977                            param_min = float(item[5].GetValue())
1978                            if not self._validate_qrange(item[5],item[2]):
1979                                if numpy.isfinite(param_min):
1980                                    item[2].SetValue(format_number(param_min))
1981                           
1982                            item[5].SetBackgroundColour(wx.WHITE)
1983                            item[2].SetBackgroundColour(wx.WHITE)
1984                                           
1985                        except:
1986                            msg = "Wrong Fit parameter range entered "
1987                            wx.PostEvent(self.parent.parent, 
1988                                         StatusEvent(status = msg))
1989                            raise ValueError, msg
1990                        is_modified = True
1991                ## check maximum value
1992                if item[6]!= None and item[6]!= "":
1993                    if item[6].GetValue().lstrip().rstrip()!="":
1994                        try:                         
1995                            param_max = float(item[6].GetValue())
1996                            if not self._validate_qrange(item[2],item[6]):
1997                                if numpy.isfinite(param_max):
1998                                    item[2].SetValue(format_number(param_max)) 
1999                           
2000                            item[6].SetBackgroundColour(wx.WHITE)
2001                            item[2].SetBackgroundColour(wx.WHITE)
2002                        except:
2003                            msg = "Wrong Fit parameter range entered "
2004                            wx.PostEvent(self.parent.parent, 
2005                                         StatusEvent(status = msg))
2006                            raise ValueError, msg
2007                        is_modified = True
2008               
2009
2010                if param_min != None and param_max !=None:
2011                    if not self._validate_qrange(item[5], item[6]):
2012                        msg= "Wrong Fit range entered for parameter "
2013                        msg+= "name %s of model %s "%(name, self.model.name)
2014                        wx.PostEvent(self.parent.parent, 
2015                                     StatusEvent(status = msg))
2016               
2017                if name in self.model.details.keys():   
2018                        self.model.details[name][1:3] = param_min, param_max
2019                        is_modified = True
2020             
2021                else:
2022                        self.model.details [name] = ["", param_min, param_max] 
2023                        is_modified = True
2024            try:   
2025                # Check if the textctr is enabled
2026                if item[2].IsEnabled():
2027                    value= float(item[2].GetValue())
2028                    item[2].SetBackgroundColour("white")
2029                    # If the value of the parameter has changed,
2030                    # +update the model and set the is_modified flag
2031                    if value != self.model.getParam(name) and \
2032                                                numpy.isfinite(value):
2033                        self.model.setParam(name, value)
2034                       
2035            except:
2036                item[2].SetBackgroundColour("pink")
2037                msg = "Wrong Fit parameter value entered "
2038                wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
2039               
2040        return is_modified
2041       
2042 
2043    def _set_dipers_Param(self, event):
2044        """
2045        respond to self.enable_disp and self.disable_disp radio box.
2046        The dispersity object is reset inside the model into Gaussian.
2047        When the user select yes , this method display a combo box for more selection
2048        when the user selects No,the combo box disappears.
2049        Redraw the model with the default dispersity (Gaussian)
2050        """
2051        #if self.check_invalid_panel():
2052        #    return
2053        ## On selction if no model exists.
2054        if self.model ==None:
2055            self.disable_disp.SetValue(True)
2056            msg="Please select a Model first..."
2057            wx.MessageBox(msg, 'Info')
2058            wx.PostEvent(self._manager.parent, StatusEvent(status=\
2059                            "Polydispersion: %s"%msg))
2060            return
2061
2062        self._reset_dispersity()
2063   
2064        if self.model ==None:
2065            self.model_disp.Hide()
2066            self.sizer4_4.Clear(True)
2067            return
2068
2069        if self.enable_disp.GetValue():
2070            ## layout for model containing no dispersity parameters
2071           
2072            self.disp_list= self.model.getDispParamList()
2073             
2074            if len(self.disp_list)==0 and len(self.disp_cb_dict)==0:
2075                self._layout_sizer_noDipers() 
2076            else:
2077                ## set gaussian sizer
2078                self._on_select_Disp(event=None)
2079        else:
2080            self.sizer4_4.Clear(True)
2081           
2082        ## post state to fit panel
2083        self.save_current_state()
2084        if event !=None:
2085            #self._undo.Enable(True)
2086            event = PageInfoEvent(page = self)
2087            wx.PostEvent(self.parent, event)
2088        #draw the model with the current dispersity
2089        self._draw_model()
2090        self.sizer4_4.Layout()
2091        self.sizer5.Layout()
2092        self.Layout()
2093        self.Refresh()     
2094         
2095       
2096    def _layout_sizer_noDipers(self):
2097        """
2098        Draw a sizer with no dispersity info
2099        """
2100        ix=0
2101        iy=1
2102        self.fittable_param=[]
2103        self.fixed_param=[]
2104        self.orientation_params_disp=[]
2105       
2106        self.sizer4_4.Clear(True)
2107        text = "No polydispersity available for this model"
2108        text = "No polydispersity available for this model"
2109        model_disp = wx.StaticText(self, -1, text)
2110        self.sizer4_4.Add(model_disp,( iy, ix),(1,1), 
2111                            wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 10)
2112        self.sizer4_4.Layout()
2113        self.sizer4.Layout()
2114   
2115    def _reset_dispersity(self):
2116        """
2117        put gaussian dispersity into current model
2118        """
2119        if len(self.param_toFit)>0:
2120            for item in self.fittable_param:
2121                if item in self.param_toFit:
2122                    self.param_toFit.remove(item)
2123
2124            for item in self.orientation_params_disp:
2125                if item in self.param_toFit:
2126                    self.param_toFit.remove(item)
2127         
2128        self.fittable_param=[]
2129        self.fixed_param=[]
2130        self.orientation_params_disp=[]
2131        self.values={}
2132        self.weights={}
2133     
2134        from sans.models.dispersion_models import GaussianDispersion, ArrayDispersion
2135        if len(self.disp_cb_dict)==0:
2136            self.save_current_state()
2137            self.sizer4_4.Clear(True)
2138            self.Layout()
2139 
2140            return 
2141        if (len(self.disp_cb_dict)>0) :
2142            for p in self.disp_cb_dict:
2143                # The parameter was un-selected. Go back to Gaussian model (with 0 pts)                   
2144                disp_model = GaussianDispersion()
2145               
2146                self._disp_obj_dict[p] = disp_model
2147                # Set the new model as the dispersion object for the selected parameter
2148                try:
2149                   self.model.set_dispersion(p, disp_model)
2150                except:
2151
2152                    pass
2153
2154        ## save state into
2155        self.save_current_state()
2156        self.Layout() 
2157        self.Refresh()
2158                 
2159    def _on_select_Disp(self,event):
2160        """
2161        allow selecting different dispersion
2162        self.disp_list should change type later .now only gaussian
2163        """
2164        self._set_sizer_dispersion()
2165
2166        ## Redraw the model
2167        self._draw_model() 
2168        #self._undo.Enable(True)
2169        event = PageInfoEvent(page = self)
2170        wx.PostEvent(self.parent, event)
2171       
2172        self.sizer4_4.Layout()
2173        self.sizer4.Layout()
2174        self.SetupScrolling()
2175   
2176    def _on_disp_func(self, event=None): 
2177        """
2178        Select a distribution function for the polydispersion
2179       
2180        :Param event: ComboBox event
2181        """
2182        # get ready for new event
2183        if event != None:
2184            event.Skip()
2185        # Get event object
2186        disp_box =  event.GetEventObject() 
2187
2188        # Try to select a Distr. function
2189        try:   
2190            disp_box.SetBackgroundColour("white")
2191            selection = disp_box.GetCurrentSelection()
2192            param_name = disp_box.Name.split('.')[0]
2193            disp_name = disp_box.GetValue()
2194            dispersity= disp_box.GetClientData(selection)
2195   
2196            #disp_model =  GaussianDispersion()
2197            disp_model = dispersity()
2198            # Get param names to reset the values of the param
2199            name1 = param_name + ".width"
2200            name2 = param_name + ".npts"
2201            name3 = param_name + ".nsigmas"
2202            # Check Disp. function whether or not it is 'array'
2203            if disp_name.lower() == "array":
2204                value2= ""
2205                value3= ""
2206                value1 = self._set_array_disp(name=name1, disp=disp_model)
2207            else:
2208                self._del_array_values(name1)
2209                #self._reset_array_disp(param_name)
2210                self._disp_obj_dict[name1] = disp_model
2211                self.model.set_dispersion(param_name, disp_model)
2212                self.state._disp_obj_dict[name1]= disp_model
2213 
2214                value1= str(format_number(self.model.getParam(name1), True))
2215                value2= str(format_number(self.model.getParam(name2)))
2216                value3= str(format_number(self.model.getParam(name3)))
2217            # Reset fittable polydispersin parameter value
2218            for item in self.fittable_param:
2219                 if item[1] == name1:
2220                    item[2].SetValue(value1) 
2221                    item[5].SetValue("")
2222                    item[6].SetValue("")
2223                    # Disable for array
2224                    if disp_name.lower() == "array":
2225                        item[0].SetValue(False)
2226                        item[0].Disable()
2227                        item[2].Disable()
2228                        item[5].Disable()
2229                        item[6].Disable()
2230                    else:
2231                        item[0].Enable()
2232                        item[2].Enable()
2233                        item[5].Enable()
2234                        item[6].Enable()                       
2235                    break
2236            # Reset fixed polydispersion params
2237            for item in self.fixed_param:
2238                if item[1] == name2:
2239                    item[2].SetValue(value2) 
2240                    # Disable Npts for array
2241                    if disp_name.lower() == "array":
2242                        item[2].Disable()
2243                    else:
2244                        item[2].Enable()
2245                if item[1] == name3:
2246                    item[2].SetValue(value3) 
2247                    # Disable Nsigs for array
2248                    if disp_name.lower() == "array":
2249                        item[2].Disable()
2250                    else:
2251                        item[2].Enable()
2252               
2253            # Make sure the check box updated when all checked
2254            if self.cb1.GetValue():
2255                self.select_all_param(None)
2256
2257            # update params
2258            self._update_paramv_on_fit() 
2259            # draw
2260            self._draw_model()
2261            self.Refresh()
2262        except:
2263            # Error msg
2264            msg = "Error occurred:"
2265            msg += " Could not select the distribution function..."
2266            msg += " Please select another distribution function."
2267            disp_box.SetBackgroundColour("pink")
2268            # Focus on Fit button so that users can see the pinky box
2269            self.btFit.SetFocus()
2270            wx.PostEvent(self.parent.parent, 
2271                         StatusEvent(status=msg, info="error"))
2272       
2273       
2274    def _set_array_disp(self, name=None, disp=None):
2275        """
2276        Set array dispersion
2277       
2278        :param name: name of the parameter for the dispersion to be set
2279        :param disp: the polydisperion object
2280        """
2281        # The user wants this parameter to be averaged.
2282        # Pop up the file selection dialog.
2283        path = self._selectDlg()
2284        # Array data
2285        values = []
2286        weights = []
2287        # If nothing was selected, just return
2288        if path is None:
2289            self.disp_cb_dict[name].SetValue(False)
2290            #self.noDisper_rbox.SetValue(True)
2291            return
2292        self._default_save_location = os.path.dirname(path)
2293
2294        basename  = os.path.basename(path)
2295        values,weights = self.read_file(path)
2296       
2297        # If any of the two arrays is empty, notify the user that we won't
2298        # proceed
2299        if len(self.param_toFit)>0:
2300            if name in self.param_toFit:
2301                self.param_toFit.remove(name)
2302
2303        # Tell the user that we are about to apply the distribution
2304        msg = "Applying loaded %s distribution: %s" % (name, path)
2305        wx.PostEvent(self.parent.parent, StatusEvent(status=msg)) 
2306
2307        disp.set_weights(values, weights)
2308        self._disp_obj_dict[name] = disp
2309        self.model.set_dispersion(name.split('.')[0], disp)
2310        self.state._disp_obj_dict[name]= disp
2311        self.values[name] = values
2312        self.weights[name] = weights
2313        # Store the object to make it persist outside the
2314        # scope of this method
2315        #TODO: refactor model to clean this up?
2316        self.state.values = {}
2317        self.state.weights = {}
2318        self.state.values = copy.deepcopy(self.values)
2319        self.state.weights = copy.deepcopy(self.weights)
2320
2321        # Set the new model as the dispersion object for the
2322        #selected parameter
2323        #self.model.set_dispersion(p, disp_model)
2324        # Store a reference to the weights in the model object
2325        #so that
2326        # it's not lost when we use the model within another thread.
2327        #TODO: total hack - fix this
2328        self.state.model= self.model.clone()
2329        self.model._persistency_dict[name.split('.')[0]] = \
2330                                        [values, weights]
2331        self.state.model._persistency_dict[name.split('.')[0]] = \
2332                                        [values,weights]
2333        return basename
2334   
2335    def _del_array_values(self, name=None): 
2336        """
2337        Reset array dispersion
2338       
2339        :param name: name of the parameter for the dispersion to be set
2340        """
2341        # Try to delete values and weight of the names array dic if exists
2342        try:
2343            del self.values[name]
2344            del self.weights[name]
2345            # delete all other dic
2346            del self.state.values[name]
2347            del self.state.weights[name]
2348            del self.model._persistency_dict[name.split('.')[0]] 
2349            del self.state.model._persistency_dict[name.split('.')[0]]
2350        except:
2351            pass
2352                                           
2353    def _lay_out(self):
2354        """
2355        returns self.Layout
2356       
2357        :Note: Mac seems to like this better when self.
2358            Layout is called after fitting.
2359        """
2360        self._sleep4sec()
2361        self.Layout()
2362        return 
2363   
2364    def _sleep4sec(self):
2365        """
2366            sleep for 1 sec only applied on Mac
2367            Note: This 1sec helps for Mac not to crash on self.:ayout after self._draw_model
2368        """
2369        if ON_MAC == True:
2370            time.sleep(1)
2371           
2372    def _find_polyfunc_selection(self, disp_func = None):
2373        """
2374        FInd Comboox selection from disp_func
2375       
2376        :param disp_function: dispersion distr. function
2377        """
2378        # List of the poly_model name in the combobox
2379        list = ["RectangleDispersion", "ArrayDispersion", 
2380                    "LogNormalDispersion", "GaussianDispersion", 
2381                    "SchulzDispersion"]
2382
2383        # Find the selection
2384        try:
2385            selection = list.index(disp_func.__class__.__name__)
2386            return selection
2387        except:
2388             return 3
2389                           
2390    def on_reset_clicked(self,event):
2391        """
2392        On 'Reset' button  for Q range clicked
2393        """
2394        flag = True
2395        #if self.check_invalid_panel():
2396        #    return
2397        ##For 3 different cases: Data2D, Data1D, and theory
2398        if self.model == None:
2399            msg="Please select a model first..."
2400            wx.MessageBox(msg, 'Info')
2401            flag = False
2402            return
2403           
2404        elif self.data.__class__.__name__ == "Data2D":
2405            data_min= 0
2406            x= max(math.fabs(self.data.xmin), math.fabs(self.data.xmax)) 
2407            y= max(math.fabs(self.data.ymin), math.fabs(self.data.ymax))
2408            self.qmin_x = data_min
2409            self.qmax_x = math.sqrt(x*x + y*y)
2410            # check smearing
2411            if not self.disable_smearer.GetValue():
2412                temp_smearer= self.current_smearer
2413                ## set smearing value whether or not the data contain the smearing info
2414                if self.pinhole_smearer.GetValue():
2415                    flag = self.update_pinhole_smear()
2416                else:
2417                    flag = True
2418                   
2419        elif self.data == None:
2420            self.qmin_x = _QMIN_DEFAULT
2421            self.qmax_x = _QMAX_DEFAULT
2422            self.num_points = _NPTS_DEFAULT           
2423            self.state.npts = self.num_points
2424           
2425        elif self.data.__class__.__name__ != "Data2D":
2426            self.qmin_x = min(self.data.x)
2427            self.qmax_x = max(self.data.x)
2428            # check smearing
2429            if not self.disable_smearer.GetValue():
2430                temp_smearer= self.current_smearer
2431                ## set smearing value whether or not the data contain the smearing info
2432                if self.slit_smearer.GetValue():
2433                    flag = self.update_slit_smear()
2434                elif self.pinhole_smearer.GetValue():
2435                    flag = self.update_pinhole_smear()
2436                else:
2437                    flag = True
2438        else:
2439            flag = False
2440           
2441        if flag == False:
2442            msg= "Cannot Plot :Must enter a number!!!  "
2443            wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
2444        else:
2445            # set relative text ctrs.
2446            self.qmin.SetValue(str(self.qmin_x))
2447            self.qmax.SetValue(str(self.qmax_x))
2448            self.set_npts2fit()
2449            # At this point, some button and variables satatus (disabled?) should be checked
2450            # such as color that should be reset to white in case that it was pink.
2451            self._onparamEnter_helper()
2452
2453        self.save_current_state()
2454        self.state.qmin = self.qmin_x
2455        self.state.qmax = self.qmax_x
2456       
2457        #reset the q range values
2458        self._reset_plotting_range(self.state)
2459        #self.compute_chisqr(smearer=self.current_smearer)
2460        #Re draw plot
2461        self._draw_model()
2462       
2463    def get_images(self):
2464        """
2465        Get the images of the plots corresponding this panel for report
2466       
2467        : return graphs: list of figures
2468        : TODO: Move to guiframe
2469        """
2470        # set list of graphs
2471        graphs = []
2472        canvases = []
2473        # call gui_manager
2474        gui_manager = self.parent.parent
2475        # loops through the panels [dic]
2476        for item1, item2 in gui_manager.panels.iteritems():
2477             data_title = self.data.group_id
2478             data_name = str(self.data.name).split(" [")[0]
2479            # try to get all plots belonging to this control panel
2480             try:
2481                 title = ''
2482                 # check titles (main plot)
2483                 if hasattr(item2,"data2D"):
2484                     title = item2.data2D.title
2485                 # and data_names (model plot[2D], and residuals)
2486                 if item2.group_id == data_title or \
2487                                title.count(data_name) or \
2488                                item2.window_name.count(data_name) or \
2489                                item2.window_caption.count(data_name):
2490                     #panel = gui_manager._mgr.GetPane(item2.window_name)
2491                     # append to the list
2492                     graphs.append(item2.figure) 
2493                     canvases.append(item2.canvas)     
2494             except:
2495                 # Not for control panels
2496                 pass
2497        # return the list of graphs
2498        return graphs, canvases
2499
2500    def on_model_help_clicked(self,event):
2501        """
2502        on 'More details' button
2503        """
2504        from help_panel import  HelpWindow
2505        import sans.models as models 
2506       
2507        # Get models help model_function path
2508        path = models.get_data_path(media='media')
2509        model_path = os.path.join(path,"model_functions.html")
2510        if self.model == None:
2511            name = 'FuncHelp'
2512        else:
2513            name = self.model.__class__.__name__
2514
2515        frame = HelpWindow(None, -1,  pageToOpen=model_path)   
2516        frame.Show(True)
2517        if frame.rhelp.HasAnchor(name):
2518            frame.rhelp.ScrollToAnchor(name)
2519        else:
2520           msg= "Model does not contains an available description "
2521           msg +="Please try searching in the Help window"
2522           wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))     
2523   
2524    def on_pd_help_clicked(self, event):
2525        """
2526        Button event for PD help
2527        """
2528        from help_panel import  HelpWindow
2529        import sans.models as models 
2530       
2531        # Get models help model_function path
2532        path = models.get_data_path(media='media')
2533        pd_path = os.path.join(path,"pd_help.html")
2534
2535        frame = HelpWindow(None, -1,  pageToOpen=pd_path)   
2536        frame.Show(True)
2537       
2538    def on_left_down(self, event):
2539        """
2540        Get key stroke event
2541        """
2542        # Figuring out key combo: Cmd for copy, Alt for paste
2543        if event.CmdDown() and event.ShiftDown():
2544            flag = self.get_paste()
2545        elif event.CmdDown():
2546            flag = self.get_copy()
2547        else:
2548            event.Skip()
2549            return
2550        # make event free
2551        event.Skip()
2552        # messages depending on the flag
2553        if flag == None:
2554            msg = " Parameter values are copied to the clipboard..."
2555            infor = 'warning'
2556        elif flag:
2557            msg = " Parameter values are pasted from the clipboad..."
2558            infor = "warning"
2559        else:
2560            msg = "Error was occured during pasting the parameter values "
2561            msg += "from the clipboard..."
2562            infor = "error"
2563        # inform msg to wx
2564        wx.PostEvent( self.parent.parent, 
2565                      StatusEvent(status= msg, info=infor))
2566       
2567           
2568    def get_copy(self): 
2569        """
2570        Get the string copies of the param names and values in the tap
2571        """ 
2572        content = 'sansview_parameter_values:'
2573        # Do it if params exist       
2574        if  self.parameters !=[]:
2575           
2576            # go through the parameters
2577            string = self._get_copy_helper(self.parameters, 
2578                                           self.orientation_params)
2579            content += string
2580           
2581            # go through the fittables
2582            string = self._get_copy_helper(self.fittable_param, 
2583                                           self.orientation_params_disp)
2584            content += string
2585
2586            # go through the fixed params
2587            string = self._get_copy_helper(self.fixed_param, 
2588                                           self.orientation_params_disp)
2589            content += string
2590               
2591            # go through the str params
2592            string = self._get_copy_helper(self.str_parameters, 
2593                                           self.orientation_params)
2594            content += string
2595
2596        if wx.TheClipboard.Open():
2597            wx.TheClipboard.SetData(wx.TextDataObject(str(content)))
2598            data = wx.TextDataObject()
2599            success = wx.TheClipboard.GetData(data)
2600            text = data.GetText()
2601            wx.TheClipboard.Close()
2602           
2603        return None
2604   
2605    def _get_copy_helper(self, param, orient_param):
2606        """
2607        Helping get value and name of the params
2608       
2609        : param param:  parameters
2610        : param orient_param = oritational params
2611        : return content: strings [list] [name,value:....]
2612        """
2613        content = ''
2614        # go through the str params
2615        for item in param: 
2616            # 2D
2617            if self.data.__class__.__name__== "Data2D":
2618                name = item[1]
2619                value = item[2].GetValue()
2620            # 1D
2621            else:
2622                ## for 1D all parameters except orientation
2623                if not item[1] in orient_param:
2624                    name = item[1]
2625                    value = item[2].GetValue()
2626            # add to the content
2627            content += name + ',' + value + ':'
2628           
2629        return content
2630   
2631    def get_paste(self): 
2632        """
2633        Get the string copies of the param names and values in the tap
2634        """ 
2635        context = {}
2636        text = ""
2637       
2638        # Get text from the clip board       
2639        if wx.TheClipboard.Open():
2640            if wx.TheClipboard.IsSupported(wx.DataFormat(wx.DF_TEXT)):
2641                data = wx.TextDataObject()
2642                # get wx dataobject
2643                success = wx.TheClipboard.GetData(data)
2644                # get text
2645                text = data.GetText()
2646            # close clipboard
2647            wx.TheClipboard.Close()
2648           
2649        # put the text into dictionary   
2650        lines = text.split(':')
2651        if lines[0] != 'sansview_parameter_values':
2652            return False
2653        for line in lines[1:-1]:
2654            if len(line) != 0:
2655                item =line.split(',')
2656                name = item[0]
2657                value = item[1]
2658                # Transfer the text to content[dictionary]
2659                context[name] = value
2660       
2661        # Do it if params exist       
2662        if  self.parameters != []:
2663            # go through the parameters 
2664            self._get_paste_helper(self.parameters, 
2665                                   self.orientation_params, context)
2666
2667            # go through the fittables
2668            self._get_paste_helper(self.fittable_param, 
2669                                   self.orientation_params_disp, context)
2670
2671            # go through the fixed params
2672            self._get_paste_helper(self.fixed_param, 
2673                                   self.orientation_params_disp, context)
2674           
2675            # go through the str params
2676            self._get_paste_helper(self.str_parameters, 
2677                                   self.orientation_params, context)
2678           
2679        return True
2680   
2681    def _get_paste_helper(self, param, orient_param, content):
2682        """
2683        Helping set values of the params
2684       
2685        : param param:  parameters
2686        : param orient_param: oritational params
2687        : param content: dictionary [ name, value: name1.value1,...]
2688        """
2689        # go through the str params
2690        for item in param: 
2691            # 2D
2692            if self.data.__class__.__name__== "Data2D":
2693                name = item[1]
2694                if name in content.keys():
2695                    item[2].SetValue(content[name])
2696            # 1D
2697            else:
2698                ## for 1D all parameters except orientation
2699                if not item[1] in orient_param:
2700                    name = item[1]
2701                    if name in content.keys():
2702                        # Avoid changing combox content which needs special care
2703                        if item[2].__class__.__name__ != "ComboBox":
2704                            item[2].SetValue(content[name])
2705                           
2706               
Note: See TracBrowser for help on using the repository browser.