objects.JetAnalyzer Namespace Reference

Classes
class	JetAnalyzer

Functions
def	cleanJetsAndLeptons (jets, leptons, deltaR, arbitration)
def	cleanNearestJetOnly (jets, leptons, deltaR)
def	shiftJERfactor (JERShift, aeta)

Function Documentation

def objects.JetAnalyzer.cleanJetsAndLeptons	(		jets,
			leptons,
			deltaR,
			arbitration
	)

Definition at line 24 of file JetAnalyzer.py.

References reco.deltaR2().

def cleanJetsAndLeptons(jets,leptons,deltaR,arbitration):
    dr2 = deltaR**2
    goodjet = [ True for j in jets ]
    goodlep = [ True for l in leptons ]
    for il, l in enumerate(leptons):
        ibest, d2m = -1, dr2
        for i,j in enumerate(jets):
            d2i = deltaR2(l.eta(),l.phi(), j.eta(),j.phi())
            if d2i < dr2:
                choice = arbitration(j,l)
                if choice == j:
                   # if the two match, and we prefer the jet, then drop the lepton and be done
                   goodlep[il] = False
                   break 
                elif choice == (j,l) or choice == (l,j):
                   # asked to keep both, so we don't consider this match
                   continue
            if d2i < d2m:
                ibest, d2m = i, d2i
        # this lepton has been killed by a jet, then we clean the jet that best matches it
        if not goodlep[il]: continue 
        if ibest != -1: goodjet[ibest] = False
    return ( [ j for (i ,j) in enumerate(jets)    if goodjet[i ] == True ], 
             [ l for (il,l) in enumerate(leptons) if goodlep[il] == True ] )

def objects.JetAnalyzer.cleanNearestJetOnly	(		jets,
			leptons,
			deltaR
	)

Definition at line 12 of file JetAnalyzer.py.

References reco.deltaR2().

def cleanNearestJetOnly(jets,leptons,deltaR):
    dr2 = deltaR**2
    good = [ True for j in jets ]
    for l in leptons:
        ibest, d2m = -1, dr2
        for i,j in enumerate(jets):
            d2i = deltaR2(l.eta(),l.phi(), j.eta(),j.phi())
            if d2i < d2m:
                ibest, d2m = i, d2i
        if ibest != -1: good[ibest] = False
    return [ j for (i,j) in enumerate(jets) if good[i] == True ] 

def objects.JetAnalyzer.shiftJERfactor	(		JERShift,
			aeta
	)

Definition at line 51 of file JetAnalyzer.py.

def shiftJERfactor(JERShift, aeta):
        factor = 1.079 + JERShift*0.026
        if   aeta > 3.2: factor = 1.056 + JERShift * 0.191
        elif aeta > 2.8: factor = 1.395 + JERShift * 0.063
        elif aeta > 2.3: factor = 1.254 + JERShift * 0.062
        elif aeta > 1.7: factor = 1.208 + JERShift * 0.046
        elif aeta > 1.1: factor = 1.121 + JERShift * 0.029
        elif aeta > 0.5: factor = 1.099 + JERShift * 0.028
        return factor