fitpage.py @ ba64348

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

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

SasView

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

Download in other formats: