ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/PatAna_cff.py
Revision: 1.11
Committed: Wed May 9 13:53:56 2012 UTC (12 years, 11 months ago) by yilmaz
Content type: text/x-python
Branch: MAIN
Changes since 1.10: +1 -2 lines
Log Message:
up

File Contents

# User Rev Content
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