fitpage.py @ e4c9030

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

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

SasView

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

Download in other formats: