def elementFactory(self, **kwds):
datapath = kwds.get("idf-data-path")
hdfpath = kwds.get("histogram-hdf-path")
if datapath:
from mccomponents.sample.idf import readSQE
sqe = readSQE(datapath)
pass
elif hdfpath:
from histogram.hdf import load
try:
sqe = load(hdfpath)
except:
import os, traceback
f = os.path.dirname(hdfpath)
e = os.path.basename(hdfpath)
t = traceback.format_exc()
raise IOError, "unable to load histogram from hdf5 file %s, entry %s. Original traceback:\n%s" % (
f,
e,
t,
)
pass
else:
raise ValueError, "GridSQE needs path to " "idf data files or " "histogram hdf5 file "
auto_normalization = kwds.get("auto-normalization")
if auto_normalization:
auto_normalization = bool(auto_normalization)
norm = _calcNorm(sqe)
if abs(norm - 1) > 0.2:
if auto_normalization:
sqe.I /= norm
else:
raise RuntimeError, "S(Q,E) should average to ~1, got %s" % norm
from mccomponents.sample import gridsqe
return gridsqe(sqe)
def elementFactory(self, **kwds):
datapath = kwds.get('idf-data-path')
hdfpath = kwds.get('histogram-hdf-path')
if datapath:
from mccomponents.sample.idf import readSQE
sqe = readSQE(datapath)
pass
elif hdfpath:
from histogram.hdf import load
try:
sqe = load(hdfpath)
except:
import os, traceback
f = os.path.dirname(hdfpath)
e = os.path.basename(hdfpath)
t = traceback.format_exc()
raise IOError, "unable to load histogram from hdf5 file %s, entry %s. Original traceback:\n%s" % (
f, e, t)
pass
else:
raise ValueError, "GridSQE needs path to "\
"idf data files or "\
"histogram hdf5 file "
auto_normalization = kwds.get('auto-normalization')
if auto_normalization:
auto_normalization = bool(auto_normalization)
norm = _calcNorm(sqe)
if abs(norm - 1) > 0.2:
if auto_normalization:
sqe.I /= norm
else:
raise RuntimeError, "S(Q,E) should average to ~1, got %s" % norm
from mccomponents.sample import gridsqe
return gridsqe(sqe)
from mccomponents.sample.idf import readSQE
sqe = readSQE('SQE-examples')
import histogram.hdf as hh
hh.dump(sqe, 'sqe.h5', '/', 'c')
def test1(self):
from mccomponents.sample.idf import readSQE
sqe = readSQE(datapath)
import pickle
pickle.dump(sqe, open('sqehist.pkl', 'w'))
return
from mccomponents.sample.idf import readSQE sqe = readSQE( 'SQE-examples' ) import histogram.hdf as hh hh.dump( sqe, 'sqe.h5', '/', 'c' )
def test1(self):
from mccomponents.sample.idf import readSQE
sqe = readSQE( datapath )
import pickle
pickle.dump(sqe, open('sqehist.pkl','w') )
return
#!/usr/bin/env python import sys path = sys.argv[1] from mccomponents.sample.idf import readSQE sqe = readSQE( path ) from histogram.plotter import defaultPlotter defaultPlotter.plot( sqe )
