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

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

fixed a bug that crashes sansview on_selete model in win7 (within python env)

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