ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitProd/TreeFiller/python/JetsMCFlavourMatching_cfi.py
Revision: 1.8
Committed: Wed Jul 22 04:41:00 2009 UTC (15 years, 9 months ago) by bendavid
Content type: text/x-python
Branch: MAIN
CVS Tags: Mit_011a, Mit_011, Mit_010a, Mit_010
Changes since 1.7: +101 -3 lines
Log Message:
Add antikt jets

File Contents

# User Rev Content
1 bendavid 1.8 # $Id: JetsMCFlavourMatching_cfi.py,v 1.7 2009/03/22 08:48:59 loizides Exp $
2 loizides 1.3
3 sixie 1.1 import FWCore.ParameterSet.Config as cms
4    
5     # select the partons for Jet MC Flavour
6     myPartons = cms.EDFilter("PartonSelector",
7     withLeptons = cms.bool(False)
8     )
9    
10     #Iterative Cone
11     # Flavour byReference
12     IC5byRef = cms.EDFilter("JetPartonMatcher",
13 loizides 1.7 jets = cms.InputTag("iterativeCone5CaloJets"),
14 sixie 1.1 coneSizeToAssociate = cms.double(0.3),
15 loizides 1.7 partons = cms.InputTag("myPartons")
16 sixie 1.1 )
17    
18     # Flavour byValue PhysDef
19     IC5byValPhys = cms.EDFilter("JetFlavourIdentifier",
20 loizides 1.7 srcByReference = cms.InputTag("IC5byRef"),
21 sixie 1.1 physicsDefinition = cms.bool(True)
22     )
23    
24     # Flavour byValue AlgoDef
25     IC5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
26 loizides 1.7 srcByReference = cms.InputTag("IC5byRef"),
27 sixie 1.1 physicsDefinition = cms.bool(False)
28     )
29    
30     #SisCone Jets
31     # Flavour byReference
32     SC5byRef = cms.EDFilter("JetPartonMatcher",
33 loizides 1.7 jets = cms.InputTag("sisCone5CaloJets"),
34 sixie 1.1 coneSizeToAssociate = cms.double(0.3),
35 loizides 1.7 partons = cms.InputTag("myPartons")
36 sixie 1.1 )
37    
38     # Flavour byValue PhysDef
39     SC5byValPhys = cms.EDFilter("JetFlavourIdentifier",
40 loizides 1.7 srcByReference = cms.InputTag("SC5byRef"),
41 sixie 1.1 physicsDefinition = cms.bool(True)
42     )
43    
44     # Flavour byValue AlgoDef
45     SC5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
46 loizides 1.7 srcByReference = cms.InputTag("SC5byRef"),
47 sixie 1.1 physicsDefinition = cms.bool(False)
48     )
49    
50     # Flavour byReference
51     SC7byRef = cms.EDFilter("JetPartonMatcher",
52 loizides 1.7 jets = cms.InputTag("sisCone7CaloJets"),
53 sixie 1.1 coneSizeToAssociate = cms.double(0.3),
54 loizides 1.7 partons = cms.InputTag("myPartons")
55 sixie 1.1 )
56    
57     # Flavour byValue PhysDef
58     SC7byValPhys = cms.EDFilter("JetFlavourIdentifier",
59 loizides 1.7 srcByReference = cms.InputTag("SC7byRef"),
60 sixie 1.1 physicsDefinition = cms.bool(True)
61     )
62    
63     # Flavour byValue AlgoDef
64     SC7byValAlgo = cms.EDFilter("JetFlavourIdentifier",
65 loizides 1.7 srcByReference = cms.InputTag("SC7byRef"),
66 sixie 1.1 physicsDefinition = cms.bool(False)
67     )
68    
69     #Kt Jets
70     # Flavour byReference
71     KT4byRef = cms.EDFilter("JetPartonMatcher",
72 loizides 1.7 jets = cms.InputTag("kt4CaloJets"),
73 sixie 1.1 coneSizeToAssociate = cms.double(0.3),
74 loizides 1.7 partons = cms.InputTag("myPartons")
75 sixie 1.1 )
76    
77     # Flavour byValue PhysDef
78     KT4byValPhys = cms.EDFilter("JetFlavourIdentifier",
79 loizides 1.7 srcByReference = cms.InputTag("KT4byRef"),
80 sixie 1.1 physicsDefinition = cms.bool(True)
81     )
82    
83     # Flavour byValue AlgoDef
84     KT4byValAlgo = cms.EDFilter("JetFlavourIdentifier",
85 loizides 1.7 srcByReference = cms.InputTag("KT4byRef"),
86 sixie 1.1 physicsDefinition = cms.bool(False)
87     )
88    
89    
90     # Flavour byReference
91     KT6byRef = cms.EDFilter("JetPartonMatcher",
92 loizides 1.7 jets = cms.InputTag("kt6CaloJets"),
93 sixie 1.1 coneSizeToAssociate = cms.double(0.3),
94 loizides 1.7 partons = cms.InputTag("myPartons")
95 sixie 1.1 )
96    
97     # Flavour byValue PhysDef
98     KT6byValPhys = cms.EDFilter("JetFlavourIdentifier",
99 loizides 1.7 srcByReference = cms.InputTag("KT6byRef"),
100 sixie 1.1 physicsDefinition = cms.bool(True)
101     )
102    
103     # Flavour byValue AlgoDef
104     KT6byValAlgo = cms.EDFilter("JetFlavourIdentifier",
105 loizides 1.7 srcByReference = cms.InputTag("KT6byRef"),
106 sixie 1.1 physicsDefinition = cms.bool(False)
107     )
108    
109 bendavid 1.8 #AntiKt Jets
110     # Flavour byReference
111     AKT5byRef = cms.EDFilter("JetPartonMatcher",
112     jets = cms.InputTag("antikt5CaloJets"),
113     coneSizeToAssociate = cms.double(0.3),
114     partons = cms.InputTag("myPartons")
115     )
116    
117     # Flavour byValue PhysDef
118     AKT5byValPhys = cms.EDFilter("JetFlavourIdentifier",
119     srcByReference = cms.InputTag("AKT5byRef"),
120     physicsDefinition = cms.bool(True)
121     )
122    
123     # Flavour byValue AlgoDef
124     AKT5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
125     srcByReference = cms.InputTag("AKT5byRef"),
126     physicsDefinition = cms.bool(False)
127     )
128    
129     # Flavour byReference
130     AKT7byRef = cms.EDFilter("JetPartonMatcher",
131     jets = cms.InputTag("antikt7CaloJets"),
132     coneSizeToAssociate = cms.double(0.3),
133     partons = cms.InputTag("myPartons")
134     )
135    
136     # Flavour byValue PhysDef
137     AKT7byValPhys = cms.EDFilter("JetFlavourIdentifier",
138     srcByReference = cms.InputTag("AKT7byRef"),
139     physicsDefinition = cms.bool(True)
140     )
141    
142     # Flavour byValue AlgoDef
143     AKT7byValAlgo = cms.EDFilter("JetFlavourIdentifier",
144     srcByReference = cms.InputTag("AKT7byRef"),
145     physicsDefinition = cms.bool(False)
146     )
147    
148 bendavid 1.5 #IC5PF Jets
149     # Flavour byReference
150     IC5PFbyRef = cms.EDFilter("JetPartonMatcher",
151 loizides 1.7 jets = cms.InputTag("iterativeCone5PFJets"),
152 bendavid 1.5 coneSizeToAssociate = cms.double(0.3),
153 loizides 1.7 partons = cms.InputTag("myPartons")
154 bendavid 1.5 )
155    
156     # Flavour byValue PhysDef
157     IC5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
158 loizides 1.7 srcByReference = cms.InputTag("IC5PFbyRef"),
159 bendavid 1.5 physicsDefinition = cms.bool(True)
160     )
161    
162     # Flavour byValue AlgoDef
163     IC5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
164 loizides 1.7 srcByReference = cms.InputTag("IC5PFbyRef"),
165 bendavid 1.5 physicsDefinition = cms.bool(False)
166     )
167    
168     #SisCone Jets
169     # Flavour byReference
170     SC5PFbyRef = cms.EDFilter("JetPartonMatcher",
171 loizides 1.7 jets = cms.InputTag("sisCone5PFJets"),
172 bendavid 1.5 coneSizeToAssociate = cms.double(0.3),
173 loizides 1.7 partons = cms.InputTag("myPartons")
174 bendavid 1.5 )
175    
176     # Flavour byValue PhysDef
177     SC5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
178 loizides 1.7 srcByReference = cms.InputTag("SC5PFbyRef"),
179 bendavid 1.5 physicsDefinition = cms.bool(True)
180     )
181    
182     # Flavour byValue AlgoDef
183     SC5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
184 loizides 1.7 srcByReference = cms.InputTag("SC5PFbyRef"),
185 bendavid 1.5 physicsDefinition = cms.bool(False)
186     )
187    
188     # Flavour byReference
189     SC7PFbyRef = cms.EDFilter("JetPartonMatcher",
190 loizides 1.7 jets = cms.InputTag("sisCone7PFJets"),
191 bendavid 1.5 coneSizeToAssociate = cms.double(0.3),
192 loizides 1.7 partons = cms.InputTag("myPartons")
193 bendavid 1.5 )
194    
195     # Flavour byValue PhysDef
196     SC7PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
197 loizides 1.7 srcByReference = cms.InputTag("SC7PFbyRef"),
198 bendavid 1.5 physicsDefinition = cms.bool(True)
199     )
200    
201     # Flavour byValue AlgoDef
202     SC7PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
203 loizides 1.7 srcByReference = cms.InputTag("SC7PFbyRef"),
204 bendavid 1.5 physicsDefinition = cms.bool(False)
205     )
206    
207     #Kt Jets
208     # Flavour byReference
209     KT4PFbyRef = cms.EDFilter("JetPartonMatcher",
210 loizides 1.7 jets = cms.InputTag("kt4PFJets"),
211 bendavid 1.5 coneSizeToAssociate = cms.double(0.3),
212 loizides 1.7 partons = cms.InputTag("myPartons")
213 bendavid 1.5 )
214    
215     # Flavour byValue PhysDef
216     KT4PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
217 loizides 1.7 srcByReference = cms.InputTag("KT4PFbyRef"),
218 bendavid 1.5 physicsDefinition = cms.bool(True)
219     )
220    
221     # Flavour byValue AlgoDef
222     KT4PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
223 loizides 1.7 srcByReference = cms.InputTag("KT4PFbyRef"),
224 bendavid 1.5 physicsDefinition = cms.bool(False)
225     )
226    
227    
228     # Flavour byReference
229     KT6PFbyRef = cms.EDFilter("JetPartonMatcher",
230 loizides 1.7 jets = cms.InputTag("kt6PFJets"),
231 bendavid 1.5 coneSizeToAssociate = cms.double(0.3),
232 loizides 1.7 partons = cms.InputTag("myPartons")
233 bendavid 1.5 )
234    
235     # Flavour byValue PhysDef
236     KT6PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
237 loizides 1.7 srcByReference = cms.InputTag("KT6PFbyRef"),
238 bendavid 1.5 physicsDefinition = cms.bool(True)
239     )
240    
241     # Flavour byValue AlgoDef
242     KT6PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
243 loizides 1.7 srcByReference = cms.InputTag("KT6PFbyRef"),
244 bendavid 1.5 physicsDefinition = cms.bool(False)
245     )
246    
247 bendavid 1.8 #AntiKt Jets
248     # Flavour byReference
249     AKT5PFbyRef = cms.EDFilter("JetPartonMatcher",
250     jets = cms.InputTag("antikt5PFJets"),
251     coneSizeToAssociate = cms.double(0.3),
252     partons = cms.InputTag("myPartons")
253     )
254    
255     # Flavour byValue PhysDef
256     AKT5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
257     srcByReference = cms.InputTag("AKT5PFbyRef"),
258     physicsDefinition = cms.bool(True)
259     )
260    
261     # Flavour byValue AlgoDef
262     AKT5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
263     srcByReference = cms.InputTag("AKT5PFbyRef"),
264     physicsDefinition = cms.bool(False)
265     )
266    
267     # Flavour byReference
268     AKT7PFbyRef = cms.EDFilter("JetPartonMatcher",
269     jets = cms.InputTag("antikt7PFJets"),
270     coneSizeToAssociate = cms.double(0.3),
271     partons = cms.InputTag("myPartons")
272     )
273    
274     # Flavour byValue PhysDef
275     AKT7PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
276     srcByReference = cms.InputTag("AKT7PFbyRef"),
277     physicsDefinition = cms.bool(True)
278     )
279    
280     # Flavour byValue AlgoDef
281     AKT7PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
282     srcByReference = cms.InputTag("AKT7PFbyRef"),
283     physicsDefinition = cms.bool(False)
284     )
285    
286 sixie 1.2 #Flavour by Reference
287     GenIC5byRef = cms.EDFilter("JetPartonMatcher",
288 loizides 1.7 jets = cms.InputTag("iterativeCone5GenJets"),
289 sixie 1.2 coneSizeToAssociate = cms.double(0.3),
290 loizides 1.7 partons = cms.InputTag("myPartons")
291 sixie 1.2 )
292    
293     #Flavour by Reference
294     GenSC5byRef = cms.EDFilter("JetPartonMatcher",
295 loizides 1.7 jets = cms.InputTag("sisCone5GenJets"),
296 sixie 1.2 coneSizeToAssociate = cms.double(0.3),
297 loizides 1.7 partons = cms.InputTag("myPartons")
298 sixie 1.2 )
299 sixie 1.1
300 sixie 1.4 #Flavour by Reference
301     GenSC7byRef = cms.EDFilter("JetPartonMatcher",
302 loizides 1.7 jets = cms.InputTag("sisCone7GenJets"),
303 sixie 1.4 coneSizeToAssociate = cms.double(0.3),
304 loizides 1.7 partons = cms.InputTag("myPartons")
305 sixie 1.4 )
306    
307     #Flavour by Reference
308     GenKT4byRef = cms.EDFilter("JetPartonMatcher",
309 loizides 1.7 jets = cms.InputTag("kt4GenJets"),
310 sixie 1.4 coneSizeToAssociate = cms.double(0.3),
311 loizides 1.7 partons = cms.InputTag("myPartons")
312 sixie 1.4 )
313    
314     #Flavour by Reference
315     GenKT6byRef = cms.EDFilter("JetPartonMatcher",
316 loizides 1.7 jets = cms.InputTag("kt6GenJets"),
317 sixie 1.4 coneSizeToAssociate = cms.double(0.3),
318 loizides 1.7 partons = cms.InputTag("myPartons")
319 sixie 1.4 )
320    
321 bendavid 1.8 #Flavour by Reference
322     GenAKT5byRef = cms.EDFilter("JetPartonMatcher",
323     jets = cms.InputTag("antikt5GenJets"),
324     coneSizeToAssociate = cms.double(0.3),
325     partons = cms.InputTag("myPartons")
326     )
327    
328     #Flavour by Reference
329     GenAKT7byRef = cms.EDFilter("JetPartonMatcher",
330     jets = cms.InputTag("antikt7GenJets"),
331     coneSizeToAssociate = cms.double(0.3),
332     partons = cms.InputTag("myPartons")
333     )
334    
335 sixie 1.1 iterativeCone5Flavour = cms.Sequence(IC5byRef*IC5byValPhys*IC5byValAlgo)
336     sisCone5Flavour = cms.Sequence(SC5byRef*SC5byValPhys*SC5byValAlgo)
337     sisCone7Flavour = cms.Sequence(SC7byRef*SC7byValPhys*SC7byValAlgo)
338     KT4Flavour = cms.Sequence(KT4byRef*KT4byValPhys*KT4byValAlgo)
339     KT6Flavour = cms.Sequence(KT6byRef*KT6byValPhys*KT6byValAlgo)
340 bendavid 1.8 AKT5Flavour = cms.Sequence(AKT5byRef*AKT5byValPhys*AKT5byValAlgo)
341     AKT7Flavour = cms.Sequence(AKT7byRef*AKT7byValPhys*AKT7byValAlgo)
342 bendavid 1.5 iterativeCone5PFFlavour = cms.Sequence(IC5PFbyRef*IC5PFbyValPhys*IC5PFbyValAlgo)
343     sisCone5PFFlavour = cms.Sequence(SC5PFbyRef*SC5PFbyValPhys*SC5PFbyValAlgo)
344     sisCone7PFFlavour = cms.Sequence(SC7PFbyRef*SC7PFbyValPhys*SC7PFbyValAlgo)
345     KT4PFFlavour = cms.Sequence(KT4PFbyRef*KT4PFbyValPhys*KT4PFbyValAlgo)
346     KT6PFFlavour = cms.Sequence(KT6PFbyRef*KT6PFbyValPhys*KT6PFbyValAlgo)
347 bendavid 1.8 AKT5PFFlavour = cms.Sequence(AKT5PFbyRef*AKT5PFbyValPhys*AKT5PFbyValAlgo)
348     AKT7PFFlavour = cms.Sequence(AKT7PFbyRef*AKT7PFbyValPhys*AKT7PFbyValAlgo)
349 sixie 1.4 GenJetFlavour = cms.Sequence(GenIC5byRef*GenSC5byRef*GenSC7byRef*GenKT4byRef*GenKT6byRef)
350 loizides 1.7 caloJetMCFlavour = cms.Sequence(myPartons*iterativeCone5Flavour
351     +sisCone5Flavour
352     +sisCone7Flavour
353     +KT4Flavour
354     +KT6Flavour
355 bendavid 1.6 +GenJetFlavour)
356 bendavid 1.8
357     pfJetMCFlavour = cms.Sequence(iterativeCone5PFFlavour+sisCone5PFFlavour+sisCone7PFFlavour+KT4PFFlavour+KT6PFFlavour)
358    
359     antiKtJetMCFlavour = cms.Sequence(GenAKT5byRef*GenAKT7byRef*AKT5Flavour*AKT7Flavour*AKT5PFFlavour*AKT7PFFlavour)