def run(eventfile, tofbinsize, nxsfile, type, instrument): print(eventfile, nxsfile) from mccomponents.detector.event_utils import readEvents events = readEvents(eventfile) prefix = 'mcvine.instruments.%s.nxs' % instrument.upper() mod = '%s.%s' % (prefix, type) mod = __import__(mod, fromlist=['']) mod.write(events, tofbinsize, nxsfile) return
def run(eventfile, nxsfile, tofbinsize, type): print (eventfile, nxsfile) from mccomponents.detector.event_utils import readEvents events = readEvents(eventfile) prefix = 'mcvine.instruments.SEQUOIA.nxs' mod = '%s.%s' % (prefix, type) mod = __import__(mod, fromlist = ['']) mod.write(events, tofbinsize, nxsfile) return
def test2(self): "ARCS detector system. A neutron missing all pixels." # cmd = './sd --source.velocity=1000,0,-2000 --ncount=1 --output-dir=out2' cmd = './sd --source.position=-0.00875,0.00462,0.005 --source.velocity=5713.19,-765.203,9068.39 --ncount=1 --output-dir=out2' import os if os.system(cmd): raise RuntimeError, "%s failed" % cmd eventdata = 'out2/events.dat' from mccomponents.detector.event_utils import readEvents events = readEvents(eventdata) print events # self.assertEqual(len(events), 0) return
def run(eventfile, nxsfile, tofbinsize=0.1, type="processed", Ei = None): """tofbinsize: in microsecond type: processed or raw. processed is obsolete but kept for backward compatibility Ei: nominal incident energy in meV """ print (eventfile, nxsfile) from mccomponents.detector.event_utils import readEvents events = readEvents(eventfile) prefix = 'mcvine.instruments.ARCS.nxs' mod = '%s.%s' % (prefix, type) mod = __import__(mod, fromlist = ['']) mod.write(events, tofbinsize, nxsfile, Ei=Ei) return
def run(eventfile, nxsfile, tofbinsize=0.1, type="processed", Ei = None): """tofbinsize: in microsecond type: processed or raw. processed is obsolete but kept for backward compatibility Ei: nominal incident energy in meV """ print((eventfile, nxsfile)) from mccomponents.detector.event_utils import readEvents events = readEvents(eventfile) prefix = 'mcvine.instruments.ARCS.nxs' mod = '%s.%s' % (prefix, type) mod = __import__(mod, fromlist = ['']) mod.write(events, tofbinsize, nxsfile, Ei=Ei) return
def test1a(self): "ARCS detector system. Neutrons shotting at one single pixel." cmd = './sd --mpirun.nodes=2 --output-dir=out1a' import os if os.system(cmd): raise RuntimeError, "%s failed" % cmd eventdata = 'out1a/events.dat' from mccomponents.detector.event_utils import readEvents events = readEvents(eventdata) # print events # print events['p'].sum() p = events['p'].sum() print p self.assert_(p > 0.9 and p < 1) return
def test1a(self): "ARCS detector system. Neutrons shotting at one single pixel." cmd = './sd --mpirun.nodes=2 --output-dir=out1a' import os if os.system(cmd): raise RuntimeError, "%s failed" % cmd eventdata = 'out1a/events.dat' from mccomponents.detector.event_utils import readEvents events = readEvents(eventdata) # print events # print events['p'].sum() p = events['p'].sum() print p self.assert_(p>0.9 and p<1) return
def run(eventfile, nxsfile, tofbinsize=0.1, type=None, Ei=None): """tofbinsize: in microsecond """ assert type == 'processed', "Unsupported type %r" % type print((eventfile, nxsfile)) # read events from mccomponents.detector.event_utils import readEvents events = readEvents(eventfile) # write nxs write(events, tofbinsize, nxsfile) # if Ei: from .nxs import setEnergyRequest setEnergyRequest(nxsfile, Ei) return
def test(self): "ARCS detector system. Neutrons shotting at one single pixel." cmd = "MCVINE_MPI_LAUNCHER=serial ./sd" import os if os.system(cmd): raise RuntimeError, "%s failed" % cmd eventdata = "out/events.dat" from mccomponents.detector.event_utils import readEvents events = readEvents(eventdata) # print events # print events['p'].sum() p = events["p"].sum() print p self.assert_(p > 0.9 and p < 1) return
def getNumberOfEvts(filename, pixelID): from mccomponents.detector import event_utils evts = event_utils.readEvents(filename) _ = lambda evt: evt['pixelID'] == pixelID return len(list(filter(_, evts)))
def getNumberOfEvts(filename, pixelID): from mccomponents.detector import event_utils evts = event_utils.readEvents(filename) _ = lambda evt: evt['pixelID'] == pixelID return len(filter(_, evts))
#!/usr/bin/env python from mccomponents.detector import event_utils as evu f = 'out/events.dat' print print "* reading events from %s" % f events = evu.readEvents(f) print print "* # of events", len(events) print print "* events (pixelID, tofID, weight):" print events
#!/usr/bin/env python from mccomponents.detector import event_utils as evu f = 'out/events.dat' print() print(("* reading events from %s" % f)) events = evu.readEvents(f) print() print(("* # of events", len(events))) print() print("* events (pixelID, tofID, weight):") print(events)