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

	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,

	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

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)
         self.doJEC = self.recalibrateJets or (self.shiftJEC != 0)
         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 = 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)
         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=""