iCol1 = obsFile.nCol spt = obsFile.tickDuration masks = {} pickle.dump(offsets,pickleFile) for offset in offsets: print "begin offset=",offset masks[offset] = interval() for sec in range(beginTime,beginTime+expTime): tl = sec+obsFile.tickDuration*(offset-2) tr = sec+obsFile.tickDuration*(offset+3) masks[offset] = masks[offset] | interval([sec,tl]) | interval([tr,sec+1]) obsFile.cosmicMask = masks[offset] obsFile.switchOnCosmicTimeMask() nPhotonSum = 0 rows = np.zeros(0,dtype=np.int) cols = np.zeros(0,dtype=np.int) dts = np.zeros(0,dtype=np.double) secs = np.zeros(0,dtype=np.int) phs = np.zeros(0,dtype=np.double) for iRow in range(iRow0, iRow1): for iCol in range(iCol0, iCol1): #tpl = obsFile.getTimedPacketList(iRow, iCol, firstSec=beginTime, integrationTime=expTime) tpl = obsFile.getPackets(iRow, iCol, firstSec=beginTime, integrationTime=expTime, fields=['peakHeights']) nPhoton = len(tpl['timestamps']) if nPhoton>0: print "offset=",offset,"row=",iRow, "iCol=",iCol, "nPhotons=",nPhoton