TreeFiller/python/JetsMCFlavourMatching_cfi.py

# $Id: JetsMCFlavourMatching_cfi.py,v 1.3 2008/09/28 02:27:52 loizides 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)
)

#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)
GenJetFlavour = cms.Sequence(GenIC5byRef*GenSC5byRef*GenSC7byRef*GenKT4byRef*GenKT6byRef)
caloJetMCFlavour = cms.Sequence(myPartons*iterativeCone5Flavour+sisCone5Flavour+sisCone7Flavour+KT4Flavour+KT6Flavour+GenJetFlavour)
Revision:	1.4
Committed:	Mon Sep 29 17:06:29 2008 UTC (16 years, 7 months ago) by sixie
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	Mit_008pre1, Mit_006b, Mit_006a, Mit_006, Mit_005, Mit_004
Changes since 1.3:	+23 -2 lines
Log Message:	add more types of jet collections
#	User	Rev	Content
1	sixie	1.4	# $Id: JetsMCFlavourMatching_cfi.py,v 1.3 2008/09/28 02:27:52 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			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	sixie	1.2	#Flavour by Reference
110			GenIC5byRef = cms.EDFilter("JetPartonMatcher",
111			jets = cms.InputTag("iterativeCone5GenJets"),
112			coneSizeToAssociate = cms.double(0.3),
113			partons = cms.InputTag("myPartons")
114			)
115
116			#Flavour by Reference
117			GenSC5byRef = cms.EDFilter("JetPartonMatcher",
118			jets = cms.InputTag("sisCone5GenJets"),
119			coneSizeToAssociate = cms.double(0.3),
120			partons = cms.InputTag("myPartons")
121			)
122	sixie	1.1
123	sixie	1.4	#Flavour by Reference
124			GenSC7byRef = cms.EDFilter("JetPartonMatcher",
125			jets = cms.InputTag("sisCone7GenJets"),
126			coneSizeToAssociate = cms.double(0.3),
127			partons = cms.InputTag("myPartons")
128			)
129
130			#Flavour by Reference
131			GenKT4byRef = cms.EDFilter("JetPartonMatcher",
132			jets = cms.InputTag("kt4GenJets"),
133			coneSizeToAssociate = cms.double(0.3),
134			partons = cms.InputTag("myPartons")
135			)
136
137			#Flavour by Reference
138			GenKT6byRef = cms.EDFilter("JetPartonMatcher",
139			jets = cms.InputTag("kt6GenJets"),
140			coneSizeToAssociate = cms.double(0.3),
141			partons = cms.InputTag("myPartons")
142			)
143
144	sixie	1.1	iterativeCone5Flavour = cms.Sequence(IC5byRefIC5byValPhysIC5byValAlgo)
145			sisCone5Flavour = cms.Sequence(SC5byRefSC5byValPhysSC5byValAlgo)
146			sisCone7Flavour = cms.Sequence(SC7byRefSC7byValPhysSC7byValAlgo)
147			KT4Flavour = cms.Sequence(KT4byRefKT4byValPhysKT4byValAlgo)
148			KT6Flavour = cms.Sequence(KT6byRefKT6byValPhysKT6byValAlgo)
149	sixie	1.4	GenJetFlavour = cms.Sequence(GenIC5byRefGenSC5byRefGenSC7byRefGenKT4byRefGenKT6byRef)
150	sixie	1.2	caloJetMCFlavour = cms.Sequence(myPartons*iterativeCone5Flavour+sisCone5Flavour+sisCone7Flavour+KT4Flavour+KT6Flavour+GenJetFlavour)