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 = 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 = True,
91  filljetHighLevel = False,
92  CleaningParameters = cleaningParameters.clone(
93  bypassAllPVChecks = True,
94  ),
95  jetAnalysis=jetDQMParameters.clone(
96  ptThreshold = 20.,
97  asymmetryThirdJetCut = 30,
98  balanceThirdJetCut = 0.2,
99  )
100 )
101 
102 jetDQMAnalyzerAk4PFUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
103  CleaningParameters = cleaningParameters.clone(
104  bypassAllPVChecks = False,
105  ),
106  #for PFJets: LOOSE,TIGHT
107  JetIDQuality = cms.string("TIGHT"),
108  #options for Calo and JPT: PURE09,DQM09,CRAFT08
109  #for PFJets: RUN2ULCHS for 11_1_X onwards
110  JetIDVersion = cms.string("RUN2ULCHS"),
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 = True,
126  filljetHighLevel = False,
127  jetAnalysis=jetDQMParameters.clone(
128  ptThreshold = 20.,
129  asymmetryThirdJetCut = 30,
130  balanceThirdJetCut = 0.2,
131  ),
132  METCollectionLabel = "pfMet"
133 )
134 
135 jetDQMAnalyzerAk4PFCHSCleaned=jetDQMAnalyzerAk4PFCleaned.clone(
136  filljetHighLevel = True,
137  JetCorrections = "dqmAk4PFCHSL1FastL2L3ResidualCorrector",
138  jetsrc = "ak4PFJetsCHS",
139  METCollectionLabel = "pfMETT1",
140  InputMVAPUIDDiscriminant = "pileupJetIdEvaluatorCHSDQM:fullDiscriminant",
141  InputCutPUIDDiscriminant = "pileupJetIdEvaluatorCHSDQM:cutbasedDiscriminant",
142  InputMVAPUIDValue = "pileupJetIdEvaluatorCHSDQM:fullId",
143  InputCutPUIDValue = "pileupJetIdEvaluatorCHSDQM:cutbasedId",
144  fillCHShistos = True
145 )
146 
147 jetDQMAnalyzerAk4PFCHSUncleanedMiniAOD=jetDQMAnalyzerAk4PFUncleaned.clone(
148  filljetHighLevel = True,
149  CleaningParameters = cleaningParameters.clone(
150  vertexCollection = "goodOfflinePrimaryVerticesDQMforMiniAOD" ,
151  ),
152  JetType = 'miniaod',#pf, calo or jpt
153  jetsrc = "slimmedJets",
154  METCollectionLabel = "slimmedMETs"
155 )
156 
157 jetDQMAnalyzerAk4PFCHSCleanedMiniAOD=jetDQMAnalyzerAk4PFCleaned.clone(
158  CleaningParameters = cleaningParameters.clone(
159  vertexCollection = "goodOfflinePrimaryVerticesDQMforMiniAOD"
160  ),
161  JetType = 'miniaod',#pf, calo or jpt
162  jetsrc = "slimmedJets"
163 )
164 
165 jetDQMAnalyzerAk8PFPUPPICleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
166  jetsrc = "slimmedJetsAK8",
167  #for PUPPI jets: TIGHT
168  JetIDQuality = "TIGHT",
169  #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
170  JetIDVersion = "RUN2ULPUPPI",
171  fillsubstructure =True
172 )
173 
174 jetDQMAnalyzerAk4PFCHSPuppiCleanedMiniAOD=jetDQMAnalyzerAk4PFCHSCleanedMiniAOD.clone(
175  JetType = 'miniaod',#pf, calo or jpt
176  jetsrc = "slimmedJetsPuppi",
177  #for PUPPI jets: TIGHT
178  JetIDQuality = "TIGHT",
179  #for PUPPI jets: RUN2ULPUPPI from 11_1_X onwards
180  JetIDVersion = "RUN2ULPUPPI"
181 )
182 
183 jetDQMAnalyzerIC5CaloHIUncleaned=jetDQMAnalyzerAk4CaloUncleaned.clone(
184  filljetHighLevel = True,
185  CleaningParameters = cleaningParameters.clone(
186  bypassAllPVChecks = False,
187  vertexCollection = "hiSelectedVertex",
188  ),
189  JetType = 'calo',#pf, calo or jpt
190  JetCorrections = "",# no jet correction available yet?
191  jetsrc = "iterativeConePu5CaloJets",
192  JetCleaningFlag = False,
193  runcosmics = True,
194  DCSFilterForJetMonitoring = dict(
195  DetectorTypes = "ecal:hbhe:hf:pixel:sistrip:es:muon",
196  #DebugOn = True,
197  alwaysPass = False
198  )
199 )
200 
201 jetDQMAnalyzerAkPU3Calo = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
202  JetType = cms.untracked.string('calo'),
203  UEAlgo = cms.untracked.string('Pu'),
204  OutputFile = cms.untracked.string(''),
205  src = cms.InputTag("akPu3CaloJets"),
206  PFcands = cms.InputTag("particleFlowTmp"),
207  Background = cms.InputTag("voronoiBackgroundCalo"),
208  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
209  centralitycollection = cms.InputTag("hiCentrality"),
210  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
211  JetCorrections = cms.string(""),
212  recoJetPtThreshold = cms.double(10),
213  RThreshold = cms.double(0.3),
214  reverseEnergyFractionThreshold = cms.double(0.5)
215 )
216 jetDQMAnalyzerAkPU4Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu4CaloJets")
217 jetDQMAnalyzerAkPU5Calo = jetDQMAnalyzerAkPU3Calo.clone(src = "akPu5CaloJets")
218 
219 jetDQMAnalyzerAkPU3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
220  JetType = cms.untracked.string('pf'),
221  UEAlgo = cms.untracked.string('Pu'),
222  OutputFile = cms.untracked.string(''),
223  src = cms.InputTag("akPu3PFJets"),
224  PFcands = cms.InputTag("particleFlowTmp"),
225  Background = cms.InputTag("voronoiBackgroundPF"),
226  centralitycollection = cms.InputTag("hiCentrality"),
227  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
228  JetCorrections = cms.string(""),
229  recoJetPtThreshold = cms.double(10),
230  RThreshold = cms.double(0.3),
231  reverseEnergyFractionThreshold = cms.double(0.5)
232 )
233 jetDQMAnalyzerAkPU4PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu4PFJets")
234 jetDQMAnalyzerAkPU5PF = jetDQMAnalyzerAkPU3PF.clone(src = "akPu5PFJets")
235 
236 jetDQMAnalyzerAkCs3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons',
237  JetType = cms.untracked.string('pf'),
238  UEAlgo = cms.untracked.string('Cs'),
239  OutputFile = cms.untracked.string(''),
240  src = cms.InputTag("akCs3PFJets"),
241  CScands = cms.InputTag("akCs3PFJets","pfParticlesCs"),
242  PFcands = cms.InputTag("particleFlowTmp"),
243  centralitycollection = cms.InputTag("hiCentrality"),
244  #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
245  JetCorrections = cms.string(""),
246  recoJetPtThreshold = cms.double(10),
247  RThreshold = cms.double(0.3),
248  reverseEnergyFractionThreshold = cms.double(0.5),
249  etaMap = cms.InputTag('hiFJRhoProducer','mapEtaEdges'),
250  rho = cms.InputTag('hiFJRhoProducer','mapToRho'),
251  rhom = cms.InputTag('hiFJRhoProducer','mapToRhoM')
252 )
253 
254 jetDQMAnalyzerAkCs4PF=jetDQMAnalyzerAkCs3PF.clone(src = "akCs4PFJets",
255  CScands = "akCs4PFJets:pfParticlesCs"
256 )
257 
258 
259 jetDQMMatchAkPu3CaloAkPu3PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
260  src_Jet1 = cms.InputTag("akPu3CaloJets"),
261  src_Jet2 = cms.InputTag("akPu3PFJets"),
262  Jet1 = cms.untracked.string("PuCalo"),
263  Jet2 = cms.untracked.string("PuPF"),
264  recoJetPtThreshold = cms.double(20.),
265  recoDelRMatch = cms.double(0.2),
266  recoJetEtaCut = cms.double(2.0)
267 )
268 
269 jetDQMMatchAkPu4CaloAkPu4PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
270  src_Jet1 = cms.InputTag("akPu4CaloJets"),
271  src_Jet2 = cms.InputTag("akPu4PFJets"),
272  Jet1 = cms.untracked.string("PuCalo"),
273  Jet2 = cms.untracked.string("PuPF"),
274  recoJetPtThreshold = cms.double(20.),
275  recoDelRMatch = cms.double(0.2),
276  recoJetEtaCut = cms.double(2.0)
277 )
278 
279 jetDQMMatchAkPu5CaloAkPu5PF = DQMEDAnalyzer('JetAnalyzer_HeavyIons_matching',
280  src_Jet1 = cms.InputTag("akPu5CaloJets"),
281  src_Jet2 = cms.InputTag("akPu5PFJets"),
282  Jet1 = cms.untracked.string("PuCalo"),
283  Jet2 = cms.untracked.string("PuPF"),
284  recoJetPtThreshold = cms.double(20.),
285  recoDelRMatch = cms.double(0.2),
286  recoJetEtaCut = cms.double(2.0)
287 )
288