jet_maps = pickle.load(infile)
print 'build jet by jet map from flat tree'
flat_map = {}
inputFile = io.root_open('trees/CombinedSV_ALL.root')
flat_tree = inputFile.tree
for entry in flat_tree:
evtid = (entry.run, entry.lumi, entry.evt)
if evtid not in jet_maps:
continue
if (entry.jetPt, entry.jetEta) not in jet_maps[evtid]:
continue
if evtid not in flat_map: flat_map[evtid] = {}
flat_map[evtid][(entry.jetPt, entry.jetEta)] = Struct.from_entry(entry)
handle = Handle('std::vector<pat::Jet>')
vtx_handle = Handle('vector<reco::Vertex>')
print 'analyzing pat output'
same_mva, same_inputs = 0, 0
n_analyzed_jets = 0
different_jets = {}
for evt in events:
evtid = (evt.eventAuxiliary().run(), evt.eventAuxiliary().luminosityBlock(), evt.eventAuxiliary().event())
evt.getByLabel('selectedPatJets', handle)
jets = handle.product()
ext_jets = []
for jet in jets:
