pagestate.py @ 3152a02

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalcmagnetic_scattrelease-4.2.2ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload

Last change on this file since 3152a02 was 959eb01, checked in by ajj, 8 years ago
normalising line endings
Property mode set to `100644`
File size: 56.7 KB

Line
1	"""
2	Class that holds a fit page state
3	"""
4	# TODO: Refactor code so we don't need to use getattr/setattr
5	################################################################################
6	# This software was developed by the University of Tennessee as part of the
7	# Distributed Data Analysis of Neutron Scattering Experiments (DANSE)
8	# project funded by the US National Science Foundation.
9	#
10	# See the license text in license.txt
11	#
12	# copyright 2009, University of Tennessee
13	################################################################################
14	import time
15	import os
16	import sys
17	import wx
18	import copy
19	import logging
20	import numpy as np
21	import traceback
22
23	import xml.dom.minidom
24	from xml.dom.minidom import parseString
25	from lxml import etree
26
27	from sasmodels import convert
28	import sasmodels.weights
29
30	import sas.sascalc.dataloader
31	from sas.sascalc.dataloader.readers.cansas_reader import Reader as CansasReader
32	from sas.sascalc.dataloader.readers.cansas_reader import get_content, write_node
33	from sas.sascalc.dataloader.data_info import Data2D, Collimation, Detector
34	from sas.sascalc.dataloader.data_info import Process, Aperture
35
36	logger = logging.getLogger(__name__)
37
38	# Information to read/write state as xml
39	FITTING_NODE_NAME = 'fitting_plug_in'
40	CANSAS_NS = "cansas1d/1.0"
41
42	CUSTOM_MODEL = 'Plugin Models'
43	CUSTOM_MODEL_OLD = 'Customized Models'
44
45	LIST_OF_DATA_ATTRIBUTES = [["is_data", "is_data", "bool"],
46	["group_id", "data_group_id", "string"],
47	["data_name", "data_name", "string"],
48	["data_id", "data_id", "string"],
49	["name", "name", "string"],
50	["data_name", "data_name", "string"]]
51	LIST_OF_STATE_ATTRIBUTES = [["qmin", "qmin", "float"],
52	["qmax", "qmax", "float"],
53	["npts", "npts", "float"],
54	["categorycombobox", "categorycombobox", "string"],
55	["formfactorcombobox", "formfactorcombobox",
56	"string"],
57	["structurecombobox", "structurecombobox",
58	"string"],
59	["multi_factor", "multi_factor", "float"],
60	["magnetic_on", "magnetic_on", "bool"],
61	["enable_smearer", "enable_smearer", "bool"],
62	["disable_smearer", "disable_smearer", "bool"],
63	["pinhole_smearer", "pinhole_smearer", "bool"],
64	["slit_smearer", "slit_smearer", "bool"],
65	["enable_disp", "enable_disp", "bool"],
66	["disable_disp", "disable_disp", "bool"],
67	["dI_noweight", "dI_noweight", "bool"],
68	["dI_didata", "dI_didata", "bool"],
69	["dI_sqrdata", "dI_sqrdata", "bool"],
70	["dI_idata", "dI_idata", "bool"],
71	["enable2D", "enable2D", "bool"],
72	["cb1", "cb1", "bool"],
73	["tcChi", "tcChi", "float"],
74	["smearer", "smearer", "float"],
75	["smear_type", "smear_type", "string"],
76	["dq_l", "dq_l", "float"],
77	["dq_r", "dq_r", "float"],
78	["dx_percent", "dx_percent", "float"],
79	["dxl", "dxl", "float"],
80	["dxw", "dxw", "float"]]
81
82	LIST_OF_MODEL_ATTRIBUTES = [["values", "values"],
83	["weights", "weights"]]
84
85	DISPERSION_LIST = [["disp_obj_dict", "_disp_obj_dict", "string"]]
86
87	LIST_OF_STATE_PARAMETERS = [["parameters", "parameters"],
88	["str_parameters", "str_parameters"],
89	["orientation_parameters", "orientation_params"],
90	["dispersity_parameters",
91	"orientation_params_disp"],
92	["fixed_param", "fixed_param"],
93	["fittable_param", "fittable_param"]]
94	LIST_OF_DATA_2D_ATTR = [["xmin", "xmin", "float"],
95	["xmax", "xmax", "float"],
96	["ymin", "ymin", "float"],
97	["ymax", "ymax", "float"],
98	["_xaxis", "_xaxis", "string"],
99	["_xunit", "_xunit", "string"],
100	["_yaxis", "_yaxis", "string"],
101	["_yunit", "_yunit", "string"],
102	["_zaxis", "_zaxis", "string"],
103	["_zunit", "_zunit", "string"]]
104	LIST_OF_DATA_2D_VALUES = [["qx_data", "qx_data", "float"],
105	["qy_data", "qy_data", "float"],
106	["dqx_data", "dqx_data", "float"],
107	["dqy_data", "dqy_data", "float"],
108	["data", "data", "float"],
109	["q_data", "q_data", "float"],
110	["err_data", "err_data", "float"],
111	["mask", "mask", "bool"]]
112
113
114	def parse_entry_helper(node, item):
115	"""
116	Create a numpy list from value extrated from the node
117
118	:param node: node from each the value is stored
119	:param item: list name of three strings.the two first are name of data
120	attribute and the third one is the type of the value of that
121	attribute. type can be string, float, bool, etc.
122
123	: return: numpy array
124	"""
125	if node is not None:
126	if item[2] == "string":
127	return str(node.get(item[0]).strip())
128	elif item[2] == "bool":
129	try:
130	return node.get(item[0]).strip() == "True"
131	except Exception:
132	return None
133	else:
134	try:
135	return float(node.get(item[0]))
136	except Exception:
137	return None
138
139
140	class PageState(object):
141	"""
142	Contains information to reconstruct a page of the fitpanel.
143	"""
144	def __init__(self, parent=None, model=None, data=None):
145	"""
146	Initialize the current state
147
148	:param model: a selected model within a page
149	:param data:
150
151	"""
152	self.file = None
153	# Time of state creation
154	self.timestamp = time.time()
155	# Data member to store the dispersion object created
156	self._disp_obj_dict = {}
157	# ------------------------
158	# Data used for fitting
159	self.data = data
160	# model data
161	self.theory_data = None
162	# Is 2D
163	self.is_2D = False
164	self.images = None
165
166	# save additional information on data that dataloader.reader
167	# does not read
168	self.is_data = None
169	self.data_name = ""
170
171	if self.data is not None:
172	self.data_name = self.data.name
173	self.data_id = None
174	if self.data is not None and hasattr(self.data, "id"):
175	self.data_id = self.data.id
176	self.data_group_id = None
177	if self.data is not None and hasattr(self.data, "group_id"):
178	self.data_group_id = self.data.group_id
179
180	# reset True change the state of existing button
181	self.reset = False
182
183	# flag to allow data2D plot
184	self.enable2D = False
185	# model on which the fit would be performed
186	self.model = model
187	self.m_name = None
188	# list of process done to model
189	self.process = []
190	# fit page manager
191	self.manager = None
192	# Store the parent of this panel parent
193	# For this application fitpanel is the parent
194	self.parent = parent
195	# Event_owner is the owner of model event
196	self.event_owner = None
197	# page name
198	self.page_name = ""
199	# Contains link between model, its parameters, and panel organization
200	self.parameters = []
201	# String parameter list that can not be fitted
202	self.str_parameters = []
203	# Contains list of parameters that cannot be fitted and reference to
204	# panel objects
205	self.fixed_param = []
206	# Contains list of parameters with dispersity and reference to
207	# panel objects
208	self.fittable_param = []
209	# orientation parameters
210	self.orientation_params = []
211	# orientation parameters for gaussian dispersity
212	self.orientation_params_disp = []
213	# smearer info
214	self.smearer = None
215	self.smear_type = None
216	self.dq_l = None
217	self.dq_r = None
218	self.dx_percent = None
219	self.dx_old = False
220	self.dxl = None
221	self.dxw = None
222	# list of dispersion parameters
223	self.disp_list = []
224	if self.model is not None:
225	self.disp_list = self.model.getDispParamList()
226
227	self.disp_cb_dict = {}
228	self.values = {}
229	self.weights = {}
230
231	# contains link between a model and selected parameters to fit
232	self.param_toFit = []
233	# dictionary of model type and model class
234	self.model_list_box = None
235	# save the state of the context menu
236	self.saved_states = {}
237	# save selection of combobox
238	self.formfactorcombobox = None
239	self.categorycombobox = None
240	self.structurecombobox = None
241
242	# radio box to select type of model
243	# self.shape_rbutton = False
244	# self.shape_indep_rbutton = False
245	# self.struct_rbutton = False
246	# self.plugin_rbutton = False
247	# the indice of the current selection
248	self.disp_box = 0
249	# Qrange
250	# Q range
251	self.qmin = 0.001
252	self.qmax = 0.1
253	# reset data range
254	self.qmax_x = None
255	self.qmin_x = None
256
257	self.npts = None
258	self.name = ""
259	self.multi_factor = None
260	self.magnetic_on = False
261	# enable smearering state
262	self.enable_smearer = False
263	self.disable_smearer = True
264	self.pinhole_smearer = False
265	self.slit_smearer = False
266	# weighting options
267	self.dI_noweight = False
268	self.dI_didata = True
269	self.dI_sqrdata = False
270	self.dI_idata = False
271	# disperity selection
272	self.enable_disp = False
273	self.disable_disp = True
274
275	# state of selected all check button
276	self.cb1 = False
277	# store value of chisqr
278	self.tcChi = None
279	self.version = (1,0,0)
280
281	def clone(self):
282	"""
283	Create a new copy of the current object
284	"""
285	model = None
286	if self.model is not None:
287	model = self.model.clone()
288	model.name = self.model.name
289	obj = PageState(self.parent, model=model)
290	obj.file = copy.deepcopy(self.file)
291	obj.data = copy.deepcopy(self.data)
292	if self.data is not None:
293	self.data_name = self.data.name
294	obj.data_name = self.data_name
295	obj.is_data = self.is_data
296	obj.model_list_box = copy.deepcopy(self.model_list_box)
297
298	obj.categorycombobox = self.categorycombobox
299	obj.formfactorcombobox = self.formfactorcombobox
300	obj.structurecombobox = self.structurecombobox
301
302	# obj.shape_rbutton = self.shape_rbutton
303	# obj.shape_indep_rbutton = self.shape_indep_rbutton
304	# obj.struct_rbutton = self.struct_rbutton
305	# obj.plugin_rbutton = self.plugin_rbutton
306
307	obj.manager = self.manager
308	obj.event_owner = self.event_owner
309	obj.disp_list = copy.deepcopy(self.disp_list)
310
311	obj.enable2D = copy.deepcopy(self.enable2D)
312	obj.parameters = copy.deepcopy(self.parameters)
313	obj.str_parameters = copy.deepcopy(self.str_parameters)
314	obj.fixed_param = copy.deepcopy(self.fixed_param)
315	obj.fittable_param = copy.deepcopy(self.fittable_param)
316	obj.orientation_params = copy.deepcopy(self.orientation_params)
317	obj.orientation_params_disp = \
318	copy.deepcopy(self.orientation_params_disp)
319	obj.enable_disp = copy.deepcopy(self.enable_disp)
320	obj.disable_disp = copy.deepcopy(self.disable_disp)
321	obj.tcChi = self.tcChi
322
323	if len(self._disp_obj_dict) > 0:
324	for k, v in self._disp_obj_dict.iteritems():
325	obj._disp_obj_dict[k] = v
326	if len(self.disp_cb_dict) > 0:
327	for k, v in self.disp_cb_dict.iteritems():
328	obj.disp_cb_dict[k] = v
329	if len(self.values) > 0:
330	for k, v in self.values.iteritems():
331	obj.values[k] = v
332	if len(self.weights) > 0:
333	for k, v in self.weights.iteritems():
334	obj.weights[k] = v
335	obj.enable_smearer = copy.deepcopy(self.enable_smearer)
336	obj.disable_smearer = copy.deepcopy(self.disable_smearer)
337	obj.pinhole_smearer = copy.deepcopy(self.pinhole_smearer)
338	obj.slit_smearer = copy.deepcopy(self.slit_smearer)
339	obj.smear_type = copy.deepcopy(self.smear_type)
340	obj.dI_noweight = copy.deepcopy(self.dI_noweight)
341	obj.dI_didata = copy.deepcopy(self.dI_didata)
342	obj.dI_sqrdata = copy.deepcopy(self.dI_sqrdata)
343	obj.dI_idata = copy.deepcopy(self.dI_idata)
344	obj.dq_l = copy.deepcopy(self.dq_l)
345	obj.dq_r = copy.deepcopy(self.dq_r)
346	obj.dx_percent = copy.deepcopy(self.dx_percent)
347	obj.dx_old = copy.deepcopy(self.dx_old)
348	obj.dxl = copy.deepcopy(self.dxl)
349	obj.dxw = copy.deepcopy(self.dxw)
350	obj.disp_box = copy.deepcopy(self.disp_box)
351	obj.qmin = copy.deepcopy(self.qmin)
352	obj.qmax = copy.deepcopy(self.qmax)
353	obj.multi_factor = self.multi_factor
354	obj.magnetic_on = self.magnetic_on
355	obj.npts = copy.deepcopy(self.npts)
356	obj.cb1 = copy.deepcopy(self.cb1)
357	obj.smearer = copy.deepcopy(self.smearer)
358	obj.version = copy.deepcopy(self.version)
359
360	for name, state in self.saved_states.iteritems():
361	copy_name = copy.deepcopy(name)
362	copy_state = state.clone()
363	obj.saved_states[copy_name] = copy_state
364	return obj
365
366	def _old_first_model(self):
367	"""
368	Handle save states from 4.0.1 and before where the first item in the
369	selection boxes of category, formfactor and structurefactor were not
370	saved.
371	:return: None
372	"""
373	if self.categorycombobox == CUSTOM_MODEL_OLD:
374	self.categorycombobox = CUSTOM_MODEL
375	if self.formfactorcombobox == '':
376	FIRST_FORM = {
377	'Shapes' : 'BCCrystalModel',
378	'Uncategorized' : 'LineModel',
379	'StructureFactor' : 'HardsphereStructure',
380	'Ellipsoid' : 'core_shell_ellipsoid',
381	'Lamellae' : 'lamellar',
382	'Paracrystal' : 'bcc_paracrystal',
383	'Parallelepiped' : 'core_shell_parallelepiped',
384	'Shape Independent' : 'be_polyelectrolyte',
385	'Sphere' : 'adsorbed_layer',
386	'Structure Factor' : 'hardsphere',
387	CUSTOM_MODEL : ''
388	}
389	if self.categorycombobox == '':
390	if len(self.parameters) == 3:
391	self.categorycombobox = "Shape-Independent"
392	self.formfactorcombobox = 'PowerLawAbsModel'
393	elif len(self.parameters) == 9:
394	self.categorycombobox = 'Cylinder'
395	self.formfactorcombobox = 'barbell'
396	else:
397	msg = "Save state does not have enough information to load"
398	msg += " the all of the data."
399	logger.warning(msg=msg)
400	else:
401	self.formfactorcombobox = FIRST_FORM[self.categorycombobox]
402
403	@staticmethod
404	def param_remap_to_sasmodels_convert(params, is_string=False):
405	"""
406	Remaps the parameters for sasmodels conversion
407
408	:param params: list of parameters (likely self.parameters)
409	:return: remapped dictionary of parameters
410	"""
411	p = dict()
412	for fittable, name, value, _, uncert, lower, upper, units in params:
413	if not value:
414	value = np.nan
415	if not uncert or uncert[1] == '' or uncert[1] == 'None':
416	uncert[0] = False
417	uncert[1] = np.nan
418	if not upper or upper[1] == '' or upper[1] == 'None':
419	upper[0] = False
420	upper[1] = np.nan
421	if not lower or lower[1] == '' or lower[1] == 'None':
422	lower[0] = False
423	lower[1] = np.nan
424	if is_string:
425	p[name] = str(value)
426	else:
427	p[name] = float(value)
428	p[name + ".fittable"] = bool(fittable)
429	p[name + ".std"] = float(uncert[1])
430	p[name + ".upper"] = float(upper[1])
431	p[name + ".lower"] = float(lower[1])
432	p[name + ".units"] = units
433	return p
434
435	@staticmethod
436	def param_remap_from_sasmodels_convert(params):
437	"""
438	Converts {name : value} map back to [] param list
439	:param params: parameter map returned from sasmodels
440	:return: None
441	"""
442	p_map = []
443	for name, info in params.iteritems():
444	if ".fittable" in name or ".std" in name or ".upper" in name or \
445	".lower" in name or ".units" in name:
446	pass
447	else:
448	fittable = params.get(name + ".fittable", True)
449	std = params.get(name + ".std", '0.0')
450	upper = params.get(name + ".upper", 'inf')
451	lower = params.get(name + ".lower", '-inf')
452	units = params.get(name + ".units")
453	if std is not None and std is not np.nan:
454	std = [True, str(std)]
455	else:
456	std = [False, '']
457	if lower is not None and lower is not np.nan:
458	lower = [True, str(lower)]
459	else:
460	lower = [True, '-inf']
461	if upper is not None and upper is not np.nan:
462	upper = [True, str(upper)]
463	else:
464	upper = [True, 'inf']
465	param_list = [bool(fittable), str(name), str(info),
466	"+/-", std, lower, upper, str(units)]
467	p_map.append(param_list)
468	return p_map
469
470	def _convert_to_sasmodels(self):
471	"""
472	Convert parameters to a form usable by sasmodels converter
473
474	:return: None
475	"""
476	# Create conversion dictionary to send to sasmodels
477	self._old_first_model()
478	p = self.param_remap_to_sasmodels_convert(self.parameters)
479	structurefactor, params = convert.convert_model(self.structurecombobox,
480	p, False, self.version)
481	formfactor, params = convert.convert_model(self.formfactorcombobox,
482	params, False, self.version)
483	if len(self.str_parameters) > 0:
484	str_pars = self.param_remap_to_sasmodels_convert(
485	self.str_parameters, True)
486	formfactor, str_params = convert.convert_model(
487	self.formfactorcombobox, str_pars, False, self.version)
488	for key, value in str_params.iteritems():
489	params[key] = value
490
491	if self.formfactorcombobox == 'SphericalSLDModel':
492	self.multi_factor += 1
493	self.formfactorcombobox = formfactor
494	self.structurecombobox = structurefactor
495	self.parameters = []
496	self.parameters = self.param_remap_from_sasmodels_convert(params)
497
498	def _repr_helper(self, list, rep):
499	"""
500	Helper method to print a state
501	"""
502	for item in list:
503	rep += "parameter name: %s \n" % str(item[1])
504	rep += "value: %s\n" % str(item[2])
505	rep += "selected: %s\n" % str(item[0])
506	rep += "error displayed : %s \n" % str(item[4][0])
507	rep += "error value:%s \n" % str(item[4][1])
508	rep += "minimum displayed : %s \n" % str(item[5][0])
509	rep += "minimum value : %s \n" % str(item[5][1])
510	rep += "maximum displayed : %s \n" % str(item[6][0])
511	rep += "maximum value : %s \n" % str(item[6][1])
512	rep += "parameter unit: %s\n\n" % str(item[7])
513	return rep
514
515	def __repr__(self):
516	"""
517	output string for printing
518	"""
519	rep = "\nState name: %s\n" % self.file
520	t = time.localtime(self.timestamp)
521	time_str = time.strftime("%b %d %Y %H;%M;%S ", t)
522
523	rep += "State created: %s\n" % time_str
524	rep += "State form factor combobox selection: %s\n" % \
525	self.formfactorcombobox
526	rep += "State structure factor combobox selection: %s\n" % \
527	self.structurecombobox
528	rep += "is data : %s\n" % self.is_data
529	rep += "data's name : %s\n" % self.data_name
530	rep += "data's id : %s\n" % self.data_id
531	if self.model is not None:
532	m_name = self.model.__class__.__name__
533	if m_name == 'Model':
534	m_name = self.m_name
535	rep += "model name : %s\n" % m_name
536	else:
537	rep += "model name : None\n"
538	rep += "multi_factor : %s\n" % str(self.multi_factor)
539	rep += "magnetic_on : %s\n" % str(self.magnetic_on)
540	rep += "model type (Category) selected: %s\n" % self.categorycombobox
541	rep += "data : %s\n" % str(self.data)
542	rep += "Plotting Range: min: %s, max: %s, steps: %s\n" % \
543	(str(self.qmin), str(self.qmax), str(self.npts))
544	rep += "Dispersion selection : %s\n" % str(self.disp_box)
545	rep += "Smearing enable : %s\n" % str(self.enable_smearer)
546	rep += "Smearing disable : %s\n" % str(self.disable_smearer)
547	rep += "Pinhole smearer enable : %s\n" % str(self.pinhole_smearer)
548	rep += "Slit smearer enable : %s\n" % str(self.slit_smearer)
549	rep += "Dispersity enable : %s\n" % str(self.enable_disp)
550	rep += "Dispersity disable : %s\n" % str(self.disable_disp)
551	rep += "Slit smearer enable: %s\n" % str(self.slit_smearer)
552
553	rep += "dI_noweight : %s\n" % str(self.dI_noweight)
554	rep += "dI_didata : %s\n" % str(self.dI_didata)
555	rep += "dI_sqrdata : %s\n" % str(self.dI_sqrdata)
556	rep += "dI_idata : %s\n" % str(self.dI_idata)
557
558	rep += "2D enable : %s\n" % str(self.enable2D)
559	rep += "All parameters checkbox selected: %s\n" % self.cb1
560	rep += "Value of Chisqr : %s\n" % str(self.tcChi)
561	rep += "Smear object : %s\n" % self.smearer
562	rep += "Smear type : %s\n" % self.smear_type
563	rep += "dq_l : %s\n" % self.dq_l
564	rep += "dq_r : %s\n" % self.dq_r
565	rep += "dx_percent : %s\n" % str(self.dx_percent)
566	rep += "dxl : %s\n" % str(self.dxl)
567	rep += "dxw : %s\n" % str(self.dxw)
568	rep += "model : %s\n\n" % str(self.model)
569	temp_parameters = []
570	temp_fittable_param = []
571	if self.data.__class__.__name__ == "Data2D":
572	self.is_2D = True
573	else:
574	self.is_2D = False
575	if self.data is not None:
576	if not self.is_2D:
577	for item in self.parameters:
578	if item not in self.orientation_params:
579	temp_parameters.append(item)
580	for item in self.fittable_param:
581	if item not in self.orientation_params_disp:
582	temp_fittable_param.append(item)
583	else:
584	temp_parameters = self.parameters
585	temp_fittable_param = self.fittable_param
586
587	rep += "number parameters(self.parameters): %s\n" % \
588	len(temp_parameters)
589	rep = self._repr_helper(list=temp_parameters, rep=rep)
590	rep += "number str_parameters(self.str_parameters): %s\n" % \
591	len(self.str_parameters)
592	rep = self._repr_helper(list=self.str_parameters, rep=rep)
593	rep += "number fittable_param(self.fittable_param): %s\n" % \
594	len(temp_fittable_param)
595	rep = self._repr_helper(list=temp_fittable_param, rep=rep)
596	return rep
597
598	def set_report_string(self):
599	"""
600	Get the values (strings) from __str__ for report
601	"""
602	# Dictionary of the report strings
603	repo_time = ""
604	model_name = ""
605	title = ""
606	title_name = ""
607	file_name = ""
608	param_string = ""
609	paramval_string = ""
610	chi2_string = ""
611	q_range = ""
612	strings = self.__repr__()
613	lines = strings.split('\n')
614
615	# get all string values from __str__()
616	for line in lines:
617	value = ""
618	content = line.split(":")
619	name = content[0]
620	try:
621	value = content[1]
622	except Exception:
623	msg = "Report string expected 'name: value' but got %r" % line
624	logger.error(msg)
625	if name.count("State created"):
626	repo_time = "" + value
627	if name.count("parameter name"):
628	val_name = value.split(".")
629	if len(val_name) > 1:
630	if val_name[1].count("width"):
631	param_string += value + ','
632	else:
633	continue
634	else:
635	param_string += value + ','
636	if name == "value":
637	param_string += value + ','
638	fixed_parameter = False
639	if name == "selected":
640	if value == u' False':
641	fixed_parameter = True
642	if name == "error value":
643	if fixed_parameter:
644	param_string += '(fixed),'
645	else:
646	param_string += value + ','
647	if name == "parameter unit":
648	param_string += value + ':'
649	if name == "Value of Chisqr ":
650	chi2 = ("Chi2/Npts = " + value)
651	chi2_string = CENTRE % chi2
652	if name == "Title":
653	if len(value.strip()) == 0:
654	continue
655	title = value + " [" + repo_time + "]"
656	title_name = HEADER % title
657	if name == "data ":
658	try:
659	file_value = ("File name:" + content[2])
660	file_name = CENTRE % file_value
661	if len(title) == 0:
662	title = content[2] + " [" + repo_time + "]"
663	title_name = HEADER % title
664	except Exception:
665	msg = "While parsing 'data: ...'\n"
666	logger.error(msg + traceback.format_exc())
667	if name == "model name ":
668	try:
669	modelname = "Model name:" + content[1]
670	except:
671	modelname = "Model name:" + " NAN"
672	model_name = CENTRE % modelname
673
674	if name == "Plotting Range":
675	try:
676	q_range = content[1] + " = " + content[2] \
677	+ " = " + content[3].split(",")[0]
678	q_name = ("Q Range: " + q_range)
679	q_range = CENTRE % q_name
680	except Exception:
681	msg = "While parsing 'Plotting Range: ...'\n"
682	logger.error(msg + traceback.format_exc())
683	paramval = ""
684	for lines in param_string.split(":"):
685	line = lines.split(",")
686	if len(lines) > 0:
687	param = line[0]
688	param += " = " + line[1]
689	if len(line[2].split()) > 0 and not line[2].count("None"):
690	param += " +- " + line[2]
691	if len(line[3].split()) > 0 and not line[3].count("None"):
692	param += " " + line[3]
693	if not paramval.count(param):
694	paramval += param + "\n"
695	paramval_string += CENTRE % param + "\n"
696
697	text_string = "\n\n\n%s\n\n%s\n%s\n%s\n\n%s" % \
698	(title, file, q_name, chi2, paramval)
699
700	title_name = self._check_html_format(title_name)
701	file_name = self._check_html_format(file_name)
702	title = self._check_html_format(title)
703
704	html_string = title_name + "\n" + file_name + \
705	"\n" + model_name + \
706	"\n" + q_range + \
707	"\n" + chi2_string + \
708	"\n" + ELINE + \
709	"\n" + paramval_string + \
710	"\n" + ELINE + \
711	"\n" + FEET_1 % title + \
712	"\n" + FEET_2
713
714	return html_string, text_string, title
715
716	def _check_html_format(self, name):
717	"""
718	Check string '%' for html format
719	"""
720	if name.count('%'):
721	name = name.replace('%', '&#37')
722
723	return name
724
725	def report(self, figs=None, canvases=None):
726	"""
727	Invoke report dialog panel
728
729	: param figs: list of pylab figures [list]
730	"""
731	from sas.sasgui.perspectives.fitting.report_dialog import ReportDialog
732	# get the strings for report
733	html_str, text_str, title = self.set_report_string()
734	# Allow 2 figures to append
735	if len(figs) == 1:
736	add_str = FEET_3
737	elif len(figs) == 2:
738	add_str = ELINE
739	add_str += FEET_2 % ("%s")
740	add_str += ELINE
741	add_str += FEET_3
742	elif len(figs) > 2:
743	add_str = ELINE
744	add_str += FEET_2 % ("%s")
745	add_str += ELINE
746	add_str += FEET_2 % ("%s")
747	add_str += ELINE
748	add_str += FEET_3
749	else:
750	add_str = ""
751
752	# final report html strings
753	report_str = html_str % ("%s") + add_str
754
755	# make plot image
756	images = self.set_plot_state(figs, canvases)
757	report_list = [report_str, text_str, images]
758	dialog = ReportDialog(report_list, None, wx.ID_ANY, "")
759	dialog.Show()
760
761	def _to_xml_helper(self, thelist, element, newdoc):
762	"""
763	Helper method to create xml file for saving state
764	"""
765	for item in thelist:
766	sub_element = newdoc.createElement('parameter')
767	sub_element.setAttribute('name', str(item[1]))
768	sub_element.setAttribute('value', str(item[2]))
769	sub_element.setAttribute('selected_to_fit', str(item[0]))
770	sub_element.setAttribute('error_displayed', str(item[4][0]))
771	sub_element.setAttribute('error_value', str(item[4][1]))
772	sub_element.setAttribute('minimum_displayed', str(item[5][0]))
773	sub_element.setAttribute('minimum_value', str(item[5][1]))
774	sub_element.setAttribute('maximum_displayed', str(item[6][0]))
775	sub_element.setAttribute('maximum_value', str(item[6][1]))
776	sub_element.setAttribute('unit', str(item[7]))
777	element.appendChild(sub_element)
778
779	def to_xml(self, file="fitting_state.fitv", doc=None,
780	entry_node=None, batch_fit_state=None):
781	"""
782	Writes the state of the fit panel to file, as XML.
783
784	Compatible with standalone writing, or appending to an
785	already existing XML document. In that case, the XML document is
786	required. An optional entry node in the XML document may also be given.
787
788	:param file: file to write to
789	:param doc: XML document object [optional]
790	:param entry_node: XML node within the XML document at which we
791	will append the data [optional]
792	:param batch_fit_state: simultaneous fit state
793	"""
794	from xml.dom.minidom import getDOMImplementation
795
796	# Check whether we have to write a standalone XML file
797	if doc is None:
798	impl = getDOMImplementation()
799	doc_type = impl.createDocumentType(FITTING_NODE_NAME, "1.0", "1.0")
800	newdoc = impl.createDocument(None, FITTING_NODE_NAME, doc_type)
801	top_element = newdoc.documentElement
802	else:
803	# We are appending to an existing document
804	newdoc = doc
805	try:
806	top_element = newdoc.createElement(FITTING_NODE_NAME)
807	except:
808	string = etree.tostring(doc, pretty_print=True)
809	newdoc = parseString(string)
810	top_element = newdoc.createElement(FITTING_NODE_NAME)
811	if entry_node is None:
812	newdoc.documentElement.appendChild(top_element)
813	else:
814	try:
815	entry_node.appendChild(top_element)
816	except:
817	node_name = entry_node.tag
818	node_list = newdoc.getElementsByTagName(node_name)
819	entry_node = node_list.item(0)
820	entry_node.appendChild(top_element)
821
822	attr = newdoc.createAttribute("version")
823	from sas import sasview
824	attr.nodeValue = sasview.__version__
825	# attr.nodeValue = '1.0'
826	top_element.setAttributeNode(attr)
827
828	# File name
829	element = newdoc.createElement("filename")
830	if self.file is not None:
831	element.appendChild(newdoc.createTextNode(str(self.file)))
832	else:
833	element.appendChild(newdoc.createTextNode(str(file)))
834	top_element.appendChild(element)
835
836	element = newdoc.createElement("timestamp")
837	element.appendChild(newdoc.createTextNode(time.ctime(self.timestamp)))
838	attr = newdoc.createAttribute("epoch")
839	attr.nodeValue = str(self.timestamp)
840	element.setAttributeNode(attr)
841	top_element.appendChild(element)
842
843	# Inputs
844	inputs = newdoc.createElement("Attributes")
845	top_element.appendChild(inputs)
846
847	if self.data is not None and hasattr(self.data, "group_id"):
848	self.data_group_id = self.data.group_id
849	if self.data is not None and hasattr(self.data, "is_data"):
850	self.is_data = self.data.is_data
851	if self.data is not None:
852	self.data_name = self.data.name
853	if self.data is not None and hasattr(self.data, "id"):
854	self.data_id = self.data.id
855
856	for item in LIST_OF_DATA_ATTRIBUTES:
857	element = newdoc.createElement(item[0])
858	element.setAttribute(item[0], str(getattr(self, item[1])))
859	inputs.appendChild(element)
860
861	for item in LIST_OF_STATE_ATTRIBUTES:
862	element = newdoc.createElement(item[0])
863	element.setAttribute(item[0], str(getattr(self, item[1])))
864	inputs.appendChild(element)
865
866	# For self.values ={ disp_param_name: [vals,...],...}
867	# and for self.weights ={ disp_param_name: [weights,...],...}
868	for item in LIST_OF_MODEL_ATTRIBUTES:
869	element = newdoc.createElement(item[0])
870	value_list = getattr(self, item[1])
871	for key, value in value_list.iteritems():
872	sub_element = newdoc.createElement(key)
873	sub_element.setAttribute('name', str(key))
874	for val in value:
875	sub_element.appendChild(newdoc.createTextNode(str(val)))
876
877	element.appendChild(sub_element)
878	inputs.appendChild(element)
879
880	# Create doc for the dictionary of self._disp_obj_dic
881	for tagname, varname, tagtype in DISPERSION_LIST:
882	element = newdoc.createElement(tagname)
883	value_list = getattr(self, varname)
884	for key, value in value_list.iteritems():
885	sub_element = newdoc.createElement(key)
886	sub_element.setAttribute('name', str(key))
887	sub_element.setAttribute('value', str(value))
888	element.appendChild(sub_element)
889	inputs.appendChild(element)
890
891	for item in LIST_OF_STATE_PARAMETERS:
892	element = newdoc.createElement(item[0])
893	self._to_xml_helper(thelist=getattr(self, item[1]),
894	element=element, newdoc=newdoc)
895	inputs.appendChild(element)
896
897	# Combined and Simultaneous Fit Parameters
898	if batch_fit_state is not None:
899	batch_combo = newdoc.createElement('simultaneous_fit')
900	top_element.appendChild(batch_combo)
901
902	# Simultaneous Fit Number For Linking Later
903	element = newdoc.createElement('sim_fit_number')
904	element.setAttribute('fit_number', str(batch_fit_state.fit_page_no))
905	batch_combo.appendChild(element)
906
907	# Save constraints
908	constraints = newdoc.createElement('constraints')
909	batch_combo.appendChild(constraints)
910	for constraint in batch_fit_state.constraints_list:
911	if constraint.model_cbox.GetValue() != "":
912	# model_cbox, param_cbox, egal_txt, constraint,
913	# btRemove, sizer
914	doc_cons = newdoc.createElement('constraint')
915	doc_cons.setAttribute('model_cbox',
916	str(constraint.model_cbox.GetValue()))
917	doc_cons.setAttribute('param_cbox',
918	str(constraint.param_cbox.GetValue()))
919	doc_cons.setAttribute('egal_txt',
920	str(constraint.egal_txt.GetLabel()))
921	doc_cons.setAttribute('constraint',
922	str(constraint.constraint.GetValue()))
923	constraints.appendChild(doc_cons)
924
925	# Save all models
926	models = newdoc.createElement('model_list')
927	batch_combo.appendChild(models)
928	for model in batch_fit_state.model_list:
929	doc_model = newdoc.createElement('model_list_item')
930	doc_model.setAttribute('checked', str(model[0].GetValue()))
931	keys = model[1].keys()
932	doc_model.setAttribute('name', str(keys[0]))
933	values = model[1].get(keys[0])
934	doc_model.setAttribute('fit_number', str(model[2]))
935	doc_model.setAttribute('fit_page_source', str(model[3]))
936	doc_model.setAttribute('model_name', str(values.model.id))
937	models.appendChild(doc_model)
938
939	# Select All Checkbox
940	element = newdoc.createElement('select_all')
941	if batch_fit_state.select_all:
942	element.setAttribute('checked', 'True')
943	else:
944	element.setAttribute('checked', 'False')
945	batch_combo.appendChild(element)
946
947	# Save the file
948	if doc is None:
949	fd = open(file, 'w')
950	fd.write(newdoc.toprettyxml())
951	fd.close()
952	return None
953	else:
954	return newdoc
955
956	def _from_xml_helper(self, node, list):
957	"""
958	Helper function to write state to xml
959	"""
960	for item in node:
961	try:
962	name = item.get('name')
963	except:
964	name = None
965	try:
966	value = item.get('value')
967	except:
968	value = None
969	try:
970	selected_to_fit = (item.get('selected_to_fit') == "True")
971	except:
972	selected_to_fit = None
973	try:
974	error_displayed = (item.get('error_displayed') == "True")
975	except:
976	error_displayed = None
977	try:
978	error_value = item.get('error_value')
979	except:
980	error_value = None
981	try:
982	minimum_displayed = (item.get('minimum_displayed') == "True")
983	except:
984	minimum_displayed = None
985	try:
986	minimum_value = item.get('minimum_value')
987	except:
988	minimum_value = None
989	try:
990	maximum_displayed = (item.get('maximum_displayed') == "True")
991	except:
992	maximum_displayed = None
993	try:
994	maximum_value = item.get('maximum_value')
995	except:
996	maximum_value = None
997	try:
998	unit = item.get('unit')
999	except:
1000	unit = None
1001	list.append([selected_to_fit, name, value, "+/-",
1002	[error_displayed, error_value],
1003	[minimum_displayed, minimum_value],
1004	[maximum_displayed, maximum_value], unit])
1005
1006	def from_xml(self, file=None, node=None):
1007	"""
1008	Load fitting state from a file
1009
1010	:param file: .fitv file
1011	:param node: node of a XML document to read from
1012	"""
1013	if file is not None:
1014	msg = "PageState no longer supports non-CanSAS"
1015	msg += " format for fitting files"
1016	raise RuntimeError, msg
1017
1018	if node.get('version'):
1019	# Get the version for model conversion purposes
1020	self.version = tuple(int(e) for e in
1021	str.split(node.get('version'), "."))
1022	# The tuple must be at least 3 items long
1023	while len(self.version) < 3:
1024	ver_list = list(self.version)
1025	ver_list.append(0)
1026	self.version = tuple(ver_list)
1027
1028	# Get file name
1029	entry = get_content('ns:filename', node)
1030	if entry is not None:
1031	self.file = entry.text.strip()
1032
1033	# Get time stamp
1034	entry = get_content('ns:timestamp', node)
1035	if entry is not None and entry.get('epoch'):
1036	try:
1037	self.timestamp = float(entry.get('epoch'))
1038	except:
1039	msg = "PageState.fromXML: Could not"
1040	msg += " read timestamp\n %s" % sys.exc_value
1041	logger.error(msg)
1042
1043	if entry is not None:
1044	# Parse fitting attributes
1045	entry = get_content('ns:Attributes', node)
1046	for item in LIST_OF_DATA_ATTRIBUTES:
1047	node = get_content('ns:%s' % item[0], entry)
1048	setattr(self, item[0], parse_entry_helper(node, item))
1049
1050	dx_old_node = get_content('ns:%s' % 'dx_min', entry)
1051	for item in LIST_OF_STATE_ATTRIBUTES:
1052	if item[0] == "dx_percent" and dx_old_node is not None:
1053	dxmin = ["dx_min", "dx_min", "float"]
1054	setattr(self, item[0], parse_entry_helper(dx_old_node,
1055	dxmin))
1056	self.dx_old = True
1057	else:
1058	node = get_content('ns:%s' % item[0], entry)
1059	setattr(self, item[0], parse_entry_helper(node, item))
1060
1061	for item in LIST_OF_STATE_PARAMETERS:
1062	node = get_content("ns:%s" % item[0], entry)
1063	self._from_xml_helper(node=node,
1064	list=getattr(self, item[1]))
1065
1066	# Recover _disp_obj_dict from xml file
1067	self._disp_obj_dict = {}
1068	for tagname, varname, tagtype in DISPERSION_LIST:
1069	node = get_content("ns:%s" % tagname, entry)
1070	for attr in node:
1071	parameter = str(attr.get('name'))
1072	value = attr.get('value')
1073	if value.startswith("<"):
1074	try:
1075	# <path.to.NamedDistribution object/instance...>
1076	cls_name = value[1:].split()[0].split('.')[-1]
1077	cls = getattr(sasmodels.weights, cls_name)
1078	value = cls.type
1079	except Exception:
1080	base = "unable to load distribution %r for %s"
1081	logger.error(base % (value, parameter))
1082	continue
1083	_disp_obj_dict = getattr(self, varname)
1084	_disp_obj_dict[parameter] = value
1085
1086	# get self.values and self.weights dic. if exists
1087	for tagname, varname in LIST_OF_MODEL_ATTRIBUTES:
1088	node = get_content("ns:%s" % tagname, entry)
1089	dic = {}
1090	value_list = []
1091	for par in node:
1092	name = par.get('name')
1093	values = par.text.split()
1094	# Get lines only with numbers
1095	for line in values:
1096	try:
1097	val = float(line)
1098	value_list.append(val)
1099	except Exception:
1100	# pass if line is empty (it happens)
1101	msg = ("Error reading %r from %s %s\n"
1102	% (line, tagname, name))
1103	logger.error(msg + traceback.format_exc())
1104	dic[name] = np.array(value_list)
1105	setattr(self, varname, dic)
1106
1107	def set_plot_state(self, figs, canvases):
1108	"""
1109	Build image state that wx.html understand
1110	by plotting, putting it into wx.FileSystem image object
1111
1112	"""
1113	images = []
1114
1115	# Reset memory
1116	self.imgRAM = None
1117	wx.MemoryFSHandler()
1118
1119	# For no figures in the list, prepare empty plot
1120	if figs is None or len(figs) == 0:
1121	figs = [None]
1122
1123	# Loop over the list of figures
1124	# use wx.MemoryFSHandler
1125	self.imgRAM = wx.MemoryFSHandler()
1126	for fig in figs:
1127	if fig is not None:
1128	ind = figs.index(fig)
1129	canvas = canvases[ind]
1130
1131	# store the image in wx.FileSystem Object
1132	wx.FileSystem.AddHandler(wx.MemoryFSHandler())
1133
1134	# index of the fig
1135	ind = figs.index(fig)
1136
1137	# AddFile, image can be retrieved with 'memory:filename'
1138	self.imgRAM.AddFile('img_fit%s.png' % ind,
1139	canvas.bitmap, wx.BITMAP_TYPE_PNG)
1140
1141	# append figs
1142	images.append(fig)
1143
1144	return images
1145
1146
1147	class Reader(CansasReader):
1148	"""
1149	Class to load a .fitv fitting file
1150	"""
1151	# File type
1152	type_name = "Fitting"
1153
1154	# Wildcards
1155	type = ["Fitting files (.fitv)\|.fitv"
1156	"SASView file (.svs)\|.svs"]
1157	# List of allowed extensions
1158	ext = ['.fitv', '.FITV', '.svs', 'SVS']
1159
1160	def __init__(self, call_back=None, cansas=True):
1161	CansasReader.__init__(self)
1162	"""
1163	Initialize the call-back method to be called
1164	after we load a file
1165
1166	:param call_back: call-back method
1167	:param cansas: True = files will be written/read in CanSAS format
1168	False = write CanSAS format
1169
1170	"""
1171	# Call back method to be executed after a file is read
1172	self.call_back = call_back
1173	# CanSAS format flag
1174	self.cansas = cansas
1175	self.state = None
1176	# batch fitting params for saving
1177	self.batchfit_params = []
1178
1179	def get_state(self):
1180	return self.state
1181
1182	def read(self, path):
1183	"""
1184	Load a new P(r) inversion state from file
1185
1186	:param path: file path
1187
1188	"""
1189	if self.cansas:
1190	return self._read_cansas(path)
1191
1192	def _parse_state(self, entry):
1193	"""
1194	Read a fit result from an XML node
1195
1196	:param entry: XML node to read from
1197	:return: PageState object
1198	"""
1199	# Create an empty state
1200	state = None
1201	# Locate the P(r) node
1202	try:
1203	nodes = entry.xpath('ns:%s' % FITTING_NODE_NAME,
1204	namespaces={'ns': CANSAS_NS})
1205	if nodes:
1206	# Create an empty state
1207	state = PageState()
1208	state.from_xml(node=nodes[0])
1209
1210	except:
1211	logger.info("XML document does not contain fitting information.\n"
1212	+ traceback.format_exc())
1213
1214	return state
1215
1216	def _parse_simfit_state(self, entry):
1217	"""
1218	Parses the saved data for a simultaneous fit
1219	:param entry: XML object to read from
1220	:return: XML object for a simultaneous fit or None
1221	"""
1222	nodes = entry.xpath('ns:%s' % FITTING_NODE_NAME,
1223	namespaces={'ns': CANSAS_NS})
1224	if nodes:
1225	simfitstate = nodes[0].xpath('ns:simultaneous_fit',
1226	namespaces={'ns': CANSAS_NS})
1227	if simfitstate:
1228	from simfitpage import SimFitPageState
1229	sim_fit_state = SimFitPageState()
1230	simfitstate_0 = simfitstate[0]
1231	all = simfitstate_0.xpath('ns:select_all',
1232	namespaces={'ns': CANSAS_NS})
1233	atts = all[0].attrib
1234	checked = atts.get('checked')
1235	sim_fit_state.select_all = bool(checked)
1236	model_list = simfitstate_0.xpath('ns:model_list',
1237	namespaces={'ns': CANSAS_NS})
1238	model_list_items = model_list[0].xpath('ns:model_list_item',
1239	namespaces={'ns':
1240	CANSAS_NS})
1241	for model in model_list_items:
1242	attrs = model.attrib
1243	sim_fit_state.model_list.append(attrs)
1244
1245	constraints = simfitstate_0.xpath('ns:constraints',
1246	namespaces={'ns': CANSAS_NS})
1247	constraint_list = constraints[0].xpath('ns:constraint',
1248	namespaces={'ns': CANSAS_NS})
1249	for constraint in constraint_list:
1250	attrs = constraint.attrib
1251	sim_fit_state.constraints_list.append(attrs)
1252
1253	return sim_fit_state
1254	else:
1255	return None
1256
1257	def _parse_save_state_entry(self, dom):
1258	"""
1259	Parse a SASentry
1260
1261	:param node: SASentry node
1262
1263	:return: Data1D/Data2D object
1264
1265	"""
1266	node = dom.xpath('ns:data_class', namespaces={'ns': CANSAS_NS})
1267	return_value, _ = self._parse_entry(dom)
1268	return return_value, _
1269
1270	def _read_cansas(self, path):
1271	"""
1272	Load data and fitting information from a CanSAS XML file.
1273
1274	:param path: file path
1275	:return: Data1D object if a single SASentry was found,
1276	or a list of Data1D objects if multiple entries were found,
1277	or None of nothing was found
1278	:raise RuntimeError: when the file can't be opened
1279	:raise ValueError: when the length of the data vectors are inconsistent
1280	"""
1281	output = []
1282	simfitstate = None
1283	basename = os.path.basename(path)
1284	root, extension = os.path.splitext(basename)
1285	ext = extension.lower()
1286	try:
1287	if os.path.isfile(path):
1288	if ext in self.ext or ext == '.xml':
1289	tree = etree.parse(path, parser=etree.ETCompatXMLParser())
1290	# Check the format version number
1291	# Specifying the namespace will take care of the file
1292	# format version
1293	root = tree.getroot()
1294	entry_list = root.xpath('ns:SASentry',
1295	namespaces={'ns': CANSAS_NS})
1296	for entry in entry_list:
1297	try:
1298	sas_entry, _ = self._parse_save_state_entry(entry)
1299	except:
1300	raise
1301	fitstate = self._parse_state(entry)
1302
1303	# state could be None when .svs file is loaded
1304	# in this case, skip appending to output
1305	if fitstate is not None:
1306	sas_entry.meta_data['fitstate'] = fitstate
1307	sas_entry.filename = fitstate.file
1308	output.append(sas_entry)
1309
1310	else:
1311	self.call_back(format=ext)
1312	raise RuntimeError, "%s is not a file" % path
1313
1314	# Return output consistent with the loader's api
1315	if len(output) == 0:
1316	self.call_back(state=None, datainfo=None, format=ext)
1317	return None
1318	else:
1319	for ind in range(len(output)):
1320	# Call back to post the new state
1321	state = output[ind].meta_data['fitstate']
1322	t = time.localtime(state.timestamp)
1323	time_str = time.strftime("%b %d %H:%M", t)
1324	# Check that no time stamp is already appended
1325	max_char = state.file.find("[")
1326	if max_char < 0:
1327	max_char = len(state.file)
1328	original_fname = state.file[0:max_char]
1329	state.file = original_fname + ' [' + time_str + ']'
1330
1331	if state is not None and state.is_data is not None:
1332	output[ind].is_data = state.is_data
1333
1334	output[ind].filename = state.file
1335	state.data = output[ind]
1336	state.data.name = output[ind].filename # state.data_name
1337	state.data.id = state.data_id
1338	if state.is_data is not None:
1339	state.data.is_data = state.is_data
1340	if output[ind].run_name is not None\
1341	and len(output[ind].run_name) != 0:
1342	if isinstance(output[ind].run_name, dict):
1343	name = output[ind].run_name.keys()[0]
1344	else:
1345	name = output[ind].run_name
1346	else:
1347	name = original_fname
1348	state.data.group_id = name
1349	state.version = fitstate.version
1350	# store state in fitting
1351	self.call_back(state=state,
1352	datainfo=output[ind], format=ext)
1353	self.state = state
1354	simfitstate = self._parse_simfit_state(entry)
1355	if simfitstate is not None:
1356	self.call_back(state=simfitstate)
1357
1358	return output
1359	except:
1360	self.call_back(format=ext)
1361	raise
1362
1363	def write(self, filename, datainfo=None, fitstate=None):
1364	"""
1365	Write the content of a Data1D as a CanSAS XML file only for standalone
1366
1367	:param filename: name of the file to write
1368	:param datainfo: Data1D object
1369	:param fitstate: PageState object
1370
1371	"""
1372	# Sanity check
1373	if self.cansas:
1374	# Add fitting information to the XML document
1375	doc = self.write_toXML(datainfo, fitstate)
1376	# Write the XML document
1377	else:
1378	doc = fitstate.to_xml(file=filename)
1379
1380	# Save the document no matter the type
1381	fd = open(filename, 'w')
1382	fd.write(doc.toprettyxml())
1383	fd.close()
1384
1385	def write_toXML(self, datainfo=None, state=None, batchfit=None):
1386	"""
1387	Write toXML, a helper for write(),
1388	could be used by guimanager._on_save()
1389
1390	: return: xml doc
1391	"""
1392
1393	self.batchfit_params = batchfit
1394	if state.data is None or not state.data.is_data:
1395	return None
1396	# make sure title and data run are filled.
1397	if state.data.title is None or state.data.title == '':
1398	state.data.title = state.data.name
1399	if state.data.run_name is None or state.data.run_name == {}:
1400	state.data.run = [str(state.data.name)]
1401	state.data.run_name[0] = state.data.name
1402
1403	data = state.data
1404	doc, sasentry = self._to_xml_doc(data)
1405
1406	if state is not None:
1407	doc = state.to_xml(doc=doc, file=data.filename, entry_node=sasentry,
1408	batch_fit_state=self.batchfit_params)
1409
1410	return doc
1411
1412	# Simple html report templet
1413	HEADER = "<html>\n"
1414	HEADER += "<head>\n"
1415	HEADER += "<meta http-equiv=Content-Type content='text/html; "
1416	HEADER += "charset=windows-1252'> \n"
1417	HEADER += "<meta name=Generator >\n"
1418	HEADER += "</head>\n"
1419	HEADER += "<body lang=EN-US>\n"
1420	HEADER += "<div class=WordSection1>\n"
1421	HEADER += "<p class=MsoNormal><b><span ><center><font size='4' >"
1422	HEADER += "%s</font></center></span></center></b></p>"
1423	HEADER += "<p class=MsoNormal> </p>"
1424	PARA = "<p class=MsoNormal><font size='4' > %s \n"
1425	PARA += "</font></p>"
1426	CENTRE = "<p class=MsoNormal><center><font size='4' > %s \n"
1427	CENTRE += "</font></center></p>"
1428	FEET_1 = \
1429	"""
1430	<p class=MsoNormal> </p>
1431	<br>
1432	<p class=MsoNormal><b><span ><center> <font size='4' > Graph
1433	</font></span></center></b></p>
1434	<p class=MsoNormal> </p>
1435	<center>
1436	<br><font size='4' >Model Computation</font>
1437	<br><font size='4' >Data: "%s"</font><br>
1438	"""
1439	FEET_2 = \
1440	"""
1441	<img src="%s" >
1442	</img>
1443	"""
1444	FEET_3 = \
1445	"""
1446	</center>
1447	</div>
1448	</body>
1449	</html>
1450	"""
1451	ELINE = "<p class=MsoNormal> </p>"

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

SasView

source: sasview/src/sas/sasgui/perspectives/fitting/pagestate.py @ 3152a02

Download in other formats: