def test2(self): "read events and check them" from mccomponents.detector.reduction_utils import readevents evts = readevents(eventfilename) self.assertEqual(len(evts), 1) evt = evts[0] longpixelID, tofchannel, probability = evt self.assertEqual(longpixelID, pixelID + npixels * (tubeID + ntubes * (packID))) self.assertEqual(tofchannel, int((neutronT - tofmin) / tofstep)) self.assertEqual(probability, prob) return
def main(): from sim_params import instrument, tofparams, eventsdat, Idpt_filename as filename from mccomponents.detector.reduction_utils import readevents events = readevents(eventsdat) Idpt = events2Idpt(events, instrument, tofparams) #save to file import os if os.path.exists(filename): os.remove(filename) import histogram.hdf as hh hh.dump(Idpt, filename, '/', 'c') return
def main(): from sim_params import instrument, tofparams, eventsdat, Idpt_filename as filename from mccomponents.detector.reduction_utils import readevents events = readevents( eventsdat ) Idpt = events2Idpt( events, instrument, tofparams ) #save to file import os if os.path.exists(filename): os.remove( filename ) import histogram.hdf as hh hh.dump( Idpt, filename, '/', 'c' ) return
def test1a(self): from mccomponents.detector.reduction_utils import readevents events = readevents( outfilename ) n = len(events) print "number of cases where absorption happen: ", n self.assert_( abs(n-(nevents*absorption_weight)) < 3*N.sqrt(n) ) p = sum([ e[2] for e in events ] ) print "absorbed neutrons: ", p self.assert_( abs( p-(nevents*0.91) ) < 3*N.sqrt(p) ) t = L/vi tchannel = int( (t-tmin)/tstep ) for e in events: # should hit center of tube self.assert_( abs(e[0] - (npixelsperdet*((packindexat0+0.5)*ndetsperpack)) ) <= 1 ) # tof channel self.assert_( abs(e[1] - tchannel) <= 1 ) continue return
def test1a(self): "simple. one detector -- verify" from mccomponents.detector.reduction_utils import readevents events = readevents( 'events.dat' ) n = len(events) self.assert_( abs(n-(nevents*absorption_weight)) < 3*N.sqrt(n) ) p = sum([ e[2] for e in events ] ) # print "%s should be almost equalt to %s" % (p, nevents*absorption_probability) self.assert_( abs( p-(nevents*absorption_probability) ) < 3*N.sqrt(p) ) # self.assert_( abs( p-(nevents*0.908484) ) < 3*N.sqrt(p) ) t = L/vi tchannel = int( (t-tmin)/tstep ) for e in events: # should hit center of tube self.assert_( abs(e[0] - (npixels/2)) <= 1 ) self.assertEqual( e[1], tchannel ) # continue return
def test1a(self): from mccomponents.detector.reduction_utils import readevents events = readevents(outfilename) n = len(events) print "number of cases where absorption happen: ", n self.assert_(abs(n - (nevents * absorption_weight)) < 3 * N.sqrt(n)) p = sum([e[2] for e in events]) print "absorbed neutrons: ", p self.assert_(abs(p - (nevents * 0.91)) < 3 * N.sqrt(p)) t = L / vi tchannel = int((t - tmin) / tstep) for e in events: # should hit center of tube self.assert_( abs(e[0] - (npixelsperdet * ((packindexat0 + 0.5) * ndetsperpack))) <= 1) # tof channel self.assert_(abs(e[1] - tchannel) <= 1) continue return