source: sasview/fittingview/src/sans/perspectives/fitting/fitpage.py @ d40038e

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 d40038e was d40038e, checked in by Jae Cho <jhjcho@…>, 12 years ago

fix: ticket #10

  • Property mode set to 100644
File size: 135.2 KB
Line 
1
2
3import sys
4import wx
5import wx.lib.newevent
6import numpy
7import copy
8import math
9import time
10from sans.models.dispersion_models import ArrayDispersion, GaussianDispersion
11#from sans.dataloader.data_info import Data1D
12from sans.guiframe.events import StatusEvent
13from sans.guiframe.events import NewPlotEvent 
14from sans.guiframe.dataFitting import check_data_validity
15from sans.guiframe.utils import format_number
16from sans.guiframe.utils import check_float
17
18(Chi2UpdateEvent, EVT_CHI2_UPDATE)   = wx.lib.newevent.NewEvent()
19_BOX_WIDTH = 76
20_DATA_BOX_WIDTH = 300
21SMEAR_SIZE_L = 0.00
22SMEAR_SIZE_H = 0.00
23
24import basepage
25from basepage import BasicPage
26from basepage import PageInfoEvent
27from sans.models.qsmearing import smear_selection
28
29class FitPage(BasicPage):
30    """
31    FitPanel class contains fields allowing to display results when
32    fitting  a model and one data
33   
34    :note: For Fit to be performed the user should check at least one parameter
35        on fit Panel window.
36    """
37   
38    def __init__(self,parent, color='rand'):
39        """
40        Initialization of the Panel
41        """
42        BasicPage.__init__(self, parent, color=color)
43       
44        ## draw sizer
45        self._fill_data_sizer()
46        self.is_2D = None
47        self.fit_started = False
48        self.weightbt_string = None
49        # get smear info from data
50        self._get_smear_info()
51        self._fill_model_sizer( self.sizer1)
52        self._get_defult_custom_smear()
53        self._fill_range_sizer() 
54        self._set_smear(self.data)
55        self.Bind(EVT_CHI2_UPDATE, self.on_complete_chisqr)
56        # bind key event
57        self.Bind(wx.EVT_RIGHT_DOWN, self.on_right_down)
58        self._set_bookmark_flag(False)
59        self._set_save_flag(False)
60        self._set_preview_flag(False)
61        self._set_copy_flag(False)
62        self._set_paste_flag(False)
63        self.btFit.SetFocus()
64        self.enable_fit_button()
65        self.fill_data_combobox(data_list=self.data_list)
66        #create a default data for an empty panel
67        self.create_default_data()
68   
69    def enable_fit_button(self):
70        """
71        Enable fit button if data is valid and model is valid
72        """
73        flag = check_data_validity(self.data) & (self.model is not None)
74        self.btFit.Enable(flag)
75       
76    def _fill_data_sizer(self):
77        """
78        fill sizer 0 with data info
79        """
80        self.data_box_description = wx.StaticBox(self, -1, 'I(q) Data Source')
81        if check_data_validity(self.data):
82            dname_color = wx.BLUE
83        else:
84            dname_color = wx.RED
85        self.data_box_description.SetForegroundColour(dname_color)
86        boxsizer1 = wx.StaticBoxSizer(self.data_box_description, wx.VERTICAL)
87        #----------------------------------------------------------
88        sizer_data = wx.BoxSizer(wx.HORIZONTAL)
89        self.dataSource = wx.ComboBox(self, -1, style=wx.CB_READONLY)
90        wx.EVT_COMBOBOX(self.dataSource, -1, self.on_select_data)
91        self.dataSource.SetMinSize((_DATA_BOX_WIDTH, -1))
92        sizer_data.Add(wx.StaticText(self, -1, 'Name : '))
93        sizer_data.Add(self.dataSource)
94        sizer_data.Add( (0,5) )
95        boxsizer1.Add(sizer_data,0, wx.ALL, 10)
96        self.sizer0.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
97        self.sizer0.Layout()
98       
99    def enable_datasource(self):
100        """
101        Enable or disable data source control depending on existing data
102        """
103        if not self.data_list:
104            self.dataSource.Disable()
105        else:
106            self.dataSource.Enable()
107           
108    def fill_data_combobox(self, data_list):
109        """
110        Get a list of data and fill the corresponding combobox
111        """
112        self.dataSource.Clear()
113        self.data_list = data_list
114        self.enable_datasource()
115        if len(data_list) > 0:
116            #find the maximum range covering all data
117            qmin, qmax, npts = self.compute_data_set_range(data_list)
118            self.qmin_data_set = qmin
119            self.qmax_data_set = qmax
120            self.npts_data_set = npts
121
122            self.qmin.SetValue(str(self.qmin_data_set))
123            self.qmax.SetValue(str(self.qmax_data_set))
124            self.qmin.SetBackgroundColour("white")
125            self.qmax.SetBackgroundColour("white")
126            self.qmin_x = self.qmin_data_set
127            self.qmax_x =  self.qmax_data_set
128            self.state.qmin = self.qmin_x
129            self.state.qmax = self.qmax_x
130        is_data = False   
131        for data in self.data_list:
132            if data is not None:
133                self.dataSource.Append(str(data.name), clientData=data)
134                if not is_data:
135                    is_data = check_data_validity(data)
136        self.dataSource.SetSelection(0)
137        self.on_select_data(event=None)
138               
139    def on_select_data(self, event=None):
140        """
141        """
142        if event is None and self.dataSource.GetCount() > 0:
143            data = self.dataSource.GetClientData(0)
144            self.set_data(data)
145        elif self.dataSource.GetCount() > 0:
146            pos = self.dataSource.GetSelection()
147            data = self.dataSource.GetClientData(pos)
148            self.set_data(data)
149   
150   
151       
152    def _on_fit_complete(self):
153        """
154        When fit is complete ,reset the fit button label.
155        """
156        self.btFit.SetLabel("Fit")
157       
158    def _is_2D(self):
159        """
160        Check if data_name is Data2D
161       
162        :return: True or False
163       
164        """
165        if self.data.__class__.__name__ == "Data2D" or \
166                        self.enable2D:
167            return True
168        return False
169           
170    def _on_engine_change(self, name):
171        """
172        get the current name of the fit engine type
173         and update the panel accordingly
174        """
175       
176        self.engine_type = str(name)
177        self.state.engine_type = self.engine_type
178        if not self.is_mac:
179            if len(self.parameters) == 0:
180                self.Layout()
181                return
182            self.Layout()
183            self.Refresh()
184       
185 
186       
187    def _fill_range_sizer(self):
188        """
189        Fill the sizer containing the plotting range
190        add  access to npts
191        """
192        is_2Ddata = False
193       
194        # Check if data is 2D
195        if self.data.__class__.__name__ ==  "Data2D" or \
196                        self.enable2D:
197            is_2Ddata = True
198           
199        title = "Fitting"
200        #smear messages & titles
201        smear_message_none  =  "No smearing is selected..."
202        smear_message_dqdata  =  "The dQ data is being used for smearing..."
203        smear_message_2d  =  \
204              "Higher accuracy is very time-expensive. Use it with care..."
205        smear_message_new_ssmear  = \
206              "Please enter only the value of interest to customize smearing..."
207        smear_message_new_psmear  = \
208              "Please enter both; the dQ will be generated by interpolation..."
209        smear_message_2d_x_title = "<dQp>[1/A]:"
210        smear_message_2d_y_title = "<dQs>[1/A]:"       
211        smear_message_pinhole_min_title = "dQ_low[1/A]:"
212        smear_message_pinhole_max_title = "dQ_high[1/A]:"
213        smear_message_slit_height_title = "Slit height[1/A]:"
214        smear_message_slit_width_title = "Slit width[1/A]:"
215       
216        self._get_smear_info()
217       
218        #Sizers
219        box_description_range = wx.StaticBox(self, -1,str(title))
220        box_description_range.SetForegroundColour(wx.BLUE)
221        boxsizer_range = wx.StaticBoxSizer(box_description_range, wx.VERTICAL)     
222        self.sizer_set_smearer = wx.BoxSizer(wx.VERTICAL)
223        sizer_smearer = wx.BoxSizer(wx.HORIZONTAL)
224        self.sizer_new_smear= wx.BoxSizer(wx.HORIZONTAL)
225        self.sizer_set_masking = wx.BoxSizer(wx.HORIZONTAL)
226        sizer_chi2 = wx.BoxSizer(wx.VERTICAL)
227        smear_set_box= wx.StaticBox(self, -1,'Set Instrumental Smearing')
228        sizer_smearer_box = wx.StaticBoxSizer(smear_set_box, wx.HORIZONTAL)
229        sizer_smearer_box.SetMinSize((_DATA_BOX_WIDTH, 60))
230       
231        weighting_set_box= wx.StaticBox(self, -1,\
232                                'Set Weighting by Selecting dI Source')
233        weighting_box = wx.StaticBoxSizer(weighting_set_box, wx.HORIZONTAL)
234        sizer_weighting = wx.BoxSizer(wx.HORIZONTAL)
235        weighting_box.SetMinSize((_DATA_BOX_WIDTH, 40))
236        #Filling the sizer containing weighting info.
237        self.dI_noweight = wx.RadioButton(self, -1, 'No Weighting', 
238                                          style=wx.RB_GROUP)
239        self.dI_didata = wx.RadioButton(self, -1, 'Use dI Data')
240        self.dI_sqrdata = wx.RadioButton(self, -1, 'Use |sqrt(I Data)|')
241        self.dI_idata = wx.RadioButton(self, -1, 'Use |I Data|')
242        self.Bind(wx.EVT_RADIOBUTTON, self.onWeighting, 
243                  id=self.dI_noweight.GetId())
244        self.Bind(wx.EVT_RADIOBUTTON, self.onWeighting, 
245                  id=self.dI_didata.GetId())
246        self.Bind(wx.EVT_RADIOBUTTON, self.onWeighting, 
247                  id=self.dI_sqrdata.GetId())
248        self.Bind(wx.EVT_RADIOBUTTON, self.onWeighting, 
249                  id=self.dI_idata.GetId())
250        self.dI_didata.SetValue(True)
251        # add 4 types of weighting to the sizer
252        sizer_weighting.Add( self.dI_noweight,0, wx.LEFT, 10)
253        sizer_weighting.Add((14,10))
254        sizer_weighting.Add( self.dI_didata)
255        sizer_weighting.Add((14,10))
256        sizer_weighting.Add( self.dI_sqrdata) 
257        sizer_weighting.Add((14,10))
258        sizer_weighting.Add( self.dI_idata) 
259        sizer_weighting.Add((10,10)) 
260        self.dI_noweight.Enable(False)
261        self.dI_didata.Enable(False)
262        self.dI_sqrdata.Enable(False)
263        self.dI_idata.Enable(False)
264        weighting_box.Add(sizer_weighting)
265       
266        sizer_fit = wx.GridSizer(2, 4, 2, 6)
267       
268        # combobox for smear2d accuracy selection
269        self.smear_accuracy = wx.ComboBox(self, -1,size=(50,-1),
270                                          style=wx.CB_READONLY)
271        self._set_accuracy_list()
272        self.smear_accuracy.SetValue(self.smear2d_accuracy)
273        self.smear_accuracy.SetSelection(0)
274        self.smear_accuracy.SetToolTipString(\
275            "'Higher' uses more Gaussian points for smearing computation.")
276                   
277        wx.EVT_COMBOBOX(self.smear_accuracy,-1, self._on_select_accuracy)
278
279        #Fit button
280        self.btFit = wx.Button(self,wx.NewId(),'Fit', size=(88,25))
281        self.default_bt_colour =  self.btFit.GetDefaultAttributes()
282        self.btFit.Bind(wx.EVT_BUTTON, self._onFit,id= self.btFit.GetId())
283        self.btFit.SetToolTipString("Start fitting.")
284       
285        #textcntrl for custom resolution
286        self.smear_pinhole_max = self.ModelTextCtrl(self, -1,
287                            size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
288                            text_enter_callback = self.onPinholeSmear)
289        self.smear_pinhole_min = self.ModelTextCtrl(self, -1,
290                            size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
291                            text_enter_callback = self.onPinholeSmear)
292        self.smear_slit_height= self.ModelTextCtrl(self, -1,
293                            size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
294                            text_enter_callback = self.onSlitSmear)
295        self.smear_slit_width = self.ModelTextCtrl(self, -1,
296                            size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
297                            text_enter_callback = self.onSlitSmear)
298
299        ## smear
300        self.smear_data_left= BGTextCtrl(self, -1, 
301                                        size=(_BOX_WIDTH-25,20), style=0)
302        self.smear_data_left.SetValue(str(self.dq_l))
303        self.smear_data_right = BGTextCtrl(self, -1, 
304                                        size=(_BOX_WIDTH-25,20), style=0)
305        self.smear_data_right.SetValue(str(self.dq_r))
306
307        #set default values for smear
308        self.smear_pinhole_max.SetValue(str(self.dx_max))
309        self.smear_pinhole_min.SetValue(str(self.dx_min))
310        self.smear_slit_height.SetValue(str(self.dxl))
311        self.smear_slit_width.SetValue(str(self.dxw))
312
313        #Filling the sizer containing instruments smearing info.
314        self.disable_smearer = wx.RadioButton(self, -1, 
315                                              'None', style=wx.RB_GROUP)
316        self.enable_smearer = wx.RadioButton(self, -1, 
317                                             'Use dQ Data')
318        #self.enable_smearer.SetToolTipString(
319        #"Click to use the loaded dQ data for smearing.")
320        self.pinhole_smearer = wx.RadioButton(self, -1, 
321                                              'Custom Pinhole Smear')
322        #self.pinhole_smearer.SetToolTipString
323        #("Click to input custom resolution for pinhole smearing.")
324        self.slit_smearer = wx.RadioButton(self, -1, 'Custom Slit Smear')
325        #self.slit_smearer.SetToolTipString
326        #("Click to input custom resolution for slit smearing.")
327        self.Bind(wx.EVT_RADIOBUTTON, self.onSmear, 
328                  id=self.disable_smearer.GetId())
329        self.Bind(wx.EVT_RADIOBUTTON, self.onSmear, 
330                  id=self.enable_smearer.GetId())
331        self.Bind(wx.EVT_RADIOBUTTON, self.onPinholeSmear, 
332                  id=self.pinhole_smearer.GetId())
333        self.Bind(wx.EVT_RADIOBUTTON, self.onSlitSmear, 
334                  id=self.slit_smearer.GetId())
335        self.disable_smearer.SetValue(True)
336       
337        # add 4 types of smearing to the sizer
338        sizer_smearer.Add( self.disable_smearer,0, wx.LEFT, 10)
339        sizer_smearer.Add((10,10))
340        sizer_smearer.Add( self.enable_smearer)
341        sizer_smearer.Add((10,10))
342        sizer_smearer.Add( self.pinhole_smearer ) 
343        sizer_smearer.Add((10,10))
344        sizer_smearer.Add( self.slit_smearer ) 
345        sizer_smearer.Add((10,10))       
346       
347        # StaticText for chi2, N(for fitting), Npts
348        self.tcChi    =  BGTextCtrl(self, -1, "-", size=(75,20), style=0)
349        self.tcChi.SetToolTipString("Chi2/Npts(Fit)")
350        self.Npts_fit    =  BGTextCtrl(self, -1, "-", size=(75,20), style=0)
351        self.Npts_fit.SetToolTipString(\
352                            " Npts : number of points selected for fitting")
353        self.Npts_total  =  self.ModelTextCtrl(self, -1, 
354                        size=(_BOX_WIDTH, 20), 
355                        style=wx.TE_PROCESS_ENTER, 
356                        text_enter_callback=self._onQrangeEnter)
357        self.Npts_total.SetValue(format_number(self.npts_x))
358        self.Npts_total.SetToolTipString(\
359                                " Total Npts : total number of data points")
360       
361        # Update and Draw button
362        self.draw_button = wx.Button(self,wx.NewId(),'Compute', size=(88,24))
363        self.draw_button.Bind(wx.EVT_BUTTON, \
364                              self._onDraw,id= self.draw_button.GetId())
365        self.draw_button.SetToolTipString("Compute and Draw.")
366       
367        box_description_1= wx.StaticText(self, -1,'   Chi2/Npts')
368        box_description_2= wx.StaticText(self, -1,'Npts(Fit)')
369        box_description_3= wx.StaticText(self, -1,'Total Npts')
370        box_description_3.SetToolTipString( \
371                                " Total Npts : total number of data points")
372        #box_description_4= wx.StaticText(self, -1,' ')
373       
374       
375        sizer_fit.Add(box_description_1,0,0)
376        sizer_fit.Add(box_description_2,0,0)
377        sizer_fit.Add(box_description_3,0,0)       
378        sizer_fit.Add(self.draw_button,0,0)
379        sizer_fit.Add(self.tcChi,0,0)
380        sizer_fit.Add(self.Npts_fit ,0,0)
381        sizer_fit.Add(self.Npts_total,0,0)
382        sizer_fit.Add(self.btFit,0,0) 
383
384        # StaticText for smear
385        #self.tcChi    =  wx.StaticText(self, -1, "-", style=wx.ALIGN_LEFT)
386        self.smear_description_none =  wx.StaticText(self, -1, 
387                                    smear_message_none , style=wx.ALIGN_LEFT)
388        self.smear_description_dqdata =  wx.StaticText(self, 
389                                -1, smear_message_dqdata , style=wx.ALIGN_LEFT)
390        self.smear_description_type =  wx.StaticText(self,
391                                     -1, "Type:" , style=wx.ALIGN_LEFT)
392        self.smear_description_accuracy_type =  wx.StaticText(self, -1, 
393                                        "Accuracy:" , style=wx.ALIGN_LEFT)
394        self.smear_description_smear_type =  BGTextCtrl(self, -1, 
395                                                        size=(57,20), style=0)
396        self.smear_description_smear_type.SetValue(str(self.dq_l))
397        self.SetBackgroundColour(self.GetParent().GetBackgroundColour())
398        self.smear_description_2d =  wx.StaticText(self, -1, 
399                                    smear_message_2d  , style=wx.ALIGN_LEFT)
400        self.smear_message_new_s = wx.StaticText(self, -1,
401                         smear_message_new_ssmear, style=wx.ALIGN_LEFT)
402        self.smear_message_new_p = wx.StaticText(self, -1,
403                            smear_message_new_psmear , style=wx.ALIGN_LEFT)
404        self.smear_description_2d_x     =  wx.StaticText(self, -1, 
405                            smear_message_2d_x_title  , style=wx.ALIGN_LEFT)
406        self.smear_description_2d_x.SetToolTipString(\
407                                        "  dQp(parallel) in q_r direction.")
408        self.smear_description_2d_y     =  wx.StaticText(self, -1, 
409                            smear_message_2d_y_title  , style=wx.ALIGN_LEFT)
410        self.smear_description_2d_y.SetToolTipString(\
411                                    " dQs(perpendicular) in q_phi direction.")
412        self.smear_description_pin_min     =  wx.StaticText(self, -1, 
413                        smear_message_pinhole_min_title  , style=wx.ALIGN_LEFT)
414        self.smear_description_pin_max     =  wx.StaticText(self, -1, 
415                        smear_message_pinhole_max_title  , style=wx.ALIGN_LEFT)
416        self.smear_description_slit_height    =  wx.StaticText(self, -1, 
417                        smear_message_slit_height_title   , style=wx.ALIGN_LEFT)
418        self.smear_description_slit_width    =  wx.StaticText(self, -1, 
419                        smear_message_slit_width_title   , style=wx.ALIGN_LEFT)
420       
421        #arrange sizers
422        #boxsizer1.Add( self.tcChi ) 
423        self.sizer_set_smearer.Add(sizer_smearer )
424        self.sizer_set_smearer.Add((10,10))
425        self.sizer_set_smearer.Add( self.smear_description_none,
426                                    0, wx.CENTER, 10 ) 
427        self.sizer_set_smearer.Add( self.smear_description_dqdata,
428                                    0, wx.CENTER, 10 )
429        self.sizer_set_smearer.Add( self.smear_description_2d,
430                                    0, wx.CENTER, 10 )
431        self.sizer_new_smear.Add( self.smear_description_type,
432                                  0, wx.CENTER, 10 )
433        self.sizer_new_smear.Add( self.smear_description_accuracy_type,
434                                  0, wx.CENTER, 10 )
435        self.sizer_new_smear.Add( self.smear_accuracy )
436        self.sizer_new_smear.Add( self.smear_description_smear_type,
437                                  0, wx.CENTER, 10 )
438        self.sizer_new_smear.Add((15,-1))
439        self.sizer_new_smear.Add( self.smear_description_2d_x,
440                                  0, wx.CENTER, 10 )
441        self.sizer_new_smear.Add( self.smear_description_pin_min,
442                                  0, wx.CENTER, 10 )
443        self.sizer_new_smear.Add( self.smear_description_slit_height,
444                                  0, wx.CENTER, 10 )
445
446        self.sizer_new_smear.Add( self.smear_pinhole_min,
447                                  0, wx.CENTER, 10 )
448        self.sizer_new_smear.Add( self.smear_slit_height,
449                                  0, wx.CENTER, 10 )
450        self.sizer_new_smear.Add( self.smear_data_left,
451                                  0, wx.CENTER, 10 )
452        self.sizer_new_smear.Add((20,-1))
453        self.sizer_new_smear.Add( self.smear_description_2d_y,
454                                  0, wx.CENTER, 10 )
455        self.sizer_new_smear.Add( self.smear_description_pin_max,
456                                  0, wx.CENTER, 10 )
457        self.sizer_new_smear.Add( self.smear_description_slit_width,
458                                  0, wx.CENTER, 10 )
459
460        self.sizer_new_smear.Add( self.smear_pinhole_max,0, wx.CENTER, 10 )
461        self.sizer_new_smear.Add( self.smear_slit_width,0, wx.CENTER, 10 )
462        self.sizer_new_smear.Add( self.smear_data_right,0, wx.CENTER, 10 )
463           
464        self.sizer_set_smearer.Add( self.smear_message_new_s,0, wx.CENTER, 10)
465        self.sizer_set_smearer.Add( self.smear_message_new_p,0, wx.CENTER, 10)
466        self.sizer_set_smearer.Add((5,2))
467        self.sizer_set_smearer.Add( self.sizer_new_smear,0, wx.CENTER, 10 )
468       
469        # add all to chi2 sizer
470        sizer_smearer_box.Add(self.sizer_set_smearer)   
471        sizer_chi2.Add(sizer_smearer_box)
472        sizer_chi2.Add((-1,5))
473        sizer_chi2.Add(weighting_box) 
474        sizer_chi2.Add((-1,5))
475       
476        # hide all smear messages and textctrl
477        self._hide_all_smear_info()
478       
479        # get smear_selection
480        self.current_smearer= smear_selection( self.data, self.model )
481
482        # Show only the relevant smear messages, etc
483        if self.current_smearer == None:
484            if not is_2Ddata:
485                self.smear_description_none.Show(True)
486                self.enable_smearer.Disable() 
487            else:
488                self.smear_description_none.Show(True)
489                #self.smear_description_2d.Show(True)
490                #self.pinhole_smearer.Disable()
491                self.slit_smearer.Disable()   
492                #self.enable_smearer.Disable()
493            if self.data == None:
494                self.slit_smearer.Disable() 
495                self.pinhole_smearer.Disable() 
496                self.enable_smearer.Disable() 
497        else: self._show_smear_sizer()
498        boxsizer_range.Add(self.sizer_set_masking)
499         #2D data? default
500        is_2Ddata = False
501       
502        #check if it is 2D data
503        if self.data.__class__.__name__ ==  "Data2D" or \
504                        self.enable2D:
505            is_2Ddata = True
506           
507        self.sizer5.Clear(True)
508     
509        self.qmin  = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH,20),
510                                          style=wx.TE_PROCESS_ENTER,
511                                    text_enter_callback = self._onQrangeEnter)
512        self.qmin.SetValue(str(self.qmin_x))
513        self.qmin.SetToolTipString("Minimun value of Q in linear scale.")
514     
515        self.qmax  = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH,20),
516                                          style=wx.TE_PROCESS_ENTER,
517                                        text_enter_callback=self._onQrangeEnter)
518        self.qmax.SetValue(str(self.qmax_x))
519        self.qmax.SetToolTipString("Maximum value of Q in linear scale.")
520        """
521        self.theory_npts_tcrtl  = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH, 20),
522                        style=wx.TE_PROCESS_ENTER,
523                        text_enter_callback=self._onQrangeEnter)
524        self.theory_npts_tcrtl.SetValue(format_number(self.npts_x))
525        self.theory_npts_tcrtl.SetToolTipString("Number of point to plot.")
526        """
527        id = wx.NewId()
528        self.reset_qrange =wx.Button(self,id,'Reset',size=(77,20))
529     
530        self.reset_qrange.Bind(wx.EVT_BUTTON, self.on_reset_clicked,id=id)
531        self.reset_qrange.SetToolTipString("Reset Q range to the default values")
532     
533        sizer_horizontal=wx.BoxSizer(wx.HORIZONTAL)
534        sizer= wx.GridSizer(2, 4, 2, 6)
535
536        self.btEditMask = wx.Button(self,wx.NewId(),'Editor', size=(88,23))
537        self.btEditMask.Bind(wx.EVT_BUTTON, self._onMask,id= self.btEditMask.GetId())
538        self.btEditMask.SetToolTipString("Edit Mask.")
539        self.EditMask_title = wx.StaticText(self, -1, ' Masking(2D)')
540
541        sizer.Add(wx.StaticText(self, -1, '   Q range'))     
542        sizer.Add(wx.StaticText(self, -1, ' Min[1/A]'))
543        sizer.Add(wx.StaticText(self, -1, ' Max[1/A]'))
544        sizer.Add(self.EditMask_title)
545        #sizer.Add(wx.StaticText(self, -1, ''))
546        sizer.Add(self.reset_qrange)   
547        sizer.Add(self.qmin)
548        sizer.Add(self.qmax)
549        #sizer.Add(self.theory_npts_tcrtl)
550        sizer.Add(self.btEditMask)
551        boxsizer_range.Add(sizer_chi2) 
552        boxsizer_range.Add((10,10))
553        boxsizer_range.Add(sizer)
554       
555        boxsizer_range.Add((10,15))
556        boxsizer_range.Add(sizer_fit)
557        if is_2Ddata:
558            self.btEditMask.Enable() 
559            self.EditMask_title.Enable() 
560        else:
561            self.btEditMask.Disable() 
562            self.EditMask_title.Disable()
563        ## save state
564        self.save_current_state()
565        self.sizer5.Add(boxsizer_range,0, wx.EXPAND | wx.ALL, 10)
566        self.sizer5.Layout()
567
568       
569    def _fill_model_sizer(self, sizer):
570        """
571        fill sizer containing model info
572        """
573        ##Add model function Details button in fitpanel.
574        ##The following 3 lines are for Mac. Let JHC know before modifying...
575        title = "Model"
576        self.formfactorbox = None
577        self.multifactorbox = None
578        self.mbox_description= wx.StaticBox(self, -1,str(title))
579        boxsizer1 = wx.StaticBoxSizer(self.mbox_description, wx.VERTICAL)
580        self.mbox_description.SetForegroundColour(wx.RED)
581        id = wx.NewId()
582        self.model_help =wx.Button(self,id,'Details', size=(80,23))
583        self.model_help.Bind(wx.EVT_BUTTON, self.on_model_help_clicked,id=id)
584        self.model_help.SetToolTipString("Model Function Help")
585        id = wx.NewId()
586        self.model_view = wx.Button(self, id,"Show 2D", size=(80, 23))
587        self.model_view.Bind(wx.EVT_BUTTON, self._onModel2D, id=id)
588        hint = "toggle view of model from 1D to 2D  or 2D to 1D"
589        self.model_view.SetToolTipString(hint)
590     
591        self.shape_rbutton = wx.RadioButton(self, -1, 'Shapes',
592                                             style=wx.RB_GROUP)
593        self.shape_indep_rbutton = wx.RadioButton(self, -1, "Shape-Independent")
594        self.struct_rbutton = wx.RadioButton(self, -1, "Structure Factor ")
595        self.plugin_rbutton = wx.RadioButton(self, -1, "Customized Models")
596               
597        self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
598                            id= self.shape_rbutton.GetId()) 
599        self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
600                            id= self.shape_indep_rbutton.GetId()) 
601        self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
602                            id= self.struct_rbutton.GetId()) 
603        self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
604                            id= self.plugin_rbutton.GetId()) 
605        #MAC needs SetValue
606        self.shape_rbutton.SetValue(True)
607     
608        sizer_radiobutton = wx.GridSizer(2, 3, 5, 5)
609        sizer_radiobutton.Add(self.shape_rbutton)
610        sizer_radiobutton.Add(self.shape_indep_rbutton)
611        #sizer_radiobutton.Add((5, 5))
612        sizer_radiobutton.Add(self.model_view,1, wx.LEFT, 20)
613        sizer_radiobutton.Add(self.plugin_rbutton)
614        sizer_radiobutton.Add(self.struct_rbutton)
615        #sizer_radiobutton.Add((5, 5))
616        sizer_radiobutton.Add(self.model_help,1, wx.LEFT, 20)
617       
618        sizer_selection = wx.BoxSizer(wx.HORIZONTAL)
619        mutifactor_selection = wx.BoxSizer(wx.HORIZONTAL)
620       
621        self.text1 = wx.StaticText(self,-1,"" )
622        self.text2 = wx.StaticText(self,-1,"P(Q)*S(Q)" )
623        self.mutifactor_text = wx.StaticText( self,-1,"No. of Shells: ")
624        self.mutifactor_text1 = wx.StaticText( self,-1,"" )
625        self.show_sld_button = wx.Button( self,-1,"Show SLD Profile" )
626        self.show_sld_button.Bind(wx.EVT_BUTTON,self._on_show_sld)
627
628        self.formfactorbox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
629        if self.model!= None:
630            self.formfactorbox.SetValue(self.model.name)
631        self.structurebox = wx.ComboBox(self, -1, style=wx.CB_READONLY)
632        self.multifactorbox = wx.ComboBox(self, -1, style=wx.CB_READONLY)
633        self.initialize_combox()
634        wx.EVT_COMBOBOX(self.formfactorbox, -1, self._on_select_model)
635
636        wx.EVT_COMBOBOX(self.structurebox, -1, self._on_select_model)
637        wx.EVT_COMBOBOX(self.multifactorbox, -1, self._on_select_model)
638        ## check model type to show sizer
639        if self.model !=None:
640            self._set_model_sizer_selection(self.model)
641       
642        sizer_selection.Add(self.text1)
643        sizer_selection.Add((5, 5))
644        sizer_selection.Add(self.formfactorbox)
645        sizer_selection.Add((5, 5))
646        sizer_selection.Add(self.text2)
647        sizer_selection.Add((5, 5))
648        sizer_selection.Add(self.structurebox)
649       
650        mutifactor_selection.Add((10,5))
651        mutifactor_selection.Add(self.mutifactor_text)
652        mutifactor_selection.Add(self.multifactorbox)
653        mutifactor_selection.Add((5, 5))
654        mutifactor_selection.Add(self.mutifactor_text1)
655        mutifactor_selection.Add((10, 5))
656        mutifactor_selection.Add(self.show_sld_button)
657
658       
659        boxsizer1.Add(sizer_radiobutton)
660        boxsizer1.Add((10, 10))
661        boxsizer1.Add(sizer_selection)
662        boxsizer1.Add((10, 10))
663        boxsizer1.Add(mutifactor_selection)
664       
665        self._set_multfactor_combobox()
666        self.multifactorbox.SetSelection(1)
667        self.show_sld_button.Hide()
668        sizer.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
669        sizer.Layout()
670       
671    def _set_sizer_dispersion(self):
672        """
673        draw sizer with gaussian dispersity parameters
674        """
675        self.fittable_param=[]
676        self.fixed_param=[]
677        self.orientation_params_disp=[]
678
679        self.sizer4_4.Clear(True)
680        if self.model==None:
681            ##no model is selected
682            return
683        if not self.enable_disp.GetValue():
684            ## the user didn't select dispersity display
685            return 
686           
687        self._reset_dispersity()
688       
689        ## fill a sizer with the combobox to select dispersion type
690        #sizer_select_dispers = wx.BoxSizer(wx.HORIZONTAL) 
691        model_disp = wx.StaticText(self, -1, 'Function')
692        CHECK_STATE = self.cb1.GetValue()     
693        import sans.models.dispersion_models 
694        self.polydisp= sans.models.dispersion_models.models
695
696        ix = 0
697        iy = 0
698        disp = wx.StaticText(self, -1, ' ')
699        self.sizer4_4.Add(disp,( iy, ix),(1,1), 
700                           wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
701        ix += 1 
702        values = wx.StaticText(self, -1, 'PD[ratio]')
703        polytext = "Polydispersity (= STD/mean); "
704        polytext +=  "the standard deviation over the mean value."
705        values.SetToolTipString(polytext)
706
707        self.sizer4_4.Add(values,( iy, ix),(1,1), wx.EXPAND|wx.ADJUST_MINSIZE, 
708                          0)
709        ix +=2 
710        if self.is_mac:
711            err_text = 'Error'
712        else:
713            err_text = ''
714        self.text_disp_1 = wx.StaticText(self, -1, err_text)
715        self.sizer4_4.Add( self.text_disp_1,(iy, ix),(1,1),\
716                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
717       
718       
719        ix +=1 
720        self.text_disp_min = wx.StaticText(self, -1, 'Min')
721        self.sizer4_4.Add(self.text_disp_min,(iy, ix),(1,1),\
722                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
723
724        ix +=1 
725        self.text_disp_max = wx.StaticText(self, -1, 'Max')
726        self.sizer4_4.Add(self.text_disp_max,(iy, ix),(1,1),\
727                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
728       
729        ix += 1 
730        npts = wx.StaticText(self, -1, 'Npts')
731        npts.SetToolTipString("Number of sampling points for the numerical\n\
732        integration over the distribution function.")
733        self.sizer4_4.Add(npts,( iy, ix),(1,1), wx.EXPAND|wx.ADJUST_MINSIZE, 0)
734        ix += 1 
735        nsigmas = wx.StaticText(self, -1, 'Nsigs')
736        nsigmas.SetToolTipString("   Number of sigmas between which the range\n\
737         of the distribution function will be used for weighting. \n\
738        The value '3' covers 99.5% for Gaussian distribution \n\
739        function. Note: Not recommended to change this value.")
740        self.sizer4_4.Add(nsigmas,( iy, ix),(1,1), wx.EXPAND|wx.ADJUST_MINSIZE,
741                           0)
742        ix +=1 
743        self.sizer4_4.Add(model_disp,(iy, ix),(1,1),\
744                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
745       
746        if self.engine_type=="park":
747            self.text_disp_max.Show(True)
748            self.text_disp_min.Show(True)
749
750        for item in self.model.dispersion.keys():
751            if not item in self.model.orientation_params:
752                if not self.disp_cb_dict.has_key(item):
753                    self.disp_cb_dict[item]= None
754                name0="Distribution of " + item
755                name1=item+".width"
756                name2=item+".npts"
757                name3=item+".nsigmas"
758                if not self.model.details.has_key(name1):
759                    self.model.details [name1] = ["",None,None] 
760
761                iy += 1
762                for p in self.model.dispersion[item].keys(): 
763       
764                    if p=="width":
765                        ix = 0
766                        cb = wx.CheckBox(self, -1, name0, (10, 10))
767                        cb.SetValue(CHECK_STATE)
768                        cb.SetToolTipString("Check mark to fit")
769                        wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
770                        self.sizer4_4.Add( cb,( iy, ix),(1,1), 
771                                           wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 
772                                           5)
773                        ix = 1
774                        value= self.model.getParam(name1)
775                        ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.3,20)
776                                                  ,style=wx.TE_PROCESS_ENTER)
777                        ctl1.SetLabel('PD[ratio]')
778                        poly_text = "Polydispersity (STD/mean) of %s\n" % item
779                        poly_text += "STD: the standard deviation"
780                        poly_text += " from the mean value."
781                        ctl1.SetToolTipString(poly_text)
782                        ctl1.SetValue(str (format_number(value, True)))
783                        self.sizer4_4.Add(ctl1, (iy,ix),(1,1),wx.EXPAND)
784                        ## text to show error sign
785                        ix = 2
786                        text2=wx.StaticText(self, -1, '+/-')
787                        self.sizer4_4.Add(text2,(iy, ix),(1,1),
788                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
789                        if not self.is_mac:
790                            text2.Hide() 
791
792                        ix = 3
793                        ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.3,20), 
794                                           style=0)
795                 
796                        self.sizer4_4.Add(ctl2, (iy,ix),(1,1), 
797                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
798                        if not self.is_mac:
799                            ctl2.Hide()
800
801                        ix = 4
802                        ctl3 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/2,
803                                            20), style=wx.TE_PROCESS_ENTER,
804                                text_enter_callback = self._onparamRangeEnter)
805                       
806                        self.sizer4_4.Add(ctl3, (iy,ix),(1,1), 
807                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
808                       
809               
810                        ix = 5
811                        ctl4 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/2,
812                                            20), style=wx.TE_PROCESS_ENTER,
813                            text_enter_callback = self._onparamRangeEnter)
814                       
815                        self.sizer4_4.Add(ctl4, (iy,ix),(1,1), 
816                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
817
818                       
819                        if self.engine_type=="park":
820                            ctl3.Show(True)
821                            ctl4.Show(True)
822                                                             
823                    elif p=="npts":
824                            ix = 6
825                            value= self.model.getParam(name2)
826                            Tctl = self.ModelTextCtrl(self, -1, 
827                                                       size=(_BOX_WIDTH/2.2,20),
828                                                style=wx.TE_PROCESS_ENTER)
829                           
830                            Tctl.SetValue(str (format_number(value)))
831                            self.sizer4_4.Add(Tctl, (iy,ix),(1,1),
832                                               wx.EXPAND|wx.ADJUST_MINSIZE, 0)
833                            self.fixed_param.append([None,name2, Tctl,None,None,
834                                                      None, None,None])
835                    elif p=="nsigmas":
836                            ix = 7
837                            value= self.model.getParam(name3)
838                            Tct2 = self.ModelTextCtrl(self, -1, 
839                                                      size=(_BOX_WIDTH/2.2,20),
840                                                style=wx.TE_PROCESS_ENTER)
841                           
842                            Tct2.SetValue(str (format_number(value)))
843                            self.sizer4_4.Add(Tct2, (iy,ix),(1,1),
844                                               wx.EXPAND|wx.ADJUST_MINSIZE, 0)
845                            #ix +=1
846                            #self.sizer4_4.Add((20,20), (iy,ix),(1,1),
847                            #                   wx.EXPAND|wx.ADJUST_MINSIZE, 0)
848                           
849                            self.fixed_param.append([None, name3, Tct2,
850                                                     None, None, None, 
851                                                     None, None])
852
853
854                ix = 8     
855                disp_box = wx.ComboBox(self, -1,size=(65,-1),
856                                style=wx.CB_READONLY, name = '%s'% name1)
857                for key, value in self.polydisp.iteritems():
858                    name_disp = str(key)
859                    disp_box.Append(name_disp,value)
860                    disp_box.SetStringSelection("gaussian") 
861                wx.EVT_COMBOBOX(disp_box,-1, self._on_disp_func)     
862                self.sizer4_4.Add(disp_box,(iy,ix),(1,1), wx.EXPAND)
863                self.fittable_param.append([cb,name1,ctl1,text2,
864                                                    ctl2, ctl3, ctl4, disp_box])                     
865                           
866        ix =0
867        iy +=1 
868        self.sizer4_4.Add((20,20),(iy,ix),(1,1), 
869                          wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15) 
870        first_orient  = True
871        for item in self.model.dispersion.keys():
872            if  item in self.model.orientation_params:
873                if not self.disp_cb_dict.has_key(item):
874                    self.disp_cb_dict[item]= None
875                name0="Distribution of " + item
876                name1=item+".width"
877                name2=item+".npts"
878                name3=item+".nsigmas"
879               
880                if not self.model.details.has_key(name1):
881                    self.model.details [name1] = ["",None,None]                 
882 
883
884                iy += 1
885                for p in self.model.dispersion[item].keys(): 
886       
887                    if p=="width":
888                        ix = 0
889                        cb = wx.CheckBox(self, -1, name0, (10, 10))
890                        cb.SetValue(CHECK_STATE)
891                        cb.SetToolTipString("Check mark to fit")
892                        wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
893                        self.sizer4_4.Add( cb,( iy, ix),(1,1), 
894                                           wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 
895                                           5)
896                        if self.data.__class__.__name__ ==  "Data2D" or \
897                                    self.enable2D:
898                            cb.Show(True)
899                        elif cb.IsShown():
900                            cb.Hide()
901                        ix = 1
902                        value= self.model.getParam(name1)
903                        ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.3,
904                                        20),style=wx.TE_PROCESS_ENTER)
905                        poly_tip = "Absolute Sigma for %s." % item
906                        ctl1.SetToolTipString(poly_tip)
907                        ctl1.SetValue(str (format_number(value, True)))
908                        if self.data.__class__.__name__ == "Data2D" or \
909                                    self.enable2D:
910                            if first_orient:
911                                values.SetLabel('PD[ratio], Sig[deg]')
912                                poly_text = "PD(polydispersity for lengths):\n"
913                                poly_text +=  "It should be a value between"
914                                poly_text +=  "0 and 1\n"
915                                poly_text += "Sigma for angles: \n"
916                                poly_text += "It is the STD (ratio*mean)"
917                                poly_text += " of the distribution.\n "
918                           
919                                values.SetToolTipString(poly_text)
920                                first_orient = False
921                            ctl1.Show(True)
922                        elif ctl1.IsShown():
923                            ctl1.Hide()
924                       
925                        self.sizer4_4.Add(ctl1, (iy,ix),(1,1),wx.EXPAND)
926                        ## text to show error sign
927                        ix = 2
928                        text2=wx.StaticText(self, -1, '+/-')
929                        self.sizer4_4.Add(text2,(iy, ix),(1,1),
930                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
931
932                        text2.Hide() 
933
934                        ix = 3
935                        ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.3,20), 
936                                           style=0)
937                   
938                        self.sizer4_4.Add(ctl2, (iy,ix),(1,1), 
939                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
940
941                        ctl2.Hide()
942                        if self.data.__class__.__name__ ==  "Data2D" or \
943                                self.enable2D:
944                            if self.is_mac: 
945                                text2.Show(True)
946                                ctl2.Show(True) 
947                           
948                        ix = 4
949                        ctl3 = self.ModelTextCtrl(self, -1, 
950                                                  size=(_BOX_WIDTH/2,20), 
951                                                  style=wx.TE_PROCESS_ENTER,
952                                text_enter_callback = self._onparamRangeEnter)
953
954                        self.sizer4_4.Add(ctl3, (iy,ix),(1,1), 
955                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
956
957                        ctl3.Hide()
958               
959                        ix = 5
960                        ctl4 = self.ModelTextCtrl(self, -1, 
961                            size=(_BOX_WIDTH/2,20), style=wx.TE_PROCESS_ENTER,
962                                text_enter_callback = self._onparamRangeEnter)
963                        self.sizer4_4.Add(ctl4, (iy,ix),(1,1), 
964                                          wx.EXPAND|wx.ADJUST_MINSIZE, 0)
965                        ctl4.Hide()
966                       
967                        if self.data.__class__.__name__ ==  "Data2D" or \
968                                self.enable2D:
969                            ctl3.Show(True)
970                            ctl4.Show(True) 
971                             
972                    elif p=="npts":
973                            ix = 6
974                            value= self.model.getParam(name2)
975                            Tctl = self.ModelTextCtrl(self, -1, 
976                                                     size=(_BOX_WIDTH/2.2, 20),
977                                                style=wx.TE_PROCESS_ENTER)
978                           
979                            Tctl.SetValue(str (format_number(value)))
980                            if self.data.__class__.__name__ ==  "Data2D" or \
981                                    self.enable2D:
982                                Tctl.Show(True)
983                            else:
984                                Tctl.Hide()
985                            self.sizer4_4.Add(Tctl, (iy,ix),(1,1),
986                                               wx.EXPAND|wx.ADJUST_MINSIZE, 0)
987                            self.fixed_param.append([None,name2, Tctl,None,None,
988                                                      None, None,None])
989                            self.orientation_params_disp.append([None,name2, 
990                                                    Tctl,None,None,
991                                                      None, None,None])
992                    elif p=="nsigmas":
993                            ix = 7
994                            value= self.model.getParam(name3)
995                            Tct2 = self.ModelTextCtrl(self, -1, 
996                                                      size=(_BOX_WIDTH/2.2,20),
997                                                style=wx.TE_PROCESS_ENTER)
998                           
999                            Tct2.SetValue(str (format_number(value)))
1000                            if self.data.__class__.__name__ ==  "Data2D" or \
1001                                    self.enable2D:
1002                                Tct2.Show(True)
1003                            else:
1004                                Tct2.Hide()
1005                            self.sizer4_4.Add(Tct2, (iy,ix),(1,1),
1006                                               wx.EXPAND|wx.ADJUST_MINSIZE, 0)
1007
1008
1009                            self.fixed_param.append([None,name3, Tct2
1010                                                 ,None,None, None, None,None])   
1011                                                       
1012                            self.orientation_params_disp.append([None,name3, 
1013                                            Tct2 ,None,None, None, None,None])
1014       
1015
1016                ix = 8     
1017                disp_box = wx.ComboBox(self, -1,size=(65,-1),
1018                                style=wx.CB_READONLY, name = '%s'% name1)
1019                for key, value in self.polydisp.iteritems():
1020                    name_disp = str(key)
1021                    disp_box.Append(name_disp,value)
1022                    disp_box.SetStringSelection("gaussian") 
1023                wx.EVT_COMBOBOX(disp_box,-1, self._on_disp_func)     
1024                self.sizer4_4.Add(disp_box,(iy,ix),(1,1), wx.EXPAND)
1025                self.fittable_param.append([cb,name1,ctl1,text2,
1026                                            ctl2, ctl3, ctl4, disp_box])
1027                self.orientation_params_disp.append([cb,name1,ctl1,
1028                                            text2, ctl2, ctl3, ctl4, disp_box])
1029                       
1030                if self.data.__class__.__name__ == "Data2D" or \
1031                                self.enable2D:
1032                    disp_box.Show(True)
1033                else:
1034                    disp_box.Hide()
1035       
1036
1037        self.state.disp_cb_dict = copy.deepcopy(self.disp_cb_dict) 
1038         
1039        self.state.model = self.model.clone() 
1040         ## save state into
1041        self.state.cb1 = self.cb1.GetValue()
1042        self._copy_parameters_state(self.parameters, self.state.parameters)
1043        self._copy_parameters_state(self.orientation_params_disp,
1044                                     self.state.orientation_params_disp)
1045        self._copy_parameters_state(self.fittable_param, 
1046                                    self.state.fittable_param)
1047        self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
1048
1049
1050        wx.PostEvent(self.parent, StatusEvent(status=\
1051                        " Selected Distribution: Gaussian"))   
1052        #Fill the list of fittable parameters
1053        #self.select_all_param(event=None)
1054        self.get_all_checked_params()
1055        self.Layout()
1056
1057   
1058    def _onDraw(self, event):
1059        """
1060        Update and Draw the model
1061        """ 
1062        if self.model ==None:
1063            msg="Please select a Model first..."
1064            wx.MessageBox(msg, 'Info')
1065
1066            return
1067        """
1068        if not self.data.is_data:
1069            self.npts_x = self.Npts_total.GetValue()
1070            self.Npts_fit.SetValue(self.npts_x)
1071            self.create_default_data()
1072        """
1073        flag = self._update_paramv_on_fit() 
1074   
1075        wx.CallAfter(self._onparamEnter_helper)
1076        if not flag:
1077            msg= "The parameters are invalid"
1078            wx.PostEvent(self.parent.parent, StatusEvent(status= msg ))
1079            return 
1080       
1081
1082    def _onFit(self, event):     
1083        """
1084        Allow to fit
1085        """
1086        if self.fit_started:
1087            self._StopFit()
1088           
1089        if event != None:
1090            event.Skip() 
1091        if len(self.parent._manager.fit_thread_list)>0 and\
1092                    self.parent._manager._fit_engine != "park" and\
1093                    self._manager.sim_page != None and \
1094                    self._manager.sim_page.uid == self.uid: 
1095            msg = "The FitEnging will be set to 'ParkMC'\n"
1096            msg += " to fit with more than one data set..."
1097            wx.MessageBox(msg, 'Info')
1098            #wx.PostEvent(self._manager.parent, StatusEvent(status=\
1099            #                "Fitting: %s"%msg))
1100           
1101        if self.data is None:
1102            msg = "Please get Data first..."
1103            wx.MessageBox(msg, 'Info')
1104            wx.PostEvent(self._manager.parent, StatusEvent(status=\
1105                            "Fit: %s" % msg))
1106            return
1107        if self.model is None:
1108            msg = "Please select a Model first..."
1109            wx.MessageBox(msg, 'Info')
1110            wx.PostEvent(self._manager.parent, StatusEvent(status=\
1111                            "Fit: %s"%msg, type="stop"))
1112            return
1113       
1114        if len(self.param_toFit) <= 0:
1115            msg= "Select at least one parameter to fit"
1116            wx.MessageBox(msg, 'Info')
1117            wx.PostEvent(self._manager.parent, StatusEvent(status= msg, 
1118                                                         type="stop" ))
1119            return 
1120       
1121        flag = self._update_paramv_on_fit() 
1122       
1123        if self.batch_on and not self._is_2D():
1124            if not self._validate_Npts_1D():
1125                return
1126               
1127        if not flag:
1128            msg= "Fitting range or parameters are invalid"
1129            wx.PostEvent(self.parent.parent, StatusEvent(status= msg, 
1130                                                         type="stop"))
1131            return 
1132             
1133        self.select_param(event =None)
1134       
1135        #Clear errors if exist from previous fitting
1136        #self._clear_Err_on_Fit()
1137
1138        # Remove or do not allow fitting on the Q=0 point, especially
1139        # when y(q=0)=None at x[0].         
1140        self.qmin_x = float(self.qmin.GetValue())
1141        self.qmax_x = float(self.qmax.GetValue())
1142        self._manager._reset_schedule_problem(value=0, uid=self.uid)
1143        self._manager.schedule_for_fit(uid=self.uid,value=1) 
1144        self._manager.set_fit_range(uid=self.uid,qmin=self.qmin_x, 
1145                                   qmax=self.qmax_x)
1146
1147        #single fit
1148        self._manager.onFit(uid=self.uid)
1149        self.fit_started = True
1150        self.btFit.SetLabel("Stop")
1151           
1152    def get_weight_flag(self):
1153        """
1154        Get flag corresponding to a given weighting dI data.
1155        """
1156        button_list = [self.dI_noweight,
1157                       self.dI_didata,
1158                       self.dI_sqrdata,
1159                       self.dI_idata]
1160        flag = 1
1161        for item in button_list:
1162            if item.GetValue():
1163                if button_list.index(item) == 0:
1164                    flag = 0 #dy = numpy.ones_like(dy_data)
1165                elif button_list.index(item) == 1:
1166                    flag = 1 #dy = dy_data
1167                elif button_list.index(item) == 2:
1168                    flag = 2 #dy = numpy.sqrt(numpy.abs(data))
1169                elif button_list.index(item) == 3:
1170                    flag = 3 # dy = numpy.abs(data)
1171                break
1172        return flag
1173               
1174    def bind_fit_button(self):
1175        """
1176        bind the fit button to either fit handler or stop fit handler
1177        """
1178        self.btFit.Unbind(event=wx.EVT_BUTTON, id= self.btFit.GetId())
1179        if self.btFit.GetLabel().lower() == "stop":
1180            self.fit_started = True
1181            self.btFit.SetForegroundColour('red')
1182            self.btFit.Bind(event=wx.EVT_BUTTON, handler=self._StopFit,
1183                             id=self.btFit.GetId())
1184        elif self.btFit.GetLabel().lower() == "fit":
1185            self.fit_started = False
1186            self.btFit.SetDefault()
1187            self.btFit.SetForegroundColour('black')
1188            #self.btFit.SetBackgroundColour(self.default_bt_colour)
1189            self.btFit.Bind(event=wx.EVT_BUTTON, handler=self._onFit, 
1190                            id=self.btFit.GetId())
1191        else:
1192            msg = "FitPage: fit button has unknown label"
1193            raise ValuerError, msg
1194        self._manager._reset_schedule_problem(value=0)
1195         
1196    def is_fitting(self):
1197        if self.fit_started:
1198            self._StopFit(event=None)
1199           
1200    def _StopFit(self, event=None):
1201        """
1202        Stop fit
1203        """
1204        #time.sleep(0.1)
1205        if event != None:
1206            event.Skip()
1207        #if self.engine_type=="scipy":
1208        self._manager.stop_fit(self.uid)
1209        self._manager._reset_schedule_problem(value=0)
1210        self._on_fit_complete()
1211         
1212    def rename_model(self):
1213        """
1214        find a short name for model
1215        """
1216        if self.model is not None:
1217            self.model.name = "M" + str(self.index_model)
1218   
1219
1220    def _on_select_model(self, event=None): 
1221        """
1222        call back for model selection
1223        """ 
1224        self.Show(False) 
1225        copy_flag = False 
1226        is_poly_enabled = None 
1227        if event != None:
1228            if (event.GetEventObject() == self.formfactorbox\
1229                        and self.structurebox.GetValue() != 'None')\
1230                        or event.GetEventObject() == self.structurebox\
1231                        or event.GetEventObject() == self.multifactorbox:
1232                copy_flag = self.get_copy_params()
1233                is_poly_enabled = self.enable_disp.GetValue() 
1234
1235        self._on_select_model_helper() 
1236        self.set_model_param_sizer(self.model)                   
1237        if self.model is None:
1238            self._set_bookmark_flag(False)
1239            self._keep.Enable(False)
1240            self._set_save_flag(False)
1241        self.enable_disp.SetValue(False)
1242        self.disable_disp.SetValue(True)
1243        try:
1244            self.set_dispers_sizer()
1245        except:
1246            pass
1247        #self.btFit.SetFocus()
1248        self.state.enable_disp = self.enable_disp.GetValue()
1249        self.state.disable_disp = self.disable_disp.GetValue()
1250        self.state.pinhole_smearer = self.pinhole_smearer.GetValue()
1251        self.state.slit_smearer = self.slit_smearer.GetValue()
1252   
1253        self.state.structurecombobox = self.structurebox.GetValue()#.GetCurrentSelection()
1254        self.state.formfactorcombobox = self.formfactorbox.GetValue()#.GetCurrentSelection()
1255        self.enable_fit_button()
1256        if self.model != None:
1257            self.rename_model()
1258            self._set_copy_flag(True)
1259            self._set_paste_flag(True)
1260            if self.data is not None:
1261                is_data = check_data_validity(self.data) 
1262                if is_data:
1263                    self._set_bookmark_flag(True)
1264                    self._keep.Enable(True)
1265                    self._set_save_flag(True)
1266            # Reset smearer, model and data
1267            if not copy_flag:
1268                self.disable_smearer.SetValue(True)
1269                self.enable_smearer.SetValue(False)
1270   
1271            # more disables for 2D
1272            self._set_smear_buttons()
1273           
1274            try:
1275                # update smearer sizer
1276                #if not self.enable_smearer.GetValue():
1277                #    self.disable_smearer.SetValue(True)
1278                self.onSmear(None)
1279                temp_smear = None
1280                if not self.disable_smearer.GetValue():
1281                    # Set the smearer environments
1282                    temp_smear = self.current_smearer
1283            except:
1284                raise
1285                ## error occured on chisqr computation
1286                #pass
1287            ## event to post model to fit to fitting plugins
1288            (ModelEventbox, EVT_MODEL_BOX) = wx.lib.newevent.NewEvent()
1289         
1290            ## set smearing value whether or not
1291            #    the data contain the smearing info
1292            evt = ModelEventbox(model=self.model, 
1293                                        smearer=temp_smear, 
1294                             enable_smearer=not self.disable_smearer.GetValue(),
1295                                        qmin=float(self.qmin_x),
1296                                        uid=self.uid,
1297                                        caption=self.window_caption,
1298                                     qmax=float(self.qmax_x)) 
1299   
1300            self._manager._on_model_panel(evt=evt)
1301            self.mbox_description.SetLabel("Model [ %s ]" % str(self.model.name))
1302            self.mbox_description.SetForegroundColour(wx.BLUE)
1303            self.state.model = self.model.clone()
1304            self.state.model.name = self.model.name
1305
1306           
1307        if event != None:
1308            ## post state to fit panel
1309            new_event = PageInfoEvent(page = self)
1310            wx.PostEvent(self.parent, new_event) 
1311            #update list of plugins if new plugin is available
1312            if self.plugin_rbutton.GetValue():
1313                temp = self.parent.update_model_list()
1314                if temp:
1315                    self.model_list_box = temp
1316                    current_val = self.formfactorbox.GetValue()
1317                    pos = self.formfactorbox.GetSelection()
1318                    self._show_combox_helper()
1319                    self.formfactorbox.SetSelection(pos)
1320                    self.formfactorbox.SetValue(current_val)
1321            # when select a model only from guictr/button
1322            if is_poly_enabled != None:
1323                self.enable_disp.SetValue(is_poly_enabled)
1324                self.disable_disp.SetValue(not is_poly_enabled)
1325                self._set_dipers_Param(event=None)
1326                self.state.enable_disp = self.enable_disp.GetValue()
1327                self.state.disable_disp = self.disable_disp.GetValue()
1328
1329            # Keep the previous param values
1330            if copy_flag:
1331                self.get_paste_params(copy_flag)
1332                wx.CallAfter(self._onDraw, None)
1333               
1334        else:
1335            self._draw_model()
1336           
1337        if self.batch_on:
1338            self.slit_smearer.Enable(False)
1339            self.pinhole_smearer.Enable(False)
1340            self.btEditMask.Disable() 
1341            self.EditMask_title.Disable()
1342           
1343        self.Show(True)     
1344        self.SetupScrolling()
1345         
1346    def _onparamEnter(self,event):
1347        """
1348        when enter value on panel redraw model according to changed
1349        """
1350        if self.model ==None:
1351            msg="Please select a Model first..."
1352            wx.MessageBox(msg, 'Info')
1353            #wx.PostEvent(self._manager.parent, StatusEvent(status=\
1354            #                "Parameters: %s"%msg))
1355            return
1356
1357        #default flag
1358        flag = False
1359        self.fitrange = True
1360        #get event object
1361        tcrtl = event.GetEventObject()
1362       
1363        #wx.PostEvent(self._manager.parent, StatusEvent(status=" \
1364        #                        updating ... ",type="update"))
1365        #Clear msg if previously shown.
1366        msg= ""
1367        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1368
1369        if check_float(tcrtl):
1370            flag = self._onparamEnter_helper() 
1371            self.set_npts2fit() 
1372            if self.fitrange:             
1373                temp_smearer = None
1374                if not self.disable_smearer.GetValue():
1375                    temp_smearer= self.current_smearer
1376                    ## set smearing value whether or not
1377                    #        the data contain the smearing info
1378                    if self.slit_smearer.GetValue():
1379                        flag1 = self.update_slit_smear()
1380                        flag = flag or flag1
1381                    elif self.pinhole_smearer.GetValue():
1382                        flag1 = self.update_pinhole_smear()
1383                        flag = flag or flag1
1384                elif self.data.__class__.__name__ !=  "Data2D" and \
1385                        not self.enable2D:
1386                    self._manager.set_smearer(smearer=temp_smearer, 
1387                                              fid=self.data.id,
1388                                              uid=self.uid,
1389                                             qmin= float(self.qmin_x),
1390                                            qmax= float(self.qmax_x),
1391                            enable_smearer=not self.disable_smearer.GetValue(),
1392                                            draw=True) 
1393                if flag:   
1394                    #self.compute_chisqr(smearer= temp_smearer)
1395       
1396                    ## new state posted
1397                    if self.state_change:
1398                        #self._undo.Enable(True)
1399                        event = PageInfoEvent(page = self)
1400                        wx.PostEvent(self.parent, event)
1401                    self.state_change= False 
1402            else: 
1403                # invalid fit range: do nothing here:
1404                # msg already displayed in validate
1405                return     
1406        else:
1407            self.save_current_state()
1408            msg= "Cannot Plot :Must enter a number!!!  "
1409            wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1410             
1411        self.save_current_state() 
1412        return 
1413   
1414    def _onparamRangeEnter(self, event):
1415        """
1416        Check validity of value enter in the parameters range field
1417        """
1418        #if self.check_invalid_panel():
1419        #    return
1420        tcrtl= event.GetEventObject()
1421        #Clear msg if previously shown.
1422        msg= ""
1423        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1424        # Flag to register when a parameter has changed.
1425        is_modified = False
1426        if tcrtl.GetValue().lstrip().rstrip()!="":
1427            try:
1428                value = float(tcrtl.GetValue())
1429                tcrtl.SetBackgroundColour(wx.WHITE)
1430                self._check_value_enter(self.fittable_param ,is_modified)
1431                self._check_value_enter(self.parameters ,is_modified) 
1432            except:
1433                tcrtl.SetBackgroundColour("pink")
1434                msg= "Model Error:wrong value entered : %s"% sys.exc_value
1435                wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1436                return 
1437        else:
1438           tcrtl.SetBackgroundColour(wx.WHITE)
1439           
1440        #self._undo.Enable(True)
1441        self.save_current_state()
1442        event = PageInfoEvent(page = self)
1443        wx.PostEvent(self.parent, event)
1444        self.state_change= False
1445
1446                   
1447    def _onQrangeEnter(self, event):
1448        """
1449        Check validity of value enter in the Q range field
1450        """
1451        tcrtl = event.GetEventObject()
1452        #Clear msg if previously shown.
1453        msg= ""
1454        wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1455        # For theory mode
1456        if not self.data.is_data:
1457            self.npts_x = self.Npts_total.GetValue()
1458            self.Npts_fit.SetValue(self.npts_x)
1459            self.create_default_data()
1460        # Flag to register when a parameter has changed.
1461        is_modified = False
1462        if tcrtl.GetValue().lstrip().rstrip()!="":
1463            try:
1464                value = float(tcrtl.GetValue())
1465                tcrtl.SetBackgroundColour(wx.WHITE)
1466
1467                # If qmin and qmax have been modified, update qmin and qmax
1468                if self._validate_qrange( self.qmin, self.qmax):
1469                    tempmin = float(self.qmin.GetValue())
1470                    if tempmin != self.qmin_x:
1471                        self.qmin_x = tempmin
1472                    tempmax = float(self.qmax.GetValue())
1473                    if tempmax != self.qmax_x:
1474                        self.qmax_x = tempmax
1475                else:
1476                    tcrtl.SetBackgroundColour("pink")
1477                    msg= "Model Error:wrong value entered : %s"% sys.exc_value
1478                    wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1479                    return 
1480               
1481            except:
1482                tcrtl.SetBackgroundColour("pink")
1483                msg= "Model Error:wrong value entered : %s"% sys.exc_value
1484                wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1485                return 
1486            #Check if # of points for theory model are valid(>0).
1487            # check for 2d
1488            if self.data.__class__.__name__ == "Data2D" or \
1489                    self.enable2D:
1490                # set mask   
1491                radius= numpy.sqrt(self.data.qx_data*self.data.qx_data + 
1492                                    self.data.qy_data*self.data.qy_data )
1493                index_data = ((self.qmin_x <= radius)& \
1494                                (radius<= self.qmax_x))
1495                index_data = (index_data)&(self.data.mask)
1496                index_data = (index_data)&(numpy.isfinite(self.data.data))
1497                if len(index_data[index_data]) < 10:
1498                    msg = "Cannot Plot :No or too little npts in"
1499                    msg += " that data range!!!  "
1500                    wx.PostEvent(self.parent.parent, 
1501                                 StatusEvent(status=msg))
1502                    return
1503                else:
1504                    #self.data.mask = index_data
1505                    #self.Npts_fit.SetValue(str(len(self.data.mask)))
1506                    self.set_npts2fit() 
1507            else:
1508                index_data = ((self.qmin_x <= self.data.x)& \
1509                              (self.data.x <= self.qmax_x))
1510                self.Npts_fit.SetValue(str(len(self.data.x[index_data])))
1511           
1512            self.npts_x = self.Npts_total.GetValue()
1513            self.create_default_data()
1514            self._save_plotting_range()
1515        else:
1516           tcrtl.SetBackgroundColour("pink")
1517           msg= "Model Error:wrong value entered!!!"
1518           wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1519       
1520        self._draw_model()
1521        self.save_current_state()
1522        event = PageInfoEvent(page = self)
1523        wx.PostEvent(self.parent, event)
1524        self.state_change= False
1525        return
1526   
1527    def _clear_Err_on_Fit(self):
1528        """
1529        hide the error text control shown
1530        after fitting
1531        """
1532        if self.is_mac:
1533            return
1534        if hasattr(self,"text2_3"):
1535            self.text2_3.Hide()
1536
1537        if len(self.parameters)>0:
1538            for item in self.parameters:
1539                #Skip t ifhe angle parameters if 1D data
1540                if self.data.__class__.__name__ !=  "Data2D" and \
1541                        not self.enable2D:
1542                    if item in self.orientation_params:
1543                        continue
1544                if item in self.param_toFit:
1545                    continue
1546                ## hide statictext +/-   
1547                if len(item) < 4 :
1548                    continue
1549                if item[3]!=None and item[3].IsShown():
1550                    item[3].Hide()
1551                ## hide textcrtl  for error after fit
1552                if item[4]!=None and item[4].IsShown():                   
1553                    item[4].Hide()
1554 
1555        if len(self.fittable_param)>0:
1556            for item in self.fittable_param:
1557                #Skip t ifhe angle parameters if 1D data
1558                if self.data.__class__.__name__ !=  "Data2D" and \
1559                        not self.enable2D:
1560                    if item in self.orientation_params:
1561                        continue
1562                if item in self.param_toFit:
1563                    continue
1564                if len(item) < 4 :
1565                    continue
1566                ## hide statictext +/-   
1567                if item[3]!=None and item[3].IsShown():
1568                    item[3].Hide()
1569                ## hide textcrtl  for error after fit
1570                if item[4]!=None and item[4].IsShown():
1571                    item[4].Hide()
1572        return       
1573               
1574    def _get_defult_custom_smear(self):
1575        """
1576        Get the defult values for custum smearing.
1577        """
1578        # get the default values
1579        if self.dxl == None: self.dxl = 0.0
1580        if self.dxw == None: self.dxw = ""
1581        if self.dx_min == None: self.dx_min = SMEAR_SIZE_L
1582        if self.dx_max == None: self.dx_max = SMEAR_SIZE_H
1583       
1584    def _get_smear_info(self):
1585        """
1586        Get the smear info from data.
1587       
1588        :return: self.smear_type, self.dq_l and self.dq_r,
1589            respectively the type of the smear, dq_min and
1590            dq_max for pinhole smear data
1591            while dxl and dxw for slit smear
1592           
1593        """
1594
1595        # default
1596        self.smear_type = None
1597        self.dq_l = None
1598        self.dq_r = None
1599        data = self.data
1600        if self.data is None:
1601            return
1602        elif self.data.__class__.__name__ ==  "Data2D" or \
1603                        self.enable2D: 
1604            if data.dqx_data == None or  data.dqy_data ==None: 
1605                return
1606            elif self.current_smearer != None and (data.dqx_data.any()!=0) and \
1607                            (data.dqx_data.any()!=0): 
1608                self.smear_type = "Pinhole2d"
1609                self.dq_l = format_number(numpy.average(data.dqx_data)) 
1610                self.dq_r = format_number(numpy.average(data.dqy_data)) 
1611                return 
1612            else: 
1613                return
1614        # check if it is pinhole smear and get min max if it is.
1615        if data.dx != None and all(data.dx !=0): 
1616            self.smear_type = "Pinhole" 
1617            self.dq_l = data.dx[0]
1618            self.dq_r = data.dx[-1]
1619           
1620        # check if it is slit smear and get min max if it is.
1621        elif data.dxl != None or data.dxw != None: 
1622            self.smear_type = "Slit" 
1623            if data.dxl != None and all(data.dxl !=0):
1624                self.dq_l = data.dxl[0]
1625            if data.dxw != None and all(data.dxw !=0): 
1626                self.dq_r = data.dxw[0]   
1627        #return self.smear_type,self.dq_l,self.dq_r       
1628   
1629    def _show_smear_sizer(self):   
1630        """
1631        Show only the sizers depending on smear selection
1632        """
1633        # smear disabled
1634        if self.disable_smearer.GetValue():
1635            self.smear_description_none.Show(True)
1636        # 2Dsmear
1637        elif self._is_2D():
1638            self.smear_description_accuracy_type.Show(True)
1639            self.smear_accuracy.Show(True)
1640            self.smear_description_accuracy_type.Show(True)
1641            self.smear_description_2d.Show(True)
1642            self.smear_description_2d_x.Show(True)
1643            self.smear_description_2d_y.Show(True)
1644            if self.pinhole_smearer.GetValue():
1645                self.smear_pinhole_min.Show(True)
1646                self.smear_pinhole_max.Show(True)
1647        # smear from data
1648        elif self.enable_smearer.GetValue():
1649
1650            self.smear_description_dqdata.Show(True)
1651            if self.smear_type != None:
1652                self.smear_description_smear_type.Show(True)
1653                if self.smear_type == 'Slit':
1654                    self.smear_description_slit_height.Show(True)
1655                    self.smear_description_slit_width.Show(True)               
1656                elif self.smear_type == 'Pinhole':
1657                    self.smear_description_pin_min.Show(True)
1658                    self.smear_description_pin_max.Show(True)
1659                self.smear_description_smear_type.Show(True)
1660                self.smear_description_type.Show(True)
1661                self.smear_data_left.Show(True)
1662                self.smear_data_right.Show(True)
1663        # custom pinhole smear
1664        elif self.pinhole_smearer.GetValue():
1665            if self.smear_type == 'Pinhole':
1666                self.smear_message_new_p.Show(True)
1667                self.smear_description_pin_min.Show(True)
1668                self.smear_description_pin_max.Show(True)
1669
1670            self.smear_pinhole_min.Show(True)
1671            self.smear_pinhole_max.Show(True)
1672        # custom slit smear
1673        elif self.slit_smearer.GetValue():
1674            self.smear_message_new_s.Show(True)
1675            self.smear_description_slit_height.Show(True)
1676            self.smear_slit_height.Show(True)
1677            self.smear_description_slit_width.Show(True)
1678            self.smear_slit_width.Show(True)
1679
1680    def _hide_all_smear_info(self):
1681        """
1682        Hide all smearing messages in the set_smearer sizer
1683        """
1684        self.smear_description_none.Hide()
1685        self.smear_description_dqdata.Hide()
1686        self.smear_description_type.Hide()
1687        self.smear_description_smear_type.Hide()
1688        self.smear_description_accuracy_type.Hide()
1689        self.smear_description_2d_x.Hide()
1690        self.smear_description_2d_y.Hide()
1691        self.smear_description_2d.Hide()
1692       
1693        self.smear_accuracy.Hide()
1694        self.smear_data_left.Hide()
1695        self.smear_data_right.Hide()
1696        self.smear_description_pin_min.Hide()
1697        self.smear_pinhole_min.Hide()
1698        self.smear_description_pin_max.Hide()
1699        self.smear_pinhole_max.Hide()
1700        self.smear_description_slit_height.Hide()
1701        self.smear_slit_height.Hide()
1702        self.smear_description_slit_width.Hide()
1703        self.smear_slit_width.Hide()
1704        self.smear_message_new_p.Hide()
1705        self.smear_message_new_s.Hide()
1706   
1707    def _set_accuracy_list(self):
1708        """
1709        Set the list of an accuracy in 2D custum smear: Xhigh, High, Med, or Low
1710        """
1711        # list of accuracy choices
1712        list = ['Low','Med','High','Xhigh']
1713        for idx in range(len(list)):
1714            self.smear_accuracy.Append(list[idx],idx)
1715           
1716    def _set_fun_box_list(self,fun_box):
1717        """
1718        Set the list of func for multifunctional models
1719        """
1720        # Check if it is multi_functional model
1721        if self.model.__class__ not in self.model_list_box["Multi-Functions"] \
1722                and not self.temp_multi_functional:
1723            return None
1724        # Get the func name list
1725        list = self.model.fun_list
1726        if len(list) == 0:
1727            return None
1728        # build function (combo)box
1729        ind = 0
1730        while(ind < len(list)):
1731            for key, val in list.iteritems():
1732                if (val == ind):
1733                    fun_box.Append(key,val)
1734                    break
1735            ind += 1
1736       
1737    def _on_select_accuracy(self,event):
1738        """
1739        Select an accuracy in 2D custom smear: Xhigh, High, Med, or Low
1740        """
1741        #event.Skip()
1742        # Check if the accuracy is same as before
1743        #self.smear2d_accuracy = event.GetEventObject().GetValue()
1744        self.smear2d_accuracy = self.smear_accuracy.GetValue()
1745        if self.pinhole_smearer.GetValue():
1746            self.onPinholeSmear(event=None)
1747        else:   
1748            self.onSmear(event=None)
1749            if self.current_smearer != None:
1750                self.current_smearer.set_accuracy(accuracy = 
1751                                                  self.smear2d_accuracy) 
1752        event.Skip()
1753
1754    def _on_fun_box(self,event):
1755        """
1756        Select an func: Erf,Rparabola,LParabola
1757        """
1758        fun_val = None
1759        fun_box = event.GetEventObject()
1760        name = fun_box.Name
1761        value = fun_box.GetValue()
1762        if self.model.fun_list.has_key(value):
1763            fun_val = self.model.fun_list[value]
1764
1765        self.model.setParam(name,fun_val)
1766        # save state
1767        self._copy_parameters_state(self.str_parameters, 
1768                                    self.state.str_parameters)
1769        # update params
1770        self._update_paramv_on_fit() 
1771        # draw
1772        self._draw_model()
1773        self.Refresh()
1774        # get ready for new event
1775        event.Skip()
1776       
1777    def _onMask(self, event):     
1778        """
1779        Build a panel to allow to edit Mask
1780        """
1781       
1782        from sans.guiframe.local_perspectives.plotting.masking \
1783        import MaskPanel as MaskDialog
1784       
1785        self.panel = MaskDialog(base=self, data=self.data, id=wx.NewId())
1786        #self.panel.Bind(wx.EVT_CLOSE, self._draw_masked_model)
1787        self.panel.ShowModal()
1788        #wx.PostEvent(self.parent, event)
1789       
1790    def _draw_masked_model(self, event):
1791        """
1792        Draw model image w/mask
1793        """
1794        #event.Skip()
1795        is_valid_qrange = self._update_paramv_on_fit()
1796
1797        if is_valid_qrange and self.model != None:
1798            #self.panel.Show(0)
1799            #self.panel.Destroy() # frame
1800            self.panel.MakeModal(False)
1801            event.Skip() 
1802            # try re draw the model plot if it exists
1803            self._draw_model()
1804            self.set_npts2fit()
1805        elif self.model == None:
1806            #self.panel.Show(0)
1807            #self.panel.Destroy() # frame
1808            self.panel.MakeModal(False)
1809            event.Skip()
1810            self.set_npts2fit()
1811            msg= "No model is found on updating MASK in the model plot... "
1812            wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1813        else:
1814            event.Skip()
1815            msg = ' Please consider your Q range, too.'
1816            self.panel.ShowMessage(msg)
1817
1818       
1819    def _set_smear(self, data):
1820        """
1821        """
1822        if data is None:
1823            return
1824        self.current_smearer = smear_selection(data, self.model)
1825        flag =  self.disable_smearer.GetValue()
1826        self.disable_smearer.SetValue(flag)
1827        if self.current_smearer == None:
1828            self.enable_smearer.Disable()
1829        else:
1830            self.enable_smearer.Enable()
1831        if not flag:
1832            self.onSmear(None)
1833
1834    def _mac_sleep(self, sec=0.2):
1835        """
1836        Give sleep to MAC
1837        """
1838        if self.is_mac:
1839            time.sleep(sec)
1840
1841    def get_view_mode(self):
1842        """
1843        return True if the panel allow 2D or False if 1D
1844        """
1845        return self.enable2D
1846   
1847    def compute_data_set_range(self, data_list):
1848        """
1849        find the range that include all data  in the set
1850        return the minimum and the maximum values
1851        """
1852        if data_list is not None and data_list != []:
1853            for data in data_list:
1854                qmin, qmax, npts = self.compute_data_range(data)
1855                self.qmin_data_set = min(self.qmin_data_set, qmin)
1856                self.qmax_data_set = max(self.qmax_data_set, qmax)
1857                self.npts_data_set += npts
1858        return self.qmin_data_set, self.qmax_data_set, self.npts_data_set
1859       
1860    def compute_data_range(self, data):
1861        """
1862        compute the minimum and the maximum range of the data
1863        return the npts contains in data
1864        :param data:
1865        """
1866        qmin, qmax, npts = None, None, None
1867        if data is not None:
1868            if not hasattr(data,"data"):
1869                qmin = min(data.x)
1870                # Maximum value of data 
1871                qmax = max(data.x)
1872                npts = len(data.x)
1873            else:
1874                qmin = 0
1875                x = max(math.fabs(data.xmin), math.fabs(data.xmax)) 
1876                y = max(math.fabs(data.ymin), math.fabs(data.ymax))
1877                ## Maximum value of data 
1878                qmax = math.sqrt(x*x + y*y)
1879                npts = len(data.data)
1880        return qmin, qmax, npts
1881           
1882   
1883    def set_data(self, data):
1884        """
1885        reset the current data
1886        """
1887        from sans.guiframe.dataFitting import check_data_validity
1888        id = None
1889        flag = False
1890        is_data = False
1891        if self.data is not None:
1892            is_data = check_data_validity(self.data) 
1893        if not is_data and data is not None:
1894                flag = True
1895        if data is not None:
1896            id = data.id
1897            if is_data:
1898                self.graph_id = self.data.group_id
1899                flag = (data.id != self.data.id)
1900        self.data = data
1901        if check_data_validity(data):
1902            self.graph_id = data.group_id
1903        self.data.group_id = self.graph_id
1904       
1905        if self.data is None:
1906            data_name = ""
1907            self._set_bookmark_flag(False)
1908            self._keep.Enable(False)
1909            self._set_save_flag(False)
1910        else:
1911            if self.model != None:
1912                self._set_bookmark_flag(True)
1913                self._keep.Enable(True)
1914               
1915            self._set_save_flag(True)
1916            self._set_preview_flag(True)
1917
1918            self._set_smear(data)
1919            # more disables for 2D
1920            if self.data.__class__.__name__ ==  "Data2D" or \
1921                        self.enable2D:
1922                self.slit_smearer.Disable()
1923                self.pinhole_smearer.Enable(True) 
1924                self.default_mask = copy.deepcopy(self.data.mask)
1925                if self.data.err_data == None or\
1926                        (self.data.err_data == 1).all() or\
1927                        (self.data.err_data == 0).all():
1928                    self.dI_didata.Enable(False)
1929                    self.dI_noweight.SetValue(True)
1930                    self.weightbt_string = self.dI_noweight.GetLabelText()
1931                else:
1932                    self.dI_didata.Enable(True)
1933                    self.dI_didata.SetValue(True)
1934                    self.weightbt_string = self.dI_didata.GetLabelText()
1935            else:
1936                self.slit_smearer.Enable(True) 
1937                self.pinhole_smearer.Enable(True) 
1938                if self.data.dy == None or\
1939                     (self.data.dy == 1).all() or\
1940                     (self.data.dy == 0).all():
1941                    self.dI_didata.Enable(False)
1942                    self.dI_noweight.SetValue(True)
1943                    self.weightbt_string = self.dI_noweight.GetLabelText()
1944                else:
1945                    self.dI_didata.Enable(True)
1946                    self.dI_didata.SetValue(True)
1947                    self.weightbt_string = self.dI_didata.GetLabelText()
1948            # Enable weighting radio uttons     
1949            self.dI_noweight.Enable(True)   
1950            self.dI_sqrdata.Enable(True)
1951            self.dI_idata.Enable(True)
1952             
1953            self.formfactorbox.Enable()
1954            self.structurebox.Enable()
1955            data_name = self.data.name
1956            _, _, npts =  self.compute_data_range(self.data)
1957            #set maximum range for x in linear scale
1958            if not hasattr(self.data, "data"): #Display only for 1D data fit
1959                self.btEditMask.Disable() 
1960                self.EditMask_title.Disable()
1961            else:
1962                self.btEditMask.Enable() 
1963                self.EditMask_title.Enable() 
1964   
1965        self.Npts_total.SetValue(str(npts))
1966        #default:number of data points selected to fit
1967        self.Npts_fit.SetValue(str(npts))
1968        self.Npts_total.SetEditable(False)
1969        self.Npts_total.SetBackgroundColour(\
1970                                    self.GetParent().GetBackgroundColour())
1971       
1972        self.Npts_total.Bind(wx.EVT_MOUSE_EVENTS, self._npts_click)
1973        self.dataSource.SetValue(data_name)
1974        self.state.data = data
1975        self.enable_fit_button()
1976        # send graph_id to page_finder
1977        self._manager.set_graph_id(uid=self.uid, graph_id=self.graph_id)
1978        #focus the page
1979        if check_data_validity(data):
1980            self.data_box_description.SetForegroundColour(wx.BLUE)
1981       
1982        if self.batch_on:
1983            self.slit_smearer.Enable(False)
1984            self.pinhole_smearer.Enable(False)
1985            self.btEditMask.Disable() 
1986            self.EditMask_title.Disable()
1987
1988       
1989        self.on_set_focus(None)
1990        self.Refresh()
1991        #update model plot with new data information
1992        if flag:
1993            #set model view button
1994            if not self.enable_smearer.GetValue():
1995                    self.disable_smearer.SetValue(True)
1996            self.onSmear(None)
1997            #self._set_smear(self.data)
1998
1999            if self.data.__class__.__name__ == "Data2D":
2000                self.enable2D = True
2001                self.model_view.SetLabel("2D Mode")
2002            else:
2003                self.enable2D = False
2004                self.model_view.SetLabel("1D Mode")
2005            self.model_view.Disable()
2006            #replace data plot on combo box selection
2007            #by removing the previous selected data
2008            wx.PostEvent(self._manager.parent, NewPlotEvent(action="remove",
2009                                                    group_id=self.graph_id, id=id))
2010            #plot the current selected data
2011            wx.PostEvent(self._manager.parent, NewPlotEvent(plot=self.data, 
2012                                                title=str(self.data.title)))
2013            #self._manager.store_data(uid=self.uid, data_list=self.data_list,
2014            #                          caption=self.window_name)
2015            self._draw_model()
2016   
2017    def _npts_click(self, event):
2018        """
2019        Prevent further handling of the mouse event on Npts_total
2020        by not calling Skip().
2021        """ 
2022        pass
2023   
2024    def reset_page(self, state,first=False):
2025        """
2026        reset the state
2027        """
2028        self.reset_page_helper(state)
2029
2030        self.select_param(event = None) 
2031        #Save state_fit
2032        self.save_current_state_fit()
2033        self._lay_out()
2034        self.Refresh()
2035       
2036    def get_range(self):
2037        """
2038        return the fitting range
2039        """
2040        return float(self.qmin_x) , float(self.qmax_x)
2041   
2042    def get_npts2fit(self):
2043        """
2044        return numbers of data points within qrange
2045       
2046        :Note: This is for Park where chi2 is not normalized by Npts of fit
2047       
2048        """
2049        if self.data is None:
2050            return
2051        npts2fit = 0
2052        qmin,qmax = self.get_range()
2053        if self.data.__class__.__name__ == "Data2D" or \
2054                        self.enable2D:
2055            radius= numpy.sqrt( self.data.qx_data*self.data.qx_data + 
2056                                self.data.qy_data*self.data.qy_data )
2057            index_data = (self.qmin_x <= radius)&(radius<= self.qmax_x)
2058            index_data= (index_data)&(self.data.mask)
2059            index_data = (index_data)&(numpy.isfinite(self.data.data))
2060            npts2fit = len(self.data.data[index_data])
2061        else:
2062            for qx in self.data.x:
2063                   if qx >= qmin and qx <= qmax:
2064                       npts2fit += 1
2065        return npts2fit
2066
2067    def set_npts2fit(self):
2068        """
2069        setValue Npts for fitting
2070        """
2071        self.Npts_fit.SetValue(str(self.get_npts2fit()))
2072       
2073    def get_chi2(self):
2074        """
2075        return the current chi2
2076        """
2077        return self.tcChi.GetValue()
2078       
2079    def get_param_list(self):
2080        """
2081        :return self.param_toFit: list containing  references to TextCtrl
2082            checked.Theses TextCtrl will allow reference to parameters to fit.
2083       
2084        :raise: if return an empty list of parameter fit will nnote work
2085            properly so raise ValueError,"missing parameter to fit"
2086        """
2087        if self.param_toFit !=[]:
2088            return self.param_toFit
2089        else:
2090            msg = "missing parameters to fit" 
2091            wx.MessageBox(msg, 'warning')
2092            return False
2093            #raise ValueError,"missing parameters to fit"   
2094     
2095    def onsetValues(self, chisqr, p_name, out, cov):
2096        """
2097        Build the panel from the fit result
2098       
2099        :param chisqr: Value of the goodness of fit metric
2100        :param p_name: the name of parameters
2101        :param out: list of parameter with the best value found during fitting
2102        :param cov: Covariance matrix
2103   
2104        """
2105        # make sure stop button to fit button all the time
2106        self._on_fit_complete()
2107        if out == None or not numpy.isfinite(chisqr):
2108            raise ValueError,"Fit error occured..." 
2109       
2110        is_modified = False
2111        has_error = False 
2112        dispersity = ''   
2113       
2114        #Hide textctrl boxes of errors.
2115        self._clear_Err_on_Fit()   
2116       
2117        #Check if chi2 is finite
2118        if chisqr != None and numpy.isfinite(chisqr):
2119        #format chi2 
2120            if self.engine_type == "park": 
2121                npt_fit = float(self.get_npts2fit()) 
2122                # Taking care this in parkemgine
2123                #if npt_fit > 0:
2124                #    chisqr =chisqr/npt_fit   
2125            chi2 = format_number(chisqr, True)   
2126            self.tcChi.SetValue(chi2)   
2127            self.tcChi.Refresh()   
2128        else:
2129            self.tcChi.SetValue("-")
2130       
2131        #Hide error title
2132        if self.text2_3.IsShown() and not self.is_mac:
2133            self.text2_3.Hide()
2134     
2135        try:
2136            if self.enable_disp.GetValue():
2137                if hasattr(self,"text_disp_1" ):
2138                    if self.text_disp_1 != None and not self.is_mac:
2139                        self.text_disp_1.Hide()
2140        except:
2141            dispersity = None
2142            pass
2143     
2144           
2145        i = 0
2146        #Set the panel when fit result are list
2147        for item in self.param_toFit:     
2148            if len(item)>5 and item != None:     
2149                ## reset error value to initial state
2150                if not self.is_mac:
2151                    item[3].Hide()
2152                    item[4].Hide()
2153                for ind in range(len(out)):
2154                    if item[1] == p_name[ind]:
2155                        break       
2156                if len(out)<=len(self.param_toFit) and out[ind] !=None:   
2157                    val_out = format_number(out[ind], True)                 
2158                    item[2].SetValue(val_out)
2159
2160                if(cov !=None and len(cov) == len(out)):
2161                    try:
2162                        if dispersity !=None:
2163                            if self.enable_disp.GetValue():
2164                                if hasattr(self,"text_disp_1" ):
2165                                    if self.text_disp_1!=None:
2166                                        if not self.text_disp_1.IsShown()\
2167                                            and not self.is_mac:
2168                                            self.text_disp_1.Show(True)
2169                    except:
2170                        pass   
2171                   
2172                    if cov[ind]!=None :
2173                        if numpy.isfinite(float(cov[ind])):
2174                            val_err = format_number(cov[ind], True)
2175                            if not self.is_mac:
2176                                item[3].Show(True)
2177                                item[4].Show(True)
2178                            item[4].SetValue(val_err)
2179                            has_error = True
2180                i += 1         
2181        #Show error title when any errors displayed
2182        if has_error: 
2183            if not self.text2_3.IsShown():
2184                self.text2_3.Show(True)   
2185        ## save current state 
2186        self.save_current_state()         
2187       
2188        #self._lay_out()
2189        if not self.is_mac:
2190            self.Layout() 
2191            self.Refresh() 
2192        self._mac_sleep(0.1) 
2193        #plot model ( when drawing, do not update chisqr value again)
2194        self._draw_model(update_chisqr=False, source = 'fit')   
2195        #PostStatusEvent     
2196        #msg = "Fit completed!dddd "
2197        #wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
2198   
2199    def onWeighting(self, event):
2200        """
2201        On Weighting radio button event, sets the weightbt_string
2202        """
2203        self.weightbt_string = event.GetEventObject().GetLabelText()
2204        self._set_weight()
2205       
2206    def _set_weight(self, is_2D=None):
2207        """
2208        Set weight in fit problem
2209        """
2210        # compute weight for the current data
2211        from .utils import get_weight
2212        flag_weight = self.get_weight_flag()
2213        if is_2D == None:
2214            is_2D = self._is_2D()
2215        weight = get_weight(data=self.data, 
2216                            is2d=is_2D, 
2217                            flag=flag_weight)
2218        self._manager.set_fit_weight(uid=self.uid, 
2219                                     flag=flag_weight, 
2220                                     is2d=is_2D, 
2221                                     fid=None)
2222   
2223    def onPinholeSmear(self, event):
2224        """
2225        Create a custom pinhole smear object that will change the way residuals
2226        are compute when fitting
2227       
2228        :Note: accuracy is given by strings'High','Med', 'Low' FOR 2d,
2229                     None for 1D
2230                     
2231        """
2232        if self.model == None:
2233            self.disable_smearer.SetValue(True)
2234            if event == None:
2235                return
2236            msg="Please select a Model first..."
2237            wx.MessageBox(msg, 'Info')
2238            wx.PostEvent(self._manager.parent, StatusEvent(status=\
2239                            "Smear: %s"%msg))
2240            return
2241
2242        # Need update param values
2243        self._update_paramv_on_fit()     
2244
2245        # msg default
2246        msg = None
2247        if event != None:
2248            tcrtl= event.GetEventObject()
2249            # event case of radio button
2250            if tcrtl.GetValue()== True:
2251                self.dx_min = 0.0
2252                self.dx_max = 0.0
2253                is_new_pinhole = True
2254            else:
2255                is_new_pinhole = self._is_changed_pinhole()
2256        else:
2257            is_new_pinhole = True 
2258        # if any value is changed
2259        if is_new_pinhole:
2260            msg = self._set_pinhole_smear()
2261        # hide all silt sizer   
2262        self._hide_all_smear_info()
2263       
2264        ##Calculate chi2
2265        temp_smearer = self.current_smearer
2266        #self.compute_chisqr(smearer= temp_smearer)
2267           
2268        # show relevant slit sizers
2269        self._show_smear_sizer()
2270
2271        self.sizer_set_smearer.Layout()
2272        self.Layout()   
2273       
2274        if event != None: 
2275            event.Skip()                 
2276        #self._undo.Enable(True)
2277        self.save_current_state()
2278        event = PageInfoEvent(page = self)
2279        wx.PostEvent(self.parent, event)
2280       
2281    def _is_changed_pinhole(self): 
2282        """
2283        check if any of pinhole smear is changed
2284       
2285        :return: True or False
2286       
2287        """
2288        # get the values
2289        pin_min = self.smear_pinhole_min.GetValue()
2290        pin_max = self.smear_pinhole_max.GetValue()
2291                   
2292        # Check changes in slit width   
2293        try:
2294            dx_min = float(pin_min)
2295        except:
2296            return True
2297        if self.dx_min != dx_min:
2298            return True
2299       
2300        # Check changes in slit heigth
2301        try:
2302            dx_max = float(pin_max)
2303        except:
2304            return True
2305        if self.dx_max != dx_max:
2306            return True
2307        return False
2308   
2309    def _set_pinhole_smear(self):
2310        """
2311        Set custom pinhole smear
2312       
2313        :return: msg
2314       
2315        """
2316        # copy data
2317        data = copy.deepcopy(self.data)
2318        if self._is_2D():
2319            self.smear_type = 'Pinhole2d'
2320            len_data = len(data.data)
2321            data.dqx_data = numpy.zeros(len_data)
2322            data.dqy_data = numpy.zeros(len_data)
2323        else:
2324            self.smear_type = 'Pinhole'
2325            len_data = len(data.x)
2326            data.dx = numpy.zeros(len_data)
2327            data.dxl = None
2328            data.dxw = None
2329        msg = None
2330   
2331        get_pin_min = self.smear_pinhole_min
2332        get_pin_max = self.smear_pinhole_max       
2333
2334        if not check_float(get_pin_min):
2335            get_pin_min.SetBackgroundColour("pink")
2336            msg= "Model Error:wrong value entered!!!"
2337        elif not check_float(get_pin_max ):
2338            get_pin_max.SetBackgroundColour("pink")
2339            msg= "Model Error:wrong value entered!!!"
2340        else:
2341            if len_data < 2: len_data = 2     
2342            self.dx_min = float(get_pin_min.GetValue())
2343            self.dx_max = float(get_pin_max.GetValue())
2344            if self.dx_min < 0:
2345                get_pin_min.SetBackgroundColour("pink")
2346                msg= "Model Error:This value can not be negative!!!"
2347            elif self.dx_max <0:
2348                get_pin_max.SetBackgroundColour("pink")
2349                msg= "Model Error:This value can not be negative!!!"
2350            elif self.dx_min != None and self.dx_max != None:   
2351                if self._is_2D():
2352                    data.dqx_data[data.dqx_data==0] = self.dx_min
2353                    data.dqy_data[data.dqy_data==0] = self.dx_max         
2354                elif self.dx_min == self.dx_max:
2355                    data.dx[data.dx==0] = self.dx_min
2356                else:
2357                    step = (self.dx_max - self.dx_min)/(len_data-1)   
2358                    data.dx = numpy.arange(self.dx_min,self.dx_max+step/1.1,
2359                                           step)           
2360            elif self.dx_min != None: 
2361                if self._is_2D(): data.dqx_data[data.dqx_data==0] = self.dx_min
2362                else: data.dx[data.dx==0] = self.dx_min
2363            elif self.dx_max != None:
2364                if self._is_2D(): data.dqy_data[data.dqy_data==0] = self.dx_max
2365                else: data.dx[data.dx==0] = self.dx_max         
2366            self.current_smearer = smear_selection(data, self.model)
2367            # 2D need to set accuracy
2368            if self._is_2D(): 
2369                self.current_smearer.set_accuracy(accuracy = \
2370                                                  self.smear2d_accuracy)
2371
2372        if msg != None:
2373            wx.PostEvent(self._manager.parent, StatusEvent(status = msg ))
2374        else:
2375            get_pin_min.SetBackgroundColour("white")
2376            get_pin_max.SetBackgroundColour("white")
2377        ## set smearing value whether or not the data contain the smearing info
2378       
2379        self._manager.set_smearer(smearer=self.current_smearer,
2380                                fid=self.data.id,
2381                                 qmin=float(self.qmin_x),
2382                                 qmax= float(self.qmax_x),
2383                                 enable_smearer=not self.disable_smearer.GetValue(),
2384                                 uid=self.uid)
2385        return msg
2386       
2387    def update_pinhole_smear(self):
2388        """
2389        called by kill_focus on pinhole TextCntrl
2390        to update the changes
2391       
2392        :return: False when wrong value was entered
2393       
2394        """
2395        # msg default
2396        msg = None
2397        # check if any value is changed
2398        if self._is_changed_pinhole():
2399            msg = self._set_pinhole_smear()           
2400        #self._undo.Enable(True)
2401        wx.CallAfter(self.save_current_state)
2402
2403        if msg != None:
2404            return False   
2405        else:
2406            return True
2407                     
2408    def onSlitSmear(self, event):
2409        """
2410        Create a custom slit smear object that will change the way residuals
2411        are compute when fitting
2412        """
2413        if self.model == None:
2414            self.disable_smearer.SetValue(True)
2415            if event == None:
2416                return
2417            msg="Please select a Model first..."
2418            wx.MessageBox(msg, 'Info')
2419            wx.PostEvent(self._manager.parent, StatusEvent(status=\
2420                            "Smear: %s"%msg))
2421            return
2422
2423        # Need update param values
2424        self._update_paramv_on_fit()
2425
2426        # msg default
2427        msg = None
2428        # for event given
2429        if event != None:
2430            tcrtl= event.GetEventObject()
2431            # event case of radio button
2432            if tcrtl.GetValue():
2433                self.dxl = 0.0
2434                self.dxw = 0.0
2435                is_new_slit = True
2436            else:
2437                is_new_slit = self._is_changed_slit()
2438        else:
2439            is_new_slit = True 
2440       
2441        # if any value is changed
2442        if is_new_slit:
2443            msg = self._set_slit_smear()
2444           
2445        # hide all silt sizer
2446        self._hide_all_smear_info()       
2447        ##Calculate chi2
2448        #self.compute_chisqr(smearer= self.current_smearer) 
2449        # show relevant slit sizers       
2450        self._show_smear_sizer()
2451        self.sizer_set_smearer.Layout()
2452        self.Layout()
2453
2454        if event != None:
2455            event.Skip()     
2456        #self._undo.Enable(True)
2457        self.save_current_state()
2458        event = PageInfoEvent(page = self)
2459        wx.PostEvent(self.parent, event)
2460        if msg != None:
2461            wx.PostEvent(self._manager.parent, StatusEvent(status = msg))
2462
2463    def _is_changed_slit(self): 
2464        """
2465        check if any of slit lengths is changed
2466       
2467        :return: True or False
2468       
2469        """
2470        # get the values
2471        width = self.smear_slit_width.GetValue()
2472        height = self.smear_slit_height.GetValue()
2473       
2474        # check and change the box bg color if it was pink
2475        #    but it should be white now
2476        # because this is the case that _set_slit_smear() will not handle
2477        if height.lstrip().rstrip()=="":
2478            self.smear_slit_height.SetBackgroundColour(wx.WHITE)
2479        if width.lstrip().rstrip()=="":
2480            self.smear_slit_width.SetBackgroundColour(wx.WHITE)
2481           
2482        # Check changes in slit width   
2483        if width == "": 
2484            dxw = 0.0
2485        else:
2486            try:
2487                dxw = float(width)
2488            except:
2489                return True
2490        if self.dxw != dxw:
2491            return True
2492       
2493        # Check changes in slit heigth
2494        if height == "": 
2495            dxl = 0.0
2496        else:
2497            try:
2498                dxl = float(height)
2499            except:
2500                return True
2501        if self.dxl != dxl:
2502            return True
2503
2504        return False
2505   
2506    def _set_slit_smear(self):
2507        """
2508        Set custom slit smear
2509       
2510        :return: message to inform the user about the validity
2511            of the values entered for slit smear
2512        """
2513        if self.data.__class__.__name__ ==  "Data2D" or \
2514                        self.enable2D:
2515            return
2516        temp_smearer = None
2517        # make sure once more if it is smearer
2518        data = copy.deepcopy(self.data)
2519        data_len = len(data.x)
2520        data.dx = None
2521        data.dxl = None
2522        data.dxw = None
2523        msg = None
2524   
2525        try:
2526            self.dxl = float(self.smear_slit_height.GetValue())
2527            data.dxl = self.dxl* numpy.ones(data_len)
2528            self.smear_slit_height.SetBackgroundColour(wx.WHITE)
2529        except: 
2530            data.dxl = numpy.zeros(data_len)
2531            if self.smear_slit_height.GetValue().lstrip().rstrip()!="":
2532                self.smear_slit_height.SetBackgroundColour("pink")
2533                msg = "Wrong value entered... " 
2534            else:
2535                self.smear_slit_height.SetBackgroundColour(wx.WHITE)
2536        try:
2537            self.dxw = float(self.smear_slit_width.GetValue())
2538            self.smear_slit_width.SetBackgroundColour(wx.WHITE)
2539            data.dxw = self.dxw* numpy.ones(data_len)
2540        except: 
2541            data.dxw = numpy.zeros(data_len)
2542            if self.smear_slit_width.GetValue().lstrip().rstrip()!="":
2543                self.smear_slit_width.SetBackgroundColour("pink")
2544                msg = "Wrong Fit value entered... "
2545            else:
2546                self.smear_slit_width.SetBackgroundColour(wx.WHITE)
2547             
2548        self.current_smearer = smear_selection(data, self.model)
2549        #temp_smearer = self.current_smearer
2550        ## set smearing value whether or not the data contain the smearing info
2551        self._manager.set_smearer(smearer=self.current_smearer, 
2552                                 fid=self.data.id,
2553                                 qmin=float(self.qmin_x), 
2554                                 qmax= float(self.qmax_x),
2555                        enable_smearer=not self.disable_smearer.GetValue(),
2556                                 uid=self.uid) 
2557        return msg
2558   
2559    def update_slit_smear(self):
2560        """
2561        called by kill_focus on pinhole TextCntrl
2562        to update the changes
2563       
2564        :return: False when wrong value was entered
2565       
2566        """             
2567        # msg default
2568        msg = None
2569        # check if any value is changed
2570        if self._is_changed_slit():
2571            msg = self._set_slit_smear()           
2572        #self._undo.Enable(True)
2573        self.save_current_state()
2574
2575        if msg != None:
2576            return False   
2577        else:
2578            return True
2579                           
2580    def onSmear(self, event):
2581        """
2582        Create a smear object that will change the way residuals
2583        are compute when fitting
2584        """
2585        if event != None: 
2586            event.Skip()   
2587        if self.data is None:
2588            return
2589
2590        if self.model == None:
2591            self.disable_smearer.SetValue(True)
2592            if event == None:
2593                return
2594            msg="Please select a Model first..."
2595            wx.MessageBox(msg, 'Info')
2596            wx.PostEvent(self._manager.parent, StatusEvent(status=\
2597                            "Smear: %s"%msg))
2598            return
2599        # Need update param values
2600        self._update_paramv_on_fit()
2601        temp_smearer = self.on_smear_helper()           
2602       
2603        self.sizer_set_smearer.Layout()
2604        self.Layout()
2605        self._set_weight()
2606       
2607        ## set smearing value whether or not the data contain the smearing info
2608        wx.CallAfter(self._manager.set_smearer, uid=self.uid, smearer=temp_smearer,
2609                                  fid=self.data.id,
2610                        qmin=float(self.qmin_x),
2611                        qmax=float(self.qmax_x),
2612                        enable_smearer=not self.disable_smearer.GetValue(),
2613                         draw=True) 
2614       
2615        self.state.enable_smearer=  self.enable_smearer.GetValue()
2616        self.state.disable_smearer=self.disable_smearer.GetValue()
2617        self.state.pinhole_smearer = self.pinhole_smearer.GetValue()
2618        self.state.slit_smearer = self.slit_smearer.GetValue()
2619       
2620    def on_smear_helper(self, update=False):
2621        """
2622        Help for onSmear
2623       
2624        :param update: force or not to update
2625        """
2626        self._get_smear_info()
2627        temp_smearer = None
2628        #renew smear sizer
2629        if self.smear_type != None:
2630            self.smear_description_smear_type.SetValue(str(self.smear_type))
2631            self.smear_data_left.SetValue(str(self.dq_l))   
2632            self.smear_data_right.SetValue(str(self.dq_r)) 
2633
2634        self._hide_all_smear_info()
2635        data = copy.deepcopy(self.data)
2636       
2637        # make sure once more if it is smearer
2638        temp_smearer = smear_selection(data, self.model)
2639        if self.current_smearer != temp_smearer or update:
2640            self.current_smearer = temp_smearer
2641        if self.enable_smearer.GetValue():
2642            if hasattr(self.data,"dxl"):
2643                msg= ": Resolution smearing parameters"
2644            if hasattr(self.data,"dxw"):
2645                msg= ": Slit smearing parameters"
2646            if self.current_smearer ==None:
2647                wx.PostEvent(self._manager.parent, StatusEvent(status=\
2648                            "Data contains no smearing information"))
2649            else:
2650                wx.PostEvent(self._manager.parent, StatusEvent(status=\
2651                            "Data contains smearing information"))
2652
2653            #self.smear_description_dqdata.Show(True)
2654            self.smear_data_left.Show(True)
2655            self.smear_data_right.Show(True)
2656            temp_smearer= self.current_smearer
2657        elif self.disable_smearer.GetValue():
2658            self.smear_description_none.Show(True)
2659        elif self.pinhole_smearer.GetValue():
2660            self.onPinholeSmear(None)
2661        elif self.slit_smearer.GetValue():
2662            self.onSlitSmear(None)
2663        self._show_smear_sizer()
2664       
2665        return temp_smearer
2666   
2667    def on_complete_chisqr(self, event): 
2668        """
2669        Display result chisqr on the panel
2670        :event: activated by fitting/ complete after draw
2671        """
2672        try:
2673            if event == None:
2674                output= "-"
2675            elif not numpy.isfinite(event.output):
2676                output= "-"
2677            else:
2678                output = event.output
2679            self.tcChi.SetValue(str(format_number(output, True)))
2680            self.state.tcChi = self.tcChi.GetValue()
2681        except:
2682            pass 
2683           
2684   
2685    def get_all_checked_params(self):
2686        """
2687        Found all parameters current check and add them to list of parameters to
2688        fit
2689        """
2690        self.param_toFit = []
2691        for item in self.parameters:
2692            if item[0].GetValue() and item not in self.param_toFit:
2693                self.param_toFit.append(item)
2694        for item in self.fittable_param:
2695            if item[0].GetValue() and item not in self.param_toFit:
2696                self.param_toFit.append(item)
2697        self.save_current_state_fit() 
2698       
2699     
2700        event = PageInfoEvent(page = self)
2701        wx.PostEvent(self.parent, event) 
2702        param2fit = []
2703        for item in self.param_toFit:
2704            if item[0]:
2705                param2fit.append(item[1])     
2706        self.parent._manager.set_param2fit(self.uid, param2fit)
2707               
2708    def select_all_param(self, event): 
2709        """
2710        set to true or false all checkBox given the main checkbox value cb1
2711        """           
2712        self.param_toFit = []
2713        if  self.parameters !=[]:
2714            if  self.cb1.GetValue():
2715                for item in self.parameters:
2716                    ## for data2D select all to fit
2717                    if self.data.__class__.__name__==  "Data2D" or \
2718                            self.enable2D:
2719                        item[0].SetValue(True)
2720                        self.param_toFit.append(item)
2721                    else:
2722                        ## for 1D all parameters except orientation
2723                        if not item in self.orientation_params:
2724                            item[0].SetValue(True)
2725                            self.param_toFit.append(item )
2726                #if len(self.fittable_param)>0:
2727                for item in self.fittable_param:
2728                    if self.data.__class__.__name__== "Data2D" or \
2729                            self.enable2D:
2730                        item[0].SetValue(True)
2731                        self.param_toFit.append(item)
2732                        try:
2733                            if len(self.values[item[1]]) > 0:
2734                                item[0].SetValue(False)
2735                        except:
2736                            pass
2737
2738                    else:
2739                        ## for 1D all parameters except orientation
2740                        if not item in self.orientation_params_disp:
2741                            item[0].SetValue(True)
2742                            self.param_toFit.append(item )
2743                            try:
2744                                if len(self.values[item[1]]) > 0:
2745                                    item[0].SetValue(False)
2746                            except:
2747                                pass
2748
2749            else:
2750                for item in self.parameters:
2751                    item[0].SetValue(False)
2752                for item in self.fittable_param:
2753                    item[0].SetValue(False)
2754                self.param_toFit=[]
2755           
2756        self.save_current_state_fit() 
2757       
2758        if event !=None:
2759            #self._undo.Enable(True)
2760            ## post state to fit panel
2761            event = PageInfoEvent(page = self)
2762            wx.PostEvent(self.parent, event) 
2763        param2fit = []
2764        for item in self.param_toFit:
2765            if item[0]:
2766                param2fit.append(item[1])     
2767        self.parent._manager.set_param2fit(self.uid, param2fit)
2768
2769       
2770    def select_param(self,event):
2771        """
2772        Select TextCtrl  checked for fitting purpose and stores them
2773        in  self.param_toFit=[] list
2774        """
2775        self.param_toFit = []
2776        for item in self.parameters:
2777            #Skip t ifhe angle parameters if 1D data
2778            if self.data.__class__.__name__ != "Data2D" and\
2779                        not self.enable2D:
2780                if item in self.orientation_params:
2781                    continue
2782            #Select parameters to fit for list of primary parameters
2783            if item[0].GetValue():
2784                if not (item in self.param_toFit):
2785                    self.param_toFit.append(item ) 
2786            else:
2787                #remove parameters from the fitting list
2788                if item in self.param_toFit:
2789                    self.param_toFit.remove(item)
2790
2791        #Select parameters to fit for list of fittable parameters
2792        #        with dispersion         
2793        for item in self.fittable_param:
2794            #Skip t ifhe angle parameters if 1D data
2795            if self.data.__class__.__name__ !=  "Data2D" and\
2796                        not self.enable2D:
2797                if item in self.orientation_params:
2798                    continue
2799            if item[0].GetValue():
2800                if not (item in self.param_toFit):
2801                    self.param_toFit.append(item) 
2802            else:
2803                #remove parameters from the fitting list
2804                if item in self.param_toFit:
2805                    self.param_toFit.remove(item)
2806
2807        #Calculate num. of angle parameters
2808        if self.data.__class__.__name__ ==  "Data2D" or \
2809                       self.enable2D:
2810            len_orient_para = 0
2811        else:
2812            len_orient_para = len(self.orientation_params)  #assume even len
2813        #Total num. of angle parameters
2814        if len(self.fittable_param) > 0:
2815            len_orient_para *= 2
2816        #Set the value of checkbox that selected every checkbox or not       
2817        if len(self.parameters)+len(self.fittable_param)-len_orient_para ==\
2818                len(self.param_toFit):
2819            self.cb1.SetValue(True)
2820        else:
2821            self.cb1.SetValue(False)
2822       
2823       
2824       
2825        self.save_current_state_fit()
2826        if event !=None:
2827            #self._undo.Enable(True)
2828            ## post state to fit panel
2829            event = PageInfoEvent(page = self)
2830            wx.PostEvent(self.parent, event) 
2831       
2832        param2fit = []
2833        for item in self.param_toFit:
2834            if item[0]:
2835                param2fit.append(item[1])     
2836        self.parent._manager.set_param2fit(self.uid, param2fit)
2837       
2838    def set_model_param_sizer(self, model):
2839        """
2840        Build the panel from the model content
2841       
2842        :param model: the model selected in combo box for fitting purpose
2843       
2844        """
2845        self.sizer3.Clear(True)
2846        self.parameters = []
2847        self.str_parameters = []
2848        self.param_toFit=[]
2849        self.fittable_param=[]
2850        self.fixed_param=[]
2851        self.orientation_params=[]
2852        self.orientation_params_disp=[]
2853       
2854        if model ==None:
2855            self.sizer3.Layout()
2856            self.SetupScrolling()
2857            return
2858        ## the panel is drawn using the current value of the fit engine
2859        if self.engine_type==None and self._manager !=None:
2860            self.engine_type= self._manager._return_engine_type()
2861
2862        box_description= wx.StaticBox(self, -1,str("Model Parameters"))
2863        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
2864        sizer = wx.GridBagSizer(5,5)
2865        ## save the current model
2866        self.model = model
2867           
2868        keys = self.model.getParamList()
2869        #list of dispersion paramaters
2870        self.disp_list=self.model.getDispParamList()
2871
2872        def custom_compare(a,b):
2873            """
2874            Custom compare to order, first by alphabets then second by number.
2875            """ 
2876            # number at the last digit
2877            a_last = a[len(a)-1]
2878            b_last = b[len(b)-1]
2879            # default
2880            num_a = None
2881            num_b = None
2882            # split the names
2883            a2 = a.lower().split('_')
2884            b2 = b.lower().split('_')
2885            # check length of a2, b2
2886            len_a2 = len(a2)
2887            len_b2 = len(b2)
2888            # check if it contains a int number(<10)
2889            try: 
2890                num_a = int(a_last)
2891            except: pass
2892            try:
2893                num_b = int(b_last)
2894            except: pass
2895            # Put 'scale' near the top; happens
2896            # when numbered param name exists
2897            if a == 'scale':
2898                return -1
2899            # both have a number   
2900            if num_a != None and num_b != None:
2901                if num_a > num_b: return -1
2902                # same number
2903                elif num_a == num_b: 
2904                    # different last names
2905                    if a2[len_a2-1] != b2[len_b2-1] and num_a != 0:
2906                        return -cmp(a2[len_a2-1], b2[len_b2-1])
2907                    else: 
2908                        return cmp(a, b) 
2909                else: return 1
2910            # one of them has a number
2911            elif num_a != None: return 1
2912            elif num_b != None: return -1
2913            # no numbers
2914            else: return cmp(a.lower(), b.lower())
2915
2916                       
2917        keys.sort(custom_compare)
2918   
2919        iy = 0
2920        ix = 0
2921        select_text = "Select All"
2922        self.cb1 = wx.CheckBox(self, -1,str(select_text), (10, 10))
2923        wx.EVT_CHECKBOX(self, self.cb1.GetId(), self.select_all_param)
2924        self.cb1.SetToolTipString("To check/uncheck all the boxes below.")
2925        self.cb1.SetValue(True)
2926       
2927        sizer.Add(self.cb1,(iy, ix),(1,1),\
2928                             wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 5)
2929        ix +=1
2930        self.text2_2 = wx.StaticText(self, -1, 'Value')
2931        sizer.Add(self.text2_2,(iy, ix),(1,1),\
2932                            wx.EXPAND|wx.ADJUST_MINSIZE, 0)
2933        ix +=2 
2934        self.text2_3 = wx.StaticText(self, -1, 'Error')
2935        sizer.Add(self.text2_3,(iy, ix),(1,1),\
2936                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
2937        if not self.is_mac:
2938            self.text2_3.Hide()
2939        ix +=1 
2940        self.text2_min = wx.StaticText(self, -1, 'Min')
2941        sizer.Add(self.text2_min,(iy, ix),(1,1),\
2942                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
2943        #self.text2_min.Hide()
2944        ix +=1 
2945        self.text2_max = wx.StaticText(self, -1, 'Max')
2946        sizer.Add(self.text2_max,(iy, ix),(1,1),\
2947                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
2948        #self.text2_max.Hide()
2949        ix += 1
2950        self.text2_4 = wx.StaticText(self, -1, '[Units]')
2951        sizer.Add(self.text2_4,(iy, ix),(1,1),\
2952                            wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
2953        self.text2_4.Hide()
2954       
2955        CHECK_STATE = self.cb1.GetValue()
2956        for item in keys:
2957            if not item in self.disp_list and not item in \
2958                    self.model.orientation_params:
2959               
2960                ##prepare a spot to store errors
2961                if not self.model.details.has_key(item):
2962                    self.model.details [item] = ["",None,None] 
2963         
2964                iy += 1
2965                ix = 0
2966                if (self.model.__class__ in \
2967                    self.model_list_box["Multi-Functions"] or \
2968                    self.temp_multi_functional)\
2969                    and (item in self.model.non_fittable):
2970                    non_fittable_name = wx.StaticText(self, -1, item )
2971                    sizer.Add(non_fittable_name,(iy, ix),(1,1),\
2972                            wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 21)
2973                    ## add parameter value
2974                    ix += 1
2975                    value= self.model.getParam(item)
2976                    if len(self.model.fun_list) > 0:
2977                        num = item.split('_')[1][5:7]
2978                        fun_box = wx.ComboBox(self, -1,size=(100,-1),
2979                                    style=wx.CB_READONLY, name = '%s'% item)
2980                        self._set_fun_box_list(fun_box)
2981                        fun_box.SetSelection(0)
2982                        #self.fun_box.SetToolTipString("A function
2983                        #    describing the interface")
2984                        wx.EVT_COMBOBOX(fun_box,-1, self._on_fun_box)
2985                    else:
2986                        fun_box = self.ModelTextCtrl(self, -1, 
2987                                                     size=(_BOX_WIDTH,20),
2988                                style=wx.TE_PROCESS_ENTER, name ='%s'% item)
2989                        fun_box.SetToolTipString(\
2990                                "Hit 'Enter' after typing to update the plot.")
2991                        fun_box.SetValue(format_number(value, True))
2992                    sizer.Add(fun_box, (iy,ix),(1,1), wx.EXPAND)
2993                    self.str_parameters.append([None,item, fun_box,
2994                                                None, None, None, 
2995                                                None, None])
2996
2997                   
2998                else:
2999                    ## add parameters name with checkbox for selecting to fit
3000                    cb = wx.CheckBox(self, -1, item) 
3001                    cb.SetValue(CHECK_STATE)           
3002                    cb.SetToolTipString(" Check mark to fit.")
3003                    #cb.SetValue(True)
3004                    wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
3005                   
3006                    sizer.Add( cb,( iy, ix),(1,1),
3007                                 wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 5)
3008
3009                    ## add parameter value
3010                    ix += 1
3011                    value= self.model.getParam(item)
3012                    ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH,20),
3013                                        style=wx.TE_PROCESS_ENTER)
3014                    ctl1.SetToolTipString(\
3015                                "Hit 'Enter' after typing to update the plot.")
3016                    ctl1.SetValue(format_number(value, True))
3017                    sizer.Add(ctl1, (iy,ix),(1,1), wx.EXPAND)
3018                    ## text to show error sign
3019                    ix += 1
3020                    text2=wx.StaticText(self, -1, '+/-')
3021                    sizer.Add(text2,(iy, ix),(1,1),\
3022                                    wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
3023                    if not self.is_mac:
3024                        text2.Hide() 
3025                    ix += 1
3026                    ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.2,20), style=0)
3027                    sizer.Add(ctl2, (iy,ix),(1,1), 
3028                              wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3029                    if not self.is_mac:
3030                        ctl2.Hide()
3031                   
3032                    ix += 1
3033                    ctl3 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.9,20),
3034                                               style=wx.TE_PROCESS_ENTER,
3035                                text_enter_callback = self._onparamRangeEnter)
3036         
3037                    sizer.Add(ctl3, (iy,ix),(1,1), 
3038                              wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3039           
3040                    ix += 1
3041                    ctl4 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.9,20),
3042                                               style=wx.TE_PROCESS_ENTER,
3043                                text_enter_callback = self._onparamRangeEnter)
3044                    sizer.Add(ctl4, (iy,ix),(1,1), 
3045                              wx.EXPAND|wx.FIXED_MINSIZE, 0)
3046   
3047                    ix +=1
3048                    # Units
3049                    if self.model.details.has_key(item):
3050                        units = wx.StaticText(self, -1, 
3051                            self.model.details[item][0], style=wx.ALIGN_LEFT)
3052                    else:
3053                        units = wx.StaticText(self, -1, "", style=wx.ALIGN_LEFT)
3054                    sizer.Add(units, (iy,ix),(1,1), 
3055                               wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3056                       
3057                    ##[cb state, name, value, "+/-", error of fit, min, max , units]
3058                    self.parameters.append([cb,item, ctl1,
3059                                            text2,ctl2, ctl3, ctl4,units])
3060                                 
3061        iy+=1
3062        sizer.Add((10,10),(iy,ix),(1,1), 
3063                  wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15)
3064       
3065        # type can be either Guassian or Array
3066        if len(self.model.dispersion.values())>0:
3067            type= self.model.dispersion.values()[0]["type"]
3068        else:
3069            type = "Gaussian"
3070           
3071        iy += 1
3072        ix = 0
3073        #Add tile for orientational angle
3074        for item in keys:
3075            if item in self.model.orientation_params:       
3076                orient_angle = wx.StaticText(self, -1, '[For 2D only]:')
3077                sizer.Add(orient_angle,(iy, ix),(1,1), 
3078                          wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 15) 
3079                if not self.data.__class__.__name__ == "Data2D" and \
3080                        not self.enable2D:
3081                    orient_angle.Hide()
3082                else:
3083                    orient_angle.Show(True)
3084                break
3085     
3086        #For Gaussian only
3087        if type.lower() != "array":
3088            for item in self.model.orientation_params:
3089                if not item in self.disp_list:
3090                    ##prepare a spot to store min max
3091                    if not self.model.details.has_key(item):
3092                        self.model.details [item] = ["",None,None] 
3093                         
3094                    iy += 1
3095                    ix = 0
3096                    ## add parameters name with checkbox for selecting to fit
3097                    cb = wx.CheckBox(self, -1, item )
3098                    cb.SetValue(CHECK_STATE)
3099                    cb.SetToolTipString("Check mark to fit")
3100                    wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
3101                    if self.data.__class__.__name__ == "Data2D" or \
3102                            self.enable2D:
3103                        cb.Show(True)
3104                    else:
3105                        cb.Hide()
3106                    sizer.Add( cb,( iy, ix),(1,1),
3107                                 wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE, 5)
3108   
3109                    ## add parameter value
3110                    ix += 1
3111                    value= self.model.getParam(item)
3112                    ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH,20),
3113                                        style=wx.TE_PROCESS_ENTER)
3114                    ctl1.SetToolTipString(\
3115                                "Hit 'Enter' after typing to update the plot.")
3116                    ctl1.SetValue(format_number(value, True))
3117                    if self.data.__class__.__name__ == "Data2D" or \
3118                            self.enable2D:
3119                        ctl1.Show(True)
3120                    else:
3121                        ctl1.Hide()
3122                    sizer.Add(ctl1, (iy,ix),(1,1), wx.EXPAND)
3123                    ## text to show error sign
3124                    ix += 1
3125                    text2=wx.StaticText(self, -1, '+/-')
3126                    sizer.Add(text2,(iy, ix),(1,1),\
3127                                    wx.EXPAND|wx.ADJUST_MINSIZE, 0) 
3128
3129                    text2.Hide() 
3130                    ix += 1
3131                    ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.2,20), style=0)
3132                    sizer.Add(ctl2, (iy,ix),(1,1), 
3133                              wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3134
3135                    ctl2.Hide()
3136                   
3137                   
3138                    ix += 1
3139                    ctl3 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.8,20), 
3140                                              style=wx.TE_PROCESS_ENTER,
3141                                text_enter_callback = self._onparamRangeEnter)
3142                   
3143                    sizer.Add(ctl3, (iy,ix),(1,1), 
3144                              wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3145                    ctl3.Hide()
3146                 
3147                    ix += 1
3148                    ctl4 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.8,20),
3149                                               style=wx.TE_PROCESS_ENTER,
3150                            text_enter_callback = self._onparamRangeEnter)
3151                    sizer.Add(ctl4, (iy,ix),(1,1), 
3152                              wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3153                   
3154                    ctl4.Hide()
3155                   
3156                    if self.data.__class__.__name__ == "Data2D" or \
3157                            self.enable2D: 
3158                        if self.is_mac: 
3159                            text2.Show(True)
3160                            ctl2.Show(True)                 
3161                        ctl3.Show(True)
3162                        ctl4.Show(True)
3163                   
3164                    ix +=1
3165                    # Units
3166                    if self.model.details.has_key(item):
3167                        units = wx.StaticText(self, -1, 
3168                                              self.model.details[item][0], 
3169                                              style=wx.ALIGN_LEFT)
3170                    else:
3171                        units = wx.StaticText(self, -1, "", style=wx.ALIGN_LEFT)
3172                    if self.data.__class__.__name__ == "Data2D" or \
3173                            self.enable2D:
3174                        units.Show(True)
3175                   
3176                    else:
3177                        units.Hide()
3178                   
3179                    sizer.Add(units, (iy,ix),(1,1), 
3180                              wx.EXPAND|wx.ADJUST_MINSIZE, 0)
3181                                         
3182                    ##[cb state, name, value, "+/-", error of fit,min,max,units]
3183                    self.parameters.append([cb,item, ctl1,
3184                                            text2,ctl2, ctl3, ctl4,units])
3185                    self.orientation_params.append([cb,item, ctl1,
3186                                            text2,ctl2, ctl3, ctl4,units])
3187             
3188        iy+=1
3189       
3190        #Display units text on panel
3191        for item in keys:   
3192            if self.model.details.has_key(item):
3193                self.text2_4.Show()
3194        #Fill the list of fittable parameters
3195        self.get_all_checked_params()
3196        self.save_current_state_fit()
3197        boxsizer1.Add(sizer)
3198        self.sizer3.Add(boxsizer1,0, wx.EXPAND | wx.ALL, 10)
3199        self.sizer3.Layout()
3200        self.Layout()
3201        #self.Refresh()
3202
3203    def on_right_down(self, event):
3204        """
3205        Get key stroke event
3206        """
3207        if self.data == None:
3208            return
3209        # Figuring out key combo: Cmd for copy, Alt for paste
3210        if event.AltDown() and event.ShiftDown():
3211            self._manager.show_ftol_dialog()
3212            flag = True
3213        elif event.AltDown() or event.ShiftDown():
3214            flag = False
3215        else:
3216            return
3217        # make event free
3218        event.Skip()
3219        # messages depending on the flag
3220        if not flag:
3221            msg = " Could not open ftol dialog;"
3222            msg += " Check if the Scipy fit engine is selected in the menubar."
3223            infor = 'warning'
3224            # inform msg to wx
3225            wx.PostEvent( self.parent.parent, 
3226                          StatusEvent(status= msg, info=infor))
3227
3228       
3229    def _onModel2D(self, event):
3230        """
3231        toggle view of model from 1D to 2D  or 2D from 1D
3232        """
3233        if self.model_view.GetLabelText() == "Show 2D":
3234            self.model_view.SetLabel("Show 1D")
3235            self.enable2D = True
3236             
3237        else:
3238            self.model_view.SetLabel("Show 2D")
3239            self.enable2D = False
3240        self.Show(False)
3241        self.create_default_data()
3242        self._manager.store_data(self.uid, data_list=[self.data])
3243
3244        self.set_model_param_sizer(self.model)
3245        self._set_sizer_dispersion() 
3246        self._set_weight(is_2D= self.enable2D)
3247        self._set_smear_buttons()
3248        self.Show(True)
3249        self.SetupScrolling()
3250        self._draw_model()
3251       
3252        self.state.enable2D =  copy.deepcopy(self.enable2D)
3253   
3254    def _set_smear_buttons(self):
3255        """
3256        Set semarer radio buttons
3257        """
3258        # more disables for 2D
3259        if self.data.__class__.__name__ ==  "Data2D" or \
3260                    self.enable2D:
3261            self.slit_smearer.Disable()
3262            self.pinhole_smearer.Enable(True) 
3263            self.default_mask = copy.deepcopy(self.data.mask)
3264        else:
3265            self.slit_smearer.Enable(True) 
3266            self.pinhole_smearer.Enable(True) 
3267           
3268           
3269class BGTextCtrl(wx.TextCtrl):
3270    """
3271    Text control used to display outputs.
3272    No editing allowed. The background is
3273    grayed out. User can't select text.
3274    """
3275    def __init__(self, *args, **kwds):
3276        wx.TextCtrl.__init__(self, *args, **kwds)
3277        self.SetEditable(False)
3278        self.SetBackgroundColour(self.GetParent().parent.GetBackgroundColour())
3279       
3280        # Bind to mouse event to avoid text highlighting
3281        # The event will be skipped once the call-back
3282        # is called.
3283        self.Bind(wx.EVT_MOUSE_EVENTS, self._click)
3284       
3285    def _click(self, event):
3286        """
3287        Prevent further handling of the mouse event
3288        by not calling Skip().
3289        """ 
3290        pass
3291 
Note: See TracBrowser for help on using the repository browser.