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 jetDQMAnalizerAk4PUPPICleaned=jetDQMAnalyzerAk4PFCleaned.clone(
0148     JetType = cms.string('puppi'),
0149     jetsrc = "ak4PFJetsPuppi",
0150     METCollectionLabel = "pfMetPuppi",
0151     JetCorrections = "ak4PFPuppiL1FastL2L3ResidualCorrector",
0152     JetIDVersion = "RUN2ULPUPPI",
0153     JetIDQuality = cms.string("TIGHT"),
0154     fillCHShistos = True,
0155 )
0156
0157 jetDQMAnalyzerAk4PFCHSUncleanedMiniAOD=jetDQMAnalyzerAk4PFUncleaned.clone(
0158     filljetHighLevel = True,
0159     CleaningParameters = cleaningParameters.clone(
0160         vertexCollection  =  "goodOfflinePrimaryVerticesDQMforMiniAOD" ,
0161         ),
0162     JetType = 'miniaod',#pf, calo or jpt
0163     jetsrc = "slimmedJets",
0164     METCollectionLabel     = "slimmedMETs"
0165 )
0166
0167 jetDQMAnalyzerAk4PFCHSCleanedMiniAOD=jetDQMAnalyzerAk4PFCleaned.clone(
0168     CleaningParameters = cleaningParameters.clone(
0169         vertexCollection    =  "goodOfflinePrimaryVerticesDQMforMiniAOD"
0170         ),
0171     JetType = 'miniaod',#pf, calo or jpt
0172     jetsrc = "slimmedJets"
0173 )
0174
0175 jetDQMAnalyzerAk8PFPUPPICleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
0176     jetsrc = "slimmedJetsAK8",
0177     #for PUPPI jets: TIGHT
0178     JetIDQuality  = "TIGHT",
0179     #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
0180     JetIDVersion  = "RUN2ULPUPPI",
0181     fillsubstructure =True
0182 )
0183
0184 jetDQMAnalyzerAk4PFCHSPuppiCleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
0185     JetType = 'miniaod',#pf, calo or jpt
0186     jetsrc = "slimmedJetsPuppi",
0187     #for PUPPI jets: TIGHT
0188     JetIDQuality  = "TIGHT",
0189     #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
0190     JetIDVersion  = "RUN2ULPUPPI"
0191 )
0192
0193 jetDQMAnalyzerIC5CaloHIUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
0194     filljetHighLevel = True,
0195     CleaningParameters = cleaningParameters.clone(
0196         bypassAllPVChecks  = False,
0197         vertexCollection =  "hiSelectedVertex",
0198         ),
0199     JetType = 'calo',#pf, calo or jpt
0200     JetCorrections = "",# no jet correction available yet?
0201     jetsrc = "iterativeConePu5CaloJets",
0202     JetCleaningFlag            = False,
0203     runcosmics                 = True,
0204     DCSFilterForJetMonitoring = dict(
0205         DetectorTypes = "ecal:hbhe:hf:pixel:sistrip:es:muon",
0206         #DebugOn = True,
0207         alwaysPass = False
0208     )
0209 )
0210
0211 jetDQMAnalyzerAkPU3Calo = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0212                                          JetType = cms.untracked.string('calo'),
0213                                          UEAlgo = cms.untracked.string('Pu'),
0214                                          OutputFile = cms.untracked.string(''),
0215                                          src = cms.InputTag("akPu3CaloJets"),
0216                                          PFcands = cms.InputTag("particleFlowTmp"),
0217                                          Background = cms.InputTag("voronoiBackgroundCalo"),
0218                                          #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
0219                                          centralitycollection = cms.InputTag("hiCentrality"),
0220                                          centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0221                                          JetCorrections = cms.string(""),
0222                                          recoJetPtThreshold = cms.double(10),
0223                                          RThreshold = cms.double(0.3),
0224                                          reverseEnergyFractionThreshold = cms.double(0.5)
0225 )
0226 jetDQMAnalyzerAkPU4Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu4CaloJets")
0227 jetDQMAnalyzerAkPU5Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu5CaloJets")
0228
0229 jetDQMAnalyzerAkPU3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0230                                        JetType = cms.untracked.string('pf'),
0231                                        UEAlgo = cms.untracked.string('Pu'),
0232                                        OutputFile = cms.untracked.string(''),
0233                                        src = cms.InputTag("akPu3PFJets"),
0234                                        PFcands = cms.InputTag("particleFlowTmp"),
0235                                        Background = cms.InputTag("voronoiBackgroundPF"),
0236                                        centralitycollection = cms.InputTag("hiCentrality"),
0237                                        centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0238                                        JetCorrections = cms.string(""),
0239                                        recoJetPtThreshold = cms.double(10),
0240                                        RThreshold = cms.double(0.3),
0241                                        reverseEnergyFractionThreshold = cms.double(0.5)
0242 )
0243 jetDQMAnalyzerAkPU4PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu4PFJets")
0244 jetDQMAnalyzerAkPU5PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu5PFJets")
0245
0246 jetDQMAnalyzerAkCs3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
0247                                          JetType = cms.untracked.string('pf'),
0248                                          UEAlgo = cms.untracked.string('Cs'),
0249                                          OutputFile = cms.untracked.string(''),
0250                                          src = cms.InputTag("akCs3PFJets"),
0251                                          CScands = cms.InputTag("akCs3PFJets","pfParticlesCs"),
0252                                          PFcands = cms.InputTag("particleFlowTmp"),
0253                                          centralitycollection = cms.InputTag("hiCentrality"),
0254                                          #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
0255                                          JetCorrections = cms.string(""),
0256                                          recoJetPtThreshold = cms.double(10),
0257                                          RThreshold = cms.double(0.3),
0258                                          reverseEnergyFractionThreshold = cms.double(0.5),
0259                                          etaMap    = cms.InputTag('hiFJRhoProducer','mapEtaEdges'),
0260                                          rho       = cms.InputTag('hiFJRhoProducer','mapToRho'),
0261                                          rhom      = cms.InputTag('hiFJRhoProducer','mapToRhoM')
0262 )
0263
0264 jetDQMAnalyzerAkCs4PF=jetDQMAnalyzerAkCs3PF.clone(src = "akCs4PFJets",
0265                           CScands = "akCs4PFJets:pfParticlesCs"
0266 )
0267
0268
0269 jetDQMMatchAkPu3CaloAkPu3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0270                                              src_Jet1 = cms.InputTag("akPu3CaloJets"),
0271                                              src_Jet2 = cms.InputTag("akPu3PFJets"),
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 jetDQMMatchAkPu4CaloAkPu4PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0280                                              src_Jet1 = cms.InputTag("akPu4CaloJets"),
0281                                              src_Jet2 = cms.InputTag("akPu4PFJets"),
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
0289 jetDQMMatchAkPu5CaloAkPu5PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
0290                                              src_Jet1 = cms.InputTag("akPu5CaloJets"),
0291                                              src_Jet2 = cms.InputTag("akPu5PFJets"),
0292                                              Jet1     = cms.untracked.string("PuCalo"),
0293                                              Jet2     = cms.untracked.string("PuPF"),
0294                                              recoJetPtThreshold = cms.double(20.),
0295                                              recoDelRMatch = cms.double(0.2),
0296                                              recoJetEtaCut = cms.double(2.0)
0297 )
0298