source: sasview/sansview/perspectives/fitting/basepage.py @ fca9cbd9

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 fca9cbd9 was fca9cbd9, checked in by Gervaise Alina <gervyh@…>, 15 years ago

drawing with smearering

  • Property mode set to 100644
File size: 50.5 KB
Line 
1
2import sys, os
3import wx
4import numpy
5import time
6import copy 
7from sans.guiframe.utils import format_number,check_float
8from sans.guicomm.events import StatusEvent   
9import pagestate
10from pagestate import PageState
11(PageInfoEvent, EVT_PAGE_INFO)   = wx.lib.newevent.NewEvent()
12_BOX_WIDTH = 80
13
14class BasicPage(wx.ScrolledWindow):
15    """
16        This class provide general structure of  fitpanel page
17    """
18     ## Internal name for the AUI manager
19    window_name = "Basic Page"
20    ## Title to appear on top of the window
21    window_caption = "Basic page "
22   
23    def __init__(self,parent, page_info):
24        wx.ScrolledWindow.__init__(self, parent)
25        ##window_name
26        self.window_name = page_info.window_name
27        ##window_caption
28        self.window_caption = page_info.window_caption
29        ## parent of the page
30        self.parent = parent
31        ## manager is the fitting plugin
32        self.manager= page_info.manager
33        ## owner of the page (fitting plugin)
34        self.event_owner= page_info.event_owner
35         ## current model
36        self.model = page_info.model
37        ## data
38        self.data = page_info.data
39        ## dictionary containing list of models
40        self.model_list_box = page_info.model_list_box
41        ## Data member to store the dispersion object created
42        self._disp_obj_dict = {}
43        ## selected parameters to apply dispersion
44        self.disp_cb_dict ={}
45        ## smearer object
46        self.smearer = None
47        ##list of model parameters. each item must have same length
48        ## each item related to a given parameters
49        ##[cb state, name, value, "+/-", error of fit, min, max , units]
50        self.parameters=[]
51        ## list of parameters to fit , must be like self.parameters
52        self.param_toFit=[]
53        ## list of looking like parameters but with non fittable parameters info
54        self.fixed_param=[]
55        ## list of looking like parameters but with  fittable parameters info
56        self.fittable_param=[]
57        ##list of dispersion parameters
58        self.disp_list=[]
59        ## list of orientation parameters
60        self.orientation_params=[]
61        self.orientation_params_disp=[]
62        if self.model !=None:
63            self.disp_list= self.model.getDispParamList()
64        ##enable model 2D draw
65        self.enable2D= False
66        ## check that the fit range is correct to plot the model again
67        self.fitrange= True
68        ## Q range
69        self.qmin_x= 0.001
70        self.qmax_x= 0.13
71        self.num_points= 100
72        ## Create memento to save the current state
73       
74        self.state= PageState(parent= self.parent,model=self.model, data=self.data)
75       
76        ## retrieve saved state
77        self.number_saved_state= 0
78        ## dictionary of saved state
79        self.saved_states={}
80        self.slicerpop = wx.Menu()
81        ## Default locations
82        self._default_save_location = os.getcwd()     
83        ## save initial state on context menu
84        self.onSave(event=None)
85        self.Bind(wx.EVT_CONTEXT_MENU, self.onContextMenu)
86       
87        ## create the basic structure of the panel with empty sizer
88        self.define_page_structure()
89        ## drawing Initial dispersion parameters sizer
90        self.set_dispers_sizer()
91        self._fill_save_sizer()
92        ## layout
93        self.set_layout()
94       
95       
96    def onContextMenu(self, event): 
97        """
98            Retrieve the state selected state
99        """
100        pos = event.GetPosition()
101        pos = self.ScreenToClient(pos)
102        self.PopupMenu(self.slicerpop, pos) 
103       
104       
105       
106    def define_page_structure(self):
107        """
108            Create empty sizer for a panel
109        """
110        self.vbox  = wx.BoxSizer(wx.VERTICAL)
111        self.sizer0 = wx.BoxSizer(wx.VERTICAL)
112        self.sizer1 = wx.BoxSizer(wx.VERTICAL)
113        self.sizer2 = wx.BoxSizer(wx.VERTICAL)
114        self.sizer3 = wx.BoxSizer(wx.VERTICAL)
115        self.sizer4 = wx.BoxSizer(wx.VERTICAL)
116        self.sizer5 = wx.BoxSizer(wx.VERTICAL)
117        self.sizer6 = wx.BoxSizer(wx.VERTICAL)
118       
119        self.sizer0.SetMinSize((375,-1))
120        self.sizer1.SetMinSize((375,-1))
121        self.sizer2.SetMinSize((375,-1))
122        self.sizer3.SetMinSize((375,-1))
123        self.sizer4.SetMinSize((375,-1))
124        self.sizer5.SetMinSize((375,-1))
125        self.sizer6.SetMinSize((375,-1))
126       
127        self.vbox.Add(self.sizer0)
128        self.vbox.Add(self.sizer1)
129        self.vbox.Add(self.sizer2)
130        self.vbox.Add(self.sizer3)
131        self.vbox.Add(self.sizer4)
132        self.vbox.Add(self.sizer5)
133        self.vbox.Add(self.sizer6)
134       
135       
136    def set_layout(self):
137        """
138             layout
139        """
140        self.vbox.Layout()
141        self.vbox.Fit(self) 
142        self.SetSizer(self.vbox)
143       
144        self.set_scroll()
145        self.Centre()
146       
147       
148    def set_scroll(self):
149        self.SetScrollbars(20,20,200,100)
150        self.Layout()   
151        self.SetAutoLayout(True)
152         
153         
154    def set_owner(self,owner):
155        """
156            set owner of fitpage
157            @param owner: the class responsible of plotting
158        """
159        self.event_owner = owner   
160        self.state.event_owner = owner
161 
162    def set_manager(self, manager):
163        """
164             set panel manager
165             @param manager: instance of plugin fitting
166        """
167        self.manager = manager 
168        self.state.manager = manager
169       
170    def populate_box(self, dict):
171        """
172             Store list of model
173             @param dict: dictionary containing list of models
174        """
175        self.model_list_box = dict
176        self.state.model_list_box = self.model_list_box
177           
178   
179       
180    def set_dispers_sizer(self):
181        """
182            fill sizer containing dispersity info
183        """
184        self.sizer4.Clear(True)
185        name="Polydispersity and Orientational Distribution"
186        box_description= wx.StaticBox(self, -1,name)
187        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
188        #----------------------------------------------------
189        self.disable_disp = wx.RadioButton(self, -1, 'Off', (10, 10), style=wx.RB_GROUP)
190        self.enable_disp = wx.RadioButton(self, -1, 'On', (10, 30))
191        ## saving the state of enable dispersity button
192        self.state.enable_disp= self.enable_disp.GetValue()
193       
194        self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param, id=self.disable_disp.GetId())
195        self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param, id=self.enable_disp.GetId())
196       
197        sizer_dispersion = wx.BoxSizer(wx.HORIZONTAL)
198        sizer_dispersion.Add((20,20))
199        name=""#Polydispersity and \nOrientational Distribution "
200        sizer_dispersion.Add(wx.StaticText(self,-1,name))
201        sizer_dispersion.Add(self.enable_disp )
202        sizer_dispersion.Add((20,20))
203        sizer_dispersion.Add(self.disable_disp )
204        sizer_dispersion.Add((10,10))
205       
206        ## fill a sizer with the combobox to select dispersion type
207        sizer_select_dispers = wx.BoxSizer(wx.HORIZONTAL) 
208        self.model_disp = wx.StaticText(self, -1, 'Distribution Function ')
209           
210        import sans.models.dispersion_models 
211        self.polydisp= sans.models.dispersion_models.models
212        self.disp_box = wx.ComboBox(self, -1)
213       
214        for key in self.polydisp.iterkeys():
215            name = str(key.__name__)
216            if name=="ArrayDispersion":
217                # Remove the option until the rest of the code is ready for it
218                self.disp_box.Append("Select customized Model",key)
219               
220            else:
221                self.disp_box.Append(name,key)
222        self.disp_box.SetSelection(0) 
223       
224        wx.EVT_COMBOBOX(self.disp_box,-1, self._on_select_Disp) 
225             
226        sizer_select_dispers.Add((10,10)) 
227        sizer_select_dispers.Add(self.model_disp) 
228        sizer_select_dispers.Add(self.disp_box,0,
229                wx.TOP|wx.BOTTOM|wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE,border=5)
230        #sizer_select_dispers.Add((10,10))
231        self.model_disp.Hide()
232        self.disp_box.Hide()
233       
234        boxsizer1.Add( sizer_dispersion,0,
235                wx.TOP|wx.BOTTOM|wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE,border=5)
236        #boxsizer1.Add( (10,10) )
237        boxsizer1.Add( sizer_select_dispers )
238        self.sizer4_4 = wx.GridBagSizer(5,5)
239        boxsizer1.Add( self.sizer4_4  )
240        #-----------------------------------------------------
241        self.sizer4.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
242        self.sizer4_4.Layout()
243        self.sizer4.Layout()
244        self.Layout()
245        self.SetScrollbars(20,20,200,100)
246        self.Refresh()
247       
248   
249    def select_disp_angle(self, event): 
250        """
251            Event for when a user select a parameter to average over.
252            @param event: check box event
253        """
254       
255       
256        # Go through the list of dispersion check boxes to identify which one has changed
257        for p in self.disp_cb_dict:
258            # Catch which one of the box was just checked or unchecked.
259            if event.GetEventObject() == self.disp_cb_dict[p]:             
260
261               
262                if self.disp_cb_dict[p].GetValue() == True:
263                    # The user wants this parameter to be averaged.
264                    # Pop up the file selection dialog.
265                    path = self._selectDlg()
266                   
267                    # If nothing was selected, just return
268                    if path is None:
269                        self.disp_cb_dict[p].SetValue(False)
270                        return
271                    try:
272                        self._default_save_location = os.path.dirname(path)
273                    except:
274                        pass 
275                    try:
276                        values,weights = self.read_file(path)
277                    except:
278                        msg="Could not read input file"
279                        wx.PostEvent(self.parent.parent, StatusEvent(status= msg))
280                        return
281                   
282                    # If any of the two arrays is empty, notify the user that we won't
283                    # proceed
284                    if values is None or weights is None:
285                        wx.PostEvent(self.parent.parent, StatusEvent(status=\
286                            "The loaded %s distrubtion is corrupted or empty" % p))
287                        return
288                       
289                    # Tell the user that we are about to apply the distribution
290                    wx.PostEvent(self.parent.parent, StatusEvent(status=\
291                            "Applying loaded %s distribution: %s" % (p, path))) 
292                   
293                    # Create the dispersion objects
294                    from sans.models.dispersion_models import ArrayDispersion
295                    disp_model = ArrayDispersion()
296                    disp_model.set_weights(values, weights)
297                    # Store the object to make it persist outside the scope of this method
298                    #TODO: refactor model to clean this up?
299                    self._disp_obj_dict[p] = disp_model
300                    # Set the new model as the dispersion object for the selected parameter
301                    self.model.set_dispersion(p, disp_model)
302                         
303                else:
304                    # The parameter was un-selected. Go back to Gaussian model (with 0 pts)
305                    self._reset_dispersity()
306                   
307                ## Redraw the model
308                self._draw_model()
309        return
310   
311   
312    def onResetModel(self, event):
313        """
314            Reset model state
315        """
316        ## post help message for the selected model
317        msg = self.slicerpop.GetHelpString(event.GetId())
318        msg +=" reloaded"
319        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
320       
321        name= self.slicerpop.GetLabel(event.GetId())
322        if name in self.saved_states.keys():
323            previous_state = self.saved_states[name]
324            self.reset_page(previous_state)
325           
326   
327    def onSave(self, event):
328        """
329            save history of the data and model
330        """
331        if self.model==None:
332            return 
333        if hasattr(self,"enable_disp"):
334            self.state.enable_disp = self.enable_disp.GetValue()
335        if hasattr(self, "disp_box"):
336            self.state.disp_box = self.disp_box.GetSelection()
337       
338        self.state.model = self.model.clone()
339        new_state = self.state.clone()
340       
341        ##Add model state on context menu
342        self.number_saved_state += 1
343        name= self.model.name+"[%g]"%self.number_saved_state
344        self.saved_states[name]= new_state
345       
346        ## Add item in the context menu
347        year, month, day,hour,minute,second,tda,ty,tm_isdst= time.gmtime()
348        my_time= str(hour)+"hrs "+str(minute)+"min "+str(second)+"s"
349        date= str( month)+"|"+str(day)+"|"+str(year)
350        msg=  "Model saved at %s on %s"%(my_time, date)
351         ## post help message for the selected model
352        msg +=" Saved! right click on this page to retrieve this model"
353        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
354       
355        id = wx.NewId()
356        self.slicerpop.Append(id,name,str(msg))
357        wx.EVT_MENU(self, id, self.onResetModel)
358       
359    def onSetFocus(self, evt):
360        """
361            highlight the current textcrtl and hide the error text control shown
362            after fitting
363        """
364       
365        if hasattr(self,"text2_3"):
366            self.text2_3.Hide()
367        if len(self.parameters)>0:
368            for item in self.parameters:
369                ## hide statictext +/-   
370                if item[3]!=None:
371                    item[3].Hide()
372                ## hide textcrtl  for error after fit
373                if item[4]!=None:
374                    item[4].Clear()
375                    item[4].Hide()
376        if len(self.fittable_param)>0:
377            for item in self.fittable_param:
378                ## hide statictext +/-   
379                if item[3]!=None:
380                    item[3].Hide()
381                ## hide textcrtl  for error after fit
382                if item[4]!=None:
383                    item[4].Clear()
384                    item[4].Hide()
385        self.Layout()
386        # Get a handle to the TextCtrl
387        widget = evt.GetEventObject()
388        # Select the whole control, after this event resolves
389        wx.CallAfter(widget.SetSelection, -1,-1)
390        return
391   
392   
393    def read_file(self, path):
394        """
395            Read two columns file
396            @param path: the path to the file to read
397        """
398        try:
399            if path==None:
400                wx.PostEvent(self.parent.parent, StatusEvent(status=\
401                            " Selected Distribution was not loaded: %s"%path))
402                return None, None
403            input_f = open(path, 'r')
404            buff = input_f.read()
405            lines = buff.split('\n')
406           
407            angles = []
408            weights=[]
409            for line in lines:
410                toks = line.split()
411                if len(toks)==2:
412                    try:
413                        angle = float(toks[0])
414                        weight = float(toks[1])
415                    except:
416                        # Skip non-data lines
417                        pass
418                    angles.append(angle)
419                    weights.append(weight)
420            return numpy.array(angles), numpy.array(weights)
421        except:
422            raise 
423   
424   
425    def createMemento(self):
426        """
427            return the current state of the page
428        """
429        return self.state.clone()
430   
431   
432    def save_current_state(self):
433        """
434            Store current state
435        """
436        if self.model!= None:
437            self.state.model = self.model.clone()
438        self.state.save_data(self.data)
439   
440        if hasattr(self,"cb1"):
441            self.state.cb1= self.cb1.GetValue()
442           
443        if hasattr(self,"enable_disp"):
444            self.state.enable_disp= self.enable_disp.GetValue()
445       
446        self.state.smearer = copy.deepcopy(self.smearer)
447        if hasattr(self,"enable_smearer"):
448            self.state.enable_smearer = copy.deepcopy(self.enable_smearer.GetValue())
449           
450        if hasattr(self,"disp_box"):
451            self.state.disp_box = self.disp_box.GetCurrentSelection()
452        self._save_plotting_range()
453     
454        ## save checkbutton state and txtcrtl values
455        self.state.parameters=[]
456        self.state.fittable_param=[]
457        self.state.fixed_param=[]
458       
459        self._copy_parameters_state(self.parameters, self.state.parameters)
460        self._copy_parameters_state(self.fittable_param, self.state.fittable_param)
461        self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
462       
463        ## post state to fit panel
464        event = PageInfoEvent(page = self)
465        wx.PostEvent(self.parent, event)
466   
467   
468    def reset_page_helper(self, state):
469        """
470            Use page_state and change the state of existing page
471        """
472        self.state = state.clone()
473        self.model= self.state.model
474        self.data = self.state.data
475        self.smearer= self.state.smearer
476        self.state.enable_smearer=state.enable_smearer
477       
478        ##model parameter values restore
479        self._set_model_sizer_selection( self.model )
480        self.set_model_param_sizer(self.model)
481        if hasattr(self, "cb1"):   
482            self.cb1.SetValue(self.state.cb1)
483        ## display dispersion info layer
484        self.enable_disp.SetValue(self.state.enable_disp)
485        if hasattr(self, "disp_box"):
486            self.disp_box.SetSelection(self.state.disp_box) 
487        self._set_dipers_Param(event=None)
488        ##plotting range restore   
489        self._reset_plotting_range()
490        ## smearing info  restore
491        if hasattr(self,"enable_smearer"):
492            ## set smearing value whether or not the data contain the smearing info
493            self.enable_smearer.SetValue(state.enable_smearer)
494            self.compute_chisqr(smearer= self.smearer) 
495           
496        ## reset context menu items
497        self._reset_context_menu()
498        ## reset state of checkbox,textcrtl  and parameters value
499        self._reset_parameters_state(self.parameters,state.parameters)
500        self._reset_parameters_state(self.fittable_param,state.fittable_param)
501        self._reset_parameters_state(self.fixed_param,state.fixed_param)
502       
503        ## draw the model with previous parameters value
504        self._draw_model()
505       
506   
507   
508         
509       
510    def _selectDlg(self):
511        """
512            open a dialog file to selected the customized dispersity
513        """
514        import os
515        dlg = wx.FileDialog(self, "Choose a weight file",
516                                self._default_save_location , "", "*.*", wx.OPEN)
517        path = None
518        if dlg.ShowModal() == wx.ID_OK:
519            path = dlg.GetPath()
520        dlg.Destroy()
521        return path
522   
523   
524    def _reset_context_menu(self):
525        """
526            reset the context menu
527        """
528        for name, state in self.state.saved_states.iteritems():
529            self.number_saved_state += 1
530            ## Add item in the context menu
531            id = wx.NewId()
532            self.slicerpop.Append(id,name, 'Save model and state %g'%self.number_saved_state)
533            wx.EVT_MENU(self, id, self.onResetModel)
534   
535    def _reset_plotting_range(self):
536        """
537            Reset the plotting range to a given state
538        """
539       
540        #self.qmin.SetValue(format_number(self.state.qmin))
541        #self.qmax.SetValue(format_number(self.state.qmax))
542        self.qmin.SetValue(str(self.state.qmin))
543        self.qmax.SetValue(str(self.state.qmax)) 
544        if self.state.npts!=None:
545            self.npts.SetValue(format_number(self.state.npts)) 
546            self.num_points = float(self.state.npts)
547           
548        self.qmin_x = float(self.qmin.GetValue())
549        self.qmax_x = float(self.qmax.GetValue())
550       
551       
552    def _save_plotting_range(self ):
553        """
554            save the state of plotting range
555        """
556        self.state.qmin = self.qmin_x
557        self.state.qmax = self.qmax_x
558        if self.npts!=None:
559            self.state.npts= self.num_points
560           
561           
562    def _onparamEnter_helper(self):
563        """
564             check if values entered by the user are changed and valid to replot
565             model
566             use : _check_value_enter
567        """
568        if self.model !=None:
569            ## save current state
570            self.save_current_state()
571            # Flag to register when a parameter has changed.
572            is_modified = False
573            is_modified =self._check_value_enter( self.fittable_param ,is_modified)
574            is_modified =self._check_value_enter( self.fixed_param ,is_modified)
575            is_modified =self._check_value_enter( self.parameters ,is_modified)       
576           
577           
578            self.Layout()
579            # Here we should check whether the boundaries have been modified.
580            # If qmin and qmax have been modified, update qmin and qmax and
581            # set the is_modified flag to True
582            from sans.guiframe.utils import check_value
583            if check_value( self.qmin, self.qmax):
584                if float(self.qmin.GetValue()) != self.qmin_x:
585                    self.qmin_x = float(self.qmin.GetValue())
586                    is_modified = True
587                if float(self.qmax.GetValue()) != self.qmax_x:
588                    self.qmax_x = float(self.qmax.GetValue())
589                    is_modified = True
590                self.fitrange = True
591            else:
592                self.fitrange = False
593            if self.npts != None:
594                if check_float(self.npts):
595                    if float(self.npts.GetValue()) !=  self.num_points:
596                        self.num_points = float(self.npts.GetValue())
597                        is_modified = True
598                else:
599                    msg= "Cannot Plot :Must enter a number!!!  "
600                    wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
601                   
602           
603            ## if any value is modify draw model with new value
604            if is_modified:
605                self._draw_model() 
606               
607               
608    def _reset_parameters_state(self, listtorestore,statelist):
609        """
610            Reset the parameters at the given state
611        """
612        if len(statelist)==0 or  len(listtorestore)==0 :
613            return
614        if len(statelist)!=  len(listtorestore) :
615            return
616        for j in range(len(listtorestore)):
617            item_page = listtorestore[j]
618            item_page_info = statelist[j]
619            ##change the state of the check box for simple parameters
620            if item_page[0]!=None:
621                item_page[0].SetValue(item_page_info[0])
622   
623            if item_page[2]!=None:
624                item_page[2].SetValue(item_page_info[2])
625               
626            if item_page[3]!=None:
627                ## show or hide text +/-
628                if item_page_info[2]:
629                    item_page[3].Show(True)
630                else:
631                    item_page[3].Hide()
632            if item_page[4]!=None:
633                ## show of hide the text crtl for fitting error
634                if item_page_info[4][0]:
635                    item_page[4].Show(True)
636                    item_page[4].SetValue(item_page_info[4][1])
637                else:
638                    item_page[3].Hide()
639            if item_page[5]!=None:
640                ## show of hide the text crtl for fitting error
641                if item_page_info[5][0]:
642                    item_page[5].Show(True)
643                    item_page[5].SetValue(item_page_info[4][1])
644                else:
645                    item_page[5].Hide()
646                   
647            if item_page[6]!=None:
648                ## show of hide the text crtl for fitting error
649                if item_page_info[6][0]:
650                    item_page[6].Show(True)
651                    item_page[6].SetValue(item_page_info[6][1])
652                else:
653                    item_page[6].Hide()
654                           
655                           
656    def _copy_parameters_state(self, listtocopy, statelist):
657        """
658            copy the state of button
659            @param listtocopy: the list of check button to copy
660            @param statelist: list of state object to store the current state
661        """
662        if len(listtocopy)==0:
663            return
664       
665        for item in listtocopy:
666            checkbox_state = None
667            if item[0]!= None:
668                checkbox_state= item[0].GetValue()
669            parameter_name = item[1]
670            parameter_value = None
671            if item[2]!=None:
672                parameter_value = item[2].GetValue()
673            static_text = None
674            if item[3]!=None:
675                static_text = item[3].IsShown()
676            error_value = None
677            error_state = None
678            if item[4]!= None:
679                error_value = item[4].GetValue()
680                error_state = item[4].IsShown()
681               
682            min_value = None
683            min_state = None
684            if item[5]!= None:
685                min_value = item[5].GetValue()
686                min_state = item[5].IsShown()
687               
688            max_value = None
689            max_state = None
690            if item[6]!= None:
691                max_value = item[6].GetValue()
692                max_state = item[6].IsShown()
693               
694            statelist.append([checkbox_state, parameter_name, parameter_value,
695                              static_text ,[error_state,error_value],
696                                [min_state,min_value],[max_state , max_value],None])
697           
698       
699   
700    def _set_model_sizer_selection(self, model):
701        """
702            Display the sizer according to the type of the current model
703        """
704        if hasattr(model ,"model2"):
705           
706            class_name= model.model2.__class__
707            name= model.model2.name
708            flag= name != "NoStructure"
709            if flag and (class_name in self.model_list_box["Structure Factors"]):
710                self.structurebox.Enable()
711                items = self.structurebox.GetItems()
712                self.sizer1.Layout()
713                self.SetScrollbars(20,20,200,100)
714                for i in range(len(items)):
715                    if items[i]== str(name):
716                        self.structurebox.SetSelection(i)
717                        break
718                   
719        if hasattr(model ,"model1"):
720            class_name = model.model1.__class__
721            name = model.model1.name
722            self.formfactorbox.Clear()
723           
724            for k, list in self.model_list_box.iteritems():
725                if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
726                    self.shape_rbutton.SetValue(True)
727                    ## fill the form factor list with new model
728                    self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
729                    items = self.formfactorbox.GetItems()
730                    ## set comboxbox to the selected item
731                    for i in range(len(items)):
732                        if items[i]== str(name):
733                            self.formfactorbox.SetSelection(i)
734                            break
735                    return
736                elif k == "Shape-Independent":
737                    self.shape_indep_rbutton.SetValue(True)
738                elif k == "Structure Factors":
739                     self.struct_rbutton.SetValue(True)
740                else:
741                    self.plugin_rbutton.SetValue(True)
742               
743                if class_name in list:
744                    ## fill the form factor list with new model
745                    self._populate_box(self.formfactorbox, list)
746                    items = self.formfactorbox.GetItems()
747                    ## set comboxbox to the selected item
748                    for i in range(len(items)):
749                        if items[i]== str(name):
750                            self.formfactorbox.SetSelection(i)
751                            break
752                    break
753        else:
754            ## Select the model from the combobox
755            class_name = model.__class__
756            name = model.name
757            self.formfactorbox.Clear()
758            items = self.formfactorbox.GetItems()
759   
760            for k, list in self.model_list_box.iteritems():
761                if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
762                   
763                    if class_name in self.model_list_box["P(Q)*S(Q)"]:
764                        self.structurebox.Enable()
765                    else:
766                        self.structurebox.Disable()
767                        self.structurebox.SetSelection(0)
768                       
769                    self.shape_rbutton.SetValue(True)
770                    ## fill the form factor list with new model
771                    self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
772                    items = self.formfactorbox.GetItems()
773                    ## set comboxbox to the selected item
774                    for i in range(len(items)):
775                        if items[i]== str(name):
776                            self.formfactorbox.SetSelection(i)
777                            break
778                    return
779                elif k == "Shape-Independent":
780                    self.shape_indep_rbutton.SetValue(True)
781                elif k == "Structure Factors":
782                     self.struct_rbutton.SetValue(True)
783                else:
784                    self.plugin_rbutton.SetValue(True)
785                if class_name in list:
786                    self.structurebox.SetSelection(0)
787                    self.structurebox.Disable()
788                    ## fill the form factor list with new model
789                    self._populate_box(self.formfactorbox, list)
790                    items = self.formfactorbox.GetItems()
791                    ## set comboxbox to the selected item
792                    for i in range(len(items)):
793                        if items[i]== str(name):
794                            self.formfactorbox.SetSelection(i)
795                            break
796                    break
797                   
798       
799    def _draw_model(self):
800        """
801            Method to draw or refresh a plotted model.
802            The method will use the data member from the model page
803            to build a call to the fitting perspective manager.
804           
805            [Note to coder: This way future changes will be done in only one place.]
806        """
807        if self.model !=None:
808            temp_smear=None
809            if hasattr(self, "enable_smearer"):
810                if self.enable_smearer.GetValue():
811                    temp_smear= self.smearer
812            self.manager.draw_model(self.model, data=self.data,
813                                    smearer= temp_smearer,
814                                    qmin=float(self.qmin_x), qmax=float(self.qmax_x),
815                                    qstep= float(self.num_points),
816                                    enable2D=self.enable2D) 
817       
818    def _set_model_sizer(self, sizer, title="", object=None):
819        """
820            Use lists to fill a sizer for model info
821        """
822       
823        sizer.Clear(True)
824        box_description= wx.StaticBox(self, -1,str(title))
825        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
826        #--------------------------------------------------------
827        self.shape_rbutton = wx.RadioButton(self, -1, 'Shapes', style=wx.RB_GROUP)
828        self.shape_indep_rbutton = wx.RadioButton(self, -1, "Shape-Independent")
829        self.struct_rbutton = wx.RadioButton(self, -1, "Structure Factor ")
830        self.plugin_rbutton = wx.RadioButton(self, -1, "Customized Models")
831       
832        self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
833                            id= self.shape_rbutton.GetId() ) 
834        self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
835                            id= self.shape_indep_rbutton.GetId() ) 
836        self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
837                            id= self.struct_rbutton.GetId() ) 
838        self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
839                            id= self.plugin_rbutton.GetId() ) 
840       
841     
842        sizer_radiobutton = wx.GridSizer(2, 2,5, 5)
843        sizer_radiobutton.Add(self.shape_rbutton)
844        sizer_radiobutton.Add(self.shape_indep_rbutton)
845        sizer_radiobutton.Add(self.plugin_rbutton)
846        sizer_radiobutton.Add(self.struct_rbutton)
847       
848        sizer_selection = wx.BoxSizer(wx.HORIZONTAL)
849       
850        self.text1 = wx.StaticText( self,-1,"" )
851        self.text2 = wx.StaticText( self,-1,"P(Q)*S(Q)" )
852       
853       
854        self.formfactorbox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
855        if self.model!=None:
856            self.formfactorbox.SetValue(self.model.name)
857           
858           
859        self.structurebox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
860        wx.EVT_COMBOBOX(self.formfactorbox,-1, self._on_select_model)
861        wx.EVT_COMBOBOX(self.structurebox,-1, self._on_select_model)
862       
863       
864       
865        ## fill combox box
866        if len(self.model_list_box)>0:
867            self._populate_box( self.formfactorbox,self.model_list_box["Shapes"])
868       
869        if len(self.model_list_box)>0:
870            self._populate_box( self.structurebox,
871                                self.model_list_box["Structure Factors"])
872            self.structurebox.Insert("None", 0,None)
873            self.structurebox.SetSelection(0)
874            self.structurebox.Disable()
875 
876            if self.model.__class__ in self.model_list_box["P(Q)*S(Q)"]:
877                self.structurebox.Enable()
878           
879       
880        ## check model type to show sizer
881        if self.model !=None:
882            self._set_model_sizer_selection( self.model )
883       
884        sizer_selection.Add(self.text1)
885        sizer_selection.Add((5,5))
886        sizer_selection.Add(self.formfactorbox)
887        sizer_selection.Add((5,5))
888        sizer_selection.Add(self.text2)
889        sizer_selection.Add((5,5))
890        sizer_selection.Add(self.structurebox)
891        sizer_selection.Add((5,5))
892       
893        boxsizer1.Add( sizer_radiobutton )
894        boxsizer1.Add( (20,20))
895        boxsizer1.Add( sizer_selection )
896        if object !=None:
897            boxsizer1.Add( (-72,-72))
898            boxsizer1.Add( object,  0, wx.ALIGN_RIGHT| wx.RIGHT, 10)
899            boxsizer1.Add( (60,60))
900        #--------------------------------------------------------
901        sizer.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
902        sizer.Layout()
903        self.SetScrollbars(20,20,200,100)
904       
905       
906    def _show_combox(self, event):
907        """
908            Show combox box associate with type of model selected
909        """
910        ## Don't want to populate combo box again if the event comes from check box
911        if self.shape_rbutton.GetValue()and\
912             event.GetEventObject()==self.shape_rbutton:
913            ##fill the combobox with form factor list
914            self.structurebox.SetSelection(0)
915            self.structurebox.Disable()
916            self.formfactorbox.Clear()
917            self._populate_box( self.formfactorbox,self.model_list_box["Shapes"])
918           
919        if self.shape_indep_rbutton.GetValue()and\
920             event.GetEventObject()==self.shape_indep_rbutton:
921            ##fill the combobox with shape independent  factor list
922            self.structurebox.SetSelection(0)
923            self.structurebox.Disable()
924            self.formfactorbox.Clear()
925            self._populate_box( self.formfactorbox,
926                                self.model_list_box["Shape-Independent"])
927           
928        if self.struct_rbutton.GetValue() and\
929             event.GetEventObject()==self.struct_rbutton:
930            ##fill the combobox with structure factor list
931            self.structurebox.SetSelection(0)
932            self.structurebox.Disable()
933            self.formfactorbox.Clear()
934            self._populate_box( self.formfactorbox,
935                                self.model_list_box["Structure Factors"])
936           
937        if self.plugin_rbutton.GetValue()and\
938             event.GetEventObject()==self.plugin_rbutton:
939           
940            ##fill the combobox with form factor list
941            self.structurebox.Disable()
942            self.formfactorbox.Clear()
943            self._populate_box( self.formfactorbox,
944                                self.model_list_box["Customized Models"])
945       
946        self._on_select_model(event=None)
947        self.sizer4_4.Layout()
948        self.sizer4.Layout()
949        self.Layout()
950        self.Refresh()
951        self.SetScrollbars(20,20,200,100)
952           
953    def _populate_box(self, combobox, list):
954        """
955            fill combox box with dict item
956            @param list: contains item to fill the combox
957            item must model class
958        """
959        for models in list:
960            model= models()
961            name = model.__class__.__name__
962            if models.__name__!="NoStructure":
963                if hasattr(model, "name"):
964                    name = model.name
965                combobox.Append(name,models)
966        try:
967
968            combobox.SetSelection(0)
969           
970        except:
971            pass
972   
973        return 0
974   
975   
976    def _on_select_model_helper(self): 
977        """
978             call back for model selection
979        """
980        ## reset dictionary containing reference to dispersion
981        self._disp_obj_dict = {}
982        self.disp_cb_dict ={}
983        f_id = self.formfactorbox.GetCurrentSelection()
984        form_factor = self.formfactorbox.GetClientData( f_id )
985        if not form_factor in  self.model_list_box["multiplication"]:
986            self.structurebox.Disable()
987            self.structurebox.SetSelection(0)
988        else:
989            self.structurebox.Enable()
990           
991        s_id = self.structurebox.GetCurrentSelection()
992        struct_factor = self.structurebox.GetClientData( s_id )
993       
994        if  struct_factor !=None:
995            from sans.models.MultiplicationModel import MultiplicationModel
996            self.model= MultiplicationModel(form_factor(),struct_factor())
997        else:
998            self.model= form_factor()
999       
1000        ## post state to fit panel
1001        self.save_current_state()
1002        self.sizer4_4.Layout()
1003        self.sizer4.Layout()
1004        self.Layout()
1005        self.SetScrollbars(20,20,200,100)
1006        self.Refresh()
1007       
1008       
1009       
1010    def _onparamEnter(self,event):
1011        """
1012            when enter value on panel redraw model according to changed
1013        """
1014        tcrtl= event.GetEventObject()
1015        if check_float(tcrtl):
1016            self._onparamEnter_helper()
1017        else:
1018            msg= "Cannot Plot :Must enter a number!!!  "
1019            wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1020            return 
1021       
1022       
1023    def _check_value_enter(self, list, modified):
1024        """
1025            @param list: model parameter and panel info
1026            each item of the list should be as follow:
1027            item=[cb state, name, value, "+/-", error of fit, min, max , units]
1028        """ 
1029        is_modified =  modified
1030        if len(list)==0:
1031            return is_modified
1032        for item in list:
1033            try:
1034                name = str(item[1])
1035                if hasattr(self,"text2_3"):
1036                    self.text2_3.Hide()
1037                ## check model parameters range
1038                ## check minimun value
1039                param_min= None
1040                param_max= None
1041                if item[5]!= None:
1042                    if item[5].GetValue().lstrip().rstrip()!="":
1043                        param_min = float(item[5].GetValue())
1044                   
1045                ## check maximum value
1046                if item[6]!= None:
1047                    if item[6].GetValue().lstrip().rstrip()!="":
1048                        param_max = float(item[6].GetValue())
1049                       
1050                from sans.guiframe.utils import check_value
1051                if param_min != None and param_max !=None:
1052                    if not check_value(item[5], item[6]):
1053                        msg= "Wrong Fit range entered for parameter "
1054                        msg+= "name %s of model %s "%(name, self.model.name)
1055                        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1056                if name in self.model.details.keys():   
1057                    self.model.details[name][1:]= param_min,param_max
1058               
1059                ## hide statictext +/-   
1060                if item[3]!=None:
1061                    item[3].Hide()
1062                ## hide textcrtl  for error after fit
1063                if item[4]!=None:
1064                    item[4].Clear()
1065                    item[4].Hide()
1066                   
1067           
1068                value= float(item[2].GetValue())
1069               
1070                # If the value of the parameter has changed,
1071                # +update the model and set the is_modified flag
1072                if value != self.model.getParam(name):
1073                    self.model.setParam(name,value)
1074                    is_modified = True   
1075           
1076            except:
1077                msg= "Model Drawing  Error:wrong value entered : %s"% sys.exc_value
1078                wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1079                return 
1080       
1081        return is_modified
1082       
1083 
1084    def _set_dipers_Param(self, event):
1085        """
1086            Add more item to select user dispersity
1087        """
1088        if self.model == None:
1089            msg= " Select non - model value:%s !"%self.model
1090            wx.PostEvent(self.parent.parent, StatusEvent(status= msg))
1091            return 
1092        else:
1093           
1094            if self.enable_disp.GetValue():
1095                self.model_disp.Show(True)
1096                self.disp_box.Show(True)
1097                ## layout for model containing no dispersity parameters
1098                if len(self.disp_list)==0:
1099                    self._layout_sizer_noDipers() 
1100                else:
1101                    ## set gaussian sizer
1102                    self._on_select_Disp(event=None)
1103            else:
1104                self.model_disp.Hide()
1105                self.disp_box.Hide()
1106                self.sizer4_4.Clear(True)
1107                self._reset_dispersity()
1108               
1109         
1110            ## post state to fit panel
1111            self.save_current_state()
1112           
1113         
1114           
1115       
1116    def _layout_sizer_noDipers(self):
1117        """
1118            Draw a sizer with no dispersity info
1119        """
1120        ix=0
1121        iy=1
1122        self.fittable_param=[]
1123        self.fixed_param=[]
1124        self.model_disp.Hide()
1125        self.disp_box.Hide()
1126        self.sizer4_4.Clear(True)
1127        model_disp = wx.StaticText(self, -1, 'No PolyDispersity for this model')
1128        self.sizer4_4.Add(model_disp,( iy, ix),(1,1),  wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
1129        self.sizer4_4.Layout()
1130        self.sizer4.Layout()
1131        self.SetScrollbars(20,20,200,100)
1132     
1133           
1134    def _reset_dispersity(self):
1135        """
1136             put gaussian dispersity into current model
1137        """
1138        self.fittable_param=[]
1139        self.fixed_param=[]
1140       
1141        from sans.models.dispersion_models import GaussianDispersion
1142        if len(self.disp_cb_dict)==0:
1143            self.sizer4_4.Clear(True)
1144            self.sizer4_4.Layout()
1145            self.sizer4.Layout()
1146            self.Layout()
1147            self.Refresh()
1148            self.SetScrollbars(20,20,200,100)   
1149            return 
1150       
1151        for p in self.disp_cb_dict:
1152            # The parameter was un-selected. Go back to Gaussian model (with 0 pts)
1153            disp_model = GaussianDispersion()
1154            # Store the object to make it persist outside the scope of this method
1155            #TODO: refactor model to clean this up?
1156            self._disp_obj_dict[p] = disp_model
1157            # Set the new model as the dispersion object for the selected parameter
1158            self.model.set_dispersion(p, disp_model)
1159            # Redraw the model
1160            self._draw_model()
1161           
1162        self.sizer4_4.Layout()
1163        self.sizer4.Layout()
1164        self.Layout()
1165        self.SetScrollbars(20,20,200,100)   
1166        self.Refresh()
1167       
1168           
1169           
1170    def _on_select_Disp(self,event):
1171        """
1172             allow selecting different dispersion
1173             self.disp_list should change type later .now only gaussian
1174        """
1175       
1176        n = self.disp_box.GetCurrentSelection()
1177        dispersity= self.disp_box.GetClientData(n)
1178        name= dispersity.__name__
1179        if name == "GaussianDispersion":
1180            self._set_sizer_gaussian()
1181           
1182        if  name=="ArrayDispersion":
1183            self._set_sizer_arraydispersion()
1184           
1185        self.state.disp_box= n
1186        ## post state to fit panel
1187        event = PageInfoEvent(page = self)
1188        wx.PostEvent(self.parent, event)
1189       
1190        self.sizer4_4.Layout()
1191        self.sizer4.Layout()
1192        self.Layout()
1193        self.SetScrollbars(20,20,200,100)
1194        self.Refresh()
1195       
1196       
1197       
1198       
1199    def _set_sizer_arraydispersion(self):
1200        """
1201            draw sizer with array dispersity  parameters
1202        """
1203        self.orientation_params_disp=[]
1204        self.sizer4_4.Clear(True) 
1205        ix=0
1206        iy=1     
1207        disp1 = wx.StaticText(self, -1, 'Array Dispersion')
1208        self.sizer4_4.Add(disp1,( iy, ix),(1,1),  wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
1209       
1210        # Look for model parameters to which we can apply an ArrayDispersion model
1211        # Add a check box for each parameter.
1212        self.disp_cb_dict = {}
1213        for p in self.model.dispersion.keys():
1214            ix+=1 
1215            self.disp_cb_dict[p] = wx.CheckBox(self, -1, p, (10, 10))
1216           
1217            wx.EVT_CHECKBOX(self, self.disp_cb_dict[p].GetId(), self.select_disp_angle)
1218            self.sizer4_4.Add(self.disp_cb_dict[p], (iy, ix), (1,1), wx.EXPAND|wx.ADJUST_MINSIZE, 0)
1219       
1220        ix =0
1221        iy +=1 
1222        self.sizer4_4.Add((20,20),(iy,ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)       
1223        self.sizer4_4.Layout()
1224        self.sizer4.Layout()
1225        self.SetScrollbars(20,20,200,100)
1226       
1227       
1228    def _set_range_sizer(self, title, object1=None,object=None):
1229        """
1230            Fill the
1231        """
1232        self.sizer5.Clear(True)
1233        box_description= wx.StaticBox(self, -1,str(title))
1234        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
1235        #--------------------------------------------------------------
1236        self.qmin    = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
1237        #For qmin and qmax, do not use format_number.(If do, qmin and max could be different from what is in the data.)
1238        #self.qmin.SetValue(format_number(self.qmin_x))
1239        self.qmin.SetValue(str(self.qmin_x))
1240        self.qmin.SetToolTipString("Minimun value of Q in linear scale.")
1241        self.qmin.Bind(wx.EVT_SET_FOCUS, self.onSetFocus)
1242        self.qmin.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
1243        self.qmin.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
1244     
1245        self.qmax    = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
1246        #For qmin and qmax, do not use format_number.(If do, qmin and max could be different from what is in the data.)
1247        #self.qmax.SetValue(format_number(self.qmax_x))
1248        self.qmax.SetValue(str(self.qmax_x))
1249        self.qmax.SetToolTipString("Maximum value of Q in linear scale.")
1250        self.qmax.Bind(wx.EVT_SET_FOCUS, self.onSetFocus)
1251        self.qmax.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
1252        self.qmax.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
1253     
1254        sizer_horizontal=wx.BoxSizer(wx.HORIZONTAL)
1255        sizer= wx.GridSizer(3, 3,5, 5)
1256       
1257        sizer.Add((5,5))
1258        sizer.Add(wx.StaticText(self, -1, 'Min'))
1259        sizer.Add(wx.StaticText(self, -1, 'Max'))
1260        sizer.Add(wx.StaticText(self, -1, 'Q range'))
1261             
1262        sizer.Add(self.qmin)
1263        sizer.Add(self.qmax)
1264        sizer_horizontal.Add(sizer)
1265        if object!=None:
1266            sizer_horizontal.Add(object)
1267       
1268        if object1!=None:
1269           boxsizer1.Add(object1) 
1270           boxsizer1.Add((10,10))
1271        boxsizer1.Add(sizer_horizontal)
1272        ## save state
1273        self.save_current_state()
1274        #----------------------------------------------------------------
1275        self.sizer5.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
1276        self.sizer5.Layout()
1277        self.Layout()
1278        self.SetScrollbars(20,20,200,100)
1279   
1280   
1281    def _fill_save_sizer(self):
1282        """
1283            Draw the layout for saving option
1284        """
1285        self.sizer6.Clear(True)
1286        box_description= wx.StaticBox(self, -1,"Save Status")
1287        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
1288        sizer_save = wx.BoxSizer(wx.HORIZONTAL)
1289       
1290        self.btSave_title = wx.StaticText(self, -1, 'Save the current panel status')
1291        self.btSave = wx.Button(self,wx.NewId(),'Save')
1292        self.btSave.Bind(wx.EVT_BUTTON, self.onSave,id= self.btSave.GetId())
1293        self.btSave.SetToolTipString("Save current panel status")
1294       
1295        sizer_save.Add(self.btSave_title) 
1296        sizer_save.Add((20,20),0, wx.LEFT|wx.RIGHT|wx.EXPAND,45) 
1297             
1298        sizer_save.Add(self.btSave)     
1299       
1300        boxsizer1.Add(sizer_save)
1301        self.sizer6.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
1302        self.sizer6.Layout()
1303        self.SetScrollbars(20,20,200,100)
1304       
1305       
1306               
Note: See TracBrowser for help on using the repository browser.