fitpage.py @ 29ef718

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

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

SasView

source: sasview/sansview/perspectives/fitting/fitpage.py @ 29ef718

Download in other formats: