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     #Only for JPT and CaloJetID  -> need to define InputJetIDValueMap
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 jetDQMAnalizerAk4PUPPICleaned=jetDQMAnalyzerAk4PFCleaned.clone(
0148     JetType = cms.string('puppi'),
0149     jetsrc = "ak4PFJetsPuppi",
0150     METCollectionLabel = "pfMetPuppi",
0151     JetCorrections = cms.InputTag("dqmAk4PFPuppiL1FastL2L3ResidualCorrector"),
0152     JetIDVersion = "RUN2ULPUPPI",
0153     JetIDQuality = cms.string("TIGHT"),
0154     fillCHShistos = True,
0155 )
0156
0157 jetDQMAnalyzerAk4ScoutingUncleaned = jetDQMAnalyzerAk4CaloUncleaned.clone(
0158     JetType = cms.string('scouting'),
0159     JetCorrections = cms.InputTag("dqmAk4PFScoutingL1FastL2L3ResidualCorrector"),
0160     jetsrc = cms.InputTag("hltScoutingPFPacker"),                       ###---> this name goes to GUI and TTree under .../JetMET/Run summary/Jet
0161     srcRho = cms.InputTag("hltScoutingPFPacker","rho"),
0162     METCollectionLabel = cms.InputTag("hltScoutingPFPacker","pfMetPt"), ###---> this name goes to GUI and TTree under .../JetMET/Run summary/MET
0163     muonsrc = cms.InputTag("hltScoutingMuonPackerNoVtx","","HLT"),
0164     l1algoname = cms.string("L1Tech_BPTX_plus_AND_minus.v0"),
0165     filljetHighLevel = False,                                           ### for plots: "vertices" and "cleanup" in .../JetMET/Run summary/
0166     fillsubstructure = False,
0167
0168     highPtJetTrigger = cms.PSet(
0169         andOr         = cms.bool( False ),
0170         dbLabel        = cms.string("JetMETDQMTrigger"),
0171         hltInputTag    = cms.InputTag( "TriggerResults::HLT" ),
0172         hltPaths       = cms.vstring( 'DST_PFScouting_JetHT_v*'),
0173         andOrHlt       = cms.bool( True ),
0174         errorReplyHlt  = cms.bool( False ),
0175     ),
0176     lowPtJetTrigger = cms.PSet(
0177         andOr         = cms.bool( False ),
0178         dbLabel        = cms.string("JetMETDQMTrigger"),
0179         hltInputTag    = cms.InputTag( "TriggerResults::HLT" ),
0180         hltPaths       = cms.vstring( 'DST_PFScouting_ZeroBias_v*'),
0181         andOrHlt       = cms.bool( True ),
0182         errorReplyHlt  = cms.bool( False ),
0183     ),
0184
0185     TriggerResultsLabel        = cms.InputTag("TriggerResults::HLT"),
0186     processname                = cms.string("HLT"),
0187
0188     #
0189     # Jet-related
0190     #
0191
0192     JetCleaningFlag            = cms.untracked.bool(False),
0193
0194     runcosmics                 = False,
0195
0196     #Cleanup parameters
0197     CleaningParameters = cleaningParameters.clone(
0198         bypassAllPVChecks = False,
0199         ),
0200
0201     #for scouting PF jets: TIGHT
0202     JetIDQuality               = cms.string("TIGHT"),
0203     #for scouting PF jets: RUN3Scouting
0204     JetIDVersion               = cms.string("RUN3Scouting"),
0205     #
0206     #Pileup JetID anf quark-gluon discrimination not exist for scouting PF jets. The following: actually done only for PFJets at the moment
0207     ###InputMVAPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorDQM","fullDiscriminant"),
0208     ###InputCutPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorDQM","cutbasedDiscriminant"),
0209     ###InputMVAPUIDValue = cms.InputTag("pileupJetIdEvaluatorDQM","fullId"),
0210     ###InputCutPUIDValue = cms.InputTag("pileupJetIdEvaluatorDQM","cutbasedId"),
0211
0212     ###InputQGMultiplicity = cms.InputTag("QGTagger", "mult"),
0213     ###InputQGLikelihood = cms.InputTag("QGTagger", "qgLikelihood"),
0214     ###InputQGPtDToken = cms.InputTag("QGTagger", "ptD"),
0215     ###InputQGAxis2 = cms.InputTag("QGTagger", "axis2"),
0216
0217     fillCHShistos = False,
0218     #
0219     # For jetAnalysis
0220     #
0221     jetAnalysis = jetDQMParameters.clone(),
0222
0223     #
0224     # DCS ### -> only used in JetMETDQMFilter.cc
0225     #
0226     DCSFilterForJetMonitoring = cms.PSet(
0227       DetectorTypes = cms.untracked.string("ecal:hbhe:hf:pixel:sistrip:es:muon"),
0228       #DebugOn = cms.untracked.bool(True),
0229       alwaysPass = cms.untracked.bool(False)
0230     )
0231 )
0232
0233 jetDQMAnalyzerAk4ScoutingCleaned = jetDQMAnalyzerAk4ScoutingUncleaned.clone(
0234     JetCleaningFlag = True
0235 )
0236
0237 jetDQMAnalyzerAk4PFCHSUncleanedMiniAOD=jetDQMAnalyzerAk4PFUncleaned.clone(
0238     filljetHighLevel = True,
0239     CleaningParameters = cleaningParameters.clone(
0240         vertexCollection  =  "goodOfflinePrimaryVerticesDQMforMiniAOD" ,
0241         ),
0242     JetType = 'miniaod',#pf, calo or jpt
0243     jetsrc = "slimmedJets",
0244     METCollectionLabel     = "slimmedMETs"
0245 )
0246
0247 jetDQMAnalyzerAk4PFCHSCleanedMiniAOD=jetDQMAnalyzerAk4PFCleaned.clone(
0248     CleaningParameters = cleaningParameters.clone(
0249         vertexCollection    =  "goodOfflinePrimaryVerticesDQMforMiniAOD"
0250         ),
0251     JetType = 'miniaod',#pf, calo or jpt
0252     jetsrc = "slimmedJets"
0253 )
0254
0255 jetDQMAnalyzerAk8PFPUPPICleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
0256     jetsrc = "slimmedJetsAK8",
0257     #for PUPPI jets: TIGHT
0258     JetIDQuality  = "TIGHT",
0259     #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
0260     JetIDVersion  = "RUN2ULPUPPI",
0261     fillsubstructure =True
0262 )
0263
0264 jetDQMAnalyzerAk4PFCHSPuppiCleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
0265     JetType = 'miniaod',#pf, calo or jpt
0266     jetsrc = "slimmedJetsPuppi",
0267     #for PUPPI jets: TIGHT
0268     JetIDQuality  = "TIGHT",
0269     #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
0270     JetIDVersion  = "RUN2ULPUPPI"
0271 )
0272
0273 jetDQMAnalyzerIC5CaloHIUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
0274     filljetHighLevel = True,
0275     CleaningParameters = cleaningParameters.clone(
0276         bypassAllPVChecks  = False,
0277         vertexCollection =  "hiSelectedVertex",
0278         ),
0279     JetType = 'calo',#pf, calo or jpt
0280     JetCorrections = "",# no jet correction available yet?
0281     jetsrc = "iterativeConePu5CaloJets",
0282     JetCleaningFlag            = False,
0283     runcosmics                 = True,
0284     DCSFilterForJetMonitoring = dict(
0285         DetectorTypes = "ecal:hbhe:hf:pixel:sistrip:es:muon",
0286         #DebugOn = True,
0287         alwaysPass = False
0288     )
0289 )
0290
0291 jetDQMAnalyzerAkPU3Calo = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0292                                          JetType = cms.untracked.string('calo'),
0293                                          UEAlgo = cms.untracked.string('Pu'),
0294                                          OutputFile = cms.untracked.string(''),
0295                                          src = cms.InputTag("akPu3CaloJets"),
0296                                          PFcands = cms.InputTag("particleFlowTmp"),
0297                                          Background = cms.InputTag("voronoiBackgroundCalo"),
0298                                          #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0299                                          centralitycollection = cms.InputTag("hiCentrality"),
0300                                          centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0301                                          JetCorrections = cms.string(""),
0302                                          recoJetPtThreshold = cms.double(10),
0303                                          RThreshold = cms.double(0.3),
0304                                          reverseEnergyFractionThreshold = cms.double(0.5)
0305 )
0306 jetDQMAnalyzerAkPU4Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu4CaloJets")
0307 jetDQMAnalyzerAkPU5Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu5CaloJets")
0308
0309 jetDQMAnalyzerAkPU3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0310                                        JetType = cms.untracked.string('pf'),
0311                                        UEAlgo = cms.untracked.string('Pu'),
0312                                        OutputFile = cms.untracked.string(''),
0313                                        src = cms.InputTag("akPu3PFJets"),
0314                                        PFcands = cms.InputTag("particleFlowTmp"),
0315                                        Background = cms.InputTag("voronoiBackgroundPF"),
0316                                        centralitycollection = cms.InputTag("hiCentrality"),
0317                                        centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0318                                        JetCorrections = cms.string(""),
0319                                        recoJetPtThreshold = cms.double(10),
0320                                        RThreshold = cms.double(0.3),
0321                                        reverseEnergyFractionThreshold = cms.double(0.5)
0322 )
0323 jetDQMAnalyzerAkPU4PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu4PFJets")
0324 jetDQMAnalyzerAkPU5PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu5PFJets")
0325
0326 jetDQMAnalyzerAkCs3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0327                                          JetType = cms.untracked.string('pf'),
0328                                          UEAlgo = cms.untracked.string('Cs'),
0329                                          OutputFile = cms.untracked.string(''),
0330                                          src = cms.InputTag("akCs3PFJets"),
0331                                          CScands = cms.InputTag("akCs3PFJets","pfParticlesCs"),
0332                                          PFcands = cms.InputTag("particleFlowTmp"),
0333                                          centralitycollection = cms.InputTag("hiCentrality"),
0334                                          #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0335                                          JetCorrections = cms.string(""),
0336                                          recoJetPtThreshold = cms.double(10),
0337                                          RThreshold = cms.double(0.3),
0338                                          reverseEnergyFractionThreshold = cms.double(0.5),
0339                                          etaMap    = cms.InputTag('hiFJRhoProducer','mapEtaEdges'),
0340                                          rho       = cms.InputTag('hiFJRhoProducer','mapToRho'),
0341                                          rhom      = cms.InputTag('hiFJRhoProducer','mapToRhoM')
0342 )
0343
0344 jetDQMAnalyzerAkCs4PF=jetDQMAnalyzerAkCs3PF.clone(src = "akCs4PFJets",
0345                           CScands = "akCs4PFJets:pfParticlesCs"
0346 )
0347
0348
0349 jetDQMMatchAkPu3CaloAkPu3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0350                                              src_Jet1 = cms.InputTag("akPu3CaloJets"),
0351                                              src_Jet2 = cms.InputTag("akPu3PFJets"),
0352                                              Jet1     = cms.untracked.string("PuCalo"),
0353                                              Jet2     = cms.untracked.string("PuPF"),
0354                                              recoJetPtThreshold = cms.double(20.),
0355                                              recoDelRMatch = cms.double(0.2),
0356                                              recoJetEtaCut = cms.double(2.0)
0357 )
0358
0359 jetDQMMatchAkPu4CaloAkPu4PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0360                                              src_Jet1 = cms.InputTag("akPu4CaloJets"),
0361                                              src_Jet2 = cms.InputTag("akPu4PFJets"),
0362                                              Jet1     = cms.untracked.string("PuCalo"),
0363                                              Jet2     = cms.untracked.string("PuPF"),
0364                                              recoJetPtThreshold = cms.double(20.),
0365                                              recoDelRMatch = cms.double(0.2),
0366                                              recoJetEtaCut = cms.double(2.0)
0367 )
0368
0369 jetDQMMatchAkPu5CaloAkPu5PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0370                                              src_Jet1 = cms.InputTag("akPu5CaloJets"),
0371                                              src_Jet2 = cms.InputTag("akPu5PFJets"),
0372                                              Jet1     = cms.untracked.string("PuCalo"),
0373                                              Jet2     = cms.untracked.string("PuPF"),
0374                                              recoJetPtThreshold = cms.double(20.),
0375                                              recoDelRMatch = cms.double(0.2),
0376                                              recoJetEtaCut = cms.double(2.0)
0377 )
0378