fitpage.py @ 1b43306

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 1b43306 was 5062bbf, checked in by Gervaise Alina <gervyh@…>, 14 years ago
working on documentation
Property mode set to `100644`
File size: 99.5 KB

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

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

SasView

source: sasview/sansview/perspectives/fitting/fitpage.py @ 1b43306

Download in other formats: