TreeFiller/python/JetsMCFlavourMatching_cfi.py

# $Id: JetsMCFlavourMatching_cfi.py,v 1.4 2008/09/29 17:06:29 sixie Exp $

import FWCore.ParameterSet.Config as cms

# select the partons for Jet MC Flavour
myPartons = cms.EDFilter("PartonSelector",
    withLeptons = cms.bool(False)
)

#Iterative Cone
# Flavour byReference
IC5byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("iterativeCone5CaloJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)
 
# Flavour byValue PhysDef
IC5byValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("IC5byRef"),
    physicsDefinition = cms.bool(True)
)
 
# Flavour byValue AlgoDef
IC5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("IC5byRef"),
    physicsDefinition = cms.bool(False)
)

#SisCone Jets
# Flavour byReference
SC5byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("sisCone5CaloJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
SC5byValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC5byRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
SC5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC5byRef"),
    physicsDefinition = cms.bool(False)
)

# Flavour byReference
SC7byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("sisCone7CaloJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
SC7byValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC7byRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
SC7byValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC7byRef"),
    physicsDefinition = cms.bool(False)
)

#Kt Jets
# Flavour byReference
KT4byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("kt4CaloJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
KT4byValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT4byRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
KT4byValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT4byRef"),
    physicsDefinition = cms.bool(False)
)


# Flavour byReference
KT6byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("kt6CaloJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
KT6byValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT6byRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
KT6byValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT6byRef"),
    physicsDefinition = cms.bool(False)
)

#IC5PF Jets
# Flavour byReference
IC5PFbyRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("iterativeCone5PFJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)
 
# Flavour byValue PhysDef
IC5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("IC5PFbyRef"),
    physicsDefinition = cms.bool(True)
)
 
# Flavour byValue AlgoDef
IC5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("IC5PFbyRef"),
    physicsDefinition = cms.bool(False)
)

#SisCone Jets
# Flavour byReference
SC5PFbyRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("sisCone5PFJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
SC5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC5PFbyRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
SC5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC5PFbyRef"),
    physicsDefinition = cms.bool(False)
)

# Flavour byReference
SC7PFbyRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("sisCone7PFJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
SC7PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC7PFbyRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
SC7PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("SC7PFbyRef"),
    physicsDefinition = cms.bool(False)
)

#Kt Jets
# Flavour byReference
KT4PFbyRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("kt4PFJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
KT4PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT4PFbyRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
KT4PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT4PFbyRef"),
    physicsDefinition = cms.bool(False)
)


# Flavour byReference
KT6PFbyRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("kt6PFJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

# Flavour byValue PhysDef
KT6PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT6PFbyRef"),
    physicsDefinition = cms.bool(True)
)

# Flavour byValue AlgoDef
KT6PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
    srcByReference = cms.InputTag("KT6PFbyRef"),
    physicsDefinition = cms.bool(False)
)


#Flavour by Reference
GenIC5byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("iterativeCone5GenJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

#Flavour by Reference
GenSC5byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("sisCone5GenJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

#Flavour by Reference
GenSC7byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("sisCone7GenJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

#Flavour by Reference
GenKT4byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("kt4GenJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

#Flavour by Reference
GenKT6byRef = cms.EDFilter("JetPartonMatcher",
    jets = cms.InputTag("kt6GenJets"),
    coneSizeToAssociate = cms.double(0.3),
    partons = cms.InputTag("myPartons")
)

iterativeCone5Flavour = cms.Sequence(IC5byRef*IC5byValPhys*IC5byValAlgo)
sisCone5Flavour = cms.Sequence(SC5byRef*SC5byValPhys*SC5byValAlgo)
sisCone7Flavour = cms.Sequence(SC7byRef*SC7byValPhys*SC7byValAlgo)
KT4Flavour = cms.Sequence(KT4byRef*KT4byValPhys*KT4byValAlgo)
KT6Flavour = cms.Sequence(KT6byRef*KT6byValPhys*KT6byValAlgo)
iterativeCone5PFFlavour = cms.Sequence(IC5PFbyRef*IC5PFbyValPhys*IC5PFbyValAlgo)
sisCone5PFFlavour = cms.Sequence(SC5PFbyRef*SC5PFbyValPhys*SC5PFbyValAlgo)
sisCone7PFFlavour = cms.Sequence(SC7PFbyRef*SC7PFbyValPhys*SC7PFbyValAlgo)
KT4PFFlavour = cms.Sequence(KT4PFbyRef*KT4PFbyValPhys*KT4PFbyValAlgo)
KT6PFFlavour = cms.Sequence(KT6PFbyRef*KT6PFbyValPhys*KT6PFbyValAlgo)
GenJetFlavour = cms.Sequence(GenIC5byRef*GenSC5byRef*GenSC7byRef*GenKT4byRef*GenKT6byRef)
caloJetMCFlavour = cms.Sequence(myPartons*iterativeCone5Flavour+sisCone5Flavour+sisCone7Flavour+KT4Flavour+KT6Flavour
                                +iterativeCone5PFFlavour+sisCone5PFFlavour+sisCone7PFFlavour+KT4PFFlavour+KT6PFFlavour+GenJetFlavour)
Revision:	1.5
Committed:	Wed Mar 11 18:17:10 2009 UTC (16 years, 1 month ago) by bendavid
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	Mit_008pre2
Changes since 1.4:	+108 -2 lines
Log Message:	Added PFJets flavour matching
#	User	Rev	Content
1	bendavid	1.5	# $Id: JetsMCFlavourMatching_cfi.py,v 1.4 2008/09/29 17:06:29 sixie 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			jets = cms.InputTag("iterativeCone5CaloJets"),
14			coneSizeToAssociate = cms.double(0.3),
15			partons = cms.InputTag("myPartons")
16			)
17
18			# Flavour byValue PhysDef
19			IC5byValPhys = cms.EDFilter("JetFlavourIdentifier",
20			srcByReference = cms.InputTag("IC5byRef"),
21			physicsDefinition = cms.bool(True)
22			)
23
24			# Flavour byValue AlgoDef
25			IC5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
26			srcByReference = cms.InputTag("IC5byRef"),
27			physicsDefinition = cms.bool(False)
28			)
29
30			#SisCone Jets
31			# Flavour byReference
32			SC5byRef = cms.EDFilter("JetPartonMatcher",
33			jets = cms.InputTag("sisCone5CaloJets"),
34			coneSizeToAssociate = cms.double(0.3),
35			partons = cms.InputTag("myPartons")
36			)
37
38			# Flavour byValue PhysDef
39			SC5byValPhys = cms.EDFilter("JetFlavourIdentifier",
40			srcByReference = cms.InputTag("SC5byRef"),
41			physicsDefinition = cms.bool(True)
42			)
43
44			# Flavour byValue AlgoDef
45			SC5byValAlgo = cms.EDFilter("JetFlavourIdentifier",
46			srcByReference = cms.InputTag("SC5byRef"),
47			physicsDefinition = cms.bool(False)
48			)
49
50			# Flavour byReference
51			SC7byRef = cms.EDFilter("JetPartonMatcher",
52			jets = cms.InputTag("sisCone7CaloJets"),
53			coneSizeToAssociate = cms.double(0.3),
54			partons = cms.InputTag("myPartons")
55			)
56
57			# Flavour byValue PhysDef
58			SC7byValPhys = cms.EDFilter("JetFlavourIdentifier",
59			srcByReference = cms.InputTag("SC7byRef"),
60			physicsDefinition = cms.bool(True)
61			)
62
63			# Flavour byValue AlgoDef
64			SC7byValAlgo = cms.EDFilter("JetFlavourIdentifier",
65			srcByReference = cms.InputTag("SC7byRef"),
66			physicsDefinition = cms.bool(False)
67			)
68
69			#Kt Jets
70			# Flavour byReference
71			KT4byRef = cms.EDFilter("JetPartonMatcher",
72			jets = cms.InputTag("kt4CaloJets"),
73			coneSizeToAssociate = cms.double(0.3),
74			partons = cms.InputTag("myPartons")
75			)
76
77			# Flavour byValue PhysDef
78			KT4byValPhys = cms.EDFilter("JetFlavourIdentifier",
79			srcByReference = cms.InputTag("KT4byRef"),
80			physicsDefinition = cms.bool(True)
81			)
82
83			# Flavour byValue AlgoDef
84			KT4byValAlgo = cms.EDFilter("JetFlavourIdentifier",
85			srcByReference = cms.InputTag("KT4byRef"),
86			physicsDefinition = cms.bool(False)
87			)
88
89
90			# Flavour byReference
91			KT6byRef = cms.EDFilter("JetPartonMatcher",
92			jets = cms.InputTag("kt6CaloJets"),
93			coneSizeToAssociate = cms.double(0.3),
94			partons = cms.InputTag("myPartons")
95			)
96
97			# Flavour byValue PhysDef
98			KT6byValPhys = cms.EDFilter("JetFlavourIdentifier",
99			srcByReference = cms.InputTag("KT6byRef"),
100			physicsDefinition = cms.bool(True)
101			)
102
103			# Flavour byValue AlgoDef
104			KT6byValAlgo = cms.EDFilter("JetFlavourIdentifier",
105			srcByReference = cms.InputTag("KT6byRef"),
106			physicsDefinition = cms.bool(False)
107			)
108
109	bendavid	1.5	#IC5PF Jets
110			# Flavour byReference
111			IC5PFbyRef = cms.EDFilter("JetPartonMatcher",
112			jets = cms.InputTag("iterativeCone5PFJets"),
113			coneSizeToAssociate = cms.double(0.3),
114			partons = cms.InputTag("myPartons")
115			)
116
117			# Flavour byValue PhysDef
118			IC5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
119			srcByReference = cms.InputTag("IC5PFbyRef"),
120			physicsDefinition = cms.bool(True)
121			)
122
123			# Flavour byValue AlgoDef
124			IC5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
125			srcByReference = cms.InputTag("IC5PFbyRef"),
126			physicsDefinition = cms.bool(False)
127			)
128
129			#SisCone Jets
130			# Flavour byReference
131			SC5PFbyRef = cms.EDFilter("JetPartonMatcher",
132			jets = cms.InputTag("sisCone5PFJets"),
133			coneSizeToAssociate = cms.double(0.3),
134			partons = cms.InputTag("myPartons")
135			)
136
137			# Flavour byValue PhysDef
138			SC5PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
139			srcByReference = cms.InputTag("SC5PFbyRef"),
140			physicsDefinition = cms.bool(True)
141			)
142
143			# Flavour byValue AlgoDef
144			SC5PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
145			srcByReference = cms.InputTag("SC5PFbyRef"),
146			physicsDefinition = cms.bool(False)
147			)
148
149			# Flavour byReference
150			SC7PFbyRef = cms.EDFilter("JetPartonMatcher",
151			jets = cms.InputTag("sisCone7PFJets"),
152			coneSizeToAssociate = cms.double(0.3),
153			partons = cms.InputTag("myPartons")
154			)
155
156			# Flavour byValue PhysDef
157			SC7PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
158			srcByReference = cms.InputTag("SC7PFbyRef"),
159			physicsDefinition = cms.bool(True)
160			)
161
162			# Flavour byValue AlgoDef
163			SC7PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
164			srcByReference = cms.InputTag("SC7PFbyRef"),
165			physicsDefinition = cms.bool(False)
166			)
167
168			#Kt Jets
169			# Flavour byReference
170			KT4PFbyRef = cms.EDFilter("JetPartonMatcher",
171			jets = cms.InputTag("kt4PFJets"),
172			coneSizeToAssociate = cms.double(0.3),
173			partons = cms.InputTag("myPartons")
174			)
175
176			# Flavour byValue PhysDef
177			KT4PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
178			srcByReference = cms.InputTag("KT4PFbyRef"),
179			physicsDefinition = cms.bool(True)
180			)
181
182			# Flavour byValue AlgoDef
183			KT4PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
184			srcByReference = cms.InputTag("KT4PFbyRef"),
185			physicsDefinition = cms.bool(False)
186			)
187
188
189			# Flavour byReference
190			KT6PFbyRef = cms.EDFilter("JetPartonMatcher",
191			jets = cms.InputTag("kt6PFJets"),
192			coneSizeToAssociate = cms.double(0.3),
193			partons = cms.InputTag("myPartons")
194			)
195
196			# Flavour byValue PhysDef
197			KT6PFbyValPhys = cms.EDFilter("JetFlavourIdentifier",
198			srcByReference = cms.InputTag("KT6PFbyRef"),
199			physicsDefinition = cms.bool(True)
200			)
201
202			# Flavour byValue AlgoDef
203			KT6PFbyValAlgo = cms.EDFilter("JetFlavourIdentifier",
204			srcByReference = cms.InputTag("KT6PFbyRef"),
205			physicsDefinition = cms.bool(False)
206			)
207
208
209	sixie	1.2	#Flavour by Reference
210			GenIC5byRef = cms.EDFilter("JetPartonMatcher",
211			jets = cms.InputTag("iterativeCone5GenJets"),
212			coneSizeToAssociate = cms.double(0.3),
213			partons = cms.InputTag("myPartons")
214			)
215
216			#Flavour by Reference
217			GenSC5byRef = cms.EDFilter("JetPartonMatcher",
218			jets = cms.InputTag("sisCone5GenJets"),
219			coneSizeToAssociate = cms.double(0.3),
220			partons = cms.InputTag("myPartons")
221			)
222	sixie	1.1
223	sixie	1.4	#Flavour by Reference
224			GenSC7byRef = cms.EDFilter("JetPartonMatcher",
225			jets = cms.InputTag("sisCone7GenJets"),
226			coneSizeToAssociate = cms.double(0.3),
227			partons = cms.InputTag("myPartons")
228			)
229
230			#Flavour by Reference
231			GenKT4byRef = cms.EDFilter("JetPartonMatcher",
232			jets = cms.InputTag("kt4GenJets"),
233			coneSizeToAssociate = cms.double(0.3),
234			partons = cms.InputTag("myPartons")
235			)
236
237			#Flavour by Reference
238			GenKT6byRef = cms.EDFilter("JetPartonMatcher",
239			jets = cms.InputTag("kt6GenJets"),
240			coneSizeToAssociate = cms.double(0.3),
241			partons = cms.InputTag("myPartons")
242			)
243
244	sixie	1.1	iterativeCone5Flavour = cms.Sequence(IC5byRefIC5byValPhysIC5byValAlgo)
245			sisCone5Flavour = cms.Sequence(SC5byRefSC5byValPhysSC5byValAlgo)
246			sisCone7Flavour = cms.Sequence(SC7byRefSC7byValPhysSC7byValAlgo)
247			KT4Flavour = cms.Sequence(KT4byRefKT4byValPhysKT4byValAlgo)
248			KT6Flavour = cms.Sequence(KT6byRefKT6byValPhysKT6byValAlgo)
249	bendavid	1.5	iterativeCone5PFFlavour = cms.Sequence(IC5PFbyRefIC5PFbyValPhysIC5PFbyValAlgo)
250			sisCone5PFFlavour = cms.Sequence(SC5PFbyRefSC5PFbyValPhysSC5PFbyValAlgo)
251			sisCone7PFFlavour = cms.Sequence(SC7PFbyRefSC7PFbyValPhysSC7PFbyValAlgo)
252			KT4PFFlavour = cms.Sequence(KT4PFbyRefKT4PFbyValPhysKT4PFbyValAlgo)
253			KT6PFFlavour = cms.Sequence(KT6PFbyRefKT6PFbyValPhysKT6PFbyValAlgo)
254	sixie	1.4	GenJetFlavour = cms.Sequence(GenIC5byRefGenSC5byRefGenSC7byRefGenKT4byRefGenKT6byRef)
255	bendavid	1.5	caloJetMCFlavour = cms.Sequence(myPartons*iterativeCone5Flavour+sisCone5Flavour+sisCone7Flavour+KT4Flavour+KT6Flavour
256			+iterativeCone5PFFlavour+sisCone5PFFlavour+sisCone7PFFlavour+KT4PFFlavour+KT6PFFlavour+GenJetFlavour)