source: sasview/src/sas/sasgui/perspectives/fitting/basepage.py @ 922ef55

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 922ef55 was 922ef55, checked in by Tim Snow <tim.snow@…>, 7 years ago

Quick stash commit

Changing branches

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