source: sasview/sansview/perspectives/fitting/modelpage.py @ 32d802f

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

plotted the result from fit 2D on a separate plot

  • Property mode set to 100644
File size: 17.7 KB
Line 
1import sys
2import wx
3import wx.lib
4import numpy
5import copy
6
7from sans.guicomm.events import StatusEvent   
8(ModelEventbox, EVT_MODEL_BOX) = wx.lib.newevent.NewEvent()
9_BOX_WIDTH = 80
10
11def format_number(value, high=False):
12    """
13        Return a float in a standardized, human-readable formatted string
14    """
15    try: 
16        value = float(value)
17    except:
18        print "returning 0"
19        return "0"
20   
21    if high:
22        return "%-6.4g" % value
23    else:
24        return "%-5.3g" % value
25
26   
27class ModelPage(wx.ScrolledWindow):
28    """
29        FitPanel class contains fields allowing to display results when
30        fitting  a model and one data
31        @note: For Fit to be performed the user should check at least one parameter
32        on fit Panel window.
33 
34    """
35    ## Internal name for the AUI manager
36    window_name = "Fit page"
37    ## Title to appear on top of the window
38    window_caption = "Fit Page"
39   
40   
41    def __init__(self, parent,model,name, *args, **kwargs):
42        wx.ScrolledWindow.__init__(self, parent, *args, **kwargs)
43        """
44            Initialization of the Panel
45        """
46        #self.scroll = wx.ScrolledWindow(self)
47        self.manager = None
48        self.parent  = parent
49        self.event_owner=None
50        #panel interface
51        self.vbox  = wx.BoxSizer(wx.VERTICAL)
52        self.sizer3 = wx.GridBagSizer(5,5)
53        self.sizer1 = wx.GridBagSizer(5,5)
54        self.sizer2 = wx.GridBagSizer(5,5)
55        self.sizer4 = wx.GridBagSizer(5,5)
56        self.sizer5 = wx.GridBagSizer(5,5)
57        self.static_line_1 = wx.StaticLine(self, -1)
58        self.modelbox = wx.ComboBox(self, -1)
59        id = wx.NewId()
60        self.vbox.Add(self.sizer3)
61        self.vbox.Add(self.sizer1)
62        self.vbox.Add(self.sizer2)
63        self.vbox.Add(self.static_line_1, 0, wx.EXPAND, 0)
64        self.vbox.Add(self.sizer5)
65        self.vbox.Add(self.sizer4)
66       
67        id = wx.NewId()
68        self.btClose =wx.Button(self,id,'Close')
69        self.btClose.Bind(wx.EVT_BUTTON, self.onClose,id=id)
70        self.btClose.SetToolTipString("Close page.")
71        ix = 1
72        iy = 1 
73        self.sizer4.Add(wx.StaticText(self, -1, 'Min'),(iy, ix),(1,1),\
74                            wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
75        ix += 1
76        self.sizer4.Add(wx.StaticText(self, -1, 'Max'),(iy, ix),(1,1),\
77                            wx.EXPAND|wx.ADJUST_MINSIZE, 0)
78        ix += 1
79        self.sizer4.Add(wx.StaticText(self, -1, 'Npts'),(iy, ix),(1,1),\
80                            wx.EXPAND|wx.ADJUST_MINSIZE, 0)
81        ix = 0
82        iy += 1
83        self.sizer4.Add(wx.StaticText(self, -1, 'x range'),(iy, ix),(1,1),\
84                            wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
85        ## Q range
86        self.qmin= 0.001
87        self.qmax= 0.1
88        self.num_points= 100
89       
90        ix += 1
91        self.xmin    = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
92        self.xmin.SetValue(format_number(self.qmin))
93        self.xmin.SetToolTipString("Minimun value of x in linear scale.")
94        self.xmin.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
95        self.xmin.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
96        self.sizer4.Add(self.xmin,(iy, ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
97       
98       
99        ix += 1
100        self.xmax    = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
101        self.xmax.SetValue(format_number(self.qmax))
102        self.xmax.SetToolTipString("Maximum value of x in linear scale.")
103        self.xmax.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
104        self.xmax.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
105        self.sizer4.Add(self.xmax,(iy,ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 0)
106        ix += 1
107        self.npts    = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
108        self.npts.SetValue(format_number(self.num_points))
109        self.npts.SetToolTipString("Number of point to plot.")
110        self.npts.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
111        self.npts.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
112       
113        self.sizer4.Add(self.npts,(iy,ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 0)
114        ix = 0
115        iy += 1
116        self.sizer4.Add((20,20),(iy, ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
117        ix +=3
118        self.sizer4.Add( self.btClose,(iy,ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 0)
119        ix = 0
120        iy = 1
121        self.sizer3.Add(wx.StaticText(self,-1,'Model'),(iy,ix),(1,1)\
122                  , wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
123        ix += 1
124        self.sizer3.Add(self.modelbox,(iy,ix),(1,1),  wx.EXPAND|wx.ADJUST_MINSIZE, 0)
125        #ix = 0
126        #iy += 1
127        ix += 1
128        id = wx.NewId()
129        self.model_view =wx.Button(self,id,'View 2D')
130        self.model_view.Bind(wx.EVT_BUTTON, self.onModel2D,id=id)
131        self.model_view.SetToolTipString("View model in 2D")
132        self.sizer3.Add(self.model_view,(iy,ix),(1,1),\
133                   wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
134        # contains link between  model ,all its parameters, and panel organization
135        self.parameters=[]
136        #contains link between a model and selected parameters to fit
137        self.param_toFit=[]
138        # model on which the fit would be performed
139        self.model=model
140        try:
141            #print"init modelpage",model.name
142            self.set_panel(model)
143        except:
144            raise
145        # preview selected model name
146        self.prevmodel_name=name
147        #print "model view prev_model",name
148        self.modelbox.SetValue(self.prevmodel_name)
149        # flag to check if the user has selected a new model in the combox box
150        self.model_hasChanged=False
151        #dictionary of model name and model class
152        self.model_list_box={}
153       
154        #enable model 2D draw
155        self.enable2D= False
156        # Data1D to make a deep comparison between 2 Data1D for checking data
157        #change
158        self.vbox.Layout()
159        self.vbox.Fit(self) 
160       
161        self.SetSizer(self.vbox)
162        self.SetScrollbars(20,20,55,40)
163        self.Centre()
164       
165    def set_range(self, qmin, qmax, npts):
166        """
167            Set the range for the plotted models
168            @param qmin: minimum Q
169            @param qmax: maximum Q
170            @param npts: number of Q bins
171        """
172        # Set the data members
173        self.qmin = qmin
174        self.qmax = qmax
175        self.num_points = npts
176       
177        # Set the controls
178        self.xmin.SetValue(format_number(self.qmin))
179        self.xmax.SetValue(format_number(self.qmax))
180        self.npts.SetValue(format_number(self.num_points))
181       
182    def onClose(self,event):
183        """ close the page associated with this panel"""
184        self.GrandParent.onClose()
185       
186    def set_owner(self,owner):
187        """
188            set owner of fitpage
189            @param owner: the class responsible of plotting
190        """
191        self.event_owner=owner   
192   
193 
194    def set_manager(self, manager):
195        """
196             set panel manager
197             @param manager: instance of plugin fitting
198        """
199        self.manager = manager
200       
201    def onModel2D(self, event):
202        """
203         call manager to plot model in 2D
204        """
205        # If the 2D display is not currently enabled, plot the model in 2D
206        # and set the enable2D flag.
207        if self.enable2D==False:
208            self.enable2D=True
209            self._draw_model()
210           
211        else:
212            print "enable is true:",self.enable2D
213            #self.manager.parent. show_panel(147)
214            self.manager.show_panel2D( id=None )
215            #self.manager.menu1.Append(event_id, new_panel.window_caption,
216            #             "Show %s plot panel" % new_panel.window_caption)
217           
218    def populate_box(self, dict):
219        """
220            Populate each combox box of each page
221            @param page: the page to populate
222        """
223        id=0
224        self.model_list_box=dict
225        list_name=[]
226        for item in  self.model_list_box.itervalues():
227            name = item.__name__
228            if hasattr(item, "name"):
229                name = item.name
230            list_name.append(name)
231        list_name.sort()   
232        for name in list_name:
233            self.modelbox.Insert(name,int(id))
234            id+=1
235        wx.EVT_COMBOBOX(self.modelbox,-1, self._on_select_model) 
236        return 0
237   
238    def select_model(self, model, name):
239        """
240            Select a new model
241            @param model: model object
242        """
243        self.model= model
244        print "select_model", model.__class__
245        self.set_panel(model)
246        self._draw_model(name)
247       
248        # Select the model from the combo box
249        items = self.modelbox.GetItems()
250        for i in range(len(items)):
251            print "model name",items[i],model.name, model.__class__.__name__
252            #if items[i]==model.__class__.__name__:
253            if items[i]==name:
254                self.modelbox.SetSelection(i)
255       
256    def _on_select_model(self,event):
257        """
258            react when a model is selected from page's combo box
259            post an event to its owner to draw an appropriate theory
260        """
261        self.model_view.SetFocus()
262        for item in self.model_list_box.itervalues():
263            name = item.__name__
264            if hasattr(item, "name"):
265                name = item.name
266            if name ==event.GetString():
267                model=item()
268                self.model= model
269                self.set_panel(model)
270                self.name= name
271                #self.manager.draw_model(model, name)
272                self.enable2D=False
273                self._draw_model(name)
274           
275           
276    def get_model_box(self): 
277        """ return reference to combox box self.model"""
278        return self.modelbox
279
280   
281    def get_param_list(self):
282        """
283            @return self.param_toFit: list containing  references to TextCtrl
284            checked.Theses TextCtrl will allow reference to parameters to fit.
285            @raise: if return an empty list of parameter fit will nnote work
286            properly so raise ValueError,"missing parameter to fit"
287        """
288        if self.param_toFit !=[]:
289            return self.param_toFit
290        else:
291            raise ValueError,"missing parameter to fit"
292       
293       
294    def set_panel(self,model):
295        """
296            Build the panel from the model content
297            @param model: the model selected in combo box for fitting purpose
298        """
299       
300        self.sizer2.Clear(True)
301        self.sizer1.Clear(True)
302        self.sizer5.Clear(True)
303        self.parameters = []
304        self.param_toFit=[]
305        self.model = model
306        keys = self.model.getParamList()
307        keys.sort()
308        description=None
309        if hasattr(self.model,'description'):
310            description =model.description
311        disp_list=self.model.getDispParamList()
312        ip=0
313        iq=1
314        if len(disp_list)>0:
315            disp = wx.StaticText(self, -1, 'Dispersion')
316            self.sizer5.Add(disp,( iq, ip),(1,1),  wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
317            ip += 1 
318            values = wx.StaticText(self, -1, 'Values')
319            self.sizer5.Add(values,( iq, ip),(1,1), wx.EXPAND|wx.ADJUST_MINSIZE, 0)
320           
321        disp_list.sort()
322        #print "went here",self.model.name,model.description
323        iy = 1
324        ix = 0
325        self.cb0 = wx.StaticText(self, -1,'Model Description:')
326        self.sizer1.Add(self.cb0,(iy, ix),(1,1),\
327                          wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
328        iy += 1
329       
330        self.cb01 = wx.StaticText(self, -1,str(description),style=wx.ALIGN_LEFT)
331        self.cb01.Wrap(400) 
332        #self.cb01 = wx.StaticText(self, -1,str(description),(45, 25),style=wx.ALIGN_LEFT)
333       
334        self.sizer1.Add(self.cb01,(iy, ix),(1,1),\
335                          wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
336        ix = 0
337        iy = 1
338        self.cb1 = wx.StaticText(self, -1,'Parameters')
339        self.sizer2.Add(self.cb1,(iy, ix),(1,1),\
340                          wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
341        ix +=1
342        self.text2_2 = wx.StaticText(self, -1, 'Values')
343        self.sizer2.Add(self.text2_2,(iy, ix),(1,1),\
344                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
345        ix +=1
346        self.text2_4 = wx.StaticText(self, -1, 'Units')
347        self.sizer2.Add(self.text2_4,(iy, ix),(1,1),\
348                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
349        self.text2_4.Hide()
350        for item in keys:
351            if not item in disp_list:
352                iy += 1
353                ix = 0
354                cb=wx.StaticText(self, -1, item)
355                self.sizer2.Add( cb,( iy, ix),(1,1),  wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
356                ix += 1
357                value= self.model.getParam(item)
358                ctl1 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH,20), style=wx.TE_PROCESS_ENTER)
359                ctl1.SetValue(str (format_number(value)))
360                ctl1.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
361                ctl1.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
362                self.sizer2.Add(ctl1, (iy,ix),(1,1), wx.EXPAND|wx.ADJUST_MINSIZE, 0)
363                ix +=1
364               
365                # Units
366                try:
367                    units = wx.StaticText(self, -1, self.model.details[item][0], style=wx.ALIGN_LEFT)
368                except:
369                    units = wx.StaticText(self, -1, "", style=wx.ALIGN_LEFT)
370                self.sizer2.Add(units, (iy,ix),(1,1),  wx.EXPAND|wx.ADJUST_MINSIZE, 0)
371            else:
372                ip = 0
373                iq += 1
374                cb=wx.StaticText(self, -1, item)
375                self.sizer5.Add( cb,( iq, ip),(1,1),  wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
376                wx.EVT_CHECKBOX(self, cb.GetId(), self._on_select_model)
377               
378                ip += 1
379                value= self.model.getParam(item)
380                ctl1 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH,20), style=wx.TE_PROCESS_ENTER)
381                ctl1.SetValue(str (format_number(value)))
382                ctl1.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
383                ctl1.Bind(wx.EVT_TEXT_ENTER,self._onparamEnter)
384                self.sizer5.Add(ctl1, (iq,ip),(1,1), wx.EXPAND|wx.ADJUST_MINSIZE, 0)
385            #save data
386            self.parameters.append([cb,ctl1])
387        iy+=1
388        self.sizer2.Add((20,20),(iy,ix),(1,1), wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
389       
390        #Display units text on panel
391        for item in keys:   
392            if self.model.details[item][0]!='':
393                self.text2_4.Show()
394                break
395            else:
396                self.text2_4.Hide()
397        self.vbox.Layout()
398        self.GrandParent.GetSizer().Layout()
399       
400       
401    def _onparamEnter(self,event):
402        """
403            when enter value on panel redraw model according to changed
404        """
405        self.set_model_parameter()
406       
407    def set_model_parameter(self):
408        if len(self.parameters) !=0 and self.model !=None:
409            # Flag to register when a parameter has changed.
410            is_modified = False
411            for item in self.parameters:
412                try:
413                     name=str(item[0].GetLabelText())
414                     value= float(item[1].GetValue())
415                     # If the value of the parameter has changed,
416                     # update the model and set the is_modified flag
417                     if value != self.model.getParam(name):
418                         self.model.setParam(name,value)
419                         is_modified = True
420                except:
421                     wx.PostEvent(self.parent.GrandParent, StatusEvent(status=\
422                            "Model Drawing  Error:wrong value entered : %s"% sys.exc_value))
423           
424            # Here we should check whether the boundaries have been modified.
425            # If qmin and qmax have been modified, update qmin and qmax and
426            # set the is_modified flag to True
427            if float(self.xmin.GetValue()) != self.qmin:
428                self.qmin = float(self.xmin.GetValue())
429                is_modified = True
430            if float(self.xmax.GetValue()) != self.qmax:
431                self.qmax = float(self.xmax.GetValue())
432                is_modified = True
433           
434            if float(self.npts.GetValue()) !=  self.num_points:
435                self.num_points = float(self.npts.GetValue())
436                is_modified = True
437         
438            if is_modified:
439                self._draw_model()           
440           
441    def _draw_model(self, name=None):
442        """
443            Method to draw or refresh a plotted model.
444            The method will use the data member from the model page
445            to build a call to the fitting perspective manager.
446           
447            [Note to coder: This way future changes will be done in only one place.]
448        """
449        if name==None:
450            name= self.model.name
451        self.manager.draw_model(self.model, name, 
452                                qmin=self.qmin, qmax=self.qmax,
453                                qstep= self.num_points,
454                                enable2D=self.enable2D)
455        """
456            self.manager.draw_model(self.model, self.model.name,
457                                    qmin=self.qmin, qmax=self.qmax,
458                                    qstep= self.num_points,
459                                    enable2D=self.enable2D)
460        """
461   
462       
463             
Note: See TracBrowser for help on using the repository browser.