# $Id$
import DST
from SOM import SOM
from SOM import SO
filename_SOM1 = "stuff1.dat"
SOM1 = SOM()
SOM1.attr_list["filename"] = filename_SOM1
SOM1.setTitle("This is a test")
for i in range(2):
SO1 = SO()
for j in range(10):
SO1.id = i
SO1.axis[0].val.append(j+1)
SO1.y.append(1000+j+(20*j))
SO1.var_y.append(100+j)
SO1.axis[0].val.append(11)
SOM1.append(SO1)
file = open(filename_SOM1, "w")
gsas = DST.GsasDST(file)
gsas.writeSOM(SOM1)
gsas.release_resource()
from time import localtime, strftime, time
filename_SOM1 = "stuff1.dat"
SOM1 = SOM()
SOM1.attr_list["filename"] = filename_SOM1
SOM1.attr_list["epoch"] = time()
SOM1.attr_list["timestamp"] = DST.make_ISO8601(SOM1.attr_list["epoch"])
SOM1.attr_list["username"] = "******"
SOM1.setAllAxisLabels(["Q", "E"])
SOM1.setAllAxisUnits(["A-1", "meV"])
SOM1.setYLabel("Intensity")
SOM1.setYUnits("Counts/(meV A-1))")
SO1 = SO(2)
SO1.id = 0
SO1.axis[0].val.extend(range(5))
SO1.axis[1].val.extend(range(10))
y_len = (len(SO1.axis[0].val) - 1) * (len(SO1.axis[1].val) - 1)
y = range(y_len)
SO1.y.extend(y)
SO1.var_y.extend(y)
SOM1.append(SO1)
file = open(filename_SOM1, "w")
d2d = DST.Dave2dDST(file)
d2d.writeSOM(SOM1)
d2d.release_resource()
SOM2 = SOM()
SOM2.copyAttributes(SOM1)
SOM2.setAllAxisLabels(["Wavelength"])
SOM2.setAllAxisUnits(["Angstroms"])
SOM2.setYLabel("Counts")
SOM2.setYUnits("counts")
SOM2.setDataSetType("histogram")
else:
SOM2 = None
for i in range(2):
SO1 = SO(construct=True, withVar=options.withXvar)
if options.extraSom:
SO2 = SO(construct=True)
for j in range(10):
SO1.id = i
SO1.axis[0].val.append(j+1)
if options.withXvar:
SO1.axis[0].var.append(j+1)
SO1.y.append(1000+j+(20*j))
SO1.var_y.append(100+j)
if options.extraSom:
SO2.id = i
SO2.axis[0].val.append((j+1)*1.0579838)
SO2.y.append(1000+j+(20*j))
SO2.var_y.append(100+j)
SO1.axis[0].val.append(11)
if options.withXvar:
SO1.axis[0].var.append(11)
SOM1.setAllAxisLabels(["momentum transfer", "energy transfer"])
SOM1.setAllAxisUnits(["1/A", "meV"])
SOM1.attr_list["data-title"] = "Test S(Q,E)"
SOM1.attr_list["data-run_number"] = "1344"
DSample = Sample()
DSample.name = "Test Sample"
DSample.nature = "CoCo"
SOM1.attr_list.sample = DSample
x = NessiList()
y = NessiList()
z = NessiList()
x.extend(0,1,2,3)
y.extend(0,1,2)
z.extend(1,2,3,4,5,6)
SO1 = SO(2)
SO1.id = ("bank1", (4, 32))
SO1.y = z
SO1.var_y = z
SO1.axis[0].val = x
SO1.axis[1].val = y
SOM1.append(SO1)
rednxs = DST.RedNxsDST("test_red.nxs")
rednxs.writeSOM(SOM1, entry_name="sqe")
rednxs.release_resource()
SOM1.setAllAxisLabels(["momentum transfer", "energy transfer"])
SOM1.setAllAxisUnits(["1/A", "meV"])
SOM1.attr_list["data-title"] = "Test S(Q,E)"
SOM1.attr_list["data-run_number"] = "1344"
DSample = Sample()
DSample.name = "Test Sample"
DSample.nature = "CoCo"
SOM1.attr_list.sample = DSample
x = NessiList()
y = NessiList()
z = NessiList()
x.extend(0, 1, 2, 3)
y.extend(0, 1, 2)
z.extend(1, 2, 3, 4, 5, 6)
SO1 = SO(2)
SO1.id = ("bank1", (4, 32))
SO1.y = z
SO1.var_y = z
SO1.axis[0].val = x
SO1.axis[1].val = y
SOM1.append(SO1)
rednxs = DST.RedNxsDST("test_red.nxs")
rednxs.writeSOM(SOM1, entry_name="sqe")
rednxs.release_resource()
#
# $Id$
import DST
from SOM import SOM
from SOM import SO
filename_SOM1 = "stuff1.dat"
SOM1 = SOM()
SOM1.attr_list["filename"] = filename_SOM1
SOM1.setTitle("This is a test")
for i in range(2):
SO1 = SO()
for j in range(10):
SO1.id = i
SO1.axis[0].val.append(j + 1)
SO1.y.append(1000 + j + (20 * j))
SO1.var_y.append(100 + j)
SO1.axis[0].val.append(11)
SOM1.append(SO1)
file = open(filename_SOM1, "w")
gsas = DST.GsasDST(file)
gsas.writeSOM(SOM1)
gsas.release_resource()
from time import localtime, strftime, time
filename_SOM1 = "stuff1.dat"
SOM1 = SOM()
SOM1.attr_list["filename"] = filename_SOM1
SOM1.attr_list["epoch"] = time()
SOM1.attr_list["timestamp"] = DST.make_ISO8601(SOM1.attr_list["epoch"])
SOM1.attr_list["username"] = "******"
SOM1.setAllAxisLabels(["Q", "E"])
SOM1.setAllAxisUnits(["A-1", "meV"])
SOM1.setYLabel("Intensity")
SOM1.setYUnits("Counts/(meV A-1))")
SO1 = SO(2)
SO1.id = 0
SO1.axis[0].val.extend(range(5))
SO1.axis[1].val.extend(range(10))
y_len = (len(SO1.axis[0].val)-1) * (len(SO1.axis[1].val)-1)
y = range(y_len)
SO1.y.extend(y)
SO1.var_y.extend(y)
SOM1.append(SO1)
file = open(filename_SOM1, "w")
d2d = DST.Dave2dDST(file)
d2d.writeSOM(SOM1)
d2d.release_resource()