source: sasview/src/sas/sasgui/perspectives/fitting/basepage.py @ 367a31c

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalcmagnetic_scattrelease-4.2.2ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since 367a31c was 367a31c, checked in by Paul Kienzle <pkienzle@…>, 7 years ago

Merge branch 'master' into ticket-869

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