| 1 |
"""Python helper tools for CMS FWLite
|
| 2 |
|
| 3 |
benedikt.hegner@cern.ch
|
| 4 |
|
| 5 |
"""
|
| 6 |
import re
|
| 7 |
import ROOT
|
| 8 |
import exceptions
|
| 9 |
### define tab completion
|
| 10 |
try:
|
| 11 |
import readline, cmscompleter
|
| 12 |
readline.parse_and_bind('tab: complete')
|
| 13 |
except:
|
| 14 |
print 'WARNING: Could not load tab completion'
|
| 15 |
|
| 16 |
|
| 17 |
# for adding iterators at runtime
|
| 18 |
import iterators
|
| 19 |
|
| 20 |
|
| 21 |
### workaround iterator generators for ROOT classes
|
| 22 |
def all(container):
|
| 23 |
|
| 24 |
# loop over ROOT::TTree and similar
|
| 25 |
if hasattr(container,'GetEntries'):
|
| 26 |
try:
|
| 27 |
entries = container.GetEntries()
|
| 28 |
for entry in xrange(entries):
|
| 29 |
yield entry
|
| 30 |
except:
|
| 31 |
raise cmserror("Looping of %s failed" %container)
|
| 32 |
|
| 33 |
# loop over std::vectors and similar
|
| 34 |
elif hasattr(container, 'size'):
|
| 35 |
try:
|
| 36 |
entries = container.size()
|
| 37 |
for entry in xrange(entries):
|
| 38 |
yield container[entry]
|
| 39 |
except:
|
| 40 |
pass
|
| 41 |
|
| 42 |
# loop over containers with begin and end iterators
|
| 43 |
def loop(begin, end):
|
| 44 |
"""Convert a pair of C++ iterators into a python generator"""
|
| 45 |
while (begin != end):
|
| 46 |
yield begin.__deref__() #*b
|
| 47 |
begin.__preinc__() #++b
|
| 48 |
|
| 49 |
### auto branch types (Chris Jones)
|
| 50 |
def createBranchBuffer(branch):
|
| 51 |
reColons = re.compile(r'::')
|
| 52 |
reCloseTemplate =re.compile(r'>')
|
| 53 |
reOpenTemplate =re.compile(r'<')
|
| 54 |
branchType = ROOT.branchToClass(branch)
|
| 55 |
#buffer = eval ('ROOT.'+reColons.sub(".",reOpenTemplate.sub("(ROOT.",reCloseTemplate.sub(")",branchType.GetName())))+'()')
|
| 56 |
buffer = ROOT.MakeRootClass(branchType.GetName()) ()
|
| 57 |
if( branch.GetName()[-1] != '.'):
|
| 58 |
branch.SetAddress(buffer)
|
| 59 |
else:
|
| 60 |
branch.SetAddress(ROOT.AddressOf(buffer))
|
| 61 |
return buffer
|
| 62 |
|
| 63 |
|
| 64 |
class EventTree(object):
|
| 65 |
def __init__(self,obj):
|
| 66 |
if isinstance(obj, ROOT.TTree):
|
| 67 |
self._tree = obj
|
| 68 |
elif isinstance(obj, ROOT.TFile):
|
| 69 |
self._tree = obj.Get("Events")
|
| 70 |
elif isinstance(obj, str):
|
| 71 |
self._tree = ROOT.TFile.Open(obj).Get("Events")
|
| 72 |
else:
|
| 73 |
raise cmserror("EventTree accepts only TTrees, TFiles and filenames")
|
| 74 |
self._usedBranches = dict()
|
| 75 |
self._index = -1
|
| 76 |
self._aliases = self._tree.GetListOfAliases()
|
| 77 |
def branch(self,name):
|
| 78 |
# support for aliases
|
| 79 |
alias = self._tree.GetAlias(name)
|
| 80 |
if alias != '': name = alias
|
| 81 |
# access the branch in ttree
|
| 82 |
if name in self._usedBranches:
|
| 83 |
return self._usedBranches[name]
|
| 84 |
self._usedBranches[name]=EventBranch(self,name)
|
| 85 |
return self._usedBranches[name]
|
| 86 |
def cppCode(self, name):
|
| 87 |
"""C++ code for accessing the product inside the full framework"""
|
| 88 |
alias = self._tree.GetAlias(name)
|
| 89 |
if alias != '': name = alias
|
| 90 |
tmpBranch = self._tree.GetBranch(name)
|
| 91 |
typeString = ROOT.branchToClass(tmpBranch).GetName()
|
| 92 |
if "edm::Wrapper" in typeString:
|
| 93 |
typeString = typeString.replace("<edm::Wrapper","")
|
| 94 |
typeString = typeString.rstrip(">")
|
| 95 |
nameParts = name.split("_")
|
| 96 |
if nameParts[2] == "":
|
| 97 |
cppCode = 'edm::Handle<%s > dummy;\nevent.getByLabel("%s", dummy);'\
|
| 98 |
%(typeString, nameParts[1])
|
| 99 |
else:
|
| 100 |
cppCode = 'edm::Handle<%s > dummy;\nevent.getByLabel("%s", "%s", dummy);'\
|
| 101 |
%(typeString, nameParts[1], nameParts[2])
|
| 102 |
return cppCode
|
| 103 |
def getListOfAliases(self):
|
| 104 |
return self._aliases
|
| 105 |
def index(self):
|
| 106 |
return self._index
|
| 107 |
def tree(self):
|
| 108 |
return self._tree
|
| 109 |
def __setBranchIndicies(self):
|
| 110 |
for branch in self._usedBranches.itervalues():
|
| 111 |
branch.setIndex(self._index)
|
| 112 |
def __getattr__(self, name):
|
| 113 |
return self.branch(name)
|
| 114 |
def __getitem__(self,key):
|
| 115 |
if key <0 or key > self._tree.GetEntries():
|
| 116 |
raise IndexError
|
| 117 |
self._index = key
|
| 118 |
self.__setBranchIndicies()
|
| 119 |
return self
|
| 120 |
def __iter__(self):
|
| 121 |
# flushing/initializing the root buffers
|
| 122 |
entry = 0
|
| 123 |
self._index = entry
|
| 124 |
self.__setBranchIndicies()
|
| 125 |
self._tree.GetEntry(self._index,0)
|
| 126 |
# the real loop
|
| 127 |
for entry in xrange(self._tree.GetEntries()):
|
| 128 |
self._index = entry
|
| 129 |
self.__setBranchIndicies()
|
| 130 |
self._tree.GetEntry(self._index,0)
|
| 131 |
yield Event(self) # TODO: don't return a new object but update the old one
|
| 132 |
|
| 133 |
|
| 134 |
class Event(object):
|
| 135 |
def __init__(self, eventTree):
|
| 136 |
self._eventTree = eventTree
|
| 137 |
|
| 138 |
def getProduct(self, name):
|
| 139 |
return iterators.addIterator(self._eventTree.branch(name)())
|
| 140 |
|
| 141 |
def __getattr__(self, name):
|
| 142 |
return iterators.addIterator(self._eventTree.branch(name)())
|
| 143 |
|
| 144 |
|
| 145 |
class EventBranch(object):
|
| 146 |
def __init__(self,parent,name):
|
| 147 |
self._branch = parent.tree().GetBranch(name)
|
| 148 |
if self._branch == None:
|
| 149 |
raise cmserror("Unknown branch "+name)
|
| 150 |
self._buffer = createBranchBuffer(self._branch)
|
| 151 |
self._index = parent.index()
|
| 152 |
self._readData = False
|
| 153 |
def setIndex(self,index):
|
| 154 |
self._index = index
|
| 155 |
self._readData = False
|
| 156 |
def __readData(self):
|
| 157 |
self._branch.GetEntry(self._index)
|
| 158 |
self._readData = True
|
| 159 |
|
| 160 |
# replace this by __getattr__ to allow branch.attr instead of branch().attr
|
| 161 |
def __call__(self):
|
| 162 |
if not self._readData:
|
| 163 |
self.__readData()
|
| 164 |
return self._buffer
|
| 165 |
|
| 166 |
|
| 167 |
class cmserror(exceptions.StandardError):
|
| 168 |
def __init__(self, message):
|
| 169 |
length = len(message)+7 #7=len("ERROR: ")
|
| 170 |
print "="*length
|
| 171 |
print "ERROR:", message
|
| 172 |
print "="*length
|
