Beispiel #1
0
# $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()


Beispiel #2
0
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()
Beispiel #3
0
    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)
Beispiel #4
0
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()
Beispiel #5
0
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()
Beispiel #6
0
#

# $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()
Beispiel #7
0
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()