Ejemplo n.º 1
0
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),
Ejemplo n.º 2
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
Ejemplo n.º 3
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 )
Ejemplo n.º 4
0
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
Ejemplo n.º 5
0
# 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