basepage.py @ fc6ea43

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 fc6ea43 was fc6ea43, checked in by Gervaise Alina <gervyh@…>, 15 years ago
disable and enable dispersity change
Property mode set to `100644`
File size: 51.4 KB

Line
1
2	import sys, os
3	import wx
4	import numpy
5	import time
6	import copy
7	from sans.guiframe.utils import format_number,check_float
8	from sans.guicomm.events import StatusEvent
9	import pagestate
10	from pagestate import PageState
11	(PageInfoEvent, EVT_PAGE_INFO) = wx.lib.newevent.NewEvent()
12	_BOX_WIDTH = 80
13
14	class BasicPage(wx.ScrolledWindow):
15	"""
16	This class provide general structure of fitpanel page
17	"""
18	## Internal name for the AUI manager
19	window_name = "Basic Page"
20	## Title to appear on top of the window
21	window_caption = "Basic page "
22
23	def __init__(self,parent, page_info):
24	wx.ScrolledWindow.__init__(self, parent)
25	##window_name
26	self.window_name = page_info.window_name
27	##window_caption
28	self.window_caption = page_info.window_caption
29	## parent of the page
30	self.parent = parent
31	## manager is the fitting plugin
32	self.manager= page_info.manager
33	## owner of the page (fitting plugin)
34	self.event_owner= page_info.event_owner
35	## current model
36	self.model = page_info.model
37	## data
38	self.data = page_info.data
39	## dictionary containing list of models
40	self.model_list_box = page_info.model_list_box
41	## Data member to store the dispersion object created
42	self._disp_obj_dict = {}
43	## selected parameters to apply dispersion
44	self.disp_cb_dict ={}
45	## smearer object
46	self.smearer = None
47	##list of model parameters. each item must have same length
48	## each item related to a given parameters
49	##[cb state, name, value, "+/-", error of fit, min, max , units]
50	self.parameters=[]
51	## list of parameters to fit , must be like self.parameters
52	self.param_toFit=[]
53	## list of looking like parameters but with non fittable parameters info
54	self.fixed_param=[]
55	## list of looking like parameters but with fittable parameters info
56	self.fittable_param=[]
57	##list of dispersion parameters
58	self.disp_list=[]
59	## list of orientation parameters
60	self.orientation_params=[]
61	self.orientation_params_disp=[]
62	if self.model !=None:
63	self.disp_list= self.model.getDispParamList()
64	##enable model 2D draw
65	self.enable2D= False
66	## check that the fit range is correct to plot the model again
67	self.fitrange= True
68	## Q range
69	self.qmin_x= 0.001
70	self.qmax_x= 0.13
71	self.num_points= 100
72	## Create memento to save the current state
73
74	self.state= PageState(parent= self.parent,model=self.model, data=self.data)
75
76	## retrieve saved state
77	self.number_saved_state= 0
78	## dictionary of saved state
79	self.saved_states={}
80	self.slicerpop = wx.Menu()
81	## Default locations
82	self._default_save_location = os.getcwd()
83	## save initial state on context menu
84	self.onSave(event=None)
85	self.Bind(wx.EVT_CONTEXT_MENU, self.onContextMenu)
86
87	## create the basic structure of the panel with empty sizer
88	self.define_page_structure()
89	## drawing Initial dispersion parameters sizer
90	self.set_dispers_sizer()
91	self._fill_save_sizer()
92	## layout
93	self.set_layout()
94
95
96	def onContextMenu(self, event):
97	"""
98	Retrieve the state selected state
99	"""
100	pos = event.GetPosition()
101	pos = self.ScreenToClient(pos)
102	self.PopupMenu(self.slicerpop, pos)
103
104
105
106	def define_page_structure(self):
107	"""
108	Create empty sizer for a panel
109	"""
110	self.vbox = wx.BoxSizer(wx.VERTICAL)
111	self.sizer0 = wx.BoxSizer(wx.VERTICAL)
112	self.sizer1 = wx.BoxSizer(wx.VERTICAL)
113	self.sizer2 = wx.BoxSizer(wx.VERTICAL)
114	self.sizer3 = wx.BoxSizer(wx.VERTICAL)
115	self.sizer4 = wx.BoxSizer(wx.VERTICAL)
116	self.sizer5 = wx.BoxSizer(wx.VERTICAL)
117	self.sizer6 = wx.BoxSizer(wx.VERTICAL)
118
119	self.sizer0.SetMinSize((375,-1))
120	self.sizer1.SetMinSize((375,-1))
121	self.sizer2.SetMinSize((375,-1))
122	self.sizer3.SetMinSize((375,-1))
123	self.sizer4.SetMinSize((375,-1))
124	self.sizer5.SetMinSize((375,-1))
125	self.sizer6.SetMinSize((375,-1))
126
127	self.vbox.Add(self.sizer0)
128	self.vbox.Add(self.sizer1)
129	self.vbox.Add(self.sizer2)
130	self.vbox.Add(self.sizer3)
131	self.vbox.Add(self.sizer4)
132	self.vbox.Add(self.sizer5)
133	self.vbox.Add(self.sizer6)
134
135
136	def set_layout(self):
137	"""
138	layout
139	"""
140	self.vbox.Layout()
141	self.vbox.Fit(self)
142	self.SetSizer(self.vbox)
143
144	self.set_scroll()
145	self.Centre()
146
147
148	def set_scroll(self):
149	self.SetScrollbars(20,20,200,100)
150	self.Layout()
151	self.SetAutoLayout(True)
152
153
154	def set_owner(self,owner):
155	"""
156	set owner of fitpage
157	@param owner: the class responsible of plotting
158	"""
159	self.event_owner = owner
160	self.state.event_owner = owner
161
162	def set_manager(self, manager):
163	"""
164	set panel manager
165	@param manager: instance of plugin fitting
166	"""
167	self.manager = manager
168	self.state.manager = manager
169
170	def populate_box(self, dict):
171	"""
172	Store list of model
173	@param dict: dictionary containing list of models
174	"""
175	self.model_list_box = dict
176	self.state.model_list_box = self.model_list_box
177
178
179
180	def set_dispers_sizer(self):
181	"""
182	fill sizer containing dispersity info
183	"""
184	self.sizer4.Clear(True)
185	name="Polydispersity and Orientational Distribution"
186	box_description= wx.StaticBox(self, -1,name)
187	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
188	#----------------------------------------------------
189	self.disable_disp = wx.RadioButton(self, -1, 'Off', (10, 10), style=wx.RB_GROUP)
190	self.enable_disp = wx.RadioButton(self, -1, 'On', (10, 30))
191	## saving the state of enable dispersity button
192	self.state.enable_disp= self.enable_disp.GetValue()
193
194	self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param, id=self.disable_disp.GetId())
195	self.Bind(wx.EVT_RADIOBUTTON, self._set_dipers_Param, id=self.enable_disp.GetId())
196
197	sizer_dispersion = wx.BoxSizer(wx.HORIZONTAL)
198	sizer_dispersion.Add((20,20))
199	name=""#Polydispersity and \nOrientational Distribution "
200	sizer_dispersion.Add(wx.StaticText(self,-1,name))
201	sizer_dispersion.Add(self.enable_disp )
202	sizer_dispersion.Add((20,20))
203	sizer_dispersion.Add(self.disable_disp )
204	sizer_dispersion.Add((10,10))
205
206	## fill a sizer with the combobox to select dispersion type
207	sizer_select_dispers = wx.BoxSizer(wx.HORIZONTAL)
208	self.model_disp = wx.StaticText(self, -1, 'Distribution Function ')
209
210	import sans.models.dispersion_models
211	self.polydisp= sans.models.dispersion_models.models
212	self.disp_box = wx.ComboBox(self, -1)
213
214	for key in self.polydisp.iterkeys():
215	name = str(key.__name__)
216	if name=="ArrayDispersion":
217	# Remove the option until the rest of the code is ready for it
218	self.disp_box.Append("Select customized Model",key)
219
220	else:
221	self.disp_box.Append(name,key)
222	self.disp_box.SetSelection(0)
223
224	wx.EVT_COMBOBOX(self.disp_box,-1, self._on_select_Disp)
225
226	sizer_select_dispers.Add((10,10))
227	sizer_select_dispers.Add(self.model_disp)
228	sizer_select_dispers.Add(self.disp_box,0,
229	wx.TOP\|wx.BOTTOM\|wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE,border=5)
230	#sizer_select_dispers.Add((10,10))
231	self.model_disp.Hide()
232	self.disp_box.Hide()
233
234	boxsizer1.Add( sizer_dispersion,0,
235	wx.TOP\|wx.BOTTOM\|wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE,border=5)
236	#boxsizer1.Add( (10,10) )
237	boxsizer1.Add( sizer_select_dispers )
238	self.sizer4_4 = wx.GridBagSizer(5,5)
239	boxsizer1.Add( self.sizer4_4 )
240	#-----------------------------------------------------
241	self.sizer4.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
242	self.sizer4_4.Layout()
243	self.sizer4.Layout()
244	self.Layout()
245	self.SetScrollbars(20,20,200,100)
246	self.Refresh()
247
248
249	def select_disp_angle(self, event):
250	"""
251	Event for when a user select a parameter to average over.
252	@param event: check box event
253	"""
254
255
256	# Go through the list of dispersion check boxes to identify which one has changed
257	for p in self.disp_cb_dict:
258	# Catch which one of the box was just checked or unchecked.
259	if event.GetEventObject() == self.disp_cb_dict[p]:
260
261
262	if self.disp_cb_dict[p].GetValue() == True:
263	# The user wants this parameter to be averaged.
264	# Pop up the file selection dialog.
265	path = self._selectDlg()
266
267	# If nothing was selected, just return
268	if path is None:
269	self.disp_cb_dict[p].SetValue(False)
270	return
271	try:
272	self._default_save_location = os.path.dirname(path)
273	except:
274	pass
275	try:
276	values,weights = self.read_file(path)
277	except:
278	msg="Could not read input file"
279	wx.PostEvent(self.parent.parent, StatusEvent(status= msg))
280	return
281
282	# If any of the two arrays is empty, notify the user that we won't
283	# proceed
284	if values is None or weights is None:
285	wx.PostEvent(self.parent.parent, StatusEvent(status=\
286	"The loaded %s distrubtion is corrupted or empty" % p))
287	return
288
289	# Tell the user that we are about to apply the distribution
290	wx.PostEvent(self.parent.parent, StatusEvent(status=\
291	"Applying loaded %s distribution: %s" % (p, path)))
292
293	# Create the dispersion objects
294	from sans.models.dispersion_models import ArrayDispersion
295	disp_model = ArrayDispersion()
296	disp_model.set_weights(values, weights)
297	# Store the object to make it persist outside the scope of this method
298	#TODO: refactor model to clean this up?
299	self._disp_obj_dict[p] = disp_model
300	# Set the new model as the dispersion object for the selected parameter
301	self.model.set_dispersion(p, disp_model)
302
303	else:
304	# The parameter was un-selected. Go back to Gaussian model (with 0 pts)
305	self._reset_dispersity()
306
307	## Redraw the model
308	self._draw_model()
309	return
310
311
312	def onResetModel(self, event):
313	"""
314	Reset model state
315	"""
316	## post help message for the selected model
317	msg = self.slicerpop.GetHelpString(event.GetId())
318	msg +=" reloaded"
319	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
320
321	name= self.slicerpop.GetLabel(event.GetId())
322	if name in self.saved_states.keys():
323	previous_state = self.saved_states[name]
324	self.reset_page(previous_state)
325
326
327	def onSave(self, event):
328	"""
329	save history of the data and model
330	"""
331	if self.model==None:
332	return
333	if hasattr(self,"enable_disp"):
334	self.state.enable_disp = copy.deepcopy(self.enable_disp.GetValue())
335	if hasattr(self, "disp_box"):
336	self.state.disp_box = copy.deepcopy(self.disp_box.GetSelection())
337
338	self.state.model = self.model.clone()
339	new_state = self.state.clone()
340
341	##Add model state on context menu
342	self.number_saved_state += 1
343	name= self.model.name+"[%g]"%self.number_saved_state
344	self.saved_states[name]= new_state
345
346	## Add item in the context menu
347
348	year, month, day,hour,minute,second,tda,ty,tm_isdst= time.localtime()
349	my_time= str(hour)+" : "+str(minute)+" : "+str(second)+" "
350	date= str( month)+"\|"+str(day)+"\|"+str(year)
351	msg= "Model saved at %s on %s"%(my_time, date)
352	## post help message for the selected model
353	msg +=" Saved! right click on this page to retrieve this model"
354	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
355
356	id = wx.NewId()
357	self.slicerpop.Append(id,name,str(msg))
358	wx.EVT_MENU(self, id, self.onResetModel)
359
360	def onSetFocus(self, evt):
361	"""
362	highlight the current textcrtl and hide the error text control shown
363	after fitting
364	"""
365
366	if hasattr(self,"text2_3"):
367	self.text2_3.Hide()
368	if len(self.parameters)>0:
369	for item in self.parameters:
370	## hide statictext +/-
371	if item[3]!=None:
372	item[3].Hide()
373	## hide textcrtl for error after fit
374	if item[4]!=None:
375	item[4].Clear()
376	item[4].Hide()
377	if len(self.fittable_param)>0:
378	for item in self.fittable_param:
379	## hide statictext +/-
380	if item[3]!=None:
381	item[3].Hide()
382	## hide textcrtl for error after fit
383	if item[4]!=None:
384	item[4].Clear()
385	item[4].Hide()
386	self.Layout()
387	# Get a handle to the TextCtrl
388	widget = evt.GetEventObject()
389	# Select the whole control, after this event resolves
390	wx.CallAfter(widget.SetSelection, -1,-1)
391	return
392
393
394	def read_file(self, path):
395	"""
396	Read two columns file
397	@param path: the path to the file to read
398	"""
399	try:
400	if path==None:
401	wx.PostEvent(self.parent.parent, StatusEvent(status=\
402	" Selected Distribution was not loaded: %s"%path))
403	return None, None
404	input_f = open(path, 'r')
405	buff = input_f.read()
406	lines = buff.split('\n')
407
408	angles = []
409	weights=[]
410	for line in lines:
411	toks = line.split()
412	if len(toks)==2:
413	try:
414	angle = float(toks[0])
415	weight = float(toks[1])
416	except:
417	# Skip non-data lines
418	pass
419	angles.append(angle)
420	weights.append(weight)
421	return numpy.array(angles), numpy.array(weights)
422	except:
423	raise
424
425
426	def createMemento(self):
427	"""
428	return the current state of the page
429	"""
430	return self.state.clone()
431
432
433	def save_current_state(self):
434	"""
435	Store current state
436	"""
437	if self.model!= None:
438	self.state.model = self.model.clone()
439	self.state.save_data(self.data)
440
441	if hasattr(self,"cb1"):
442	self.state.cb1= self.cb1.GetValue()
443
444	if hasattr(self,"enable_disp"):
445	self.state.enable_disp= self.enable_disp.GetValue()
446	self.state.disable_disp = self.disable_disp.GetValue()
447
448	self.state.smearer = copy.deepcopy(self.smearer)
449	if hasattr(self,"enable_smearer"):
450	self.state.enable_smearer = copy.deepcopy(self.enable_smearer.GetValue())
451	self.state.disable_smearer = copy.deepcopy(self.disable_smearer.GetValue())
452
453
454	if hasattr(self,"disp_box"):
455	self.state.disp_box = self.disp_box.GetCurrentSelection()
456	self._save_plotting_range()
457
458	## save checkbutton state and txtcrtl values
459
460	self._copy_parameters_state(self.parameters, self.state.parameters)
461	self._copy_parameters_state(self.fittable_param, self.state.fittable_param)
462	self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
463	self._copy_parameters_state(self.orientation_params,
464	self.state.orientation_params)
465	self._copy_parameters_state(self.orientation_params_disp,
466	self.state.orientation_params_disp)
467
468	## post state to fit panel
469	event = PageInfoEvent(page = self)
470	wx.PostEvent(self.parent, event)
471
472
473	def reset_page_helper(self, state):
474	"""
475	Use page_state and change the state of existing page
476	"""
477	self.state = state.clone()
478	self.model= self.state.model
479	self.data = self.state.data
480	self.smearer= self.state.smearer
481	self.state.enable_smearer=state.enable_smearer
482
483	##model parameter values restore
484	self._set_model_sizer_selection( self.model )
485	self.set_model_param_sizer(self.model)
486	if hasattr(self, "cb1"):
487	self.cb1.SetValue(self.state.cb1)
488	## display dispersion info layer
489	self.enable_disp.SetValue(self.state.enable_disp)
490	self.disable_disp.SetValue(self.state.disable_disp)
491	if hasattr(self, "disp_box"):
492	self.disp_box.SetSelection(self.state.disp_box)
493	self._set_dipers_Param(event=None)
494	##plotting range restore
495	self._reset_plotting_range()
496	## smearing info restore
497	if hasattr(self,"enable_smearer"):
498	## set smearing value whether or not the data contain the smearing info
499	self.enable_smearer.SetValue(state.enable_smearer)
500	self.disable_smearer.SetValue(state.disable_smearer)
501	self.compute_chisqr(smearer= self.smearer)
502
503	## reset context menu items
504	self._reset_context_menu()
505	## reset state of checkbox,textcrtl and parameters value
506	self._reset_parameters_state(self.parameters,state.parameters)
507	self._reset_parameters_state(self.fittable_param,state.fittable_param)
508	self._reset_parameters_state(self.fixed_param,state.fixed_param)
509	self._reset_parameters_state(self.orientation_params_disp,
510	state.orientation_params_disp)
511	self._reset_parameters_state(self.orientation_params,
512	state.orientation_params)
513
514	## draw the model with previous parameters value
515	self._draw_model()
516
517
518
519
520
521	def _selectDlg(self):
522	"""
523	open a dialog file to selected the customized dispersity
524	"""
525	import os
526	dlg = wx.FileDialog(self, "Choose a weight file",
527	self._default_save_location , "", ".", wx.OPEN)
528	path = None
529	if dlg.ShowModal() == wx.ID_OK:
530	path = dlg.GetPath()
531	dlg.Destroy()
532	return path
533
534
535	def _reset_context_menu(self):
536	"""
537	reset the context menu
538	"""
539	for name, state in self.state.saved_states.iteritems():
540	self.number_saved_state += 1
541	## Add item in the context menu
542	id = wx.NewId()
543	self.slicerpop.Append(id,name, 'Save model and state %g'%self.number_saved_state)
544	wx.EVT_MENU(self, id, self.onResetModel)
545
546	def _reset_plotting_range(self):
547	"""
548	Reset the plotting range to a given state
549	"""
550
551	#self.qmin.SetValue(format_number(self.state.qmin))
552	#self.qmax.SetValue(format_number(self.state.qmax))
553	self.qmin.SetValue(str(self.state.qmin))
554	self.qmax.SetValue(str(self.state.qmax))
555	if self.state.npts!=None:
556	self.npts.SetValue(format_number(self.state.npts))
557	self.num_points = float(self.state.npts)
558
559	self.qmin_x = float(self.qmin.GetValue())
560	self.qmax_x = float(self.qmax.GetValue())
561
562
563	def _save_plotting_range(self ):
564	"""
565	save the state of plotting range
566	"""
567	self.state.qmin = self.qmin_x
568	self.state.qmax = self.qmax_x
569	if self.npts!=None:
570	self.state.npts= self.num_points
571
572
573	def _onparamEnter_helper(self):
574	"""
575	check if values entered by the user are changed and valid to replot
576	model
577	use : _check_value_enter
578	"""
579	if self.model !=None:
580	## save current state
581	self.save_current_state()
582	# Flag to register when a parameter has changed.
583	is_modified = False
584	is_modified =self._check_value_enter( self.fittable_param ,is_modified)
585	is_modified =self._check_value_enter( self.fixed_param ,is_modified)
586	is_modified =self._check_value_enter( self.parameters ,is_modified)
587
588
589	self.Layout()
590	# Here we should check whether the boundaries have been modified.
591	# If qmin and qmax have been modified, update qmin and qmax and
592	# set the is_modified flag to True
593	from sans.guiframe.utils import check_value
594	if check_value( self.qmin, self.qmax):
595	if float(self.qmin.GetValue()) != self.qmin_x:
596	self.qmin_x = float(self.qmin.GetValue())
597	is_modified = True
598	if float(self.qmax.GetValue()) != self.qmax_x:
599	self.qmax_x = float(self.qmax.GetValue())
600	is_modified = True
601	self.fitrange = True
602	else:
603	self.fitrange = False
604	if self.npts != None:
605	if check_float(self.npts):
606	if float(self.npts.GetValue()) != self.num_points:
607	self.num_points = float(self.npts.GetValue())
608	is_modified = True
609	else:
610	msg= "Cannot Plot :Must enter a number!!! "
611	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
612
613
614	## if any value is modify draw model with new value
615	if is_modified:
616	self._draw_model()
617
618
619	def _reset_parameters_state(self, listtorestore,statelist):
620	"""
621	Reset the parameters at the given state
622	"""
623	if len(statelist)==0 or len(listtorestore)==0 :
624	return
625	if len(statelist)!= len(listtorestore) :
626	return
627	for j in range(len(listtorestore)):
628	item_page = listtorestore[j]
629	item_page_info = statelist[j]
630	##change the state of the check box for simple parameters
631	if item_page[0]!=None:
632	item_page[0].SetValue(item_page_info[0])
633
634	if item_page[2]!=None:
635	item_page[2].SetValue(item_page_info[2])
636
637	if item_page[3]!=None:
638	## show or hide text +/-
639	if item_page_info[2]:
640	item_page[3].Show(True)
641	else:
642	item_page[3].Hide()
643	if item_page[4]!=None:
644	## show of hide the text crtl for fitting error
645	if item_page_info[4][0]:
646	item_page[4].Show(True)
647	item_page[4].SetValue(item_page_info[4][1])
648	else:
649	item_page[3].Hide()
650	if item_page[5]!=None:
651	## show of hide the text crtl for fitting error
652	if item_page_info[5][0]:
653	item_page[5].Show(True)
654	item_page[5].SetValue(item_page_info[4][1])
655	else:
656	item_page[5].Hide()
657
658	if item_page[6]!=None:
659	## show of hide the text crtl for fitting error
660	if item_page_info[6][0]:
661	item_page[6].Show(True)
662	item_page[6].SetValue(item_page_info[6][1])
663	else:
664	item_page[6].Hide()
665
666
667	def _copy_parameters_state(self, listtocopy, statelist):
668	"""
669	copy the state of button
670	@param listtocopy: the list of check button to copy
671	@param statelist: list of state object to store the current state
672	"""
673	if len(listtocopy)==0:
674	return
675
676	for item in listtocopy:
677	checkbox_state = None
678	if item[0]!= None:
679	checkbox_state= item[0].GetValue()
680	parameter_name = item[1]
681	parameter_value = None
682	if item[2]!=None:
683	parameter_value = item[2].GetValue()
684	static_text = None
685	if item[3]!=None:
686	static_text = item[3].IsShown()
687	error_value = None
688	error_state = None
689	if item[4]!= None:
690	error_value = item[4].GetValue()
691	error_state = item[4].IsShown()
692
693	min_value = None
694	min_state = None
695	if item[5]!= None:
696	min_value = item[5].GetValue()
697	min_state = item[5].IsShown()
698
699	max_value = None
700	max_state = None
701	if item[6]!= None:
702	max_value = item[6].GetValue()
703	max_state = item[6].IsShown()
704
705	statelist.append([checkbox_state, parameter_name, parameter_value,
706	static_text ,[error_state,error_value],
707	[min_state,min_value],[max_state , max_value],None])
708
709
710
711	def _set_model_sizer_selection(self, model):
712	"""
713	Display the sizer according to the type of the current model
714	"""
715	if hasattr(model ,"model2"):
716
717	class_name= model.model2.__class__
718	name= model.model2.name
719	flag= name != "NoStructure"
720	if flag and (class_name in self.model_list_box["Structure Factors"]):
721	self.structurebox.Enable()
722	items = self.structurebox.GetItems()
723	self.sizer1.Layout()
724	self.SetScrollbars(20,20,200,100)
725	for i in range(len(items)):
726	if items[i]== str(name):
727	self.structurebox.SetSelection(i)
728	break
729
730	if hasattr(model ,"model1"):
731	class_name = model.model1.__class__
732	name = model.model1.name
733	self.formfactorbox.Clear()
734
735	for k, list in self.model_list_box.iteritems():
736	if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
737	self.shape_rbutton.SetValue(True)
738	## fill the form factor list with new model
739	self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
740	items = self.formfactorbox.GetItems()
741	## set comboxbox to the selected item
742	for i in range(len(items)):
743	if items[i]== str(name):
744	self.formfactorbox.SetSelection(i)
745	break
746	return
747	elif k == "Shape-Independent":
748	self.shape_indep_rbutton.SetValue(True)
749	elif k == "Structure Factors":
750	self.struct_rbutton.SetValue(True)
751	else:
752	self.plugin_rbutton.SetValue(True)
753
754	if class_name in list:
755	## fill the form factor list with new model
756	self._populate_box(self.formfactorbox, list)
757	items = self.formfactorbox.GetItems()
758	## set comboxbox to the selected item
759	for i in range(len(items)):
760	if items[i]== str(name):
761	self.formfactorbox.SetSelection(i)
762	break
763	break
764	else:
765	## Select the model from the combobox
766	class_name = model.__class__
767	name = model.name
768	self.formfactorbox.Clear()
769	items = self.formfactorbox.GetItems()
770
771	for k, list in self.model_list_box.iteritems():
772	if k in["P(Q)*S(Q)","Shapes" ] and class_name in self.model_list_box["Shapes"]:
773
774	if class_name in self.model_list_box["P(Q)*S(Q)"]:
775	self.structurebox.Enable()
776	else:
777	self.structurebox.Disable()
778	self.structurebox.SetSelection(0)
779
780	self.shape_rbutton.SetValue(True)
781	## fill the form factor list with new model
782	self._populate_box(self.formfactorbox,self.model_list_box["Shapes"])
783	items = self.formfactorbox.GetItems()
784	## set comboxbox to the selected item
785	for i in range(len(items)):
786	if items[i]== str(name):
787	self.formfactorbox.SetSelection(i)
788	break
789	return
790	elif k == "Shape-Independent":
791	self.shape_indep_rbutton.SetValue(True)
792	elif k == "Structure Factors":
793	self.struct_rbutton.SetValue(True)
794	else:
795	self.plugin_rbutton.SetValue(True)
796	if class_name in list:
797	self.structurebox.SetSelection(0)
798	self.structurebox.Disable()
799	## fill the form factor list with new model
800	self._populate_box(self.formfactorbox, list)
801	items = self.formfactorbox.GetItems()
802	## set comboxbox to the selected item
803	for i in range(len(items)):
804	if items[i]== str(name):
805	self.formfactorbox.SetSelection(i)
806	break
807	break
808
809
810	def _draw_model(self):
811	"""
812	Method to draw or refresh a plotted model.
813	The method will use the data member from the model page
814	to build a call to the fitting perspective manager.
815
816	[Note to coder: This way future changes will be done in only one place.]
817	"""
818	if self.model !=None:
819	temp_smear=None
820	if hasattr(self, "enable_smearer"):
821	if self.enable_smearer.GetValue():
822	temp_smear= self.smearer
823	self.manager.draw_model(self.model, data=self.data,
824	smearer= temp_smear,
825	qmin=float(self.qmin_x), qmax=float(self.qmax_x),
826	qstep= float(self.num_points),
827	enable2D=self.enable2D)
828
829	def _set_model_sizer(self, sizer, title="", object=None):
830	"""
831	Use lists to fill a sizer for model info
832	"""
833
834	sizer.Clear(True)
835	box_description= wx.StaticBox(self, -1,str(title))
836	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
837	#--------------------------------------------------------
838	self.shape_rbutton = wx.RadioButton(self, -1, 'Shapes', style=wx.RB_GROUP)
839	self.shape_indep_rbutton = wx.RadioButton(self, -1, "Shape-Independent")
840	self.struct_rbutton = wx.RadioButton(self, -1, "Structure Factor ")
841	self.plugin_rbutton = wx.RadioButton(self, -1, "Customized Models")
842
843	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
844	id= self.shape_rbutton.GetId() )
845	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
846	id= self.shape_indep_rbutton.GetId() )
847	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
848	id= self.struct_rbutton.GetId() )
849	self.Bind( wx.EVT_RADIOBUTTON, self._show_combox,
850	id= self.plugin_rbutton.GetId() )
851
852
853	sizer_radiobutton = wx.GridSizer(2, 2,5, 5)
854	sizer_radiobutton.Add(self.shape_rbutton)
855	sizer_radiobutton.Add(self.shape_indep_rbutton)
856	sizer_radiobutton.Add(self.plugin_rbutton)
857	sizer_radiobutton.Add(self.struct_rbutton)
858
859	sizer_selection = wx.BoxSizer(wx.HORIZONTAL)
860
861	self.text1 = wx.StaticText( self,-1,"" )
862	self.text2 = wx.StaticText( self,-1,"P(Q)*S(Q)" )
863
864
865	self.formfactorbox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
866	if self.model!=None:
867	self.formfactorbox.SetValue(self.model.name)
868
869
870	self.structurebox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
871	wx.EVT_COMBOBOX(self.formfactorbox,-1, self._on_select_model)
872	wx.EVT_COMBOBOX(self.structurebox,-1, self._on_select_model)
873
874
875
876	## fill combox box
877	if len(self.model_list_box)>0:
878	self._populate_box( self.formfactorbox,self.model_list_box["Shapes"])
879
880	if len(self.model_list_box)>0:
881	self._populate_box( self.structurebox,
882	self.model_list_box["Structure Factors"])
883	self.structurebox.Insert("None", 0,None)
884	self.structurebox.SetSelection(0)
885	self.structurebox.Disable()
886
887	if self.model.__class__ in self.model_list_box["P(Q)*S(Q)"]:
888	self.structurebox.Enable()
889
890
891	## check model type to show sizer
892	if self.model !=None:
893	self._set_model_sizer_selection( self.model )
894
895	sizer_selection.Add(self.text1)
896	sizer_selection.Add((5,5))
897	sizer_selection.Add(self.formfactorbox)
898	sizer_selection.Add((5,5))
899	sizer_selection.Add(self.text2)
900	sizer_selection.Add((5,5))
901	sizer_selection.Add(self.structurebox)
902	sizer_selection.Add((5,5))
903
904	boxsizer1.Add( sizer_radiobutton )
905	boxsizer1.Add( (20,20))
906	boxsizer1.Add( sizer_selection )
907	if object !=None:
908	boxsizer1.Add( (-72,-72))
909	boxsizer1.Add( object, 0, wx.ALIGN_RIGHT\| wx.RIGHT, 10)
910	boxsizer1.Add( (60,60))
911	#--------------------------------------------------------
912	sizer.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
913	sizer.Layout()
914	self.SetScrollbars(20,20,200,100)
915
916
917	def _show_combox(self, event):
918	"""
919	Show combox box associate with type of model selected
920	"""
921	## Don't want to populate combo box again if the event comes from check box
922	if self.shape_rbutton.GetValue()and\
923	event.GetEventObject()==self.shape_rbutton:
924	##fill the combobox with form factor list
925	self.structurebox.SetSelection(0)
926	self.structurebox.Disable()
927	self.formfactorbox.Clear()
928	self._populate_box( self.formfactorbox,self.model_list_box["Shapes"])
929
930	if self.shape_indep_rbutton.GetValue()and\
931	event.GetEventObject()==self.shape_indep_rbutton:
932	##fill the combobox with shape independent factor list
933	self.structurebox.SetSelection(0)
934	self.structurebox.Disable()
935	self.formfactorbox.Clear()
936	self._populate_box( self.formfactorbox,
937	self.model_list_box["Shape-Independent"])
938
939	if self.struct_rbutton.GetValue() and\
940	event.GetEventObject()==self.struct_rbutton:
941	##fill the combobox with structure factor list
942	self.structurebox.SetSelection(0)
943	self.structurebox.Disable()
944	self.formfactorbox.Clear()
945	self._populate_box( self.formfactorbox,
946	self.model_list_box["Structure Factors"])
947
948	if self.plugin_rbutton.GetValue()and\
949	event.GetEventObject()==self.plugin_rbutton:
950
951	##fill the combobox with form factor list
952	self.structurebox.Disable()
953	self.formfactorbox.Clear()
954	self._populate_box( self.formfactorbox,
955	self.model_list_box["Customized Models"])
956
957	self._on_select_model(event=None)
958	self.sizer4_4.Layout()
959	self.sizer4.Layout()
960	self.Layout()
961	self.Refresh()
962	self.SetScrollbars(20,20,200,100)
963
964	def _populate_box(self, combobox, list):
965	"""
966	fill combox box with dict item
967	@param list: contains item to fill the combox
968	item must model class
969	"""
970	for models in list:
971	model= models()
972	name = model.__class__.__name__
973	if models.__name__!="NoStructure":
974	if hasattr(model, "name"):
975	name = model.name
976	combobox.Append(name,models)
977	try:
978
979	combobox.SetSelection(0)
980
981	except:
982	pass
983
984	return 0
985
986
987	def _on_select_model_helper(self):
988	"""
989	call back for model selection
990	"""
991	## reset dictionary containing reference to dispersion
992	self._disp_obj_dict = {}
993	self.disp_cb_dict ={}
994	f_id = self.formfactorbox.GetCurrentSelection()
995	form_factor = self.formfactorbox.GetClientData( f_id )
996	if not form_factor in self.model_list_box["multiplication"]:
997	self.structurebox.Disable()
998	self.structurebox.SetSelection(0)
999	else:
1000	self.structurebox.Enable()
1001
1002	s_id = self.structurebox.GetCurrentSelection()
1003	struct_factor = self.structurebox.GetClientData( s_id )
1004
1005	if struct_factor !=None:
1006	from sans.models.MultiplicationModel import MultiplicationModel
1007	self.model= MultiplicationModel(form_factor(),struct_factor())
1008	else:
1009	self.model= form_factor()
1010
1011	## post state to fit panel
1012	self.state.model =self.model
1013	## post state to fit panel
1014	event = PageInfoEvent(page = self)
1015	wx.PostEvent(self.parent, event)
1016
1017	self.sizer4_4.Layout()
1018	self.sizer4.Layout()
1019	self.Layout()
1020	self.SetScrollbars(20,20,200,100)
1021	self.Refresh()
1022
1023
1024
1025	def _onparamEnter(self,event):
1026	"""
1027	when enter value on panel redraw model according to changed
1028	"""
1029	tcrtl= event.GetEventObject()
1030	if check_float(tcrtl):
1031	self._onparamEnter_helper()
1032	else:
1033	msg= "Cannot Plot :Must enter a number!!! "
1034	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1035	return
1036
1037
1038	def _check_value_enter(self, list, modified):
1039	"""
1040	@param list: model parameter and panel info
1041	each item of the list should be as follow:
1042	item=[cb state, name, value, "+/-", error of fit, min, max , units]
1043	"""
1044	is_modified = modified
1045	if len(list)==0:
1046	return is_modified
1047	for item in list:
1048	try:
1049	name = str(item[1])
1050	if hasattr(self,"text2_3"):
1051	self.text2_3.Hide()
1052	## check model parameters range
1053	## check minimun value
1054	param_min= None
1055	param_max= None
1056	if item[5]!= None:
1057	if item[5].GetValue().lstrip().rstrip()!="":
1058	param_min = float(item[5].GetValue())
1059
1060	## check maximum value
1061	if item[6]!= None:
1062	if item[6].GetValue().lstrip().rstrip()!="":
1063	param_max = float(item[6].GetValue())
1064
1065	from sans.guiframe.utils import check_value
1066	if param_min != None and param_max !=None:
1067	if not check_value(item[5], item[6]):
1068	msg= "Wrong Fit range entered for parameter "
1069	msg+= "name %s of model %s "%(name, self.model.name)
1070	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1071	if name in self.model.details.keys():
1072	self.model.details[name][1:]= param_min,param_max
1073
1074	## hide statictext +/-
1075	if item[3]!=None:
1076	item[3].Hide()
1077	## hide textcrtl for error after fit
1078	if item[4]!=None:
1079	item[4].Clear()
1080	item[4].Hide()
1081
1082
1083	value= float(item[2].GetValue())
1084
1085	# If the value of the parameter has changed,
1086	# +update the model and set the is_modified flag
1087	if value != self.model.getParam(name):
1088	self.model.setParam(name,value)
1089	is_modified = True
1090
1091	except:
1092	msg= "Model Drawing Error:wrong value entered : %s"% sys.exc_value
1093	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1094	return
1095
1096	return is_modified
1097
1098
1099	def _set_dipers_Param(self, event):
1100	"""
1101	Add more item to select user dispersity
1102	"""
1103	if self.model == None:
1104	msg= " Select non - model value:%s !"%self.model
1105	wx.PostEvent(self.parent.parent, StatusEvent(status= msg))
1106	return
1107	else:
1108
1109	if self.enable_disp.GetValue():
1110	self.model_disp.Show(True)
1111	self.disp_box.Show(True)
1112	## layout for model containing no dispersity parameters
1113	if len(self.disp_list)==0:
1114	self._layout_sizer_noDipers()
1115	else:
1116	## set gaussian sizer
1117	self._on_select_Disp(event=None)
1118	else:
1119	self.model_disp.Hide()
1120	self.disp_box.Hide()
1121	self.sizer4_4.Clear(True)
1122	self._reset_dispersity()
1123
1124
1125	## post state to fit panel
1126	self.save_current_state()
1127
1128
1129
1130
1131	def _layout_sizer_noDipers(self):
1132	"""
1133	Draw a sizer with no dispersity info
1134	"""
1135	ix=0
1136	iy=1
1137	self.fittable_param=[]
1138	self.fixed_param=[]
1139	self.model_disp.Hide()
1140	self.disp_box.Hide()
1141	self.sizer4_4.Clear(True)
1142	model_disp = wx.StaticText(self, -1, 'No PolyDispersity for this model')
1143	self.sizer4_4.Add(model_disp,( iy, ix),(1,1), wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
1144	self.sizer4_4.Layout()
1145	self.sizer4.Layout()
1146	self.SetScrollbars(20,20,200,100)
1147
1148
1149	def _reset_dispersity(self):
1150	"""
1151	put gaussian dispersity into current model
1152	"""
1153	self.fittable_param=[]
1154	self.fixed_param=[]
1155
1156	from sans.models.dispersion_models import GaussianDispersion
1157	if len(self.disp_cb_dict)==0:
1158	self.sizer4_4.Clear(True)
1159	self.sizer4_4.Layout()
1160	self.sizer4.Layout()
1161	self.Layout()
1162	self.Refresh()
1163	self.SetScrollbars(20,20,200,100)
1164	return
1165
1166	for p in self.disp_cb_dict:
1167	# The parameter was un-selected. Go back to Gaussian model (with 0 pts)
1168	disp_model = GaussianDispersion()
1169	# Store the object to make it persist outside the scope of this method
1170	#TODO: refactor model to clean this up?
1171	self._disp_obj_dict[p] = disp_model
1172	# Set the new model as the dispersion object for the selected parameter
1173	self.model.set_dispersion(p, disp_model)
1174	# Redraw the model
1175	self._draw_model()
1176
1177	self.sizer4_4.Layout()
1178	self.sizer4.Layout()
1179	self.Layout()
1180	self.SetScrollbars(20,20,200,100)
1181	self.Refresh()
1182
1183
1184
1185	def _on_select_Disp(self,event):
1186	"""
1187	allow selecting different dispersion
1188	self.disp_list should change type later .now only gaussian
1189	"""
1190
1191	n = self.disp_box.GetCurrentSelection()
1192	dispersity= self.disp_box.GetClientData(n)
1193	name= dispersity.__name__
1194	if name == "GaussianDispersion":
1195	self._set_sizer_gaussian()
1196
1197	if name=="ArrayDispersion":
1198	self._set_sizer_arraydispersion()
1199
1200	self.state.disp_box= n
1201	## post state to fit panel
1202	event = PageInfoEvent(page = self)
1203	wx.PostEvent(self.parent, event)
1204
1205	self.sizer4_4.Layout()
1206	self.sizer4.Layout()
1207	self.Layout()
1208	self.SetScrollbars(20,20,200,100)
1209	self.Refresh()
1210
1211
1212
1213
1214	def _set_sizer_arraydispersion(self):
1215	"""
1216	draw sizer with array dispersity parameters
1217	"""
1218	self.orientation_params_disp=[]
1219	self.sizer4_4.Clear(True)
1220	ix=0
1221	iy=1
1222	disp1 = wx.StaticText(self, -1, 'Array Dispersion')
1223	self.sizer4_4.Add(disp1,( iy, ix),(1,1), wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
1224
1225	# Look for model parameters to which we can apply an ArrayDispersion model
1226	# Add a check box for each parameter.
1227	self.disp_cb_dict = {}
1228	for p in self.model.dispersion.keys():
1229	ix+=1
1230	self.disp_cb_dict[p] = wx.CheckBox(self, -1, p, (10, 10))
1231
1232	wx.EVT_CHECKBOX(self, self.disp_cb_dict[p].GetId(), self.select_disp_angle)
1233	self.sizer4_4.Add(self.disp_cb_dict[p], (iy, ix), (1,1), wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
1234
1235	ix =0
1236	iy +=1
1237	self.sizer4_4.Add((20,20),(iy,ix),(1,1), wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
1238	self.sizer4_4.Layout()
1239	self.sizer4.Layout()
1240	self.SetScrollbars(20,20,200,100)
1241
1242
1243	def _set_range_sizer(self, title, object1=None,object=None):
1244	"""
1245	Fill the
1246	"""
1247	self.sizer5.Clear(True)
1248	box_description= wx.StaticBox(self, -1,str(title))
1249	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
1250	#--------------------------------------------------------------
1251	self.qmin = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
1252	#For qmin and qmax, do not use format_number.(If do, qmin and max could be different from what is in the data.)
1253	#self.qmin.SetValue(format_number(self.qmin_x))
1254	self.qmin.SetValue(str(self.qmin_x))
1255	self.qmin.SetToolTipString("Minimun value of Q in linear scale.")
1256	self.qmin.Bind(wx.EVT_SET_FOCUS, self.onSetFocus)
1257	self.qmin.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
1258	self.qmin.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
1259
1260	self.qmax = wx.TextCtrl(self, -1,size=(_BOX_WIDTH,20))
1261	#For qmin and qmax, do not use format_number.(If do, qmin and max could be different from what is in the data.)
1262	#self.qmax.SetValue(format_number(self.qmax_x))
1263	self.qmax.SetValue(str(self.qmax_x))
1264	self.qmax.SetToolTipString("Maximum value of Q in linear scale.")
1265	self.qmax.Bind(wx.EVT_SET_FOCUS, self.onSetFocus)
1266	self.qmax.Bind(wx.EVT_KILL_FOCUS, self._onparamEnter)
1267	self.qmax.Bind(wx.EVT_TEXT_ENTER, self._onparamEnter)
1268
1269	sizer_horizontal=wx.BoxSizer(wx.HORIZONTAL)
1270	sizer= wx.GridSizer(3, 3,5, 5)
1271
1272	sizer.Add((5,5))
1273	sizer.Add(wx.StaticText(self, -1, 'Min'))
1274	sizer.Add(wx.StaticText(self, -1, 'Max'))
1275	sizer.Add(wx.StaticText(self, -1, 'Q range'))
1276
1277	sizer.Add(self.qmin)
1278	sizer.Add(self.qmax)
1279	sizer_horizontal.Add(sizer)
1280	if object!=None:
1281	sizer_horizontal.Add(object)
1282
1283	if object1!=None:
1284	boxsizer1.Add(object1)
1285	boxsizer1.Add((10,10))
1286	boxsizer1.Add(sizer_horizontal)
1287	## save state
1288	self.save_current_state()
1289	#----------------------------------------------------------------
1290	self.sizer5.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
1291	self.sizer5.Layout()
1292	self.Layout()
1293	self.SetScrollbars(20,20,200,100)
1294
1295
1296	def _fill_save_sizer(self):
1297	"""
1298	Draw the layout for saving option
1299	"""
1300	self.sizer6.Clear(True)
1301	box_description= wx.StaticBox(self, -1,"Save Status")
1302	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
1303	sizer_save = wx.BoxSizer(wx.HORIZONTAL)
1304
1305	self.btSave_title = wx.StaticText(self, -1, 'Save the current panel status')
1306	self.btSave = wx.Button(self,wx.NewId(),'Save')
1307	self.btSave.Bind(wx.EVT_BUTTON, self.onSave,id= self.btSave.GetId())
1308	self.btSave.SetToolTipString("Save current panel status")
1309
1310	sizer_save.Add(self.btSave_title)
1311	sizer_save.Add((20,20),0, wx.LEFT\|wx.RIGHT\|wx.EXPAND,45)
1312
1313	sizer_save.Add(self.btSave)
1314
1315	boxsizer1.Add(sizer_save)
1316	self.sizer6.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
1317	self.sizer6.Layout()
1318	self.SetScrollbars(20,20,200,100)
1319
1320
1321

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

SasView

source: sasview/sansview/perspectives/fitting/basepage.py @ fc6ea43

Download in other formats: