basepage.py @ 883e5f5

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload

Last change on this file since 883e5f5 was 6bbeacd4, checked in by Gervaise Alina <gervyh@…>, 14 years ago
remove other type of data
Property mode set to `100644`
File size: 103.2 KB

Line
1
2	import sys
3	import os
4	import wx
5	import numpy
6	import time
7	import copy
8	import math
9	import string
10	from sans.guiframe.panel_base import PanelBase
11	from wx.lib.scrolledpanel import ScrolledPanel
12	from sans.guiframe.utils import format_number,check_float
13	from sans.guiframe.events import PanelOnFocusEvent
14	from sans.guiframe.events import StatusEvent
15	from sans.guiframe.events import AppendBookmarkEvent
16	import pagestate
17	from pagestate import PageState
18	(PageInfoEvent, EVT_PAGE_INFO) = wx.lib.newevent.NewEvent()
19	(PreviousStateEvent, EVT_PREVIOUS_STATE) = wx.lib.newevent.NewEvent()
20	(NextStateEvent, EVT_NEXT_STATE) = wx.lib.newevent.NewEvent()
21
22	_BOX_WIDTH = 76
23	_QMIN_DEFAULT = 0.0005
24	_QMAX_DEFAULT = 0.5
25	_NPTS_DEFAULT = 50
26	#Control panel width
27	if sys.platform.count("darwin")==0:
28	PANEL_WIDTH = 450
29	FONT_VARIANT = 0
30	ON_MAC = False
31	else:
32	PANEL_WIDTH = 500
33	FONT_VARIANT = 1
34	ON_MAC = True
35
36
37
38	class BasicPage(ScrolledPanel, PanelBase):
39	"""
40	This class provide general structure of fitpanel page
41	"""
42	## Internal name for the AUI manager
43	window_name = "Fit Page"
44	## Title to appear on top of the window
45	window_caption = "Fit Page "
46
47	def __init__(self, parent,color='blue', **kwargs):
48	"""
49	"""
50	ScrolledPanel.__init__(self, parent, **kwargs)
51	PanelBase.__init__(self)
52	self.SetupScrolling()
53	#Set window's font size
54	self.SetWindowVariant(variant=FONT_VARIANT)
55
56	self.SetBackgroundColour(color)
57	## parent of the page
58	self.parent = parent
59	## manager is the fitting plugin
60	## owner of the page (fitting plugin)
61	self.event_owner = None
62	## current model
63	self.model = None
64	## data
65	self.data = None
66	self.mask = None
67	self.id = None
68	## Q range
69	self.qmax_x = _QMAX_DEFAULT
70	self.qmin_x = _QMIN_DEFAULT
71	self.npts_x = _NPTS_DEFAULT
72	## total number of point: float
73	self.npts = None
74	## default fitengine type
75	self.engine_type = 'scipy'
76	## smear default
77	self.smearer = None
78	self.current_smearer = None
79	## 2D smear accuracy default
80	self.smear2d_accuracy = 'Low'
81	## slit smear:
82	self.dxl = None
83	self.dxw = None
84	## pinhole smear
85	self.dx_min = None
86	self.dx_max = None
87
88	self.disp_cb_dict = {}
89
90	self.state = PageState(parent=parent)
91	## dictionary containing list of models
92	self.model_list_box = {}
93
94	## Data member to store the dispersion object created
95	self._disp_obj_dict = {}
96	## selected parameters to apply dispersion
97	self.disp_cb_dict ={}
98
99	## smearer object
100	self.smearer = None
101
102	##list of model parameters. each item must have same length
103	## each item related to a given parameters
104	##[cb state, name, value, "+/-", error of fit, min, max , units]
105	self.parameters = []
106	# non-fittable parameter whose value is astring
107	self.str_parameters = []
108	## list of parameters to fit , must be like self.parameters
109	self.param_toFit = []
110	## list of looking like parameters but with non fittable parameters info
111	self.fixed_param = []
112	## list of looking like parameters but with fittable parameters info
113	self.fittable_param = []
114	##list of dispersion parameters
115	self.disp_list = []
116	self.disp_name = ""
117
118	## list of orientation parameters
119	self.orientation_params = []
120	self.orientation_params_disp = []
121	if self.model != None:
122	self.disp_list = self.model.getDispParamList()
123
124	##enable model 2D draw
125	self.enable2D = False
126	## check that the fit range is correct to plot the model again
127	self.fitrange = True
128	## Create memento to save the current state
129	self.state = PageState(parent=self.parent,
130	model=self.model, data=self.data)
131	## flag to determine if state has change
132	self.state_change = False
133	## save customized array
134	self.values = []
135	self.weights = []
136	## retrieve saved state
137	self.number_saved_state = 0
138	## dictionary of saved state
139	self.saved_states = {}
140	## Create context menu for page
141	self.popUpMenu = wx.Menu()
142
143	id = wx.NewId()
144	self._keep = wx.MenuItem(self.popUpMenu,id,"BookMark",
145	" Keep the panel status to recall it later")
146	self.popUpMenu.AppendItem(self._keep)
147	self._keep.Enable(True)
148	self._set_bookmark_flag(False)
149	self._set_save_flag(False)
150	wx.EVT_MENU(self, id, self.on_bookmark)
151	self.popUpMenu.AppendSeparator()
152
153	## Default locations
154	self._default_save_location = os.getcwd()
155	## save initial state on context menu
156	#self.onSave(event=None)
157	self.Bind(wx.EVT_CONTEXT_MENU, self.onContextMenu)
158
159	## create the basic structure of the panel with empty sizer
160	self.define_page_structure()
161	## drawing Initial dispersion parameters sizer
162	self.set_dispers_sizer()
163
164	## layout
165	self.set_layout()
166
167
168
169	def on_set_focus(self, event):
170	"""
171	"""
172	if self._manager is not None:
173	wx.PostEvent(self._manager.parent, PanelOnFocusEvent(panel=self))
174
175	class ModelTextCtrl(wx.TextCtrl):
176	"""
177	Text control for model and fit parameters.
178	Binds the appropriate events for user interactions.
179	Default callback methods can be overwritten on initialization
180
181	:param kill_focus_callback: callback method for EVT_KILL_FOCUS event
182	:param set_focus_callback: callback method for EVT_SET_FOCUS event
183	:param mouse_up_callback: callback method for EVT_LEFT_UP event
184	:param text_enter_callback: callback method for EVT_TEXT_ENTER event
185
186	"""
187	## Set to True when the mouse is clicked while the whole string is selected
188	full_selection = False
189	## Call back for EVT_SET_FOCUS events
190	_on_set_focus_callback = None
191
192	def __init__(self, parent, id=-1,
193	value=wx.EmptyString,
194	pos=wx.DefaultPosition,
195	size=wx.DefaultSize,
196	style=0,
197	validator=wx.DefaultValidator,
198	name=wx.TextCtrlNameStr,
199	kill_focus_callback=None,
200	set_focus_callback=None,
201	mouse_up_callback=None,
202	text_enter_callback = None):
203
204	wx.TextCtrl.__init__(self, parent, id, value, pos,
205	size, style, validator, name)
206
207	# Bind appropriate events
208	self._on_set_focus_callback = parent.onSetFocus \
209	if set_focus_callback is None else set_focus_callback
210	self.Bind(wx.EVT_SET_FOCUS, self._on_set_focus)
211	self.Bind(wx.EVT_KILL_FOCUS, self._silent_kill_focus \
212	if kill_focus_callback is None else kill_focus_callback)
213	self.Bind(wx.EVT_TEXT_ENTER, parent._onparamEnter \
214	if text_enter_callback is None else text_enter_callback)
215	if not ON_MAC :
216	self.Bind(wx.EVT_LEFT_UP, self._highlight_text \
217	if mouse_up_callback is None else mouse_up_callback)
218
219	def _on_set_focus(self, event):
220	"""
221	Catch when the text control is set in focus to highlight the whole
222	text if necessary
223
224	:param event: mouse event
225
226	"""
227	event.Skip()
228	self.full_selection = True
229	return self._on_set_focus_callback(event)
230
231
232
233	def _highlight_text(self, event):
234	"""
235	Highlight text of a TextCtrl only of no text has be selected
236
237	:param event: mouse event
238
239	"""
240	# Make sure the mouse event is available to other listeners
241	event.Skip()
242	control = event.GetEventObject()
243	if self.full_selection:
244	self.full_selection = False
245	# Check that we have a TextCtrl
246	if issubclass(control.__class__, wx.TextCtrl):
247	# Check whether text has been selected,
248	# if not, select the whole string
249	(start, end) = control.GetSelection()
250	if start==end:
251	control.SetSelection(-1,-1)
252
253	def _silent_kill_focus(self,event):
254	"""
255	Save the state of the page
256	"""
257
258	event.Skip()
259	pass
260
261	def set_page_info(self, page_info):
262	"""
263	set some page important information at once
264	"""
265	##window_name
266	self.window_name = page_info.window_name
267	##window_caption
268	self.window_caption = page_info.window_caption
269	## manager is the fitting plugin
270	self._manager= page_info.manager
271	## owner of the page (fitting plugin)
272	self.event_owner= page_info.event_owner
273	## current model
274	self.model = page_info.model
275	## data
276	self.data = page_info.data
277	## dictionary containing list of models
278	self.model_list_box = page_info.model_list_box
279	## Data member to store the dispersion object created
280	self.populate_box(dict=self.model_list_box)
281
282	def onContextMenu(self, event):
283	"""
284	Retrieve the state selected state
285	"""
286	# Skipping the save state functionality for release 0.9.0
287	#return
288
289	pos = event.GetPosition()
290	pos = self.ScreenToClient(pos)
291
292	self.PopupMenu(self.popUpMenu, pos)
293
294
295	def onUndo(self, event):
296	"""
297	Cancel the previous action
298	"""
299	event = PreviousStateEvent(page = self)
300	wx.PostEvent(self.parent, event)
301
302	def onRedo(self, event):
303	"""
304	Restore the previous action cancelled
305	"""
306	event = NextStateEvent(page= self)
307	wx.PostEvent(self.parent, event)
308
309	def define_page_structure(self):
310	"""
311	Create empty sizer for a panel
312	"""
313	self.vbox = wx.BoxSizer(wx.VERTICAL)
314	self.sizer0 = wx.BoxSizer(wx.VERTICAL)
315	self.sizer1 = wx.BoxSizer(wx.VERTICAL)
316	self.sizer2 = wx.BoxSizer(wx.VERTICAL)
317	self.sizer3 = wx.BoxSizer(wx.VERTICAL)
318	self.sizer4 = wx.BoxSizer(wx.VERTICAL)
319	self.sizer5 = wx.BoxSizer(wx.VERTICAL)
320	self.sizer6 = wx.BoxSizer(wx.VERTICAL)
321
322	self.sizer0.SetMinSize((PANEL_WIDTH,-1))
323	self.sizer1.SetMinSize((PANEL_WIDTH,-1))
324	self.sizer2.SetMinSize((PANEL_WIDTH,-1))
325	self.sizer3.SetMinSize((PANEL_WIDTH,-1))
326	self.sizer4.SetMinSize((PANEL_WIDTH,-1))
327	self.sizer5.SetMinSize((PANEL_WIDTH,-1))
328	self.sizer6.SetMinSize((PANEL_WIDTH,-1))
329
330	self.vbox.Add(self.sizer0)
331	self.vbox.Add(self.sizer1)
332	self.vbox.Add(self.sizer2)
333	self.vbox.Add(self.sizer3)
334	self.vbox.Add(self.sizer4)
335	self.vbox.Add(self.sizer5)
336	self.vbox.Add(self.sizer6)
337
338	def set_layout(self):
339	"""
340	layout
341	"""
342	self.vbox.Layout()
343	self.vbox.Fit(self)
344	self.SetSizer(self.vbox)
345	self.Centre()
346
347	def set_owner(self,owner):
348	"""
349	set owner of fitpage
350
351	:param owner: the class responsible of plotting
352
353	"""
354	self.event_owner = owner
355	self.state.event_owner = owner
356
357	def get_state(self):
358	"""
359	"""
360	return self.state
361	def get_data(self):
362	"""
363	return the current data
364	"""
365	return self.data
366
367	def set_manager(self, manager):
368	"""
369	set panel manager
370
371	:param manager: instance of plugin fitting
372
373	"""
374	self._manager = manager
375	self.state.manager = manager
376
377	def populate_box(self, dict):
378	"""
379	Store list of model
380
381	:param dict: dictionary containing list of models
382
383	"""
384	self.model_list_box = dict
385	self.state.model_list_box = self.model_list_box
386	self.initialize_combox()
387
388	def initialize_combox(self):
389	"""
390	put default value in the combobox
391	"""
392	## fill combox box
393	if self.model_list_box is None:
394	return
395	if len(self.model_list_box) > 0:
396	self._populate_box(self.formfactorbox,
397	self.model_list_box["Shapes"])
398
399	if len(self.model_list_box) > 0:
400	self._populate_box(self.structurebox,
401	self.model_list_box["Structure Factors"])
402	self.structurebox.Insert("None", 0, None)
403	self.structurebox.SetSelection(0)
404	self.structurebox.Hide()
405	self.text2.Hide()
406	self.structurebox.Disable()
407	self.text2.Disable()
408
409	if self.model.__class__ in self.model_list_box["P(Q)*S(Q)"]:
410	self.structurebox.Show()
411	self.text2.Show()
412	self.structurebox.Enable()
413	self.text2.Enable()
414
415	def set_dispers_sizer(self):
416	"""
417	fill sizer containing dispersity info
418	"""
419	self.sizer4.Clear(True)
420	name="Polydispersity and Orientational Distribution"
421	box_description= wx.StaticBox(self, -1,name)
422	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
423	#----------------------------------------------------
424	self.disable_disp = wx.RadioButton(self, -1, 'Off', (10, 10),
425	style=wx.RB_GROUP)
426	self.enable_disp = wx.RadioButton(self, -1, 'On', (10, 30))
427
428
429	self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param,
430	id=self.disable_disp.GetId())
431	self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param,
432	id=self.enable_disp.GetId())
433	#MAC needs SetValue
434	self.disable_disp.SetValue(True)
435	sizer_dispersion = wx.BoxSizer(wx.HORIZONTAL)
436	sizer_dispersion.Add((20,20))
437	name=""#Polydispersity and \nOrientational Distribution "
438	sizer_dispersion.Add(wx.StaticText(self,-1,name))
439	sizer_dispersion.Add(self.enable_disp )
440	sizer_dispersion.Add((20,20))
441	sizer_dispersion.Add(self.disable_disp )
442	sizer_dispersion.Add((10,10))
443
444	## fill a sizer with the combobox to select dispersion type
445	sizer_select_dispers = wx.BoxSizer(wx.HORIZONTAL)
446	self.model_disp = wx.StaticText(self, -1, 'Distribution Function ')
447
448	import sans.models.dispersion_models
449	self.polydisp= sans.models.dispersion_models.models
450	self.disp_box = wx.ComboBox(self, -1)
451
452	for key, value in self.polydisp.iteritems():
453	name = str(key)
454	self.disp_box.Append(name,value)
455	self.disp_box.SetStringSelection("gaussian")
456	wx.EVT_COMBOBOX(self.disp_box,-1, self._on_select_Disp)
457
458	sizer_select_dispers.Add((10,10))
459	sizer_select_dispers.Add(self.model_disp)
460	sizer_select_dispers.Add(self.disp_box,0,
461	wx.TOP\|wx.BOTTOM\|wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE,border=5)
462
463	self.model_disp.Hide()
464	self.disp_box.Hide()
465
466	boxsizer1.Add( sizer_dispersion,0,
467	wx.TOP\|wx.BOTTOM\|wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE,border=5)
468	#boxsizer1.Add( (10,10) )
469	boxsizer1.Add( sizer_select_dispers )
470	self.sizer4_4 = wx.GridBagSizer(5,5)
471	boxsizer1.Add( self.sizer4_4 )
472	#-----------------------------------------------------
473	self.sizer4.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
474	self.sizer4_4.Layout()
475	self.sizer4.Layout()
476	self.Layout()
477
478	self.Refresh()
479	## saving the state of enable dispersity button
480	self.state.enable_disp= self.enable_disp.GetValue()
481	self.state.disable_disp= self.disable_disp.GetValue()
482
483	def select_disp_angle(self, event):
484	"""
485	Event for when a user select a parameter to average over.
486
487	:param event: radiobutton event
488
489	"""
490	self.values=[]
491	self.weights=[]
492	if event.GetEventObject()==self.noDisper_rbox:
493	if self.noDisper_rbox.GetValue():
494	#No array dispersity apply yet
495	self._reset_dispersity()
496	## Redraw the model ???
497	self._draw_model()
498	# Go through the list of dispersion check boxes to identify
499	# which one has changed
500	for p in self.disp_cb_dict:
501	self.state.disp_cb_dict[p]= self.disp_cb_dict[p].GetValue()
502	# Catch which one of the box was just checked or unchecked.
503	if event.GetEventObject() == self.disp_cb_dict[p]:
504	if self.disp_cb_dict[p].GetValue() == True:
505
506	##Temp. FIX for V1.0 regarding changing checkbox
507	#to radiobutton.
508	##This (self._reset_dispersity) should be removed
509	#when the array dispersion is fixed.
510	self._reset_dispersity()
511
512	# The user wants this parameter to be averaged.
513	# Pop up the file selection dialog.
514	path = self._selectDlg()
515
516	# If nothing was selected, just return
517	if path is None:
518	self.disp_cb_dict[p].SetValue(False)
519	self.noDisper_rbox.SetValue(True)
520	return
521	try:
522	self._default_save_location = os.path.dirname(path)
523	except:
524	pass
525	try:
526	self.values,self.weights = self.read_file(path)
527	except:
528	msg="Could not read input file"
529	wx.PostEvent(self.parent.parent,
530	StatusEvent(status=msg))
531	return
532
533	# If any of the two arrays is empty, notify the user that we won't
534	# proceed
535	if self.values is None or self.weights is None or \
536	self.values ==[] or self.weights ==[]:
537	msg = "The loaded %s distrubtion is"
538	msg + " corrupted or empty" % p
539	wx.PostEvent(self.parent.parent,
540	StatusEvent(status=msg))
541	return
542
543	# Tell the user that we are about to apply the distribution
544	msg = "Applying loaded %s distribution: %s" % (p, path)
545	wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
546
547	# Create the dispersion objects
548	from sans.models.dispersion_models import ArrayDispersion
549	disp_model = ArrayDispersion()
550	disp_model.set_weights(self.values, self.weights)
551
552	# Store the object to make it persist outside the
553	# scope of this method
554	#TODO: refactor model to clean this up?
555	self._disp_obj_dict[p] = disp_model
556	self.state._disp_obj_dict [p]= disp_model
557	self.state.values = []
558	self.state.weights = []
559	self.state.values = copy.deepcopy(self.values)
560	self.state.weights = copy.deepcopy(self.weights)
561	# Set the new model as the dispersion object for the
562	#selected parameter
563	self.model.set_dispersion(p, disp_model)
564	# Store a reference to the weights in the model object
565	#so that
566	# it's not lost when we use the model within another thread.
567	#TODO: total hack - fix this
568	self.state.model= self.model.clone()
569
570	self.model._persistency_dict = {}
571	self.model._persistency_dict[p] = \
572	[self.values, self.weights]
573	self.state.model._persistency_dict[p] = \
574	[self.values,self.weights]
575	else:
576	self._reset_dispersity()
577
578	## Redraw the model
579	self._draw_model()
580
581	## post state to fit panel
582	event = PageInfoEvent(page = self)
583	wx.PostEvent(self.parent, event)
584
585
586	def onResetModel(self, event):
587	"""
588	Reset model state
589	"""
590	menu = event.GetEventObject()
591	## post help message for the selected model
592	msg = menu.GetHelpString(event.GetId())
593	msg +=" reloaded"
594	wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
595
596	name = menu.GetLabel(event.GetId())
597	self._on_select_model_helper()
598
599	if name in self.saved_states.keys():
600	previous_state = self.saved_states[name]
601	## reset state of checkbox,textcrtl and regular parameters value
602	self.reset_page(previous_state)
603
604	def on_save(self, event):
605	"""
606	Save the current state into file
607	"""
608	self.save_current_state()
609	new_state = self.state.clone()
610	# Ask the user the location of the file to write to.
611	path = None
612	dlg = wx.FileDialog(self, "Choose a file", self._default_save_location,
613	"", "*.fitv", wx.SAVE)
614	if dlg.ShowModal() == wx.ID_OK:
615	path = dlg.GetPath()
616	self._default_save_location = os.path.dirname(path)
617	else:
618	return None
619	#the manager write the state into file
620	self._manager.save_fit_state(filepath=path, fitstate=new_state)
621	return new_state
622
623	def _get_time_stamp(self):
624	"""
625	return time and date stings
626	"""
627	# date and time
628	year, month, day,hour,minute,second,tda,ty,tm_isdst= time.localtime()
629	current_time= str(hour)+":"+str(minute)+":"+str(second)
630	current_date= str( month)+"/"+str(day)+"/"+str(year)
631	return current_time, current_date
632
633	def on_bookmark(self, event):
634	"""
635	save history of the data and model
636	"""
637	if self.model==None:
638	msg="Can not bookmark; Please select Data and Model first..."
639	wx.MessageBox(msg, 'Info')
640	return
641	self.save_current_state()
642	new_state = self.state.clone()
643	##Add model state on context menu
644	self.number_saved_state += 1
645	current_time, current_date = self._get_time_stamp()
646	#name= self.model.name+"[%g]"%self.number_saved_state
647	name = "Fitting: %g]" % self.number_saved_state
648	name += self.model.__class__.__name__
649	name += "bookmarked at %s on %s" % (current_time, current_date)
650	self.saved_states[name]= new_state
651
652	## Add item in the context menu
653	msg = "Model saved at %s on %s"%(current_time, current_date)
654	## post help message for the selected model
655	msg +=" Saved! right click on this page to retrieve this model"
656	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
657
658	#id = wx.NewId()
659	#self.popUpMenu.Append(id,name,str(msg))
660	#wx.EVT_MENU(self, id, self.onResetModel)
661	wx.PostEvent(self.parent.parent,
662	AppendBookmarkEvent(title=name,
663	hint=str(msg), handler=self.onResetModel))
664
665	def old_on_bookmark(self, event):
666	"""
667	save history of the data and model
668	"""
669	if self.model==None:
670	msg="Can not bookmark; Please select Data and Model first..."
671	wx.MessageBox(msg, 'Info')
672	return
673	if hasattr(self,"enable_disp"):
674	self.state.enable_disp = copy.deepcopy(self.enable_disp.GetValue())
675	if hasattr(self, "disp_box"):
676	self.state.disp_box = copy.deepcopy(self.disp_box.GetSelection())
677
678	self.state.model.name= self.model.name
679
680	#Remember fit engine_type for fit panel
681	if self.engine_type == None:
682	self.engine_type = "scipy"
683	if self._manager !=None:
684	self._manager._on_change_engine(engine=self.engine_type)
685
686	self.state.engine_type = self.engine_type
687
688	new_state = self.state.clone()
689	new_state.model.name = self.state.model.name
690
691	new_state.enable2D = copy.deepcopy(self.enable2D)
692	##Add model state on context menu
693	self.number_saved_state += 1
694	#name= self.model.name+"[%g]"%self.number_saved_state
695	name= self.model.__class__.__name__+"[%g]"%self.number_saved_state
696	self.saved_states[name]= new_state
697
698	## Add item in the context menu
699
700	year, month, day,hour,minute,second,tda,ty,tm_isdst= time.localtime()
701	my_time= str(hour)+" : "+str(minute)+" : "+str(second)+" "
702	date= str( month)+"\|"+str(day)+"\|"+str(year)
703	msg= "Model saved at %s on %s"%(my_time, date)
704	## post help message for the selected model
705	msg +=" Saved! right click on this page to retrieve this model"
706	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
707
708	id = wx.NewId()
709	self.popUpMenu.Append(id,name,str(msg))
710	wx.EVT_MENU(self, id, self.onResetModel)
711
712	def onSetFocus(self, evt):
713	"""
714	highlight the current textcrtl and hide the error text control shown
715	after fitting
716	"""
717	return
718
719	def read_file(self, path):
720	"""
721	Read two columns file
722
723	:param path: the path to the file to read
724
725	"""
726	try:
727	if path==None:
728	wx.PostEvent(self.parent.parent, StatusEvent(status=\
729	" Selected Distribution was not loaded: %s"%path))
730	return None, None
731	input_f = open(path, 'r')
732	buff = input_f.read()
733	lines = buff.split('\n')
734
735	angles = []
736	weights=[]
737	for line in lines:
738	toks = line.split()
739	try:
740	angle = float(toks[0])
741	weight = float(toks[1])
742	except:
743	# Skip non-data lines
744	pass
745	angles.append(angle)
746	weights.append(weight)
747	return numpy.array(angles), numpy.array(weights)
748	except:
749	raise
750
751	def createMemento(self):
752	"""
753	return the current state of the page
754	"""
755	return self.state.clone()
756
757
758	def save_current_state(self):
759	"""
760	Store current state
761	"""
762	self.state.engine_type = copy.deepcopy(self.engine_type)
763	## save model option
764	if self.model!= None:
765	self.disp_list= self.model.getDispParamList()
766	self.state.disp_list= copy.deepcopy(self.disp_list)
767	self.state.model = self.model.clone()
768	#save radiobutton state for model selection
769	self.state.shape_rbutton = self.shape_rbutton.GetValue()
770	self.state.shape_indep_rbutton = self.shape_indep_rbutton.GetValue()
771	self.state.struct_rbutton = self.struct_rbutton.GetValue()
772	self.state.plugin_rbutton = self.plugin_rbutton.GetValue()
773	#model combobox
774	self.state.structurebox = self.structurebox.GetSelection()
775	self.state.formfactorbox = self.formfactorbox.GetSelection()
776
777	self.state.enable2D = copy.deepcopy(self.enable2D)
778	self.state.values= copy.deepcopy(self.values)
779	self.state.weights = copy.deepcopy( self.weights)
780	## save data
781	self.state.data= copy.deepcopy(self.data)
782	self.state.qmax_x = self.qmax_x
783	self.state.qmin_x = self.qmin_x
784	try:
785	n = self.disp_box.GetCurrentSelection()
786	dispersity= self.disp_box.GetClientData(n)
787	name= dispersity.__name__
788	self.disp_name = name
789	if name == "GaussianDispersion" :
790	if hasattr(self,"cb1"):
791	self.state.cb1= self.cb1.GetValue()
792	except:
793	pass
794
795	if hasattr(self,"enable_disp"):
796	self.state.enable_disp= self.enable_disp.GetValue()
797	self.state.disable_disp = self.disable_disp.GetValue()
798
799	self.state.smearer = copy.deepcopy(self.smearer)
800	if hasattr(self,"enable_smearer"):
801	self.state.enable_smearer = \
802	copy.deepcopy(self.enable_smearer.GetValue())
803	self.state.disable_smearer = \
804	copy.deepcopy(self.disable_smearer.GetValue())
805
806	self.state.pinhole_smearer = \
807	copy.deepcopy(self.pinhole_smearer.GetValue())
808	self.state.slit_smearer = copy.deepcopy(self.slit_smearer.GetValue())
809
810	if hasattr(self,"disp_box"):
811	self.state.disp_box = self.disp_box.GetSelection()
812
813	if len(self.disp_cb_dict)>0:
814	for k , v in self.disp_cb_dict.iteritems():
815
816	if v ==None :
817	self.state.disp_cb_dict[k]= v
818	else:
819	try:
820	self.state.disp_cb_dict[k]=v.GetValue()
821	except:
822	self.state.disp_cb_dict[k]= None
823
824	if len(self._disp_obj_dict)>0:
825	for k , v in self._disp_obj_dict.iteritems():
826
827	self.state._disp_obj_dict[k]= v
828
829
830	self.state.values = copy.deepcopy(self.values)
831	self.state.weights = copy.deepcopy(self.weights)
832	## save plotting range
833	self._save_plotting_range()
834
835	self.state.orientation_params =[]
836	self.state.orientation_params_disp =[]
837	self.state.parameters =[]
838	self.state.fittable_param =[]
839	self.state.fixed_param =[]
840
841
842	## save checkbutton state and txtcrtl values
843	self._copy_parameters_state(self.orientation_params,
844	self.state.orientation_params)
845	self._copy_parameters_state(self.orientation_params_disp,
846	self.state.orientation_params_disp)
847
848	self._copy_parameters_state(self.parameters, self.state.parameters)
849	self._copy_parameters_state(self.fittable_param,
850	self.state.fittable_param)
851	self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
852	#save chisqr
853	self.state.tcChi = self.tcChi.GetValue()
854
855	def save_current_state_fit(self):
856	"""
857	Store current state for fit_page
858	"""
859	## save model option
860	if self.model!= None:
861	self.disp_list= self.model.getDispParamList()
862	self.state.disp_list= copy.deepcopy(self.disp_list)
863	self.state.model = self.model.clone()
864	if hasattr(self, "engine_type"):
865	self.state.engine_type = copy.deepcopy(self.engine_type)
866
867	self.state.enable2D = copy.deepcopy(self.enable2D)
868	self.state.values= copy.deepcopy(self.values)
869	self.state.weights = copy.deepcopy( self.weights)
870	## save data
871	self.state.data= copy.deepcopy(self.data)
872	try:
873	n = self.disp_box.GetCurrentSelection()
874	dispersity= self.disp_box.GetClientData(n)
875	name= dispersity.__name__
876	self.disp_name = name
877	if name == "GaussianDispersion" :
878	if hasattr(self,"cb1"):
879	self.state.cb1= self.cb1.GetValue()
880
881	except:
882	pass
883
884	if hasattr(self,"enable_disp"):
885	self.state.enable_disp= self.enable_disp.GetValue()
886	self.state.disable_disp = self.disable_disp.GetValue()
887
888	self.state.smearer = copy.deepcopy(self.smearer)
889	if hasattr(self,"enable_smearer"):
890	self.state.enable_smearer = \
891	copy.deepcopy(self.enable_smearer.GetValue())
892	self.state.disable_smearer = \
893	copy.deepcopy(self.disable_smearer.GetValue())
894
895	self.state.pinhole_smearer = \
896	copy.deepcopy(self.pinhole_smearer.GetValue())
897	self.state.slit_smearer = copy.deepcopy(self.slit_smearer.GetValue())
898
899	if hasattr(self,"disp_box"):
900	self.state.disp_box = self.disp_box.GetCurrentSelection()
901
902	if len(self.disp_cb_dict) > 0:
903	for k, v in self.disp_cb_dict.iteritems():
904
905	if v == None :
906	self.state.disp_cb_dict[k] = v
907	else:
908	try:
909	self.state.disp_cb_dict[k] = v.GetValue()
910	except:
911	self.state.disp_cb_dict[k] = None
912
913	if len(self._disp_obj_dict) > 0:
914	for k , v in self._disp_obj_dict.iteritems():
915
916	self.state._disp_obj_dict[k] = v
917
918
919	self.state.values = copy.deepcopy(self.values)
920	self.state.weights = copy.deepcopy(self.weights)
921
922	## save plotting range
923	self._save_plotting_range()
924
925	## save checkbutton state and txtcrtl values
926	self._copy_parameters_state(self.orientation_params,
927	self.state.orientation_params)
928	self._copy_parameters_state(self.orientation_params_disp,
929	self.state.orientation_params_disp)
930	self._copy_parameters_state(self.parameters, self.state.parameters)
931	self._copy_parameters_state(self.fittable_param,
932	self.state.fittable_param)
933	self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
934
935
936	def check_invalid_panel(self):
937	"""
938	check if the user can already perform some action with this panel
939	"""
940	flag = False
941	if self.data is None:
942	self.disable_smearer.SetValue(True)
943	self.disable_disp.SetValue(True)
944	msg = "Please load Data and select Model to start..."
945	wx.MessageBox(msg, 'Info')
946	return True
947
948
949	def reset_page_helper(self, state):
950	"""
951	Use page_state and change the state of existing page
952
953	:precondition: the page is already drawn or created
954
955	:postcondition: the state of the underlying data change as well as the
956	state of the graphic interface
957	"""
958	if state == None:
959	#self._undo.Enable(False)
960	return
961
962	self.set_data(state.data)
963	self.enable2D= state.enable2D
964	self.engine_type = state.engine_type
965
966	self.disp_cb_dict = state.disp_cb_dict
967	self.disp_list = state.disp_list
968
969	## set the state of the radio box
970	self.shape_rbutton.SetValue(state.shape_rbutton )
971	self.shape_indep_rbutton.SetValue(state.shape_indep_rbutton)
972	self.struct_rbutton.SetValue(state.struct_rbutton)
973	self.plugin_rbutton.SetValue(state.plugin_rbutton)
974
975	## fill model combobox
976	self._show_combox_helper()
977	#select the current model
978	self.formfactorbox.Select(int(state.formfactorcombobox))
979	self.structurebox.SetSelection(state.structurecombobox )
980	if state.multi_factor != None:
981	self.multifactorbox.SetSelection(state.multi_factor)
982
983	#reset the fitting engine type
984	self.engine_type = state.engine_type
985	#draw the pnael according to the new model parameter
986	self._on_select_model(event=None)
987
988	if self._manager !=None:
989	self._manager._on_change_engine(engine=self.engine_type)
990	## set the select all check box to the a given state
991	self.cb1.SetValue(state.cb1)
992
993	## reset state of checkbox,textcrtl and regular parameters value
994	self._reset_parameters_state(self.orientation_params_disp,
995	state.orientation_params_disp)
996	self._reset_parameters_state(self.orientation_params,
997	state.orientation_params)
998	self._reset_parameters_state(self.str_parameters,
999	state.str_parameters)
1000	self._reset_parameters_state(self.parameters,state.parameters)
1001	## display dispersion info layer
1002	self.enable_disp.SetValue(state.enable_disp)
1003	self.disable_disp.SetValue(state.disable_disp)
1004
1005	if hasattr(self, "disp_box"):
1006
1007	self.disp_box.SetSelection(state.disp_box)
1008	n= self.disp_box.GetCurrentSelection()
1009	dispersity= self.disp_box.GetClientData(n)
1010	name = dispersity.__name__
1011
1012	self._set_dipers_Param(event=None)
1013
1014	if name == "ArrayDispersion":
1015
1016	for item in self.disp_cb_dict.keys():
1017
1018	if hasattr(self.disp_cb_dict[item], "SetValue") :
1019	self.disp_cb_dict[item].SetValue(\
1020	state.disp_cb_dict[item])
1021	# Create the dispersion objects
1022	from sans.models.dispersion_models import ArrayDispersion
1023	disp_model = ArrayDispersion()
1024	if hasattr(state,"values")and\
1025	self.disp_cb_dict[item].GetValue() == True:
1026	if len(state.values)>0:
1027	self.values=state.values
1028	self.weights=state.weights
1029	disp_model.set_weights(self.values,
1030	state.weights)
1031	else:
1032	self._reset_dispersity()
1033
1034	self._disp_obj_dict[item] = disp_model
1035	# Set the new model as the dispersion object
1036	#for the selected parameter
1037	self.model.set_dispersion(item, disp_model)
1038
1039	self.model._persistency_dict[item] = \
1040	[state.values, state.weights]
1041
1042	else:
1043	keys = self.model.getParamList()
1044	for item in keys:
1045	if item in self.disp_list and \
1046	not self.model.details.has_key(item):
1047	self.model.details[item] = ["", None, None]
1048	for k,v in self.state.disp_cb_dict.iteritems():
1049	self.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1050	self.state.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1051
1052	##plotting range restore
1053	self._reset_plotting_range(state)
1054	## smearing info restore
1055	if hasattr(self, "enable_smearer"):
1056	## set smearing value whether or not the data
1057	#contain the smearing info
1058	self.enable_smearer.SetValue(state.enable_smearer)
1059	self.disable_smearer.SetValue(state.disable_smearer)
1060	self.onSmear(event=None)
1061	self.pinhole_smearer.SetValue(state.pinhole_smearer)
1062	self.slit_smearer.SetValue(state.slit_smearer)
1063	## we have two more options for smearing
1064	if self.pinhole_smearer.GetValue(): self.onPinholeSmear(event=None)
1065	elif self.slit_smearer.GetValue(): self.onSlitSmear(event=None)
1066
1067	## reset state of checkbox,textcrtl and dispersity parameters value
1068	self._reset_parameters_state(self.fittable_param,state.fittable_param)
1069	self._reset_parameters_state(self.fixed_param,state.fixed_param)
1070
1071	## draw the model with previous parameters value
1072	self._onparamEnter_helper()
1073	#reset the value of chisqr when not consistent with the value computed
1074	self.tcChi.SetValue(str(self.state.tcChi))
1075	## reset context menu items
1076	self._reset_context_menu()
1077
1078	## set the value of the current state to the state given as parameter
1079	self.state = state.clone()
1080
1081
1082	def old_reset_page_helper(self, state):
1083	"""
1084	Use page_state and change the state of existing page
1085
1086	:precondition: the page is already drawn or created
1087
1088	:postcondition: the state of the underlying data change as well as the
1089	state of the graphic interface
1090	"""
1091	if state ==None:
1092	#self._undo.Enable(False)
1093	return
1094
1095	self.model= state.model
1096	self.data = state.data
1097	if self.data !=None:
1098	from DataLoader.qsmearing import smear_selection
1099	self.smearer= smear_selection(self.data, self.model)
1100	self.enable2D= state.enable2D
1101	self.engine_type = state.engine_type
1102
1103	self.disp_cb_dict = state.disp_cb_dict
1104	self.disp_list = state.disp_list
1105
1106	## set the state of the radio box
1107	self.shape_rbutton.SetValue(state.shape_rbutton )
1108	self.shape_indep_rbutton.SetValue(state.shape_indep_rbutton)
1109	self.struct_rbutton.SetValue(state.struct_rbutton )
1110	self.plugin_rbutton.SetValue(state.plugin_rbutton)
1111	##draw sizer containing model parameters value for the current model
1112	self._set_model_sizer_selection( self.model )
1113	self.set_model_param_sizer(self.model)
1114
1115	## reset value of combox box
1116	self.structurebox.SetSelection(state.structurecombobox )
1117	self.formfactorbox.SetSelection(state.formfactorcombobox)
1118
1119
1120	## enable the view 2d button if this is a modelpage type
1121	if hasattr(self,"model_view"):
1122	if self.enable2D:
1123	self.model_view.Disable()
1124	else:
1125	self.model_view.Enable()
1126	## set the select all check box to the a given state
1127	if hasattr(self, "cb1"):
1128	self.cb1.SetValue(state.cb1)
1129
1130	## reset state of checkbox,textcrtl and regular parameters value
1131
1132	self._reset_parameters_state(self.orientation_params_disp,
1133	state.orientation_params_disp)
1134	self._reset_parameters_state(self.orientation_params,
1135	state.orientation_params)
1136	self._reset_parameters_state(self.parameters,state.parameters)
1137	## display dispersion info layer
1138	self.enable_disp.SetValue(state.enable_disp)
1139	self.disable_disp.SetValue(state.disable_disp)
1140
1141	if hasattr(self, "disp_box"):
1142
1143	self.disp_box.SetSelection(state.disp_box)
1144	n= self.disp_box.GetCurrentSelection()
1145	dispersity= self.disp_box.GetClientData(n)
1146	name= dispersity.__name__
1147
1148	self._set_dipers_Param(event=None)
1149
1150	if name=="ArrayDispersion":
1151
1152	for item in self.disp_cb_dict.keys():
1153
1154	if hasattr(self.disp_cb_dict[item],"SetValue") :
1155	self.disp_cb_dict[item].SetValue(state.disp_cb_dict[item])
1156	# Create the dispersion objects
1157	from sans.models.dispersion_models import ArrayDispersion
1158	disp_model = ArrayDispersion()
1159	if hasattr(state,"values") and\
1160	self.disp_cb_dict[item].GetValue()==True:
1161	if len(state.values) > 0:
1162	self.values = state.values
1163	self.weights = state.weights
1164	disp_model.set_weights(self.values,
1165	state.weights)
1166	else:
1167	self._reset_dispersity()
1168
1169	self._disp_obj_dict[item] = disp_model
1170	# Set the new model as the dispersion
1171	#object for the selected parameter
1172	self.model.set_dispersion(item, disp_model)
1173
1174	self.model._persistency_dict[item] = [state.values,
1175	state.weights]
1176
1177	else:
1178	keys = self.model.getParamList()
1179	for item in keys:
1180	if item in self.disp_list and \
1181	not self.model.details.has_key(item):
1182	self.model.details[item]=["",None,None]
1183	for k,v in self.state.disp_cb_dict.iteritems():
1184	self.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1185	self.state.disp_cb_dict = copy.deepcopy(state.disp_cb_dict)
1186
1187	##plotting range restore
1188	self._reset_plotting_range(state)
1189
1190	## smearing info restore
1191	if hasattr(self,"enable_smearer"):
1192	## set smearing value whether or not the data
1193	#contain the smearing info
1194	self.enable_smearer.SetValue(state.enable_smearer)
1195	self.disable_smearer.SetValue(state.disable_smearer)
1196	self.onSmear(event=None)
1197	self.pinhole_smearer.SetValue(state.pinhole_smearer)
1198	self.slit_smearer.SetValue(state.slit_smearer)
1199	## we have two more options for smearing
1200	if self.pinhole_smearer.GetValue(): self.onPinholeSmear(event=None)
1201	elif self.slit_smearer.GetValue(): self.onSlitSmear(event=None)
1202
1203	## reset state of checkbox,textcrtl and dispersity parameters value
1204	self._reset_parameters_state(self.fittable_param,state.fittable_param)
1205	self._reset_parameters_state(self.fixed_param,state.fixed_param)
1206
1207	## draw the model with previous parameters value
1208	self._onparamEnter_helper()
1209
1210	## reset context menu items
1211	self._reset_context_menu()
1212
1213	## set the value of the current state to the state given as parameter
1214	self.state = state.clone()
1215	self._draw_model()
1216
1217	def _selectDlg(self):
1218	"""
1219	open a dialog file to selected the customized dispersity
1220	"""
1221	import os
1222	dlg = wx.FileDialog(self, "Choose a weight file",
1223	self._default_save_location , "",
1224	".", wx.OPEN)
1225	path = None
1226	if dlg.ShowModal() == wx.ID_OK:
1227	path = dlg.GetPath()
1228	dlg.Destroy()
1229	return path
1230
1231	def _reset_context_menu(self):
1232	"""
1233	reset the context menu
1234	"""
1235	for name, state in self.state.saved_states.iteritems():
1236	self.number_saved_state += 1
1237	## Add item in the context menu
1238	id = wx.NewId()
1239	msg = 'Save model and state %g' % self.number_saved_state
1240	self.popUpMenu.Append(id, name, msg)
1241	wx.EVT_MENU(self, id, self.onResetModel)
1242
1243	def _reset_plotting_range(self, state):
1244	"""
1245	Reset the plotting range to a given state
1246	"""
1247	# if self.check_invalid_panel():
1248	# return
1249	self.qmin_tctrl.SetValue(str(state.qmin))
1250	self.qmax.SetValue(str(state.qmax))
1251
1252	def _save_typeOfmodel(self):
1253	"""
1254	save radiobutton containing the type model that can be selected
1255	"""
1256	self.state.shape_rbutton = self.shape_rbutton.GetValue()
1257	self.state.shape_indep_rbutton = self.shape_indep_rbutton.GetValue()
1258	self.state.struct_rbutton = self.struct_rbutton.GetValue()
1259	self.state.plugin_rbutton = self.plugin_rbutton.GetValue()
1260	self.state.structurebox= self.structurebox.GetCurrentSelection()
1261	self.state.formfactorbox = self.formfactorbox.GetCurrentSelection()
1262
1263	#self._undo.Enable(True)
1264	## post state to fit panel
1265	event = PageInfoEvent(page = self)
1266	wx.PostEvent(self.parent, event)
1267
1268	def _save_plotting_range(self ):
1269	"""
1270	save the state of plotting range
1271	"""
1272	self.state.qmin = self.qmin_x
1273	self.state.qmax = self.qmax_x
1274	self.state.npts = self.npts_x
1275
1276	def _onparamEnter_helper(self):
1277	"""
1278	check if values entered by the user are changed and valid to replot
1279	model
1280	"""
1281	# Flag to register when a parameter has changed.
1282	is_modified = False
1283	self.fitrange = True
1284	is_2Ddata = False
1285	#self._undo.Enable(True)
1286	# check if 2d data
1287	if self.data.__class__.__name__ == "Data2D":
1288	is_2Ddata = True
1289	if self.model !=None:
1290	try:
1291	is_modified = self._check_value_enter(self.fittable_param,
1292	is_modified)
1293	is_modified = self._check_value_enter(self.fixed_param,
1294	is_modified)
1295	is_modified = self._check_value_enter(self.parameters,
1296	is_modified)
1297	except:
1298	pass
1299	#if is_modified:
1300
1301	# Here we should check whether the boundaries have been modified.
1302	# If qmin and qmax have been modified, update qmin and qmax and
1303	# set the is_modified flag to True
1304	if self._validate_qrange(self.qmin_tcrl, self.qmax):
1305	tempmin = float(self.qmin_tcrl.GetValue())
1306	if tempmin != self.qmin_x:
1307	self.qmin_x = tempmin
1308	is_modified = True
1309	tempmax = float(self.qmax.GetValue())
1310	if tempmax != self.qmax_x:
1311	self.qmax_x = tempmax
1312	is_modified = True
1313
1314	if is_2Ddata:
1315	# set mask
1316	is_modified = self._validate_Npts()
1317
1318	else:
1319	self.fitrange = False
1320
1321	## if any value is modify draw model with new value
1322	if not self.fitrange:
1323	#self.btFit.Disable()
1324	if is_2Ddata: self.btEditMask.Disable()
1325	else:
1326	#self.btFit.Enable(True)
1327	if is_2Ddata: self.btEditMask.Enable(True)
1328
1329	if is_modified and self.fitrange:
1330	self.state_change= True
1331	self._draw_model()
1332	self.Refresh()
1333	return is_modified
1334
1335	def _update_paramv_on_fit(self):
1336	"""
1337	make sure that update param values just before the fitting
1338	"""
1339	#flag for qmin qmax check values
1340	flag = True
1341	self.fitrange = True
1342	is_modified = False
1343
1344	wx.PostEvent(self._manager.parent, StatusEvent(status=" \
1345	updating ... ",type="update"))
1346
1347	##So make sure that update param values on_Fit.
1348	#self._undo.Enable(True)
1349	if self.model !=None:
1350	##Check the values
1351	self._check_value_enter( self.fittable_param ,is_modified)
1352	self._check_value_enter( self.fixed_param ,is_modified)
1353	self._check_value_enter( self.parameters ,is_modified)
1354
1355	# If qmin and qmax have been modified, update qmin and qmax and
1356	# Here we should check whether the boundaries have been modified.
1357	# If qmin and qmax have been modified, update qmin and qmax and
1358	# set the is_modified flag to True
1359	self.fitrange = self._validate_qrange(self.qmin_tcrl, self.qmax)
1360	if self.fitrange:
1361	tempmin = float(self.qmin_tcrl.GetValue())
1362	if tempmin != self.qmin_x:
1363	self.qmin_x = tempmin
1364	tempmax = float(self.qmax.GetValue())
1365	if tempmax != self.qmax_x:
1366	self.qmax_x = tempmax
1367	if tempmax == tempmin:
1368	flag = False
1369	temp_smearer = None
1370	if not self.disable_smearer.GetValue():
1371	temp_smearer= self.current_smearer
1372	if self.slit_smearer.GetValue():
1373	flag = self.update_slit_smear()
1374	elif self.pinhole_smearer.GetValue():
1375	flag = self.update_pinhole_smear()
1376	else:
1377	self._manager.set_smearer(smearer=temp_smearer,
1378	id=self.id,
1379	qmin=float(self.qmin_x),
1380	qmax=float(self.qmax_x),
1381	draw=False)
1382	elif not self._is_2D():
1383	self._manager.set_smearer(smearer=temp_smearer,
1384	qmin=float(self.qmin_x),
1385	id=self.id,
1386	qmax= float(self.qmax_x))
1387	index_data = ((self.qmin_x <= self.data.x)&\
1388	(self.data.x <= self.qmax_x))
1389	val = str(len(self.data.x[index_data==True]))
1390	self.Npts_fit.SetValue(val)
1391	flag = True
1392	if self._is_2D():
1393	# only 2D case set mask
1394	flag = self._validate_Npts()
1395	if not flag:
1396	return flag
1397	else: flag = False
1398	else:
1399	flag = False
1400
1401	#For invalid q range, disable the mask editor and fit button, vs.
1402	if not self.fitrange:
1403	#self.btFit.Disable()
1404	if self._is_2D():self.btEditMask.Disable()
1405	else:
1406	#self.btFit.Enable(True)
1407	if self._is_2D():self.btEditMask.Enable(True)
1408
1409	if not flag:
1410	msg = "Cannot Plot or Fit :Must select a "
1411	msg += " model or Fitting range is not valid!!! "
1412	wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
1413
1414	self.save_current_state()
1415
1416	return flag
1417
1418	def _is_modified(self, is_modified):
1419	"""
1420	return to self._is_modified
1421	"""
1422	return is_modified
1423
1424	def _reset_parameters_state(self, listtorestore, statelist):
1425	"""
1426	Reset the parameters at the given state
1427	"""
1428	if len(statelist) == 0 or len(listtorestore) == 0:
1429	return
1430	if len(statelist) != len(listtorestore):
1431	return
1432
1433	for j in range(len(listtorestore)):
1434	item_page = listtorestore[j]
1435	item_page_info = statelist[j]
1436	##change the state of the check box for simple parameters
1437	if item_page[0]!=None:
1438	item_page[0].SetValue(item_page_info[0])
1439	if item_page[2]!=None:
1440	item_page[2].SetValue(item_page_info[2])
1441	if item_page[3]!=None:
1442	## show or hide text +/-
1443	if item_page_info[2]:
1444	item_page[3].Show(True)
1445	else:
1446	item_page[3].Hide()
1447	if item_page[4]!=None:
1448	## show of hide the text crtl for fitting error
1449	if item_page_info[4][0]:
1450	item_page[4].Show(True)
1451	item_page[4].SetValue(item_page_info[4][1])
1452	else:
1453	item_page[3].Hide()
1454	if item_page[5]!=None:
1455	## show of hide the text crtl for fitting error
1456	if item_page_info[5][0]:
1457	item_page[5].Show(True)
1458	item_page[5].SetValue(item_page_info[5][1])
1459	else:
1460	item_page[5].Hide()
1461
1462	if item_page[6]!=None:
1463	## show of hide the text crtl for fitting error
1464	if item_page_info[6][0]:
1465	item_page[6].Show(True)
1466	item_page[6].SetValue(item_page_info[6][1])
1467	else:
1468	item_page[6].Hide()
1469
1470	def _copy_parameters_state(self, listtocopy, statelist):
1471	"""
1472	copy the state of button
1473
1474	:param listtocopy: the list of check button to copy
1475	:param statelist: list of state object to store the current state
1476
1477	"""
1478	if len(listtocopy)==0:
1479	return
1480
1481	for item in listtocopy:
1482
1483	checkbox_state = None
1484	if item[0]!= None:
1485	checkbox_state= item[0].GetValue()
1486	parameter_name = item[1]
1487	parameter_value = None
1488	if item[2]!=None:
1489	parameter_value = item[2].GetValue()
1490	static_text = None
1491	if item[3]!=None:
1492	static_text = item[3].IsShown()
1493	error_value = None
1494	error_state = None
1495	if item[4]!= None:
1496	error_value = item[4].GetValue()
1497	error_state = item[4].IsShown()
1498
1499	min_value = None
1500	min_state = None
1501	if item[5]!= None:
1502	min_value = item[5].GetValue()
1503	min_state = item[5].IsShown()
1504
1505	max_value = None
1506	max_state = None
1507	if item[6]!= None:
1508	max_value = item[6].GetValue()
1509	max_state = item[6].IsShown()
1510	unit=None
1511	if item[7]!=None:
1512	unit = item[7].GetLabel()
1513
1514	statelist.append([checkbox_state, parameter_name, parameter_value,
1515	static_text ,[error_state, error_value],
1516	[min_state, min_value],
1517	[max_state, max_value], unit])
1518
1519	def _set_model_sizer_selection(self, model):
1520	"""
1521	Display the sizer according to the type of the current model
1522	"""
1523	if model == None:
1524	return
1525	if hasattr(model ,"s_model"):
1526
1527	class_name = model.s_model.__class__
1528	name = model.s_model.name
1529	flag = (name != "NoStructure")
1530	if flag and \
1531	(class_name in self.model_list_box["Structure Factors"]):
1532	self.structurebox.Show()
1533	self.text2.Show()
1534	self.structurebox.Enable()
1535	self.text2.Enable()
1536	items = self.structurebox.GetItems()
1537	self.sizer1.Layout()
1538
1539	for i in range(len(items)):
1540	if items[i]== str(name):
1541	self.structurebox.SetSelection(i)
1542	break
1543
1544	if hasattr(model ,"p_model"):
1545	class_name = model.p_model.__class__
1546	name = model.p_model.name
1547	self.formfactorbox.Clear()
1548
1549	for k, list in self.model_list_box.iteritems():
1550	if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
1551	self.shape_rbutton.SetValue(True)
1552	## fill the form factor list with new model
1553	self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
1554	items = self.formfactorbox.GetItems()
1555	## set comboxbox to the selected item
1556	for i in range(len(items)):
1557	if items[i]== str(name):
1558	self.formfactorbox.SetSelection(i)
1559	break
1560	return
1561	elif k == "Shape-Independent":
1562	self.shape_indep_rbutton.SetValue(True)
1563	elif k == "Structure Factors":
1564	self.struct_rbutton.SetValue(True)
1565	elif k == "Multi-Functions":
1566	continue
1567	else:
1568	self.plugin_rbutton.SetValue(True)
1569
1570	if class_name in list:
1571	## fill the form factor list with new model
1572	self._populate_box(self.formfactorbox, list)
1573	items = self.formfactorbox.GetItems()
1574	## set comboxbox to the selected item
1575	for i in range(len(items)):
1576	if items[i]== str(name):
1577	self.formfactorbox.SetSelection(i)
1578	break
1579	break
1580	else:
1581
1582	## Select the model from the menu
1583	class_name = model.__class__
1584	name = model.name
1585	self.formfactorbox.Clear()
1586	items = self.formfactorbox.GetItems()
1587
1588	for k, list in self.model_list_box.iteritems():
1589	if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
1590	if class_name in self.model_list_box["P(Q)*S(Q)"]:
1591	self.structurebox.Show()
1592	self.text2.Show()
1593	self.structurebox.Enable()
1594	self.structurebox.SetSelection(0)
1595	self.text2.Enable()
1596	else:
1597	self.structurebox.Hide()
1598	self.text2.Hide()
1599	self.structurebox.Disable()
1600	self.structurebox.SetSelection(0)
1601	self.text2.Disable()
1602
1603	self.shape_rbutton.SetValue(True)
1604	## fill the form factor list with new model
1605	self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
1606	items = self.formfactorbox.GetItems()
1607	## set comboxbox to the selected item
1608	for i in range(len(items)):
1609	if items[i]== str(name):
1610	self.formfactorbox.SetSelection(i)
1611	break
1612	return
1613	elif k == "Shape-Independent":
1614	self.shape_indep_rbutton.SetValue(True)
1615	elif k == "Structure Factors":
1616	self.struct_rbutton.SetValue(True)
1617	elif k == "Multi-Functions":
1618	continue
1619	else:
1620	self.plugin_rbutton.SetValue(True)
1621	if class_name in list:
1622	self.structurebox.SetSelection(0)
1623	self.structurebox.Disable()
1624	self.text2.Disable()
1625	## fill the form factor list with new model
1626	self._populate_box(self.formfactorbox, list)
1627	items = self.formfactorbox.GetItems()
1628	## set comboxbox to the selected item
1629	for i in range(len(items)):
1630	if items[i]== str(name):
1631	self.formfactorbox.SetSelection(i)
1632	break
1633	break
1634
1635	def _draw_model(self):
1636	"""
1637	Method to draw or refresh a plotted model.
1638	The method will use the data member from the model page
1639	to build a call to the fitting perspective manager.
1640	"""
1641	#if self.check_invalid_panel():
1642	# return
1643	if self.model !=None:
1644	temp_smear=None
1645	if hasattr(self, "enable_smearer"):
1646	if not self.disable_smearer.GetValue():
1647	temp_smear= self.current_smearer
1648	self._manager.draw_model(self.model,
1649	data=self.data,
1650	smearer= temp_smear,
1651	qmin=float(self.qmin_x),
1652	qmax=float(self.qmax_x),
1653	qstep= float(self.npts_x),
1654	id=self.id,
1655	enable2D=self.enable2D)
1656
1657	def _set_model_sizer(self,sizer, box_sizer, title="", object=None):
1658	"""
1659	Use lists to fill a sizer for model info
1660	"""
1661
1662	sizer.Clear(True)
1663	##For MAC, this should defined here.
1664	if box_sizer == None:
1665	box_description= wx.StaticBox(self, -1,str(title))
1666	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
1667	else:
1668	boxsizer1 = box_sizer
1669	sizer_buttons = wx.BoxSizer(wx.HORIZONTAL)
1670	#--------------------------------------------------------
1671	self.shape_rbutton = wx.RadioButton(self, -1, 'Shapes', style=wx.RB_GROUP)
1672	self.shape_indep_rbutton = wx.RadioButton(self, -1, "Shape-Independent")
1673	self.struct_rbutton = wx.RadioButton(self, -1, "Structure Factor ")
1674	self.plugin_rbutton = wx.RadioButton(self, -1, "Customized Models")
1675
1676	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
1677	id= self.shape_rbutton.GetId() )
1678	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
1679	id= self.shape_indep_rbutton.GetId() )
1680	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
1681	id= self.struct_rbutton.GetId() )
1682	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
1683	id= self.plugin_rbutton.GetId() )
1684	#MAC needs SetValue
1685	self.shape_rbutton.SetValue(True)
1686
1687	sizer_radiobutton = wx.GridSizer(2, 2,5, 5)
1688	sizer_radiobutton.Add(self.shape_rbutton)
1689	sizer_radiobutton.Add(self.shape_indep_rbutton)
1690
1691	sizer_radiobutton.Add(self.plugin_rbutton)
1692	sizer_radiobutton.Add(self.struct_rbutton)
1693	sizer_buttons.Add(sizer_radiobutton)
1694	# detail button
1695	if object !=None:
1696	sizer_buttons.Add((50,0))
1697	sizer_buttons.Add(object)
1698
1699	sizer_selection = wx.BoxSizer(wx.HORIZONTAL)
1700	mutifactor_selection = wx.BoxSizer(wx.HORIZONTAL)
1701
1702	self.text1 = wx.StaticText( self,-1,"" )
1703	self.text2 = wx.StaticText( self,-1,"P(Q)*S(Q)" )
1704	self.mutifactor_text = wx.StaticText( self,-1,"No. of Shells: " )
1705	self.mutifactor_text1 = wx.StaticText( self,-1,"" )
1706	self.show_sld_button = wx.Button( self,-1,"Show SLD Profile" )
1707	self.show_sld_button.Bind(wx.EVT_BUTTON,self._on_show_sld)
1708
1709	self.formfactorbox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
1710	if self.model!= None:
1711	self.formfactorbox.SetValue(self.model.name)
1712
1713	self.structurebox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
1714	self.multifactorbox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
1715	self.initialize_combox()
1716
1717	wx.EVT_COMBOBOX(self.formfactorbox,-1, self._on_select_model)
1718	wx.EVT_COMBOBOX(self.structurebox,-1, self._on_select_model)
1719	wx.EVT_COMBOBOX(self.multifactorbox,-1, self._on_select_model)
1720
1721	## check model type to show sizer
1722	if self.model !=None:
1723	self._set_model_sizer_selection( self.model )
1724
1725	sizer_selection.Add(self.text1)
1726	sizer_selection.Add((5,5))
1727	sizer_selection.Add(self.formfactorbox)
1728	sizer_selection.Add((5,5))
1729	sizer_selection.Add(self.text2)
1730	sizer_selection.Add((5,5))
1731	sizer_selection.Add(self.structurebox)
1732	#sizer_selection.Add((5,5))
1733	mutifactor_selection.Add((10,5))
1734	mutifactor_selection.Add(self.mutifactor_text)
1735	mutifactor_selection.Add(self.multifactorbox)
1736	mutifactor_selection.Add((5,5))
1737	mutifactor_selection.Add(self.mutifactor_text1)
1738	mutifactor_selection.Add((10,5))
1739	mutifactor_selection.Add(self.show_sld_button)
1740
1741	boxsizer1.Add( sizer_buttons )
1742	boxsizer1.Add( (15,15))
1743	boxsizer1.Add( sizer_selection )
1744	boxsizer1.Add( (10,10))
1745	boxsizer1.Add(mutifactor_selection)
1746
1747	self._set_multfactor_combobox()
1748	self.multifactorbox.SetSelection(1)
1749	self.show_sld_button.Hide()
1750	#--------------------------------------------------------
1751	sizer.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
1752	sizer.Layout()
1753
1754	def _on_show_sld(self, event=None):
1755	"""
1756	Plot SLD profile
1757	"""
1758	# get profile data
1759	x,y=self.model.getProfile()
1760
1761	from danse.common.plottools import Data1D
1762	#from sans.perspectives.theory.profile_dialog import SLDPanel
1763	from sans.guiframe.local_perspectives.plotting.profile_dialog \
1764	import SLDPanel
1765	sld_data = Data1D(x,y)
1766	sld_data.name = 'SLD'
1767	sld_data.axes = self.sld_axes
1768	self.panel = SLDPanel(self, data=sld_data,axes =self.sld_axes,id =-1 )
1769	self.panel.ShowModal()
1770
1771	def _set_multfactor_combobox(self, multiplicity=10):
1772	"""
1773	Set comboBox for muitfactor of CoreMultiShellModel
1774	:param multiplicit: no. of multi-functionality
1775	"""
1776	# build content of the combobox
1777	for idx in range(0,multiplicity):
1778	self.multifactorbox.Append(str(idx),int(idx))
1779	#self.multifactorbox.SetSelection(1)
1780	self._hide_multfactor_combobox()
1781
1782	def _show_multfactor_combobox(self):
1783	"""
1784	Show the comboBox of muitfactor of CoreMultiShellModel
1785	"""
1786	if not self.mutifactor_text.IsShown():
1787	self.mutifactor_text.Show(True)
1788	self.mutifactor_text1.Show(True)
1789	if not self.multifactorbox.IsShown():
1790	self.multifactorbox.Show(True)
1791
1792	def _hide_multfactor_combobox(self):
1793	"""
1794	Hide the comboBox of muitfactor of CoreMultiShellModel
1795	"""
1796	if self.mutifactor_text.IsShown():
1797	self.mutifactor_text.Hide()
1798	self.mutifactor_text1.Hide()
1799	if self.multifactorbox.IsShown():
1800	self.multifactorbox.Hide()
1801
1802
1803	def _show_combox_helper(self):
1804	"""
1805	Fill panel's combo box according to the type of model selected
1806	"""
1807	if self.shape_rbutton.GetValue():
1808	##fill the combobox with form factor list
1809	self.structurebox.SetSelection(0)
1810	self.structurebox.Disable()
1811	self.formfactorbox.Clear()
1812	self._populate_box( self.formfactorbox,self.model_list_box["Shapes"])
1813	if self.shape_indep_rbutton.GetValue():
1814	##fill the combobox with shape independent factor list
1815	self.structurebox.SetSelection(0)
1816	self.structurebox.Disable()
1817	self.formfactorbox.Clear()
1818	self._populate_box( self.formfactorbox,
1819	self.model_list_box["Shape-Independent"])
1820	if self.struct_rbutton.GetValue():
1821	##fill the combobox with structure factor list
1822	self.structurebox.SetSelection(0)
1823	self.structurebox.Disable()
1824	self.formfactorbox.Clear()
1825	self._populate_box( self.formfactorbox,
1826	self.model_list_box["Structure Factors"])
1827	if self.plugin_rbutton.GetValue():
1828	##fill the combobox with form factor list
1829	self.structurebox.Disable()
1830	self.formfactorbox.Clear()
1831	self._populate_box( self.formfactorbox,
1832	self.model_list_box["Customized Models"])
1833
1834	def _show_combox(self, event=None):
1835	"""
1836	Show combox box associate with type of model selected
1837	"""
1838	#if self.check_invalid_panel():
1839	# self.shape_rbutton.SetValue(True)
1840	# return
1841
1842	self._show_combox_helper()
1843	self._on_select_model(event=None)
1844	self._save_typeOfmodel()
1845	self.sizer4_4.Layout()
1846	self.sizer4.Layout()
1847	self.Layout()
1848	self.Refresh()
1849
1850	def _populate_box(self, combobox, list):
1851	"""
1852	fill combox box with dict item
1853
1854	:param list: contains item to fill the combox
1855	item must model class
1856	"""
1857	st = time.time()
1858	for models in list:
1859	model= models()
1860	name = model.__class__.__name__
1861	if models.__name__!="NoStructure":
1862	if hasattr(model, "name"):
1863	name = model.name
1864	combobox.Append(name,models)
1865	return 0
1866
1867	def _onQrangeEnter(self, event):
1868	"""
1869	Check validity of value enter in the Q range field
1870
1871	"""
1872	tcrtl = event.GetEventObject()
1873	#Clear msg if previously shown.
1874	msg = ""
1875	wx.PostEvent(self.parent, StatusEvent(status=msg))
1876	# Flag to register when a parameter has changed.
1877	is_modified = False
1878	if tcrtl.GetValue().lstrip().rstrip() != "":
1879	try:
1880	value = float(tcrtl.GetValue())
1881	tcrtl.SetBackgroundColour(wx.WHITE)
1882	# If qmin and qmax have been modified, update qmin and qmax
1883	if self._validate_qrange(self.qmin_tcrl, self.qmax):
1884	tempmin = float(self.qmin_tcrl.GetValue())
1885	if tempmin != self.qmin_x:
1886	self.qmin_x = tempmin
1887	tempmax = float(self.qmax.GetValue())
1888	if tempmax != self.qmax_x:
1889	self.qmax_x = tempmax
1890	else:
1891	tcrtl.SetBackgroundColour("pink")
1892	msg = "Model Error:wrong value entered : %s" % sys.exc_value
1893	wx.PostEvent(self.parent, StatusEvent(status=msg))
1894	return
1895	except:
1896	tcrtl.SetBackgroundColour("pink")
1897	msg = "Model Error:wrong value entered : %s" % sys.exc_value
1898	wx.PostEvent(self.parent, StatusEvent(status=msg))
1899	return
1900	#Check if # of points for theory model are valid(>0).
1901	if self.npts != None:
1902	if check_float(self.npts):
1903	temp_npts = float(self.npts.GetValue())
1904	if temp_npts != self.num_points:
1905	self.num_points = temp_npts
1906	is_modified = True
1907	else:
1908	msg = "Cannot Plot :No npts in that Qrange!!! "
1909	wx.PostEvent(self.parent, StatusEvent(status=msg))
1910	else:
1911	tcrtl.SetBackgroundColour("pink")
1912	msg = "Model Error:wrong value entered!!!"
1913	wx.PostEvent(self.parent, StatusEvent(status=msg))
1914	#self._undo.Enable(True)
1915	self.save_current_state()
1916	event = PageInfoEvent(page=self)
1917	wx.PostEvent(self.parent, event)
1918	self.state_change = False
1919	#Draw the model for a different range
1920	self._draw_model()
1921
1922	def _theory_qrange_enter(self, event):
1923	"""
1924	Check validity of value enter in the Q range field
1925	"""
1926
1927	tcrtl= event.GetEventObject()
1928	#Clear msg if previously shown.
1929	msg= ""
1930	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1931	# Flag to register when a parameter has changed.
1932	is_modified = False
1933	if tcrtl.GetValue().lstrip().rstrip()!="":
1934	try:
1935	value = float(tcrtl.GetValue())
1936	tcrtl.SetBackgroundColour(wx.WHITE)
1937
1938	# If qmin and qmax have been modified, update qmin and qmax
1939	if self._validate_qrange(self.theory_qmin, self.theory_qmax):
1940	tempmin = float(self.theory_qmin.GetValue())
1941	if tempmin != self.theory_qmin_x:
1942	self.theory_qmin_x = tempmin
1943	tempmax = float(self.theory_qmax.GetValue())
1944	if tempmax != self.qmax_x:
1945	self.theory_qmax_x = tempmax
1946	else:
1947	tcrtl.SetBackgroundColour("pink")
1948	msg= "Model Error:wrong value entered : %s"% sys.exc_value
1949	wx.PostEvent(self._manager.parent, StatusEvent(status = msg ))
1950	return
1951	except:
1952	tcrtl.SetBackgroundColour("pink")
1953	msg= "Model Error:wrong value entered : %s"% sys.exc_value
1954	wx.PostEvent(self._manager.parent, StatusEvent(status = msg ))
1955	return
1956	#Check if # of points for theory model are valid(>0).
1957	if self.theory_npts != None:
1958	if check_float(self.theory_npts):
1959	temp_npts = float(self.theory_npts.GetValue())
1960	if temp_npts != self.num_points:
1961	self.num_points = temp_npts
1962	is_modified = True
1963	else:
1964	msg= "Cannot Plot :No npts in that Qrange!!! "
1965	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1966	else:
1967	tcrtl.SetBackgroundColour("pink")
1968	msg = "Model Error:wrong value entered!!!"
1969	wx.PostEvent(self.parent.parent, StatusEvent(status=msg))
1970	#self._undo.Enable(True)
1971	self.save_current_state()
1972	event = PageInfoEvent(page = self)
1973	wx.PostEvent(self.parent, event)
1974	self.state_change= False
1975	#Draw the model for a different range
1976	self._draw_model()
1977
1978	def _on_select_model_helper(self):
1979	"""
1980	call back for model selection
1981	"""
1982	## reset dictionary containing reference to dispersion
1983	self._disp_obj_dict = {}
1984	self.disp_cb_dict ={}
1985
1986	f_id = self.formfactorbox.GetCurrentSelection()
1987	#For MAC
1988	form_factor = None
1989	if f_id >= 0:
1990	form_factor = self.formfactorbox.GetClientData(f_id)
1991
1992	if not form_factor in self.model_list_box["multiplication"]:
1993	self.structurebox.Hide()
1994	self.text2.Hide()
1995	self.structurebox.Disable()
1996	self.structurebox.SetSelection(0)
1997	self.text2.Disable()
1998	else:
1999	self.structurebox.Show()
2000	self.text2.Show()
2001	self.structurebox.Enable()
2002	self.text2.Enable()
2003
2004	if form_factor != None:
2005	# set multifactor for Mutifunctional models
2006	if form_factor().__class__ in self.model_list_box["Multi-Functions"]:
2007	m_id = self.multifactorbox.GetCurrentSelection()
2008	multiplicity = form_factor().multiplicity_info[0]
2009	self.multifactorbox.Clear()
2010	#self.mutifactor_text.SetLabel(form_factor().details[])
2011	self._set_multfactor_combobox(multiplicity)
2012	self._show_multfactor_combobox()
2013	#ToDo: this info should be called directly from the model
2014	text = form_factor().multiplicity_info[1]#'No. of Shells: '
2015
2016	#self.mutifactor_text.Clear()
2017	self.mutifactor_text.SetLabel(text)
2018	if m_id > multiplicity -1:
2019	# default value
2020	m_id = 1
2021
2022	self.multi_factor = self.multifactorbox.GetClientData(m_id)
2023	if self.multi_factor == None: self.multi_factor =0
2024	form_factor = form_factor(int(self.multi_factor))
2025	self.multifactorbox.SetSelection(m_id)
2026	# Check len of the text1 and max_multiplicity
2027	text = ''
2028	if form_factor.multiplicity_info[0] == len(form_factor.multiplicity_info[2]):
2029	text = form_factor.multiplicity_info[2][self.multi_factor]
2030	self.mutifactor_text1.SetLabel(text)
2031	# Check if model has get sld profile.
2032	if len(form_factor.multiplicity_info[3]) > 0:
2033	self.sld_axes = form_factor.multiplicity_info[3]
2034	self.show_sld_button.Show(True)
2035	else:
2036	self.sld_axes = ""
2037
2038	else:
2039	self._hide_multfactor_combobox()
2040	self.show_sld_button.Hide()
2041	form_factor = form_factor()
2042	self.multi_factor = None
2043	else:
2044	self._hide_multfactor_combobox()
2045	self.show_sld_button.Hide()
2046	self.multi_factor = None
2047
2048	s_id = self.structurebox.GetCurrentSelection()
2049	struct_factor = self.structurebox.GetClientData( s_id )
2050
2051	if struct_factor !=None:
2052	from sans.models.MultiplicationModel import MultiplicationModel
2053	self.model= MultiplicationModel(form_factor,struct_factor())
2054
2055	else:
2056	if form_factor != None:
2057	self.model= form_factor
2058	else:
2059	self.model = None
2060	return self.model
2061
2062
2063	## post state to fit panel
2064	self.state.parameters =[]
2065	self.state.model =self.model
2066	self.state.qmin = self.qmin_x
2067	self.state.multi_factor = self.multi_factor
2068	self.disp_list =self.model.getDispParamList()
2069	self.state.disp_list = self.disp_list
2070	self.Layout()
2071
2072	def _validate_qrange(self, qmin_ctrl, qmax_ctrl):
2073	"""
2074	Verify that the Q range controls have valid values
2075	and that Qmin < Qmax.
2076
2077	:param qmin_ctrl: text control for Qmin
2078	:param qmax_ctrl: text control for Qmax
2079
2080	:return: True is the Q range is value, False otherwise
2081
2082	"""
2083	qmin_validity = check_float(qmin_ctrl)
2084	qmax_validity = check_float(qmax_ctrl)
2085	if not (qmin_validity and qmax_validity):
2086	return False
2087	else:
2088	qmin = float(qmin_ctrl.GetValue())
2089	qmax = float(qmax_ctrl.GetValue())
2090	if qmin <= qmax:
2091	#Make sure to set both colours white.
2092	qmin_ctrl.SetBackgroundColour(wx.WHITE)
2093	qmin_ctrl.Refresh()
2094	qmax_ctrl.SetBackgroundColour(wx.WHITE)
2095	qmax_ctrl.Refresh()
2096	else:
2097	qmin_ctrl.SetBackgroundColour("pink")
2098	qmin_ctrl.Refresh()
2099	qmax_ctrl.SetBackgroundColour("pink")
2100	qmax_ctrl.Refresh()
2101	msg= "Invalid Q range: Q min must be smaller than Q max"
2102	wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
2103	return False
2104	return True
2105
2106	def _validate_Npts(self):
2107	"""
2108	Validate the number of points for fitting is more than 10 points.
2109	If valid, setvalues Npts_fit otherwise post msg.
2110	"""
2111	#default flag
2112	flag = True
2113
2114	# q value from qx and qy
2115	radius= numpy.sqrt( self.data.qx_data * self.data.qx_data +
2116	self.data.qy_data * self.data.qy_data )
2117	#get unmasked index
2118	index_data = (float(self.qmin_tcrl.GetValue()) <= radius) & \
2119	(radius <= float(self.qmax.GetValue()))
2120	index_data = (index_data) & (self.data.mask)
2121	index_data = (index_data) & (numpy.isfinite(self.data.data))
2122
2123	if len(index_data[index_data]) < 10:
2124	# change the color pink.
2125	self.qmin_tcrl.SetBackgroundColour("pink")
2126	self.qmin_tcrl.Refresh()
2127	self.qmax.SetBackgroundColour("pink")
2128	self.qmax.Refresh()
2129	msg= "Cannot Plot :No or too little npts in that data range!!! "
2130	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
2131	self.fitrange = False
2132	flag = False
2133	else:
2134	self.Npts_fit.SetValue(str(len(self.data.mask[index_data==True])))
2135	self.fitrange = True
2136
2137	return flag
2138
2139	def _check_value_enter(self, list, modified):
2140	"""
2141	:param list: model parameter and panel info
2142	:Note: each item of the list should be as follow:
2143	item=[check button state, parameter's name,
2144	paramater's value, string="+/-",
2145	parameter's error of fit,
2146	parameter's minimum value,
2147	parrameter's maximum value ,
2148	parameter's units]
2149	"""
2150	is_modified = modified
2151	if len(list)==0:
2152	return is_modified
2153	for item in list:
2154	#skip angle parameters for 1D
2155	if self.data.__class__.__name__ !="Data2D":
2156	if item in self.orientation_params:
2157	continue
2158	#try:
2159	name = str(item[1])
2160
2161	if string.find(name,".npts") ==-1 and string.find(name,".nsigmas")==-1:
2162	## check model parameters range
2163	param_min= None
2164	param_max= None
2165
2166	## check minimun value
2167	if item[5]!= None and item[5]!= "":
2168	if item[5].GetValue().lstrip().rstrip()!="":
2169	try:
2170
2171	param_min = float(item[5].GetValue())
2172	if not self._validate_qrange(item[5],item[2]):
2173	if numpy.isfinite(param_min):
2174	item[2].SetValue(format_number(param_min))
2175
2176	item[5].SetBackgroundColour(wx.WHITE)
2177	item[2].SetBackgroundColour(wx.WHITE)
2178
2179	except:
2180	msg = "Wrong Fit parameter range entered "
2181	wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
2182	raise ValueError, msg
2183	is_modified = True
2184	## check maximum value
2185	if item[6]!= None and item[6]!= "":
2186	if item[6].GetValue().lstrip().rstrip()!="":
2187	try:
2188	param_max = float(item[6].GetValue())
2189	if not self._validate_qrange(item[2],item[6]):
2190	if numpy.isfinite(param_max):
2191	item[2].SetValue(format_number(param_max))
2192
2193	item[6].SetBackgroundColour(wx.WHITE)
2194	item[2].SetBackgroundColour(wx.WHITE)
2195	except:
2196	msg = "Wrong Fit parameter range entered "
2197	wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
2198	raise ValueError, msg
2199	is_modified = True
2200
2201
2202	if param_min != None and param_max !=None:
2203	if not self._validate_qrange(item[5], item[6]):
2204	msg= "Wrong Fit range entered for parameter "
2205	msg+= "name %s of model %s "%(name, self.model.name)
2206	wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
2207
2208	if name in self.model.details.keys():
2209	self.model.details[name][1:3]= param_min,param_max
2210	is_modified = True
2211
2212	else:
2213	self.model.details [name] = ["",param_min,param_max]
2214	is_modified = True
2215	try:
2216	value= float(item[2].GetValue())
2217
2218	# If the value of the parameter has changed,
2219	# +update the model and set the is_modified flag
2220	if value != self.model.getParam(name) and numpy.isfinite(value):
2221	self.model.setParam(name,value)
2222	is_modified = True
2223	except:
2224	msg = "Wrong Fit parameter value entered "
2225	wx.PostEvent(self.parent.parent, StatusEvent(status = msg))
2226
2227	return is_modified
2228
2229
2230	def _set_dipers_Param(self, event):
2231	"""
2232	respond to self.enable_disp and self.disable_disp radio box.
2233	The dispersity object is reset inside the model into Gaussian.
2234	When the user select yes , this method display a combo box for more selection
2235	when the user selects No,the combo box disappears.
2236	Redraw the model with the default dispersity (Gaussian)
2237	"""
2238	#if self.check_invalid_panel():
2239	# return
2240	## On selction if no model exists.
2241	if self.model ==None:
2242	self.disable_disp.SetValue(True)
2243	msg="Please select a Model first..."
2244	wx.MessageBox(msg, 'Info')
2245	wx.PostEvent(self._manager.parent, StatusEvent(status=\
2246	"Polydispersion: %s"%msg))
2247	return
2248
2249	self._reset_dispersity()
2250
2251	if self.model ==None:
2252	self.model_disp.Hide()
2253	self.disp_box.Hide()
2254	self.sizer4_4.Clear(True)
2255	return
2256
2257	if self.enable_disp.GetValue():
2258	self.model_disp.Show(True)
2259	self.disp_box.Show(True)
2260	## layout for model containing no dispersity parameters
2261
2262	self.disp_list= self.model.getDispParamList()
2263
2264	if len(self.disp_list)==0 and len(self.disp_cb_dict)==0:
2265	self._layout_sizer_noDipers()
2266	else:
2267	## set gaussian sizer
2268	self._on_select_Disp(event=None)
2269	else:
2270	self.model_disp.Hide()
2271	self.disp_box.Hide()
2272	self.disp_box.SetSelection(0)
2273	self.sizer4_4.Clear(True)
2274
2275	## post state to fit panel
2276	self.save_current_state()
2277	if event !=None:
2278	#self._undo.Enable(True)
2279	event = PageInfoEvent(page = self)
2280	wx.PostEvent(self.parent, event)
2281	#draw the model with the current dispersity
2282	self._draw_model()
2283	self.sizer4_4.Layout()
2284	self.sizer5.Layout()
2285	self.Layout()
2286	self.Refresh()
2287
2288
2289	def _layout_sizer_noDipers(self):
2290	"""
2291	Draw a sizer with no dispersity info
2292	"""
2293	ix=0
2294	iy=1
2295	self.fittable_param=[]
2296	self.fixed_param=[]
2297	self.orientation_params_disp=[]
2298
2299	self.model_disp.Hide()
2300	self.disp_box.Hide()
2301	self.sizer4_4.Clear(True)
2302	text = "No polydispersity available for this model"
2303	model_disp = wx.StaticText(self, -1, text)
2304	self.sizer4_4.Add(model_disp,( iy, ix),(1,1), wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
2305	self.sizer4_4.Layout()
2306	self.sizer4.Layout()
2307
2308	def _reset_dispersity(self):
2309	"""
2310	put gaussian dispersity into current model
2311	"""
2312	if len(self.param_toFit)>0:
2313	for item in self.fittable_param:
2314	if item in self.param_toFit:
2315	self.param_toFit.remove(item)
2316
2317	for item in self.orientation_params_disp:
2318	if item in self.param_toFit:
2319	self.param_toFit.remove(item)
2320
2321	self.fittable_param=[]
2322	self.fixed_param=[]
2323	self.orientation_params_disp=[]
2324	self.values=[]
2325	self.weights=[]
2326
2327	from sans.models.dispersion_models import GaussianDispersion, ArrayDispersion
2328	if len(self.disp_cb_dict)==0:
2329	self.save_current_state()
2330	self.sizer4_4.Clear(True)
2331	self.Layout()
2332
2333	return
2334	if (len(self.disp_cb_dict)>0) :
2335	for p in self.disp_cb_dict:
2336	# The parameter was un-selected. Go back to Gaussian model (with 0 pts)
2337	disp_model = GaussianDispersion()
2338
2339	self._disp_obj_dict[p] = disp_model
2340	# Set the new model as the dispersion object for the selected parameter
2341	try:
2342	self.model.set_dispersion(p, disp_model)
2343	except:
2344
2345	pass
2346
2347	## save state into
2348	self.save_current_state()
2349	self.Layout()
2350	self.Refresh()
2351
2352	def _on_select_Disp(self,event):
2353	"""
2354	allow selecting different dispersion
2355	self.disp_list should change type later .now only gaussian
2356	"""
2357	n = self.disp_box.GetCurrentSelection()
2358	name = self.disp_box.GetValue()
2359	dispersity= self.disp_box.GetClientData(n)
2360	self.disp_name = name
2361
2362	if name.lower() == "array":
2363	self._set_sizer_arraydispersion()
2364	else:
2365	self._set_sizer_dispersion(dispersity= dispersity)
2366
2367	self.state.disp_box= n
2368	## Redraw the model
2369	self._draw_model()
2370	#self._undo.Enable(True)
2371	event = PageInfoEvent(page = self)
2372	wx.PostEvent(self.parent, event)
2373
2374	self.sizer4_4.Layout()
2375	self.sizer4.Layout()
2376
2377	def _set_sizer_arraydispersion(self):
2378	"""
2379	draw sizer with array dispersity parameters
2380	"""
2381
2382	if len(self.param_toFit)>0:
2383	for item in self.fittable_param:
2384	if item in self.param_toFit:
2385	self.param_toFit.remove(item)
2386	for item in self.orientation_params_disp:
2387	if item in self.param_toFit:
2388	self.param_toFit.remove(item)
2389	for item in self.model.details.keys():
2390	if item in self.model.fixed:
2391	del self.model.details [item]
2392
2393	self.fittable_param=[]
2394	self.fixed_param=[]
2395	self.orientation_params_disp=[]
2396	self.sizer4_4.Clear(True)
2397	self._reset_dispersity()
2398	ix=0
2399	iy=1
2400	disp1 = wx.StaticText(self, -1, 'Array Dispersion')
2401	self.sizer4_4.Add(disp1,( iy, ix),(1,1), wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
2402
2403	# Look for model parameters to which we can apply an ArrayDispersion model
2404	# Add a check box for each parameter.
2405	self.disp_cb_dict = {}
2406	ix+=1
2407	self.noDisper_rbox = wx.RadioButton(self, -1,"None", (10, 10),style= wx.RB_GROUP)
2408	self.Bind(wx.EVT_RADIOBUTTON,self.select_disp_angle , id=self.noDisper_rbox.GetId())
2409	#MAC needs SetValue
2410	self.noDisper_rbox.SetValue(True)
2411	self.sizer4_4.Add(self.noDisper_rbox, (iy, ix),
2412	(1,1), wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2413
2414	for p in self.model.dispersion.keys():
2415	if not p in self.model.orientation_params:
2416	ix+=1
2417	self.disp_cb_dict[p] = wx.RadioButton(self, -1, p, (10, 10))
2418	self.state.disp_cb_dict[p]= self.disp_cb_dict[p].GetValue()
2419	self.Bind(wx.EVT_RADIOBUTTON, self.select_disp_angle, id=self.disp_cb_dict[p].GetId())
2420	self.sizer4_4.Add(self.disp_cb_dict[p], (iy, ix), (1,1), wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2421
2422	for p in self.model.dispersion.keys():
2423	if p in self.model.orientation_params:
2424	ix+=1
2425	self.disp_cb_dict[p] = wx.RadioButton(self, -1, p, (10, 10))
2426	self.state.disp_cb_dict[p]= self.disp_cb_dict[p].GetValue()
2427	if not (self.enable2D or self.data.__class__.__name__ =="Data2D"):
2428	self.disp_cb_dict[p].Hide()
2429	else:
2430	self.disp_cb_dict[p].Show(True)
2431	self.Bind(wx.EVT_RADIOBUTTON, self.select_disp_angle, id=self.disp_cb_dict[p].GetId())
2432	self.sizer4_4.Add(self.disp_cb_dict[p], (iy, ix), (1,1), wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2433
2434
2435	ix =0
2436	iy +=1
2437	self.sizer4_4.Add((20,20),(iy,ix),(1,1), wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
2438	self.Layout()
2439
2440	self.state.orientation_params =[]
2441	self.state.orientation_params_disp =[]
2442	self.state.parameters =[]
2443	self.state.fittable_param =[]
2444	self.state.fixed_param =[]
2445
2446	## save checkbutton state and txtcrtl values
2447
2448	self._copy_parameters_state(self.orientation_params,
2449	self.state.orientation_params)
2450
2451	self._copy_parameters_state(self.orientation_params_disp,
2452	self.state.orientation_params_disp)
2453
2454	self._copy_parameters_state(self.parameters, self.state.parameters)
2455	self._copy_parameters_state(self.fittable_param, self.state.fittable_param)
2456	self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
2457
2458
2459	## post state to fit panel
2460	event = PageInfoEvent(page = self)
2461	wx.PostEvent(self.parent, event)
2462
2463	def _lay_out(self):
2464	"""
2465	returns self.Layout
2466
2467	:Note: Mac seems to like this better when self.
2468	Layout is called after fitting.
2469	"""
2470	self._sleep4sec()
2471	self.Layout()
2472	return
2473
2474	def _sleep4sec(self):
2475	"""
2476	sleep for 1 sec only applied on Mac
2477	Note: This 1sec helps for Mac not to crash on self.:ayout after self._draw_model
2478	"""
2479	if ON_MAC == True:
2480	time.sleep(1)
2481
2482	def on_reset_clicked(self,event):
2483	"""
2484	On 'Reset' button for Q range clicked
2485	"""
2486	flag = True
2487	#if self.check_invalid_panel():
2488	# return
2489	##For 3 different cases: Data2D, Data1D, and theory
2490	if self.data.__class__.__name__ == "Data2D":
2491	data_min= 0
2492	x= max(math.fabs(self.data.xmin), math.fabs(self.data.xmax))
2493	y= max(math.fabs(self.data.ymin), math.fabs(self.data.ymax))
2494	self.qmin_x = data_min
2495	self.qmax_x = math.sqrt(xx + yy)
2496	# check smearing
2497	if not self.disable_smearer.GetValue():
2498	temp_smearer= self.current_smearer
2499	## set smearing value whether or not the data contain the smearing info
2500	if self.pinhole_smearer.GetValue():
2501	flag = self.update_pinhole_smear()
2502	else:
2503	flag = True
2504	elif self.data.__class__.__name__ != "Data2D":
2505	self.qmin_x = min(self.data.x)
2506	self.qmax_x = max(self.data.x)
2507	# check smearing
2508	if not self.disable_smearer.GetValue():
2509	temp_smearer= self.current_smearer
2510	## set smearing value whether or not the data contain the smearing info
2511	if self.slit_smearer.GetValue():
2512	flag = self.update_slit_smear()
2513	elif self.pinhole_smearer.GetValue():
2514	flag = self.update_pinhole_smear()
2515	else:
2516	flag = True
2517	else:
2518	self.qmin_x = _QMIN_DEFAULT
2519	self.qmax_x = _QMAX_DEFAULT
2520	self.num_points = _NPTS_DEFAULT
2521	self.state.npts = self.num_points
2522
2523	if flag == False:
2524	msg= "Cannot Plot :Must enter a number!!! "
2525	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
2526	else:
2527	# set relative text ctrs.
2528	self.qmin_tcrl.SetValue(str(self.qmin_x))
2529	self.qmax.SetValue(str(self.qmax_x))
2530	self.set_npts2fit()
2531	# At this point, some button and variables satatus (disabled?) should be checked
2532	# such as color that should be reset to white in case that it was pink.
2533	self._onparamEnter_helper()
2534
2535	self.save_current_state()
2536	self.state.qmin = self.qmin_x
2537	self.state.qmax = self.qmax_x
2538
2539	#reset the q range values
2540	self._reset_plotting_range(self.state)
2541	#self.compute_chisqr(smearer=self.current_smearer)
2542	#Re draw plot
2543	self._draw_model()
2544
2545	def on_model_help_clicked(self,event):
2546	"""
2547	on 'More details' button
2548	"""
2549	from help_panel import HelpWindow
2550
2551	if self.model == None:
2552	name = 'FuncHelp'
2553	else:
2554	name = self.model.origin_name
2555
2556	frame = HelpWindow(None, -1, pageToOpen="media/model_functions.html")
2557	frame.Show(True)
2558	if frame.rhelp.HasAnchor(name):
2559	frame.rhelp.ScrollToAnchor(name)
2560	else:
2561	msg= "Model does not contains an available description "
2562	msg +="Please try searching in the Help window"
2563	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
2564

Note: See TracBrowser for help on using the repository browser.

SasView

source: sasview/sansview/perspectives/fitting/basepage.py @ 883e5f5

Download in other formats: