1 |
frankma |
1.1 |
import FWCore.ParameterSet.Config as cms
|
2 |
|
|
|
3 |
|
|
# Pat Jets
|
4 |
|
|
from PhysicsTools.PatAlgos.patHeavyIonSequences_cff import *
|
5 |
|
|
patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets")
|
6 |
|
|
patJets.addBTagInfo = False
|
7 |
|
|
patJets.addTagInfos = False
|
8 |
|
|
patJets.addDiscriminators = False
|
9 |
|
|
patJets.addAssociatedTracks = False
|
10 |
|
|
patJets.addJetCharge = False
|
11 |
|
|
patJets.addJetID = False
|
12 |
|
|
patJets.getJetMCFlavour = False
|
13 |
|
|
patJets.addGenPartonMatch = True
|
14 |
|
|
patJets.addGenJetMatch = True
|
15 |
|
|
patJets.embedGenJetMatch = True
|
16 |
|
|
patJets.embedGenPartonMatch = True
|
17 |
|
|
patJets.embedCaloTowers = False
|
18 |
|
|
|
19 |
mnguyen |
1.8 |
|
20 |
yilmaz |
1.5 |
patJetCorrFactors.useNPV = False
|
21 |
frankma |
1.1 |
# full reco
|
22 |
|
|
icPu5corr = patJetCorrFactors.clone(
|
23 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
24 |
|
|
levels = cms.vstring('L2Relative','L3Absolute'),
|
25 |
yilmaz |
1.7 |
payload = cms.string('IC5Calo_2760GeV')
|
26 |
frankma |
1.1 |
)
|
27 |
yilmaz |
1.6 |
|
28 |
mnguyen |
1.8 |
#icPu5patJets = patJets.clone(
|
29 |
|
|
# jetSource = cms.InputTag("iterativeConePu5CaloJets"),
|
30 |
|
|
# genJetMatch = cms.InputTag("icPu5match"),
|
31 |
|
|
# genPartonMatch = cms.InputTag("icPu5parton"),
|
32 |
|
|
# jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr"))
|
33 |
|
|
# )
|
34 |
frankma |
1.1 |
|
35 |
|
|
akPu5PFcorr = icPu5corr.clone(
|
36 |
|
|
src = cms.InputTag("akPu5PFJets"),
|
37 |
yilmaz |
1.7 |
payload = cms.string('AK5PFTowers_hiGoodTightTracks')
|
38 |
frankma |
1.1 |
)
|
39 |
|
|
akPu5PFpatJets = patJets.clone(
|
40 |
|
|
jetSource = cms.InputTag("akPu5PFJets"),
|
41 |
|
|
genJetMatch = cms.InputTag("akPu5PFmatch"),
|
42 |
|
|
genPartonMatch = cms.InputTag("akPu5PFparton"),
|
43 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr"))
|
44 |
|
|
)
|
45 |
|
|
|
46 |
|
|
akPu3PFcorr = icPu5corr.clone(
|
47 |
|
|
src = cms.InputTag("akPu3PFJets"),
|
48 |
yilmaz |
1.7 |
payload = cms.string('AK3PFTowers_hiGoodTightTracks')
|
49 |
frankma |
1.1 |
)
|
50 |
yilmaz |
1.6 |
|
51 |
mnguyen |
1.8 |
#akPu3PFpatJets = patJets.clone(
|
52 |
|
|
# jetSource = cms.InputTag("akPu3PFJets"),
|
53 |
|
|
# genJetMatch = cms.InputTag("akPu3PFmatch"),
|
54 |
|
|
# genPartonMatch = cms.InputTag("akPu3PFparton"),
|
55 |
|
|
# jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr"))
|
56 |
|
|
# )
|
57 |
frankma |
1.1 |
|
58 |
yilmaz |
1.6 |
akPu5corr = icPu5corr.clone(
|
59 |
|
|
src = cms.InputTag("akPu5CaloJets"),
|
60 |
|
|
payload = cms.string('AK5Calo')
|
61 |
|
|
)
|
62 |
|
|
|
63 |
|
|
akPu5patJets = patJets.clone(
|
64 |
|
|
jetSource = cms.InputTag("akPu5CaloJets"),
|
65 |
|
|
genJetMatch = cms.InputTag("akPu5match"),
|
66 |
|
|
genPartonMatch = cms.InputTag("akPu5parton"),
|
67 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr"))
|
68 |
|
|
)
|
69 |
|
|
|
70 |
|
|
akPu3corr = icPu5corr.clone(
|
71 |
|
|
src = cms.InputTag("akPu3CaloJets"),
|
72 |
yilmaz |
1.7 |
payload = cms.string('AK5Calo')
|
73 |
yilmaz |
1.6 |
)
|
74 |
|
|
|
75 |
|
|
akPu3patJets = patJets.clone(
|
76 |
|
|
jetSource = cms.InputTag("akPu3CaloJets"),
|
77 |
|
|
genJetMatch = cms.InputTag("akPu3match"),
|
78 |
|
|
genPartonMatch = cms.InputTag("akPu3parton"),
|
79 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3corr"))
|
80 |
|
|
)
|
81 |
|
|
|
82 |
|
|
|
83 |
frankma |
1.1 |
|
84 |
|
|
# mc matching
|
85 |
|
|
patJetPartonMatch.matched = cms.InputTag("hiPartons")
|
86 |
|
|
|
87 |
|
|
icPu5clean = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets') )
|
88 |
|
|
icPu5match = patJetGenJetMatch.clone(
|
89 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets"),
|
90 |
|
|
matched = cms.InputTag("icPu5clean")
|
91 |
|
|
)
|
92 |
|
|
icPu5parton = patJetPartonMatch.clone(
|
93 |
|
|
src = cms.InputTag("iterativeConePu5CaloJets")
|
94 |
|
|
)
|
95 |
|
|
|
96 |
yilmaz |
1.3 |
|
97 |
frankma |
1.1 |
akPu5PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets') )
|
98 |
|
|
akPu5PFmatch = patJetGenJetMatch.clone(
|
99 |
|
|
src = cms.InputTag("akPu5PFJets"),
|
100 |
|
|
matched = cms.InputTag("akPu5PFclean")
|
101 |
|
|
)
|
102 |
|
|
akPu5PFparton = patJetPartonMatch.clone(
|
103 |
|
|
src = cms.InputTag("akPu5PFJets")
|
104 |
|
|
)
|
105 |
|
|
|
106 |
|
|
akPu3PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets') )
|
107 |
|
|
akPu3PFmatch = patJetGenJetMatch.clone(
|
108 |
|
|
src = cms.InputTag("akPu3PFJets"),
|
109 |
|
|
matched = cms.InputTag("akPu3PFclean")
|
110 |
|
|
)
|
111 |
|
|
akPu3PFparton = patJetPartonMatch.clone(
|
112 |
|
|
src = cms.InputTag("akPu3PFJets")
|
113 |
|
|
)
|
114 |
|
|
|
115 |
yilmaz |
1.6 |
akPu5match = patJetGenJetMatch.clone(
|
116 |
|
|
src = cms.InputTag("akPu5CaloJets"),
|
117 |
|
|
matched = cms.InputTag("akPu5PFclean")
|
118 |
|
|
)
|
119 |
|
|
akPu5parton = patJetPartonMatch.clone(
|
120 |
|
|
src = cms.InputTag("akPu5CaloJets")
|
121 |
|
|
)
|
122 |
|
|
|
123 |
|
|
akPu3match = patJetGenJetMatch.clone(
|
124 |
|
|
src = cms.InputTag("akPu3CaloJets"),
|
125 |
|
|
matched = cms.InputTag("akPu3PFclean")
|
126 |
|
|
)
|
127 |
|
|
akPu3parton = patJetPartonMatch.clone(
|
128 |
|
|
src = cms.InputTag("akPu3CaloJets")
|
129 |
|
|
)
|
130 |
|
|
|
131 |
mnguyen |
1.8 |
from RecoJets.JetAssociationProducers.ak5JTA_cff import *
|
132 |
|
|
from RecoBTag.Configuration.RecoBTag_cff import *
|
133 |
|
|
# @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ calo jet
|
134 |
|
|
#### b-tagging:
|
135 |
|
|
# b-tagging general configuration
|
136 |
|
|
# @@@@@@@@@ calo jet
|
137 |
|
|
#### b-tagging:
|
138 |
|
|
# b-tagging general configuration
|
139 |
|
|
icPu5JetTracksAssociatorAtVertex = ak5JetTracksAssociatorAtVertex.clone()
|
140 |
|
|
icPu5JetTracksAssociatorAtVertex.jets = cms.InputTag("iterativeConePu5CaloJets")
|
141 |
|
|
icPu5JetTracksAssociatorAtVertex.tracks = cms.InputTag("hiSelectedTracks")
|
142 |
|
|
|
143 |
|
|
# impact parameter b-tag
|
144 |
|
|
icPu5ImpactParameterTagInfos = impactParameterTagInfos.clone()
|
145 |
|
|
icPu5ImpactParameterTagInfos.jetTracks = cms.InputTag("icPu5JetTracksAssociatorAtVertex")
|
146 |
|
|
icPu5ImpactParameterTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
|
147 |
|
|
|
148 |
|
|
icPu5TrackCountingHighEffBJetTags = trackCountingHighEffBJetTags.clone()
|
149 |
|
|
icPu5TrackCountingHighEffBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
|
150 |
|
|
icPu5TrackCountingHighPurBJetTags = trackCountingHighPurBJetTags.clone()
|
151 |
|
|
icPu5TrackCountingHighPurBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
|
152 |
|
|
icPu5JetProbabilityBJetTags = jetProbabilityBJetTags.clone()
|
153 |
|
|
icPu5JetProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
|
154 |
|
|
icPu5JetBProbabilityBJetTags = jetBProbabilityBJetTags.clone()
|
155 |
|
|
icPu5JetBProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
|
156 |
|
|
|
157 |
|
|
# secondary vertex b-tag
|
158 |
|
|
icPu5SecondaryVertexTagInfos = secondaryVertexTagInfos.clone()
|
159 |
|
|
icPu5SecondaryVertexTagInfos.trackIPTagInfos = cms.InputTag("icPu5ImpactParameterTagInfos")
|
160 |
|
|
icPu5SimpleSecondaryVertexBJetTags = simpleSecondaryVertexBJetTags.clone()
|
161 |
|
|
icPu5SimpleSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SecondaryVertexTagInfos"))
|
162 |
|
|
icPu5CombinedSecondaryVertexBJetTags = combinedSecondaryVertexBJetTags.clone()
|
163 |
|
|
icPu5CombinedSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
|
164 |
|
|
cms.InputTag("icPu5SecondaryVertexTagInfos"))
|
165 |
|
|
icPu5CombinedSecondaryVertexMVABJetTags = combinedSecondaryVertexMVABJetTags.clone()
|
166 |
|
|
icPu5CombinedSecondaryVertexMVABJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
|
167 |
|
|
cms.InputTag("icPu5SecondaryVertexTagInfos"))
|
168 |
|
|
|
169 |
|
|
# soft muon b-tag
|
170 |
|
|
icPu5SoftMuonTagInfos = softMuonTagInfos.clone()
|
171 |
|
|
icPu5SoftMuonTagInfos.jets = cms.InputTag("iterativeConePu5CaloJets")
|
172 |
|
|
icPu5SoftMuonTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
|
173 |
|
|
icPu5SoftMuonBJetTags = softMuonBJetTags.clone()
|
174 |
|
|
icPu5SoftMuonBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
|
175 |
|
|
icPu5SoftMuonByIP3dBJetTags = softMuonByIP3dBJetTags.clone()
|
176 |
|
|
icPu5SoftMuonByIP3dBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
|
177 |
|
|
icPu5SoftMuonByPtBJetTags = softMuonByPtBJetTags.clone()
|
178 |
|
|
icPu5SoftMuonByPtBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
|
179 |
|
|
|
180 |
|
|
# ghost tracks
|
181 |
|
|
icPu5GhostTrackVertexTagInfos = ghostTrackVertexTagInfos.clone()
|
182 |
|
|
icPu5GhostTrackVertexTagInfos.trackIPTagInfos = cms.InputTag("icPu5ImpactParameterTagInfos")
|
183 |
|
|
icPu5GhostTrackBJetTags = ghostTrackBJetTags.clone()
|
184 |
|
|
icPu5GhostTrackBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
|
185 |
|
|
cms.InputTag("icPu5GhostTrackVertexTagInfos"))
|
186 |
|
|
# prepare a path running the new modules
|
187 |
|
|
icPu5JetTracksAssociator = cms.Sequence(icPu5JetTracksAssociatorAtVertex)
|
188 |
|
|
icPu5JetBtaggingIP = cms.Sequence(icPu5ImpactParameterTagInfos * (icPu5TrackCountingHighEffBJetTags +
|
189 |
|
|
icPu5TrackCountingHighPurBJetTags +
|
190 |
|
|
icPu5JetProbabilityBJetTags +
|
191 |
|
|
icPu5JetBProbabilityBJetTags
|
192 |
|
|
)
|
193 |
|
|
)
|
194 |
|
|
|
195 |
|
|
icPu5JetBtaggingSV = cms.Sequence(icPu5ImpactParameterTagInfos *
|
196 |
|
|
icPu5SecondaryVertexTagInfos * (icPu5SimpleSecondaryVertexBJetTags +
|
197 |
|
|
icPu5CombinedSecondaryVertexBJetTags +
|
198 |
|
|
icPu5CombinedSecondaryVertexMVABJetTags
|
199 |
|
|
)
|
200 |
|
|
+icPu5GhostTrackVertexTagInfos
|
201 |
|
|
*icPu5GhostTrackBJetTags
|
202 |
|
|
)
|
203 |
|
|
|
204 |
|
|
|
205 |
|
|
icPu5JetBtaggingMu = cms.Sequence(icPu5SoftMuonTagInfos * (icPu5SoftMuonBJetTags +
|
206 |
|
|
icPu5SoftMuonByIP3dBJetTags +
|
207 |
|
|
icPu5SoftMuonByPtBJetTags
|
208 |
|
|
)
|
209 |
|
|
)
|
210 |
|
|
|
211 |
|
|
icPu5JetBtagging = cms.Sequence(icPu5JetBtaggingIP
|
212 |
|
|
*icPu5JetBtaggingSV
|
213 |
|
|
*icPu5JetBtaggingMu
|
214 |
|
|
)
|
215 |
|
|
#----------------------
|
216 |
|
|
|
217 |
|
|
|
218 |
|
|
#----------------------
|
219 |
|
|
icPu5corr = patJetCorrFactors.clone(src = cms.InputTag("iterativeConePu5CaloJets"),
|
220 |
|
|
levels = cms.vstring('L2Relative','L3Absolute'),
|
221 |
|
|
payload = cms.string('IC5Calo_2760GeV'))
|
222 |
|
|
icPu5clean = heavyIonCleanedGenJets.clone(src = cms.InputTag('iterativeCone5HiGenJets')) # cleans the jets, but NOT the partons
|
223 |
|
|
icPu5match = patJetGenJetMatch.clone(src = cms.InputTag("iterativeConePu5CaloJets"),
|
224 |
|
|
matched = cms.InputTag("icPu5clean"))
|
225 |
|
|
|
226 |
|
|
icPu5parton = patJetPartonMatch.clone(src = cms.InputTag("iterativeConePu5CaloJets"))
|
227 |
|
|
|
228 |
|
|
# ----- flavour bit
|
229 |
|
|
icPu5PatJetPartonAssociation = patJetPartonAssociation.clone(jets = cms.InputTag("iterativeConePu5CaloJets"),
|
230 |
|
|
partons = cms.InputTag("genPartons"),
|
231 |
|
|
coneSizeToAssociate = cms.double(0.4))
|
232 |
|
|
icPu5PatJetFlavourAssociation = patJetFlavourAssociation.clone(srcByReference = cms.InputTag("icPu5PatJetPartonAssociation"))
|
233 |
|
|
|
234 |
|
|
icPu5PatJetFlavourId = cms.Sequence(icPu5PatJetPartonAssociation*icPu5PatJetFlavourAssociation)
|
235 |
|
|
|
236 |
|
|
#-------
|
237 |
|
|
|
238 |
|
|
icPu5patJets = patJets.clone(jetSource = cms.InputTag("iterativeConePu5CaloJets"),
|
239 |
|
|
genJetMatch = cms.InputTag("icPu5match"),
|
240 |
|
|
genPartonMatch = cms.InputTag("icPu5parton"),
|
241 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")),
|
242 |
|
|
JetPartonMapSource = cms.InputTag("icPu5PatJetFlavourAssociation"),
|
243 |
|
|
trackAssociationSource = cms.InputTag("icPu5JetTracksAssociatorAtVertex"),
|
244 |
|
|
discriminatorSources = cms.VInputTag(cms.InputTag("icPu5CombinedSecondaryVertexBJetTags"),
|
245 |
|
|
cms.InputTag("icPu5CombinedSecondaryVertexMVABJetTags"),
|
246 |
|
|
cms.InputTag("icPu5JetBProbabilityBJetTags"),
|
247 |
|
|
cms.InputTag("icPu5JetProbabilityBJetTags"),
|
248 |
|
|
cms.InputTag("icPu5SoftMuonByPtBJetTags"),
|
249 |
|
|
cms.InputTag("icPu5SoftMuonByIP3dBJetTags"),
|
250 |
|
|
cms.InputTag("icPu5TrackCountingHighEffBJetTags"),
|
251 |
|
|
cms.InputTag("icPu5TrackCountingHighPurBJetTags"),
|
252 |
|
|
),
|
253 |
|
|
)
|
254 |
|
|
|
255 |
|
|
|
256 |
|
|
#### B-tagging for this bit:
|
257 |
|
|
# b-tagging general configuration
|
258 |
|
|
akPu3PFJetTracksAssociatorAtVertex = ak5JetTracksAssociatorAtVertex.clone()
|
259 |
|
|
akPu3PFJetTracksAssociatorAtVertex.jets = cms.InputTag("akPu3PFJets")
|
260 |
|
|
akPu3PFJetTracksAssociatorAtVertex.tracks = cms.InputTag("hiSelectedTracks")
|
261 |
|
|
# impact parameter b-tag
|
262 |
|
|
akPu3PFImpactParameterTagInfos = impactParameterTagInfos.clone()
|
263 |
|
|
akPu3PFImpactParameterTagInfos.jetTracks = cms.InputTag("akPu3PFJetTracksAssociatorAtVertex")
|
264 |
|
|
akPu3PFImpactParameterTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
|
265 |
|
|
akPu3PFTrackCountingHighEffBJetTags = trackCountingHighEffBJetTags.clone()
|
266 |
|
|
akPu3PFTrackCountingHighEffBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
|
267 |
|
|
akPu3PFTrackCountingHighPurBJetTags = trackCountingHighPurBJetTags.clone()
|
268 |
|
|
akPu3PFTrackCountingHighPurBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
|
269 |
|
|
akPu3PFJetProbabilityBJetTags = jetProbabilityBJetTags.clone()
|
270 |
|
|
akPu3PFJetProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
|
271 |
|
|
akPu3PFJetBProbabilityBJetTags = jetBProbabilityBJetTags.clone()
|
272 |
|
|
akPu3PFJetBProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
|
273 |
|
|
|
274 |
|
|
# secondary vertex b-tag
|
275 |
|
|
akPu3PFSecondaryVertexTagInfos = secondaryVertexTagInfos.clone()
|
276 |
|
|
akPu3PFSecondaryVertexTagInfos.trackIPTagInfos = cms.InputTag("akPu3PFImpactParameterTagInfos")
|
277 |
|
|
akPu3PFSimpleSecondaryVertexBJetTags = simpleSecondaryVertexBJetTags.clone()
|
278 |
|
|
akPu3PFSimpleSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
|
279 |
|
|
akPu3PFCombinedSecondaryVertexBJetTags = combinedSecondaryVertexBJetTags.clone()
|
280 |
|
|
akPu3PFCombinedSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
|
281 |
|
|
cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
|
282 |
|
|
akPu3PFCombinedSecondaryVertexMVABJetTags = combinedSecondaryVertexMVABJetTags.clone()
|
283 |
|
|
akPu3PFCombinedSecondaryVertexMVABJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
|
284 |
|
|
cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
|
285 |
|
|
|
286 |
|
|
# soft muon b-tag
|
287 |
|
|
akPu3PFSoftMuonTagInfos = softMuonTagInfos.clone()
|
288 |
|
|
akPu3PFSoftMuonTagInfos.jets = cms.InputTag("akPu3PFJets")
|
289 |
|
|
akPu3PFSoftMuonTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
|
290 |
|
|
akPu3PFSoftMuonBJetTags = softMuonBJetTags.clone()
|
291 |
|
|
akPu3PFSoftMuonBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
|
292 |
|
|
akPu3PFSoftMuonByIP3dBJetTags = softMuonByIP3dBJetTags.clone()
|
293 |
|
|
akPu3PFSoftMuonByIP3dBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
|
294 |
|
|
akPu3PFSoftMuonByPtBJetTags = softMuonByPtBJetTags.clone()
|
295 |
|
|
akPu3PFSoftMuonByPtBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
|
296 |
|
|
|
297 |
|
|
# ghost tracks
|
298 |
|
|
akPu3PFGhostTrackVertexTagInfos = ghostTrackVertexTagInfos.clone()
|
299 |
|
|
akPu3PFGhostTrackVertexTagInfos.trackIPTagInfos = cms.InputTag("akPu3PFImpactParameterTagInfos")
|
300 |
|
|
akPu3PFGhostTrackBJetTags = ghostTrackBJetTags.clone()
|
301 |
|
|
akPu3PFGhostTrackBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
|
302 |
|
|
cms.InputTag("akPu3PFGhostTrackVertexTagInfos"))
|
303 |
|
|
# prepare a path running the new modules
|
304 |
|
|
akPu3PFJetTracksAssociator = cms.Sequence(akPu3PFJetTracksAssociatorAtVertex)
|
305 |
|
|
akPu3PFJetBtaggingIP = cms.Sequence(akPu3PFImpactParameterTagInfos * (akPu3PFTrackCountingHighEffBJetTags +
|
306 |
|
|
akPu3PFTrackCountingHighPurBJetTags +
|
307 |
|
|
akPu3PFJetProbabilityBJetTags +
|
308 |
|
|
akPu3PFJetBProbabilityBJetTags
|
309 |
|
|
)
|
310 |
|
|
)
|
311 |
|
|
|
312 |
|
|
akPu3PFJetBtaggingSV = cms.Sequence(akPu3PFImpactParameterTagInfos *
|
313 |
|
|
akPu3PFSecondaryVertexTagInfos * (akPu3PFSimpleSecondaryVertexBJetTags +
|
314 |
|
|
akPu3PFCombinedSecondaryVertexBJetTags +
|
315 |
|
|
akPu3PFCombinedSecondaryVertexMVABJetTags
|
316 |
|
|
)
|
317 |
|
|
+akPu3PFGhostTrackVertexTagInfos
|
318 |
|
|
*akPu3PFGhostTrackBJetTags
|
319 |
|
|
)
|
320 |
|
|
|
321 |
|
|
|
322 |
|
|
akPu3PFJetBtaggingMu = cms.Sequence(akPu3PFSoftMuonTagInfos * (akPu3PFSoftMuonBJetTags +
|
323 |
|
|
akPu3PFSoftMuonByIP3dBJetTags +
|
324 |
|
|
akPu3PFSoftMuonByPtBJetTags
|
325 |
|
|
)
|
326 |
|
|
)
|
327 |
|
|
|
328 |
|
|
akPu3PFJetBtagging = cms.Sequence(akPu3PFJetBtaggingIP
|
329 |
|
|
*akPu3PFJetBtaggingSV
|
330 |
|
|
*akPu3PFJetBtaggingMu
|
331 |
|
|
)
|
332 |
|
|
|
333 |
|
|
#__________________________________________________________
|
334 |
|
|
# ----- flavour bit
|
335 |
|
|
akPu3PFPatJetPartonAssociation = patJetPartonAssociation.clone(jets = cms.InputTag("akPu3PFJets"),
|
336 |
|
|
partons = cms.InputTag("genPartons"),
|
337 |
|
|
coneSizeToAssociate = cms.double(0.4))
|
338 |
|
|
akPu3PFPatJetFlavourAssociation = patJetFlavourAssociation.clone(srcByReference = cms.InputTag("akPu3PFPatJetPartonAssociation"))
|
339 |
|
|
|
340 |
|
|
akPu3PFPatJetFlavourId = cms.Sequence(akPu3PFPatJetPartonAssociation*akPu3PFPatJetFlavourAssociation)
|
341 |
|
|
#
|
342 |
|
|
#-------
|
343 |
|
|
akPu3PFcorr = patJetCorrFactors.clone(src = cms.InputTag("akPu3PFJets"),
|
344 |
|
|
levels = cms.vstring('L2Relative','L3Absolute'),
|
345 |
|
|
payload = cms.string('AK3PFTowers_hiGoodTightTracks')
|
346 |
|
|
)
|
347 |
|
|
akPu3PFclean = heavyIonCleanedGenJets.clone(src = cms.InputTag('ak3HiGenJets'))
|
348 |
|
|
akPu3PFmatch = patJetGenJetMatch.clone(src = cms.InputTag("akPu3PFJets"),
|
349 |
|
|
matched = cms.InputTag("akPu3PFclean"))
|
350 |
|
|
akPu3PFparton = patJetPartonMatch.clone(src = cms.InputTag("akPu3PFJets"))
|
351 |
|
|
akPu3PFpatJets = patJets.clone(jetSource = cms.InputTag("akPu3PFJets"),
|
352 |
|
|
genJetMatch = cms.InputTag("akPu3PFmatch"),
|
353 |
|
|
genPartonMatch = cms.InputTag("akPu3PFparton"),
|
354 |
|
|
jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr")),
|
355 |
|
|
JetPartonMapSource = cms.InputTag("akPu3PFPatJetFlavourAssociation"),
|
356 |
|
|
trackAssociationSource = cms.InputTag("akPu3PFJetTracksAssociatorAtVertex"),
|
357 |
|
|
discriminatorSources = cms.VInputTag(cms.InputTag("akPu3PFCombinedSecondaryVertexBJetTags"),
|
358 |
|
|
cms.InputTag("akPu3PFCombinedSecondaryVertexMVABJetTags"),
|
359 |
|
|
cms.InputTag("akPu3PFJetBProbabilityBJetTags"),
|
360 |
|
|
cms.InputTag("akPu3PFJetProbabilityBJetTags"),
|
361 |
|
|
cms.InputTag("akPu3PFSoftMuonByPtBJetTags"),
|
362 |
|
|
cms.InputTag("akPu3PFSoftMuonByIP3dBJetTags"),
|
363 |
|
|
cms.InputTag("akPu3PFTrackCountingHighEffBJetTags"),
|
364 |
|
|
cms.InputTag("akPu3PFTrackCountingHighPurBJetTags"),
|
365 |
|
|
),
|
366 |
|
|
)
|
367 |
|
|
|
368 |
|
|
# end batagging
|
369 |
|
|
|
370 |
|
|
|
371 |
frankma |
1.1 |
# === mc sequences ===
|
372 |
mnguyen |
1.8 |
icPu5patSequence = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5patJets)
|
373 |
frankma |
1.1 |
akPu5PFpatSequence = cms.Sequence(akPu5PFcorr * akPu5PFclean * akPu5PFmatch * akPu5PFparton * akPu5PFpatJets)
|
374 |
|
|
akPu3PFpatSequence = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFpatJets)
|
375 |
yilmaz |
1.6 |
akPu5patSequence = cms.Sequence(akPu5corr * akPu5PFclean * akPu5match * akPu5parton * akPu5patJets)
|
376 |
|
|
akPu3patSequence = cms.Sequence(akPu3corr * akPu3PFclean * akPu3match * akPu3parton * akPu3patJets)
|
377 |
mnguyen |
1.8 |
|
378 |
|
|
icPu5patSequence_withBtagging = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5PatJetFlavourId * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5patJets)
|
379 |
|
|
akPu3PFpatSequence_withBtagging = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFPatJetFlavourId * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFpatJets)
|
380 |
|
|
|
381 |
|
|
|
382 |
|
|
# === data sequences ===
|
383 |
|
|
# Note still need to use enableData function in cfg to remove mc dep of patjet
|
384 |
|
|
icPu5patSequence_data = cms.Sequence( icPu5corr * icPu5patJets )
|
385 |
|
|
akPu5PFpatSequence_data = cms.Sequence(akPu5PFcorr * akPu5PFpatJets )
|
386 |
|
|
akPu3PFpatSequence_data = cms.Sequence( akPu3PFcorr * akPu3PFpatJets )
|
387 |
|
|
akPu5patSequence_data = cms.Sequence( akPu5corr * akPu5patJets )
|
388 |
|
|
akPu3patSequence_data = cms.Sequence( akPu3corr * akPu3patJets )
|
389 |
|
|
|
390 |
|
|
icPu5patSequence_withBtagging_data = cms.Sequence( icPu5corr * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5patJets )
|
391 |
|
|
akPu3PFpatSequence_withBtagging_data = cms.Sequence( akPu3PFcorr * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFpatJets )
|
392 |
yilmaz |
1.9 |
|
393 |
|
|
|
394 |
|
|
# AK Algorithms
|
395 |
|
|
|
396 |
|
|
|
397 |
|
|
akPu6PFJets = akPu3PFJets.clone(rParam = 0.6)
|
398 |
|
|
akPu6CaloJets = akPu3CaloJets.clone(rParam = 0.6)
|
399 |
|
|
akPu4PFJets = akPu3PFJets.clone(rParam = 0.4)
|
400 |
|
|
akPu4CaloJets = akPu3CaloJets.clone(rParam = 0.4)
|
401 |
|
|
akPu2PFJets = akPu3PFJets.clone(rParam = 0.2)
|
402 |
|
|
akPu2CaloJets = akPu3CaloJets.clone(rParam = 0.2)
|
403 |
|
|
akPu1PFJets = akPu3PFJets.clone(rParam = 0.1)
|
404 |
|
|
akPu1CaloJets = akPu3CaloJets.clone(rParam = 0.1)
|
405 |
|
|
|
406 |
|
|
iterativeCone5CaloJets = iterativeConePu5CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
407 |
|
|
ak6PFJets = akPu6PFJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
408 |
|
|
ak6CaloJets = akPu6CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
409 |
|
|
ak5PFJets = akPu5PFJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
410 |
|
|
ak5CaloJets = akPu5CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
411 |
|
|
ak4PFJets = akPu4PFJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
412 |
|
|
ak4CaloJets = akPu4CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
413 |
|
|
ak3PFJets = akPu3PFJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
414 |
|
|
ak3CaloJets = akPu3CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
415 |
|
|
ak2PFJets = akPu2PFJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
416 |
|
|
ak2CaloJets = akPu2CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
417 |
|
|
ak1PFJets = akPu1PFJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
418 |
|
|
ak1CaloJets = akPu1CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 1)
|
419 |
|
|
|
420 |
|
|
|
421 |
yilmaz |
1.10 |
recoAk1to6 = cms.Sequence( akPu1PFJets * akPu2PFJets *akPu3PFJets * akPu4PFJets * akPu5PFJets * akPu6PFJets *
|
422 |
yilmaz |
1.9 |
ak1PFJets * ak2PFJets *ak3PFJets * ak4PFJets * ak5PFJets * ak6PFJets *
|
423 |
|
|
akPu1CaloJets * akPu2CaloJets *akPu3CaloJets * akPu4CaloJets * akPu5CaloJets * akPu6CaloJets *
|
424 |
yilmaz |
1.11 |
ak1CaloJets * ak2CaloJets *ak3CaloJets * ak4CaloJets * ak5CaloJets * ak6CaloJets
|
425 |
yilmaz |
1.9 |
)
|
426 |
|
|
|
427 |
|
|
|
428 |
|
|
|
429 |
|
|
# All corrections
|
430 |
|
|
|
431 |
|
|
|
432 |
|
|
akPu1PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu1PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
433 |
|
|
akPu2PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu2PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
434 |
|
|
akPu3PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu3PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
435 |
|
|
akPu4PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu4PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
436 |
|
|
akPu5PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu5PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
437 |
|
|
akPu6PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu6PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
438 |
|
|
|
439 |
|
|
ak1PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak1PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
440 |
|
|
ak2PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak2PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
441 |
|
|
ak3PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak3PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
442 |
|
|
ak4PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak4PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
443 |
|
|
ak5PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak5PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
444 |
|
|
ak6PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak6PFJets"),payload = cms.string('AK3PFTowers_hiGoodTightTracks'))
|
445 |
|
|
|
446 |
|
|
akPu1Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu1CaloJets"),payload = cms.string('AK5Calo'))
|
447 |
|
|
akPu2Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu2CaloJets"),payload = cms.string('AK5Calo'))
|
448 |
|
|
akPu3Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu3CaloJets"),payload = cms.string('AK5Calo'))
|
449 |
|
|
akPu4Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu4CaloJets"),payload = cms.string('AK5Calo'))
|
450 |
|
|
akPu5Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu5CaloJets"),payload = cms.string('AK5Calo'))
|
451 |
|
|
akPu6Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu6CaloJets"),payload = cms.string('AK5Calo'))
|
452 |
|
|
|
453 |
|
|
ak1Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak1CaloJets"),payload = cms.string('AK5Calo'))
|
454 |
|
|
ak2Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak2CaloJets"),payload = cms.string('AK5Calo'))
|
455 |
|
|
ak3Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak3CaloJets"),payload = cms.string('AK5Calo'))
|
456 |
|
|
ak4Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak4CaloJets"),payload = cms.string('AK5Calo'))
|
457 |
|
|
ak5Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak5CaloJets"),payload = cms.string('AK5Calo'))
|
458 |
|
|
ak6Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak6CaloJets"),payload = cms.string('AK5Calo'))
|
459 |
|
|
|
460 |
|
|
# Gen stuff
|
461 |
|
|
|
462 |
|
|
ak1clean = akPu3PFclean.clone()
|
463 |
|
|
ak2clean = akPu3PFclean.clone()
|
464 |
|
|
ak3clean = akPu3PFclean.clone()
|
465 |
|
|
ak4clean = akPu3PFclean.clone()
|
466 |
|
|
ak5clean = akPu3PFclean.clone()
|
467 |
|
|
ak6clean = akPu3PFclean.clone()
|
468 |
|
|
|
469 |
|
|
|
470 |
|
|
akPu1PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu1PFJets"), matched = cms.InputTag("ak1clean"))
|
471 |
|
|
akPu2PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu2PFJets"), matched = cms.InputTag("ak2clean"))
|
472 |
|
|
akPu3PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu3PFJets"), matched = cms.InputTag("ak3clean"))
|
473 |
|
|
akPu4PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu4PFJets"), matched = cms.InputTag("ak4clean"))
|
474 |
|
|
akPu5PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu5PFJets"), matched = cms.InputTag("ak5clean"))
|
475 |
|
|
akPu6PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu6PFJets"), matched = cms.InputTag("ak6clean"))
|
476 |
|
|
akPu1Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu1CaloJets"), matched = cms.InputTag("ak1clean"))
|
477 |
|
|
akPu2Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu2CaloJets"), matched = cms.InputTag("ak2clean"))
|
478 |
|
|
akPu3Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu3CaloJets"), matched = cms.InputTag("ak3clean"))
|
479 |
|
|
akPu4Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu4CaloJets"), matched = cms.InputTag("ak4clean"))
|
480 |
|
|
akPu5Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu5CaloJets"), matched = cms.InputTag("ak5clean"))
|
481 |
|
|
akPu6Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu6CaloJets"), matched = cms.InputTag("ak6clean"))
|
482 |
|
|
ak1PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak1PFJets"), matched = cms.InputTag("ak1clean"))
|
483 |
|
|
ak2PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak2PFJets"), matched = cms.InputTag("ak2clean"))
|
484 |
|
|
ak3PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak3PFJets"), matched = cms.InputTag("ak3clean"))
|
485 |
|
|
ak4PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak4PFJets"), matched = cms.InputTag("ak4clean"))
|
486 |
|
|
ak5PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak5PFJets"), matched = cms.InputTag("ak5clean"))
|
487 |
|
|
ak6PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak6PFJets"), matched = cms.InputTag("ak6clean"))
|
488 |
|
|
ak1Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak1CaloJets"), matched = cms.InputTag("ak1clean"))
|
489 |
|
|
ak2Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak2CaloJets"), matched = cms.InputTag("ak2clean"))
|
490 |
|
|
ak3Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak3CaloJets"), matched = cms.InputTag("ak3clean"))
|
491 |
|
|
ak4Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak4CaloJets"), matched = cms.InputTag("ak4clean"))
|
492 |
|
|
ak5Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak5CaloJets"), matched = cms.InputTag("ak5clean"))
|
493 |
|
|
ak6Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak6CaloJets"), matched = cms.InputTag("ak6clean"))
|
494 |
|
|
akPu1PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu1PFJets"))
|
495 |
|
|
akPu2PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu2PFJets"))
|
496 |
|
|
akPu3PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu3PFJets"))
|
497 |
|
|
akPu4PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu4PFJets"))
|
498 |
|
|
akPu5PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu5PFJets"))
|
499 |
|
|
akPu6PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu6PFJets"))
|
500 |
|
|
akPu1Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu1CaloJets"))
|
501 |
|
|
akPu2Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu2CaloJets"))
|
502 |
|
|
akPu3Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu3CaloJets"))
|
503 |
|
|
akPu4Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu4CaloJets"))
|
504 |
|
|
akPu5Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu5CaloJets"))
|
505 |
|
|
akPu6Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu6CaloJets"))
|
506 |
|
|
ak1PFparton = akPu3PFparton.clone(src = cms.InputTag("ak1PFJets"))
|
507 |
|
|
ak2PFparton = akPu3PFparton.clone(src = cms.InputTag("ak2PFJets"))
|
508 |
|
|
ak3PFparton = akPu3PFparton.clone(src = cms.InputTag("ak3PFJets"))
|
509 |
|
|
ak4PFparton = akPu3PFparton.clone(src = cms.InputTag("ak4PFJets"))
|
510 |
|
|
ak5PFparton = akPu3PFparton.clone(src = cms.InputTag("ak5PFJets"))
|
511 |
|
|
ak6PFparton = akPu3PFparton.clone(src = cms.InputTag("ak6PFJets"))
|
512 |
|
|
ak1Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak1CaloJets"))
|
513 |
|
|
ak2Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak2CaloJets"))
|
514 |
|
|
ak3Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak3CaloJets"))
|
515 |
|
|
ak4Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak4CaloJets"))
|
516 |
|
|
ak5Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak5CaloJets"))
|
517 |
|
|
ak6Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak6CaloJets"))
|
518 |
|
|
|
519 |
|
|
|
520 |
|
|
# PAT Maker
|
521 |
|
|
|
522 |
|
|
akPu1PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu1PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu1PFcorr")), genJetMatch = cms.InputTag("akPu1PFmatch"), genPartonMatch = cms.InputTag("akPu1PFparton"))
|
523 |
|
|
akPu2PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu2PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu2PFcorr")), genJetMatch = cms.InputTag("akPu2PFmatch"), genPartonMatch = cms.InputTag("akPu2PFparton"))
|
524 |
|
|
akPu3PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu3PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr")), genJetMatch = cms.InputTag("akPu3PFmatch"), genPartonMatch = cms.InputTag("akPu3PFparton"))
|
525 |
|
|
akPu4PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu4PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu4PFcorr")), genJetMatch = cms.InputTag("akPu4PFmatch"), genPartonMatch = cms.InputTag("akPu4PFparton"))
|
526 |
|
|
akPu5PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu5PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr")), genJetMatch = cms.InputTag("akPu5PFmatch"), genPartonMatch = cms.InputTag("akPu5PFparton"))
|
527 |
|
|
akPu6PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu6PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu6PFcorr")), genJetMatch = cms.InputTag("akPu6PFmatch"), genPartonMatch = cms.InputTag("akPu6PFparton"))
|
528 |
|
|
akPu1CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu1CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu1Calocorr")), genJetMatch = cms.InputTag("akPu1Calomatch"), genPartonMatch = cms.InputTag("akPu1Caloparton"))
|
529 |
|
|
akPu2CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu2CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu2Calocorr")), genJetMatch = cms.InputTag("akPu2Calomatch"), genPartonMatch = cms.InputTag("akPu2Caloparton"))
|
530 |
|
|
akPu3CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu3CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3Calocorr")), genJetMatch = cms.InputTag("akPu3Calomatch"), genPartonMatch = cms.InputTag("akPu3Caloparton"))
|
531 |
|
|
akPu4CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu4CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu4Calocorr")), genJetMatch = cms.InputTag("akPu4Calomatch"), genPartonMatch = cms.InputTag("akPu4Caloparton"))
|
532 |
|
|
akPu5CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu5CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5Calocorr")), genJetMatch = cms.InputTag("akPu5Calomatch"), genPartonMatch = cms.InputTag("akPu5Caloparton"))
|
533 |
|
|
akPu6CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu6CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu6Calocorr")), genJetMatch = cms.InputTag("akPu6Calomatch"), genPartonMatch = cms.InputTag("akPu6Caloparton"))
|
534 |
|
|
ak1PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak1PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak1PFcorr")), genJetMatch = cms.InputTag("ak1PFmatch"), genPartonMatch = cms.InputTag("ak1PFparton"))
|
535 |
|
|
ak2PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak2PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak2PFcorr")), genJetMatch = cms.InputTag("ak2PFmatch"), genPartonMatch = cms.InputTag("ak2PFparton"))
|
536 |
|
|
ak3PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak3PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak3PFcorr")), genJetMatch = cms.InputTag("ak3PFmatch"), genPartonMatch = cms.InputTag("ak3PFparton"))
|
537 |
|
|
ak4PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak4PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak4PFcorr")), genJetMatch = cms.InputTag("ak4PFmatch"), genPartonMatch = cms.InputTag("ak4PFparton"))
|
538 |
|
|
ak5PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak5PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5PFcorr")), genJetMatch = cms.InputTag("ak5PFmatch"), genPartonMatch = cms.InputTag("ak5PFparton"))
|
539 |
|
|
ak6PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak6PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak6PFcorr")), genJetMatch = cms.InputTag("ak6PFmatch"), genPartonMatch = cms.InputTag("ak6PFparton"))
|
540 |
|
|
ak1CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak1CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak1Calocorr")), genJetMatch = cms.InputTag("ak1Calomatch"), genPartonMatch = cms.InputTag("ak1Caloparton"))
|
541 |
|
|
ak2CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak2CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak2Calocorr")), genJetMatch = cms.InputTag("ak2Calomatch"), genPartonMatch = cms.InputTag("ak2Caloparton"))
|
542 |
|
|
ak3CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak3CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak3Calocorr")), genJetMatch = cms.InputTag("ak3Calomatch"), genPartonMatch = cms.InputTag("ak3Caloparton"))
|
543 |
|
|
ak4CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak4CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak4Calocorr")), genJetMatch = cms.InputTag("ak4Calomatch"), genPartonMatch = cms.InputTag("ak4Caloparton"))
|
544 |
|
|
ak5CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak5CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5Calocorr")), genJetMatch = cms.InputTag("ak5Calomatch"), genPartonMatch = cms.InputTag("ak5Caloparton"))
|
545 |
|
|
ak6CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak6CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak6Calocorr")), genJetMatch = cms.InputTag("ak6Calomatch"), genPartonMatch = cms.InputTag("ak6Caloparton"))
|
546 |
|
|
|
547 |
|
|
|
548 |
|
|
akPu1PFpatSequence = cms.Sequence(akPu1PFcorr+ak1clean+akPu1PFmatch+akPu1PFparton+akPu1PFpatJets)
|
549 |
|
|
akPu2PFpatSequence = cms.Sequence(akPu2PFcorr+ak2clean+akPu2PFmatch+akPu2PFparton+akPu2PFpatJets)
|
550 |
|
|
akPu3PFpatSequence = cms.Sequence(akPu3PFcorr+ak3clean+akPu3PFmatch+akPu3PFparton+akPu3PFpatJets)
|
551 |
|
|
akPu4PFpatSequence = cms.Sequence(akPu4PFcorr+ak4clean+akPu4PFmatch+akPu4PFparton+akPu4PFpatJets)
|
552 |
|
|
akPu5PFpatSequence = cms.Sequence(akPu5PFcorr+ak5clean+akPu5PFmatch+akPu5PFparton+akPu5PFpatJets)
|
553 |
|
|
akPu6PFpatSequence = cms.Sequence(akPu6PFcorr+ak6clean+akPu6PFmatch+akPu6PFparton+akPu6PFpatJets)
|
554 |
|
|
|
555 |
|
|
akPu1CalopatSequence = cms.Sequence(akPu1Calocorr+ak1clean+akPu1Calomatch+akPu1Caloparton+akPu1CalopatJets)
|
556 |
|
|
akPu2CalopatSequence = cms.Sequence(akPu2Calocorr+ak2clean+akPu2Calomatch+akPu2Caloparton+akPu2CalopatJets)
|
557 |
|
|
akPu3CalopatSequence = cms.Sequence(akPu3Calocorr+ak3clean+akPu3Calomatch+akPu3Caloparton+akPu3CalopatJets)
|
558 |
|
|
akPu4CalopatSequence = cms.Sequence(akPu4Calocorr+ak4clean+akPu4Calomatch+akPu4Caloparton+akPu4CalopatJets)
|
559 |
|
|
akPu5CalopatSequence = cms.Sequence(akPu5Calocorr+ak5clean+akPu5Calomatch+akPu5Caloparton+akPu5CalopatJets)
|
560 |
|
|
akPu6CalopatSequence = cms.Sequence(akPu6Calocorr+ak6clean+akPu6Calomatch+akPu6Caloparton+akPu6CalopatJets)
|
561 |
|
|
|
562 |
|
|
ak1PFpatSequence = cms.Sequence(ak1PFcorr+ak1clean+ak1PFmatch+ak1PFparton+ak1PFpatJets)
|
563 |
|
|
ak2PFpatSequence = cms.Sequence(ak2PFcorr+ak2clean+ak2PFmatch+ak2PFparton+ak2PFpatJets)
|
564 |
|
|
ak3PFpatSequence = cms.Sequence(ak3PFcorr+ak3clean+ak3PFmatch+ak3PFparton+ak3PFpatJets)
|
565 |
|
|
ak4PFpatSequence = cms.Sequence(ak4PFcorr+ak4clean+ak4PFmatch+ak4PFparton+ak4PFpatJets)
|
566 |
|
|
ak5PFpatSequence = cms.Sequence(ak5PFcorr+ak5clean+ak5PFmatch+ak5PFparton+ak5PFpatJets)
|
567 |
|
|
ak6PFpatSequence = cms.Sequence(ak6PFcorr+ak6clean+ak6PFmatch+ak6PFparton+ak6PFpatJets)
|
568 |
|
|
|
569 |
|
|
ak1CalopatSequence = cms.Sequence(ak1Calocorr+ak1clean+ak1Calomatch+ak1Caloparton+ak1CalopatJets)
|
570 |
|
|
ak2CalopatSequence = cms.Sequence(ak2Calocorr+ak2clean+ak2Calomatch+ak2Caloparton+ak2CalopatJets)
|
571 |
|
|
ak3CalopatSequence = cms.Sequence(ak3Calocorr+ak3clean+ak3Calomatch+ak3Caloparton+ak3CalopatJets)
|
572 |
|
|
ak4CalopatSequence = cms.Sequence(ak4Calocorr+ak4clean+ak4Calomatch+ak4Caloparton+ak4CalopatJets)
|
573 |
|
|
ak5CalopatSequence = cms.Sequence(ak5Calocorr+ak5clean+ak5Calomatch+ak5Caloparton+ak5CalopatJets)
|
574 |
|
|
ak6CalopatSequence = cms.Sequence(ak6Calocorr+ak6clean+ak6Calomatch+ak6Caloparton+ak6CalopatJets)
|
575 |
|
|
|
576 |
|
|
makeHeavyIonJets = cms.Sequence(icPu5patSequence +
|
577 |
|
|
|
578 |
|
|
akPu1PFpatSequence +
|
579 |
|
|
akPu2PFpatSequence +
|
580 |
|
|
akPu3PFpatSequence +
|
581 |
|
|
akPu4PFpatSequence +
|
582 |
|
|
akPu5PFpatSequence +
|
583 |
|
|
akPu6PFpatSequence +
|
584 |
|
|
|
585 |
|
|
akPu1CalopatSequence +
|
586 |
|
|
akPu2CalopatSequence +
|
587 |
|
|
akPu3CalopatSequence +
|
588 |
|
|
akPu4CalopatSequence +
|
589 |
|
|
akPu5CalopatSequence +
|
590 |
|
|
akPu6CalopatSequence +
|
591 |
|
|
|
592 |
|
|
ak1PFpatSequence +
|
593 |
|
|
ak2PFpatSequence +
|
594 |
|
|
ak3PFpatSequence +
|
595 |
|
|
ak4PFpatSequence +
|
596 |
|
|
ak5PFpatSequence +
|
597 |
|
|
ak6PFpatSequence +
|
598 |
|
|
|
599 |
|
|
ak1CalopatSequence +
|
600 |
|
|
ak2CalopatSequence +
|
601 |
|
|
ak3CalopatSequence +
|
602 |
|
|
ak4CalopatSequence +
|
603 |
|
|
ak5CalopatSequence +
|
604 |
|
|
ak6CalopatSequence
|
605 |
|
|
|
606 |
|
|
|
607 |
|
|
)
|
608 |
|
|
|
609 |
|
|
|
610 |
|
|
ic5JetAnalyzer = icPu5JetAnalyzer.clone(jetTag = cms.InputTag("iterativeCone5CaloJets"))
|
611 |
|
|
|
612 |
|
|
akPu1PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak1PFpatJets"))
|
613 |
|
|
akPu2PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak2PFpatJets"))
|
614 |
|
|
akPu3PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak3PFpatJets"))
|
615 |
|
|
akPu4PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak4PFpatJets"))
|
616 |
|
|
akPu5PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak5PFpatJets"))
|
617 |
|
|
akPu6PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak6PFpatJets"))
|
618 |
|
|
|
619 |
|
|
akPu3CaloJetAnalyzer = akPu3PFJetAnalyzer.clone()
|
620 |
|
|
|
621 |
|
|
akPu1CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak1CalopatJets"))
|
622 |
|
|
akPu2CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak2CalopatJets"))
|
623 |
|
|
akPu3CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak3CalopatJets"))
|
624 |
|
|
akPu4CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak4CalopatJets"))
|
625 |
|
|
akPu5CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak5CalopatJets"))
|
626 |
|
|
akPu6CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak6CalopatJets"))
|
627 |
|
|
|
628 |
|
|
ak1PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak1PFpatJets"))
|
629 |
|
|
ak2PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak2PFpatJets"))
|
630 |
|
|
ak3PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak3PFpatJets"))
|
631 |
|
|
ak4PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak4PFpatJets"))
|
632 |
|
|
ak5PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak5PFpatJets"))
|
633 |
|
|
ak6PFJetAnalyzer = akPu3PFJetAnalyzer.clone(jetTag = cms.InputTag("ak6PFpatJets"))
|
634 |
|
|
|
635 |
|
|
ak3CaloJetAnalyzer = akPu3PFJetAnalyzer.clone()
|
636 |
|
|
|
637 |
|
|
ak1CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak1CalopatJets"))
|
638 |
|
|
ak2CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak2CalopatJets"))
|
639 |
|
|
ak3CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak3CalopatJets"))
|
640 |
|
|
ak4CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak4CalopatJets"))
|
641 |
|
|
ak5CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak5CalopatJets"))
|
642 |
|
|
ak6CaloJetAnalyzer = akPu3CaloJetAnalyzer.clone(jetTag = cms.InputTag("ak6CalopatJets"))
|
643 |
|
|
|
644 |
|
|
jetAnalyzers = cms.Sequence(
|
645 |
yilmaz |
1.10 |
akPu1PFJetAnalyzer +
|
646 |
|
|
akPu2PFJetAnalyzer +
|
647 |
|
|
akPu3PFJetAnalyzer +
|
648 |
|
|
akPu4PFJetAnalyzer +
|
649 |
|
|
akPu5PFJetAnalyzer +
|
650 |
|
|
akPu6PFJetAnalyzer +
|
651 |
|
|
|
652 |
|
|
akPu1CaloJetAnalyzer +
|
653 |
|
|
akPu2CaloJetAnalyzer +
|
654 |
|
|
akPu3CaloJetAnalyzer +
|
655 |
|
|
akPu4CaloJetAnalyzer +
|
656 |
|
|
akPu5CaloJetAnalyzer +
|
657 |
|
|
akPu6CaloJetAnalyzer +
|
658 |
|
|
|
659 |
|
|
ak1PFJetAnalyzer +
|
660 |
|
|
ak2PFJetAnalyzer +
|
661 |
|
|
ak3PFJetAnalyzer +
|
662 |
|
|
ak4PFJetAnalyzer +
|
663 |
|
|
ak5PFJetAnalyzer +
|
664 |
|
|
ak6PFJetAnalyzer +
|
665 |
|
|
|
666 |
|
|
ak1CaloJetAnalyzer +
|
667 |
|
|
ak2CaloJetAnalyzer +
|
668 |
|
|
ak3CaloJetAnalyzer +
|
669 |
|
|
ak4CaloJetAnalyzer +
|
670 |
|
|
ak5CaloJetAnalyzer +
|
671 |
|
|
ak6CaloJetAnalyzer
|
672 |
|
|
|
673 |
|
|
)
|
674 |
yilmaz |
1.9 |
|
675 |
|
|
|
676 |
|
|
|
677 |
|
|
|
678 |
|
|
|
679 |
|
|
|
680 |
|
|
|
681 |
|
|
|
682 |
|
|
|