JetMET/python/jetMETDQMOfflineSource_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from DQMOffline.JetMET.metDQMConfig_cff     import *
0004 from DQMOffline.JetMET.jetAnalyzer_cff   import *
0005 from DQMOffline.JetMET.pfCandidateDQMConfig_cfi   import *
0006 from DQMOffline.JetMET.SUSYDQMAnalyzer_cfi  import *
0007 from DQMOffline.JetMET.goodOfflinePrimaryVerticesDQM_cfi import *
0008 from RecoJets.JetProducers.PileupJetID_cfi  import *
0009 from RecoJets.JetProducers.QGTagger_cfi  import *
0010 from RecoMET.METFilters.metFilters_cff import*
0011
0012 HcalStripHaloFilterDQM=HcalStripHaloFilter.clone(taggingMode = True)
0013 CSCTightHaloFilterDQM=CSCTightHaloFilter.clone(taggingMode = True)
0014 CSCTightHalo2015FilterDQM=CSCTightHalo2015Filter.clone(taggingMode = True)
0015 eeBadScFilterDQM=eeBadScFilter.clone(taggingMode = True)
0016 EcalDeadCellTriggerPrimitiveFilterDQM=EcalDeadCellTriggerPrimitiveFilter.clone(taggingMode = True)
0017 EcalDeadCellBoundaryEnergyFilterDQM=EcalDeadCellBoundaryEnergyFilter.clone(taggingMode = True)
0018
0019 pileupJetIdCalculatorDQM=pileupJetIdCalculator.clone(
0020     jets = "ak4PFJets",
0021     jec = "AK4PF",
0022     applyJec = True,
0023     inputIsCorrected = False
0024 )
0025
0026 pileupJetIdEvaluatorDQM=pileupJetIdEvaluator.clone(
0027     jets = "ak4PFJets",
0028     jetids = "pileupJetIdCalculatorDQM",
0029     jec = "AK4PF",
0030     applyJec = True,
0031     inputIsCorrected = False
0032 )
0033
0034 pileupJetIdCalculatorCHSDQM=pileupJetIdCalculator.clone(
0035     applyJec = True,
0036     inputIsCorrected = False
0037 )
0038
0039 pileupJetIdEvaluatorCHSDQM=pileupJetIdEvaluator.clone(
0040     jetids = "pileupJetIdCalculatorCHSDQM",
0041     applyJec = True,
0042     inputIsCorrected = False
0043     )
0044
0045 pileupJetIdCalculatorPUPPIDQM=pileupJetIdCalculator.clone(
0046     jets = "ak4PFJetsPuppi",
0047     jec = "AK4PFPuppi",
0048     applyJec = True,
0049     inputIsCorrected = False
0050 )
0051
0052 pileupJetIdEvaluatorPUPPIDQM=pileupJetIdEvaluator.clone(
0053     jets = "ak4PFJetsPuppi",
0054     jetids = "pileupJetIdCalculatorPUPPIDQM",
0055     jec = "AK4PFPuppi",
0056     applyJec = True,
0057     inputIsCorrected = False
0058 )
0059
0060 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4CaloL2L3ResidualCorrectorChain,ak4CaloL2L3ResidualCorrector,ak4CaloResidualCorrector,ak4CaloL2L3Corrector,ak4CaloL3AbsoluteCorrector,ak4CaloL2RelativeCorrector
0061
0062 dqmAk4CaloL2L3ResidualCorrector = ak4CaloL2L3ResidualCorrector.clone()
0063 dqmAk4CaloL2L3ResidualCorrectorChain = cms.Sequence(
0064     #ak4CaloL2RelativeCorrector*ak4CaloL3AbsoluteCorrector*ak4CaloResidualCorrector*
0065     dqmAk4CaloL2L3ResidualCorrector
0066 )
0067
0068 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4PFL1FastL2L3ResidualCorrectorChain,ak4PFL1FastL2L3ResidualCorrector,ak4PFCHSL1FastL2L3Corrector,ak4PFResidualCorrector,ak4PFL3AbsoluteCorrector,ak4PFL2RelativeCorrector,ak4PFL1FastjetCorrector
0069
0070 dqmAk4PFCHSL1FastL2L3Corrector = ak4PFCHSL1FastL2L3Corrector.clone()
0071 dqmAk4PFCHSL1FastL2L3CorrectorChain = cms.Sequence(
0072     #ak4CaloL2RelativeCorrector*ak4CaloL3AbsoluteCorrector*ak4CaloResidualCorrector*
0073     dqmAk4PFCHSL1FastL2L3Corrector
0074 )
0075
0076 dqmAk4PFL1FastL2L3ResidualCorrector = ak4PFL1FastL2L3ResidualCorrector.clone()
0077 dqmAk4PFL1FastL2L3ResidualCorrectorChain = cms.Sequence(
0078     #ak4PFL1FastjetCorrector*ak4PFL2RelativeCorrector*ak4PFL3AbsoluteCorrector*ak4PFResidualCorrector*
0079     dqmAk4PFL1FastL2L3ResidualCorrector
0080 )
0081
0082
0083 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4PFCHSL1FastL2L3ResidualCorrectorChain,ak4PFCHSL1FastL2L3ResidualCorrector,ak4PFCHSL1FastL2L3Corrector,ak4PFCHSResidualCorrector,ak4PFCHSL3AbsoluteCorrector,ak4PFCHSL2RelativeCorrector,ak4PFCHSL1FastjetCorrector
0084
0085 dqmAk4PFCHSL1FastL2L3ResidualCorrector = ak4PFCHSL1FastL2L3ResidualCorrector.clone()
0086 dqmAk4PFCHSL1FastL2L3ResidualCorrectorChain = cms.Sequence(
0087     #ak4PFCHSL1FastjetCorrector*ak4PFCHSL2RelativeCorrector*ak4PFCHSL3AbsoluteCorrector*ak4PFCHSResidualCorrector
0088     dqmAk4PFCHSL1FastL2L3ResidualCorrector
0089 )
0090
0091 dqmAk4PFCHSL1FastL2L3Corrector = ak4PFCHSL1FastL2L3Corrector.clone()
0092 dqmAk4PFCHSL1FastL2L3CorrectorChain = cms.Sequence(
0093     dqmAk4PFCHSL1FastL2L3Corrector
0094 )
0095
0096 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4PFPuppiL1FastL2L3ResidualCorrectorChain,ak4PFPuppiL1FastL2L3ResidualCorrector,ak4PFPuppiL1FastL2L3Corrector,ak4PFPuppiResidualCorrector,ak4PFPuppiL3AbsoluteCorrector,ak4PFPuppiL2RelativeCorrector,ak4PFPuppiL1FastjetCorrector
0097
0098 dqmAk4PFPuppiL1FastL2L3ResidualCorrector = ak4PFPuppiL1FastL2L3ResidualCorrector.clone()
0099 dqmAk4PFPuppiL1FastL2L3ResidualCorrectorChain = cms.Sequence(
0100     dqmAk4PFPuppiL1FastL2L3ResidualCorrector
0101 )
0102
0103 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4PFL1FastL2L3ResidualCorrectorChain,ak4PFL1FastL2L3ResidualCorrector,ak4PFCHSL1FastL2L3Corrector,ak4PFResidualCorrector,ak4PFL3AbsoluteCorrector,ak4PFL2RelativeCorrector,ak4PFL1FastjetCorrector
0104
0105 ak4PFScoutingL1FastjetCorrector = ak4PFL1FastjetCorrector.clone(
0106     algorithm   = cms.string('AK4PFHLT'),
0107     srcRho = cms.InputTag("hltScoutingPFPacker","rho")
0108     )
0109
0110 ak4PFScoutingL2RelativeCorrector = ak4PFL2RelativeCorrector.clone(
0111     algorithm = cms.string('AK4PFHLT'),
0112     )
0113
0114 ak4PFScoutingL3AbsoluteCorrector = ak4PFL3AbsoluteCorrector.clone(
0115     algorithm = cms.string('AK4PFHLT'),
0116     )
0117
0118 ak4PFScoutingResidualCorrector  = ak4PFResidualCorrector.clone(
0119     algorithm = cms.string('AK4PFHLT'),
0120     )
0121
0122 ak4PFScoutingL1FastL2L3ResidualCorrector = cms.EDProducer(
0123     'ChainedJetCorrectorProducer',
0124     correctors = cms.VInputTag('ak4PFScoutingL1FastjetCorrector','ak4PFScoutingL2RelativeCorrector','ak4PFScoutingL3AbsoluteCorrector','ak4PFScoutingResidualCorrector')
0125     )
0126
0127 ak4PFScoutingL1FastL2L3ResidualCorrectorTask = cms.Task(
0128     ak4PFScoutingL1FastjetCorrector, ak4PFScoutingL2RelativeCorrector, ak4PFScoutingL3AbsoluteCorrector, ak4PFScoutingResidualCorrector, ak4PFScoutingL1FastL2L3ResidualCorrector
0129 )
0130 ak4PFScoutingL1FastL2L3ResidualCorrectorChain = cms.Sequence(ak4PFScoutingL1FastL2L3ResidualCorrectorTask)
0131
0132 dqmAk4PFScoutingL1FastL2L3ResidualCorrector = ak4PFScoutingL1FastL2L3ResidualCorrector.clone()
0133 dqmAk4PFScoutingL1FastL2L3ResidualCorrectorChain = cms.Sequence(
0134     dqmAk4PFScoutingL1FastL2L3ResidualCorrector
0135 )
0136
0137
0138 HBHENoiseFilterResultProducerDQM=HBHENoiseFilterResultProducer.clone()
0139
0140 jetPreDQMTask = cms.Task(ak4CaloL2RelativeCorrector,
0141                          ak4CaloL3AbsoluteCorrector,
0142                          ak4CaloResidualCorrector,
0143                          ak4PFL1FastjetCorrector,
0144                          ak4PFL2RelativeCorrector,
0145                          ak4PFL3AbsoluteCorrector,
0146                          ak4PFResidualCorrector,
0147                          ak4PFCHSL1FastjetCorrector,
0148                          ak4PFCHSL2RelativeCorrector,
0149                          ak4PFCHSL3AbsoluteCorrector,
0150                          ak4PFCHSResidualCorrector,
0151 )
0152
0153 _jetPreDQMTaskWithPUPPI = cms.Task(ak4CaloL2RelativeCorrector,
0154                          ak4CaloL3AbsoluteCorrector,
0155                          ak4CaloResidualCorrector,
0156                          ak4PFL1FastjetCorrector,
0157                          ak4PFL2RelativeCorrector,
0158                          ak4PFL3AbsoluteCorrector,
0159                          ak4PFResidualCorrector,
0160                          ak4PFCHSL1FastjetCorrector,
0161                          ak4PFCHSL2RelativeCorrector,
0162                          ak4PFCHSL3AbsoluteCorrector,
0163                          ak4PFCHSResidualCorrector,
0164                          ak4PFPuppiL1FastjetCorrector,
0165                          ak4PFPuppiL2RelativeCorrector,
0166                          ak4PFPuppiL3AbsoluteCorrector,
0167                          ak4PFPuppiResidualCorrector,
0168 )
0169
0170 jetPreDQMTaskScouting = cms.Task(ak4PFScoutingL1FastjetCorrector,
0171                          ak4PFScoutingL2RelativeCorrector,
0172                          ak4PFScoutingL3AbsoluteCorrector,
0173                          ak4PFScoutingResidualCorrector,
0174                          )
0175
0176 jetPreDQMSeq=cms.Sequence(jetPreDQMTask)
0177 _jetPreDQMSeqWithPUPPI=cms.Sequence(_jetPreDQMTaskWithPUPPI)
0178 jetPreDQMSeqScouting=cms.Sequence(jetPreDQMTaskScouting)
0179
0180 from JetMETCorrections.Type1MET.correctedMet_cff import pfMetT1
0181 from JetMETCorrections.Type1MET.correctionTermsPfMetType0PFCandidate_cff import *
0182 from JetMETCorrections.Type1MET.correctionTermsPfMetType1Type2_cff import corrPfMetType1
0183
0184 dqmCorrPfMetType1=corrPfMetType1.clone(jetCorrLabel = cms.InputTag('dqmAk4PFCHSL1FastL2L3Corrector'),
0185                                        jetCorrLabelRes = cms.InputTag('dqmAk4PFCHSL1FastL2L3ResidualCorrector')
0186                                        )
0187 pfMETT1=pfMetT1.clone(srcCorrections = (
0188         'dqmCorrPfMetType1:type1',
0189         ))
0190
0191 jetMETDQMOfflineSource = cms.Sequence(AnalyzeSUSYDQM*QGTagger*
0192                                       pileupJetIdCalculatorCHSDQM*pileupJetIdEvaluatorCHSDQM*
0193                                       pileupJetIdCalculatorDQM*pileupJetIdEvaluatorDQM*
0194                                       jetPreDQMSeq*
0195                                       dqmAk4CaloL2L3ResidualCorrectorChain*dqmAk4PFL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3CorrectorChain*
0196                                       cms.ignore(goodOfflinePrimaryVerticesDQM)*
0197                                       dqmCorrPfMetType1*pfMETT1*jetDQMAnalyzerSequence*HBHENoiseFilterResultProducer*
0198                                       cms.ignore(CSCTightHaloFilterDQM)*cms.ignore(CSCTightHalo2015FilterDQM)*cms.ignore(eeBadScFilterDQM)*cms.ignore(EcalDeadCellTriggerPrimitiveFilterDQM)*cms.ignore(EcalDeadCellBoundaryEnergyFilterDQM)*cms.ignore(HcalStripHaloFilterDQM)
0199                                       *METDQMAnalyzerSequence
0200                                       *pfCandidateDQMAnalyzer)
0201
0202 _jetMETDQMOfflineSourceWithPUPPI = cms.Sequence(AnalyzeSUSYDQM*QGTagger*
0203                                       pileupJetIdCalculatorCHSDQM*pileupJetIdEvaluatorCHSDQM*
0204                                       pileupJetIdCalculatorDQM*pileupJetIdEvaluatorDQM*
0205                                       pileupJetIdCalculatorPUPPIDQM*pileupJetIdEvaluatorPUPPIDQM*
0206                                       _jetPreDQMSeqWithPUPPI*
0207                                       dqmAk4CaloL2L3ResidualCorrectorChain*dqmAk4PFL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3CorrectorChain*dqmAk4PFPuppiL1FastL2L3ResidualCorrectorChain*
0208                                       cms.ignore(goodOfflinePrimaryVerticesDQM)*
0209                                       dqmCorrPfMetType1*pfMETT1*jetDQMAnalyzerSequence*HBHENoiseFilterResultProducer*
0210                                       cms.ignore(CSCTightHaloFilterDQM)*cms.ignore(CSCTightHalo2015FilterDQM)*cms.ignore(eeBadScFilterDQM)*cms.ignore(EcalDeadCellTriggerPrimitiveFilterDQM)*cms.ignore(EcalDeadCellBoundaryEnergyFilterDQM)*cms.ignore(HcalStripHaloFilterDQM)
0211                                       *METDQMAnalyzerSequence
0212                                       *pfCandidateDQMAnalyzer)
0213
0214 jetMETDQMOfflineSourceScouting = cms.Sequence(jetPreDQMSeqScouting*
0215                                       dqmAk4PFScoutingL1FastL2L3ResidualCorrectorChain*
0216                                       jetDQMAnalyzerSequenceScouting)
0217
0218 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0219 (~pp_on_AA).toReplaceWith(jetMETDQMOfflineSource, _jetMETDQMOfflineSourceWithPUPPI)
0220
0221 from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
0222
0223 jetMETDQMOfflineRedoProductsMiniAOD = cms.Sequence(cms.ignore(goodOfflinePrimaryVerticesDQMforMiniAOD))
0224
0225 jetMETDQMOfflineSourceMiniAOD = cms.Sequence(jetDQMAnalyzerSequenceMiniAOD*METDQMAnalyzerSequenceMiniAOD*packedCandidateDQMAnalyzerMiniAOD)