import FWCore.ParameterSet.Config as cms
from DQMOffline.JetMET.jetMETDQMCleanup_cff import *
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5CaloL2L3, ak5CaloL2Relative, ak5CaloL3Absolute
newAk5CaloL2L3 = ak5CaloL2L3.clone()
from JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff import ak7CaloL2L3, ak7CaloL2Relative, ak7CaloL3Absolute
newAk7CaloL2L3 = ak7CaloL2L3.clone()
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5PFL1FastL2L3, ak5PFL1Fastjet, ak5PFL2Relative, ak5PFL3Absolute
newAk5PFL1FastL2L3 = ak5PFL1FastL2L3.clone()
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5JPTL1FastL2L3, ak5JPTL1Fastjet, ak5JPTL2Relative, ak5JPTL3Absolute
newAk5JPTL1FastL2L3 = ak5JPTL1FastL2L3.clone()
tcMetDQMAnalyzer = cms.EDAnalyzer(
"METAnalyzer",
OutputMEsInRootFile=cms.bool(False),
OutputFile=cms.string('jetMETMonitoring.root'),
METType=cms.untracked.string('tc'),
METCollectionLabel=cms.InputTag("tcMet"),
JetCollectionLabel=cms.InputTag("JetPlusTrackZSPCorJetAntiKt5"),
JetCorrections=cms.string("newAk5JPTL1FastL2L3"),
InputJetIDValueMap=cms.InputTag("ak5JetID"),
ptThreshold=cms.double(30),
FolderName=cms.untracked.string("JetMET/MET/"),
CleaningParameters=cleaningParameters.clone(),
TriggerResultsLabel=cms.InputTag("TriggerResults::HLT"),
runcosmics=cms.untracked.bool(False),
LSBegin=cms.int32(0),
ak5PFL1FastjetNoPU = ak5PFL1Fastjet.clone(algorithm="AK5PFchs",
srcRho=cms.InputTag(
'kt6PFJetsNoPU', 'rho'))
# L2 (relative eta-conformity) Correction Services
ak5PFL2RelativeNoPU = ak5CaloL2Relative.clone(algorithm='AK5PFchs')
# L3 (absolute) Correction Services
ak5PFL3AbsoluteNoPU = ak5CaloL3Absolute.clone(algorithm='AK5PFchs')
# Multiple corrections
ak5PFL2L3NoPU = ak5PFL2L3.clone(correctors=('ak5PFL2RelativeNoPU',
'ak5PFL3AbsoluteNoPU'))
ak5PFL1FastL2L3NoPU = ak5PFL1FastL2L3.clone(correctors=('ak5PFL1FastjetNoPU',
'ak5PFL2RelativeNoPU',
'ak5PFL3AbsoluteNoPU'))
from CommonTools.ParticleFlow.pfNoPileUp_cff import *
pfPileUp.checkClosestZVertex = cms.bool(False)
pfPileUp.Vertices = cms.InputTag('offlinePrimaryVertices')
pfPileUp.PFCandidates = cms.InputTag("particleFlow")
# Jets
from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets
ak5PFJetsNoPU = ak5PFJets.clone(src=cms.InputTag("pfNoPileUp"),
doAreaFastjet=cms.bool(True))
ak7PFJetsNoPU = ak5PFJetsNoPU.clone(rParam=0.7)
# MET
# )
ak5PFL1FastjetNoPU = ak5PFL1Fastjet.clone(
algorithm = "AK5PFchs",
srcRho = cms.InputTag('kt6PFJetsNoPU', 'rho')
)
# L2 (relative eta-conformity) Correction Services
ak5PFL2RelativeNoPU = ak5CaloL2Relative.clone( algorithm = 'AK5PFchs' )
# L3 (absolute) Correction Services
ak5PFL3AbsoluteNoPU = ak5CaloL3Absolute.clone( algorithm = 'AK5PFchs' )
# Multiple corrections
ak5PFL2L3NoPU = ak5PFL2L3.clone( correctors = ('ak5PFL2RelativeNoPU','ak5PFL3AbsoluteNoPU') )
ak5PFL1FastL2L3NoPU = ak5PFL1FastL2L3.clone( correctors = ('ak5PFL1FastjetNoPU','ak5PFL2RelativeNoPU','ak5PFL3AbsoluteNoPU') )
from CommonTools.ParticleFlow.pfNoPileUp_cff import *
pfPileUp.checkClosestZVertex = cms.bool(False)
pfPileUp.Vertices = cms.InputTag('offlinePrimaryVertices')
pfPileUp.PFCandidates = cms.InputTag("particleFlow")
# Jets
from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets
ak5PFJetsNoPU = ak5PFJets.clone(
src = cms.InputTag("pfNoPileUp"),
doAreaFastjet = cms.bool(True)
)
ak7PFJetsNoPU = ak5PFJetsNoPU.clone( rParam = 0.7 )
import FWCore.ParameterSet.Config as cms
#not useful anymore for b-tagging but used in some other sequences
#from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5PFL2L3,ak5PFL2Relative,ak5PFL3Absolute
#JEC for CHS
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5PFL1Fastjet, ak5PFL2Relative, ak5PFL3Absolute, ak5PFResidual, ak5PFL1FastL2L3, ak5PFL1FastL2L3Residual
newak5PFCHSL1Fastjet = ak5PFL1Fastjet.clone(algorithm = 'AK5PFchs')
newak5PFCHSL2Relative = ak5PFL2Relative.clone(algorithm = 'AK5PFchs')
newak5PFCHSL3Absolute = ak5PFL3Absolute.clone(algorithm = 'AK5PFchs')
newak5PFCHSResidual = ak5PFResidual.clone(algorithm = 'AK5PFchs')
newak5PFCHSL1FastL2L3 = ak5PFL1FastL2L3.clone(correctors = cms.vstring('newak5PFCHSL1Fastjet','newak5PFCHSL2Relative','newak5PFCHSL3Absolute'))
newak5PFCHSL1FastL2L3Residual = ak5PFL1FastL2L3Residual.clone(correctors = cms.vstring('newak5PFCHSL1Fastjet','newak5PFCHSL2Relative','newak5PFCHSL3Absolute','newak5PFCHSResidual'))
######### DATA ############
from DQMOffline.RecoB.bTagAnalysisData_cfi import *
bTagAnalysis.finalizePlots = False
bTagAnalysis.finalizeOnly = False
bTagAnalysis.ptRanges = cms.vdouble(0.0)
#Residual correction will be added inside the c++ code only for data (checking the presence of genParticles collection), not explicit here as this sequence also ran on MC FullSim
bTagAnalysis.doJetID = True
bTagAnalysis.doJEC = True
bTagAnalysis.JECsource = cms.string("newak5PFCHSL1FastL2L3")
bTagPlotsDATA = cms.Sequence(bTagAnalysis)
########## MC ############
#Matching
from PhysicsTools.JetMCAlgos.CaloJetsMCFlavour_cfi import *
AK5byRef.jets = cms.InputTag("ak5PFJetsCHS")
#Get gen jet collection for real jets
# Date : Mar 17 2009
# Description : cff file for DQM offline vladation
from RecoJets.Configuration.RecoJetAssociations_cff import *
from Validation.RecoJets.JetValidation_cfi import *
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5CaloL2L3,ak5CaloL2Relative,ak5CaloL3Absolute
newAk5CaloL2L3 = ak5CaloL2L3.clone()
#from JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff import ak7CaloL2L3,ak7CaloL2Relative,ak7CaloL3Absolute
#newAk7CaloL2L3 = ak7CaloL2L3.clone()
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5PFL1FastL2L3,ak5PFL1Fastjet,ak5PFL2Relative,ak5PFL3Absolute
newAk5PFL1FastL2L3 = ak5PFL1FastL2L3.clone()
from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak5JPTL1FastL2L3,ak5JPTL1Fastjet,ak5JPTL2Relative,ak5JPTL3Absolute
newAk5JPTL1FastL2L3 = ak5JPTL1FastL2L3.clone()
from JetCorrectionServices_AK5CHS_cff import ak5PFchsL1FastL2L3,ak5PFchsL1Fastjet,ak5PFchsL2Relative,ak5PFchsL3Absolute
newAk5PFchsL1FastL2L3 = ak5PFchsL1FastL2L3.clone()
JetValidation = cms.Sequence(
# JetAnalyzerKt6PF*
# JetAnalyzerKt6Calo*
JetAnalyzerAk5Calo*
# JetAnalyzerAk7Calo*
JetAnalyzerAk5PF*
JetAnalyzerAk5JPT*
JetAnalyzerAk5PFCHS