CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JetValidationHeavyIons_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 ### genjet cleaning for improved matching in HI environment
4 
6 
7 iterativeCone5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets'))
8 #iterativeCone7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone7HiGenJets'))
9 ak2HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak2HiGenJets'))
10 ak3HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets'))
11 ak4HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak4HiGenJets'))
12 ak5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets'))
13 
14 ak7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak7HiGenJets'))
15 
16 ### jet analyzer for several radii
17 ### iterative cone with PU, anti-kt with PU, anti-kt with fastjet PU, anti-kt with Vs
18 
19 
20 JetAnalyzerICPU5Calo = cms.EDAnalyzer("JetTester_HeavyIons",
21  JetType = cms.untracked.string('calo'),
22  UEAlgo = cms.untracked.string('Pu'),
23  OutputFile = cms.untracked.string(''),
24  src = cms.InputTag("iterativeConePu5CaloJets"),
25  srcGen = cms.InputTag("iterativeCone5HiCleanedGenJets"),
26  PFcands = cms.InputTag("particleFlowTmp"),
27  Background = cms.InputTag("voronoiBackgroundCalo"),
28  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
29  centralitycollection = cms.InputTag("hiCentrality"),
30  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
31  JetCorrections = cms.string(""),
32  recoJetPtThreshold = cms.double(10),
33  genEnergyFractionThreshold = cms.double(0.05),
34  genPtThreshold = cms.double(1.0),
35  matchGenPtThreshold = cms.double(20.0),
36  RThreshold = cms.double(0.3),
37  reverseEnergyFractionThreshold = cms.double(0.5)
38 )
39 '''
40 JetAnalyzerICPU7Calo = cms.EDAnalyzer("JetTester_HeavyIons",
41  JetType = cms.untracked.string('calo'),
42  UEAlgo = cms.untracked.string('Pu'),
43  OutputFile = cms.untracked.string(''),
44  src = cms.InputTag("iterativeConePu7CaloJets"),
45  srcGen = cms.InputTag("iterativeCone7HiCleanedGenJets"),
46  #srcRho = cms.InputTag("iterativeConePu7CaloJets","rho"),
47  JetCorrections = cms.string(""),
48  recoJetPtThreshold = cms.double(10),
49  genEnergyFractionThreshold = cms.double(0.05),
50  genPtThreshold = cms.double(1.0),
51  matchGenPtThreshold = cms.double(20.0),
52  RThreshold = cms.double(0.3),
53  reverseEnergyFractionThreshold = cms.double(0.5)
54 )
55 '''
56 JetAnalyzerAkPU3Calo = cms.EDAnalyzer("JetTester_HeavyIons",
57  JetType = cms.untracked.string('calo'),
58  UEAlgo = cms.untracked.string('Pu'),
59  OutputFile = cms.untracked.string(''),
60  src = cms.InputTag("akPu3CaloJets"),
61  srcGen = cms.InputTag("ak3HiCleanedGenJets"),
62  PFcands = cms.InputTag("particleFlowTmp"),
63  Background = cms.InputTag("voronoiBackgroundCalo"),
64  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
65  centralitycollection = cms.InputTag("hiCentrality"),
66  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
67  JetCorrections = cms.string(""),
68  recoJetPtThreshold = cms.double(10),
69  genEnergyFractionThreshold = cms.double(0.05),
70  genPtThreshold = cms.double(1.0),
71  matchGenPtThreshold = cms.double(20.0),
72  RThreshold = cms.double(0.3),
73  reverseEnergyFractionThreshold = cms.double(0.5)
74 )
75 
76 JetAnalyzerAkPU4Calo = cms.EDAnalyzer("JetTester_HeavyIons",
77  JetType = cms.untracked.string('calo'),
78  UEAlgo = cms.untracked.string('Pu'),
79  OutputFile = cms.untracked.string(''),
80  src = cms.InputTag("akPu4CaloJets"),
81  srcGen = cms.InputTag("ak4HiCleanedGenJets"),
82  PFcands = cms.InputTag("particleFlowTmp"),
83  Background = cms.InputTag("voronoiBackgroundCalo"),
84  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
85  centralitycollection = cms.InputTag("hiCentrality"),
86  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
87  JetCorrections = cms.string(""),
88  recoJetPtThreshold = cms.double(10),
89  genEnergyFractionThreshold = cms.double(0.05),
90  genPtThreshold = cms.double(1.0),
91  matchGenPtThreshold = cms.double(20.0),
92  RThreshold = cms.double(0.3),
93  reverseEnergyFractionThreshold = cms.double(0.5)
94 )
95 
96 JetAnalyzerAkPU5Calo = cms.EDAnalyzer("JetTester_HeavyIons",
97  JetType = cms.untracked.string('calo'),
98  UEAlgo = cms.untracked.string('Pu'),
99  OutputFile = cms.untracked.string(''),
100  src = cms.InputTag("akPu5CaloJets"),
101  srcGen = cms.InputTag("ak5HiCleanedGenJets"),
102  PFcands = cms.InputTag("particleFlowTmp"),
103  Background = cms.InputTag("voronoiBackgroundCalo"),
104  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
105  centralitycollection = cms.InputTag("hiCentrality"),
106  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
107  JetCorrections = cms.string(""),
108  recoJetPtThreshold = cms.double(10),
109  genEnergyFractionThreshold = cms.double(0.05),
110  genPtThreshold = cms.double(1.0),
111  matchGenPtThreshold = cms.double(20.0),
112  RThreshold = cms.double(0.3),
113  reverseEnergyFractionThreshold = cms.double(0.5)
114 )
115 
116 JetAnalyzerAkPU3PF = cms.EDAnalyzer("JetTester_HeavyIons",
117  JetType = cms.untracked.string('pf'),
118  UEAlgo = cms.untracked.string('Pu'),
119  OutputFile = cms.untracked.string(''),
120  src = cms.InputTag("akPu3PFJets"),
121  srcGen = cms.InputTag("ak3HiCleanedGenJets"),
122  PFcands = cms.InputTag("particleFlowTmp"),
123  Background = cms.InputTag("voronoiBackgroundPF"),
124  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
125  centralitycollection = cms.InputTag("hiCentrality"),
126  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
127  JetCorrections = cms.string(""),
128  recoJetPtThreshold = cms.double(10),
129  genEnergyFractionThreshold = cms.double(0.05),
130  genPtThreshold = cms.double(1.0),
131  matchGenPtThreshold = cms.double(20.0),
132  RThreshold = cms.double(0.3),
133  reverseEnergyFractionThreshold = cms.double(0.5)
134 )
135 
136 JetAnalyzerAkPU4PF = cms.EDAnalyzer("JetTester_HeavyIons",
137  JetType = cms.untracked.string('pf'),
138  UEAlgo = cms.untracked.string('Pu'),
139  OutputFile = cms.untracked.string(''),
140  src = cms.InputTag("akPu4PFJets"),
141  srcGen = cms.InputTag("ak4HiCleanedGenJets"),
142  PFcands = cms.InputTag("particleFlowTmp"),
143  Background = cms.InputTag("voronoiBackgroundPF"),
144  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
145  centralitycollection = cms.InputTag("hiCentrality"),
146  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
147  JetCorrections = cms.string(""),
148  recoJetPtThreshold = cms.double(10),
149  genEnergyFractionThreshold = cms.double(0.05),
150  genPtThreshold = cms.double(1.0),
151  matchGenPtThreshold = cms.double(20.0),
152  RThreshold = cms.double(0.3),
153  reverseEnergyFractionThreshold = cms.double(0.5)
154 )
155 
156 JetAnalyzerAkPU5PF = cms.EDAnalyzer("JetTester_HeavyIons",
157  JetType = cms.untracked.string('pf'),
158  UEAlgo = cms.untracked.string('Pu'),
159  OutputFile = cms.untracked.string(''),
160  src = cms.InputTag("akPu5PFJets"),
161  srcGen = cms.InputTag("ak5HiCleanedGenJets"),
162  PFcands = cms.InputTag("particleFlowTmp"),
163  Background = cms.InputTag("voronoiBackgroundPF"),
164  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
165  centralitycollection = cms.InputTag("hiCentrality"),
166  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
167  JetCorrections = cms.string(""),
168  recoJetPtThreshold = cms.double(10),
169  genEnergyFractionThreshold = cms.double(0.05),
170  genPtThreshold = cms.double(1.0),
171  matchGenPtThreshold = cms.double(20.0),
172  RThreshold = cms.double(0.3),
173  reverseEnergyFractionThreshold = cms.double(0.5)
174 )
175 
176 JetAnalyzerAkVs2Calo = cms.EDAnalyzer("JetTester_HeavyIons",
177  JetType = cms.untracked.string('calo'),
178  UEAlgo = cms.untracked.string('Vs'),
179  OutputFile = cms.untracked.string(''),
180  src = cms.InputTag("akVs2CaloJets"),
181  srcGen = cms.InputTag("ak2HiCleanedGenJets"),
182  PFcands = cms.InputTag("particleFlowTmp"),
183  #Cands = cms.InputTag("caloTowers"),
184  Background = cms.InputTag("voronoiBackgroundCalo"),
185  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
186  centralitycollection = cms.InputTag("hiCentrality"),
187  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
188  JetCorrections = cms.string(""),
189  recoJetPtThreshold = cms.double(10),
190  genEnergyFractionThreshold = cms.double(0.05),
191  genPtThreshold = cms.double(1.0),
192  matchGenPtThreshold = cms.double(20.0),
193  RThreshold = cms.double(0.3),
194  reverseEnergyFractionThreshold = cms.double(0.5)
195 )
196 
197 JetAnalyzerAkVs3Calo = cms.EDAnalyzer("JetTester_HeavyIons",
198  JetType = cms.untracked.string('calo'),
199  UEAlgo = cms.untracked.string('Vs'),
200  OutputFile = cms.untracked.string(''),
201  src = cms.InputTag("akVs3CaloJets"),
202  srcGen = cms.InputTag("ak3HiCleanedGenJets"),
203  PFcands = cms.InputTag("particleFlowTmp"),
204  Background = cms.InputTag("voronoiBackgroundCalo"),
205  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
206  centralitycollection = cms.InputTag("hiCentrality"),
207  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
208  JetCorrections = cms.string(""),
209  recoJetPtThreshold = cms.double(10),
210  genEnergyFractionThreshold = cms.double(0.05),
211  genPtThreshold = cms.double(1.0),
212  matchGenPtThreshold = cms.double(20.0),
213  RThreshold = cms.double(0.3),
214  reverseEnergyFractionThreshold = cms.double(0.5)
215 )
216 
217 JetAnalyzerAkVs4Calo = cms.EDAnalyzer("JetTester_HeavyIons",
218  JetType = cms.untracked.string('calo'),
219  UEAlgo = cms.untracked.string('Vs'),
220  OutputFile = cms.untracked.string(''),
221  src = cms.InputTag("akVs4CaloJets"),
222  srcGen = cms.InputTag("ak4HiCleanedGenJets"),
223  PFcands = cms.InputTag("particleFlowTmp"),
224  Background = cms.InputTag("voronoiBackgroundCalo"),
225  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
226  centralitycollection = cms.InputTag("hiCentrality"),
227  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
228  JetCorrections = cms.string(""),
229  recoJetPtThreshold = cms.double(10),
230  genEnergyFractionThreshold = cms.double(0.05),
231  genPtThreshold = cms.double(1.0),
232  matchGenPtThreshold = cms.double(20.0),
233  RThreshold = cms.double(0.3),
234  reverseEnergyFractionThreshold = cms.double(0.5)
235 )
236 
237 JetAnalyzerAkVs5Calo = cms.EDAnalyzer("JetTester_HeavyIons",
238  JetType = cms.untracked.string('calo'),
239  UEAlgo = cms.untracked.string('Vs'),
240  OutputFile = cms.untracked.string(''),
241  src = cms.InputTag("akVs5CaloJets"),
242  srcGen = cms.InputTag("ak5HiCleanedGenJets"),
243  PFcands = cms.InputTag("particleFlowTmp"),
244  Background = cms.InputTag("voronoiBackgroundCalo"),
245  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
246  centralitycollection = cms.InputTag("hiCentrality"),
247  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
248  JetCorrections = cms.string(""),
249  recoJetPtThreshold = cms.double(10),
250  genEnergyFractionThreshold = cms.double(0.05),
251  genPtThreshold = cms.double(1.0),
252  matchGenPtThreshold = cms.double(20.0),
253  RThreshold = cms.double(0.3),
254  reverseEnergyFractionThreshold = cms.double(0.5)
255 )
256 
257 JetAnalyzerAkVs3PF = cms.EDAnalyzer("JetTester_HeavyIons",
258  JetType = cms.untracked.string('pf'),
259  UEAlgo = cms.untracked.string('Vs'),
260  OutputFile = cms.untracked.string(''),
261  src = cms.InputTag("akVs3PFJets"),
262  srcGen = cms.InputTag("ak3HiCleanedGenJets"),
263  PFcands = cms.InputTag("particleFlowTmp"),
264  Background = cms.InputTag("voronoiBackgroundPF"),
265  #srcRho = cms.InputTag("akVs3PFJets","rho"),
266  centralitycollection = cms.InputTag("hiCentrality"),
267  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
268  JetCorrections = cms.string(""),
269  recoJetPtThreshold = cms.double(10),
270  genEnergyFractionThreshold = cms.double(0.05),
271  genPtThreshold = cms.double(1.0),
272  matchGenPtThreshold = cms.double(20.0),
273  RThreshold = cms.double(0.3),
274  reverseEnergyFractionThreshold = cms.double(0.5)
275 )
276 
277 
278 JetAnalyzerAkVs4PF = cms.EDAnalyzer("JetTester_HeavyIons",
279  JetType = cms.untracked.string('pf'),
280  UEAlgo = cms.untracked.string('Vs'),
281  OutputFile = cms.untracked.string(''),
282  src = cms.InputTag("akVs4PFJets"),
283  srcGen = cms.InputTag("ak4HiCleanedGenJets"),
284  PFcands = cms.InputTag("particleFlowTmp"),
285  Background = cms.InputTag("voronoiBackgroundPF"),
286  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
287  centralitycollection = cms.InputTag("hiCentrality"),
288  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
289  JetCorrections = cms.string(""),
290  recoJetPtThreshold = cms.double(10),
291  genEnergyFractionThreshold = cms.double(0.05),
292  genPtThreshold = cms.double(1.0),
293  matchGenPtThreshold = cms.double(20.0),
294  RThreshold = cms.double(0.3),
295  reverseEnergyFractionThreshold = cms.double(0.5)
296 )
297 
298 JetAnalyzerAkVs5PF = cms.EDAnalyzer("JetTester_HeavyIons",
299  JetType = cms.untracked.string('pf'),
300  UEAlgo = cms.untracked.string('Vs'),
301  OutputFile = cms.untracked.string(''),
302  src = cms.InputTag("akVs5PFJets"),
303  srcGen = cms.InputTag("ak5HiCleanedGenJets"),
304  PFcands = cms.InputTag("particleFlowTmp"),
305  Background = cms.InputTag("voronoiBackgroundPF"),
306  #srcRho = cms.InputTag("iterativeConePu5CaloJets","rho"),
307  centralitycollection = cms.InputTag("hiCentrality"),
308  centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
309  JetCorrections = cms.string(""),
310  recoJetPtThreshold = cms.double(10),
311  genEnergyFractionThreshold = cms.double(0.05),
312  genPtThreshold = cms.double(1.0),
313  matchGenPtThreshold = cms.double(20.0),
314  RThreshold = cms.double(0.3),
315  reverseEnergyFractionThreshold = cms.double(0.5)
316 )
317 
318 hiJetValidation = cms.Sequence(
319  ak3HiCleanedGenJets
320  * ak4HiCleanedGenJets
321  * ak5HiCleanedGenJets
322  * JetAnalyzerAkPU3Calo
323  * JetAnalyzerAkPU4Calo
324  * JetAnalyzerAkPU5Calo
325 
326  * JetAnalyzerAkPU3PF
327  * JetAnalyzerAkPU4PF
328  * JetAnalyzerAkPU5PF
329 
330  * JetAnalyzerAkVs2Calo
331  * JetAnalyzerAkVs3Calo
332  * JetAnalyzerAkVs4Calo
333 
334  * JetAnalyzerAkVs3PF
335  * JetAnalyzerAkVs4PF
336  * JetAnalyzerAkVs5PF
337 
338 )