JetMET/python/jetAnalyzer_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from DQMOffline.JetMET.jetDQMConfig_cff import *      # parameters for all jet analyzers
0004 from DQMOffline.JetMET.jetMETDQMCleanup_cff import *  # parameters for event cleanup
0005
0006 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0007 jetDQMAnalyzerAk4CaloUncleaned = DQMEDAnalyzer('JetAnalyzer',
0008     JetType = cms.string('calo'),#pf, calo or jpt
0009     JetCorrections = cms.InputTag("dqmAk4CaloL2L3ResidualCorrector"),
0010     jetsrc = cms.InputTag("ak4CaloJets"),
0011     METCollectionLabel     = cms.InputTag("caloMet"),
0012     muonsrc = cms.InputTag("muons"),
0013     l1algoname = cms.string("L1Tech_BPTX_plus_AND_minus.v0"),
0014     filljetHighLevel =cms.bool(False),
0015     fillsubstructure =cms.bool(False),
0016     ptMinBoosted = cms.double(400.),
0017     #
0018     #
0019     #
0020     highPtJetTrigger = cms.PSet(
0021         andOr         = cms.bool( False ),
0022         dbLabel        = cms.string("JetMETDQMTrigger"),
0023         hltInputTag    = cms.InputTag( "TriggerResults::HLT" ),
0024         hltPaths       = cms.vstring( 'HLT_PFJet450_v*'),
0025         andOrHlt       = cms.bool( True ),
0026         errorReplyHlt  = cms.bool( False ),
0027     ),
0028     lowPtJetTrigger = cms.PSet(
0029         andOr         = cms.bool( False ),
0030         dbLabel        = cms.string("JetMETDQMTrigger"),
0031         hltInputTag    = cms.InputTag( "TriggerResults::HLT" ),
0032         hltPaths       = cms.vstring( 'HLT_PFJet80_v*'),
0033         andOrHlt       = cms.bool( True ),
0034         errorReplyHlt  = cms.bool( False ),
0035     ),
0036
0037     TriggerResultsLabel        = cms.InputTag("TriggerResults::HLT"),
0038     processname                = cms.string("HLT"),
0039
0040     #
0041     # Jet-related
0042     #
0043
0044     JetCleaningFlag            = cms.untracked.bool(False),
0045
0046     runcosmics                 = cms.untracked.bool(False),
0047
0048     #Cleanup parameters
0049     CleaningParameters = cleaningParameters.clone(
0050         bypassAllPVChecks = True,
0051         ),
0052
0053     #for JPT and CaloJetID
0054     InputJetIDValueMap         = cms.InputTag("ak4JetID"),
0055     #options for Calo and JPT: LOOSE,LOOSE_AOD,TIGHT,MINIMAL
0056     #for PFJets: LOOSE,TIGHT
0057     JetIDQuality               = cms.string("LOOSE"),
0058     #options for Calo and JPT: PURE09,DQM09,CRAFT08
0059     #for PFJets: FIRSTDATA
0060     JetIDVersion               = cms.string("PURE09"),
0061     #
0062     #actually done only for PFJets at the moment
0063     InputMVAPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorDQM","fullDiscriminant"),
0064     InputCutPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorDQM","cutbasedDiscriminant"),
0065     InputMVAPUIDValue = cms.InputTag("pileupJetIdEvaluatorDQM","fullId"),
0066     InputCutPUIDValue = cms.InputTag("pileupJetIdEvaluatorDQM","cutbasedId"),
0067
0068     InputQGMultiplicity = cms.InputTag("QGTagger", "mult"),
0069     InputQGLikelihood = cms.InputTag("QGTagger", "qgLikelihood"),
0070     InputQGPtDToken = cms.InputTag("QGTagger", "ptD"),
0071     InputQGAxis2 = cms.InputTag("QGTagger", "axis2"),
0072
0073     fillCHShistos =cms.bool(False),
0074     #
0075     # For jetAnalysis
0076     #
0077     jetAnalysis = jetDQMParameters.clone(),
0078
0079     #
0080     # DCS
0081     #
0082     DCSFilterForJetMonitoring = cms.PSet(
0083       DetectorTypes = cms.untracked.string("ecal:hbhe:hf"),
0084       #DebugOn = cms.untracked.bool(True),
0085       alwaysPass = cms.untracked.bool(False)
0086     )
0087 )
0088
0089 jetDQMAnalyzerAk4CaloCleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
0090     JetCleaningFlag   = True,
0091     filljetHighLevel  = False,
0092     CleaningParameters = cleaningParameters.clone(
0093         bypassAllPVChecks = True,
0094     ),
0095     jetAnalysis=jetDQMParameters.clone(
0096         ptThreshold = 20.,
0097         asymmetryThirdJetCut = 30,
0098         balanceThirdJetCut   = 0.2,
0099        )
0100 )
0101
0102 jetDQMAnalyzerAk4PFUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
0103     CleaningParameters = cleaningParameters.clone(
0104        bypassAllPVChecks  = False,
0105     ),
0106     #for PFJets: LOOSE,TIGHT
0107     JetIDQuality               = cms.string("TIGHT"),
0108     #options for Calo and JPT: PURE09,DQM09,CRAFT08
0109     #for PFJets: RUN2ULCHS for 11_1_X onwards
0110     JetIDVersion               = cms.string("RUN2ULCHS"),
0111     JetType = cms.string('pf'),#pf, calo or jpt
0112     JetCorrections = cms.InputTag("dqmAk4PFL1FastL2L3ResidualCorrector"),
0113     jetsrc = cms.InputTag("ak4PFJets"),
0114     METCollectionLabel     = cms.InputTag("pfMet"),
0115     filljetHighLevel  = cms.bool(False),
0116     DCSFilterForJetMonitoring = cms.PSet(
0117       DetectorTypes = cms.untracked.string("ecal:hbhe:hf:pixel:sistrip:es:muon"),
0118       #DebugOn = cms.untracked.bool(True),
0119       alwaysPass = cms.untracked.bool(False)
0120     )
0121 )
0122
0123
0124 jetDQMAnalyzerAk4PFCleaned=jetDQMAnalyzerAk4PFUncleaned.clone(
0125     JetCleaningFlag = True,
0126     filljetHighLevel = False,
0127     jetAnalysis=jetDQMParameters.clone(
0128         ptThreshold = 20.,
0129         asymmetryThirdJetCut = 30,
0130         balanceThirdJetCut = 0.2,
0131         ),
0132     METCollectionLabel     = "pfMet"
0133 )
0134
0135 jetDQMAnalyzerAk4PFCHSCleaned=jetDQMAnalyzerAk4PFCleaned.clone(
0136     filljetHighLevel = True,
0137     JetCorrections = "dqmAk4PFCHSL1FastL2L3ResidualCorrector",
0138     jetsrc = "ak4PFJetsCHS",
0139     METCollectionLabel     = "pfMETT1",
0140     InputMVAPUIDDiscriminant = "pileupJetIdEvaluatorCHSDQM:fullDiscriminant",
0141     InputCutPUIDDiscriminant = "pileupJetIdEvaluatorCHSDQM:cutbasedDiscriminant",
0142     InputMVAPUIDValue = "pileupJetIdEvaluatorCHSDQM:fullId",
0143     InputCutPUIDValue = "pileupJetIdEvaluatorCHSDQM:cutbasedId",
0144     fillCHShistos = True
0145 )
0146
0147 jetDQMAnalyzerAk4PFCHSUncleanedMiniAOD=jetDQMAnalyzerAk4PFUncleaned.clone(
0148     filljetHighLevel = True,
0149     CleaningParameters = cleaningParameters.clone(
0150         vertexCollection  =  "goodOfflinePrimaryVerticesDQMforMiniAOD" ,
0151         ),
0152     JetType = 'miniaod',#pf, calo or jpt
0153     jetsrc = "slimmedJets",
0154     METCollectionLabel     = "slimmedMETs"
0155 )
0156
0157 jetDQMAnalyzerAk4PFCHSCleanedMiniAOD=jetDQMAnalyzerAk4PFCleaned.clone(
0158     CleaningParameters = cleaningParameters.clone(
0159         vertexCollection    =  "goodOfflinePrimaryVerticesDQMforMiniAOD"
0160         ),
0161     JetType = 'miniaod',#pf, calo or jpt
0162     jetsrc = "slimmedJets"
0163 )
0164
0165 jetDQMAnalyzerAk8PFPUPPICleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
0166     jetsrc = "slimmedJetsAK8",
0167     #for PUPPI jets: TIGHT
0168     JetIDQuality  = "TIGHT",
0169     #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
0170     JetIDVersion  = "RUN2ULPUPPI",
0171     fillsubstructure =True
0172 )
0173
0174 jetDQMAnalyzerAk4PFCHSPuppiCleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
0175     JetType = 'miniaod',#pf, calo or jpt
0176     jetsrc = "slimmedJetsPuppi",
0177     #for PUPPI jets: TIGHT
0178     JetIDQuality  = "TIGHT",
0179     #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
0180     JetIDVersion  = "RUN2ULPUPPI"
0181 )
0182
0183 jetDQMAnalyzerIC5CaloHIUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
0184     filljetHighLevel = True,
0185     CleaningParameters = cleaningParameters.clone(
0186         bypassAllPVChecks  = False,
0187         vertexCollection =  "hiSelectedVertex",
0188         ),
0189     JetType = 'calo',#pf, calo or jpt
0190     JetCorrections = "",# no jet correction available yet?
0191     jetsrc = "iterativeConePu5CaloJets",
0192     JetCleaningFlag            = False,
0193     runcosmics                 = True,
0194     DCSFilterForJetMonitoring = dict(
0195         DetectorTypes = "ecal:hbhe:hf:pixel:sistrip:es:muon",
0196         #DebugOn = True,
0197         alwaysPass = False
0198     )
0199 )
0200
0201 jetDQMAnalyzerAkPU3Calo = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0202                                          JetType = cms.untracked.string('calo'),
0203                                          UEAlgo = cms.untracked.string('Pu'),
0204                                          OutputFile = cms.untracked.string(''),
0205                                          src = cms.InputTag("akPu3CaloJets"),
0206                                          PFcands = cms.InputTag("particleFlowTmp"),
0207                                          Background = cms.InputTag("voronoiBackgroundCalo"),
0208                                          #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0209                                          centralitycollection = cms.InputTag("hiCentrality"),
0210                                          centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0211                                          JetCorrections = cms.string(""),
0212                                          recoJetPtThreshold = cms.double(10),
0213                                          RThreshold = cms.double(0.3),
0214                                          reverseEnergyFractionThreshold = cms.double(0.5)
0215 )
0216 jetDQMAnalyzerAkPU4Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu4CaloJets")
0217 jetDQMAnalyzerAkPU5Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu5CaloJets")
0218
0219 jetDQMAnalyzerAkPU3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0220                                        JetType = cms.untracked.string('pf'),
0221                                        UEAlgo = cms.untracked.string('Pu'),
0222                                        OutputFile = cms.untracked.string(''),
0223                                        src = cms.InputTag("akPu3PFJets"),
0224                                        PFcands = cms.InputTag("particleFlowTmp"),
0225                                        Background = cms.InputTag("voronoiBackgroundPF"),
0226                                        centralitycollection = cms.InputTag("hiCentrality"),
0227                                        centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0228                                        JetCorrections = cms.string(""),
0229                                        recoJetPtThreshold = cms.double(10),
0230                                        RThreshold = cms.double(0.3),
0231                                        reverseEnergyFractionThreshold = cms.double(0.5)
0232 )
0233 jetDQMAnalyzerAkPU4PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu4PFJets")
0234 jetDQMAnalyzerAkPU5PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu5PFJets")
0235
0236 jetDQMAnalyzerAkCs3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0237                                          JetType = cms.untracked.string('pf'),
0238                                          UEAlgo = cms.untracked.string('Cs'),
0239                                          OutputFile = cms.untracked.string(''),
0240                                          src = cms.InputTag("akCs3PFJets"),
0241                                          CScands = cms.InputTag("akCs3PFJets","pfParticlesCs"),
0242                                          PFcands = cms.InputTag("particleFlowTmp"),
0243                                          centralitycollection = cms.InputTag("hiCentrality"),
0244                                          #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0245                                          JetCorrections = cms.string(""),
0246                                          recoJetPtThreshold = cms.double(10),
0247                                          RThreshold = cms.double(0.3),
0248                                          reverseEnergyFractionThreshold = cms.double(0.5),
0249                                          etaMap    = cms.InputTag('hiFJRhoProducer','mapEtaEdges'),
0250                                          rho       = cms.InputTag('hiFJRhoProducer','mapToRho'),
0251                                          rhom      = cms.InputTag('hiFJRhoProducer','mapToRhoM')
0252 )
0253
0254 jetDQMAnalyzerAkCs4PF=jetDQMAnalyzerAkCs3PF.clone(src = "akCs4PFJets",
0255                           CScands = "akCs4PFJets:pfParticlesCs"
0256 )
0257
0258
0259 jetDQMMatchAkPu3CaloAkPu3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0260                                              src_Jet1 = cms.InputTag("akPu3CaloJets"),
0261                                              src_Jet2 = cms.InputTag("akPu3PFJets"),
0262                                              Jet1     = cms.untracked.string("PuCalo"),
0263                                              Jet2     = cms.untracked.string("PuPF"),
0264                                              recoJetPtThreshold = cms.double(20.),
0265                                              recoDelRMatch = cms.double(0.2),
0266                                              recoJetEtaCut = cms.double(2.0)
0267 )
0268
0269 jetDQMMatchAkPu4CaloAkPu4PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0270                                              src_Jet1 = cms.InputTag("akPu4CaloJets"),
0271                                              src_Jet2 = cms.InputTag("akPu4PFJets"),
0272                                              Jet1     = cms.untracked.string("PuCalo"),
0273                                              Jet2     = cms.untracked.string("PuPF"),
0274                                              recoJetPtThreshold = cms.double(20.),
0275                                              recoDelRMatch = cms.double(0.2),
0276                                              recoJetEtaCut = cms.double(2.0)
0277 )
0278
0279 jetDQMMatchAkPu5CaloAkPu5PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0280                                              src_Jet1 = cms.InputTag("akPu5CaloJets"),
0281                                              src_Jet2 = cms.InputTag("akPu5PFJets"),
0282                                              Jet1     = cms.untracked.string("PuCalo"),
0283                                              Jet2     = cms.untracked.string("PuPF"),
0284                                              recoJetPtThreshold = cms.double(20.),
0285                                              recoDelRMatch = cms.double(0.2),
0286                                              recoJetEtaCut = cms.double(2.0)
0287 )
0288