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)