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 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4PFCHSL1FastL2L3ResidualCorrectorChain,ak4PFCHSL1FastL2L3ResidualCorrector,ak4PFCHSL1FastL2L3Corrector,ak4PFCHSResidualCorrector,ak4PFCHSL3AbsoluteCorrector,ak4PFCHSL2RelativeCorrector,ak4PFCHSL1FastjetCorrector
0083
0084 dqmAk4PFCHSL1FastL2L3ResidualCorrector = ak4PFCHSL1FastL2L3ResidualCorrector.clone()
0085 dqmAk4PFCHSL1FastL2L3ResidualCorrectorChain = cms.Sequence(
0086     #ak4PFCHSL1FastjetCorrector*ak4PFCHSL2RelativeCorrector*ak4PFCHSL3AbsoluteCorrector*ak4PFCHSResidualCorrector
0087     dqmAk4PFCHSL1FastL2L3ResidualCorrector
0088 )
0089
0090 dqmAk4PFCHSL1FastL2L3Corrector = ak4PFCHSL1FastL2L3Corrector.clone()
0091 dqmAk4PFCHSL1FastL2L3CorrectorChain = cms.Sequence(
0092     dqmAk4PFCHSL1FastL2L3Corrector
0093 )
0094
0095 from JetMETCorrections.Configuration.JetCorrectors_cff import ak4PFPuppiL1FastL2L3ResidualCorrectorChain,ak4PFPuppiL1FastL2L3ResidualCorrector,ak4PFPuppiL1FastL2L3Corrector,ak4PFPuppiResidualCorrector,ak4PFPuppiL3AbsoluteCorrector,ak4PFPuppiL2RelativeCorrector,ak4PFPuppiL1FastjetCorrector
0096
0097 dqmAk4PFPuppiL1FastL2L3ResidualCorrector = ak4PFPuppiL1FastL2L3ResidualCorrector.clone()
0098 dqmAk4PFPuppiL1FastL2L3ResidualCorrectorChain = cms.Sequence(
0099     dqmAk4PFPuppiL1FastL2L3ResidualCorrector
0100 )
0101
0102 HBHENoiseFilterResultProducerDQM=HBHENoiseFilterResultProducer.clone()
0103
0104 jetPreDQMTask = cms.Task(ak4CaloL2RelativeCorrector,
0105                          ak4CaloL3AbsoluteCorrector,
0106                          ak4CaloResidualCorrector,
0107                          ak4PFL1FastjetCorrector,
0108                          ak4PFL2RelativeCorrector,
0109                          ak4PFL3AbsoluteCorrector,
0110                          ak4PFResidualCorrector,
0111                          ak4PFCHSL1FastjetCorrector,
0112                          ak4PFCHSL2RelativeCorrector,
0113                          ak4PFCHSL3AbsoluteCorrector,
0114                          ak4PFCHSResidualCorrector,
0115 )
0116
0117 _jetPreDQMTaskWithPUPPI = cms.Task(ak4CaloL2RelativeCorrector,
0118                          ak4CaloL3AbsoluteCorrector,
0119                          ak4CaloResidualCorrector,
0120                          ak4PFL1FastjetCorrector,
0121                          ak4PFL2RelativeCorrector,
0122                          ak4PFL3AbsoluteCorrector,
0123                          ak4PFResidualCorrector,
0124                          ak4PFCHSL1FastjetCorrector,
0125                          ak4PFCHSL2RelativeCorrector,
0126                          ak4PFCHSL3AbsoluteCorrector,
0127                          ak4PFCHSResidualCorrector,
0128                          ak4PFPuppiL1FastjetCorrector,
0129                          ak4PFPuppiL2RelativeCorrector,
0130                          ak4PFPuppiL3AbsoluteCorrector,
0131                          ak4PFPuppiResidualCorrector,
0132 )
0133 jetPreDQMSeq=cms.Sequence(jetPreDQMTask)
0134 _jetPreDQMSeqWithPUPPI=cms.Sequence(_jetPreDQMTaskWithPUPPI)
0135
0136 from JetMETCorrections.Type1MET.correctedMet_cff import pfMetT1
0137 from JetMETCorrections.Type1MET.correctionTermsPfMetType0PFCandidate_cff import *
0138 from JetMETCorrections.Type1MET.correctionTermsPfMetType1Type2_cff import corrPfMetType1
0139
0140 dqmCorrPfMetType1=corrPfMetType1.clone(jetCorrLabel = cms.InputTag('dqmAk4PFCHSL1FastL2L3Corrector'),
0141                                        jetCorrLabelRes = cms.InputTag('dqmAk4PFCHSL1FastL2L3ResidualCorrector')
0142                                        )
0143 pfMETT1=pfMetT1.clone(srcCorrections = (
0144         'dqmCorrPfMetType1:type1',
0145         ))
0146
0147 jetMETDQMOfflineSource = cms.Sequence(AnalyzeSUSYDQM*QGTagger*
0148                                       pileupJetIdCalculatorCHSDQM*pileupJetIdEvaluatorCHSDQM*
0149                                       pileupJetIdCalculatorDQM*pileupJetIdEvaluatorDQM*
0150                                       jetPreDQMSeq*
0151                                       dqmAk4CaloL2L3ResidualCorrectorChain*dqmAk4PFL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3CorrectorChain*
0152                                       cms.ignore(goodOfflinePrimaryVerticesDQM)*
0153                                       dqmCorrPfMetType1*pfMETT1*jetDQMAnalyzerSequence*HBHENoiseFilterResultProducer*
0154                                       cms.ignore(CSCTightHaloFilterDQM)*cms.ignore(CSCTightHalo2015FilterDQM)*cms.ignore(eeBadScFilterDQM)*cms.ignore(EcalDeadCellTriggerPrimitiveFilterDQM)*cms.ignore(EcalDeadCellBoundaryEnergyFilterDQM)*cms.ignore(HcalStripHaloFilterDQM)
0155                                       *METDQMAnalyzerSequence
0156                                       *pfCandidateDQMAnalyzer)
0157
0158 _jetMETDQMOfflineSourceWithPUPPI = cms.Sequence(AnalyzeSUSYDQM*QGTagger*
0159                                       pileupJetIdCalculatorCHSDQM*pileupJetIdEvaluatorCHSDQM*
0160                                       pileupJetIdCalculatorDQM*pileupJetIdEvaluatorDQM*
0161                                       pileupJetIdCalculatorPUPPIDQM*pileupJetIdEvaluatorPUPPIDQM*
0162                                       _jetPreDQMSeqWithPUPPI*
0163                                       dqmAk4CaloL2L3ResidualCorrectorChain*dqmAk4PFL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3ResidualCorrectorChain*dqmAk4PFCHSL1FastL2L3CorrectorChain*dqmAk4PFPuppiL1FastL2L3ResidualCorrectorChain*
0164                                       cms.ignore(goodOfflinePrimaryVerticesDQM)*
0165                                       dqmCorrPfMetType1*pfMETT1*jetDQMAnalyzerSequence*HBHENoiseFilterResultProducer*
0166                                       cms.ignore(CSCTightHaloFilterDQM)*cms.ignore(CSCTightHalo2015FilterDQM)*cms.ignore(eeBadScFilterDQM)*cms.ignore(EcalDeadCellTriggerPrimitiveFilterDQM)*cms.ignore(EcalDeadCellBoundaryEnergyFilterDQM)*cms.ignore(HcalStripHaloFilterDQM)
0167                                       *METDQMAnalyzerSequence
0168                                       *pfCandidateDQMAnalyzer)
0169
0170 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0171 (~pp_on_AA).toReplaceWith(jetMETDQMOfflineSource, _jetMETDQMOfflineSourceWithPUPPI)
0172
0173 from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
0174
0175 jetMETDQMOfflineRedoProductsMiniAOD = cms.Sequence(cms.ignore(goodOfflinePrimaryVerticesDQMforMiniAOD))
0176
0177 jetMETDQMOfflineSourceMiniAOD = cms.Sequence(jetDQMAnalyzerSequenceMiniAOD*METDQMAnalyzerSequenceMiniAOD*packedCandidateDQMAnalyzerMiniAOD)