source: sasview/src/sas/sasgui/perspectives/fitting/basepage.py @ 06a4306

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.2ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since 06a4306 was 06a4306, checked in by GitHub <noreply@…>, 7 years ago

Merge pull request #16 from SasView?/ticket778

Ticket778

  • Property mode set to 100644
File size: 143.0 KB
Line 
1"""
2Base Page for fitting
3"""
4import sys
5import os
6import wx
7import numpy
8import time
9import copy
10import math
11import json
12import logging
13import traceback
14
15from collections import defaultdict
16from wx.lib.scrolledpanel import ScrolledPanel
17
18from sasmodels.weights import MODELS as POLYDISPERSITY_MODELS
19
20from sas.sasgui.guiframe.panel_base import PanelBase
21from sas.sasgui.guiframe.utils import format_number, check_float, IdList, \
22    check_int
23from sas.sasgui.guiframe.events import PanelOnFocusEvent
24from sas.sasgui.guiframe.events import StatusEvent
25from sas.sasgui.guiframe.events import AppendBookmarkEvent
26from sas.sasgui.guiframe.dataFitting import Data2D
27from sas.sasgui.guiframe.dataFitting import Data1D
28from sas.sasgui.guiframe.dataFitting import check_data_validity
29from sas.sasgui.guiframe.gui_style import GUIFRAME_ID
30from sas.sascalc.dataloader.data_info import Detector
31from sas.sascalc.dataloader.data_info import Source
32from sas.sasgui.perspectives.fitting.pagestate import PageState
33from sas.sasgui.guiframe.CategoryInstaller import CategoryInstaller
34from sas.sasgui.guiframe.documentation_window import DocumentationWindow
35
36
37(PageInfoEvent, EVT_PAGE_INFO) = wx.lib.newevent.NewEvent()
38(PreviousStateEvent, EVT_PREVIOUS_STATE) = wx.lib.newevent.NewEvent()
39(NextStateEvent, EVT_NEXT_STATE) = wx.lib.newevent.NewEvent()
40
41_BOX_WIDTH = 76
42_QMIN_DEFAULT = 0.0005
43_QMAX_DEFAULT = 0.5
44_NPTS_DEFAULT = 50
45# Control panel width
46if sys.platform.count("win32") > 0:
47    PANEL_WIDTH = 450
48    FONT_VARIANT = 0
49    ON_MAC = False
50else:
51    PANEL_WIDTH = 500
52    FONT_VARIANT = 1
53    ON_MAC = True
54
55class BasicPage(ScrolledPanel, PanelBase):
56    """
57    This class provide general structure of  fitpanel page
58    """
59    # Internal name for the AUI manager
60    window_name = "Fit Page"
61    # Title to appear on top of the window
62    window_caption = "Fit Page "
63
64    # These two buttons have specific IDs since they seem to be created more
65    # frequently than they need to.  In particular, set_dispers_sizer() is
66    # called by _on_select_model
67    ID_BOOKMARK = wx.NewId()
68    ID_DISPERSER_HELP = wx.NewId()
69    _id_pool = IdList()
70
71    def __init__(self, parent, color='blue', **kwargs):
72        """
73        """
74        ScrolledPanel.__init__(self, parent, **kwargs)
75        PanelBase.__init__(self, parent)
76        self.SetupScrolling()
77        # Set window's font size
78        self.SetWindowVariant(variant=FONT_VARIANT)
79        self.SetBackgroundColour(color)
80
81        self._ids = iter(self._id_pool)
82        # parent of the page
83        self.parent = parent
84        # manager is the fitting plugin
85        # owner of the page (fitting plugin)
86        self.event_owner = None
87        # current model
88        self.model = None
89        self.m_name = None
90        self.index_model = None
91        self.panel = None
92        # data
93        self.data = None
94        # list of available data
95        self.data_list = []
96        self.mask = None
97        self.uid = wx.NewId()
98        self.graph_id = None
99        # Q range for data set
100        self.qmin_data_set = numpy.inf
101        self.qmax_data_set = None
102        self.npts_data_set = 0
103        # Q range
104        self.qmin = None
105        self.qmax = None
106        self.qmax_x = _QMAX_DEFAULT
107        self.qmin_x = _QMIN_DEFAULT
108        self.npts_x = _NPTS_DEFAULT
109        # total number of point: float
110        self.npts = None
111        self.num_points = None
112        # smear default
113        self.current_smearer = None
114        # 2D smear accuracy default
115        self.smear2d_accuracy = 'Low'
116        # slit smear:
117        self.dxl = None
118        self.dxw = None
119        # pinhole smear
120        self.dx_min = None
121        self.dx_max = None
122        # smear attrbs
123        self.enable_smearer = None
124        self.disable_smearer = None
125        self.pinhole_smearer = None
126        self.slit_smearer = None
127        # weight attrbs
128        self.dI_noweight = None
129        self.dI_didata = None
130        self.dI_sqrdata = None
131        self.dI_idata = None
132        # other attrbs
133        self.dq_l = None
134        self.dq_r = None
135        self.tcChi = None
136        self.disp_box = None
137        self.model_disp = None
138        self.Npts_fit = None
139        self.Npts_total = None
140        self.theory_qmin = None
141        self.theory_qmax = None
142        self.theory_qmin_x = None
143        self.theory_qmax_x = None
144        self.btEditMask = None
145        self.btFit = None
146        self.sld_axes = None
147        self.multi_factor = None
148
149        self.disp_cb_dict = {}
150
151        # self.state = PageState(parent=parent)
152        # dictionary containing list of models
153        self.model_list_box = {}
154
155        # Data member to store the dispersion object created
156        self._disp_obj_dict = {}
157        # selected parameters to apply dispersion
158        self.disp_cb_dict = {}
159        # smearer object
160        self.enable2D = False
161        self._has_magnetic = False
162        self.magnetic_on = False
163        self.is_mac = ON_MAC
164        self.formfactorbox = None
165        self.structurebox = None
166        self.categorybox = None
167        # list of model parameters. each item must have same length
168        # each item related to a given parameters
169        # [cb state, name, value, "+/-", error of fit, min, max , units]
170        self.parameters = []
171        # non-fittable parameter whose value is astring
172        self.str_parameters = []
173        # list of parameters to fit , must be like self.parameters
174        self.param_toFit = []
175        # list of looking like parameters but with non fittable parameters info
176        self.fixed_param = []
177        # list of looking like parameters but with  fittable parameters info
178        self.fittable_param = []
179        # list of dispersion parameters
180        self.disp_list = []
181        self.disp_name = ""
182
183        # list of orientation parameters
184        self.orientation_params = []
185        self.orientation_params_disp = []
186#       Self.model should ALWAYS be None here.  It was set to none above in
187#       this long init setting.  no obvious function call in between setting
188#       and this - commenting out on 4/8/2014 by PDB.  Remove once clear
189#       it is pointless.
190#        if self.model is not None:
191#            self.disp_list = self.model.getDispParamList()
192        self.temp_multi_functional = False
193        # enable model 2D draw
194        self.enable2D = False
195        # check that the fit range is correct to plot the model again
196        self.fitrange = True
197        # Create memento to save the current state
198        self.state = PageState(parent=self.parent,
199                               model=self.model, data=self.data)
200        # flag to determine if state has change
201        self.state_change = False
202        # save customized array
203        self.values = {}   # type: Dict[str, List[float, ...]]
204        self.weights = {}   # type: Dict[str, List[float, ...]]
205        # retrieve saved state
206        self.number_saved_state = 0
207        # dictionary of saved state
208        self.saved_states = {}
209        # Create context menu for page
210        self.popUpMenu = wx.Menu()
211
212        wx_id = self._ids.next()
213        self._keep = wx.MenuItem(self.popUpMenu, wx_id, "Add bookmark",
214                                 " Keep the panel status to recall it later")
215        self.popUpMenu.AppendItem(self._keep)
216        self._keep.Enable(False)
217        self._set_bookmark_flag(False)
218        self._set_save_flag(False)
219        wx.EVT_MENU(self, wx_id, self.on_bookmark)
220        self.popUpMenu.AppendSeparator()
221
222        # Default locations
223        self._default_save_location = os.getcwd()
224        # save initial state on context menu
225        # self.onSave(event=None)
226        self.Bind(wx.EVT_CONTEXT_MENU, self.onContextMenu)
227
228        # bind key event
229        self.Bind(wx.EVT_LEFT_DOWN, self.on_left_down)
230
231        # create the basic structure of the panel with empty sizer
232        self.define_page_structure()
233        # drawing Initial dispersion parameters sizer
234        self.set_dispers_sizer()
235
236        # layout
237        self.set_layout()
238
239    def set_index_model(self, index):
240        """
241        Index related to this page
242        """
243        self.index_model = index
244
245    def create_default_data(self):
246        """
247        Given the user selection, creates a 1D or 2D data
248        Only when the page is on theory mode.
249        """
250        if not hasattr(self, "model_view"):
251            return
252        toggle_mode_on = self.model_view.IsEnabled()
253        if toggle_mode_on:
254            if self.enable2D and not check_data_validity(self.data):
255                self._create_default_2d_data()
256            else:
257                if self.pointsbox.GetValue():
258                    self._create_log_1d_data()
259                else:
260                    self._create_default_1d_data()
261
262            if self.model is not None:
263                if not self.data.is_data:
264                    self._manager.page_finder[self.uid].set_fit_data(
265                        data=[self.data])
266            self.on_smear_helper(update=True)
267            self.state.enable_smearer = self.enable_smearer.GetValue()
268            self.state.disable_smearer = self.disable_smearer.GetValue()
269            self.state.pinhole_smearer = self.pinhole_smearer.GetValue()
270            self.state.slit_smearer = self.slit_smearer.GetValue()
271
272    def _create_default_1d_data(self):
273        """
274        Create default data for fitting perspective
275        Only when the page is on theory mode.
276        :warning: This data is never plotted.
277
278        """
279        x = numpy.linspace(start=self.qmin_x, stop=self.qmax_x,
280                           num=self.npts_x, endpoint=True)
281        self.data = Data1D(x=x)
282        self.data.xaxis('\\rm{Q}', "A^{-1}")
283        self.data.yaxis('\\rm{Intensity}', "cm^{-1}")
284        self.data.is_data = False
285        self.data.id = str(self.uid) + " data"
286        self.data.group_id = str(self.uid) + " Model1D"
287
288    def _create_log_1d_data(self):
289        """
290        Create log-spaced data for fitting perspective
291        Only when the page is on theory mode.
292        :warning: This data is never plotted.
293
294        """
295        if self.qmin_x >= 1.e-10:
296            qmin = numpy.log10(self.qmin_x)
297        else:
298            qmin = -10.
299
300        if self.qmax_x <= 1.e10:
301            qmax = numpy.log10(self.qmax_x)
302        else:
303            qmax = 10.
304
305        x = numpy.logspace(start=qmin, stop=qmax,
306                           num=self.npts_x, endpoint=True, base=10.0)
307        self.data = Data1D(x=x)
308        self.data.xaxis('\\rm{Q}', "A^{-1}")
309        self.data.yaxis('\\rm{Intensity}', "cm^{-1}")
310        self.data.is_data = False
311        self.data.id = str(self.uid) + " data"
312        self.data.group_id = str(self.uid) + " Model1D"
313
314    def _create_default_2d_data(self):
315        """
316        Create 2D data by default
317        Only when the page is on theory mode.
318        :warning: This data is never plotted.
319        """
320        self.data = Data2D()
321        qmax = self.qmax_x / math.sqrt(2)
322        self.data.xaxis('\\rm{Q_{x}}', 'A^{-1}')
323        self.data.yaxis('\\rm{Q_{y}}', 'A^{-1}')
324        self.data.is_data = False
325        self.data.id = str(self.uid) + " data"
326        self.data.group_id = str(self.uid) + " Model2D"
327        # Default values
328        self.data.detector.append(Detector())
329        index = len(self.data.detector) - 1
330        self.data.detector[index].distance = 8000   # mm
331        self.data.source.wavelength = 6             # A
332        self.data.detector[index].pixel_size.x = 5  # mm
333        self.data.detector[index].pixel_size.y = 5  # mm
334        self.data.detector[index].beam_center.x = qmax
335        self.data.detector[index].beam_center.y = qmax
336        xmax = qmax
337        xmin = -qmax
338        ymax = qmax
339        ymin = -qmax
340        qstep = self.npts_x
341
342        x = numpy.linspace(start=xmin, stop=xmax, num=qstep, endpoint=True)
343        y = numpy.linspace(start=ymin, stop=ymax, num=qstep, endpoint=True)
344        # use data info instead
345        new_x = numpy.tile(x, (len(y), 1))
346        new_y = numpy.tile(y, (len(x), 1))
347        new_y = new_y.swapaxes(0, 1)
348        # all data reuire now in 1d array
349        qx_data = new_x.flatten()
350        qy_data = new_y.flatten()
351        q_data = numpy.sqrt(qx_data * qx_data + qy_data * qy_data)
352        # set all True (standing for unmasked) as default
353        mask = numpy.ones(len(qx_data), dtype=bool)
354        # store x and y bin centers in q space
355        x_bins = x
356        y_bins = y
357
358        self.data.source = Source()
359        self.data.data = numpy.ones(len(mask))
360        self.data.err_data = numpy.ones(len(mask))
361        self.data.qx_data = qx_data
362        self.data.qy_data = qy_data
363        self.data.q_data = q_data
364        self.data.mask = mask
365        self.data.x_bins = x_bins
366        self.data.y_bins = y_bins
367        # max and min taking account of the bin sizes
368        self.data.xmin = xmin
369        self.data.xmax = xmax
370        self.data.ymin = ymin
371        self.data.ymax = ymax
372
373    def on_set_focus(self, event):
374        """
375        On Set Focus, update guimanger and menu
376        """
377        if self._manager is not None:
378            wx.PostEvent(self._manager.parent, PanelOnFocusEvent(panel=self))
379            self.on_tap_focus()
380
381    def on_tap_focus(self):
382        """
383        Update menu1 on cliking the page tap
384        """
385        if self._manager.menu1 is not None:
386            chain_menu = self._manager.menu1.FindItemById(
387                                                   self._manager.id_reset_flag)
388            chain_menu.Enable(self.batch_on)
389            sim_menu = self._manager.menu1.FindItemById(self._manager.id_simfit)
390            flag = self.data.is_data\
391                            and (self.model is not None)
392            sim_menu.Enable(not self.batch_on and flag)
393            batch_menu = \
394                    self._manager.menu1.FindItemById(self._manager.id_batchfit)
395            batch_menu.Enable(self.batch_on and flag)
396
397    def onContextMenu(self, event):
398        """
399        Retrieve the state selected state
400        """
401        pos = event.GetPosition()
402        pos = self.ScreenToClient(pos)
403        self.PopupMenu(self.popUpMenu, pos)
404
405    def onUndo(self, event):
406        """
407        Cancel the previous action
408        """
409        event = PreviousStateEvent(page=self)
410        wx.PostEvent(self.parent, event)
411
412    def onRedo(self, event):
413        """
414        Restore the previous action cancelled
415        """
416        event = NextStateEvent(page=self)
417        wx.PostEvent(self.parent, event)
418
419    def define_page_structure(self):
420        """
421        Create empty sizer for a panel
422        """
423        self.vbox = wx.BoxSizer(wx.VERTICAL)
424        self.sizer0 = wx.BoxSizer(wx.VERTICAL)
425        self.sizer1 = wx.BoxSizer(wx.VERTICAL)
426        self.sizer2 = wx.BoxSizer(wx.VERTICAL)
427        self.sizer3 = wx.BoxSizer(wx.VERTICAL)
428        self.sizer4 = wx.BoxSizer(wx.VERTICAL)
429        self.sizer5 = wx.BoxSizer(wx.VERTICAL)
430        self.sizer6 = wx.BoxSizer(wx.VERTICAL)
431
432        self.sizer0.SetMinSize((PANEL_WIDTH, -1))
433        self.sizer1.SetMinSize((PANEL_WIDTH, -1))
434        self.sizer2.SetMinSize((PANEL_WIDTH, -1))
435        self.sizer3.SetMinSize((PANEL_WIDTH, -1))
436        self.sizer4.SetMinSize((PANEL_WIDTH, -1))
437        self.sizer5.SetMinSize((PANEL_WIDTH, -1))
438        self.sizer6.SetMinSize((PANEL_WIDTH, -1))
439
440        self.vbox.Add(self.sizer0)
441        self.vbox.Add(self.sizer1)
442        self.vbox.Add(self.sizer2)
443        self.vbox.Add(self.sizer3)
444        self.vbox.Add(self.sizer4)
445        self.vbox.Add(self.sizer5)
446        self.vbox.Add(self.sizer6)
447
448    def set_layout(self):
449        """
450        layout
451        """
452        self.vbox.Layout()
453        self.vbox.Fit(self)
454        self.SetSizer(self.vbox)
455        self.Centre()
456
457    def set_owner(self, owner):
458        """
459        set owner of fitpage
460
461        :param owner: the class responsible of plotting
462
463        """
464        self.event_owner = owner
465        self.state.event_owner = owner
466
467    def get_state(self):
468        """
469        return the current page state
470        """
471        return self.state
472
473    def get_data(self):
474        """
475        return the current data
476        """
477        return self.data
478
479    def get_data_list(self):
480        """
481        return the current data
482        """
483        return self.data_list
484
485    def set_manager(self, manager):
486        """
487        set panel manager
488
489        :param manager: instance of plugin fitting
490
491        """
492        self._manager = manager
493        self.state.manager = manager
494
495    def populate_box(self, model_dict):
496        """
497        Store list of model
498
499        :param model_dict: dictionary containing list of models
500
501        """
502        self.model_list_box = model_dict
503        self.state.model_list_box = self.model_list_box
504        self.initialize_combox()
505
506    def set_model_dictionary(self, model_dict):
507        """
508        Store a dictionary linking model name -> model object
509
510        :param model_dict: dictionary containing list of models
511        """
512        self.model_dict = model_dict
513
514    def initialize_combox(self):
515        """
516        put default value in the combo box
517        """
518        if self.model_list_box is not None and len(self.model_list_box) > 0:
519            self._populate_box(self.structurebox,
520                               self.model_list_box["Structure Factors"])
521            self.structurebox.Insert("None", 0, None)
522            self.structurebox.SetSelection(0)
523            self.structurebox.Hide()
524            self.text2.Hide()
525            self.structurebox.Disable()
526            self.text2.Disable()
527
528    def set_dispers_sizer(self):
529        """
530        fill sizer containing dispersity info
531        """
532        # print "==== entering set_dispers_sizer ==="
533        self.sizer4.Clear(True)
534        name = "Polydispersity and Orientational Distribution"
535        box_description = wx.StaticBox(self, wx.ID_ANY, name)
536        box_description.SetForegroundColour(wx.BLUE)
537        boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
538        # ----------------------------------------------------
539        self.disable_disp = wx.RadioButton(self, wx.ID_ANY, 'Off', (10, 10),
540                                           style=wx.RB_GROUP)
541        self.enable_disp = wx.RadioButton(self, wx.ID_ANY, 'On', (10, 30))
542        # best size for MAC and PC
543        if ON_MAC:
544            size_q = (30, 20)
545        else:
546            size_q = (20, 15)
547        self.disp_help_bt = wx.Button(self, self.ID_DISPERSER_HELP, '?',
548                                      style=wx.BU_EXACTFIT,
549                                      size=size_q)
550        self.disp_help_bt.Bind(wx.EVT_BUTTON, self.on_pd_help_clicked,
551                               id=self.disp_help_bt.GetId())
552        self.disp_help_bt.SetToolTipString("Help for polydispersion.")
553
554        self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param,
555                  id=self.disable_disp.GetId())
556        self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param,
557                  id=self.enable_disp.GetId())
558        # MAC needs SetValue
559        self.disable_disp.SetValue(True)
560        sizer_dispersion = wx.BoxSizer(wx.HORIZONTAL)
561        sizer_dispersion.Add((20, 20))
562        name = ""  # Polydispersity and \nOrientational Distribution "
563        sizer_dispersion.Add(wx.StaticText(self, wx.ID_ANY, name))
564        sizer_dispersion.Add(self.enable_disp)
565        sizer_dispersion.Add((20, 20))
566        sizer_dispersion.Add(self.disable_disp)
567        sizer_dispersion.Add((25, 20))
568        sizer_dispersion.Add(self.disp_help_bt)
569
570        # fill a sizer for dispersion
571        boxsizer1.Add(sizer_dispersion, 0,
572                      wx.TOP|wx.BOTTOM|wx.LEFT|wx.EXPAND|wx.ADJUST_MINSIZE,
573                      border=5)
574        self.sizer4_4 = wx.GridBagSizer(6, 5)
575
576        boxsizer1.Add(self.sizer4_4)
577        # -----------------------------------------------------
578        self.sizer4.Add(boxsizer1, 0, wx.EXPAND | wx.ALL, 10)
579        self.sizer4_4.Layout()
580        self.sizer4.Layout()
581        self.Layout()
582
583        self.Refresh()
584        # saving the state of enable dispersity button
585        self.state.enable_disp = self.enable_disp.GetValue()
586        self.state.disable_disp = self.disable_disp.GetValue()
587        self.SetupScrolling()
588
589    def onResetModel(self, event):
590        """
591        Reset model state
592        """
593        menu = event.GetEventObject()
594        # post help message for the selected model
595        msg = menu.GetHelpString(event.GetId())
596        msg += " reloaded"
597        wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
598        self.Show(False)
599        name = menu.GetLabel(event.GetId())
600        self._on_select_model_helper()
601        if self.model is not None:
602            self.m_name = self.model.name
603        if name in self.saved_states.keys():
604            previous_state = self.saved_states[name]
605            # reset state of checkbox,textcrtl  and  regular parameters value
606
607            self.reset_page(previous_state)
608        self.state.m_name = self.m_name
609        self.Show(True)
610
611    def on_preview(self, event):
612        """
613        Report the current fit results
614        """
615        # Get plot image from plotpanel
616        images, canvases = self.get_images()
617        # get the report dialog
618        self.state.report(images, canvases)
619
620    def on_save(self, event):
621        """
622        Save the current state into file
623        """
624        self.save_current_state()
625        new_state = self.state.clone()
626        # Ask the user the location of the file to write to.
627        path = None
628        if self.parent is not None:
629            self._default_save_location = \
630                        self._manager.parent._default_save_location
631        dlg = wx.FileDialog(self, "Choose a file", self._default_save_location,
632                            self.window_caption, "*.fitv", wx.SAVE)
633
634        if dlg.ShowModal() == wx.ID_OK:
635            path = dlg.GetPath()
636            self._default_save_location = os.path.dirname(path)
637            self._manager.parent._default_save_location = \
638                self._default_save_location
639        else:
640            return None
641        # MAC always needs the extension for saving
642        extens = ".fitv"
643        # Make sure the ext included in the file name
644        fName = os.path.splitext(path)[0] + extens
645        # the manager write the state into file
646        self._manager.save_fit_state(filepath=fName, fitstate=new_state)
647        return new_state
648
649    def on_copy(self, event):
650        """
651        Copy Parameter values to the clipboad
652        """
653        if event is not None:
654            event.Skip()
655        # It seems MAC needs wxCallAfter
656        if event.GetId() == GUIFRAME_ID.COPYEX_ID:
657            print "copy excel"
658            wx.CallAfter(self.get_copy_excel)
659        elif event.GetId() == GUIFRAME_ID.COPYLAT_ID:
660            print "copy latex"
661            wx.CallAfter(self.get_copy_latex)
662        else:
663            wx.CallAfter(self.get_copy)
664
665    def on_paste(self, event):
666        """
667        Paste Parameter values to the panel if possible
668        """
669        # if event is not None:
670        #    event.Skip()
671        # It seems MAC needs wxCallAfter for the setvalues
672        # for multiple textctrl items, otherwise it tends to crash once a while
673        wx.CallAfter(self.get_paste)
674        # messages depending on the flag
675        # self._copy_info(True)
676
677    def _copy_info(self, flag):
678        """
679        Send event dpemding on flag
680
681        : Param flag: flag that distinguish event
682        """
683        # messages depending on the flag
684        if flag is None:
685            msg = " Parameter values are copied to the clipboard..."
686            infor = 'warning'
687        elif flag:
688            msg = " Parameter values are pasted from the clipboard..."
689            infor = "warning"
690        else:
691            msg = "Error occurred: "
692            msg += "No valid parameter values to paste from the clipboard..."
693            infor = "warning"
694        # inform msg to wx
695        wx.PostEvent(self._manager.parent,
696                     StatusEvent(status=msg, info=infor))
697
698    def _get_time_stamp(self):
699        """
700        return time and date stings
701        """
702        # date and time
703        year, month, day, hour, minute, second, _, _, _ = time.localtime()
704        current_time = str(hour) + ":" + str(minute) + ":" + str(second)
705        current_date = str(month) + "/" + str(day) + "/" + str(year)
706        return current_time, current_date
707
708    def on_bookmark(self, event):
709        """
710        save history of the data and model
711        """
712        if self.model is None:
713            msg = "Can not bookmark; Please select Data and Model first..."
714            wx.MessageBox(msg, 'Info')
715            return
716        self.save_current_state()
717        new_state = self.state.clone()
718        # Add model state on context menu
719        self.number_saved_state += 1
720        current_time, current_date = self._get_time_stamp()
721        # name= self.model.name+"[%g]"%self.number_saved_state
722        name = "Fitting: %g]" % self.number_saved_state
723        name += self.model.__class__.__name__
724        name += "bookmarked at %s on %s" % (current_time, current_date)
725        self.saved_states[name] = new_state
726
727        # Add item in the context menu
728        msg = "Model saved at %s on %s" % (current_time, current_date)
729        # post help message for the selected model
730        msg += " Saved! right click on this page to retrieve this model"
731        wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
732
733        self.popUpMenu.Append(self.ID_BOOKMARK, name, str(msg))
734        wx.EVT_MENU(self, self.ID_BOOKMARK, self.onResetModel)
735        wx.PostEvent(self._manager.parent,
736                     AppendBookmarkEvent(title=name,
737                                         hint=str(msg),
738                                         handler=self._back_to_bookmark))
739
740    def _back_to_bookmark(self, event):
741        """
742        Back to bookmark
743        """
744        self._manager.on_perspective(event)
745        self.onResetModel(event)
746        self._draw_model()
747
748    def onSetFocus(self, evt):
749        """
750        highlight the current textcrtl and hide the error text control shown
751        after fitting
752        """
753        return
754
755    def read_file(self, path):
756        """
757        Read two columns file
758
759        :param path: the path to the file to read
760
761        """
762        try:
763            if path is None:
764                status = " Selected Distribution was not loaded: %s" % path
765                wx.PostEvent(self._manager.parent,
766                             StatusEvent(status=status))
767                return None, None
768            input_f = open(path, 'r')
769            buff = input_f.read()
770            lines = buff.split('\n')
771            input_f.close()
772            angles = []
773            weights = []
774            for line in lines:
775                toks = line.split()
776                try:
777                    angle = float(toks[0])
778                    weight = float(toks[1])
779                    angles.append(angle)
780                    weights.append(weight)
781                except Exception:
782                    # Skip non-data lines
783                    logging.error(traceback.format_exc())
784            return numpy.array(angles), numpy.array(weights)
785        except:
786            raise
787
788    def createMemento(self):
789        """
790        return the current state of the page
791        """
792        return self.state.clone()
793
794    def save_current_state(self):
795        """
796        Store current state
797        """
798        # save model option
799        if self.model is not None:
800            self.disp_list = self.model.getDispParamList()
801            self.state.disp_list = copy.deepcopy(self.disp_list)
802            self.state.model = self.model.clone()
803
804            # model combobox: complex code because of mac's silent error
805            if self.structurebox is not None:
806                if self.structurebox.IsShown():
807                    self.state.structurecombobox = 'None'
808                    s_select = self.structurebox.GetSelection()
809                    if s_select > 0:
810                        self.state.structurecombobox = \
811                            self.structurebox.GetString(s_select)
812            if self.formfactorbox is not None:
813                f_select = self.formfactorbox.GetSelection()
814                if f_select > 0:
815                    self.state.formfactorcombobox = \
816                        self.formfactorbox.GetString(f_select)
817        if self.categorybox is not None:
818            cb_select = self.categorybox.GetSelection()
819            if cb_select > 0:
820                self.state.categorycombobox = \
821                    self.categorybox.GetString(cb_select)
822
823        self.state.enable2D = copy.deepcopy(self.enable2D)
824        self.state.values = copy.deepcopy(self.values)
825        self.state.weights = copy.deepcopy(self.weights)
826        # save data
827        self.state.data = copy.deepcopy(self.data)
828        self.state.qmax_x = self.qmax_x
829        self.state.qmin_x = self.qmin_x
830        self.state.dI_noweight = copy.deepcopy(self.dI_noweight.GetValue())
831        self.state.dI_didata = copy.deepcopy(self.dI_didata.GetValue())
832        self.state.dI_sqrdata = copy.deepcopy(self.dI_sqrdata.GetValue())
833        self.state.dI_idata = copy.deepcopy(self.dI_idata.GetValue())
834        self.state.dq_l = self.dq_l
835        self.state.dq_r = self.dq_r
836        if hasattr(self, "enable_disp"):
837            self.state.enable_disp = self.enable_disp.GetValue()
838            self.state.disable_disp = self.disable_disp.GetValue()
839
840        self.state.smearer = copy.deepcopy(self.current_smearer)
841        if hasattr(self, "enable_smearer"):
842            self.state.enable_smearer = \
843                                copy.deepcopy(self.enable_smearer.GetValue())
844            self.state.disable_smearer = \
845                                copy.deepcopy(self.disable_smearer.GetValue())
846
847        self.state.pinhole_smearer = \
848                                copy.deepcopy(self.pinhole_smearer.GetValue())
849        self.state.dx_max = copy.deepcopy(self.dx_max)
850        self.state.dx_min = copy.deepcopy(self.dx_min)
851        self.state.dxl = copy.deepcopy(self.dxl)
852        self.state.dxw = copy.deepcopy(self.dxw)
853        self.state.slit_smearer = copy.deepcopy(self.slit_smearer.GetValue())
854
855        if len(self._disp_obj_dict) > 0:
856            for k, v in self._disp_obj_dict.iteritems():
857                self.state._disp_obj_dict[k] = v.type
858
859            self.state.values = copy.deepcopy(self.values)
860            self.state.weights = copy.deepcopy(self.weights)
861        # save plotting range
862        self._save_plotting_range()
863
864        self.state.orientation_params = []
865        self.state.orientation_params_disp = []
866        self.state.parameters = []
867        self.state.fittable_param = []
868        self.state.fixed_param = []
869        self.state.str_parameters = []
870
871        # save checkbutton state and txtcrtl values
872        self._copy_parameters_state(self.str_parameters,
873                                    self.state.str_parameters)
874        self._copy_parameters_state(self.orientation_params,
875                                     self.state.orientation_params)
876        self._copy_parameters_state(self.orientation_params_disp,
877                                    self.state.orientation_params_disp)
878
879        self._copy_parameters_state(self.parameters, self.state.parameters)
880        self._copy_parameters_state(self.fittable_param,
881                                    self.state.fittable_param)
882        self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
883        # save chisqr
884        self.state.tcChi = self.tcChi.GetValue()
885
886    def save_current_state_fit(self):
887        """
888        Store current state for fit_page
889        """
890        # save model option
891        if self.model is not None:
892            self.disp_list = self.model.getDispParamList()
893            self.state.disp_list = copy.deepcopy(self.disp_list)
894            self.state.model = self.model.clone()
895
896        self.state.enable2D = copy.deepcopy(self.enable2D)
897        self.state.values = copy.deepcopy(self.values)
898        self.state.weights = copy.deepcopy(self.weights)
899        # save data
900        self.state.data = copy.deepcopy(self.data)
901
902        if hasattr(self, "enable_disp"):
903            self.state.enable_disp = self.enable_disp.GetValue()
904            self.state.disable_disp = self.disable_disp.GetValue()
905
906        self.state.smearer = copy.deepcopy(self.current_smearer)
907        if hasattr(self, "enable_smearer"):
908            self.state.enable_smearer = \
909                                copy.deepcopy(self.enable_smearer.GetValue())
910            self.state.disable_smearer = \
911                                copy.deepcopy(self.disable_smearer.GetValue())
912
913        self.state.pinhole_smearer = \
914                                copy.deepcopy(self.pinhole_smearer.GetValue())
915        self.state.slit_smearer = copy.deepcopy(self.slit_smearer.GetValue())
916        self.state.dI_noweight = copy.deepcopy(self.dI_noweight.GetValue())
917        self.state.dI_didata = copy.deepcopy(self.dI_didata.GetValue())
918        self.state.dI_sqrdata = copy.deepcopy(self.dI_sqrdata.GetValue())
919        self.state.dI_idata = copy.deepcopy(self.dI_idata.GetValue())
920        if hasattr(self, "disp_box") and self.disp_box is not None:
921            self.state.disp_box = self.disp_box.GetCurrentSelection()
922
923            if len(self.disp_cb_dict) > 0:
924                for k, v in self.disp_cb_dict.iteritems():
925                    if v is None:
926                        self.state.disp_cb_dict[k] = v
927                    else:
928                        try:
929                            self.state.disp_cb_dict[k] = v.GetValue()
930                        except:
931                            self.state.disp_cb_dict[k] = None
932            if len(self._disp_obj_dict) > 0:
933                for k, v in self._disp_obj_dict.iteritems():
934                    self.state._disp_obj_dict[k] = v.type
935
936            self.state.values = copy.deepcopy(self.values)
937            self.state.weights = copy.deepcopy(self.weights)
938
939        # save plotting range
940        self._save_plotting_range()
941
942        # save checkbutton state and txtcrtl values
943        self._copy_parameters_state(self.orientation_params,
944                                    self.state.orientation_params)
945        self._copy_parameters_state(self.orientation_params_disp,
946                                    self.state.orientation_params_disp)
947        self._copy_parameters_state(self.parameters, self.state.parameters)
948        self._copy_parameters_state(self.fittable_param,
949                                    self.state.fittable_param)
950        self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
951
952    def check_invalid_panel(self):
953        """
954        check if the user can already perform some action with this panel
955        """
956        if self.data is None:
957            self.disable_smearer.SetValue(True)
958            self.disable_disp.SetValue(True)
959            msg = "Please load Data and select Model to start..."
960            wx.MessageBox(msg, 'Info')
961            return True
962
963    def set_model_state(self, state):
964        """
965        reset page given a model state
966        """
967        self.disp_cb_dict = state.disp_cb_dict
968        self.disp_list = state.disp_list
969
970        # fill model combobox
971        self._show_combox_helper()
972        # select the current model
973        try:
974            # to support older version
975            category_pos = int(state.categorycombobox)
976        except:
977            category_pos = 0
978            for ind_cat in range(self.categorybox.GetCount()):
979                if self.categorycombobox.GetString(ind_cat) == \
980                                        state.categorycombobox:
981                    category_pos = int(ind_cat)
982                    break
983
984        self.categorybox.Select(category_pos)
985        try:
986            # to support older version
987            formfactor_pos = int(state.formfactorcombobox)
988        except:
989            formfactor_pos = 0
990            for ind_form in range(self.formfactorbox.GetCount()):
991                if self.formfactorbox.GetString(ind_form) == \
992                                        state.formfactorcombobox:
993                    formfactor_pos = int(ind_form)
994                    break
995
996        self.formfactorbox.Select(formfactor_pos)
997
998        try:
999            # to support older version
1000            structfactor_pos = int(state.structurecombobox)
1001        except:
1002            structfactor_pos = 0
1003            for ind_struct in range(self.structurebox.GetCount()):
1004                if self.structurebox.GetString(ind_struct) == \
1005                                        state.structurecombobox:
1006                    structfactor_pos = int(ind_struct)
1007                    break
1008
1009        self.structurebox.SetSelection(structfactor_pos)
1010
1011        if state.multi_factor is not None:
1012            self.multifactorbox.SetSelection(state.multi_factor)
1013
1014        # reset state of checkbox,textcrtl  and  regular parameters value
1015        self._reset_parameters_state(self.orientation_params_disp,
1016                                     state.orientation_params_disp)
1017        self._reset_parameters_state(self.orientation_params,
1018                                     state.orientation_params)
1019        self._reset_parameters_state(self.str_parameters,
1020                                     state.str_parameters)
1021        self._reset_parameters_state(self.parameters, state.parameters)
1022        # display dispersion info layer
1023        self.enable_disp.SetValue(state.enable_disp)
1024        self.disable_disp.SetValue(state.disable_disp)
1025
1026        if hasattr(self, "disp_box") and self.disp_box is not None:
1027            self.disp_box.SetSelection(state.disp_box)
1028            n = self.disp_box.GetCurrentSelection()
1029            dispersity = self.disp_box.GetClientData(n)
1030            name = dispersity.__name__
1031            self._set_dipers_Param(event=None)
1032
1033            if name == "ArrayDispersion":
1034
1035                for item in self.disp_cb_dict.keys():
1036
1037                    if hasattr(self.disp_cb_dict[item], "SetValue"):
1038                        self.disp_cb_dict[item].SetValue(
1039                                                    state.disp_cb_dict[item])
1040                        # Create the dispersion objects
1041                        disp_model = POLYDISPERSITY_MODELS['array']()
1042                        if hasattr(state, "values") and \
1043                                 self.disp_cb_dict[item].GetValue() is True:
1044                            if len(state.values) > 0:
1045                                self.values = state.values
1046                                self.weights = state.weights
1047                                disp_model.set_weights(self.values,
1048                                                       state.weights)
1049                            else:
1050                                self._reset_dispersity()
1051
1052                        self._disp_obj_dict[item] = disp_model
1053                        # Set the new model as the dispersion object
1054                        # for the selected parameter
1055                        self.model.set_dispersion(item, disp_model)
1056
1057                        self.model._persistency_dict[item] = \
1058                                                [state.values, state.weights]
1059
1060            else:
1061                keys = self.model.getParamList()
1062                for item in keys:
1063                    if item in self.disp_list and \
1064                            item not in self.model.details:
1065                        self.model.details[item] = ["", None, None]
1066                self.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1067                self.state.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1068        # smearing info  restore
1069        if hasattr(self, "enable_smearer"):
1070            # set smearing value whether or not the data
1071            # contain the smearing info
1072            self.enable_smearer.SetValue(state.enable_smearer)
1073            self.disable_smearer.SetValue(state.disable_smearer)
1074            self.onSmear(event=None)
1075        self.pinhole_smearer.SetValue(state.pinhole_smearer)
1076        self.slit_smearer.SetValue(state.slit_smearer)
1077
1078        self.dI_noweight.SetValue(state.dI_noweight)
1079        self.dI_didata.SetValue(state.dI_didata)
1080        self.dI_sqrdata.SetValue(state.dI_sqrdata)
1081        self.dI_idata.SetValue(state.dI_idata)
1082
1083        # we have two more options for smearing
1084        if self.pinhole_smearer.GetValue():
1085            self.onPinholeSmear(event=None)
1086        elif self.slit_smearer.GetValue():
1087            self.onSlitSmear(event=None)
1088
1089        # reset state of checkbox,textcrtl  and dispersity parameters value
1090        self._reset_parameters_state(self.fittable_param, state.fittable_param)
1091        self._reset_parameters_state(self.fixed_param, state.fixed_param)
1092
1093        # draw the model with previous parameters value
1094        self._onparamEnter_helper()
1095        self.select_param(event=None)
1096        # Save state_fit
1097        self.save_current_state_fit()
1098        self._lay_out()
1099        self.Refresh()
1100
1101    def get_cat_combo_box_pos(self, state):
1102        """
1103        Iterate through the categories to find the structurefactor
1104        :return: combo_box_position
1105        """
1106        for key, value in self.master_category_dict.iteritems():
1107            for list_item in value:
1108                if state.formfactorcombobox in list_item:
1109                    return self.categorybox.Items.index(key)
1110        return 0
1111
1112    def reset_page_helper(self, state):
1113        """
1114        Use page_state and change the state of existing page
1115
1116        :precondition: the page is already drawn or created
1117
1118        :postcondition: the state of the underlying data change as well as the
1119            state of the graphic interface
1120        """
1121        if state is None:
1122            return
1123        # set data, etc. from the state
1124        # reset page between theory and fitting from bookmarking
1125        data = state.data
1126
1127        if data is None:
1128            data_min = state.qmin
1129            data_max = state.qmax
1130            self.qmin_x = data_min
1131            self.qmax_x = data_max
1132            self.qmin.SetValue(str(data_min))
1133            self.qmax.SetValue(str(data_max))
1134
1135            self.state.data = data
1136            self.state.qmin = self.qmin_x
1137            self.state.qmax = self.qmax_x
1138        else:
1139            self.set_data(data)
1140
1141        self.enable2D = state.enable2D
1142        try:
1143            self.magnetic_on = state.magnetic_on
1144        except:
1145            # Backward compatibility (for older state files)
1146            self.magnetic_on = False
1147
1148        self.disp_cb_dict = state.disp_cb_dict
1149        self.disp_list = state.disp_list
1150
1151        # fill model combobox
1152        self._show_combox_helper()
1153        # select the current model
1154        try:
1155            # to support older version
1156            category_pos = int(state.categorycombobox)
1157        except:
1158            state.formfactorcombobox = state.formfactorcombobox.lower()
1159            state.formfactorcombobox = \
1160                state.formfactorcombobox.replace('model', '')
1161            state.formfactorcombobox = unicode(state.formfactorcombobox)
1162            state.categorycombobox = unicode(state.categorycombobox)
1163            if state.categorycombobox in self.categorybox.Items:
1164                category_pos = self.categorybox.Items.index(
1165                    state.categorycombobox)
1166            else:
1167                # Look in master list for model name (model.lower)
1168                category_pos = self.get_cat_combo_box_pos(state)
1169
1170        self.categorybox.Select(category_pos)
1171        self._show_combox(None)
1172        try:
1173            # to support older version
1174            formfactor_pos = int(state.formfactorcombobox)
1175        except:
1176            formfactor_pos = 0
1177            for ind_form in range(self.formfactorbox.GetCount()):
1178                if self.formfactorbox.GetString(ind_form) == \
1179                                                    (state.formfactorcombobox):
1180                    formfactor_pos = int(ind_form)
1181                    break
1182
1183        self.formfactorbox.Select(formfactor_pos)
1184
1185        structfactor_pos = 0
1186        try:
1187            # to support older version
1188            structfactor_pos = int(state.structurecombobox)
1189        except:
1190            if state.structurecombobox is not None:
1191                state.structurecombobox = unicode(state.structurecombobox)
1192                for ind_struct in range(self.structurebox.GetCount()):
1193                    if self.structurebox.GetString(ind_struct) == \
1194                                                    (state.structurecombobox):
1195                        structfactor_pos = int(ind_struct)
1196                        break
1197
1198        self.structurebox.SetSelection(structfactor_pos)
1199
1200        if state.multi_factor is not None:
1201            self.multifactorbox.SetSelection(state.multi_factor)
1202
1203        # draw the panel according to the new model parameter
1204        self._on_select_model(event=None)
1205
1206        # take care of 2D button
1207        if data is None and self.model_view.IsEnabled():
1208            if self.enable2D:
1209                self.model_view.SetLabel("2D Mode")
1210            else:
1211                self.model_view.SetLabel("1D Mode")
1212
1213        # reset state of checkbox,textcrtl  and  regular parameters value
1214        self._reset_parameters_state(self.orientation_params_disp,
1215                                     state.orientation_params_disp)
1216        self._reset_parameters_state(self.orientation_params,
1217                                     state.orientation_params)
1218        self._reset_parameters_state(self.str_parameters,
1219                                     state.str_parameters)
1220        self._reset_parameters_state(self.parameters, state.parameters)
1221        # display dispersion info layer
1222        self.enable_disp.SetValue(state.enable_disp)
1223        self.disable_disp.SetValue(state.disable_disp)
1224        # If the polydispersion is ON
1225        if state.enable_disp:
1226            # reset dispersion according the state
1227            self._set_dipers_Param(event=None)
1228            self._reset_page_disp_helper(state)
1229        # plotting range restore
1230        self._reset_plotting_range(state)
1231        # smearing info  restore
1232        if hasattr(self, "enable_smearer"):
1233            # set smearing value whether or not the data
1234            # contain the smearing info
1235            self.enable_smearer.SetValue(state.enable_smearer)
1236            self.disable_smearer.SetValue(state.disable_smearer)
1237            self.onSmear(event=None)
1238        self.pinhole_smearer.SetValue(state.pinhole_smearer)
1239        self.slit_smearer.SetValue(state.slit_smearer)
1240        try:
1241            self.dI_noweight.SetValue(state.dI_noweight)
1242            self.dI_didata.SetValue(state.dI_didata)
1243            self.dI_sqrdata.SetValue(state.dI_sqrdata)
1244            self.dI_idata.SetValue(state.dI_idata)
1245        except:
1246            # to support older state file formats
1247            self.dI_noweight.SetValue(False)
1248            self.dI_didata.SetValue(True)
1249            self.dI_sqrdata.SetValue(False)
1250            self.dI_idata.SetValue(False)
1251
1252        # we have two more options for smearing
1253        if self.pinhole_smearer.GetValue():
1254            self.dx_min = state.dx_min
1255            self.dx_max = state.dx_max
1256            if self.dx_min is not None:
1257                self.smear_pinhole_min.SetValue(str(self.dx_min))
1258            if self.dx_max is not None:
1259                self.smear_pinhole_max.SetValue(str(self.dx_max))
1260            self.onPinholeSmear(event=None)
1261        elif self.slit_smearer.GetValue():
1262            self.dxl = state.dxl
1263            self.dxw = state.dxw
1264            if self.dxl is not None:
1265                self.smear_slit_height.SetValue(str(self.dxl))
1266            if self.dxw is not None:
1267                self.smear_slit_width.SetValue(str(self.dxw))
1268            else:
1269                self.smear_slit_width.SetValue('')
1270            self.onSlitSmear(event=None)
1271
1272        # reset state of checkbox,textcrtl  and dispersity parameters value
1273        self._reset_parameters_state(self.fittable_param, state.fittable_param)
1274        self._reset_parameters_state(self.fixed_param, state.fixed_param)
1275
1276        # draw the model with previous parameters value
1277        self._onparamEnter_helper()
1278        # reset the value of chisqr when not consistent with the value computed
1279        self.tcChi.SetValue(str(self.state.tcChi))
1280        # reset context menu items
1281        self._reset_context_menu()
1282
1283        # set the value of the current state to the state given as parameter
1284        self.state = state.clone()
1285        self.state.m_name = self.m_name
1286
1287    def _reset_page_disp_helper(self, state):
1288        """
1289        Help to rest page for dispersions
1290        """
1291        keys = self.model.getParamList()
1292        for item in keys:
1293            if item in self.disp_list and \
1294                            item not in self.model.details:
1295                self.model.details[item] = ["", None, None]
1296        # for k,v in self.state.disp_cb_dict.iteritems():
1297        self.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1298        self.state.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1299        self.values = copy.deepcopy(state.values)
1300        self.weights = copy.deepcopy(state.weights)
1301
1302        for key, disp_type in state._disp_obj_dict.iteritems():
1303            # disp_model = disp
1304            disp_model = POLYDISPERSITY_MODELS[disp_type]()
1305            self._disp_obj_dict[key] = disp_model
1306            param_name = key.split('.')[0]
1307            # Try to set dispersion only when available
1308            # for eg., pass the orient. angles for 1D Cal
1309            try:
1310                self.model.set_dispersion(param_name, disp_model)
1311                self.model._persistency_dict[key] = \
1312                    [state.values, state.weights]
1313            except Exception:
1314                logging.error(traceback.format_exc())
1315            selection = self._find_polyfunc_selection(disp_model)
1316            for list in self.fittable_param:
1317                if list[1] == key and list[7] is not None:
1318                    list[7].SetSelection(selection)
1319                    # For the array disp_model, set the values and weights
1320                    if selection == 1:
1321                        disp_model.set_weights(self.values[key],
1322                                               self.weights[key])
1323                        try:
1324                            # Diables all fittable params for array
1325                            list[0].SetValue(False)
1326                            list[0].Disable()
1327                            list[2].Disable()
1328                            list[5].Disable()
1329                            list[6].Disable()
1330                        except Exception:
1331                            logging.error(traceback.format_exc())
1332            # For array, disable all fixed params
1333            if selection == 1:
1334                for item in self.fixed_param:
1335                    if item[1].split(".")[0] == key.split(".")[0]:
1336                        # try it and pass it for the orientation for 1D
1337                        try:
1338                            item[2].Disable()
1339                        except Exception:
1340                            logging.error(traceback.format_exc())
1341
1342    def _selectDlg(self):
1343        """
1344        open a dialog file to selected the customized dispersity
1345        """
1346        if self.parent is not None:
1347            self._default_save_location = \
1348                        self._manager.parent.get_save_location()
1349        dlg = wx.FileDialog(self, "Choose a weight file",
1350                            self._default_save_location, "",
1351                            "*.*", wx.OPEN)
1352        path = None
1353        if dlg.ShowModal() == wx.ID_OK:
1354            path = dlg.GetPath()
1355        dlg.Destroy()
1356        return path
1357
1358    def _reset_context_menu(self):
1359        """
1360        reset the context menu
1361        """
1362        ids = iter(self._id_pool)  # Reusing ids for context menu
1363        for name, _ in self.state.saved_states.iteritems():
1364            self.number_saved_state += 1
1365            # Add item in the context menu
1366            wx_id = ids.next()
1367            msg = 'Save model and state %g' % self.number_saved_state
1368            self.popUpMenu.Append(wx_id, name, msg)
1369            wx.EVT_MENU(self, wx_id, self.onResetModel)
1370
1371    def _reset_plotting_range(self, state):
1372        """
1373        Reset the plotting range to a given state
1374        """
1375        self.qmin.SetValue(str(state.qmin))
1376        self.qmax.SetValue(str(state.qmax))
1377
1378    def _save_typeOfmodel(self):
1379        """
1380        save radiobutton containing the type model that can be selected
1381        """
1382        # self.state.shape_rbutton = self.shape_rbutton.GetValue()
1383        # self.state.shape_indep_rbutton = self.shape_indep_rbutton.GetValue()
1384        # self.state.struct_rbutton = self.struct_rbutton.GetValue()
1385        # self.state.plugin_rbutton = self.plugin_rbutton.GetValue()
1386        self.state.structurecombobox = self.structurebox.GetLabel()
1387        self.state.formfactorcombobox = self.formfactorbox.GetLabel()
1388        self.state.categorycombobox = self.categorybox.GetLabel()
1389
1390        # post state to fit panel
1391        event = PageInfoEvent(page=self)
1392        wx.PostEvent(self.parent, event)
1393
1394    def _save_plotting_range(self):
1395        """
1396        save the state of plotting range
1397        """
1398        self.state.qmin = self.qmin_x
1399        self.state.qmax = self.qmax_x
1400        self.state.npts = self.npts_x
1401
1402    def _onparamEnter_helper(self, is_modified=False):
1403        """
1404        check if values entered by the user are changed and valid to replot
1405        model
1406        """
1407        # Flag to register when a parameter has changed.
1408        # is_modified = False
1409        self.fitrange = True
1410        is_2Ddata = False
1411        # self._undo.Enable(True)
1412        # check if 2d data
1413        if self.data.__class__.__name__ == "Data2D":
1414            is_2Ddata = True
1415        if self.model is not None:
1416            # Either we get a is_modified = True passed in because
1417            # _update_paramv_on_fit() has been called already or
1418            # we need to check here ourselves.
1419            if not is_modified:
1420                is_modified = (self._check_value_enter(self.fittable_param)
1421                               or self._check_value_enter(self.fixed_param)
1422                               or self._check_value_enter(self.parameters))
1423
1424            # Here we should check whether the boundaries have been modified.
1425            # If qmin and qmax have been modified, update qmin and qmax and
1426            # set the is_modified flag to True
1427            if self._validate_qrange(self.qmin, self.qmax):
1428                tempmin = float(self.qmin.GetValue())
1429                if tempmin != self.qmin_x:
1430                    self.qmin_x = tempmin
1431                    is_modified = True
1432                tempmax = float(self.qmax.GetValue())
1433                if tempmax != self.qmax_x:
1434                    self.qmax_x = tempmax
1435                    is_modified = True
1436                if is_2Ddata:
1437                    is_modified = self._validate_Npts()
1438                else:
1439                    is_modified = self._validate_Npts_1D()
1440            else:
1441                self.fitrange = False
1442
1443            # if any value is modify draw model with new value
1444            if not self.fitrange:
1445                # self.btFit.Disable()
1446                if is_2Ddata:
1447                    self.btEditMask.Disable()
1448            else:
1449                if is_2Ddata and self.data.is_data and not self.batch_on:
1450                    self.btEditMask.Enable(True)
1451            if is_modified and self.fitrange:
1452                # Theory case: need to get npts value to draw
1453                self.npts_x = float(self.Npts_total.GetValue())
1454                self.Npts_fit.SetValue(str(self.Npts_total.GetValue()))
1455                self._save_plotting_range()
1456                self.create_default_data()
1457                self.state_change = True
1458                self._draw_model()
1459                # Time delay has been introduced to prevent _handle error
1460                # on Windows
1461                # This part of code is executed when model is selected and
1462                # it's parameters are changed (with respect to previously
1463                # selected model). There are two Iq evaluations occuring one
1464                # after another and therefore there may be compilation error
1465                # if model is calculated for the first time.
1466                # This seems to be Windows only issue - haven't tested on Linux
1467                # though.The proper solution (other than time delay) requires
1468                # more fundemental code refatoring
1469                # Wojtek P. Nov 7, 2016
1470                if not ON_MAC:
1471                    time.sleep(0.1)
1472                self.Refresh()
1473
1474        # logging.info("is_modified flag set to %g",is_modified)
1475        return is_modified
1476
1477    def _update_paramv_on_fit(self):
1478        """
1479        make sure that update param values just before the fitting
1480        """
1481        # flag for qmin qmax check values
1482        flag = True
1483        self.fitrange = True
1484        is_modified = False
1485
1486        # wx.PostEvent(self._manager.parent, StatusEvent(status=" \
1487        # updating ... ",type="update"))
1488
1489        # So make sure that update param values on_Fit.
1490        # self._undo.Enable(True)
1491        if self.model is not None:
1492            if self.Npts_total.GetValue() != self.Npts_fit.GetValue():
1493                if not self.data.is_data:
1494                    self._manager.page_finder[self.uid].set_fit_data(
1495                        data=[self.data])
1496            # Check the values
1497            is_modified = (self._check_value_enter(self.fittable_param)
1498                           or self._check_value_enter(self.fixed_param)
1499                           or self._check_value_enter(self.parameters))
1500
1501            # If qmin and qmax have been modified, update qmin and qmax and
1502            # Here we should check whether the boundaries have been modified.
1503            # If qmin and qmax have been modified, update qmin and qmax and
1504            # set the is_modified flag to True
1505            self.fitrange = self._validate_qrange(self.qmin, self.qmax)
1506            if self.fitrange:
1507                tempmin = float(self.qmin.GetValue())
1508                if tempmin != self.qmin_x:
1509                    self.qmin_x = tempmin
1510                tempmax = float(self.qmax.GetValue())
1511                if tempmax != self.qmax_x:
1512                    self.qmax_x = tempmax
1513                if tempmax == tempmin:
1514                    flag = False
1515                temp_smearer = None
1516                if not self.disable_smearer.GetValue():
1517                    temp_smearer = self.current_smearer
1518                    if self.slit_smearer.GetValue():
1519                        flag = self.update_slit_smear()
1520                    elif self.pinhole_smearer.GetValue():
1521                        flag = self.update_pinhole_smear()
1522                    else:
1523                        enable_smearer = not self.disable_smearer.GetValue()
1524                        self._manager.set_smearer(smearer=temp_smearer,
1525                                                  uid=self.uid,
1526                                                  fid=self.data.id,
1527                                                  qmin=float(self.qmin_x),
1528                                                  qmax=float(self.qmax_x),
1529                                                  enable_smearer=enable_smearer,
1530                                                  draw=False)
1531                elif not self._is_2D():
1532                    enable_smearer = not self.disable_smearer.GetValue()
1533                    self._manager.set_smearer(smearer=temp_smearer,
1534                                              qmin=float(self.qmin_x),
1535                                              uid=self.uid,
1536                                              fid=self.data.id,
1537                                              qmax=float(self.qmax_x),
1538                                              enable_smearer=enable_smearer,
1539                                              draw=False)
1540                    if self.data is not None:
1541                        index_data = ((self.qmin_x <= self.data.x) &
1542                                      (self.data.x <= self.qmax_x))
1543                        val = self.data.x[index_data is True]
1544                        val = len(val) if isinstance(val, list) else 1
1545                        self.Npts_fit.SetValue(str(val))
1546
1547                    else:
1548                        # No data in the panel
1549                        try:
1550                            self.npts_x = float(self.Npts_total.GetValue())
1551                        except:
1552                            flag = False
1553                            return flag
1554                    flag = True
1555                if self._is_2D():
1556                    # only 2D case set mask
1557                    flag = self._validate_Npts()
1558                    if not flag:
1559                        return flag
1560            else:
1561                flag = False
1562        else:
1563            flag = False
1564
1565        # For invalid q range, disable the mask editor and fit button, vs.
1566        if not self.fitrange:
1567            if self._is_2D():
1568                self.btEditMask.Disable()
1569        else:
1570            if self._is_2D() and self.data.is_data and not self.batch_on:
1571                self.btEditMask.Enable(True)
1572
1573        if not flag:
1574            msg = "Cannot Plot or Fit :Must select a "
1575            msg += " model or Fitting range is not valid!!!  "
1576            wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1577
1578        try:
1579            self.save_current_state()
1580        except Exception:
1581            logging.error(traceback.format_exc())
1582
1583        return flag, is_modified
1584
1585    def _reset_parameters_state(self, listtorestore, statelist):
1586        """
1587        Reset the parameters at the given state
1588        """
1589        if len(statelist) == 0 or len(listtorestore) == 0:
1590            return
1591        if len(statelist) != len(listtorestore):
1592            return
1593
1594        for j in range(len(listtorestore)):
1595            item_page = listtorestore[j]
1596            item_page_info = statelist[j]
1597            # change the state of the check box for simple parameters
1598            if item_page[0] is not None:
1599                item_page[0].SetValue(item_page_info[0])
1600            if item_page[2] is not None:
1601                item_page[2].SetValue(item_page_info[2])
1602                if item_page[2].__class__.__name__ == "ComboBox":
1603                    if item_page_info[2] in self.model.fun_list:
1604                        fun_val = self.model.fun_list[item_page_info[2]]
1605                        self.model.setParam(item_page_info[1], fun_val)
1606            if item_page[3] is not None:
1607                # show or hide text +/-
1608                if item_page_info[2]:
1609                    item_page[3].Show(True)
1610                else:
1611                    item_page[3].Hide()
1612            if item_page[4] is not None:
1613                # show of hide the text crtl for fitting error
1614                if item_page_info[4][0]:
1615                    item_page[4].Show(True)
1616                    item_page[4].SetValue(item_page_info[4][1])
1617                else:
1618                    item_page[3].Hide()
1619            if item_page[5] is not None:
1620                # show of hide the text crtl for fitting error
1621                item_page[5].Show(item_page_info[5][0])
1622                item_page[5].SetValue(item_page_info[5][1])
1623
1624            if item_page[6] is not None:
1625                # show of hide the text crtl for fitting error
1626                item_page[6].Show(item_page_info[6][0])
1627                item_page[6].SetValue(item_page_info[6][1])
1628
1629    def _reset_strparam_state(self, listtorestore, statelist):
1630        """
1631        Reset the string parameters at the given state
1632        """
1633        if len(statelist) == 0:
1634            return
1635
1636        listtorestore = copy.deepcopy(statelist)
1637
1638        for j in range(len(listtorestore)):
1639            item_page = listtorestore[j]
1640            item_page_info = statelist[j]
1641            # change the state of the check box for simple parameters
1642
1643            if item_page[0] is not None:
1644                item_page[0].SetValue(format_number(item_page_info[0], True))
1645
1646            if item_page[2] is not None:
1647                param_name = item_page_info[1]
1648                value = item_page_info[2]
1649                selection = value
1650                if value in self.model.fun_list:
1651                    selection = self.model.fun_list[value]
1652                item_page[2].SetValue(selection)
1653                self.model.setParam(param_name, selection)
1654
1655    def _copy_parameters_state(self, listtocopy, statelist):
1656        """
1657        copy the state of button
1658
1659        :param listtocopy: the list of check button to copy
1660        :param statelist: list of state object to store the current state
1661
1662        """
1663        if len(listtocopy) == 0:
1664            return
1665
1666        for item in listtocopy:
1667
1668            checkbox_state = None
1669            if item[0] is not None:
1670                checkbox_state = item[0].GetValue()
1671            parameter_name = item[1]
1672            parameter_value = None
1673            if item[2] is not None:
1674                parameter_value = item[2].GetValue()
1675            static_text = None
1676            if item[3] is not None:
1677                static_text = item[3].IsShown()
1678            error_value = None
1679            error_state = None
1680            if item[4] is not None:
1681                error_value = item[4].GetValue()
1682                error_state = item[4].IsShown()
1683
1684            min_value = None
1685            min_state = None
1686            if item[5] is not None:
1687                min_value = item[5].GetValue()
1688                min_state = item[5].IsShown()
1689
1690            max_value = None
1691            max_state = None
1692            if item[6] is not None:
1693                max_value = item[6].GetValue()
1694                max_state = item[6].IsShown()
1695            unit = None
1696            if item[7] is not None:
1697                unit = item[7].GetLabel()
1698
1699            statelist.append([checkbox_state, parameter_name, parameter_value,
1700                              static_text, [error_state, error_value],
1701                              [min_state, min_value],
1702                              [max_state, max_value], unit])
1703
1704    def _draw_model(self, update_chisqr=True, source='model'):
1705        """
1706        Method to draw or refresh a plotted model.
1707        The method will use the data member from the model page
1708        to build a call to the fitting perspective manager.
1709
1710        :param chisqr: update chisqr value [bool]
1711        """
1712        wx.CallAfter(self._draw_model_after, update_chisqr, source)
1713
1714    def _draw_model_after(self, update_chisqr=True, source='model'):
1715        """
1716        Method to draw or refresh a plotted model.
1717        The method will use the data member from the model page
1718        to build a call to the fitting perspective manager.
1719
1720        :param chisqr: update chisqr value [bool]
1721        """
1722        # if self.check_invalid_panel():
1723        #    return
1724        if self.model is not None:
1725            temp_smear = None
1726            if hasattr(self, "enable_smearer"):
1727                if not self.disable_smearer.GetValue():
1728                    temp_smear = self.current_smearer
1729            # compute weight for the current data
1730            from sas.sasgui.perspectives.fitting.utils import get_weight
1731            flag = self.get_weight_flag()
1732            weight = get_weight(data=self.data, is2d=self._is_2D(), flag=flag)
1733            toggle_mode_on = self.model_view.IsEnabled()
1734            is_2d = self._is_2D()
1735            self._manager.draw_model(self.model,
1736                                     data=self.data,
1737                                     smearer=temp_smear,
1738                                     qmin=float(self.qmin_x),
1739                                     qmax=float(self.qmax_x),
1740                                     page_id=self.uid,
1741                                     toggle_mode_on=toggle_mode_on,
1742                                     state=self.state,
1743                                     enable2D=is_2d,
1744                                     update_chisqr=update_chisqr,
1745                                     source='model',
1746                                     weight=weight)
1747
1748    def _on_show_sld(self, event=None):
1749        """
1750        Plot SLD profile
1751        """
1752        # get profile data
1753        x, y = self.model.getProfile()
1754
1755        from sas.sasgui.plottools import Data1D as pf_data1d
1756        # from sas.sasgui.perspectives.theory.profile_dialog import SLDPanel
1757        from sas.sasgui.guiframe.local_perspectives.plotting.profile_dialog \
1758            import SLDPanel
1759        sld_data = pf_data1d(x, y)
1760        sld_data.name = 'SLD'
1761        sld_data.axes = self.sld_axes
1762        self.panel = SLDPanel(self, data=sld_data, axes=self.sld_axes,
1763                              id=wx.ID_ANY)
1764        self.panel.ShowModal()
1765
1766    def _set_multfactor_combobox(self, multiplicity=10):
1767        """
1768        Set comboBox for muitfactor of CoreMultiShellModel
1769        :param multiplicit: no. of multi-functionality
1770        """
1771        # build content of the combobox
1772        for idx in range(0, multiplicity):
1773            self.multifactorbox.Append(str(idx), int(idx))
1774        self._hide_multfactor_combobox()
1775
1776    def _show_multfactor_combobox(self):
1777        """
1778        Show the comboBox of muitfactor of CoreMultiShellModel
1779        """
1780        if not self.mutifactor_text.IsShown():
1781            self.mutifactor_text.Show(True)
1782            self.mutifactor_text1.Show(True)
1783        if not self.multifactorbox.IsShown():
1784            self.multifactorbox.Show(True)
1785
1786    def _hide_multfactor_combobox(self):
1787        """
1788        Hide the comboBox of muitfactor of CoreMultiShellModel
1789        """
1790        if self.mutifactor_text.IsShown():
1791            self.mutifactor_text.Hide()
1792            self.mutifactor_text1.Hide()
1793        if self.multifactorbox.IsShown():
1794            self.multifactorbox.Hide()
1795
1796    def formfactor_combo_init(self):
1797        """
1798        First time calls _show_combox_helper
1799        """
1800        self._show_combox(None)
1801
1802    def _show_combox_helper(self):
1803        """
1804        Fill panel's combo box according to the type of model selected
1805        """
1806        custom_model = 'Customized Models'
1807        mod_cat = self.categorybox.GetStringSelection()
1808        self.structurebox.SetSelection(0)
1809        self.structurebox.Disable()
1810        self.formfactorbox.Clear()
1811        if mod_cat is None:
1812            return
1813        m_list = []
1814        try:
1815            if mod_cat == custom_model:
1816                for model in self.model_list_box[mod_cat]:
1817                    m_list.append(self.model_dict[model.name])
1818            else:
1819                cat_dic = self.master_category_dict[mod_cat]
1820                for (model, enabled) in cat_dic:
1821                    if enabled:
1822                        m_list.append(self.model_dict[model])
1823        except Exception:
1824            msg = traceback.format_exc()
1825            wx.PostEvent(self._manager.parent,
1826                         StatusEvent(status=msg, info="error"))
1827        self._populate_box(self.formfactorbox, m_list)
1828
1829    def _on_modify_cat(self, event=None):
1830        """
1831        Called when category manager is opened
1832        """
1833        self._manager.parent.on_category_panel(event)
1834
1835    def _show_combox(self, event=None):
1836        """
1837        Show combox box associate with type of model selected
1838        """
1839        self.Show(False)
1840        self._show_combox_helper()
1841        self._on_select_model(event=None)
1842        self.Show(True)
1843        self._save_typeOfmodel()
1844        self.sizer4_4.Layout()
1845        self.sizer4.Layout()
1846        self.Layout()
1847        self.Refresh()
1848
1849    def _populate_box(self, combobox, list):
1850        """
1851        fill combox box with dict item
1852
1853        :param list: contains item to fill the combox
1854            item must model class
1855        """
1856        mlist = []
1857        for models in list:
1858            if models.name != "NoStructure":
1859                mlist.append((models.name, models))
1860
1861        # Sort the models
1862        mlist_sorted = sorted(mlist)
1863        for item in mlist_sorted:
1864            combobox.Append(item[0], item[1])
1865        return 0
1866
1867    def _onQrangeEnter(self, event):
1868        """
1869        Check validity of value enter in the Q range field
1870
1871        """
1872        tcrtl = event.GetEventObject()
1873        # Clear msg if previously shown.
1874        msg = ""
1875        wx.PostEvent(self.parent, StatusEvent(status=msg))
1876        # Flag to register when a parameter has changed.
1877        if tcrtl.GetValue().lstrip().rstrip() != "":
1878            try:
1879                float(tcrtl.GetValue())
1880                tcrtl.SetBackgroundColour(wx.WHITE)
1881                # If qmin and qmax have been modified, update qmin and qmax
1882                if self._validate_qrange(self.qmin, self.qmax):
1883                    tempmin = float(self.qmin.GetValue())
1884                    if tempmin != self.qmin_x:
1885                        self.qmin_x = tempmin
1886                    tempmax = float(self.qmax.GetValue())
1887                    if tempmax != self.qmax_x:
1888                        self.qmax_x = tempmax
1889                else:
1890                    tcrtl.SetBackgroundColour("pink")
1891                    msg = "Model Error: wrong value entered: %s" % \
1892                          sys.exc_info()[1]
1893                    wx.PostEvent(self.parent, StatusEvent(status=msg))
1894                    return
1895            except:
1896                tcrtl.SetBackgroundColour("pink")
1897                msg = "Model Error: wrong value entered: %s" % sys.exc_info()[1]
1898                wx.PostEvent(self.parent, StatusEvent(status=msg))
1899                return
1900            # Check if # of points for theory model are valid(>0).
1901            if self.npts is not None:
1902                if check_float(self.npts):
1903                    temp_npts = float(self.npts.GetValue())
1904                    if temp_npts != self.num_points:
1905                        self.num_points = temp_npts
1906                else:
1907                    msg = "Cannot plot: No points in Q range!!!  "
1908                    wx.PostEvent(self.parent, StatusEvent(status=msg))
1909        else:
1910            tcrtl.SetBackgroundColour("pink")
1911            msg = "Model Error: wrong value entered!!!"
1912            wx.PostEvent(self.parent, StatusEvent(status=msg))
1913        self.save_current_state()
1914        event = PageInfoEvent(page=self)
1915        wx.PostEvent(self.parent, event)
1916        self.state_change = False
1917        # Draw the model for a different range
1918        if not self.data.is_data:
1919            self.create_default_data()
1920        self._draw_model()
1921
1922    def _theory_qrange_enter(self, event):
1923        """
1924        Check validity of value enter in the Q range field
1925        """
1926
1927        tcrtl = event.GetEventObject()
1928        # Clear msg if previously shown.
1929        msg = ""
1930        wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1931        # Flag to register when a parameter has changed.
1932        is_modified = False
1933        if tcrtl.GetValue().lstrip().rstrip() != "":
1934            try:
1935                value = float(tcrtl.GetValue())
1936                tcrtl.SetBackgroundColour(wx.WHITE)
1937
1938                # If qmin and qmax have been modified, update qmin and qmax
1939                if self._validate_qrange(self.theory_qmin, self.theory_qmax):
1940                    tempmin = float(self.theory_qmin.GetValue())
1941                    if tempmin != self.theory_qmin_x:
1942                        self.theory_qmin_x = tempmin
1943                    tempmax = float(self.theory_qmax.GetValue())
1944                    if tempmax != self.qmax_x:
1945                        self.theory_qmax_x = tempmax
1946                else:
1947                    tcrtl.SetBackgroundColour("pink")
1948                    msg = "Model Error: wrong value entered: %s" % \
1949                          sys.exc_info()[1]
1950                    wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1951                    return
1952            except:
1953                tcrtl.SetBackgroundColour("pink")
1954                msg = "Model Error: wrong value entered: %s" % sys.exc_info()[1]
1955                wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1956                return
1957            # Check if # of points for theory model are valid(>0).
1958            if self.Npts_total.IsEditable():
1959                if check_float(self.Npts_total):
1960                    temp_npts = float(self.Npts_total.GetValue())
1961                    if temp_npts != self.num_points:
1962                        self.num_points = temp_npts
1963                        is_modified = True
1964                else:
1965                    msg = "Cannot Plot: No points in Q range!!!  "
1966                    wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1967        else:
1968            tcrtl.SetBackgroundColour("pink")
1969            msg = "Model Error: wrong value entered!!!"
1970            wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1971        self.save_current_state()
1972        event = PageInfoEvent(page=self)
1973        wx.PostEvent(self.parent, event)
1974        self.state_change = False
1975        # Draw the model for a different range
1976        self.create_default_data()
1977        self._draw_model()
1978
1979    def _on_select_model_helper(self):
1980        """
1981        call back for model selection
1982        """
1983        # reset dictionary containing reference to dispersion
1984        self._disp_obj_dict = {}
1985        self.disp_cb_dict = {}
1986        self.temp_multi_functional = False
1987        f_id = self.formfactorbox.GetCurrentSelection()
1988        # For MAC
1989        form_factor = None
1990        if f_id >= 0:
1991            form_factor = self.formfactorbox.GetClientData(f_id)
1992
1993        if form_factor is None or \
1994            not hasattr(form_factor, 'is_form_factor') or \
1995                not form_factor.is_form_factor:
1996            self.structurebox.Hide()
1997            self.text2.Hide()
1998            self.structurebox.Disable()
1999            self.structurebox.SetSelection(0)
2000            self.text2.Disable()
2001        else:
2002            self.structurebox.Show()
2003            self.text2.Show()
2004            self.structurebox.Enable()
2005            self.text2.Enable()
2006
2007        if form_factor is not None:
2008            # set multifactor for Mutifunctional models
2009            if form_factor.is_multiplicity_model:
2010                m_id = self.multifactorbox.GetCurrentSelection()
2011                multiplicity = form_factor.multiplicity_info[0]
2012                self.multifactorbox.Clear()
2013                self._set_multfactor_combobox(multiplicity)
2014                self._show_multfactor_combobox()
2015                # ToDo: this info should be called directly from the model
2016                text = form_factor.multiplicity_info[1]  # 'No. of Shells: '
2017
2018                self.mutifactor_text.SetLabel(text)
2019                if m_id > multiplicity - 1:
2020                    # default value
2021                    m_id = 1
2022
2023                self.multi_factor = self.multifactorbox.GetClientData(m_id)
2024                if self.multi_factor is None:
2025                    self.multi_factor = 0
2026                self.multifactorbox.SetSelection(m_id)
2027                # Check len of the text1 and max_multiplicity
2028                text = ''
2029                if form_factor.multiplicity_info[0] == \
2030                        len(form_factor.multiplicity_info[2]):
2031                    text = form_factor.multiplicity_info[2][self.multi_factor]
2032                self.mutifactor_text1.SetLabel(text)
2033                # Check if model has  get sld profile.
2034                if len(form_factor.multiplicity_info[3]) > 0:
2035                    self.sld_axes = form_factor.multiplicity_info[3]
2036                    self.show_sld_button.Show(True)
2037                else:
2038                    self.sld_axes = ""
2039            else:
2040                self._hide_multfactor_combobox()
2041                self.show_sld_button.Hide()
2042                self.multi_factor = None
2043        else:
2044            self._hide_multfactor_combobox()
2045            self.show_sld_button.Hide()
2046            self.multi_factor = None
2047
2048        s_id = self.structurebox.GetCurrentSelection()
2049        struct_factor = self.structurebox.GetClientData(s_id)
2050
2051        if struct_factor is not None:
2052            from sasmodels.sasview_model import MultiplicationModel
2053            self.model = MultiplicationModel(form_factor(self.multi_factor),
2054                                             struct_factor())
2055            # multifunctional form factor
2056            if len(form_factor.non_fittable) > 0:
2057                self.temp_multi_functional = True
2058        elif form_factor is not None:
2059            if self.multi_factor is not None:
2060                self.model = form_factor(self.multi_factor)
2061            else:
2062                # old style plugin models do not accept a multiplicity argument
2063                self.model = form_factor()
2064        else:
2065            self.model = None
2066            return
2067
2068        # check if model has magnetic parameters
2069        if len(self.model.magnetic_params) > 0:
2070            self._has_magnetic = True
2071        else:
2072            self._has_magnetic = False
2073        # post state to fit panel
2074        self.state.parameters = []
2075        self.state.model = self.model
2076        self.state.qmin = self.qmin_x
2077        self.state.multi_factor = self.multi_factor
2078        self.disp_list = self.model.getDispParamList()
2079        self.state.disp_list = self.disp_list
2080        self.on_set_focus(None)
2081        self.Layout()
2082
2083    def _validate_qrange(self, qmin_ctrl, qmax_ctrl):
2084        """
2085        Verify that the Q range controls have valid values
2086        and that Qmin < Qmax.
2087
2088        :param qmin_ctrl: text control for Qmin
2089        :param qmax_ctrl: text control for Qmax
2090
2091        :return: True is the Q range is value, False otherwise
2092
2093        """
2094        qmin_validity = check_float(qmin_ctrl)
2095        qmax_validity = check_float(qmax_ctrl)
2096        if not (qmin_validity and qmax_validity):
2097            return False
2098        else:
2099            qmin = float(qmin_ctrl.GetValue())
2100            qmax = float(qmax_ctrl.GetValue())
2101            if qmin < qmax:
2102                # Make sure to set both colours white.
2103                qmin_ctrl.SetBackgroundColour(wx.WHITE)
2104                qmin_ctrl.Refresh()
2105                qmax_ctrl.SetBackgroundColour(wx.WHITE)
2106                qmax_ctrl.Refresh()
2107            else:
2108                qmin_ctrl.SetBackgroundColour("pink")
2109                qmin_ctrl.Refresh()
2110                qmax_ctrl.SetBackgroundColour("pink")
2111                qmax_ctrl.Refresh()
2112                msg = "Invalid Q range: Q min must be smaller than Q max"
2113                wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
2114                return False
2115        return True
2116
2117    def _validate_Npts(self):
2118        """
2119        Validate the number of points for fitting is more than 10 points.
2120        If valid, setvalues Npts_fit otherwise post msg.
2121        """
2122        # default flag
2123        flag = True
2124        # Theory
2125        if self.data is None and self.enable2D:
2126            return flag
2127        for data in self.data_list:
2128            # q value from qx and qy
2129            radius = numpy.sqrt(data.qx_data * data.qx_data +
2130                                data.qy_data * data.qy_data)
2131            # get unmasked index
2132            index_data = (float(self.qmin.GetValue()) <= radius) & \
2133                         (radius <= float(self.qmax.GetValue()))
2134            index_data = (index_data) & (data.mask)
2135            index_data = (index_data) & (numpy.isfinite(data.data))
2136
2137            if len(index_data[index_data]) < 10:
2138                # change the color pink.
2139                self.qmin.SetBackgroundColour("pink")
2140                self.qmin.Refresh()
2141                self.qmax.SetBackgroundColour("pink")
2142                self.qmax.Refresh()
2143                msg = "Data Error: "
2144                msg += "Too few points in %s." % data.name
2145                wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
2146                self.fitrange = False
2147                flag = False
2148            else:
2149                self.Npts_fit.SetValue(str(len(index_data[index_data is True])))
2150                self.fitrange = True
2151
2152        return flag
2153
2154    def _validate_Npts_1D(self):
2155        """
2156        Validate the number of points for fitting is more than 5 points.
2157        If valid, setvalues Npts_fit otherwise post msg.
2158        """
2159        # default flag
2160        flag = True
2161        # Theory
2162        if self.data is None:
2163            return flag
2164        for data in self.data_list:
2165            # q value from qx and qy
2166            radius = data.x
2167            # get unmasked index
2168            index_data = (float(self.qmin.GetValue()) <= radius) & \
2169                         (radius <= float(self.qmax.GetValue()))
2170            index_data = (index_data) & (numpy.isfinite(data.y))
2171
2172            if len(index_data[index_data]) < 5:
2173                # change the color pink.
2174                self.qmin.SetBackgroundColour("pink")
2175                self.qmin.Refresh()
2176                self.qmax.SetBackgroundColour("pink")
2177                self.qmax.Refresh()
2178                msg = "Data Error: "
2179                msg += "Too few points in %s." % data.name
2180                wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
2181                self.fitrange = False
2182                flag = False
2183            else:
2184                val = index_data[index_data is True]
2185                val = len(val) if isinstance(val, list) else 1
2186                self.Npts_fit.SetValue(str(val))
2187                self.fitrange = True
2188
2189        return flag
2190
2191    def _check_value_enter(self, list):
2192        """
2193        :param list: model parameter and panel info
2194        :Note: each item of the list should be as follow:
2195            item=[check button state, parameter's name,
2196                paramater's value, string="+/-",
2197                parameter's error of fit,
2198                parameter's minimum value,
2199                parameter's maximum value ,
2200                parameter's units]
2201
2202        Returns True if the model parameters have changed.
2203        """
2204        is_modified = False
2205        for item in list:
2206            # skip angle parameters for 1D
2207            if not self.enable2D and item in self.orientation_params:
2208                continue
2209
2210            value_ctrl = item[2]
2211            if not value_ctrl.IsEnabled():
2212                # ArrayDispersion disables PD, Min, Max, Npts, Nsigs
2213                continue
2214
2215            name = item[1]
2216            value_str = value_ctrl.GetValue().strip()
2217            if name.endswith(".npts"):
2218                validity = check_int(value_ctrl)
2219                if not validity:
2220                    continue
2221                value = int(value_str)
2222
2223            elif name.endswith(".nsigmas"):
2224                validity = check_float(value_ctrl)
2225                if not validity:
2226                    continue
2227                value = float(value_str)
2228
2229            else:  # value or polydispersity
2230
2231                # Check that min, max and value are floats
2232                min_ctrl, max_ctrl = item[5], item[6]
2233                min_str = min_ctrl.GetValue().strip()
2234                max_str = max_ctrl.GetValue().strip()
2235                validity = check_float(value_ctrl)
2236                if min_str != "":
2237                    validity = validity and check_float(min_ctrl)
2238                if max_str != "":
2239                    validity = validity and check_float(max_ctrl)
2240                if not validity:
2241                    continue
2242
2243                # Check that min is less than max
2244                low = -numpy.inf if min_str == "" else float(min_str)
2245                high = numpy.inf if max_str == "" else float(max_str)
2246                if high < low:
2247                    min_ctrl.SetBackgroundColour("pink")
2248                    min_ctrl.Refresh()
2249                    max_ctrl.SetBackgroundColour("pink")
2250                    max_ctrl.Refresh()
2251                    # msg = "Invalid fit range for %s: min must be smaller
2252                    # than max"%name
2253                    # wx.PostEvent(self._manager.parent,
2254                    # StatusEvent(status=msg))
2255                    continue
2256
2257                # Force value between min and max
2258                value = float(value_str)
2259                if value < low:
2260                    value = low
2261                    value_ctrl.SetValue(format_number(value))
2262                elif value > high:
2263                    value = high
2264                    value_ctrl.SetValue(format_number(value))
2265
2266                if name not in self.model.details.keys():
2267                    self.model.details[name] = ["", None, None]
2268                old_low, old_high = self.model.details[name][1:3]
2269                if old_low != low or old_high != high:
2270                    # The configuration has changed but it won't change the
2271                    # computed curve so no need to set is_modified to True
2272                    # is_modified = True
2273                    self.model.details[name][1:3] = low, high
2274
2275            # Update value in model if it has changed
2276            if value != self.model.getParam(name):
2277                self.model.setParam(name, value)
2278                is_modified = True
2279
2280        return is_modified
2281
2282    def _set_dipers_Param(self, event):
2283        """
2284        respond to self.enable_disp and self.disable_disp radio box.
2285        The dispersity object is reset inside the model into Gaussian.
2286        When the user select yes , this method display a combo box for
2287        more selection when the user selects No,the combo box disappears.
2288        Redraw the model with the default dispersity (Gaussian)
2289        """
2290        # On selction if no model exists.
2291        if self.model is None:
2292            self.disable_disp.SetValue(True)
2293            msg = "Please select a Model first..."
2294            wx.MessageBox(msg, 'Info')
2295            wx.PostEvent(self._manager.parent,
2296                         StatusEvent(status="Polydispersion: %s" % msg))
2297            return
2298
2299        self._reset_dispersity()
2300
2301        if self.model is None:
2302            self.model_disp.Hide()
2303            self.sizer4_4.Clear(True)
2304            return
2305
2306        if self.enable_disp.GetValue():
2307            # layout for model containing no dispersity parameters
2308
2309            self.disp_list = self.model.getDispParamList()
2310
2311            if len(self.disp_list) == 0 and len(self.disp_cb_dict) == 0:
2312                self._layout_sizer_noDipers()
2313            else:
2314                # set gaussian sizer
2315                self._on_select_Disp(event=None)
2316        else:
2317            self.sizer4_4.Clear(True)
2318
2319        # post state to fit panel
2320        self.save_current_state()
2321        if event is not None:
2322            event = PageInfoEvent(page=self)
2323            wx.PostEvent(self.parent, event)
2324        # draw the model with the current dispersity
2325
2326        # Wojtek P, Oct 8, 2016: Calling draw_model seems to be unessecary.
2327        # By comenting it we save an extra Iq calculation
2328        # self._draw_model()
2329
2330        # Need to use FitInside again here to replace the next four lines.
2331        # Otherwised polydispersity off does not resize the scrollwindow.
2332        # PDB Nov 28, 2015
2333        self.FitInside()
2334#        self.sizer4_4.Layout()
2335#        self.sizer5.Layout()
2336#        self.Layout()
2337#        self.Refresh()
2338
2339    def _layout_sizer_noDipers(self):
2340        """
2341        Draw a sizer with no dispersity info
2342        """
2343        ix = 0
2344        iy = 1
2345        self.fittable_param = []
2346        self.fixed_param = []
2347        self.orientation_params_disp = []
2348
2349        self.sizer4_4.Clear(True)
2350        text = "No polydispersity available for this model"
2351        model_disp = wx.StaticText(self, wx.ID_ANY, text)
2352        self.sizer4_4.Add(model_disp, (iy, ix), (1, 1),
2353                          wx.LEFT | wx.EXPAND | wx.ADJUST_MINSIZE, 10)
2354        self.sizer4_4.Layout()
2355        self.sizer4.Layout()
2356
2357    def _reset_dispersity(self):
2358        """
2359        put gaussian dispersity into current model
2360        """
2361        if len(self.param_toFit) > 0:
2362            for item in self.fittable_param:
2363                if item in self.param_toFit:
2364                    self.param_toFit.remove(item)
2365
2366            for item in self.orientation_params_disp:
2367                if item in self.param_toFit:
2368                    self.param_toFit.remove(item)
2369
2370        self.fittable_param = []
2371        self.fixed_param = []
2372        self.orientation_params_disp = []
2373        self.values = {}
2374        self.weights = {}
2375
2376        # from sas.models.dispersion_models import GaussianDispersion
2377        from sasmodels.weights import GaussianDispersion
2378        if len(self.disp_cb_dict) == 0:
2379            self.save_current_state()
2380            self.sizer4_4.Clear(True)
2381            self.Layout()
2382            return
2383        if (len(self.disp_cb_dict) > 0):
2384            for p in self.disp_cb_dict:
2385                # The parameter was un-selected.
2386                # Go back to Gaussian model (with 0 pts)
2387                disp_model = GaussianDispersion()
2388
2389                self._disp_obj_dict[p] = disp_model
2390                # Set the new model as the dispersion object
2391                # for the selected parameter
2392                try:
2393                    self.model.set_dispersion(p, disp_model)
2394                except Exception:
2395                    logging.error(traceback.format_exc())
2396
2397        # save state into
2398        self.save_current_state()
2399        self.Layout()
2400        self.Refresh()
2401
2402    def _on_select_Disp(self, event):
2403        """
2404        allow selecting different dispersion
2405        self.disp_list should change type later .now only gaussian
2406        """
2407        self._set_sizer_dispersion()
2408
2409        # Redraw the model
2410        #  Wojtek P. Nov 7, 2016: Redrawing seems to be unnecessary here
2411        # self._draw_model()
2412        # self._undo.Enable(True)
2413        event = PageInfoEvent(page=self)
2414        wx.PostEvent(self.parent, event)
2415
2416        self.sizer4_4.Layout()
2417        self.sizer4.Layout()
2418        self.SetupScrolling()
2419
2420    def _on_disp_func(self, event=None):
2421        """
2422        Select a distribution function for the polydispersion
2423
2424        :Param event: ComboBox event
2425        """
2426        # get ready for new event
2427        if event is not None:
2428            event.Skip()
2429        # Get event object
2430        disp_box = event.GetEventObject()
2431
2432        # Try to select a Distr. function
2433        try:
2434            disp_box.SetBackgroundColour("white")
2435            selection = disp_box.GetCurrentSelection()
2436            param_name = disp_box.Name.split('.')[0]
2437            disp_name = disp_box.GetValue()
2438            dispersity = disp_box.GetClientData(selection)
2439
2440            # disp_model =  GaussianDispersion()
2441            disp_model = dispersity()
2442            # Get param names to reset the values of the param
2443            name1 = param_name + ".width"
2444            name2 = param_name + ".npts"
2445            name3 = param_name + ".nsigmas"
2446            # Check Disp. function whether or not it is 'array'
2447            if disp_name.lower() == "array":
2448                value2 = ""
2449                value3 = ""
2450                value1 = self._set_array_disp(name=name1, disp=disp_model)
2451            else:
2452                self._del_array_values(name1)
2453                # self._reset_array_disp(param_name)
2454                self._disp_obj_dict[name1] = disp_model
2455                self.model.set_dispersion(param_name, disp_model)
2456                self.state._disp_obj_dict[name1] = disp_model.type
2457
2458                value1 = str(format_number(self.model.getParam(name1), True))
2459                value2 = str(format_number(self.model.getParam(name2)))
2460                value3 = str(format_number(self.model.getParam(name3)))
2461            # Reset fittable polydispersin parameter value
2462            for item in self.fittable_param:
2463                if item[1] == name1:
2464                    item[2].SetValue(value1)
2465                    item[5].SetValue("")
2466                    item[6].SetValue("")
2467                    # Disable for array
2468                    if disp_name.lower() == "array":
2469                        item[0].SetValue(False)
2470                        item[0].Disable()
2471                        item[2].Disable()
2472                        item[3].Show(False)
2473                        item[4].Show(False)
2474                        item[5].Disable()
2475                        item[6].Disable()
2476                    else:
2477                        item[0].Enable()
2478                        item[2].Enable()
2479                        item[3].Show(True)
2480                        item[4].Show(True)
2481                        item[5].Enable()
2482                        item[6].Enable()
2483                    break
2484            # Reset fixed polydispersion params
2485            for item in self.fixed_param:
2486                if item[1] == name2:
2487                    item[2].SetValue(value2)
2488                    # Disable Npts for array
2489                    if disp_name.lower() == "array":
2490                        item[2].Disable()
2491                    else:
2492                        item[2].Enable()
2493                if item[1] == name3:
2494                    item[2].SetValue(value3)
2495                    # Disable Nsigs for array
2496                    if disp_name.lower() == "array":
2497                        item[2].Disable()
2498                    else:
2499                        item[2].Enable()
2500
2501            # Make sure the check box updated
2502            self.get_all_checked_params()
2503
2504            # update params
2505            self._update_paramv_on_fit()
2506            # draw
2507            self._draw_model()
2508            self.Refresh()
2509        except Exception:
2510            logging.error(traceback.format_exc())
2511            # Error msg
2512            msg = "Error occurred:"
2513            msg += " Could not select the distribution function..."
2514            msg += " Please select another distribution function."
2515            disp_box.SetBackgroundColour("pink")
2516            # Focus on Fit button so that users can see the pinky box
2517            self.btFit.SetFocus()
2518            wx.PostEvent(self._manager.parent,
2519                         StatusEvent(status=msg, info="error"))
2520
2521    def _set_array_disp(self, name=None, disp=None):
2522        """
2523        Set array dispersion
2524
2525        :param name: name of the parameter for the dispersion to be set
2526        :param disp: the polydisperion object
2527        """
2528        # The user wants this parameter to be averaged.
2529        # Pop up the file selection dialog.
2530        path = self._selectDlg()
2531        # Array data
2532        values = []
2533        weights = []
2534        # If nothing was selected, just return
2535        if path is None:
2536            self.disp_cb_dict[name].SetValue(False)
2537            # self.noDisper_rbox.SetValue(True)
2538            return
2539        self._default_save_location = os.path.dirname(path)
2540        if self._manager is not None:
2541            self._manager.parent._default_save_location = \
2542                             self._default_save_location
2543
2544        basename = os.path.basename(path)
2545        values, weights = self.read_file(path)
2546
2547        # If any of the two arrays is empty, notify the user that we won't
2548        # proceed
2549        if len(self.param_toFit) > 0:
2550            if name in self.param_toFit:
2551                self.param_toFit.remove(name)
2552
2553        # Tell the user that we are about to apply the distribution
2554        msg = "Applying loaded %s distribution: %s" % (name, path)
2555        wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
2556        self._set_array_disp_model(name=name, disp=disp,
2557                                   values=values, weights=weights)
2558        return basename
2559
2560    def _set_array_disp_model(self, name=None, disp=None,
2561                              values=[], weights=[]):
2562        """
2563        Set array dispersion model
2564
2565        :param name: name of the parameter for the dispersion to be set
2566        :param disp: the polydisperion object
2567        """
2568        disp.set_weights(values, weights)
2569        self._disp_obj_dict[name] = disp
2570        self.model.set_dispersion(name.split('.')[0], disp)
2571        self.state._disp_obj_dict[name] = disp.type
2572        self.values[name] = values
2573        self.weights[name] = weights
2574        # Store the object to make it persist outside the
2575        # scope of this method
2576        # TODO: refactor model to clean this up?
2577        self.state.values = {}
2578        self.state.weights = {}
2579        self.state.values = copy.deepcopy(self.values)
2580        self.state.weights = copy.deepcopy(self.weights)
2581
2582        # Set the new model as the dispersion object for the
2583        # selected parameter
2584        # self.model.set_dispersion(p, disp_model)
2585        # Store a reference to the weights in the model object
2586        # so that
2587        # it's not lost when we use the model within another thread.
2588        self.state.model = self.model.clone()
2589        self.model._persistency_dict[name.split('.')[0]] = \
2590            [values, weights]
2591        self.state.model._persistency_dict[name.split('.')[0]] = \
2592            [values, weights]
2593
2594    def _del_array_values(self, name=None):
2595        """
2596        Reset array dispersion
2597
2598        :param name: name of the parameter for the dispersion to be set
2599        """
2600        # Try to delete values and weight of the names array dic if exists
2601        try:
2602            if name in self.values:
2603                del self.values[name]
2604                del self.weights[name]
2605                # delete all other dic
2606                del self.state.values[name]
2607                del self.state.weights[name]
2608                del self.model._persistency_dict[name.split('.')[0]]
2609                del self.state.model._persistency_dict[name.split('.')[0]]
2610        except Exception:
2611            logging.error(traceback.format_exc())
2612
2613    def _lay_out(self):
2614        """
2615        returns self.Layout
2616
2617        :Note: Mac seems to like this better when self.
2618            Layout is called after fitting.
2619        """
2620        self._sleep4sec()
2621        self.Layout()
2622        return
2623
2624    def _sleep4sec(self):
2625        """
2626            sleep for 1 sec only applied on Mac
2627            Note: This 1sec helps for Mac not to crash on self.
2628            Layout after self._draw_model
2629        """
2630        if ON_MAC is True:
2631            time.sleep(1)
2632
2633    def _find_polyfunc_selection(self, disp_func=None):
2634        """
2635        FInd Comboox selection from disp_func
2636
2637        :param disp_function: dispersion distr. function
2638        """
2639        # Find the selection
2640        if disp_func is not None:
2641            try:
2642                return POLYDISPERSITY_MODELS.values().index(disp_func.__class__)
2643            except ValueError:
2644                pass  # Fall through to default class
2645        return POLYDISPERSITY_MODELS.keys().index('gaussian')
2646
2647    def on_reset_clicked(self, event):
2648        """
2649        On 'Reset' button  for Q range clicked
2650        """
2651        flag = True
2652        # For 3 different cases: Data2D, Data1D, and theory
2653        if self.model is None:
2654            msg = "Please select a model first..."
2655            wx.MessageBox(msg, 'Info')
2656            flag = False
2657            return
2658
2659        elif self.data.__class__.__name__ == "Data2D":
2660            data_min = 0
2661            x = max(math.fabs(self.data.xmin), math.fabs(self.data.xmax))
2662            y = max(math.fabs(self.data.ymin), math.fabs(self.data.ymax))
2663            self.qmin_x = data_min
2664            self.qmax_x = math.sqrt(x * x + y * y)
2665            # self.data.mask = numpy.ones(len(self.data.data),dtype=bool)
2666            # check smearing
2667            if not self.disable_smearer.GetValue():
2668                # set smearing value whether or
2669                # not the data contain the smearing info
2670                if self.pinhole_smearer.GetValue():
2671                    flag = self.update_pinhole_smear()
2672                else:
2673                    flag = True
2674
2675        elif self.data is None:
2676            self.qmin_x = _QMIN_DEFAULT
2677            self.qmax_x = _QMAX_DEFAULT
2678            self.num_points = _NPTS_DEFAULT
2679            self.state.npts = self.num_points
2680
2681        elif self.data.__class__.__name__ != "Data2D":
2682            self.qmin_x = min(self.data.x)
2683            self.qmax_x = max(self.data.x)
2684            # check smearing
2685            if not self.disable_smearer.GetValue():
2686                # set smearing value whether or
2687                # not the data contain the smearing info
2688                if self.slit_smearer.GetValue():
2689                    flag = self.update_slit_smear()
2690                elif self.pinhole_smearer.GetValue():
2691                    flag = self.update_pinhole_smear()
2692                else:
2693                    flag = True
2694        else:
2695            flag = False
2696
2697        if flag is False:
2698            msg = "Cannot Plot :Must enter a number!!!  "
2699            wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
2700        else:
2701            # set relative text ctrs.
2702            self.qmin.SetValue(str(self.qmin_x))
2703            self.qmax.SetValue(str(self.qmax_x))
2704            self.show_npts2fit()
2705            # At this point, some button and variables satatus (disabled?)
2706            # should be checked such as color that should be reset to
2707            # white in case that it was pink.
2708            self._onparamEnter_helper()
2709
2710        self.save_current_state()
2711        self.state.qmin = self.qmin_x
2712        self.state.qmax = self.qmax_x
2713
2714        # reset the q range values
2715        self._reset_plotting_range(self.state)
2716        self._draw_model()
2717
2718    def select_log(self, event):
2719        """
2720        Log checked to generate log spaced points for theory model
2721        """
2722
2723    def get_images(self):
2724        """
2725        Get the images of the plots corresponding this panel for report
2726
2727        : return graphs: list of figures
2728        : Need Move to guiframe
2729        """
2730        # set list of graphs
2731        graphs = []
2732        canvases = []
2733        res_item = None
2734        # call gui_manager
2735        gui_manager = self._manager.parent
2736        # loops through the panels [dic]
2737        for _, item2 in gui_manager.plot_panels.iteritems():
2738            data_title = self.data.group_id
2739            # try to get all plots belonging to this control panel
2740            try:
2741                g_id = item2.group_id
2742                if g_id == data_title or \
2743                        str(g_id).count("res" + str(self.graph_id)) or \
2744                        str(g_id).count(str(self.uid)) > 0:
2745                    if str(g_id).count("res" + str(self.graph_id)) > 0:
2746                        res_item = [item2.figure, item2.canvas]
2747                    else:
2748                        # append to the list
2749                        graphs.append(item2.figure)
2750                        canvases.append(item2.canvas)
2751            except Exception:
2752                # Not for control panels
2753                logging.error(traceback.format_exc())
2754        # Make sure the resduals plot goes to the last
2755        if res_item is not None:
2756            graphs.append(res_item[0])
2757            canvases.append(res_item[1])
2758        # return the list of graphs
2759        return graphs, canvases
2760
2761    def on_function_help_clicked(self, event):
2762        """
2763        Function called when 'Help' button is pressed next to model
2764        of interest.  This calls DocumentationWindow from
2765        documentation_window.py. It will load the top level of the model
2766        help documenation sphinx generated html if no model is presented.
2767        If a model IS present then if documention for that model exists
2768        it will load to that  point otherwise again it will go to the top.
2769        For Wx2.8 and below is used (i.e. non-released through installer)
2770        a browser is loaded and the top of the model documentation only is
2771        accessible because webbrowser module does not pass anything after
2772        the # to the browser.
2773
2774        :param event: on Help Button pressed event
2775        """
2776
2777        if self.model is not None:
2778            name = self.formfactorbox.GetValue()
2779            _TreeLocation = 'user/models/' + name.lower()+'.html'
2780            _doc_viewer = DocumentationWindow(self, wx.ID_ANY, _TreeLocation,
2781                                              "", name + " Help")
2782        else:
2783            _TreeLocation = 'user/index.html'
2784            _doc_viewer = DocumentationWindow(self, wx.ID_ANY, _TreeLocation,
2785                                              "", "General Model Help")
2786
2787    def on_model_help_clicked(self, event):
2788        """
2789        Function called when 'Description' button is pressed next to model
2790        of interest.  This calls the Description embedded in the model. This
2791        should work with either Wx2.8 and lower or higher. If no model is
2792        selected it will give the message that a model must be chosen first
2793        in the box that would normally contain the description.  If a badly
2794        behaved model is encountered which has no description then it will
2795        give the message that none is available.
2796
2797        :param event: on Description Button pressed event
2798        """
2799
2800        if self.model is None:
2801            name = 'index.html'
2802        else:
2803            name = self.formfactorbox.GetValue()
2804
2805        msg = 'Model description:\n'
2806        info = "Info"
2807        if self.model is not None:
2808            # frame.Destroy()
2809            if str(self.model.description).rstrip().lstrip() == '':
2810                msg += "Sorry, no information is available for this model."
2811            else:
2812                msg += self.model.description + '\n'
2813            wx.MessageBox(msg, info)
2814        else:
2815            msg += "You must select a model to get information on this"
2816            wx.MessageBox(msg, info)
2817
2818    def _on_mag_angle_help(self, event):
2819        """
2820        Bring up Magnetic Angle definition bmp image whenever the ? button
2821        is clicked. Calls DocumentationWindow with the path of the location
2822        within the documentation tree (after /doc/ ....". When using old
2823        versions of Wx (i.e. before 2.9 and therefore not part of release
2824        versions distributed via installer) it brings up an image viewer
2825        box which allows the user to click through the rest of the images in
2826        the directory.  Not ideal but probably better than alternative which
2827        would bring up the entire discussion of how magnetic models work?
2828        Specially since it is not likely to be accessed.  The normal release
2829        versions bring up the normal image box.
2830
2831        :param evt: Triggers on clicking ? in Magnetic Angles? box
2832        """
2833
2834        _TreeLocation = "_images/M_angles_pic.bmp"
2835        _doc_viewer = DocumentationWindow(self, wx.ID_ANY, _TreeLocation, "",
2836                                          "Magnetic Angle Defintions")
2837
2838    def _on_mag_help(self, event):
2839        """
2840        Bring up Magnetic Angle definition bmp image whenever the ? button
2841        is clicked. Calls DocumentationWindow with the path of the location
2842        within the documentation tree (after /doc/ ....". When using old
2843        versions of Wx (i.e. before 2.9 and therefore not part of release
2844        versions distributed via installer) it brings up an image viewer
2845        box which allows the user to click through the rest of the images in
2846        the directory.  Not ideal but probably better than alternative which
2847        would bring up the entire discussion of how magnetic models work?
2848        Specially since it is not likely to be accessed.  The normal release
2849        versions bring up the normal image box.
2850
2851        :param evt: Triggers on clicking ? in Magnetic Angles? box
2852        """
2853
2854        _TreeLocation = "user/magnetism.html"
2855        _doc_viewer = DocumentationWindow(self, wx.ID_ANY, _TreeLocation, "",
2856                                          "Polarized Beam/Magnetc Help")
2857
2858    def _on_mag_on(self, event):
2859        """
2860        Magnetic Parameters ON/OFF
2861        """
2862        button = event.GetEventObject()
2863
2864        if button.GetLabel().count('ON') > 0:
2865            self.magnetic_on = True
2866            button.SetLabel("Magnetic OFF")
2867            m_value = 1.0e-06
2868            for key in self.model.magnetic_params:
2869                if key.count('M0') > 0:
2870                    self.model.setParam(key, m_value)
2871                    m_value += 0.5e-06
2872        else:
2873            self.magnetic_on = False
2874            button.SetLabel("Magnetic ON")
2875            for key in self.model.magnetic_params:
2876                if key.count('M0') > 0:
2877                    # reset mag value to zero fo safety
2878                    self.model.setParam(key, 0.0)
2879
2880        self.Show(False)
2881        self.set_model_param_sizer(self.model)
2882        # self._set_sizer_dispersion()
2883        self.state.magnetic_on = self.magnetic_on
2884        self.SetupScrolling()
2885        self.Show(True)
2886
2887    def on_pd_help_clicked(self, event):
2888        """
2889        Bring up Polydispersity Documentation whenever the ? button is clicked.
2890        Calls DocumentationWindow with the path of the location within the
2891        documentation tree (after /doc/ ....".  Note that when using old
2892        versions of Wx (before 2.9) and thus not the release version of
2893        istallers, the help comes up at the top level of the file as
2894        webbrowser does not pass anything past the # to the browser when it is
2895        running "file:///...."
2896
2897        :param event: Triggers on clicking ? in polydispersity box
2898        """
2899
2900        _TreeLocation = "user/sasgui/perspectives/fitting/pd_help.html"
2901        _PageAnchor = ""
2902        _doc_viewer = DocumentationWindow(self, wx.ID_ANY, _TreeLocation,
2903                                          _PageAnchor, "Polydispersity Help")
2904
2905    def on_left_down(self, event):
2906        """
2907        Get key stroke event
2908        """
2909        # Figuring out key combo: Cmd for copy, Alt for paste
2910        if event.CmdDown() and event.ShiftDown():
2911            self.get_paste()
2912        elif event.CmdDown():
2913            self.get_copy()
2914        else:
2915            event.Skip()
2916            return
2917        # make event free
2918        event.Skip()
2919
2920    def get_copy(self):
2921        """
2922        Get copy params to clipboard
2923        """
2924        content = self.get_copy_params()
2925        flag = self.set_clipboard(content)
2926        self._copy_info(flag)
2927        return flag
2928
2929    def get_copy_params(self):
2930        """
2931        Get the string copies of the param names and values in the tap
2932        """
2933        content = 'sasview_parameter_values:'
2934        # Do it if params exist
2935        if self.parameters:
2936
2937            # go through the parameters
2938            strings = self._get_copy_helper(self.parameters,
2939                                           self.orientation_params)
2940            content += strings
2941
2942            # go through the fittables
2943            strings = self._get_copy_helper(self.fittable_param,
2944                                           self.orientation_params_disp)
2945            content += strings
2946
2947            # go through the fixed params
2948            strings = self._get_copy_helper(self.fixed_param,
2949                                           self.orientation_params_disp)
2950            content += strings
2951
2952            # go through the str params
2953            strings = self._get_copy_helper(self.str_parameters,
2954                                           self.orientation_params)
2955            content += strings
2956            return content
2957        else:
2958            return False
2959
2960    def get_copy_excel(self):
2961        """
2962        Get copy params to clipboard
2963        """
2964        content = self.get_copy_params_excel()
2965        flag = self.set_clipboard(content)
2966        self._copy_info(flag)
2967        return flag
2968
2969    def get_copy_params_excel(self):
2970        """
2971        Get the string copies of the param names and values in the tap
2972        """
2973        content = ''
2974
2975        crlf = chr(13) + chr(10)
2976        tab = chr(9)
2977
2978        # Do it if params exist
2979        if self.parameters:
2980
2981            for param in self.parameters:
2982                content += param[1]  # parameter name
2983                content += tab
2984                content += param[1] + "_err"
2985                content += tab
2986
2987            content += crlf
2988
2989            # row of values and errors...
2990            for param in self.parameters:
2991                content += param[2].GetValue()  # value
2992                content += tab
2993                content += param[4].GetValue()  # error
2994                content += tab
2995
2996            return content
2997        else:
2998            return False
2999
3000    def get_copy_latex(self):
3001        """
3002        Get copy params to clipboard
3003        """
3004        content = self.get_copy_params_latex()
3005        flag = self.set_clipboard(content)
3006        self._copy_info(flag)
3007        return flag
3008
3009    def get_copy_params_latex(self):
3010        """
3011        Get the string copies of the param names and values in the tap
3012        """
3013        content = '\\begin{table}'
3014        content += '\\begin{tabular}[h]'
3015
3016        crlf = chr(13) + chr(10)
3017        tab = chr(9)
3018
3019        # Do it if params exist
3020        if self.parameters:
3021
3022            content += '{|'
3023            for param in self.parameters:
3024                content += 'l|l|'
3025            content += '}\hline'
3026            content += crlf
3027
3028            for index, param in enumerate(self.parameters):
3029                content += param[1].replace('_', '\_')  # parameter name
3030                content += ' & '
3031                content += param[1].replace('_', '\_') + "\_err"
3032                if index < len(self.parameters) - 1:
3033                    content += ' & '
3034            content += '\\\\ \\hline'
3035            content += crlf
3036
3037            # row of values and errors...
3038            for index, param in enumerate(self.parameters):
3039                content += param[2].GetValue()  # parameter value
3040                content += ' & '
3041                content += param[4].GetValue()  # parameter error
3042                if index < len(self.parameters) - 1:
3043                    content += ' & '
3044            content += '\\\\ \\hline'
3045            content += crlf
3046
3047            content += '\\end{tabular}'
3048            content += '\\end{table}'
3049            return content
3050        else:
3051            return False
3052
3053    def set_clipboard(self, content=None):
3054        """
3055        Put the string to the clipboard
3056        """
3057        if not content:
3058            return False
3059        if wx.TheClipboard.Open():
3060            wx.TheClipboard.SetData(wx.TextDataObject(str(content)))
3061            wx.TheClipboard.Close()
3062            return True
3063        return None
3064
3065    def _get_copy_helper(self, param, orient_param):
3066        """
3067        Helping get value and name of the params
3068
3069        : param param:  parameters
3070        : param orient_param = oritational params
3071        : return content: strings [list] [name,value:....]
3072        """
3073        content = ''
3074        bound_hi = ''
3075        bound_lo = ''
3076        # go through the str params
3077        for item in param:
3078            # copy only the params shown
3079            if not item[2].IsShown():
3080                continue
3081            disfunc = ''
3082            try:
3083                if item[7].__class__.__name__ == 'ComboBox':
3084                    disfunc = str(item[7].GetValue())
3085            except Exception:
3086                logging.error(traceback.format_exc())
3087
3088            # 2D
3089            if self.data.__class__.__name__ == "Data2D":
3090                try:
3091                    check = item[0].GetValue()
3092                except Exception:
3093                    check = None
3094                name = item[1]
3095                value = item[2].GetValue()
3096            # 1D
3097            else:
3098                # for 1D all parameters except orientation
3099                if not item[1] in orient_param:
3100                    try:
3101                        check = item[0].GetValue()
3102                    except:
3103                        check = None
3104                    name = item[1]
3105                    value = item[2].GetValue()
3106
3107            # Bounds
3108            try:
3109                bound_lo = item[5].GetValue()
3110                bound_hi = item[6].GetValue()
3111            except Exception:
3112                # harmless - need to just pass
3113                pass
3114
3115            # add to the content
3116            if disfunc != '':
3117
3118                disfunc = ',' + disfunc
3119            # Need to support array func for copy/paste
3120            try:
3121                if disfunc.count('array') > 0:
3122                    disfunc += ','
3123                    for val in self.values[name]:
3124                        disfunc += ' ' + str(val)
3125                    disfunc += ','
3126                    for weight in self.weights[name]:
3127                        disfunc += ' ' + str(weight)
3128            except Exception:
3129                logging.error(traceback.format_exc())
3130            content += name + ',' + str(check) + ',' + value + disfunc + ',' + \
3131                       bound_lo + ',' + bound_hi + ':'
3132
3133        return content
3134
3135    def get_clipboard(self):
3136        """
3137        Get strings in the clipboard
3138        """
3139        text = ""
3140        # Get text from the clip board
3141        if wx.TheClipboard.Open():
3142            if wx.TheClipboard.IsSupported(wx.DataFormat(wx.DF_TEXT)):
3143                data = wx.TextDataObject()
3144                # get wx dataobject
3145                success = wx.TheClipboard.GetData(data)
3146                # get text
3147                if success:
3148                    text = data.GetText()
3149                else:
3150                    text = ''
3151            # close clipboard
3152            wx.TheClipboard.Close()
3153        return text
3154
3155    def get_paste(self):
3156        """
3157        Paste params from the clipboard
3158        """
3159        text = self.get_clipboard()
3160        flag = self.get_paste_params(text)
3161        self._copy_info(flag)
3162        return flag
3163
3164    def get_paste_params(self, text=''):
3165        """
3166        Get the string copies of the param names and values in the tap
3167        """
3168        context = {}
3169        # put the text into dictionary
3170        lines = text.split(':')
3171        if lines[0] != 'sasview_parameter_values':
3172            self._copy_info(False)
3173            return False
3174        for line in lines[1:-1]:
3175            if len(line) != 0:
3176                item = line.split(',')
3177                check = item[1]
3178                name = item[0]
3179                value = item[2]
3180                # Transfer the text to content[dictionary]
3181                context[name] = [check, value]
3182
3183                # limits
3184                limit_lo = item[3]
3185                context[name].append(limit_lo)
3186                limit_hi = item[4]
3187                context[name].append(limit_hi)
3188
3189            # ToDo: PlugIn this poly disp function for pasting
3190            try:
3191                poly_func = item[5]
3192                context[name].append(poly_func)
3193                try:
3194                    # take the vals and weights for  array
3195                    array_values = item[6].split(' ')
3196                    array_weights = item[7].split(' ')
3197                    val = [float(a_val) for a_val in array_values[1:]]
3198                    weit = [float(a_weit) for a_weit in array_weights[1:]]
3199
3200                    context[name].append(val)
3201                    context[name].append(weit)
3202                except:
3203                    raise
3204            except:
3205                poly_func = ''
3206                context[name].append(poly_func)
3207
3208        # Do it if params exist
3209        if self.parameters:
3210            # go through the parameters
3211            self._get_paste_helper(self.parameters,
3212                                   self.orientation_params, context)
3213
3214            # go through the fittables
3215            self._get_paste_helper(self.fittable_param,
3216                                   self.orientation_params_disp,
3217                                   context)
3218
3219            # go through the fixed params
3220            self._get_paste_helper(self.fixed_param,
3221                                   self.orientation_params_disp, context)
3222
3223            # go through the str params
3224            self._get_paste_helper(self.str_parameters,
3225                                   self.orientation_params, context)
3226
3227            return True
3228        return None
3229
3230    def _get_paste_helper(self, param, orient_param, content):
3231        """
3232        Helping set values of the params
3233
3234        : param param:  parameters
3235        : param orient_param: oritational params
3236        : param content: dictionary [ name, value: name1.value1,...]
3237        """
3238        # go through the str params
3239        for item in param:
3240            # 2D
3241            if self.data.__class__.__name__ == "Data2D":
3242                name = item[1]
3243                if name in content.keys():
3244                    values = content[name]
3245                    check = values[0]
3246                    pd = values[1]
3247
3248                    if name.count('.') > 0:
3249                        # If this is parameter.width, then pd may be a floating
3250                        # point value or it may be an array distribution.
3251                        # Nothing to do for parameter.npts or parameter.nsigmas.
3252                        try:
3253                            float(pd)
3254                            if name.endswith('.npts'):
3255                                pd = int(pd)
3256                        except Exception:
3257                            # continue
3258                            if not pd and pd != '':
3259                                continue
3260                    item[2].SetValue(str(pd))
3261                    if item in self.fixed_param and pd == '':
3262                        # Only array func has pd == '' case.
3263                        item[2].Enable(False)
3264                    else:
3265                        item[2].Enable(True)
3266                    if item[2].__class__.__name__ == "ComboBox":
3267                        if content[name][1] in self.model.fun_list:
3268                            fun_val = self.model.fun_list[content[name][1]]
3269                            self.model.setParam(name, fun_val)
3270                    try:
3271                        item[5].SetValue(str(values[-3]))
3272                        item[6].SetValue(str(values[-2]))
3273                    except Exception:
3274                        # passing as harmless non-update
3275                        pass
3276
3277                    value = content[name][1:]
3278                    self._paste_poly_help(item, value)
3279                    if check == 'True':
3280                        is_true = True
3281                    elif check == 'False':
3282                        is_true = False
3283                    else:
3284                        is_true = None
3285                    if is_true is not None:
3286                        item[0].SetValue(is_true)
3287            # 1D
3288            else:
3289                # for 1D all parameters except orientation
3290                if not item[1] in orient_param:
3291                    name = item[1]
3292                    if name in content.keys():
3293                        check = content[name][0]
3294                        # Avoid changing combox content
3295                        value = content[name][1:]
3296                        pd = value[0]
3297                        if name.count('.') > 0:
3298                            # If this is parameter.width, then pd may be a
3299                            # floating point value or it may be an array
3300                            # distribution. Nothing to do for parameter.npts or
3301                            # parameter.nsigmas.
3302                            try:
3303                                pd = float(pd)
3304                                if name.endswith('.npts'):
3305                                    pd = int(pd)
3306                            except:
3307                                # continue
3308                                if not pd and pd != '':
3309                                    continue
3310                        item[2].SetValue(str(pd))
3311                        if item in self.fixed_param and pd == '':
3312                            # Only array func has pd == '' case.
3313                            item[2].Enable(False)
3314                        else:
3315                            item[2].Enable(True)
3316                        if item[2].__class__.__name__ == "ComboBox":
3317                            if value[0] in self.model.fun_list:
3318                                fun_val = self.model.fun_list[value[0]]
3319                                self.model.setParam(name, fun_val)
3320                                # save state
3321                        try:
3322                            item[5].SetValue(str(value[-3]))
3323                            item[6].SetValue(str(value[-2]))
3324                        except Exception:
3325                            # passing as harmless non-update
3326                            pass
3327
3328                        self._paste_poly_help(item, value)
3329                        if check == 'True':
3330                            is_true = True
3331                        elif check == 'False':
3332                            is_true = False
3333                        else:
3334                            is_true = None
3335                        if is_true is not None:
3336                            item[0].SetValue(is_true)
3337
3338        self.select_param(event=None)
3339        self.Refresh()
3340
3341    def _paste_poly_help(self, item, value):
3342        """
3343        Helps get paste for poly function
3344
3345        *item* is the parameter name
3346
3347        *value* depends on which parameter is being processed, and whether it
3348        has array polydispersity.
3349
3350        For parameters without array polydispersity:
3351
3352            parameter => ['FLOAT', '']
3353            parameter.width => ['FLOAT', 'DISTRIBUTION', '']
3354            parameter.npts => ['FLOAT', '']
3355            parameter.nsigmas => ['FLOAT', '']
3356
3357        For parameters with array polydispersity:
3358
3359            parameter => ['FLOAT', '']
3360            parameter.width => ['FILENAME', 'array', [x1, ...], [w1, ...]]
3361            parameter.npts => ['FLOAT', '']
3362            parameter.nsigmas => ['FLOAT', '']
3363        """
3364        # Do nothing if not setting polydispersity
3365        if len(value[3]) == 0:
3366            return
3367
3368        try:
3369            name = item[7].Name
3370            param_name = name.split('.')[0]
3371            item[7].SetValue(value[1])
3372            selection = item[7].GetCurrentSelection()
3373            dispersity = item[7].GetClientData(selection)
3374            disp_model = dispersity()
3375
3376            if value[1] == 'array':
3377                pd_vals = numpy.array(value[2])
3378                pd_weights = numpy.array(value[3])
3379                if len(pd_vals) == 0 or len(pd_vals) != len(pd_weights):
3380                    msg = ("bad array distribution parameters for %s"
3381                           % param_name)
3382                    raise ValueError(msg)
3383                self._set_disp_cb(True, item=item)
3384                self._set_array_disp_model(name=name,
3385                                           disp=disp_model,
3386                                           values=pd_vals,
3387                                           weights=pd_weights)
3388            else:
3389                self._set_disp_cb(False, item=item)
3390                self._disp_obj_dict[name] = disp_model
3391                self.model.set_dispersion(param_name, disp_model)
3392                self.state._disp_obj_dict[name] = disp_model.type
3393                # TODO: It's not an array, why update values and weights?
3394                self.model._persistency_dict[param_name] = \
3395                    [self.values, self.weights]
3396                self.state.values = self.values
3397                self.state.weights = self.weights
3398
3399        except Exception:
3400            logging.error(traceback.format_exc())
3401            print "Error in BasePage._paste_poly_help: %s" % \
3402                  sys.exc_info()[1]
3403
3404    def _set_disp_cb(self, isarray, item):
3405        """
3406        Set cb for array disp
3407        """
3408        if isarray:
3409            item[0].SetValue(False)
3410            item[0].Enable(False)
3411            item[2].Enable(False)
3412            item[3].Show(False)
3413            item[4].Show(False)
3414            item[5].SetValue('')
3415            item[5].Enable(False)
3416            item[6].SetValue('')
3417            item[6].Enable(False)
3418        else:
3419            item[0].Enable()
3420            item[2].Enable()
3421            item[3].Show(True)
3422            item[4].Show(True)
3423            item[5].Enable()
3424            item[6].Enable()
3425
3426    def update_pinhole_smear(self):
3427        """
3428            Method to be called by sub-classes
3429            Moveit; This method doesn't belong here
3430        """
3431        print "BasicPage.update_pinhole_smear was called: skipping"
3432        return
3433
3434    def _read_category_info(self):
3435        """
3436        Reads the categories in from file
3437        """
3438        # # ILL mod starts here - July 2012 kieranrcampbell@gmail.com
3439        self.master_category_dict = defaultdict(list)
3440        self.by_model_dict = defaultdict(list)
3441        self.model_enabled_dict = defaultdict(bool)
3442        categorization_file = CategoryInstaller.get_user_file()
3443        with open(categorization_file, 'rb') as f:
3444            self.master_category_dict = json.load(f)
3445        self._regenerate_model_dict()
3446
3447    def _regenerate_model_dict(self):
3448        """
3449        regenerates self.by_model_dict which has each model name as the
3450        key and the list of categories belonging to that model
3451        along with the enabled mapping
3452        """
3453        self.by_model_dict = defaultdict(list)
3454        for category in self.master_category_dict:
3455            for (model, enabled) in self.master_category_dict[category]:
3456                self.by_model_dict[model].append(category)
3457                self.model_enabled_dict[model] = enabled
3458
3459    def _populate_listbox(self):
3460        """
3461        fills out the category list box
3462        """
3463        uncat_str = 'Customized Models'
3464        self._read_category_info()
3465
3466        self.categorybox.Clear()
3467        cat_list = sorted(self.master_category_dict.keys())
3468        if uncat_str not in cat_list:
3469            cat_list.append(uncat_str)
3470
3471        for category in cat_list:
3472            if category != '':
3473                self.categorybox.Append(category)
3474
3475        if self.categorybox.GetSelection() == wx.NOT_FOUND:
3476            self.categorybox.SetSelection(0)
3477        else:
3478            self.categorybox.SetSelection(
3479                self.categorybox.GetSelection())
3480        # self._on_change_cat(None)
3481
3482    def _on_change_cat(self, event):
3483        """
3484        Callback for category change action
3485        """
3486        self.model_name = None
3487        category = self.categorybox.GetStringSelection()
3488        if category is None:
3489            return
3490        self.model_box.Clear()
3491
3492        if category == 'Customized Models':
3493            for model in self.model_list_box[category]:
3494                str_m = str(model).split(".")[0]
3495                self.model_box.Append(str_m)
3496
3497        else:
3498            for (model, enabled) in sorted(self.master_category_dict[category],
3499                                           key=lambda name: name[0]):
3500                if(enabled):
3501                    self.model_box.Append(model)
3502
3503    def _fill_model_sizer(self, sizer):
3504        """
3505        fill sizer containing model info
3506        """
3507        # This should only be called once per fit tab
3508        # print "==== Entering _fill_model_sizer"
3509        # Add model function Details button in fitpanel.
3510        # The following 3 lines are for Mac. Let JHC know before modifying...
3511        title = "Model"
3512        self.formfactorbox = None
3513        self.multifactorbox = None
3514        self.mbox_description = wx.StaticBox(self, wx.ID_ANY, str(title))
3515        boxsizer1 = wx.StaticBoxSizer(self.mbox_description, wx.VERTICAL)
3516        sizer_cat = wx.BoxSizer(wx.HORIZONTAL)
3517        self.mbox_description.SetForegroundColour(wx.RED)
3518        wx_id = self._ids.next()
3519        self.model_func = wx.Button(self, wx_id, 'Help', size=(80, 23))
3520        self.model_func.Bind(wx.EVT_BUTTON, self.on_function_help_clicked,
3521                             id=wx_id)
3522        self.model_func.SetToolTipString("Full Model Function Help")
3523        wx_id = self._ids.next()
3524        self.model_help = wx.Button(self, wx_id, 'Description', size=(80, 23))
3525        self.model_help.Bind(wx.EVT_BUTTON, self.on_model_help_clicked,
3526                             id=wx_id)
3527        self.model_help.SetToolTipString("Short Model Function Description")
3528        wx_id = self._ids.next()
3529        self.model_view = wx.Button(self, wx_id, "Show 2D", size=(80, 23))
3530        self.model_view.Bind(wx.EVT_BUTTON, self._onModel2D, id=wx_id)
3531        hint = "toggle view of model from 1D to 2D  or 2D to 1D"
3532        self.model_view.SetToolTipString(hint)
3533
3534        cat_set_box = wx.StaticBox(self, wx.ID_ANY, 'Category')
3535        sizer_cat_box = wx.StaticBoxSizer(cat_set_box, wx.HORIZONTAL)
3536        sizer_cat_box.SetMinSize((200, 50))
3537        self.categorybox = wx.ComboBox(self, wx.ID_ANY,
3538                                       style=wx.CB_READONLY)
3539        self.categorybox.SetToolTip(wx.ToolTip("Select a Category/Type"))
3540        self._populate_listbox()
3541        wx.EVT_COMBOBOX(self.categorybox, wx.ID_ANY, self._show_combox)
3542        # self.shape_rbutton = wx.RadioButton(self, wx.ID_ANY, 'Shapes',
3543        #                                     style=wx.RB_GROUP)
3544        # self.shape_indep_rbutton = wx.RadioButton(self, wx.ID_ANY,
3545        #                                          "Shape-Independent")
3546        # self.struct_rbutton = wx.RadioButton(self, wx.ID_ANY,
3547        #                                     "Structure Factor ")
3548        # self.plugin_rbutton = wx.RadioButton(self, wx.ID_ANY,
3549        #                                     "Uncategorized")
3550
3551        # self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
3552        #                   id=self.shape_rbutton.GetId())
3553        # self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
3554        #                    id=self.shape_indep_rbutton.GetId())
3555        # self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
3556        #                    id=self.struct_rbutton.GetId())
3557        # self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
3558        #                    id=self.plugin_rbutton.GetId())
3559        # MAC needs SetValue
3560
3561        show_cat_button = wx.Button(self, wx.ID_ANY, "Modify")
3562        cat_tip = "Modify model categories \n"
3563        cat_tip += "(also accessible from the menu bar)."
3564        show_cat_button.SetToolTip(wx.ToolTip(cat_tip))
3565        show_cat_button.Bind(wx.EVT_BUTTON, self._on_modify_cat)
3566        sizer_cat_box.Add(self.categorybox, 1, wx.RIGHT, 3)
3567        sizer_cat_box.Add((10, 10))
3568        sizer_cat_box.Add(show_cat_button)
3569        # self.shape_rbutton.SetValue(True)
3570
3571        sizer_radiobutton = wx.GridSizer(2, 2, 5, 5)
3572        # sizer_radiobutton.Add(self.shape_rbutton)
3573        # sizer_radiobutton.Add(self.shape_indep_rbutton)
3574        sizer_radiobutton.Add((5, 5))
3575        sizer_radiobutton.Add(self.model_view, 1, wx.RIGHT, 5)
3576        # sizer_radiobutton.Add(self.plugin_rbutton)
3577        # sizer_radiobutton.Add(self.struct_rbutton)
3578        # sizer_radiobutton.Add((5,5))
3579        sizer_radiobutton.Add(self.model_help, 1, wx.RIGHT | wx.LEFT, 5)
3580        # sizer_radiobutton.Add((5,5))
3581        sizer_radiobutton.Add(self.model_func, 1, wx.RIGHT, 5)
3582        sizer_cat.Add(sizer_cat_box, 1, wx.LEFT, 2.5)
3583        sizer_cat.Add(sizer_radiobutton)
3584        sizer_selection = wx.BoxSizer(wx.HORIZONTAL)
3585        mutifactor_selection = wx.BoxSizer(wx.HORIZONTAL)
3586
3587        self.text1 = wx.StaticText(self, wx.ID_ANY, "")
3588        self.text2 = wx.StaticText(self, wx.ID_ANY, "P(Q)*S(Q)")
3589        self.mutifactor_text = wx.StaticText(self, wx.ID_ANY, "No. of Shells: ")
3590        self.mutifactor_text1 = wx.StaticText(self, wx.ID_ANY, "")
3591        self.show_sld_button = wx.Button(self, wx.ID_ANY, "Show SLD Profile")
3592        self.show_sld_button.Bind(wx.EVT_BUTTON, self._on_show_sld)
3593
3594        self.formfactorbox = wx.ComboBox(self, wx.ID_ANY, style=wx.CB_READONLY)
3595        self.formfactorbox.SetToolTip(wx.ToolTip("Select a Model"))
3596        if self.model is not None:
3597            self.formfactorbox.SetValue(self.model.name)
3598        self.structurebox = wx.ComboBox(self, wx.ID_ANY, style=wx.CB_READONLY)
3599        self.multifactorbox = wx.ComboBox(self, wx.ID_ANY, style=wx.CB_READONLY)
3600        self.initialize_combox()
3601        wx.EVT_COMBOBOX(self.formfactorbox, wx.ID_ANY, self._on_select_model)
3602
3603        wx.EVT_COMBOBOX(self.structurebox, wx.ID_ANY, self._on_select_model)
3604        wx.EVT_COMBOBOX(self.multifactorbox, wx.ID_ANY, self._on_select_model)
3605        # check model type to show sizer
3606        if self.model is not None:
3607            print "_set_model_sizer_selection: disabled."
3608            # self._set_model_sizer_selection(self.model)
3609
3610        sizer_selection.Add(self.text1)
3611        sizer_selection.Add((10, 5))
3612        sizer_selection.Add(self.formfactorbox)
3613        sizer_selection.Add((5, 5))
3614        sizer_selection.Add(self.text2)
3615        sizer_selection.Add((5, 5))
3616        sizer_selection.Add(self.structurebox)
3617
3618        mutifactor_selection.Add((13, 5))
3619        mutifactor_selection.Add(self.mutifactor_text)
3620        mutifactor_selection.Add(self.multifactorbox)
3621        mutifactor_selection.Add((5, 5))
3622        mutifactor_selection.Add(self.mutifactor_text1)
3623        mutifactor_selection.Add((10, 5))
3624        mutifactor_selection.Add(self.show_sld_button)
3625
3626        boxsizer1.Add(sizer_cat)
3627        boxsizer1.Add((10, 10))
3628        boxsizer1.Add(sizer_selection)
3629        boxsizer1.Add((10, 10))
3630        boxsizer1.Add(mutifactor_selection)
3631
3632        self._set_multfactor_combobox()
3633        self.multifactorbox.SetSelection(1)
3634        self.show_sld_button.Hide()
3635        sizer.Add(boxsizer1, 0, wx.EXPAND | wx.ALL, 10)
3636        sizer.Layout()
3637
3638    def on_smear_helper(self, update=False):
3639        """
3640        Help for onSmear if implemented
3641
3642        :param update: force or not to update
3643        """
3644    def reset_page(self, state, first=False):
3645        """
3646        reset the state  if implemented
3647        """
3648    def onSmear(self, event):
3649        """
3650        Create a smear object if implemented
3651        """
3652    def onPinholeSmear(self, event):
3653        """
3654        Create a custom pinhole smear object if implemented
3655        """
3656    def onSlitSmear(self, event):
3657        """
3658        Create a custom slit smear object if implemented
3659        """
3660    def update_slit_smear(self):
3661        """
3662        called by kill_focus on pinhole TextCntrl
3663        to update the changes if implemented
3664        """
3665    def select_param(self, event):
3666        """
3667        Select TextCtrl  checked if implemented
3668        """
3669    def set_data(self, data=None):
3670        """
3671        Sets data if implemented
3672        """
3673    def _is_2D(self):
3674        """
3675        Check if data_name is Data2D if implemented
3676        """
3677    def _on_select_model(self, event=None):
3678        """
3679        call back for model selection if implemented
3680        """
3681    def get_weight_flag(self):
3682        """
3683        Get flag corresponding to a given weighting dI data if implemented
3684        """
3685    def _set_sizer_dispersion(self):
3686        """
3687        draw sizer for dispersity if implemented
3688        """
3689    def get_all_checked_params(self):
3690        """
3691        Found all parameters current check and add them to list of parameters
3692        to fit if implemented
3693        """
3694    def show_npts2fit(self):
3695        """
3696        setValue Npts for fitting if implemented
3697        """
3698    def _onModel2D(self, event):
3699        """
3700        toggle view of model from 1D to 2D  or 2D from 1D if implemented
3701        """
3702
3703
3704class ModelTextCtrl(wx.TextCtrl):
3705    """
3706    Text control for model and fit parameters.
3707    Binds the appropriate events for user interactions.
3708    Default callback methods can be overwritten on initialization
3709
3710    :param kill_focus_callback: callback method for EVT_KILL_FOCUS event
3711    :param set_focus_callback:  callback method for EVT_SET_FOCUS event
3712    :param mouse_up_callback:   callback method for EVT_LEFT_UP event
3713    :param text_enter_callback: callback method for EVT_TEXT_ENTER event
3714
3715    """
3716    # Set to True when the mouse is clicked while whole string is selected
3717    full_selection = False
3718    # Call back for EVT_SET_FOCUS events
3719    _on_set_focus_callback = None
3720
3721    def __init__(self, parent, id=-1,
3722                 value=wx.EmptyString,
3723                 pos=wx.DefaultPosition,
3724                 size=wx.DefaultSize,
3725                 style=0,
3726                 validator=wx.DefaultValidator,
3727                 name=wx.TextCtrlNameStr,
3728                 kill_focus_callback=None,
3729                 set_focus_callback=None,
3730                 mouse_up_callback=None,
3731                 text_enter_callback=None):
3732
3733        wx.TextCtrl.__init__(self, parent, id, value, pos,
3734                             size, style, validator, name)
3735
3736        # Bind appropriate events
3737        self._on_set_focus_callback = parent.onSetFocus \
3738            if set_focus_callback is None else set_focus_callback
3739        self.Bind(wx.EVT_SET_FOCUS, self._on_set_focus)
3740        self.Bind(wx.EVT_KILL_FOCUS, self._silent_kill_focus
3741        if kill_focus_callback is None else kill_focus_callback)
3742        self.Bind(wx.EVT_TEXT_ENTER, parent._onparamEnter
3743        if text_enter_callback is None else text_enter_callback)
3744        if not ON_MAC:
3745            self.Bind(wx.EVT_LEFT_UP, self._highlight_text
3746            if mouse_up_callback is None else mouse_up_callback)
3747
3748    def _on_set_focus(self, event):
3749        """
3750        Catch when the text control is set in focus to highlight the whole
3751        text if necessary
3752
3753        :param event: mouse event
3754
3755        """
3756        event.Skip()
3757        self.full_selection = True
3758        return self._on_set_focus_callback(event)
3759
3760    def _highlight_text(self, event):
3761        """
3762        Highlight text of a TextCtrl only of no text has be selected
3763
3764        :param event: mouse event
3765
3766        """
3767        # Make sure the mouse event is available to other listeners
3768        event.Skip()
3769        control = event.GetEventObject()
3770        if self.full_selection:
3771            self.full_selection = False
3772            # Check that we have a TextCtrl
3773            if issubclass(control.__class__, wx.TextCtrl):
3774                # Check whether text has been selected,
3775                # if not, select the whole string
3776                (start, end) = control.GetSelection()
3777                if start == end:
3778                    control.SetSelection(-1, -1)
3779
3780    def _silent_kill_focus(self, event):
3781        """
3782        Save the state of the page
3783        """
3784
3785        event.Skip()
3786        # pass
Note: See TracBrowser for help on using the repository browser.