ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/PatAna_MC_cff.py
Revision: 1.4
Committed: Tue May 29 12:20:19 2012 UTC (12 years, 11 months ago) by yilmaz
Content type: text/x-python
Branch: MAIN
Changes since 1.3: +1 -1 lines
Log Message:
fix

File Contents

# User Rev Content
1 yilmaz 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    
20     patJetCorrFactors.useNPV = False
21     # full reco
22     icPu5corr = patJetCorrFactors.clone(
23     src = cms.InputTag("iterativeConePu5CaloJets"),
24     levels = cms.vstring('L2Relative','L3Absolute'),
25 yilmaz 1.3 payload = cms.string('IC5Calo_HI')
26 yilmaz 1.1 )
27    
28     #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    
35     akPu5PFcorr = icPu5corr.clone(
36     src = cms.InputTag("akPu5PFJets"),
37 yilmaz 1.2 payload = cms.string('AK5PF_hiIterativeTracks')
38 yilmaz 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.2 payload = cms.string('AK3PF_hiIterativeTracks')
49 yilmaz 1.1 )
50    
51     #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    
58     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     payload = cms.string('AK5Calo')
73     )
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    
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    
97     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     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     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 yilmaz 1.4 payload = cms.string('IC5Calo_HI'))
222 yilmaz 1.1 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    
369     # === data sequences ===
370     # Note still need to use enableData function in cfg to remove mc dep of patjet
371    
372     # All corrections
373    
374    
375 yilmaz 1.2 akPu1PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu1PFJets"),payload = cms.string('AK1PF_hiIterativeTracks'))
376     akPu2PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu2PFJets"),payload = cms.string('AK2PF_hiIterativeTracks'))
377     akPu3PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu3PFJets"),payload = cms.string('AK3PF_hiIterativeTracks'))
378     akPu4PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu4PFJets"),payload = cms.string('AK4PF_hiIterativeTracks'))
379     akPu5PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu5PFJets"),payload = cms.string('AK5PF_hiIterativeTracks'))
380     akPu6PFcorr = akPu3PFcorr.clone(src = cms.InputTag("akPu6PFJets"),payload = cms.string('AK6PF_hiIterativeTracks'))
381    
382     ak1PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak1PFJets"),payload = cms.string('AK1PF_hiIterativeTracks'))
383     ak2PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak2PFJets"),payload = cms.string('AK2PF_hiIterativeTracks'))
384     ak3PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak3PFJets"),payload = cms.string('AK3PF_hiIterativeTracks'))
385     ak4PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak4PFJets"),payload = cms.string('AK4PF_hiIterativeTracks'))
386     ak5PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak5PFJets"),payload = cms.string('AK5PF_hiIterativeTracks'))
387     ak6PFcorr = akPu3PFcorr.clone(src = cms.InputTag("ak6PFJets"),payload = cms.string('AK6PF_hiIterativeTracks'))
388    
389 yilmaz 1.3 akPu1Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu1CaloJets"),payload = cms.string('AK1Calo_HI'))
390     akPu2Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu2CaloJets"),payload = cms.string('AK2Calo_HI'))
391     akPu3Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu3CaloJets"),payload = cms.string('AK3Calo_HI'))
392     akPu4Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu4CaloJets"),payload = cms.string('AK4Calo_HI'))
393     akPu5Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu5CaloJets"),payload = cms.string('AK5Calo_HI'))
394     akPu6Calocorr = akPu3PFcorr.clone(src = cms.InputTag("akPu6CaloJets"),payload = cms.string('AK6Calo_HI'))
395    
396     ak1Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak1CaloJets"),payload = cms.string('AK1Calo_HI'))
397     ak2Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak2CaloJets"),payload = cms.string('AK2Calo_HI'))
398     ak3Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak3CaloJets"),payload = cms.string('AK3Calo_HI'))
399     ak4Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak4CaloJets"),payload = cms.string('AK4Calo_HI'))
400     ak5Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak5CaloJets"),payload = cms.string('AK5Calo_HI'))
401     ak6Calocorr = akPu3PFcorr.clone(src = cms.InputTag("ak6CaloJets"),payload = cms.string('AK6Calo_HI'))
402 yilmaz 1.1
403     # Gen stuff
404    
405     ak1clean = akPu3PFclean.clone(src = cms.InputTag("ak1HiGenJets"))
406     ak2clean = akPu3PFclean.clone(src = cms.InputTag("ak2HiGenJets"))
407     ak3clean = akPu3PFclean.clone(src = cms.InputTag("ak3HiGenJets"))
408     ak4clean = akPu3PFclean.clone(src = cms.InputTag("ak4HiGenJets"))
409     ak5clean = akPu3PFclean.clone(src = cms.InputTag("ak5HiGenJets"))
410     ak6clean = akPu3PFclean.clone(src = cms.InputTag("ak6HiGenJets"))
411    
412    
413     akPu1PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu1PFJets"), matched = cms.InputTag("ak1clean"))
414     akPu2PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu2PFJets"), matched = cms.InputTag("ak2clean"))
415     akPu3PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu3PFJets"), matched = cms.InputTag("ak3clean"))
416     akPu4PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu4PFJets"), matched = cms.InputTag("ak4clean"))
417     akPu5PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu5PFJets"), matched = cms.InputTag("ak5clean"))
418     akPu6PFmatch = akPu3PFmatch.clone(src = cms.InputTag("akPu6PFJets"), matched = cms.InputTag("ak6clean"))
419     akPu1Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu1CaloJets"), matched = cms.InputTag("ak1clean"))
420     akPu2Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu2CaloJets"), matched = cms.InputTag("ak2clean"))
421     akPu3Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu3CaloJets"), matched = cms.InputTag("ak3clean"))
422     akPu4Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu4CaloJets"), matched = cms.InputTag("ak4clean"))
423     akPu5Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu5CaloJets"), matched = cms.InputTag("ak5clean"))
424     akPu6Calomatch = akPu3PFmatch.clone(src = cms.InputTag("akPu6CaloJets"), matched = cms.InputTag("ak6clean"))
425     ak1PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak1PFJets"), matched = cms.InputTag("ak1clean"))
426     ak2PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak2PFJets"), matched = cms.InputTag("ak2clean"))
427     ak3PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak3PFJets"), matched = cms.InputTag("ak3clean"))
428     ak4PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak4PFJets"), matched = cms.InputTag("ak4clean"))
429     ak5PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak5PFJets"), matched = cms.InputTag("ak5clean"))
430     ak6PFmatch = akPu3PFmatch.clone(src = cms.InputTag("ak6PFJets"), matched = cms.InputTag("ak6clean"))
431     ak1Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak1CaloJets"), matched = cms.InputTag("ak1clean"))
432     ak2Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak2CaloJets"), matched = cms.InputTag("ak2clean"))
433     ak3Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak3CaloJets"), matched = cms.InputTag("ak3clean"))
434     ak4Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak4CaloJets"), matched = cms.InputTag("ak4clean"))
435     ak5Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak5CaloJets"), matched = cms.InputTag("ak5clean"))
436     ak6Calomatch = akPu3PFmatch.clone(src = cms.InputTag("ak6CaloJets"), matched = cms.InputTag("ak6clean"))
437     akPu1PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu1PFJets"))
438     akPu2PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu2PFJets"))
439     akPu3PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu3PFJets"))
440     akPu4PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu4PFJets"))
441     akPu5PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu5PFJets"))
442     akPu6PFparton = akPu3PFparton.clone(src = cms.InputTag("akPu6PFJets"))
443     akPu1Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu1CaloJets"))
444     akPu2Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu2CaloJets"))
445     akPu3Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu3CaloJets"))
446     akPu4Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu4CaloJets"))
447     akPu5Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu5CaloJets"))
448     akPu6Caloparton = akPu3PFparton.clone(src = cms.InputTag("akPu6CaloJets"))
449     ak1PFparton = akPu3PFparton.clone(src = cms.InputTag("ak1PFJets"))
450     ak2PFparton = akPu3PFparton.clone(src = cms.InputTag("ak2PFJets"))
451     ak3PFparton = akPu3PFparton.clone(src = cms.InputTag("ak3PFJets"))
452     ak4PFparton = akPu3PFparton.clone(src = cms.InputTag("ak4PFJets"))
453     ak5PFparton = akPu3PFparton.clone(src = cms.InputTag("ak5PFJets"))
454     ak6PFparton = akPu3PFparton.clone(src = cms.InputTag("ak6PFJets"))
455     ak1Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak1CaloJets"))
456     ak2Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak2CaloJets"))
457     ak3Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak3CaloJets"))
458     ak4Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak4CaloJets"))
459     ak5Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak5CaloJets"))
460     ak6Caloparton = akPu3PFparton.clone(src = cms.InputTag("ak6CaloJets"))
461    
462    
463     # PAT Maker
464    
465     akPu1PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu1PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu1PFcorr")), genJetMatch = cms.InputTag("akPu1PFmatch"), genPartonMatch = cms.InputTag("akPu1PFparton"))
466     akPu2PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu2PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu2PFcorr")), genJetMatch = cms.InputTag("akPu2PFmatch"), genPartonMatch = cms.InputTag("akPu2PFparton"))
467     akPu3PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu3PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr")), genJetMatch = cms.InputTag("akPu3PFmatch"), genPartonMatch = cms.InputTag("akPu3PFparton"))
468     akPu4PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu4PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu4PFcorr")), genJetMatch = cms.InputTag("akPu4PFmatch"), genPartonMatch = cms.InputTag("akPu4PFparton"))
469     akPu5PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu5PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr")), genJetMatch = cms.InputTag("akPu5PFmatch"), genPartonMatch = cms.InputTag("akPu5PFparton"))
470     akPu6PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu6PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu6PFcorr")), genJetMatch = cms.InputTag("akPu6PFmatch"), genPartonMatch = cms.InputTag("akPu6PFparton"))
471     akPu1CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu1CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu1Calocorr")), genJetMatch = cms.InputTag("akPu1Calomatch"), genPartonMatch = cms.InputTag("akPu1Caloparton"))
472     akPu2CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu2CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu2Calocorr")), genJetMatch = cms.InputTag("akPu2Calomatch"), genPartonMatch = cms.InputTag("akPu2Caloparton"))
473     akPu3CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu3CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3Calocorr")), genJetMatch = cms.InputTag("akPu3Calomatch"), genPartonMatch = cms.InputTag("akPu3Caloparton"))
474     akPu4CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu4CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu4Calocorr")), genJetMatch = cms.InputTag("akPu4Calomatch"), genPartonMatch = cms.InputTag("akPu4Caloparton"))
475     akPu5CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu5CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5Calocorr")), genJetMatch = cms.InputTag("akPu5Calomatch"), genPartonMatch = cms.InputTag("akPu5Caloparton"))
476     akPu6CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("akPu6CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu6Calocorr")), genJetMatch = cms.InputTag("akPu6Calomatch"), genPartonMatch = cms.InputTag("akPu6Caloparton"))
477     ak1PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak1PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak1PFcorr")), genJetMatch = cms.InputTag("ak1PFmatch"), genPartonMatch = cms.InputTag("ak1PFparton"))
478     ak2PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak2PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak2PFcorr")), genJetMatch = cms.InputTag("ak2PFmatch"), genPartonMatch = cms.InputTag("ak2PFparton"))
479     ak3PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak3PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak3PFcorr")), genJetMatch = cms.InputTag("ak3PFmatch"), genPartonMatch = cms.InputTag("ak3PFparton"))
480     ak4PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak4PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak4PFcorr")), genJetMatch = cms.InputTag("ak4PFmatch"), genPartonMatch = cms.InputTag("ak4PFparton"))
481     ak5PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak5PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5PFcorr")), genJetMatch = cms.InputTag("ak5PFmatch"), genPartonMatch = cms.InputTag("ak5PFparton"))
482     ak6PFpatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak6PFJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak6PFcorr")), genJetMatch = cms.InputTag("ak6PFmatch"), genPartonMatch = cms.InputTag("ak6PFparton"))
483     ak1CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak1CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak1Calocorr")), genJetMatch = cms.InputTag("ak1Calomatch"), genPartonMatch = cms.InputTag("ak1Caloparton"))
484     ak2CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak2CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak2Calocorr")), genJetMatch = cms.InputTag("ak2Calomatch"), genPartonMatch = cms.InputTag("ak2Caloparton"))
485     ak3CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak3CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak3Calocorr")), genJetMatch = cms.InputTag("ak3Calomatch"), genPartonMatch = cms.InputTag("ak3Caloparton"))
486     ak4CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak4CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak4Calocorr")), genJetMatch = cms.InputTag("ak4Calomatch"), genPartonMatch = cms.InputTag("ak4Caloparton"))
487     ak5CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak5CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5Calocorr")), genJetMatch = cms.InputTag("ak5Calomatch"), genPartonMatch = cms.InputTag("ak5Caloparton"))
488     ak6CalopatJets = akPu3PFpatJets.clone(jetSource = cms.InputTag("ak6CaloJets"), jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak6Calocorr")), genJetMatch = cms.InputTag("ak6Calomatch"), genPartonMatch = cms.InputTag("ak6Caloparton"))
489    
490    
491    
492     icPu5patSequence = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5patJets)
493    
494     akPu1PFpatSequence = cms.Sequence(akPu1PFcorr+ak1clean+akPu1PFmatch+akPu1PFparton+akPu1PFpatJets)
495     akPu2PFpatSequence = cms.Sequence(akPu2PFcorr+ak2clean+akPu2PFmatch+akPu2PFparton+akPu2PFpatJets)
496     akPu3PFpatSequence = cms.Sequence(akPu3PFcorr+ak3clean+akPu3PFmatch+akPu3PFparton+akPu3PFpatJets)
497     akPu4PFpatSequence = cms.Sequence(akPu4PFcorr+ak4clean+akPu4PFmatch+akPu4PFparton+akPu4PFpatJets)
498     akPu5PFpatSequence = cms.Sequence(akPu5PFcorr+ak5clean+akPu5PFmatch+akPu5PFparton+akPu5PFpatJets)
499     akPu6PFpatSequence = cms.Sequence(akPu6PFcorr+ak6clean+akPu6PFmatch+akPu6PFparton+akPu6PFpatJets)
500    
501     akPu1CalopatSequence = cms.Sequence(akPu1Calocorr+ak1clean+akPu1Calomatch+akPu1Caloparton+akPu1CalopatJets)
502     akPu2CalopatSequence = cms.Sequence(akPu2Calocorr+ak2clean+akPu2Calomatch+akPu2Caloparton+akPu2CalopatJets)
503     akPu3CalopatSequence = cms.Sequence(akPu3Calocorr+ak3clean+akPu3Calomatch+akPu3Caloparton+akPu3CalopatJets)
504     akPu4CalopatSequence = cms.Sequence(akPu4Calocorr+ak4clean+akPu4Calomatch+akPu4Caloparton+akPu4CalopatJets)
505     akPu5CalopatSequence = cms.Sequence(akPu5Calocorr+ak5clean+akPu5Calomatch+akPu5Caloparton+akPu5CalopatJets)
506     akPu6CalopatSequence = cms.Sequence(akPu6Calocorr+ak6clean+akPu6Calomatch+akPu6Caloparton+akPu6CalopatJets)
507    
508     ak1PFpatSequence = cms.Sequence(ak1PFcorr+ak1clean+ak1PFmatch+ak1PFparton+ak1PFpatJets)
509     ak2PFpatSequence = cms.Sequence(ak2PFcorr+ak2clean+ak2PFmatch+ak2PFparton+ak2PFpatJets)
510     ak3PFpatSequence = cms.Sequence(ak3PFcorr+ak3clean+ak3PFmatch+ak3PFparton+ak3PFpatJets)
511     ak4PFpatSequence = cms.Sequence(ak4PFcorr+ak4clean+ak4PFmatch+ak4PFparton+ak4PFpatJets)
512     ak5PFpatSequence = cms.Sequence(ak5PFcorr+ak5clean+ak5PFmatch+ak5PFparton+ak5PFpatJets)
513     ak6PFpatSequence = cms.Sequence(ak6PFcorr+ak6clean+ak6PFmatch+ak6PFparton+ak6PFpatJets)
514    
515     ak1CalopatSequence = cms.Sequence(ak1Calocorr+ak1clean+ak1Calomatch+ak1Caloparton+ak1CalopatJets)
516     ak2CalopatSequence = cms.Sequence(ak2Calocorr+ak2clean+ak2Calomatch+ak2Caloparton+ak2CalopatJets)
517     ak3CalopatSequence = cms.Sequence(ak3Calocorr+ak3clean+ak3Calomatch+ak3Caloparton+ak3CalopatJets)
518     ak4CalopatSequence = cms.Sequence(ak4Calocorr+ak4clean+ak4Calomatch+ak4Caloparton+ak4CalopatJets)
519     ak5CalopatSequence = cms.Sequence(ak5Calocorr+ak5clean+ak5Calomatch+ak5Caloparton+ak5CalopatJets)
520     ak6CalopatSequence = cms.Sequence(ak6Calocorr+ak6clean+ak6Calomatch+ak6Caloparton+ak6CalopatJets)
521    
522     makeHeavyIonJets = cms.Sequence(icPu5patSequence +
523    
524     akPu1PFpatSequence +
525     akPu2PFpatSequence +
526     akPu3PFpatSequence +
527     akPu4PFpatSequence +
528     akPu5PFpatSequence +
529     akPu6PFpatSequence +
530    
531     akPu1CalopatSequence +
532     akPu2CalopatSequence +
533     akPu3CalopatSequence +
534     akPu4CalopatSequence +
535     akPu5CalopatSequence +
536     akPu6CalopatSequence +
537    
538     ak1PFpatSequence +
539     ak2PFpatSequence +
540     ak3PFpatSequence +
541     ak4PFpatSequence +
542     ak5PFpatSequence +
543     ak6PFpatSequence +
544    
545     ak1CalopatSequence +
546     ak2CalopatSequence +
547     ak3CalopatSequence +
548     ak4CalopatSequence +
549     ak5CalopatSequence +
550     ak6CalopatSequence
551    
552     )
553    
554     akPu3PFpatSequence_withBtagging = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFPatJetFlavourId * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFpatJets)
555    
556    
557    
558    
559    
560    
561