def test1(self): "multiphonon.ms" mpsqe = hh.load(os.path.join(datadir, "V-S2..5.h5")) mssqe = ms.sqe(mpsqe, Ei=110.) if interactive: H.plot(mssqe, min=0) return
def test(): from histogram import histogram, axis, arange, plot xaxis = axis('x', arange(5), unit='meter') yaxis = axis('y', arange(7), unit='cm') axes = [xaxis, yaxis] h = histogram("intensity", axes, fromfunction=lambda x, y: x**2 + y**2) print(h) plot(h) help(h) slice = h[3, ()]
def test1b(self): iqehist = hh.load(os.path.join(datadir, "V-iqe.h5")) from multiphonon.sqe import interp newiqe = interp(iqehist, newE=np.arange(-50, 50, 1.0)) DOS = sqe2dos.singlephonon_sqe2dos(newiqe, T=300, Ecutoff=65.0, elastic_E_cutoff=(-20.0, 6.7), M=50.94) # plot if interactive: H.plot(DOS) return
def test(): from histogram import histogram, axis, arange, plot xaxis = axis('x', arange(5), unit='meter') yaxis = axis('y', arange(7), unit='cm') axes = [xaxis, yaxis] h = histogram( "intensity", axes, fromfunction=lambda x,y: x**2+y**2) print h plot(h) help(h) slice = h[3, ()]
def plotHistogram(roller): cases = list(itertools.islice(newStats(roller), 1000)) def toBuckets(x): "groups into 3 buckets (x <= 8) (8 < x < 13) (x >= 13)" if x <= 8: return 1 if x >= 13: return 3 return 2 histogram.plot(cases) histogram.plot([toBuckets(x) for x in cases])
def test2(self): "multiphonon.forward.phonon.sqe: UN N dos" E, g= readdos() dE = E[1] - E[0] # from multiphonon.forward.phonon import sqe q,e,i = sqe(E, g, T=5, M=14, N=7, Qmax=45.) from histogram import plot, histogram axes = [('Q', q, 'angstrom**-1'), ('E', e, 'meV')] iqe = histogram('iqe', axes, i) if interactive: plot(iqe) return
def test2(self): "multiphonon.sqe: UN N dos" dos = readdos() E = dos.energy dE = E[1] - E[0] g = dos.I # from mccomponents.sample.phonon.multiphonon import sqe q, e, i = sqe(E, g, T=5, M=14, N=7, Qmax=45.) from histogram import plot, histogram axes = [('Q', q, 'angstrom**-1'), ('E', e, 'meV')] iqe = histogram('iqe', axes, i) plot(iqe) return
def test2(self): "multiphonon.sqe: UN N dos" dos = readdos() E = dos.energy dE = E[1] - E[0] g = dos.I # from mccomponents.sample.phonon.multiphonon import sqe q,e,i = sqe(E, g, T=5, M=14, N=7, Qmax=45.) from histogram import plot, histogram axes = [('Q', q, 'angstrom**-1'), ('E', e, 'meV')] iqe = histogram('iqe', axes, i) plot(iqe) return
def test1b(self): "singlephonon_sqe2dos: exp vanadium SQE -> DOS" iqehist = hh.load(os.path.join(datadir, "V-iqe.h5")) from multiphonon.sqe import interp newiqe = interp(iqehist, newE=np.arange(-50, 50, 1.)) DOS = sqe2dos.singlephonon_sqe2dos(newiqe, T=300, Ecutoff=65., elastic_E_cutoff=(-20., 6.7), M=50.94) path = os.path.join(here, 'expected_results', 'test1b-dos.h5') expected = hh.load(path) self.assert_(np.allclose(DOS.I, expected.I)) self.assert_(np.allclose(DOS.E2, expected.E2)) # plot if interactive: H.plot(DOS) return