fitpage.py @ 1180528

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 1180528 was 1180528, checked in by Jae Cho <jhjcho@…>, 13 years ago
fixed a minor bug in last commit
Property mode set to `100644`
File size: 120.7 KB

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

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

SasView

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

Download in other formats: