t_runInfo = rh.getFromFile(rootFile,"Garfield/RunInfo")
t_runInfo.GetEntry(0) #Only one entry
#Dump some print
gh.dumpRunInfo(t_runInfo)
#
# Event loop
#
#Get event tree
t_event = rh.getFromFile(rootFile,"Garfield/Events")
#Get events numbers to process
eventNums = gh.getEventNumsToProcess(t_event.GetEntries(),args.maxNumEvents,args.firstEvent,args.eventStep)
#Plot raw signal for this event #TODO REMOVE
#t_event.Draw("rawSignalCurrent:rawSignalTime")
#if args.pauseForPlots : raw_input("Press Enter to continue...")
#Loop over events
for i_evt in eventNums :
if i_evt % 100 == 0 : print "Event %i" % (i_evt) #TODO % done instead
#Step tree to current event
t_event.GetEntry(i_evt)
#Check the track
gh.checkTrack(t_event.trackTime, t_event.trackOrigin, t_event.trackDirection)
gh.dumpRunInfo(t_runInfo)
# Get num events the event display info was dumped for
nStoredEventDisplays = t_runInfo.nStoredEventDisplays
print "Run stored event display info for %i events" % (nStoredEventDisplays)
print ""
#
# Event loop
#
# Get event tree
t_event = rh.getFromFile(rootFile, "Garfield/Events")
# Get events numbers to process
eventNums = gh.getEventNumsToProcess(nStoredEventDisplays, args.maxNumEvents, args.firstEvent, args.eventStep)
print eventNums
# Loop over events
for i_evt in eventNums:
print "\n\n---------------------------------------------"
print "Event %i :" % (i_evt)
# Step tree to current event
t_event.GetEntry(i_evt)
# Check the track
gh.checkTrack(t_event.trackTime, t_event.trackOrigin, t_event.trackDirection)
