source: svn/trunk/python/SimpleEventSelection.py

#very simple EventSelection class aimed at demonstrating the 
#typical implementation of an EventSelection class

#On purpose, the categories are not in any "logical" order, and there is some redundant check when testing the category.
#This is mostly for illustration.

# requirements:
#   event.muons
#   event.electrons
#   event.jets

# the list of category names
categoryNames = [ "MuonChannel/SingleMuon", "ElectronChannel/SingleElectron", "ElectronChannel/Jet", "MuonChannel/DoubleMuon", "ElectronChannel/DoubleElectron" ]

def eventCategory(event):
  """Check analysis requirements for various steps
     and return a tuple of data used to decide 
     to what category an event belong """
  categoryData = [ ]
  # 0: number of muons
  categoryData.append(event.muons.GetEntries())
  # 1: number of electrons
  categoryData.append(event.electrons.GetEntries())
  # 2: number of jets
  categoryData.append(event.jets.GetEntries())
  # 3: Pt of the leading muon > 2GeV
  if event.muons.GetEntries():
    categoryData.append(event.muons[0].PT>1.)
  else:
    categoryData.append(False)
  # 4: Pt of the leading electron > 2 GeV
  if event.electrons.GetEntries():
    categoryData.append(event.electrons[0].PT>1.)
  else:
    categoryData.append(False)
  # DONE
  return categoryData

def isInCategory(category, categoryData):
  """Check if the event enters category X, given the tuple computed by eventCategory."""
  if category==0:
    return categoryData[0]>0 
  elif category==1:
    return categoryData[1]>0
  elif category==2:
    return categoryData[2]>0
  elif category==3:
    return isInCategory(0,categoryData) and categoryData[3]==True
  elif category==4:
    return isInCategory(1,categoryData) and categoryData[4]==True
  else:
    return False