fitpage.py @ 75daa79

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 75daa79 was 75daa79, checked in by Gervaise Alina <gervyh@…>, 13 years ago
refactore compute data range
Property mode set to `100644`
File size: 121.9 KB

Line
1
2
3	import sys
4	import wx
5	import wx.lib.newevent
6	import numpy
7	import copy
8	import math
9	import time
10	from sans.models.dispersion_models import ArrayDispersion, GaussianDispersion
11	from DataLoader.data_info import Data1D
12	from sans.guiframe.events import StatusEvent
13	from sans.guiframe.events import NewPlotEvent
14	from sans.guiframe.utils import format_number,check_float
15
16	(Chi2UpdateEvent, EVT_CHI2_UPDATE) = wx.lib.newevent.NewEvent()
17	_BOX_WIDTH = 76
18	_DATA_BOX_WIDTH = 300
19	SMEAR_SIZE_L = 0.00
20	SMEAR_SIZE_H = 0.00
21
22	import basepage
23	from basepage import BasicPage
24	from basepage import PageInfoEvent
25	from sans.models.qsmearing import smear_selection
26
27
28	class FitPage(BasicPage):
29	"""
30	FitPanel class contains fields allowing to display results when
31	fitting a model and one data
32
33	:note: For Fit to be performed the user should check at least one parameter
34	on fit Panel window.
35	"""
36
37	def __init__(self,parent, color='rand'):
38	"""
39	Initialization of the Panel
40	"""
41	BasicPage.__init__(self, parent, color=color)
42
43	## draw sizer
44	self._fill_data_sizer()
45	self.is_2D = None
46	self.fit_started = False
47	# get smear info from data
48	self._get_smear_info()
49	self._fill_model_sizer( self.sizer1)
50	self._get_defult_custom_smear()
51	self._fill_range_sizer()
52	self._set_smear(self.data)
53	self.Bind(EVT_CHI2_UPDATE, self.on_complete_chisqr)
54	# bind key event
55	self.Bind(wx.EVT_RIGHT_DOWN, self.on_right_down)
56	self._set_bookmark_flag(False)
57	self._set_save_flag(False)
58	self._set_preview_flag(False)
59	self._set_copy_flag(False)
60	self._set_paste_flag(False)
61	self.btFit.SetFocus()
62	self.fill_data_combobox(data_list=self.data_list)
63
64
65	def _fill_data_sizer(self):
66	"""
67	fill sizer 0 with data info
68	"""
69	box_description= wx.StaticBox(self, -1, 'I(q) Data Source')
70	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
71	#----------------------------------------------------------
72	sizer_data = wx.BoxSizer(wx.HORIZONTAL)
73	self.dataSource = wx.ComboBox(self, -1, style=wx.CB_READONLY)
74	wx.EVT_COMBOBOX(self.dataSource, -1, self.on_select_data)
75	self.dataSource.SetMinSize((_DATA_BOX_WIDTH, -1))
76	sizer_data.Add(wx.StaticText(self, -1, 'Name : '))
77	sizer_data.Add(self.dataSource)
78	sizer_data.Add( (0,5) )
79	boxsizer1.Add(sizer_data,0, wx.ALL, 10)
80	self.sizer0.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
81	self.sizer0.Layout()
82
83	def enable_datasource(self):
84	"""
85	Enable or disable data source control depending on existing data
86	"""
87	if not self.data_list:
88	self.dataSource.Disable()
89	else:
90	self.dataSource.Enable()
91
92	def fill_data_combobox(self, data_list):
93	"""
94	Get a list of data and fill the corresponding combobox
95	"""
96	self.dataSource.Clear()
97	self.data_list = data_list
98	self.enable_datasource()
99	for data in self.data_list:
100	if data is not None:
101	self.dataSource.Append(str(data.name), clientData=data)
102	self.dataSource.SetSelection(0)
103	self.on_select_data(event=None)
104
105	def on_select_data(self, event=None):
106	"""
107	"""
108	if event is None and self.dataSource.GetCount() > 0:
109	data = self.dataSource.GetClientData(0)
110	self.set_data(data)
111	elif self.dataSource.GetCount() > 0:
112	pos = self.dataSource.GetSelection()
113	data = self.dataSource.GetClientData(pos)
114	self.set_data(data)
115
116
117
118	def _on_fit_complete(self):
119	"""
120	When fit is complete ,reset the fit button label.
121	"""
122	self.btFit.SetLabel("Fit")
123	self.bind_fit_button()
124
125	def _is_2D(self):
126	"""
127	Check if data_name is Data2D
128
129	:return: True or False
130
131	"""
132	if self.data.__class__.__name__ == "Data2D" or \
133	self.enable2D:
134	return True
135	return False
136
137	def _on_engine_change(self, name):
138	"""
139	get the current name of the fit engine type
140	and update the panel accordingly
141	"""
142
143	self.engine_type = str(name)
144	self.state.engine_type = self.engine_type
145	if not self.is_mac:
146	if len(self.parameters) == 0:
147	self.Layout()
148	return
149	self.Layout()
150	self.Refresh()
151
152
153
154	def _fill_range_sizer(self):
155	"""
156	Fill the sizer containing the plotting range
157	add access to npts
158	"""
159	is_2Ddata = False
160
161	# Check if data is 2D
162	if self.data.__class__.__name__ == "Data2D" or \
163	self.enable2D:
164	is_2Ddata = True
165
166	title = "Fitting"
167	#smear messages & titles
168	smear_message_none = "No smearing is selected..."
169	smear_message_dqdata = "The dQ data is being used for smearing..."
170	smear_message_2d = "Higher accuracy is very time-expensive. Use it with care..."
171	smear_message_new_ssmear = "Please enter only the value of interest to customize smearing..."
172	smear_message_new_psmear = "Please enter both; the dQ will be generated by interpolation..."
173	smear_message_2d_x_title = "<dQp>[1/A]:"
174	smear_message_2d_y_title = "<dQs>[1/A]:"
175	smear_message_pinhole_min_title = "dQ_low[1/A]:"
176	smear_message_pinhole_max_title = "dQ_high[1/A]:"
177	smear_message_slit_height_title = "Slit height[1/A]:"
178	smear_message_slit_width_title = "Slit width[1/A]:"
179
180	self._get_smear_info()
181
182	#Sizers
183	box_description_range = wx.StaticBox(self, -1,str(title))
184	boxsizer_range = wx.StaticBoxSizer(box_description_range, wx.VERTICAL)
185	self.sizer_set_smearer = wx.BoxSizer(wx.VERTICAL)
186	sizer_smearer = wx.BoxSizer(wx.HORIZONTAL)
187	self.sizer_new_smear= wx.BoxSizer(wx.HORIZONTAL)
188	self.sizer_set_masking = wx.BoxSizer(wx.HORIZONTAL)
189	sizer_chi2 = wx.BoxSizer(wx.VERTICAL)
190	smear_set_box= wx.StaticBox(self, -1,'Set Instrumental Smearing')
191	sizer_smearer_box = wx.StaticBoxSizer(smear_set_box, wx.HORIZONTAL)
192	sizer_smearer_box.SetMinSize((_DATA_BOX_WIDTH,85))
193	sizer_fit = wx.GridSizer(2, 4, 2, 6)
194
195	# combobox for smear2d accuracy selection
196	self.smear_accuracy = wx.ComboBox(self, -1,size=(50,-1),style=wx.CB_READONLY)
197	self._set_accuracy_list()
198	self.smear_accuracy.SetValue(self.smear2d_accuracy)
199	self.smear_accuracy.SetSelection(0)
200	self.smear_accuracy.SetToolTipString("'Higher' uses more Gaussian points for smearing computation.")
201
202	wx.EVT_COMBOBOX(self.smear_accuracy,-1, self._on_select_accuracy)
203
204	#Fit button
205	self.btFit = wx.Button(self,wx.NewId(),'Fit', size=(88,25))
206	self.default_bt_colour = self.btFit.GetDefaultAttributes()
207	self.btFit.Bind(wx.EVT_BUTTON, self._onFit,id= self.btFit.GetId())
208	self.btFit.SetToolTipString("Start fitting.")
209
210	#textcntrl for custom resolution
211	self.smear_pinhole_max = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
212	text_enter_callback = self.onPinholeSmear)
213	self.smear_pinhole_min = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
214	text_enter_callback = self.onPinholeSmear)
215	self.smear_slit_height= self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
216	text_enter_callback = self.onSlitSmear)
217	self.smear_slit_width = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH-25,20),style=wx.TE_PROCESS_ENTER,
218	text_enter_callback = self.onSlitSmear)
219
220	## smear
221	self.smear_data_left= BGTextCtrl(self, -1, size=(_BOX_WIDTH-25,20), style=0)
222	self.smear_data_left.SetValue(str(self.dq_l))
223	self.smear_data_right = BGTextCtrl(self, -1, size=(_BOX_WIDTH-25,20), style=0)
224	self.smear_data_right.SetValue(str(self.dq_r))
225
226	#set default values for smear
227	self.smear_pinhole_max.SetValue(str(self.dx_max))
228	self.smear_pinhole_min.SetValue(str(self.dx_min))
229	self.smear_slit_height.SetValue(str(self.dxl))
230	self.smear_slit_width.SetValue(str(self.dxw))
231
232	#Filling the sizer containing instruments smearing info.
233	self.disable_smearer = wx.RadioButton(self, -1, 'None', style=wx.RB_GROUP)
234	self.enable_smearer = wx.RadioButton(self, -1, 'Use dQ Data')
235	#self.enable_smearer.SetToolTipString("Click to use the loaded dQ data for smearing.")
236	self.pinhole_smearer = wx.RadioButton(self, -1, 'Custom Pinhole Smear')
237	#self.pinhole_smearer.SetToolTipString("Click to input custom resolution for pinhole smearing.")
238	self.slit_smearer = wx.RadioButton(self, -1, 'Custom Slit Smear')
239	#self.slit_smearer.SetToolTipString("Click to input custom resolution for slit smearing.")
240	self.Bind(wx.EVT_RADIOBUTTON, self.onSmear, id=self.disable_smearer.GetId())
241	self.Bind(wx.EVT_RADIOBUTTON, self.onSmear, id=self.enable_smearer.GetId())
242	self.Bind(wx.EVT_RADIOBUTTON, self.onPinholeSmear, id=self.pinhole_smearer.GetId())
243	self.Bind(wx.EVT_RADIOBUTTON, self.onSlitSmear, id=self.slit_smearer.GetId())
244	self.disable_smearer.SetValue(True)
245
246	# add 4 types of smearing to the sizer
247	sizer_smearer.Add( self.disable_smearer,0, wx.LEFT, 10)
248	sizer_smearer.Add((10,10))
249	sizer_smearer.Add( self.enable_smearer)
250	sizer_smearer.Add((10,10))
251	sizer_smearer.Add( self.pinhole_smearer )
252	sizer_smearer.Add((10,10))
253	sizer_smearer.Add( self.slit_smearer )
254	sizer_smearer.Add((10,10))
255
256	# StaticText for chi2, N(for fitting), Npts
257	self.tcChi = BGTextCtrl(self, -1, "-", size=(75,20), style=0)
258	self.tcChi.SetToolTipString("Chi2/Npts(Fit)")
259	self.Npts_fit = BGTextCtrl(self, -1, "-", size=(75,20), style=0)
260	self.Npts_fit.SetToolTipString(\
261	" Npts : number of points selected for fitting")
262	self.Npts_total = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH, 20),
263	style=wx.TE_PROCESS_ENTER,
264	text_enter_callback=self._onQrangeEnter)
265	self.Npts_total.SetValue(format_number(self.npts_x))
266	self.Npts_total.SetToolTipString(\
267	" Total Npts : total number of data points")
268
269	# Update and Draw button
270	self.draw_button = wx.Button(self,wx.NewId(),'Compute', size=(88,24))
271	self.draw_button.Bind(wx.EVT_BUTTON, \
272	self._onDraw,id= self.draw_button.GetId())
273	self.draw_button.SetToolTipString("Compute and Draw.")
274
275	box_description_1= wx.StaticText(self, -1,' Chi2/Npts')
276	box_description_2= wx.StaticText(self, -1,'Npts(Fit)')
277	box_description_3= wx.StaticText(self, -1,'Total Npts')
278	box_description_3.SetToolTipString( \
279	" Total Npts : total number of data points")
280	#box_description_4= wx.StaticText(self, -1,' ')
281
282
283	sizer_fit.Add(box_description_1,0,0)
284	sizer_fit.Add(box_description_2,0,0)
285	sizer_fit.Add(box_description_3,0,0)
286	sizer_fit.Add(self.draw_button,0,0)
287	sizer_fit.Add(self.tcChi,0,0)
288	sizer_fit.Add(self.Npts_fit ,0,0)
289	sizer_fit.Add(self.Npts_total,0,0)
290	sizer_fit.Add(self.btFit,0,0)
291
292	# StaticText for smear
293	#self.tcChi = wx.StaticText(self, -1, "-", style=wx.ALIGN_LEFT)
294	self.smear_description_none = wx.StaticText(self, -1, smear_message_none , style=wx.ALIGN_LEFT)
295	self.smear_description_dqdata = wx.StaticText(self, -1, smear_message_dqdata , style=wx.ALIGN_LEFT)
296	self.smear_description_type = wx.StaticText(self, -1, "Type:" , style=wx.ALIGN_LEFT)
297	self.smear_description_accuracy_type = wx.StaticText(self, -1, "Accuracy:" , style=wx.ALIGN_LEFT)
298	self.smear_description_smear_type = BGTextCtrl(self, -1, size=(57,20), style=0)
299	self.smear_description_smear_type.SetValue(str(self.dq_l))
300	self.SetBackgroundColour(self.GetParent().GetBackgroundColour())
301	self.smear_description_2d = wx.StaticText(self, -1, smear_message_2d , style=wx.ALIGN_LEFT)
302	self.smear_message_new_s = wx.StaticText(self, -1, smear_message_new_ssmear , style=wx.ALIGN_LEFT)
303	self.smear_message_new_p = wx.StaticText(self, -1, smear_message_new_psmear , style=wx.ALIGN_LEFT)
304	self.smear_description_2d_x = wx.StaticText(self, -1, smear_message_2d_x_title , style=wx.ALIGN_LEFT)
305	self.smear_description_2d_x.SetToolTipString(" dQp(parallel) in q_r direction.")
306	self.smear_description_2d_y = wx.StaticText(self, -1, smear_message_2d_y_title , style=wx.ALIGN_LEFT)
307	self.smear_description_2d_y.SetToolTipString(" dQs(perpendicular) in q_phi direction.")
308	self.smear_description_pin_min = wx.StaticText(self, -1, smear_message_pinhole_min_title , style=wx.ALIGN_LEFT)
309	self.smear_description_pin_max = wx.StaticText(self, -1, smear_message_pinhole_max_title , style=wx.ALIGN_LEFT)
310	self.smear_description_slit_height = wx.StaticText(self, -1, smear_message_slit_height_title , style=wx.ALIGN_LEFT)
311	self.smear_description_slit_width = wx.StaticText(self, -1, smear_message_slit_width_title , style=wx.ALIGN_LEFT)
312
313	#arrange sizers
314	#boxsizer1.Add( self.tcChi )
315	self.sizer_set_smearer.Add(sizer_smearer )
316	self.sizer_set_smearer.Add((10,10))
317	self.sizer_set_smearer.Add( self.smear_description_none,0, wx.CENTER, 10 )
318	self.sizer_set_smearer.Add( self.smear_description_dqdata,0, wx.CENTER, 10 )
319	self.sizer_set_smearer.Add( self.smear_description_2d,0, wx.CENTER, 10 )
320	self.sizer_new_smear.Add( self.smear_description_type,0, wx.CENTER, 10 )
321	self.sizer_new_smear.Add( self.smear_description_accuracy_type,0, wx.CENTER, 10 )
322	self.sizer_new_smear.Add( self.smear_accuracy )
323	self.sizer_new_smear.Add( self.smear_description_smear_type,0, wx.CENTER, 10 )
324	self.sizer_new_smear.Add((15,-1))
325	self.sizer_new_smear.Add( self.smear_description_2d_x,0, wx.CENTER, 10 )
326	self.sizer_new_smear.Add( self.smear_description_pin_min,0, wx.CENTER, 10 )
327	self.sizer_new_smear.Add( self.smear_description_slit_height,0, wx.CENTER, 10 )
328
329	self.sizer_new_smear.Add( self.smear_pinhole_min,0, wx.CENTER, 10 )
330	self.sizer_new_smear.Add( self.smear_slit_height,0, wx.CENTER, 10 )
331	self.sizer_new_smear.Add( self.smear_data_left,0, wx.CENTER, 10 )
332	self.sizer_new_smear.Add((20,-1))
333	self.sizer_new_smear.Add( self.smear_description_2d_y,0, wx.CENTER, 10 )
334	self.sizer_new_smear.Add( self.smear_description_pin_max,0, wx.CENTER, 10 )
335	self.sizer_new_smear.Add( self.smear_description_slit_width,0, wx.CENTER, 10 )
336
337	self.sizer_new_smear.Add( self.smear_pinhole_max,0, wx.CENTER, 10 )
338	self.sizer_new_smear.Add( self.smear_slit_width,0, wx.CENTER, 10 )
339	self.sizer_new_smear.Add( self.smear_data_right,0, wx.CENTER, 10 )
340
341	self.sizer_set_smearer.Add( self.smear_message_new_s,0, wx.CENTER, 10)
342	self.sizer_set_smearer.Add( self.smear_message_new_p,0, wx.CENTER, 10)
343	self.sizer_set_smearer.Add((5,2))
344	self.sizer_set_smearer.Add( self.sizer_new_smear,0, wx.CENTER, 10 )
345
346	# add all to chi2 sizer
347	sizer_smearer_box.Add(self.sizer_set_smearer)
348	sizer_chi2.Add(sizer_smearer_box)
349	sizer_chi2.Add((-1,5))
350
351	# hide all smear messages and textctrl
352	self._hide_all_smear_info()
353
354	# get smear_selection
355	self.current_smearer= smear_selection( self.data, self.model )
356
357	# Show only the relevant smear messages, etc
358	if self.current_smearer == None:
359	if not is_2Ddata:
360	self.smear_description_none.Show(True)
361	self.enable_smearer.Disable()
362	else:
363	self.smear_description_none.Show(True)
364	#self.smear_description_2d.Show(True)
365	#self.pinhole_smearer.Disable()
366	self.slit_smearer.Disable()
367	#self.enable_smearer.Disable()
368	if self.data == None:
369	self.slit_smearer.Disable()
370	self.pinhole_smearer.Disable()
371	self.enable_smearer.Disable()
372	else: self._show_smear_sizer()
373	boxsizer_range.Add(self.sizer_set_masking)
374	#2D data? default
375	is_2Ddata = False
376
377	#check if it is 2D data
378	if self.data.__class__.__name__ == "Data2D" or \
379	self.enable2D:
380	is_2Ddata = True
381
382	self.sizer5.Clear(True)
383
384	self.qmin = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH,20),
385	style=wx.TE_PROCESS_ENTER,
386	text_enter_callback = self._onQrangeEnter)
387	self.qmin.SetValue(str(self.qmin_x))
388	self.qmin.SetToolTipString("Minimun value of Q in linear scale.")
389
390	self.qmax = self.ModelTextCtrl(self, -1,size=(_BOX_WIDTH,20),
391	style=wx.TE_PROCESS_ENTER,
392	text_enter_callback=self._onQrangeEnter)
393	self.qmax.SetValue(str(self.qmax_x))
394	self.qmax.SetToolTipString("Maximum value of Q in linear scale.")
395	"""
396	self.theory_npts_tcrtl = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH, 20),
397	style=wx.TE_PROCESS_ENTER,
398	text_enter_callback=self._onQrangeEnter)
399	self.theory_npts_tcrtl.SetValue(format_number(self.npts_x))
400	self.theory_npts_tcrtl.SetToolTipString("Number of point to plot.")
401	"""
402	id = wx.NewId()
403	self.reset_qrange =wx.Button(self,id,'Reset',size=(77,20))
404
405	self.reset_qrange.Bind(wx.EVT_BUTTON, self.on_reset_clicked,id=id)
406	self.reset_qrange.SetToolTipString("Reset Q range to the default values")
407
408	sizer_horizontal=wx.BoxSizer(wx.HORIZONTAL)
409	sizer= wx.GridSizer(2, 4, 2, 6)
410
411	self.btEditMask = wx.Button(self,wx.NewId(),'Editor', size=(88,23))
412	self.btEditMask.Bind(wx.EVT_BUTTON, self._onMask,id= self.btEditMask.GetId())
413	self.btEditMask.SetToolTipString("Edit Mask.")
414	self.EditMask_title = wx.StaticText(self, -1, ' Masking(2D)')
415
416	sizer.Add(wx.StaticText(self, -1, 'Q range'))
417	sizer.Add(wx.StaticText(self, -1, ' Min[1/A]'))
418	sizer.Add(wx.StaticText(self, -1, ' Max[1/A]'))
419	sizer.Add(self.EditMask_title)
420	#sizer.Add(wx.StaticText(self, -1, ''))
421	sizer.Add(self.reset_qrange)
422	sizer.Add(self.qmin)
423	sizer.Add(self.qmax)
424	#sizer.Add(self.theory_npts_tcrtl)
425	sizer.Add(self.btEditMask)
426	boxsizer_range.Add(sizer_chi2)
427	boxsizer_range.Add((10,10))
428	boxsizer_range.Add(sizer)
429
430	boxsizer_range.Add((10,15))
431	boxsizer_range.Add(sizer_fit)
432	if is_2Ddata:
433	self.btEditMask.Enable()
434	self.EditMask_title.Enable()
435	else:
436	self.btEditMask.Disable()
437	self.EditMask_title.Disable()
438	## save state
439	self.save_current_state()
440	self.sizer5.Add(boxsizer_range,0, wx.EXPAND \| wx.ALL, 10)
441	self.sizer5.Layout()
442
443
444	def _fill_model_sizer(self, sizer):
445	"""
446	fill sizer containing model info
447	"""
448	##Add model function Details button in fitpanel.
449	##The following 3 lines are for Mac. Let JHC know before modifying...
450	title = "Model"
451	self.formfactorbox = None
452	self.multifactorbox = None
453	self.mbox_description= wx.StaticBox(self, -1,str(title))
454	boxsizer1 = wx.StaticBoxSizer(self.mbox_description, wx.VERTICAL)
455
456	id = wx.NewId()
457	self.model_help =wx.Button(self,id,'Details', size=(80,23))
458	self.model_help.Bind(wx.EVT_BUTTON, self.on_model_help_clicked,id=id)
459	self.model_help.SetToolTipString("Model Function Help")
460	id = wx.NewId()
461	self.model_view = wx.Button(self, id,"1D Mode", size=(80, 23))
462	self.model_view.Bind(wx.EVT_BUTTON, self._onModel2D, id=id)
463	hint = "toggle view of model from 1D to 2D or 2D to 1D"
464	self.model_view.SetToolTipString(hint)
465
466	self.shape_rbutton = wx.RadioButton(self, -1, 'Shapes', style=wx.RB_GROUP)
467	self.shape_indep_rbutton = wx.RadioButton(self, -1, "Shape-Independent")
468	self.struct_rbutton = wx.RadioButton(self, -1, "Structure Factor ")
469	self.plugin_rbutton = wx.RadioButton(self, -1, "Customized Models")
470
471	self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
472	id= self.shape_rbutton.GetId())
473	self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
474	id= self.shape_indep_rbutton.GetId())
475	self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
476	id= self.struct_rbutton.GetId())
477	self.Bind(wx.EVT_RADIOBUTTON, self._show_combox,
478	id= self.plugin_rbutton.GetId())
479	#MAC needs SetValue
480	self.shape_rbutton.SetValue(True)
481
482	sizer_radiobutton = wx.GridSizer(2, 3, 5, 5)
483	sizer_radiobutton.Add(self.shape_rbutton)
484	sizer_radiobutton.Add(self.shape_indep_rbutton)
485	#sizer_radiobutton.Add((5, 5))
486	sizer_radiobutton.Add(self.model_view,1, wx.LEFT, 20)
487	sizer_radiobutton.Add(self.plugin_rbutton)
488	sizer_radiobutton.Add(self.struct_rbutton)
489	#sizer_radiobutton.Add((5, 5))
490	sizer_radiobutton.Add(self.model_help,1, wx.LEFT, 20)
491
492	sizer_selection = wx.BoxSizer(wx.HORIZONTAL)
493	mutifactor_selection = wx.BoxSizer(wx.HORIZONTAL)
494
495	self.text1 = wx.StaticText(self,-1,"" )
496	self.text2 = wx.StaticText(self,-1,"P(Q)*S(Q)" )
497	self.mutifactor_text = wx.StaticText( self,-1,"No. of Shells: ")
498	self.mutifactor_text1 = wx.StaticText( self,-1,"" )
499	self.show_sld_button = wx.Button( self,-1,"Show SLD Profile" )
500	self.show_sld_button.Bind(wx.EVT_BUTTON,self._on_show_sld)
501
502	self.formfactorbox = wx.ComboBox(self, -1,style=wx.CB_READONLY)
503	if self.model!= None:
504	self.formfactorbox.SetValue(self.model.name)
505	self.structurebox = wx.ComboBox(self, -1, style=wx.CB_READONLY)
506	self.multifactorbox = wx.ComboBox(self, -1, style=wx.CB_READONLY)
507	self.initialize_combox()
508	wx.EVT_COMBOBOX(self.formfactorbox, -1, self._on_select_model)
509
510	wx.EVT_COMBOBOX(self.structurebox, -1, self._on_select_model)
511	wx.EVT_COMBOBOX(self.multifactorbox, -1, self._on_select_model)
512	## check model type to show sizer
513	if self.model !=None:
514	self._set_model_sizer_selection(self.model)
515
516	sizer_selection.Add(self.text1)
517	sizer_selection.Add((5, 5))
518	sizer_selection.Add(self.formfactorbox)
519	sizer_selection.Add((5, 5))
520	sizer_selection.Add(self.text2)
521	sizer_selection.Add((5, 5))
522	sizer_selection.Add(self.structurebox)
523
524	mutifactor_selection.Add((10,5))
525	mutifactor_selection.Add(self.mutifactor_text)
526	mutifactor_selection.Add(self.multifactorbox)
527	mutifactor_selection.Add((5, 5))
528	mutifactor_selection.Add(self.mutifactor_text1)
529	mutifactor_selection.Add((10, 5))
530	mutifactor_selection.Add(self.show_sld_button)
531
532
533	boxsizer1.Add(sizer_radiobutton)
534	boxsizer1.Add((10, 10))
535	boxsizer1.Add(sizer_selection)
536	boxsizer1.Add((10, 10))
537	boxsizer1.Add(mutifactor_selection)
538
539	self._set_multfactor_combobox()
540	self.multifactorbox.SetSelection(1)
541	self.show_sld_button.Hide()
542	sizer.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
543	sizer.Layout()
544
545	def _set_sizer_dispersion(self):
546	"""
547	draw sizer with gaussian dispersity parameters
548	"""
549	self.fittable_param=[]
550	self.fixed_param=[]
551	self.orientation_params_disp=[]
552
553	self.sizer4_4.Clear(True)
554	if self.model==None:
555	##no model is selected
556	return
557	if not self.enable_disp.GetValue():
558	## the user didn't select dispersity display
559	return
560
561	self._reset_dispersity()
562
563	## fill a sizer with the combobox to select dispersion type
564	#sizer_select_dispers = wx.BoxSizer(wx.HORIZONTAL)
565	model_disp = wx.StaticText(self, -1, 'Function')
566
567	import sans.models.dispersion_models
568	self.polydisp= sans.models.dispersion_models.models
569
570	ix = 0
571	iy = 0
572	disp = wx.StaticText(self, -1, ' ')
573	self.sizer4_4.Add(disp,( iy, ix),(1,1),
574	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
575	ix += 1
576	values = wx.StaticText(self, -1, 'PD[ratio]')
577	polytext = "Polydispersity (= STD/mean); "
578	polytext += "the standard deviation over the mean value."
579	values.SetToolTipString(polytext)
580
581	self.sizer4_4.Add(values,( iy, ix),(1,1), wx.EXPAND\|wx.ADJUST_MINSIZE,
582	0)
583	ix +=2
584	if self.is_mac:
585	err_text = 'Error'
586	else:
587	err_text = ''
588	self.text_disp_1 = wx.StaticText(self, -1, err_text)
589	self.sizer4_4.Add( self.text_disp_1,(iy, ix),(1,1),\
590	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
591
592
593	ix +=1
594	self.text_disp_min = wx.StaticText(self, -1, 'Min')
595	self.sizer4_4.Add(self.text_disp_min,(iy, ix),(1,1),\
596	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
597
598	ix +=1
599	self.text_disp_max = wx.StaticText(self, -1, 'Max')
600	self.sizer4_4.Add(self.text_disp_max,(iy, ix),(1,1),\
601	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
602
603	ix += 1
604	npts = wx.StaticText(self, -1, 'Npts')
605	npts.SetToolTipString("Number of sampling points for the numerical\n\
606	integration over the distribution function.")
607	self.sizer4_4.Add(npts,( iy, ix),(1,1), wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
608	ix += 1
609	nsigmas = wx.StaticText(self, -1, 'Nsigs')
610	nsigmas.SetToolTipString(" Number of sigmas between which the range\n\
611	of the distribution function will be used for weighting. \n\
612	The value '3' covers 99.5% for Gaussian distribution \n\
613	function. Note: Not recommended to change this value.")
614	self.sizer4_4.Add(nsigmas,( iy, ix),(1,1), wx.EXPAND\|wx.ADJUST_MINSIZE,
615	0)
616	ix +=1
617	self.sizer4_4.Add(model_disp,(iy, ix),(1,1),\
618	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
619
620	if self.engine_type=="park":
621	self.text_disp_max.Show(True)
622	self.text_disp_min.Show(True)
623
624	for item in self.model.dispersion.keys():
625	if not item in self.model.orientation_params:
626	if not self.disp_cb_dict.has_key(item):
627	self.disp_cb_dict[item]= None
628	name0="Distribution of " + item
629	name1=item+".width"
630	name2=item+".npts"
631	name3=item+".nsigmas"
632	if not self.model.details.has_key(name1):
633	self.model.details [name1] = ["",None,None]
634
635	iy += 1
636	for p in self.model.dispersion[item].keys():
637
638	if p=="width":
639	ix = 0
640	cb = wx.CheckBox(self, -1, name0, (10, 10))
641	cb.SetToolTipString("Check mark to fit")
642	wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
643	self.sizer4_4.Add( cb,( iy, ix),(1,1),
644	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE,
645	5)
646	ix = 1
647	value= self.model.getParam(name1)
648	ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.3,20)
649	,style=wx.TE_PROCESS_ENTER)
650	ctl1.SetLabel('PD[ratio]')
651	poly_text = "Polydispersity (STD/mean) of %s\n" % item
652	poly_text += "STD: the standard deviation"
653	poly_text += " from the mean value."
654	ctl1.SetToolTipString(poly_text)
655	ctl1.SetValue(str (format_number(value, True)))
656	self.sizer4_4.Add(ctl1, (iy,ix),(1,1),wx.EXPAND)
657	## text to show error sign
658	ix = 2
659	text2=wx.StaticText(self, -1, '+/-')
660	self.sizer4_4.Add(text2,(iy, ix),(1,1),
661	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
662	if not self.is_mac:
663	text2.Hide()
664
665	ix = 3
666	ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.3,20),
667	style=0)
668
669	self.sizer4_4.Add(ctl2, (iy,ix),(1,1),
670	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
671	if not self.is_mac:
672	ctl2.Hide()
673
674	ix = 4
675	ctl3 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/2,
676	20), style=wx.TE_PROCESS_ENTER,
677	text_enter_callback = self._onparamRangeEnter)
678
679	self.sizer4_4.Add(ctl3, (iy,ix),(1,1),
680	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
681
682
683	ix = 5
684	ctl4 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/2,
685	20), style=wx.TE_PROCESS_ENTER,
686	text_enter_callback = self._onparamRangeEnter)
687
688	self.sizer4_4.Add(ctl4, (iy,ix),(1,1),
689	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
690
691
692	if self.engine_type=="park":
693	ctl3.Show(True)
694	ctl4.Show(True)
695
696	elif p=="npts":
697	ix = 6
698	value= self.model.getParam(name2)
699	Tctl = self.ModelTextCtrl(self, -1,
700	size=(_BOX_WIDTH/2.2,20),
701	style=wx.TE_PROCESS_ENTER)
702
703	Tctl.SetValue(str (format_number(value)))
704	self.sizer4_4.Add(Tctl, (iy,ix),(1,1),
705	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
706	self.fixed_param.append([None,name2, Tctl,None,None,
707	None, None,None])
708	elif p=="nsigmas":
709	ix = 7
710	value= self.model.getParam(name3)
711	Tct2 = self.ModelTextCtrl(self, -1,
712	size=(_BOX_WIDTH/2.2,20),
713	style=wx.TE_PROCESS_ENTER)
714
715	Tct2.SetValue(str (format_number(value)))
716	self.sizer4_4.Add(Tct2, (iy,ix),(1,1),
717	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
718	#ix +=1
719	#self.sizer4_4.Add((20,20), (iy,ix),(1,1),
720	# wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
721
722	self.fixed_param.append([None, name3, Tct2,
723	None, None, None,
724	None, None])
725
726
727	ix = 8
728	disp_box = wx.ComboBox(self, -1,size=(65,-1),
729	style=wx.CB_READONLY, name = '%s'% name1)
730	for key, value in self.polydisp.iteritems():
731	name_disp = str(key)
732	disp_box.Append(name_disp,value)
733	disp_box.SetStringSelection("gaussian")
734	wx.EVT_COMBOBOX(disp_box,-1, self._on_disp_func)
735	self.sizer4_4.Add(disp_box,(iy,ix),(1,1), wx.EXPAND)
736	self.fittable_param.append([cb,name1,ctl1,text2,
737	ctl2, ctl3, ctl4, disp_box])
738
739	ix =0
740	iy +=1
741	self.sizer4_4.Add((20,20),(iy,ix),(1,1),
742	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
743	first_orient = True
744	for item in self.model.dispersion.keys():
745	if item in self.model.orientation_params:
746	if not self.disp_cb_dict.has_key(item):
747	self.disp_cb_dict[item]= None
748	name0="Distribution of " + item
749	name1=item+".width"
750	name2=item+".npts"
751	name3=item+".nsigmas"
752
753	if not self.model.details.has_key(name1):
754	self.model.details [name1] = ["",None,None]
755
756
757	iy += 1
758	for p in self.model.dispersion[item].keys():
759
760	if p=="width":
761	ix = 0
762	cb = wx.CheckBox(self, -1, name0, (10, 10))
763	cb.SetToolTipString("Check mark to fit")
764	wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
765	self.sizer4_4.Add( cb,( iy, ix),(1,1),
766	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE,
767	5)
768	if self.data.__class__.__name__ == "Data2D" or \
769	self.enable2D:
770	cb.Show(True)
771	elif cb.IsShown():
772	cb.Hide()
773	ix = 1
774	value= self.model.getParam(name1)
775	ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.3,
776	20),style=wx.TE_PROCESS_ENTER)
777	poly_tip = "Absolute Sigma for %s." % item
778	ctl1.SetToolTipString(poly_tip)
779	ctl1.SetValue(str (format_number(value, True)))
780	if self.data.__class__.__name__ == "Data2D" or \
781	self.enable2D:
782	if first_orient:
783	values.SetLabel('PD[ratio], Sig[deg]')
784	poly_text = "PD(polydispersity for lengths):\n"
785	poly_text += "It should be a value between"
786	poly_text += "0 and 1\n"
787	poly_text += "Sigma for angles: \n"
788	poly_text += "It is the STD (ratio*mean)"
789	poly_text += " of the distribution.\n "
790
791	values.SetToolTipString(poly_text)
792	first_orient = False
793	ctl1.Show(True)
794	elif ctl1.IsShown():
795	ctl1.Hide()
796
797	self.sizer4_4.Add(ctl1, (iy,ix),(1,1),wx.EXPAND)
798	## text to show error sign
799	ix = 2
800	text2=wx.StaticText(self, -1, '+/-')
801	self.sizer4_4.Add(text2,(iy, ix),(1,1),
802	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
803
804	text2.Hide()
805
806	ix = 3
807	ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.3,20),
808	style=0)
809
810	self.sizer4_4.Add(ctl2, (iy,ix),(1,1),
811	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
812
813	ctl2.Hide()
814	if self.data.__class__.__name__ == "Data2D" or \
815	self.enable2D:
816	if self.is_mac:
817	text2.Show(True)
818	ctl2.Show(True)
819
820	ix = 4
821	ctl3 = self.ModelTextCtrl(self, -1,
822	size=(_BOX_WIDTH/2,20),
823	style=wx.TE_PROCESS_ENTER,
824	text_enter_callback = self._onparamRangeEnter)
825
826	self.sizer4_4.Add(ctl3, (iy,ix),(1,1),
827	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
828
829	ctl3.Hide()
830
831	ix = 5
832	ctl4 = self.ModelTextCtrl(self, -1,
833	size=(_BOX_WIDTH/2,20), style=wx.TE_PROCESS_ENTER,
834	text_enter_callback = self._onparamRangeEnter)
835	self.sizer4_4.Add(ctl4, (iy,ix),(1,1),
836	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
837	ctl4.Hide()
838
839	if self.data.__class__.__name__ == "Data2D" or \
840	self.enable2D:
841	ctl3.Show(True)
842	ctl4.Show(True)
843
844	elif p=="npts":
845	ix = 6
846	value= self.model.getParam(name2)
847	Tctl = self.ModelTextCtrl(self, -1,
848	size=(_BOX_WIDTH/2.2, 20),
849	style=wx.TE_PROCESS_ENTER)
850
851	Tctl.SetValue(str (format_number(value)))
852	if self.data.__class__.__name__ == "Data2D" or \
853	self.enable2D:
854	Tctl.Show(True)
855	else:
856	Tctl.Hide()
857	self.sizer4_4.Add(Tctl, (iy,ix),(1,1),
858	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
859	self.fixed_param.append([None,name2, Tctl,None,None,
860	None, None,None])
861	self.orientation_params_disp.append([None,name2,
862	Tctl,None,None,
863	None, None,None])
864	elif p=="nsigmas":
865	ix = 7
866	value= self.model.getParam(name3)
867	Tct2 = self.ModelTextCtrl(self, -1,
868	size=(_BOX_WIDTH/2.2,20),
869	style=wx.TE_PROCESS_ENTER)
870
871	Tct2.SetValue(str (format_number(value)))
872	if self.data.__class__.__name__ == "Data2D" or \
873	self.enable2D:
874	Tct2.Show(True)
875	else:
876	Tct2.Hide()
877	self.sizer4_4.Add(Tct2, (iy,ix),(1,1),
878	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
879
880
881	self.fixed_param.append([None,name3, Tct2
882	,None,None, None, None,None])
883
884	self.orientation_params_disp.append([None,name3,
885	Tct2 ,None,None, None, None,None])
886
887
888	ix = 8
889	disp_box = wx.ComboBox(self, -1,size=(65,-1),
890	style=wx.CB_READONLY, name = '%s'% name1)
891	for key, value in self.polydisp.iteritems():
892	name_disp = str(key)
893	disp_box.Append(name_disp,value)
894	disp_box.SetStringSelection("gaussian")
895	wx.EVT_COMBOBOX(disp_box,-1, self._on_disp_func)
896	self.sizer4_4.Add(disp_box,(iy,ix),(1,1), wx.EXPAND)
897	self.fittable_param.append([cb,name1,ctl1,text2,
898	ctl2, ctl3, ctl4, disp_box])
899	self.orientation_params_disp.append([cb,name1,ctl1,
900	text2, ctl2, ctl3, ctl4, disp_box])
901
902	if self.data.__class__.__name__ == "Data2D" or \
903	self.enable2D:
904	disp_box.Show(True)
905	else:
906	disp_box.Hide()
907
908
909	self.state.disp_cb_dict = copy.deepcopy(self.disp_cb_dict)
910
911	self.state.model = self.model.clone()
912	## save state into
913	self.state.cb1 = self.cb1.GetValue()
914	self._copy_parameters_state(self.parameters, self.state.parameters)
915	self._copy_parameters_state(self.orientation_params_disp,
916	self.state.orientation_params_disp)
917	self._copy_parameters_state(self.fittable_param,
918	self.state.fittable_param)
919	self._copy_parameters_state(self.fixed_param, self.state.fixed_param)
920
921
922	wx.PostEvent(self.parent, StatusEvent(status=\
923	" Selected Distribution: Gaussian"))
924	#Fill the list of fittable parameters
925	self.select_all_param(event=None)
926
927	self.Layout()
928
929
930	def _onDraw(self, event):
931	"""
932	Update and Draw the model
933	"""
934	if self.model ==None:
935	msg="Please select a Model first..."
936	wx.MessageBox(msg, 'Info')
937
938	return
939	flag = self._update_paramv_on_fit()
940	self._onparamEnter_helper()
941	if not flag:
942	msg= "The parameters are invalid"
943	wx.PostEvent(self.parent.parent, StatusEvent(status= msg ))
944	return
945
946
947	def _onFit(self, event):
948	"""
949	Allow to fit
950	"""
951	if event != None:
952	event.Skip()
953	if len(self.parent._manager.fit_thread_list)>0 and\
954	self.parent._manager._fit_engine != "park":
955	msg = "The FitEnging will be set to 'ParkMC'\n"
956	msg += " to fit with more than one data set..."
957	wx.MessageBox(msg, 'Info')
958	#wx.PostEvent(self._manager.parent, StatusEvent(status=\
959	# "Fitting: %s"%msg))
960
961	if self.data is None:
962	msg = "Please get Data first..."
963	wx.MessageBox(msg, 'Info')
964	wx.PostEvent(self._manager.parent, StatusEvent(status=\
965	"Fit: %s" % msg))
966	return
967	if self.model is None:
968	msg = "Please select a Model first..."
969	wx.MessageBox(msg, 'Info')
970	wx.PostEvent(self._manager.parent, StatusEvent(status=\
971	"Fit: %s"%msg, type="stop"))
972	return
973
974	if len(self.param_toFit) <= 0:
975	msg= "Select at least one parameter to fit"
976	wx.MessageBox(msg, 'Info')
977	wx.PostEvent(self._manager.parent, StatusEvent(status= msg,
978	type="stop" ))
979	return
980
981	flag = self._update_paramv_on_fit()
982
983	if not flag:
984	msg= "Fitting range or parameters are invalid"
985	wx.PostEvent(self.parent.parent, StatusEvent(status= msg,
986	type="stop"))
987	return
988
989	self.select_param(event =None)
990
991	#Clear errors if exist from previous fitting
992	#self._clear_Err_on_Fit()
993
994	# Remove or do not allow fitting on the Q=0 point, especially
995	# when y(q=0)=None at x[0].
996	self.qmin_x = float(self.qmin.GetValue())
997	self.qmax_x = float( self.qmax.GetValue())
998	self._manager._reset_schedule_problem(value=0, uid=self.uid)
999	self._manager.schedule_for_fit(uid=self.uid,value=1, fitproblem =None)
1000	self._manager.set_fit_range(uid=self.uid,qmin= self.qmin_x,
1001	qmax= self.qmax_x)
1002
1003	#single fit
1004	self._manager.onFit(uid=self.uid)
1005	self.btFit.SetLabel("Stop")
1006	self.bind_fit_button()
1007
1008
1009	def bind_fit_button(self):
1010	"""
1011	bind the fit button to either fit handler or stop fit handler
1012	"""
1013	self.btFit.Unbind(event=wx.EVT_BUTTON, id= self.btFit.GetId())
1014	if self.btFit.GetLabel().lower() == "stop":
1015	self.fit_started = True
1016	self.btFit.SetForegroundColour('red')
1017	self.btFit.Bind(event=wx.EVT_BUTTON, handler=self._StopFit,
1018	id=self.btFit.GetId())
1019	elif self.btFit.GetLabel().lower() == "fit":
1020	self.fit_started = False
1021	self.btFit.SetDefault()
1022	self.btFit.SetForegroundColour('black')
1023	#self.btFit.SetBackgroundColour(self.default_bt_colour)
1024	self.btFit.Bind(event=wx.EVT_BUTTON, handler=self._onFit,
1025	id=self.btFit.GetId())
1026	else:
1027	msg = "FitPage: fit button has unknown label"
1028	raise ValuerError, msg
1029	self._manager._reset_schedule_problem(value=0)
1030
1031	def is_fitting(self):
1032	if self.fit_started:
1033	self._StopFit(event=None)
1034
1035	def _StopFit(self, event=None):
1036	"""
1037	Stop fit
1038	"""
1039	#time.sleep(0.1)
1040	if event != None:
1041	event.Skip()
1042	#if self.engine_type=="scipy":
1043	self._manager.stop_fit(self.uid)
1044	self._manager._reset_schedule_problem(value=0)
1045	self._on_fit_complete()
1046
1047	def _on_select_model(self, event=None):
1048	"""
1049	call back for model selection
1050	"""
1051
1052	self.Show(False)
1053	self._on_select_model_helper()
1054	self.set_model_param_sizer(self.model)
1055	if self.model is None:
1056	self._set_bookmark_flag(False)
1057	self._keep.Enable(False)
1058	self._set_save_flag(False)
1059	self.enable_disp.SetValue(False)
1060	self.disable_disp.SetValue(True)
1061	try:
1062	self.set_dispers_sizer()
1063	except:
1064	pass
1065	#self.btFit.SetFocus()
1066	self.state.enable_disp = self.enable_disp.GetValue()
1067	self.state.disable_disp = self.disable_disp.GetValue()
1068	self.state.pinhole_smearer = self.pinhole_smearer.GetValue()
1069	self.state.slit_smearer = self.slit_smearer.GetValue()
1070
1071	self.state.structurecombobox = self.structurebox.GetCurrentSelection()
1072	self.state.formfactorcombobox = self.formfactorbox.GetCurrentSelection()
1073
1074	if self.model != None:
1075	self._set_copy_flag(True)
1076	self._set_paste_flag(True)
1077	if self.data != None:
1078	self._set_bookmark_flag(True)
1079	self._keep.Enable(True)
1080	#self._set_save_flag(True)
1081	# Reset smearer, model and data
1082	self._set_smear(self.data)
1083	try:
1084	# update smearer sizer
1085	self.onSmear(None)
1086	temp_smear = None
1087	if self.enable_smearer.GetValue():
1088	# Set the smearer environments
1089	temp_smear = self.smearer
1090	except:
1091	raise
1092	## error occured on chisqr computation
1093	#pass
1094	## event to post model to fit to fitting plugins
1095	(ModelEventbox, EVT_MODEL_BOX) = wx.lib.newevent.NewEvent()
1096
1097	## set smearing value whether or not
1098	# the data contain the smearing info
1099	evt = ModelEventbox(model=self.model,
1100	smearer=temp_smear,
1101	qmin=float(self.qmin_x),
1102	uid=self.uid,
1103	qmax=float(self.qmax_x))
1104
1105	self._manager._on_model_panel(evt=evt)
1106	self.mbox_description.SetLabel("Model [%s]" % str(self.model.name))
1107	self.state.model = self.model.clone()
1108	self.state.model.name = self.model.name
1109
1110
1111	if event != None:
1112	## post state to fit panel
1113	new_event = PageInfoEvent(page = self)
1114	wx.PostEvent(self.parent, new_event)
1115	#update list of plugins if new plugin is available
1116	if self.plugin_rbutton.GetValue():
1117	temp = self.parent.update_model_list()
1118	if temp:
1119	self.model_list_box = temp
1120	current_val = self.formfactorbox.GetValue()
1121	pos = self.formfactorbox.GetSelection()
1122	self._show_combox_helper()
1123	self.formfactorbox.SetSelection(pos)
1124	self.formfactorbox.SetValue(current_val)
1125	self._onDraw(event=None)
1126	else:
1127	self._draw_model()
1128	self.SetupScrolling()
1129	self.Show(True)
1130
1131	def _onparamEnter(self,event):
1132	"""
1133	when enter value on panel redraw model according to changed
1134	"""
1135	if self.model ==None:
1136	msg="Please select a Model first..."
1137	wx.MessageBox(msg, 'Info')
1138	#wx.PostEvent(self._manager.parent, StatusEvent(status=\
1139	# "Parameters: %s"%msg))
1140	return
1141
1142	#default flag
1143	flag = False
1144	self.fitrange = True
1145	#get event object
1146	tcrtl = event.GetEventObject()
1147
1148	#wx.PostEvent(self._manager.parent, StatusEvent(status=" \
1149	# updating ... ",type="update"))
1150	#Clear msg if previously shown.
1151	msg= ""
1152	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1153
1154	if check_float(tcrtl):
1155	flag = self._onparamEnter_helper()
1156	self.set_npts2fit()
1157	if self.fitrange:
1158	temp_smearer = None
1159	if not self.disable_smearer.GetValue():
1160	temp_smearer= self.current_smearer
1161	## set smearing value whether or not
1162	# the data contain the smearing info
1163	if self.slit_smearer.GetValue():
1164	flag1 = self.update_slit_smear()
1165	flag = flag or flag1
1166	elif self.pinhole_smearer.GetValue():
1167	flag1 = self.update_pinhole_smear()
1168	flag = flag or flag1
1169	elif self.data.__class__.__name__ != "Data2D" and \
1170	not self.enable2D:
1171	self._manager.set_smearer(smearer=temp_smearer,
1172	enable2D=self.enable2D,
1173	uid=self.uid,
1174	qmin= float(self.qmin_x),
1175	qmax= float(self.qmax_x),
1176	draw=True)
1177	if flag:
1178	#self.compute_chisqr(smearer= temp_smearer)
1179
1180	## new state posted
1181	if self.state_change:
1182	#self._undo.Enable(True)
1183	event = PageInfoEvent(page = self)
1184	wx.PostEvent(self.parent, event)
1185	self.state_change= False
1186	else:
1187	# invalid fit range: do nothing here:
1188	# msg already displayed in validate
1189	return
1190	else:
1191	self.save_current_state()
1192	msg= "Cannot Plot :Must enter a number!!! "
1193	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1194
1195	self.save_current_state()
1196	return
1197
1198	def _onparamRangeEnter(self, event):
1199	"""
1200	Check validity of value enter in the parameters range field
1201	"""
1202	#if self.check_invalid_panel():
1203	# return
1204	tcrtl= event.GetEventObject()
1205	#Clear msg if previously shown.
1206	msg= ""
1207	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1208	# Flag to register when a parameter has changed.
1209	is_modified = False
1210	if tcrtl.GetValue().lstrip().rstrip()!="":
1211	try:
1212	value = float(tcrtl.GetValue())
1213	tcrtl.SetBackgroundColour(wx.WHITE)
1214	self._check_value_enter(self.fittable_param ,is_modified)
1215	self._check_value_enter(self.parameters ,is_modified)
1216	except:
1217	tcrtl.SetBackgroundColour("pink")
1218	msg= "Model Error:wrong value entered : %s"% sys.exc_value
1219	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1220	return
1221	else:
1222	tcrtl.SetBackgroundColour(wx.WHITE)
1223
1224	#self._undo.Enable(True)
1225	self.save_current_state()
1226	event = PageInfoEvent(page = self)
1227	wx.PostEvent(self.parent, event)
1228	self.state_change= False
1229
1230
1231	def _onQrangeEnter(self, event):
1232	"""
1233	Check validity of value enter in the Q range field
1234	"""
1235	#if self.check_invalid_panel():
1236	# return
1237	tcrtl= event.GetEventObject()
1238	#Clear msg if previously shown.
1239	msg= ""
1240	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1241	# Flag to register when a parameter has changed.
1242	is_modified = False
1243	if tcrtl.GetValue().lstrip().rstrip()!="":
1244	try:
1245	value = float(tcrtl.GetValue())
1246	tcrtl.SetBackgroundColour(wx.WHITE)
1247
1248	# If qmin and qmax have been modified, update qmin and qmax
1249	if self._validate_qrange( self.qmin, self.qmax):
1250	tempmin = float(self.qmin.GetValue())
1251	if tempmin != self.qmin_x:
1252	self.qmin_x = tempmin
1253	tempmax = float(self.qmax.GetValue())
1254	if tempmax != self.qmax_x:
1255	self.qmax_x = tempmax
1256	else:
1257	tcrtl.SetBackgroundColour("pink")
1258	msg= "Model Error:wrong value entered : %s"% sys.exc_value
1259	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1260	return
1261
1262	except:
1263	tcrtl.SetBackgroundColour("pink")
1264	msg= "Model Error:wrong value entered : %s"% sys.exc_value
1265	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1266	return
1267	#Check if # of points for theory model are valid(>0).
1268	# check for 2d
1269	if self.data is not None:
1270	if self.data.__class__.__name__ == "Data2D" or \
1271	self.enable2D:
1272	# set mask
1273	radius= numpy.sqrt( self.data.qx_data*self.data.qx_data +
1274	self.data.qy_data*self.data.qy_data )
1275	index_data = ((self.qmin_x <= radius)& \
1276	(radius<= self.qmax_x))
1277	index_data = (index_data)&(self.data.mask)
1278	index_data = (index_data)&(numpy.isfinite(self.data.data))
1279	if len(index_data[index_data]) < 10:
1280	msg = "Cannot Plot :No or too little npts in"
1281	msg += " that data range!!! "
1282	wx.PostEvent(self.parent.parent,
1283	StatusEvent(status=msg))
1284	return
1285	else:
1286	self.data.mask = index_data
1287	#self.Npts_fit.SetValue(str(len(self.data.mask)))
1288	self.set_npts2fit()
1289	else:
1290	index_data = ((self.qmin_x <= self.data.x)& \
1291	(self.data.x <= self.qmax_x))
1292	self.Npts_fit.SetValue(str(len(self.data.x[index_data])))
1293	else:
1294	self.npts_x = self.Npts_total.GetValue()
1295	self._save_plotting_range()
1296
1297	else:
1298	tcrtl.SetBackgroundColour("pink")
1299	msg= "Model Error:wrong value entered!!!"
1300	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1301
1302	self._draw_model()
1303	##update chi2
1304	#self.compute_chisqr(smearer= temp_smearer)
1305	#self._undo.Enable(True)
1306	self.save_current_state()
1307	event = PageInfoEvent(page = self)
1308	wx.PostEvent(self.parent, event)
1309	self.state_change= False
1310	return
1311
1312	def _clear_Err_on_Fit(self):
1313	"""
1314	hide the error text control shown
1315	after fitting
1316	"""
1317	if self.is_mac:
1318	return
1319	if hasattr(self,"text2_3"):
1320	self.text2_3.Hide()
1321
1322	if len(self.parameters)>0:
1323	for item in self.parameters:
1324	#Skip t ifhe angle parameters if 1D data
1325	if self.data.__class__.__name__ != "Data2D" and \
1326	not self.enable2D:
1327	if item in self.orientation_params:
1328	continue
1329	if item in self.param_toFit:
1330	continue
1331	## hide statictext +/-
1332	if len(item) < 4 :
1333	continue
1334	if item[3]!=None and item[3].IsShown():
1335	item[3].Hide()
1336	## hide textcrtl for error after fit
1337	if item[4]!=None and item[4].IsShown():
1338	item[4].Hide()
1339
1340	if len(self.fittable_param)>0:
1341	for item in self.fittable_param:
1342	#Skip t ifhe angle parameters if 1D data
1343	if self.data.__class__.__name__ != "Data2D" and \
1344	not self.enable2D:
1345	if item in self.orientation_params:
1346	continue
1347	if item in self.param_toFit:
1348	continue
1349	if len(item) < 4 :
1350	continue
1351	## hide statictext +/-
1352	if item[3]!=None and item[3].IsShown():
1353	item[3].Hide()
1354	## hide textcrtl for error after fit
1355	if item[4]!=None and item[4].IsShown():
1356	item[4].Hide()
1357	return
1358
1359	def _get_defult_custom_smear(self):
1360	"""
1361	Get the defult values for custum smearing.
1362	"""
1363	# get the default values
1364	if self.dxl == None: self.dxl = 0.0
1365	if self.dxw == None: self.dxw = ""
1366	if self.dx_min == None: self.dx_min = SMEAR_SIZE_L
1367	if self.dx_max == None: self.dx_max = SMEAR_SIZE_H
1368
1369	def _get_smear_info(self):
1370	"""
1371	Get the smear info from data.
1372
1373	:return: self.smear_type, self.dq_l and self.dq_r,
1374	respectively the type of the smear, dq_min and
1375	dq_max for pinhole smear data
1376	while dxl and dxw for slit smear
1377
1378	"""
1379
1380	# default
1381	self.smear_type = None
1382	self.dq_l = None
1383	self.dq_r = None
1384	data = self.data
1385	if self.data is None:
1386	return
1387	elif self.data.__class__.__name__ == "Data2D" or \
1388	self.enable2D:
1389	if data.dqx_data == None or data.dqy_data ==None:
1390	return
1391	elif self.smearer != None and (data.dqx_data.any()!=0) and \
1392	(data.dqx_data.any()!=0):
1393	self.smear_type = "Pinhole2d"
1394	self.dq_l = format_number(numpy.average(data.dqx_data))
1395	self.dq_r = format_number(numpy.average(data.dqy_data))
1396	return
1397	else:
1398	return
1399	# check if it is pinhole smear and get min max if it is.
1400	if data.dx != None and all(data.dx !=0):
1401	self.smear_type = "Pinhole"
1402	self.dq_l = data.dx[0]
1403	self.dq_r = data.dx[-1]
1404
1405	# check if it is slit smear and get min max if it is.
1406	elif data.dxl != None or data.dxw != None:
1407	self.smear_type = "Slit"
1408	if data.dxl != None and all(data.dxl !=0):
1409	self.dq_l = data.dxl[0]
1410	if data.dxw != None and all(data.dxw !=0):
1411	self.dq_r = data.dxw[0]
1412	#return self.smear_type,self.dq_l,self.dq_r
1413
1414	def _show_smear_sizer(self):
1415	"""
1416	Show only the sizers depending on smear selection
1417	"""
1418	# smear disabled
1419	if self.disable_smearer.GetValue():
1420	self.smear_description_none.Show(True)
1421	# 2Dsmear
1422	elif self._is_2D():
1423	self.smear_description_accuracy_type.Show(True)
1424	self.smear_accuracy.Show(True)
1425	self.smear_description_accuracy_type.Show(True)
1426	self.smear_description_2d.Show(True)
1427	self.smear_description_2d_x.Show(True)
1428	self.smear_description_2d_y.Show(True)
1429	if self.pinhole_smearer.GetValue():
1430	self.smear_pinhole_min.Show(True)
1431	self.smear_pinhole_max.Show(True)
1432	# smear from data
1433	elif self.enable_smearer.GetValue():
1434
1435	self.smear_description_dqdata.Show(True)
1436	if self.smear_type != None:
1437	self.smear_description_smear_type.Show(True)
1438	if self.smear_type == 'Slit':
1439	self.smear_description_slit_height.Show(True)
1440	self.smear_description_slit_width.Show(True)
1441	elif self.smear_type == 'Pinhole':
1442	self.smear_description_pin_min.Show(True)
1443	self.smear_description_pin_max.Show(True)
1444	self.smear_description_smear_type.Show(True)
1445	self.smear_description_type.Show(True)
1446	self.smear_data_left.Show(True)
1447	self.smear_data_right.Show(True)
1448	# custom pinhole smear
1449	elif self.pinhole_smearer.GetValue():
1450	if self.smear_type == 'Pinhole':
1451	self.smear_message_new_p.Show(True)
1452	self.smear_description_pin_min.Show(True)
1453	self.smear_description_pin_max.Show(True)
1454
1455	self.smear_pinhole_min.Show(True)
1456	self.smear_pinhole_max.Show(True)
1457	# custom slit smear
1458	elif self.slit_smearer.GetValue():
1459	self.smear_message_new_s.Show(True)
1460	self.smear_description_slit_height.Show(True)
1461	self.smear_slit_height.Show(True)
1462	self.smear_description_slit_width.Show(True)
1463	self.smear_slit_width.Show(True)
1464
1465	def _hide_all_smear_info(self):
1466	"""
1467	Hide all smearing messages in the set_smearer sizer
1468	"""
1469	self.smear_description_none.Hide()
1470	self.smear_description_dqdata.Hide()
1471	self.smear_description_type.Hide()
1472	self.smear_description_smear_type.Hide()
1473	self.smear_description_accuracy_type.Hide()
1474	self.smear_description_2d_x.Hide()
1475	self.smear_description_2d_y.Hide()
1476	self.smear_description_2d.Hide()
1477
1478	self.smear_accuracy.Hide()
1479	self.smear_data_left.Hide()
1480	self.smear_data_right.Hide()
1481	self.smear_description_pin_min.Hide()
1482	self.smear_pinhole_min.Hide()
1483	self.smear_description_pin_max.Hide()
1484	self.smear_pinhole_max.Hide()
1485	self.smear_description_slit_height.Hide()
1486	self.smear_slit_height.Hide()
1487	self.smear_description_slit_width.Hide()
1488	self.smear_slit_width.Hide()
1489	self.smear_message_new_p.Hide()
1490	self.smear_message_new_s.Hide()
1491
1492	def _set_accuracy_list(self):
1493	"""
1494	Set the list of an accuracy in 2D custum smear: Xhigh, High, Med, or Low
1495	"""
1496	# list of accuracy choices
1497	list = ['Low','Med','High','Xhigh']
1498	for idx in range(len(list)):
1499	self.smear_accuracy.Append(list[idx],idx)
1500
1501	def _set_fun_box_list(self,fun_box):
1502	"""
1503	Set the list of func for multifunctional models
1504	"""
1505	# Check if it is multi_functional model
1506	if self.model.__class__ not in self.model_list_box["Multi-Functions"] \
1507	and not self.temp_multi_functional:
1508	return None
1509	# Get the func name list
1510	list = self.model.fun_list
1511	if len(list) == 0:
1512	return None
1513	# build function (combo)box
1514	ind = 0
1515	while(ind < len(list)):
1516	for key, val in list.iteritems():
1517	if (val == ind):
1518	fun_box.Append(key,val)
1519	break
1520	ind += 1
1521
1522	def _on_select_accuracy(self,event):
1523	"""
1524	Select an accuracy in 2D custom smear: Xhigh, High, Med, or Low
1525	"""
1526	#event.Skip()
1527	# Check if the accuracy is same as before
1528	#self.smear2d_accuracy = event.GetEventObject().GetValue()
1529	self.smear2d_accuracy = self.smear_accuracy.GetValue()
1530	if self.pinhole_smearer.GetValue():
1531	self.onPinholeSmear(event=None)
1532	else:
1533	self.onSmear(event=None)
1534	if self.current_smearer != None:
1535	self.current_smearer.set_accuracy(accuracy =
1536	self.smear2d_accuracy)
1537	event.Skip()
1538
1539	def _on_fun_box(self,event):
1540	"""
1541	Select an func: Erf,Rparabola,LParabola
1542	"""
1543	fun_val = None
1544	fun_box = event.GetEventObject()
1545	name = fun_box.Name
1546	value = fun_box.GetValue()
1547	if self.model.fun_list.has_key(value):
1548	fun_val = self.model.fun_list[value]
1549
1550	self.model.setParam(name,fun_val)
1551	# save state
1552	self._copy_parameters_state(self.str_parameters,
1553	self.state.str_parameters)
1554	# update params
1555	self._update_paramv_on_fit()
1556	# draw
1557	self._draw_model()
1558	self.Refresh()
1559	# get ready for new event
1560	event.Skip()
1561
1562	def _onMask(self, event):
1563	"""
1564	Build a panel to allow to edit Mask
1565	"""
1566
1567	from sans.guiframe.local_perspectives.plotting.masking \
1568	import MaskPanel as MaskDialog
1569
1570	self.panel = MaskDialog(base=self, data=self.data, id=wx.NewId())
1571	#self.panel.Bind(wx.EVT_CLOSE, self._draw_masked_model)
1572	self.panel.ShowModal()
1573	#wx.PostEvent(self.parent, event)
1574
1575	def _draw_masked_model(self, event):
1576	"""
1577	Draw model image w/mask
1578	"""
1579	#event.Skip()
1580	is_valid_qrange = self._update_paramv_on_fit()
1581
1582	if is_valid_qrange and self.model != None:
1583	#self.panel.Show(0)
1584	#self.panel.Destroy() # frame
1585	self.panel.MakeModal(False)
1586	event.Skip()
1587	# try re draw the model plot if it exists
1588	self._draw_model()
1589	self.set_npts2fit()
1590	elif self.model == None:
1591	#self.panel.Show(0)
1592	#self.panel.Destroy() # frame
1593	self.panel.MakeModal(False)
1594	event.Skip()
1595	self.set_npts2fit()
1596	msg= "No model is found on updating MASK in the model plot... "
1597	wx.PostEvent(self.parent.parent, StatusEvent(status = msg ))
1598	else:
1599	event.Skip()
1600	msg = ' Please consider your Q range, too.'
1601	self.panel.ShowMessage(msg)
1602
1603
1604	def _set_smear(self, data):
1605	"""
1606	"""
1607	if data is None:
1608	return
1609	self.smearer = smear_selection(data, self.model)
1610	self.disable_smearer.SetValue(True)
1611	if self.smearer == None:
1612	self.enable_smearer.Disable()
1613	else:
1614	self.enable_smearer.Enable()
1615
1616	def _mac_sleep(self, sec=0.2):
1617	"""
1618	Give sleep to MAC
1619	"""
1620	if self.is_mac:
1621	time.sleep(sec)
1622
1623	def get_view_mode(self):
1624	"""
1625	return True if the panel allow 2D or False if 1D
1626	"""
1627	return self.enable2D
1628
1629	def compute_data_range(self, data):
1630	"""
1631	compute the minimum and the maximum range of the data
1632	return the npts contains in data
1633	:param data:
1634	"""
1635	qmin, qmax, npts = None, None, None
1636	if data is not None:
1637	if not hasattr(data,"data"):
1638	qmin = min(data.x)
1639	# Maximum value of data
1640	qmax = max(data.x)
1641	npts = len(data.x)
1642	else:
1643	qmin = 0
1644	x = max(math.fabs(data.xmin), math.fabs(data.xmax))
1645	y = max(math.fabs(data.ymin), math.fabs(data.ymax))
1646	## Maximum value of data
1647	qmax = math.sqrt(xx + yy)
1648	npts = len(data.data)
1649	return qmin, qmax, npts
1650
1651
1652	def set_data(self, data):
1653	"""
1654	reset the current data
1655	"""
1656	id = None
1657	group_id = None
1658	flag = False
1659	if self.data is None and data is not None:
1660	flag = True
1661	if data is not None:
1662	id = data.id
1663	group_id = data.group_id
1664	if self.data is not None:
1665	flag = (data.id != self.data.id)
1666	self.data = data
1667	if self.data is None:
1668	data_min = ""
1669	data_max = ""
1670	data_name = ""
1671	self._set_bookmark_flag(False)
1672	self._keep.Enable(False)
1673	self._set_save_flag(False)
1674	else:
1675	if self.model != None:
1676	self._set_bookmark_flag(True)
1677	self._keep.Enable(True)
1678	self._set_save_flag(True)
1679	self._set_preview_flag(True)
1680
1681	self._set_smear(data)
1682	# more disables for 2D
1683	if self.data.__class__.__name__ == "Data2D" or \
1684	self.enable2D:
1685	self.slit_smearer.Disable()
1686	self.pinhole_smearer.Enable(True)
1687	self.default_mask = copy.deepcopy(self.data.mask)
1688	else:
1689	self.slit_smearer.Enable(True)
1690	self.pinhole_smearer.Enable(True)
1691
1692	self.formfactorbox.Enable()
1693	self.structurebox.Enable()
1694	data_name = self.data.name
1695	data_min, data_max, npts = self.compute_data_range(self.data)
1696	#set maximum range for x in linear scale
1697	if not hasattr(self.data,"data"): #Display only for 1D data fit
1698	self.btEditMask.Disable()
1699	self.EditMask_title.Disable()
1700	else:
1701
1702	self.btEditMask.Enable()
1703	self.EditMask_title.Enable()
1704	self.Npts_total.SetValue(str(npts))
1705	#default:number of data points selected to fit
1706	self.Npts_fit.SetValue(str(npts))
1707	self.Npts_total.SetEditable(False)
1708	self.Npts_total.SetBackgroundColour(\
1709	self.GetParent().GetBackgroundColour())
1710
1711	self.Npts_total.Bind(wx.EVT_MOUSE_EVENTS, self._npts_click)
1712	#self.Npts_total.Disable()
1713	self.dataSource.SetValue(data_name)
1714	self.qmin_x = data_min
1715	self.qmax_x = data_max
1716	#self.minimum_q.SetValue(str(data_min))
1717	#self.maximum_q.SetValue(str(data_max))
1718	if data_min is None:
1719	data_min = ""
1720	if data_max is None:
1721	data_max = ""
1722	self.qmin.SetValue(str(data_min))
1723	self.qmax.SetValue(str(data_max))
1724	self.qmin.SetBackgroundColour("white")
1725	self.qmax.SetBackgroundColour("white")
1726	self.state.data = data
1727	self.state.qmin = self.qmin_x
1728	self.state.qmax = self.qmax_x
1729
1730	#update model plot with new data information
1731	if flag:
1732	#set model view button
1733	if self.data.__class__.__name__ == "Data2D":
1734	self.enable2D = True
1735	self.model_view.SetLabel("2D Mode")
1736	else:
1737	self.enable2D = False
1738	self.model_view.SetLabel("1D Mode")
1739
1740	self.model_view.Disable()
1741
1742	#replace data plot on combo box selection
1743	#by removing the previous selected data
1744	wx.PostEvent(self._manager.parent, NewPlotEvent(action="remove",
1745	group_id=group_id, id=id))
1746	#plot the current selected data
1747	wx.PostEvent(self._manager.parent, NewPlotEvent(plot=self.data,
1748	title=str(self.data.title)))
1749	self._manager.store_data(uid=self.uid, data=data,
1750	data_list=self.data_list,
1751	caption=self.window_name)
1752	self._draw_model()
1753
1754	def _npts_click(self, event):
1755	"""
1756	Prevent further handling of the mouse event on Npts_total
1757	by not calling Skip().
1758	"""
1759	pass
1760
1761	def reset_page(self, state,first=False):
1762	"""
1763	reset the state
1764	"""
1765	self.reset_page_helper(state)
1766	#import sans.guiframe.gui_manager
1767	#evt = ModelEventbox(model=state.model)
1768	#wx.PostEvent(self.event_owner, evt)
1769
1770	if self.engine_type != None:
1771	self._manager._on_change_engine(engine=self.engine_type)
1772
1773	self.select_param(event = None)
1774	#Save state_fit
1775	self.save_current_state_fit()
1776	self._lay_out()
1777	self.Refresh()
1778
1779	def get_range(self):
1780	"""
1781	return the fitting range
1782	"""
1783	return float(self.qmin_x) , float(self.qmax_x)
1784
1785	def get_npts2fit(self):
1786	"""
1787	return numbers of data points within qrange
1788
1789	:Note: This is for Park where chi2 is not normalized by Npts of fit
1790
1791	"""
1792	if self.data is None:
1793	return
1794	npts2fit = 0
1795	qmin,qmax = self.get_range()
1796	if self.data.__class__.__name__ == "Data2D" or \
1797	self.enable2D:
1798	radius= numpy.sqrt( self.data.qx_data*self.data.qx_data +
1799	self.data.qy_data*self.data.qy_data )
1800	index_data = (self.qmin_x <= radius)&(radius<= self.qmax_x)
1801	index_data= (index_data)&(self.data.mask)
1802	index_data = (index_data)&(numpy.isfinite(self.data.data))
1803	npts2fit = len(self.data.data[index_data])
1804	else:
1805	for qx in self.data.x:
1806	if qx >= qmin and qx <= qmax:
1807	npts2fit += 1
1808	return npts2fit
1809
1810	def set_npts2fit(self):
1811	"""
1812	setValue Npts for fitting
1813	"""
1814	self.Npts_fit.SetValue(str(self.get_npts2fit()))
1815
1816	def get_chi2(self):
1817	"""
1818	return the current chi2
1819	"""
1820	return self.tcChi.GetValue()
1821
1822	def get_param_list(self):
1823	"""
1824	:return self.param_toFit: list containing references to TextCtrl
1825	checked.Theses TextCtrl will allow reference to parameters to fit.
1826
1827	:raise: if return an empty list of parameter fit will nnote work
1828	properly so raise ValueError,"missing parameter to fit"
1829	"""
1830	if self.param_toFit !=[]:
1831	return self.param_toFit
1832	else:
1833	msg = "missing parameters to fit"
1834	wx.MessageBox(msg, 'warning')
1835	return False
1836	#raise ValueError,"missing parameters to fit"
1837
1838	def onsetValues(self, chisqr, p_name, out, cov):
1839	"""
1840	Build the panel from the fit result
1841
1842	:param chisqr: Value of the goodness of fit metric
1843	:param p_name: the name of parameters
1844	:param out: list of parameter with the best value found during fitting
1845	:param cov: Covariance matrix
1846
1847	"""
1848	if out == None or not numpy.isfinite(chisqr):
1849	raise ValueError,"Fit error occured..."
1850
1851	is_modified = False
1852	has_error = False
1853	dispersity = ''
1854
1855	#Hide textctrl boxes of errors.
1856	self._clear_Err_on_Fit()
1857
1858	#Check if chi2 is finite
1859	if chisqr != None and numpy.isfinite(chisqr):
1860	#format chi2
1861	if self.engine_type == "park":
1862	npt_fit = float(self.get_npts2fit())
1863	if npt_fit > 0:
1864	chisqr =chisqr/npt_fit
1865	chi2 = format_number(chisqr, True)
1866	self.tcChi.SetValue(chi2)
1867	self.tcChi.Refresh()
1868	else:
1869	self.tcChi.SetValue("-")
1870
1871	#Hide error title
1872	if self.text2_3.IsShown() and not self.is_mac:
1873	self.text2_3.Hide()
1874
1875	try:
1876	if self.enable_disp.GetValue():
1877	if hasattr(self,"text_disp_1" ):
1878	if self.text_disp_1 != None and not self.is_mac:
1879	self.text_disp_1.Hide()
1880	except:
1881	dispersity = None
1882	pass
1883	#set the panel when fit result are float not list
1884	if out.__class__== numpy.float64:
1885	self.param_toFit[0][2].SetValue(format_number(out, True))
1886
1887	if self.param_toFit[0][4].IsShown and not self.is_mac:
1888	self.param_toFit[0][4].Hide()
1889	if cov !=None :
1890	self.text2_3.Show(True)
1891	try:
1892	if self.enable_disp.GetValue():
1893	if hasattr(self,"text_disp_1" ):
1894	if self.text_disp_1 !=None:
1895	self.text_disp_1.Show(True)
1896	except:
1897	pass
1898
1899	if cov[0]==None or not numpy.isfinite(cov[0]):
1900	if self.param_toFit[0][3].IsShown and not self.is_mac:
1901	self.param_toFit[0][3].Hide()
1902	else:
1903	self.param_toFit[0][3].Show(True)
1904	self.param_toFit[0][4].Show(True)
1905	self.param_toFit[0][4].SetValue(format_number(cov[0], True))
1906	has_error = True
1907	else:
1908
1909	i = 0
1910	#Set the panel when fit result are list
1911	for item in self.param_toFit:
1912	if len(item)>5 and item != None:
1913	## reset error value to initial state
1914	if not self.is_mac:
1915	item[3].Hide()
1916	item[4].Hide()
1917
1918	for ind in range(len(out)):
1919
1920	if item[1] == p_name[ind]:
1921	break
1922	if len(out)<=len(self.param_toFit) and out[ind] !=None:
1923	val_out = format_number(out[ind], True)
1924	item[2].SetValue(val_out)
1925
1926
1927	if(cov !=None):
1928
1929	try:
1930	if dispersity !=None:
1931	if self.enable_disp.GetValue():
1932	if hasattr(self,"text_disp_1" ):
1933	if self.text_disp_1!=None:
1934	if not self.text_disp_1.IsShown()\
1935	and not self.is_mac:
1936	self.text_disp_1.Show(True)
1937	except:
1938	pass
1939
1940	if cov[ind]!=None :
1941	if numpy.isfinite(float(cov[ind])):
1942	val_err = format_number(cov[ind], True)
1943	if not self.is_mac:
1944	item[3].Show(True)
1945	item[4].Show(True)
1946	item[4].SetValue(val_err)
1947
1948	has_error = True
1949	i += 1
1950	#Show error title when any errors displayed
1951	if has_error:
1952	if not self.text2_3.IsShown():
1953	self.text2_3.Show(True)
1954	## save current state
1955	self.save_current_state()
1956
1957	#self._lay_out()
1958	if not self.is_mac:
1959	self.Layout()
1960	self.Refresh()
1961	self._mac_sleep(0.1)
1962	#plot model ( when drawing, do not update chisqr value again)
1963	self._draw_model(update_chisqr=False)
1964	#PostStatusEvent
1965	#msg = "Fit completed!dddd "
1966	#wx.PostEvent(self._manager.parent, StatusEvent(status=msg))
1967
1968	def onPinholeSmear(self, event):
1969	"""
1970	Create a custom pinhole smear object that will change the way residuals
1971	are compute when fitting
1972
1973	:Note: accuracy is given by strings'High','Med', 'Low' FOR 2d,
1974	None for 1D
1975
1976	"""
1977
1978	#if self.check_invalid_panel():
1979	# return
1980	if self.model ==None:
1981	self.disable_smearer.SetValue(True)
1982	msg="Please select a Model first..."
1983	wx.MessageBox(msg, 'Info')
1984	wx.PostEvent(self._manager.parent, StatusEvent(status=\
1985	"Smear: %s"%msg))
1986	return
1987
1988	# Need update param values
1989	self._update_paramv_on_fit()
1990
1991	# msg default
1992	msg = None
1993	if event != None:
1994	tcrtl= event.GetEventObject()
1995	# event case of radio button
1996	if tcrtl.GetValue()== True:
1997	self.dx_min = 0.0
1998	self.dx_max = 0.0
1999	is_new_pinhole = True
2000	else:
2001	is_new_pinhole = self._is_changed_pinhole()
2002	else:
2003	is_new_pinhole = True
2004	# if any value is changed
2005	if is_new_pinhole:
2006	msg = self._set_pinhole_smear()
2007	# hide all silt sizer
2008	self._hide_all_smear_info()
2009
2010	##Calculate chi2
2011	temp_smearer = self.current_smearer
2012	#self.compute_chisqr(smearer= temp_smearer)
2013
2014	# show relevant slit sizers
2015	self._show_smear_sizer()
2016
2017	self.sizer_set_smearer.Layout()
2018	self.Layout()
2019
2020	if event != None:
2021	event.Skip()
2022	#self._undo.Enable(True)
2023	self.save_current_state()
2024	event = PageInfoEvent(page = self)
2025	wx.PostEvent(self.parent, event)
2026
2027	def _is_changed_pinhole(self):
2028	"""
2029	check if any of pinhole smear is changed
2030
2031	:return: True or False
2032
2033	"""
2034	# get the values
2035	pin_min = self.smear_pinhole_min.GetValue()
2036	pin_max = self.smear_pinhole_max.GetValue()
2037
2038	# Check changes in slit width
2039	try:
2040	dx_min = float(pin_min)
2041	except:
2042	return True
2043	if self.dx_min != dx_min:
2044	return True
2045
2046	# Check changes in slit heigth
2047	try:
2048	dx_max = float(pin_max)
2049	except:
2050	return True
2051	if self.dx_max != dx_max:
2052	return True
2053	return False
2054
2055	def _set_pinhole_smear(self):
2056	"""
2057	Set custom pinhole smear
2058
2059	:return: msg
2060
2061	"""
2062	# copy data
2063	data = copy.deepcopy(self.data)
2064	if self._is_2D():
2065	self.smear_type = 'Pinhole2d'
2066	len_data = len(data.data)
2067	data.dqx_data = numpy.zeros(len_data)
2068	data.dqy_data = numpy.zeros(len_data)
2069	else:
2070	self.smear_type = 'Pinhole'
2071	len_data = len(data.x)
2072	data.dx = numpy.zeros(len_data)
2073	data.dxl = None
2074	data.dxw = None
2075	msg = None
2076
2077	get_pin_min = self.smear_pinhole_min
2078	get_pin_max = self.smear_pinhole_max
2079
2080	if not check_float(get_pin_min):
2081	get_pin_min.SetBackgroundColour("pink")
2082	msg= "Model Error:wrong value entered!!!"
2083	elif not check_float(get_pin_max ):
2084	get_pin_max.SetBackgroundColour("pink")
2085	msg= "Model Error:wrong value entered!!!"
2086	else:
2087	if len_data < 2: len_data = 2
2088	self.dx_min = float(get_pin_min.GetValue())
2089	self.dx_max = float(get_pin_max.GetValue())
2090	if self.dx_min < 0:
2091	get_pin_min.SetBackgroundColour("pink")
2092	msg= "Model Error:This value can not be negative!!!"
2093	elif self.dx_max <0:
2094	get_pin_max.SetBackgroundColour("pink")
2095	msg= "Model Error:This value can not be negative!!!"
2096	elif self.dx_min != None and self.dx_max != None:
2097	if self._is_2D():
2098	data.dqx_data[data.dqx_data==0] = self.dx_min
2099	data.dqy_data[data.dqy_data==0] = self.dx_max
2100	elif self.dx_min == self.dx_max:
2101	data.dx[data.dx==0] = self.dx_min
2102	else:
2103	step = (self.dx_max - self.dx_min)/(len_data-1)
2104	data.dx = numpy.arange(self.dx_min,self.dx_max+step/1.1,
2105	step)
2106	elif self.dx_min != None:
2107	if self._is_2D(): data.dqx_data[data.dqx_data==0] = self.dx_min
2108	else: data.dx[data.dx==0] = self.dx_min
2109	elif self.dx_max != None:
2110	if self._is_2D(): data.dqy_data[data.dqy_data==0] = self.dx_max
2111	else: data.dx[data.dx==0] = self.dx_max
2112	self.current_smearer = smear_selection(data, self.model)
2113	# 2D need to set accuracy
2114	if self._is_2D():
2115	self.current_smearer.set_accuracy(accuracy = \
2116	self.smear2d_accuracy)
2117
2118	if msg != None:
2119	wx.PostEvent(self._manager.parent, StatusEvent(status = msg ))
2120	else:
2121	get_pin_min.SetBackgroundColour("white")
2122	get_pin_max.SetBackgroundColour("white")
2123	## set smearing value whether or not the data contain the smearing info
2124	self._manager.set_smearer(smearer=self.current_smearer,
2125	enable2D=self.enable2D,
2126	qmin=float(self.qmin_x),
2127	qmax= float(self.qmax_x),
2128	uid=self.uid)
2129	return msg
2130
2131	def update_pinhole_smear(self):
2132	"""
2133	called by kill_focus on pinhole TextCntrl
2134	to update the changes
2135
2136	:return: False when wrong value was entered
2137
2138	"""
2139	# msg default
2140	msg = None
2141	# check if any value is changed
2142	if self._is_changed_pinhole():
2143	msg = self._set_pinhole_smear()
2144	#self._undo.Enable(True)
2145	self.save_current_state()
2146
2147	if msg != None:
2148	return False
2149	else:
2150	return True
2151
2152	def onSlitSmear(self, event):
2153	"""
2154	Create a custom slit smear object that will change the way residuals
2155	are compute when fitting
2156	"""
2157
2158	#if self.check_invalid_panel():
2159	# return
2160
2161	if self.model ==None:
2162	self.disable_smearer.SetValue(True)
2163
2164	msg="Please select a Model first..."
2165	wx.MessageBox(msg, 'Info')
2166	wx.PostEvent(self._manager.parent, StatusEvent(status=\
2167	"Smear: %s"%msg))
2168	return
2169
2170	# Need update param values
2171	self._update_paramv_on_fit()
2172
2173	# msg default
2174	msg = None
2175	# for event given
2176	if event != None:
2177	tcrtl= event.GetEventObject()
2178	# event case of radio button
2179	if tcrtl.GetValue():
2180	self.dxl = 0.0
2181	self.dxw = 0.0
2182	is_new_slit = True
2183	else:
2184	is_new_slit = self._is_changed_slit()
2185	else:
2186	is_new_slit = True
2187
2188	# if any value is changed
2189	if is_new_slit:
2190	msg = self._set_slit_smear()
2191
2192	# hide all silt sizer
2193	self._hide_all_smear_info()
2194	##Calculate chi2
2195	#self.compute_chisqr(smearer= self.current_smearer)
2196	# show relevant slit sizers
2197	self._show_smear_sizer()
2198	self.sizer_set_smearer.Layout()
2199	self.Layout()
2200
2201	if event != None:
2202	event.Skip()
2203	#self._undo.Enable(True)
2204	self.save_current_state()
2205	event = PageInfoEvent(page = self)
2206	wx.PostEvent(self.parent, event)
2207	if msg != None:
2208	wx.PostEvent(self._manager.parent, StatusEvent(status = msg))
2209
2210	def _is_changed_slit(self):
2211	"""
2212	check if any of slit lengths is changed
2213
2214	:return: True or False
2215
2216	"""
2217	# get the values
2218	width = self.smear_slit_width.GetValue()
2219	height = self.smear_slit_height.GetValue()
2220
2221	# check and change the box bg color if it was pink
2222	# but it should be white now
2223	# because this is the case that _set_slit_smear() will not handle
2224	if height.lstrip().rstrip()=="":
2225	self.smear_slit_height.SetBackgroundColour(wx.WHITE)
2226	if width.lstrip().rstrip()=="":
2227	self.smear_slit_width.SetBackgroundColour(wx.WHITE)
2228
2229	# Check changes in slit width
2230	if width == "":
2231	dxw = 0.0
2232	else:
2233	try:
2234	dxw = float(width)
2235	except:
2236	return True
2237	if self.dxw != dxw:
2238	return True
2239
2240	# Check changes in slit heigth
2241	if height == "":
2242	dxl = 0.0
2243	else:
2244	try:
2245	dxl = float(height)
2246	except:
2247	return True
2248	if self.dxl != dxl:
2249	return True
2250
2251	return False
2252
2253	def _set_slit_smear(self):
2254	"""
2255	Set custom slit smear
2256
2257	:return: message to inform the user about the validity
2258	of the values entered for slit smear
2259	"""
2260	if self.data.__class__.__name__ == "Data2D" or \
2261	self.enable2D:
2262	return
2263	temp_smearer = None
2264	# make sure once more if it is smearer
2265	data = copy.deepcopy(self.data)
2266	data_len = len(data.x)
2267	data.dx = None
2268	data.dxl = None
2269	data.dxw = None
2270	msg = None
2271
2272	try:
2273	self.dxl = float(self.smear_slit_height.GetValue())
2274	data.dxl = self.dxl* numpy.ones(data_len)
2275	self.smear_slit_height.SetBackgroundColour(wx.WHITE)
2276	except:
2277	data.dxl = numpy.zeros(data_len)
2278	if self.smear_slit_height.GetValue().lstrip().rstrip()!="":
2279	self.smear_slit_height.SetBackgroundColour("pink")
2280	msg = "Wrong value entered... "
2281	else:
2282	self.smear_slit_height.SetBackgroundColour(wx.WHITE)
2283	try:
2284	self.dxw = float(self.smear_slit_width.GetValue())
2285	self.smear_slit_width.SetBackgroundColour(wx.WHITE)
2286	data.dxw = self.dxw* numpy.ones(data_len)
2287	except:
2288	data.dxw = numpy.zeros(data_len)
2289	if self.smear_slit_width.GetValue().lstrip().rstrip()!="":
2290	self.smear_slit_width.SetBackgroundColour("pink")
2291	msg = "Wrong Fit value entered... "
2292	else:
2293	self.smear_slit_width.SetBackgroundColour(wx.WHITE)
2294
2295	self.current_smearer = smear_selection(data, self.model)
2296	#temp_smearer = self.current_smearer
2297	## set smearing value whether or not the data contain the smearing info
2298	self._manager.set_smearer(smearer=self.current_smearer,
2299	enable2D=self.enable2D,
2300	qmin=float(self.qmin_x),
2301	qmax= float(self.qmax_x),
2302	uid=self.uid)
2303	return msg
2304
2305	def update_slit_smear(self):
2306	"""
2307	called by kill_focus on pinhole TextCntrl
2308	to update the changes
2309
2310	:return: False when wrong value was entered
2311
2312	"""
2313	# msg default
2314	msg = None
2315	# check if any value is changed
2316	if self._is_changed_slit():
2317	msg = self._set_slit_smear()
2318	#self._undo.Enable(True)
2319	self.save_current_state()
2320
2321	if msg != None:
2322	return False
2323	else:
2324	return True
2325
2326	def onSmear(self, event):
2327	"""
2328	Create a smear object that will change the way residuals
2329	are compute when fitting
2330	"""
2331	if event != None:
2332	event.Skip()
2333	if self.data is None:
2334	return
2335
2336	if self.model == None:
2337	self.disable_smearer.SetValue(True)
2338	msg="Please select a Model first..."
2339	wx.MessageBox(msg, 'Info')
2340	wx.PostEvent(self._manager.parent, StatusEvent(status=\
2341	"Smear: %s"%msg))
2342	return
2343
2344	# Need update param values
2345	self._update_paramv_on_fit()
2346
2347	temp_smearer = None
2348	self._get_smear_info()
2349
2350	#renew smear sizer
2351	if self.smear_type != None:
2352	self.smear_description_smear_type.SetValue(str(self.smear_type))
2353	self.smear_data_left.SetValue(str(self.dq_l))
2354	self.smear_data_right.SetValue(str(self.dq_r))
2355
2356	self._hide_all_smear_info()
2357
2358	data = copy.deepcopy(self.data)
2359	# make sure once more if it is smearer
2360	self.current_smearer = smear_selection(data, self.model)
2361
2362	if self.enable_smearer.GetValue():
2363	if hasattr(self.data,"dxl"):
2364
2365	msg= ": Resolution smearing parameters"
2366	if hasattr(self.data,"dxw"):
2367	msg= ": Slit smearing parameters"
2368	if self.smearer ==None:
2369	wx.PostEvent(self._manager.parent, StatusEvent(status=\
2370	"Data contains no smearing information"))
2371	else:
2372	wx.PostEvent(self._manager.parent, StatusEvent(status=\
2373	"Data contains smearing information"))
2374
2375	#self.smear_description_dqdata.Show(True)
2376	self.smear_data_left.Show(True)
2377	self.smear_data_right.Show(True)
2378	temp_smearer= self.current_smearer
2379	elif self.disable_smearer.GetValue():
2380	self.smear_description_none.Show(True)
2381
2382	self._show_smear_sizer()
2383
2384	self.sizer_set_smearer.Layout()
2385	self.Layout()
2386	## set smearing value whether or not the data contain the smearing info
2387	self._manager.set_smearer(uid=self.uid, smearer=temp_smearer,
2388	enable2D=self.enable2D,
2389	qmin= float(self.qmin_x),
2390	qmax= float(self.qmax_x), draw=True)
2391
2392	self.state.enable_smearer= self.enable_smearer.GetValue()
2393	self.state.disable_smearer=self.disable_smearer.GetValue()
2394	self.state.pinhole_smearer = self.pinhole_smearer.GetValue()
2395	self.state.slit_smearer = self.slit_smearer.GetValue()
2396
2397	def on_complete_chisqr(self, event):
2398	"""
2399	print result chisqr
2400
2401	:event: activated by fitting/ complete after draw
2402
2403	"""
2404	try:
2405	if event == None:
2406	output= "-"
2407	elif not numpy.isfinite(event.output):
2408	output= "-"
2409	else:
2410	output = event.output
2411	self.tcChi.SetValue(str(format_number(output, True)))
2412
2413	self.state.tcChi = self.tcChi.GetValue()
2414	except:
2415	pass
2416
2417
2418	def select_all_param(self,event):
2419	"""
2420	set to true or false all checkBox given the main checkbox value cb1
2421	"""
2422	self.param_toFit=[]
2423	if self.parameters !=[]:
2424	if self.cb1.GetValue():
2425	for item in self.parameters:
2426	## for data2D select all to fit
2427	if self.data.__class__.__name__== "Data2D" or \
2428	self.enable2D:
2429	item[0].SetValue(True)
2430	self.param_toFit.append(item )
2431	else:
2432	## for 1D all parameters except orientation
2433	if not item in self.orientation_params:
2434	item[0].SetValue(True)
2435	self.param_toFit.append(item )
2436	#if len(self.fittable_param)>0:
2437	for item in self.fittable_param:
2438	if self.data.__class__.__name__== "Data2D" or \
2439	self.enable2D:
2440	item[0].SetValue(True)
2441	self.param_toFit.append(item )
2442	try:
2443	if len(self.values[item[1]]) > 0:
2444	item[0].SetValue(False)
2445	except:
2446	pass
2447
2448	else:
2449	## for 1D all parameters except orientation
2450	if not item in self.orientation_params_disp:
2451	item[0].SetValue(True)
2452	self.param_toFit.append(item )
2453	try:
2454	if len(self.values[item[1]]) > 0:
2455	item[0].SetValue(False)
2456	except:
2457	pass
2458
2459	else:
2460	for item in self.parameters:
2461	item[0].SetValue(False)
2462	for item in self.fittable_param:
2463	item[0].SetValue(False)
2464	self.param_toFit=[]
2465
2466	self.save_current_state_fit()
2467
2468	if event !=None:
2469	#self._undo.Enable(True)
2470	## post state to fit panel
2471	event = PageInfoEvent(page = self)
2472	wx.PostEvent(self.parent, event)
2473
2474	def select_param(self,event):
2475	"""
2476	Select TextCtrl checked for fitting purpose and stores them
2477	in self.param_toFit=[] list
2478	"""
2479	self.param_toFit=[]
2480	for item in self.parameters:
2481	#Skip t ifhe angle parameters if 1D data
2482	if self.data.__class__.__name__ != "Data2D" and\
2483	not self.enable2D:
2484	if item in self.orientation_params:
2485	continue
2486	#Select parameters to fit for list of primary parameters
2487	if item[0].GetValue():
2488	if not (item in self.param_toFit):
2489	self.param_toFit.append(item )
2490	else:
2491	#remove parameters from the fitting list
2492	if item in self.param_toFit:
2493	self.param_toFit.remove(item)
2494
2495	#Select parameters to fit for list of fittable parameters
2496	# with dispersion
2497	for item in self.fittable_param:
2498	#Skip t ifhe angle parameters if 1D data
2499	if self.data.__class__.__name__ != "Data2D" and\
2500	not self.enable2D:
2501	if item in self.orientation_params:
2502	continue
2503	if item[0].GetValue():
2504	if not (item in self.param_toFit):
2505	self.param_toFit.append(item)
2506	else:
2507	#remove parameters from the fitting list
2508	if item in self.param_toFit:
2509	self.param_toFit.remove(item)
2510
2511	#Calculate num. of angle parameters
2512	if self.data.__class__.__name__ == "Data2D" or \
2513	self.enable2D:
2514	len_orient_para = 0
2515	else:
2516	len_orient_para = len(self.orientation_params) #assume even len
2517	#Total num. of angle parameters
2518	if len(self.fittable_param) > 0:
2519	len_orient_para *= 2
2520	#Set the value of checkbox that selected every checkbox or not
2521	if len(self.parameters)+len(self.fittable_param)-len_orient_para ==\
2522	len(self.param_toFit):
2523	self.cb1.SetValue(True)
2524	else:
2525	self.cb1.SetValue(False)
2526	self.save_current_state_fit()
2527	if event !=None:
2528	#self._undo.Enable(True)
2529	## post state to fit panel
2530	event = PageInfoEvent(page = self)
2531	wx.PostEvent(self.parent, event)
2532
2533	def set_model_param_sizer(self, model):
2534	"""
2535	Build the panel from the model content
2536
2537	:param model: the model selected in combo box for fitting purpose
2538
2539	"""
2540	self.sizer3.Clear(True)
2541	self.parameters = []
2542	self.str_parameters = []
2543	self.param_toFit=[]
2544	self.fittable_param=[]
2545	self.fixed_param=[]
2546	self.orientation_params=[]
2547	self.orientation_params_disp=[]
2548
2549	if model ==None:
2550	self.sizer3.Layout()
2551	self.SetupScrolling()
2552	return
2553	## the panel is drawn using the current value of the fit engine
2554	if self.engine_type==None and self._manager !=None:
2555	self.engine_type= self._manager._return_engine_type()
2556
2557	box_description= wx.StaticBox(self, -1,str("Model Parameters"))
2558	boxsizer1 = wx.StaticBoxSizer(box_description, wx.VERTICAL)
2559	sizer = wx.GridBagSizer(5,5)
2560	## save the current model
2561	self.model = model
2562
2563	keys = self.model.getParamList()
2564	#list of dispersion paramaters
2565	self.disp_list=self.model.getDispParamList()
2566
2567	def custom_compare(a,b):
2568	"""
2569	Custom compare to order, first by alphabets then second by number.
2570	"""
2571	# number at the last digit
2572	a_last = a[len(a)-1]
2573	b_last = b[len(b)-1]
2574	# default
2575	num_a = None
2576	num_b = None
2577	# split the names
2578	a2 = a.lower().split('_')
2579	b2 = b.lower().split('_')
2580	# check length of a2, b2
2581	len_a2 = len(a2)
2582	len_b2 = len(b2)
2583	# check if it contains a int number(<10)
2584	try:
2585	num_a = int(a_last)
2586	except: pass
2587	try:
2588	num_b = int(b_last)
2589	except: pass
2590	# Put 'scale' near the top; happens
2591	# when numbered param name exists
2592	if a == 'scale':
2593	return -1
2594	# both have a number
2595	if num_a != None and num_b != None:
2596	if num_a > num_b: return -1
2597	# same number
2598	elif num_a == num_b:
2599	# different last names
2600	if a2[len_a2-1] != b2[len_b2-1] and num_a != 0:
2601	return -cmp(a2[len_a2-1], b2[len_b2-1])
2602	else:
2603	return cmp(a, b)
2604	else: return 1
2605	# one of them has a number
2606	elif num_a != None: return 1
2607	elif num_b != None: return -1
2608	# no numbers
2609	else: return cmp(a.lower(), b.lower())
2610
2611
2612	keys.sort(custom_compare)
2613
2614	iy = 0
2615	ix = 0
2616	select_text = "Select All"
2617	self.cb1 = wx.CheckBox(self, -1,str(select_text), (10, 10))
2618	wx.EVT_CHECKBOX(self, self.cb1.GetId(), self.select_all_param)
2619	self.cb1.SetToolTipString("To check/uncheck all the boxes below.")
2620	self.cb1.SetValue(True)
2621
2622	sizer.Add(self.cb1,(iy, ix),(1,1),\
2623	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 5)
2624	ix +=1
2625	self.text2_2 = wx.StaticText(self, -1, 'Value')
2626	sizer.Add(self.text2_2,(iy, ix),(1,1),\
2627	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2628	ix +=2
2629	self.text2_3 = wx.StaticText(self, -1, 'Error')
2630	sizer.Add(self.text2_3,(iy, ix),(1,1),\
2631	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2632	if not self.is_mac:
2633	self.text2_3.Hide()
2634	ix +=1
2635	self.text2_min = wx.StaticText(self, -1, 'Min')
2636	sizer.Add(self.text2_min,(iy, ix),(1,1),\
2637	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2638	#self.text2_min.Hide()
2639	ix +=1
2640	self.text2_max = wx.StaticText(self, -1, 'Max')
2641	sizer.Add(self.text2_max,(iy, ix),(1,1),\
2642	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2643	#self.text2_max.Hide()
2644	ix += 1
2645	self.text2_4 = wx.StaticText(self, -1, '[Units]')
2646	sizer.Add(self.text2_4,(iy, ix),(1,1),\
2647	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2648	self.text2_4.Hide()
2649
2650	for item in keys:
2651	if not item in self.disp_list and not item in \
2652	self.model.orientation_params:
2653
2654	##prepare a spot to store errors
2655	if not self.model.details.has_key(item):
2656	self.model.details [item] = ["",None,None]
2657
2658	iy += 1
2659	ix = 0
2660	if (self.model.__class__ in \
2661	self.model_list_box["Multi-Functions"] or \
2662	self.temp_multi_functional)\
2663	and (item in self.model.non_fittable):
2664	non_fittable_name = wx.StaticText(self, -1, item )
2665	sizer.Add(non_fittable_name,(iy, ix),(1,1),\
2666	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 21)
2667	## add parameter value
2668	ix += 1
2669	value= self.model.getParam(item)
2670	if len(self.model.fun_list) > 0:
2671	num = item.split('_')[1][5:7]
2672	fun_box = wx.ComboBox(self, -1,size=(100,-1),
2673	style=wx.CB_READONLY, name = '%s'% item)
2674	self._set_fun_box_list(fun_box)
2675	fun_box.SetSelection(0)
2676	#self.fun_box.SetToolTipString("A function
2677	# describing the interface")
2678	wx.EVT_COMBOBOX(fun_box,-1, self._on_fun_box)
2679	else:
2680	fun_box = self.ModelTextCtrl(self, -1,
2681	size=(_BOX_WIDTH,20),
2682	style=wx.TE_PROCESS_ENTER, name ='%s'% item)
2683	fun_box.SetToolTipString("Hit 'Enter' after typing.")
2684	fun_box.SetValue(format_number(value, True))
2685	sizer.Add(fun_box, (iy,ix),(1,1), wx.EXPAND)
2686	self.str_parameters.append([None,item, fun_box,
2687	None, None, None,
2688	None, None])
2689
2690
2691	else:
2692	## add parameters name with checkbox for selecting to fit
2693	cb = wx.CheckBox(self, -1, item )
2694	cb.SetToolTipString(" Check mark to fit.")
2695	#cb.SetValue(True)
2696	wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
2697
2698	sizer.Add( cb,( iy, ix),(1,1),
2699	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 5)
2700
2701	## add parameter value
2702	ix += 1
2703	value= self.model.getParam(item)
2704	ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH,20),
2705	style=wx.TE_PROCESS_ENTER)
2706	ctl1.SetToolTipString("Hit 'Enter' after typing.")
2707	ctl1.SetValue(format_number(value, True))
2708	sizer.Add(ctl1, (iy,ix),(1,1), wx.EXPAND)
2709	## text to show error sign
2710	ix += 1
2711	text2=wx.StaticText(self, -1, '+/-')
2712	sizer.Add(text2,(iy, ix),(1,1),\
2713	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2714	if not self.is_mac:
2715	text2.Hide()
2716	ix += 1
2717	ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.2,20), style=0)
2718	sizer.Add(ctl2, (iy,ix),(1,1),
2719	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2720	if not self.is_mac:
2721	ctl2.Hide()
2722
2723	ix += 1
2724	ctl3 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.9,20),
2725	style=wx.TE_PROCESS_ENTER,
2726	text_enter_callback = self._onparamRangeEnter)
2727
2728	sizer.Add(ctl3, (iy,ix),(1,1),
2729	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2730
2731	ix += 1
2732	ctl4 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.9,20),
2733	style=wx.TE_PROCESS_ENTER,
2734	text_enter_callback = self._onparamRangeEnter)
2735	sizer.Add(ctl4, (iy,ix),(1,1),
2736	wx.EXPAND\|wx.FIXED_MINSIZE, 0)
2737
2738	ix +=1
2739	# Units
2740	if self.model.details.has_key(item):
2741	units = wx.StaticText(self, -1,
2742	self.model.details[item][0], style=wx.ALIGN_LEFT)
2743	else:
2744	units = wx.StaticText(self, -1, "", style=wx.ALIGN_LEFT)
2745	sizer.Add(units, (iy,ix),(1,1),
2746	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2747
2748	##[cb state, name, value, "+/-", error of fit, min, max , units]
2749	self.parameters.append([cb,item, ctl1,
2750	text2,ctl2, ctl3, ctl4,units])
2751
2752	iy+=1
2753	sizer.Add((10,10),(iy,ix),(1,1),
2754	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
2755
2756	# type can be either Guassian or Array
2757	if len(self.model.dispersion.values())>0:
2758	type= self.model.dispersion.values()[0]["type"]
2759	else:
2760	type = "Gaussian"
2761
2762	iy += 1
2763	ix = 0
2764	#Add tile for orientational angle
2765	for item in keys:
2766	if item in self.model.orientation_params:
2767	orient_angle = wx.StaticText(self, -1, '[For 2D only]:')
2768	sizer.Add(orient_angle,(iy, ix),(1,1),
2769	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 15)
2770	if not self.data.__class__.__name__ == "Data2D" and \
2771	not self.enable2D:
2772	orient_angle.Hide()
2773	else:
2774	orient_angle.Show(True)
2775	break
2776
2777	#For Gaussian only
2778	if type.lower() != "array":
2779	for item in self.model.orientation_params:
2780	if not item in self.disp_list:
2781	##prepare a spot to store min max
2782	if not self.model.details.has_key(item):
2783	self.model.details [item] = ["",None,None]
2784
2785	iy += 1
2786	ix = 0
2787	## add parameters name with checkbox for selecting to fit
2788	cb = wx.CheckBox(self, -1, item )
2789	cb.SetValue(False)
2790	cb.SetToolTipString("Check mark to fit")
2791	wx.EVT_CHECKBOX(self, cb.GetId(), self.select_param)
2792	if self.data.__class__.__name__ == "Data2D" or \
2793	self.enable2D:
2794	cb.Show(True)
2795	else:
2796	cb.Hide()
2797	sizer.Add( cb,( iy, ix),(1,1),
2798	wx.LEFT\|wx.EXPAND\|wx.ADJUST_MINSIZE, 5)
2799
2800	## add parameter value
2801	ix += 1
2802	value= self.model.getParam(item)
2803	ctl1 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH,20),
2804	style=wx.TE_PROCESS_ENTER)
2805	ctl1.SetToolTipString("Hit 'Enter' after typing.")
2806	ctl1.SetValue(format_number(value, True))
2807	if self.data.__class__.__name__ == "Data2D" or \
2808	self.enable2D:
2809	ctl1.Show(True)
2810	else:
2811	ctl1.Hide()
2812	sizer.Add(ctl1, (iy,ix),(1,1), wx.EXPAND)
2813	## text to show error sign
2814	ix += 1
2815	text2=wx.StaticText(self, -1, '+/-')
2816	sizer.Add(text2,(iy, ix),(1,1),\
2817	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2818
2819	text2.Hide()
2820	ix += 1
2821	ctl2 = wx.TextCtrl(self, -1, size=(_BOX_WIDTH/1.2,20), style=0)
2822	sizer.Add(ctl2, (iy,ix),(1,1),
2823	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2824
2825	ctl2.Hide()
2826
2827
2828	ix += 1
2829	ctl3 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.8,20),
2830	style=wx.TE_PROCESS_ENTER,
2831	text_enter_callback = self._onparamRangeEnter)
2832
2833	sizer.Add(ctl3, (iy,ix),(1,1),
2834	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2835	ctl3.Hide()
2836
2837	ix += 1
2838	ctl4 = self.ModelTextCtrl(self, -1, size=(_BOX_WIDTH/1.8,20),
2839	style=wx.TE_PROCESS_ENTER,
2840	text_enter_callback = self._onparamRangeEnter)
2841	sizer.Add(ctl4, (iy,ix),(1,1),
2842	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2843
2844	ctl4.Hide()
2845
2846	if self.data.__class__.__name__ == "Data2D" or \
2847	self.enable2D:
2848	if self.is_mac:
2849	text2.Show(True)
2850	ctl2.Show(True)
2851	ctl3.Show(True)
2852	ctl4.Show(True)
2853
2854	ix +=1
2855	# Units
2856	if self.model.details.has_key(item):
2857	units = wx.StaticText(self, -1,
2858	self.model.details[item][0],
2859	style=wx.ALIGN_LEFT)
2860	else:
2861	units = wx.StaticText(self, -1, "", style=wx.ALIGN_LEFT)
2862	if self.data.__class__.__name__ == "Data2D" or \
2863	self.enable2D:
2864	units.Show(True)
2865
2866	else:
2867	units.Hide()
2868
2869	sizer.Add(units, (iy,ix),(1,1),
2870	wx.EXPAND\|wx.ADJUST_MINSIZE, 0)
2871
2872	##[cb state, name, value, "+/-", error of fit,min,max,units]
2873	self.parameters.append([cb,item, ctl1,
2874	text2,ctl2, ctl3, ctl4,units])
2875	self.orientation_params.append([cb,item, ctl1,
2876	text2,ctl2, ctl3, ctl4,units])
2877
2878	iy+=1
2879
2880	#Display units text on panel
2881	for item in keys:
2882	if self.model.details.has_key(item):
2883	self.text2_4.Show()
2884	#Fill the list of fittable parameters
2885	self.select_all_param(event=None)
2886
2887	self.save_current_state_fit()
2888	boxsizer1.Add(sizer)
2889	self.sizer3.Add(boxsizer1,0, wx.EXPAND \| wx.ALL, 10)
2890	self.sizer3.Layout()
2891	self.Layout()
2892	#self.Refresh()
2893
2894	def on_right_down(self, event):
2895	"""
2896	Get key stroke event
2897	"""
2898	if self.data == None:
2899	return
2900	# Figuring out key combo: Cmd for copy, Alt for paste
2901	if event.AltDown() and event.ShiftDown():
2902	self._manager.show_ftol_dialog()
2903	flag = True
2904	elif event.AltDown() or event.ShiftDown():
2905	flag = False
2906	else:
2907	return
2908	# make event free
2909	event.Skip()
2910	# messages depending on the flag
2911	if not flag:
2912	msg = " Could not open ftol dialog;"
2913	msg += " Check if the Scipy fit engine is selected in the menubar."
2914	infor = 'warning'
2915	# inform msg to wx
2916	wx.PostEvent( self.parent.parent,
2917	StatusEvent(status= msg, info=infor))
2918
2919
2920	def _onModel2D(self, event):
2921	"""
2922	toggle view of model from 1D to 2D or 2D from 1D
2923	"""
2924	if self.model_view.GetLabelText() == "1D Mode":
2925	self.model_view.SetLabel("2D Mode")
2926	self.enable2D = True
2927
2928	else:
2929	self.model_view.SetLabel("1D Mode")
2930	self.enable2D = False
2931
2932	self.set_model_param_sizer(self.model)
2933	self._set_sizer_dispersion()
2934	self._draw_model()
2935
2936
2937	self.state.enable2D = copy.deepcopy(self.enable2D)
2938
2939	class BGTextCtrl(wx.TextCtrl):
2940	"""
2941	Text control used to display outputs.
2942	No editing allowed. The background is
2943	grayed out. User can't select text.
2944	"""
2945	def __init__(self, args, *kwds):
2946	wx.TextCtrl.__init__(self, args, *kwds)
2947	self.SetEditable(False)
2948	self.SetBackgroundColour(self.GetParent().parent.GetBackgroundColour())
2949
2950	# Bind to mouse event to avoid text highlighting
2951	# The event will be skipped once the call-back
2952	# is called.
2953	self.Bind(wx.EVT_MOUSE_EVENTS, self._click)
2954
2955	def _click(self, event):
2956	"""
2957	Prevent further handling of the mouse event
2958	by not calling Skip().
2959	"""
2960	pass
2961

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

SasView

source: sasview/sansview/perspectives/fitting/fitpage.py @ 75daa79

Download in other formats: