fitpage.py @ 4fa5f3b

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

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

SasView

source: sasview/sansview/perspectives/fitting/fitpage.py @ 4fa5f3b

Download in other formats: