CMS 3D CMS Logo

jetAnalyzer_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 from DQMOffline.JetMET.jetDQMConfig_cff import * # parameters for all jet analyzers
4 from DQMOffline.JetMET.jetMETDQMCleanup_cff import * # parameters for event cleanup
5 
6 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
7 jetDQMAnalyzerAk4CaloUncleaned = DQMEDAnalyzer('JetAnalyzer',
8  JetType = cms.string('calo'),#pf, calo or jpt
9  JetCorrections = cms.InputTag("dqmAk4CaloL2L3ResidualCorrector"),
10  jetsrc = cms.InputTag("ak4CaloJets"),
11  METCollectionLabel = cms.InputTag("caloMet"),
12  muonsrc = cms.InputTag("muons"),
13  l1algoname = cms.string("L1Tech_BPTX_plus_AND_minus.v0"),
14  filljetHighLevel =cms.bool(False),
15  fillsubstructure =cms.bool(False),
16  ptMinBoosted = cms.double(400.),
17  #
18  #
19  #
20  highPtJetTrigger = cms.PSet(
21  andOr = cms.bool( False ),
22  dbLabel = cms.string("JetMETDQMTrigger"),
23  hltInputTag = cms.InputTag( "TriggerResults::HLT" ),
24  hltPaths = cms.vstring( 'HLT_PFJet450_v*'),
25  andOrHlt = cms.bool( True ),
26  errorReplyHlt = cms.bool( False ),
27  ),
28  lowPtJetTrigger = cms.PSet(
29  andOr = cms.bool( False ),
30  dbLabel = cms.string("JetMETDQMTrigger"),
31  hltInputTag = cms.InputTag( "TriggerResults::HLT" ),
32  hltPaths = cms.vstring( 'HLT_PFJet80_v*'),
33  andOrHlt = cms.bool( True ),
34  errorReplyHlt = cms.bool( False ),
35  ),
36 
37  TriggerResultsLabel = cms.InputTag("TriggerResults::HLT"),
38  processname = cms.string("HLT"),
39 
40  #
41  # Jet-related
42  #
43 
44  JetCleaningFlag = cms.untracked.bool(False),
45 
46  runcosmics = cms.untracked.bool(False),
47 
48  #Cleanup parameters
49  CleaningParameters = cleaningParameters.clone(
50  bypassAllPVChecks = cms.bool(True),
51  ),
52 
53  #for JPT and CaloJetID
54  InputJetIDValueMap = cms.InputTag("ak4JetID"),
55  #options for Calo and JPT: LOOSE,LOOSE_AOD,TIGHT,MINIMAL
56  #for PFJets: LOOSE,TIGHT
57  JetIDQuality = cms.string("LOOSE"),
58  #options for Calo and JPT: PURE09,DQM09,CRAFT08
59  #for PFJets: FIRSTDATA
60  JetIDVersion = cms.string("PURE09"),
61  #
62  #actually done only for PFJets at the moment
63  InputMVAPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorDQM","fullDiscriminant"),
64  InputCutPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorDQM","cutbasedDiscriminant"),
65  InputMVAPUIDValue = cms.InputTag("pileupJetIdEvaluatorDQM","fullId"),
66  InputCutPUIDValue = cms.InputTag("pileupJetIdEvaluatorDQM","cutbasedId"),
67 
68  InputQGMultiplicity = cms.InputTag("QGTagger", "mult"),
69  InputQGLikelihood = cms.InputTag("QGTagger", "qgLikelihood"),
70  InputQGPtDToken = cms.InputTag("QGTagger", "ptD"),
71  InputQGAxis2 = cms.InputTag("QGTagger", "axis2"),
72 
73  fillCHShistos =cms.bool(False),
74  #
75  # For jetAnalysis
76  #
77  jetAnalysis = jetDQMParameters.clone(),
78 
79  #
80  # DCS
81  #
82  DCSFilterForJetMonitoring = cms.PSet(
83  DetectorTypes = cms.untracked.string("ecal:hbhe:hf"),
84  #DebugOn = cms.untracked.bool(True),
85  alwaysPass = cms.untracked.bool(False)
86  )
87 )
88 
89 jetDQMAnalyzerAk4CaloCleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
90  JetCleaningFlag = cms.untracked.bool(True),
91  filljetHighLevel = cms.bool(False),
92  CleaningParameters = cleaningParameters.clone(
93  bypassAllPVChecks = cms.bool(True),
94  ),
95  jetAnalysis=jetDQMParameters.clone(
96  ptThreshold = cms.double(20.),
97  asymmetryThirdJetCut = cms.double(30),
98  balanceThirdJetCut = cms.double(0.2),
99  )
100 )
101 
102 jetDQMAnalyzerAk4PFUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
103  CleaningParameters = cleaningParameters.clone(
104  bypassAllPVChecks = cms.bool(False),
105  ),
106  #for PFJets: LOOSE,TIGHT
107  JetIDQuality = cms.string("LOOSE"),
108  #options for Calo and JPT: PURE09,DQM09,CRAFT08
109  #for PFJets: FIRSTDATA or RUNIISTARTUP (suitable for RECO beyond 7_2_X) or WINTER16 (for 8_0_X onwards)
110  JetIDVersion = cms.string("WINTER16"),
111  JetType = cms.string('pf'),#pf, calo or jpt
112  JetCorrections = cms.InputTag("dqmAk4PFL1FastL2L3ResidualCorrector"),
113  jetsrc = cms.InputTag("ak4PFJets"),
114  METCollectionLabel = cms.InputTag("pfMet"),
115  filljetHighLevel = cms.bool(False),
116  DCSFilterForJetMonitoring = cms.PSet(
117  DetectorTypes = cms.untracked.string("ecal:hbhe:hf:pixel:sistrip:es:muon"),
118  #DebugOn = cms.untracked.bool(True),
119  alwaysPass = cms.untracked.bool(False)
120  )
121 )
122 
123 
124 jetDQMAnalyzerAk4PFCleaned=jetDQMAnalyzerAk4PFUncleaned.clone(
125  JetCleaningFlag = cms.untracked.bool(True),
126  filljetHighLevel = cms.bool(False),
127  jetAnalysis=jetDQMParameters.clone(
128  ptThreshold = cms.double(20.),
129  asymmetryThirdJetCut = cms.double(30),
130  balanceThirdJetCut = cms.double(0.2),
131  ),
132  METCollectionLabel = cms.InputTag("pfMet"),
133 )
134 
135 jetDQMAnalyzerAk4PFCHSCleaned=jetDQMAnalyzerAk4PFCleaned.clone(
136  filljetHighLevel =cms.bool(True),
137  JetCorrections = cms.InputTag("dqmAk4PFCHSL1FastL2L3ResidualCorrector"),
138  jetsrc = cms.InputTag("ak4PFJetsCHS"),
139  METCollectionLabel = cms.InputTag("pfMETT1"),
140  InputMVAPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorCHSDQM","fullDiscriminant"),
141  InputCutPUIDDiscriminant = cms.InputTag("pileupJetIdEvaluatorCHSDQM","cutbasedDiscriminant"),
142  InputMVAPUIDValue = cms.InputTag("pileupJetIdEvaluatorCHSDQM","fullId"),
143  InputCutPUIDValue = cms.InputTag("pileupJetIdEvaluatorCHSDQM","cutbasedId"),
144  fillCHShistos =cms.bool(True),
145 )
146 
147 jetDQMAnalyzerAk4PFCHSUncleanedMiniAOD=jetDQMAnalyzerAk4PFUncleaned.clone(
148  filljetHighLevel =cms.bool(True),
149  CleaningParameters = cleaningParameters.clone(
150  vertexCollection = cms.InputTag( "goodOfflinePrimaryVerticesDQMforMiniAOD" ),
151  ),
152  JetType = cms.string('miniaod'),#pf, calo or jpt
153  jetsrc = cms.InputTag("slimmedJets"),
154  METCollectionLabel = cms.InputTag("slimmedMETs"),
155 )
156 
157 jetDQMAnalyzerAk4PFCHSCleanedMiniAOD=jetDQMAnalyzerAk4PFCleaned.clone(
158  CleaningParameters = cleaningParameters.clone(
159  vertexCollection = cms.InputTag( "goodOfflinePrimaryVerticesDQMforMiniAOD" ),
160  ),
161  JetType = cms.string('miniaod'),#pf, calo or jpt
162  jetsrc = cms.InputTag("slimmedJets"),
163 )
164 
165 jetDQMAnalyzerAk8PFPUPPICleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
166  jetsrc = cms.InputTag("slimmedJetsAK8"),
167  fillsubstructure =cms.bool(True),
168 )
169 
170 jetDQMAnalyzerAk4PFCHSPuppiCleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
171  JetType = cms.string('miniaod'),#pf, calo or jpt
172  jetsrc = cms.InputTag("slimmedJetsPuppi"),
173 )
174 
175 jetDQMAnalyzerIC5CaloHIUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
176  filljetHighLevel =cms.bool(True),
177  CleaningParameters = cleaningParameters.clone(
178  bypassAllPVChecks = cms.bool(False),
179  vertexCollection = cms.InputTag( "hiSelectedVertex" ),
180  ),
181  JetType = cms.string('calo'),#pf, calo or jpt
182  JetCorrections = cms.InputTag(""),# no jet correction available yet?
183  jetsrc = cms.InputTag("iterativeConePu5CaloJets"),
184  JetCleaningFlag = cms.untracked.bool(False),
185  runcosmics = cms.untracked.bool(True),
186  DCSFilterForJetMonitoring = cms.PSet(
187  DetectorTypes = cms.untracked.string("ecal:hbhe:hf:pixel:sistrip:es:muon"),
188  #DebugOn = cms.untracked.bool(True),
189  alwaysPass = cms.untracked.bool(False)
190  )
191 )
192 
193 jetDQMAnalyzerAkPU3Calo = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
194  JetType = cms.untracked.string('calo'),
195  UEAlgo = cms.untracked.string('Pu'),
196  OutputFile = cms.untracked.string(''),
197  src = cms.InputTag("akPu3CaloJets"),
198  PFcands = cms.InputTag("particleFlowTmp"),
199  Background = cms.InputTag("voronoiBackgroundCalo"),
200  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
201  centralitycollection = cms.InputTag("hiCentrality"),
202  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
203  JetCorrections = cms.string(""),
204  recoJetPtThreshold = cms.double(10),
205  RThreshold = cms.double(0.3),
206  reverseEnergyFractionThreshold = cms.double(0.5)
207 )
208 jetDQMAnalyzerAkPU4Calo = jetDQMAnalyzerAkPU3Calo.clone(src = cms.InputTag("akPu4CaloJets"))
209 jetDQMAnalyzerAkPU5Calo = jetDQMAnalyzerAkPU3Calo.clone(src = cms.InputTag("akPu5CaloJets"))
210 
211 jetDQMAnalyzerAkPU3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
212  JetType = cms.untracked.string('pf'),
213  UEAlgo = cms.untracked.string('Pu'),
214  OutputFile = cms.untracked.string(''),
215  src = cms.InputTag("akPu3PFJets"),
216  PFcands = cms.InputTag("particleFlowTmp"),
217  Background = cms.InputTag("voronoiBackgroundPF"),
218  centralitycollection = cms.InputTag("hiCentrality"),
219  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
220  JetCorrections = cms.string(""),
221  recoJetPtThreshold = cms.double(10),
222  RThreshold = cms.double(0.3),
223  reverseEnergyFractionThreshold = cms.double(0.5)
224 )
225 jetDQMAnalyzerAkPU4PF = jetDQMAnalyzerAkPU3PF.clone(src = cms.InputTag("akPu4PFJets"))
226 jetDQMAnalyzerAkPU5PF = jetDQMAnalyzerAkPU3PF.clone(src = cms.InputTag("akPu5PFJets"))
227 
228 jetDQMAnalyzerAkCs3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
229  JetType = cms.untracked.string('pf'),
230  UEAlgo = cms.untracked.string('Cs'),
231  OutputFile = cms.untracked.string(''),
232  src = cms.InputTag("akCs3PFJets"),
233  CScands = cms.InputTag("akCs3PFJets","pfParticlesCs"),
234  PFcands = cms.InputTag("particleFlowTmp"),
235  centralitycollection = cms.InputTag("hiCentrality"),
236  #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
237  JetCorrections = cms.string(""),
238  recoJetPtThreshold = cms.double(10),
239  RThreshold = cms.double(0.3),
240  reverseEnergyFractionThreshold = cms.double(0.5),
241  etaMap = cms.InputTag('hiFJRhoProducer','mapEtaEdges'),
242  rho = cms.InputTag('hiFJRhoProducer','mapToRho'),
243  rhom = cms.InputTag('hiFJRhoProducer','mapToRhoM')
244 )
245 
246 jetDQMAnalyzerAkCs4PF=jetDQMAnalyzerAkCs3PF.clone(src = cms.InputTag("akCs4PFJets"),
247  CScands = cms.InputTag("akCs4PFJets","pfParticlesCs")
248 )
249 
250 
251 jetDQMMatchAkPu3CaloAkPu3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
252  src_Jet1 = cms.InputTag("akPu3CaloJets"),
253  src_Jet2 = cms.InputTag("akPu3PFJets"),
254  Jet1 = cms.untracked.string("PuCalo"),
255  Jet2 = cms.untracked.string("PuPF"),
256  recoJetPtThreshold = cms.double(20.),
257  recoDelRMatch = cms.double(0.2),
258  recoJetEtaCut = cms.double(2.0)
259 )
260 
261 jetDQMMatchAkPu4CaloAkPu4PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
262  src_Jet1 = cms.InputTag("akPu4CaloJets"),
263  src_Jet2 = cms.InputTag("akPu4PFJets"),
264  Jet1 = cms.untracked.string("PuCalo"),
265  Jet2 = cms.untracked.string("PuPF"),
266  recoJetPtThreshold = cms.double(20.),
267  recoDelRMatch = cms.double(0.2),
268  recoJetEtaCut = cms.double(2.0)
269 )
270 
271 jetDQMMatchAkPu5CaloAkPu5PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
272  src_Jet1 = cms.InputTag("akPu5CaloJets"),
273  src_Jet2 = cms.InputTag("akPu5PFJets"),
274  Jet1 = cms.untracked.string("PuCalo"),
275  Jet2 = cms.untracked.string("PuPF"),
276  recoJetPtThreshold = cms.double(20.),
277  recoDelRMatch = cms.double(0.2),
278  recoJetEtaCut = cms.double(2.0)
279 )
280