fitpage.py @ d5c2f4d

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

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

SasView

source: sasview/sansview/perspectives/fitting/fitpage.py @ d5c2f4d

Download in other formats: