objects.JetAnalyzer.JetAnalyzer Class Reference

Inheritance diagram for objects.JetAnalyzer.JetAnalyzer:

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

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

Static Public Attributes
	addJECShifts = False,
int	addJERShifts = 0
	alwaysCleanPhotons = False,
string	applyL2L3Residual = 'Data'
	attachNeutrinos = True,
	calculateSeparateCorrections = False,
	calculateType1METCorrection = False,
	checkLeptonPFOverlap = True,
	class_object = JetAnalyzer,
	cleanFromLepAndGammaSimultaneously = False,
	cleanGenJetsFromPhoton = False,
	cleanJetsFromFirstPhoton = False,
	cleanJetsFromIsoTracks = False,
	cleanJetsFromTaus = False,
	cleanSelectedLeptons = True,
string	collectionPostFix = ""
	copyJetsByValue = False,
	do_mc_match = True,
	doPuId = False,
	doQG = False,
string	genJetCol = 'slimmedGenJets'
	genNuSelection = lambdanu:True,
tuple	has_overlaps = p1.key()
tuple	incleptons = event.inclusiveLeptonsifhasattr(event, 'inclusiveLeptons')
	Jet Id, after jet/lepton cleaning. More...
string	jecPath = ""
list	jet = jlpairs[lep]
string	jetCol = 'slimmedJets'
float	jetEta = 4.7
float	jetEtaCentral = 2.4
float	jetGammaDR = 0.4
float	jetGammaLepDR = 0.4
tuple	jetLepArbitration = (lambda jet,lepton : lepton)
float	jetLepDR = 0.4
list	jetNus = [x for x in event.genParticles if abs(x.pdgId()) in [12,14,16] and self.cfg_ana.genNuSelection(x) ]
int	jetPt = 25
list	jetsEtaCut = [j for j in self.jets if abs(j.eta()) < self.cfg_ana.jetEta ]
tuple	jlpairs = matchObjectCollection(incleptons, allJets, self.jetLepDR**2)
tuple	jtaupairs = matchObjectCollection( taus, allJets, self.jetLepDR**2)
list	leps_with_overlaps = []
	Apply jet selection. More...
	lepSelCut = lambdalep:True,
list	leptons = []
	Sort Jets by pT. More...
int	minLepPt = 10
tuple	p1 = jet.sourceCandidatePtr(i)
tuple	p2 = lep.sourceCandidatePtr(j)
tuple	pairs = matchObjectCollection(jetNus, self.genJets, 0.4**2)
list	photons = []
	First cleaning, then Jet Id. 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...
dictionary	type1METParams = { 'jetPtThreshold':15., 'skipEMfractionThreshold':0.9, 'skipMuons':True }

Detailed Description

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

Definition at line 65 of file JetAnalyzer.py.

Constructor & Destructor Documentation

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

Definition at line 67 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
        self.addJECShifts = self.cfg_ana.addJECShifts if hasattr(self.cfg_ana, 'addJECShifts') else 0
        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)
       
        calculateSeparateCorrections = getattr(cfg_ana,"calculateSeparateCorrections", False);
        calculateType1METCorrection  = getattr(cfg_ana,"calculateType1METCorrection",  False);
        self.doJEC = self.recalibrateJets or (self.shiftJEC != 0) or self.addJECShifts or calculateSeparateCorrections or calculateType1METCorrection
        if self.doJEC:
          doResidual = getattr(cfg_ana, 'applyL2L3Residual', 'Data')
          if   doResidual == "MC":   doResidual = self.cfg_comp.isMC
          elif doResidual == "Data": doResidual = not self.cfg_comp.isMC
          elif doResidual not in [True,False]: raise RuntimeError("If specified, applyL2L3Residual must be any of { True, False, 'MC', 'Data'(default)}")
          GT = getattr(cfg_comp, 'jecGT', mcGT if self.cfg_comp.isMC else dataGT)
          # Now take care of the optional arguments
          kwargs = { 'calculateSeparateCorrections':calculateSeparateCorrections,
                     'calculateType1METCorrection' :calculateType1METCorrection, }
          if kwargs['calculateType1METCorrection']: kwargs['type1METParams'] = cfg_ana.type1METParams
          # instantiate the jet re-calibrator
          self.jetReCalibrator = JetReCalibrator(GT, cfg_ana.recalibrationType, doResidual, cfg_ana.jecPath, **kwargs)
        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)
        self.jetGammaLepDR =  getattr(self.cfg_ana, 'jetGammaLepDR', 0.4)
        self.cleanFromLepAndGammaSimultaneously = getattr(self.cfg_ana, 'cleanFromLepAndGammaSimultaneously', False)
        if self.cleanFromLepAndGammaSimultaneously:
            if hasattr(self.cfg_ana, 'jetGammaLepDR'):
                self.jetGammaLepDR =  self.jetGammaLepDR 
            elif (self.jetGammaDR == self.jetLepDR):
                self.jetGammaLepDR = self.jetGammaDR
            else:
                raise RuntimeError("DR for simultaneous cleaning of jets from leptons and photons is not defined, and dR(gamma, jet)!=dR(lep, jet)")
        if(self.cfg_ana.doQG):
            qgdefname="{CMSSW_BASE}/src/PhysicsTools/Heppy/data/pdfQG_AK4chs_13TeV_v2b.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=""