objects.JetAnalyzer.JetAnalyzer Class Reference

Inheritance diagram for objects.JetAnalyzer.JetAnalyzer:

Public Member Functions
def	__init__
def	beginLoop
def	computeQGvars
def	declareHandles
def	jetFlavour
def	matchJets
def	process
def	smearJets
def	testJetID
def	testJetNoID

Public Attributes
	allJetsUsedForMET
	bqObjects
	cqObjects
	deltaMetFromJEC
	Read jets, if necessary recalibrate and shift MET. More...
	doJEC
	genJets
	heaviestQCDFlavour
	jetGammaDR
	jetLepArbitration
	jetLepDR
	jetReCalibrator
	lepPtMin
	partons
	qglcalc
	rho
	shiftJEC
	shiftJER

Static Public Attributes
	alwaysCleanPhotons = False,
	class_object = JetAnalyzer,
	cleanGenJetsFromPhoton = False,
	cleanJetsFromFirstPhoton = False,
	cleanJetsFromIsoTracks = False,
	cleanJetsFromTaus = False,
	cleanSelectedLeptons = True,
string	collectionPostFix = ""
	copyJetsByValue = False,
	doPuId = False,
	doQG = False,
string	genJetCol = 'slimmedGenJets'
string	jecPath = ""
list	jet = jlpairs[lep]
string	jetCol = 'slimmedJets'
float	jetEta = 4.7
float	jetEtaCentral = 2.4
tuple	jetLepArbitration = (lambda jet,lepton : lepton)
float	jetLepDR = 0.4
int	jetPt = 25
tuple	jlpairs = matchObjectCollection( leptons, allJets, self.jetLepDR**2)
tuple	jtaupairs = matchObjectCollection( taus, allJets, self.jetLepDR**2)
	lepSelCut = lambdalep:True,
list	leptons = []
	Sort Jets by pT. More...
tuple	leptons = event.inclusiveLeptonsifhasattr(event, 'inclusiveLeptons')
	Associate jets to leptons. More...
int	minLepPt = 10
list	photons = []
	Clean Jets from photons. More...
	recalibrateJets = False,
string	recalibrationType = "AK4PFchs"
	relaxJetId = False,
tuple	rho = ('fixedGridRhoFastjetAll','','')
int	shiftJEC = 0
int	shiftJER = 0
	smearJets = True,
tuple	taus = getattr(event,'selectedTaus',[])
	Associate jets to taus. More...

Detailed Description

Taken from RootTools.JetAnalyzer, simplified, modified, added corrections    

Definition at line 49 of file JetAnalyzer.py.

Constructor & Destructor Documentation

def objects.JetAnalyzer.JetAnalyzer.__init__	(		self,
			cfg_ana,
			cfg_comp,
			looperName
	)

Definition at line 51 of file JetAnalyzer.py.

    def __init__(self, cfg_ana, cfg_comp, looperName):
        super(JetAnalyzer,self).__init__(cfg_ana, cfg_comp, looperName)
        mcGT   = cfg_ana.mcGT   if hasattr(cfg_ana,'mcGT')   else "PHYS14_25_V2"
        dataGT = cfg_ana.dataGT if hasattr(cfg_ana,'dataGT') else "GR_70_V2_AN1"
        self.shiftJEC = self.cfg_ana.shiftJEC if hasattr(self.cfg_ana, 'shiftJEC') else 0
        self.recalibrateJets = self.cfg_ana.recalibrateJets
        if   self.recalibrateJets == "MC"  : self.recalibrateJets =     self.cfg_comp.isMC
        elif self.recalibrateJets == "Data": self.recalibrateJets = not self.cfg_comp.isMC
        elif self.recalibrateJets not in [True,False]: raise RuntimeError, "recalibrateJets must be any of { True, False, 'MC', 'Data' }, while it is %r " % self.recalibrateJets
        self.doJEC = self.recalibrateJets or (self.shiftJEC != 0)
        if self.doJEC:
          if self.cfg_comp.isMC:
            self.jetReCalibrator = JetReCalibrator(mcGT,self.cfg_ana.recalibrationType, False,cfg_ana.jecPath)
          else:
            self.jetReCalibrator = JetReCalibrator(dataGT,self.cfg_ana.recalibrationType, True,cfg_ana.jecPath)
        self.doPuId = getattr(self.cfg_ana, 'doPuId', True)
        self.jetLepDR = getattr(self.cfg_ana, 'jetLepDR', 0.4)
        self.jetLepArbitration = getattr(self.cfg_ana, 'jetLepArbitration', lambda jet,lepton: lepton) 
        self.lepPtMin = getattr(self.cfg_ana, 'minLepPt', -1)
        self.lepSelCut = getattr(self.cfg_ana, 'lepSelCut', lambda lep : True)
        self.jetGammaDR =  getattr(self.cfg_ana, 'jetGammaDR', 0.4)
        if(self.cfg_ana.doQG):
            qgdefname="{CMSSW_BASE}/src/PhysicsTools/Heppy/data/pdfQG_AK4chs_antib_13TeV_v1.root"
            self.qglcalc = QGLikelihoodCalculator(getattr(self.cfg_ana,"QGpath",qgdefname).format(CMSSW_BASE= os.environ['CMSSW_BASE']))
        if not hasattr(self.cfg_ana ,"collectionPostFix"):self.cfg_ana.collectionPostFix=""