ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/HbbAnalysis/python/jetConfig_cff.py
Revision: 1.2
Committed: Fri Apr 17 14:26:35 2009 UTC (16 years ago) by amagnan
Content type: text/x-python
Branch: MAIN
CVS Tags: v00-03-00, v00-02-01, v00-02-00, v00-01-00
Changes since 1.1: +25 -24 lines
Log Message:
add Tau histos+debug btag for JPT/PF jets

File Contents

# User Rev Content
1 amagnan 1.1 ########################################
2     #### Created 07-Apr-09
3     #### by A.-M. Magnan
4     #### a.magnan@imperial.ac.uk
5     #### from TauAnalysis globalV00-00-03
6     ########################################
7    
8     import FWCore.ParameterSet.Config as cms
9     import copy
10    
11     ######################################################################################
12     #### PAT layer 0 cleaning
13     ######################################################################################
14    
15     from PhysicsTools.PatAlgos.cleaningLayer0.caloJetCleaner_cfi import *
16    
17     allLayer0JetsForTauAnalyses = copy.deepcopy(allLayer0Jets)
18     allLayer0JetsForTauAnalyses.jetSource = cms.InputTag("iterativeCone5CaloJets")
19     allLayer0JetsForTauAnalyses.removeOverlaps.electrons = cms.PSet() # disable default to remove overlap with electrons
20     #allLayer0JetsForTauAnalyses.removeOverlaps.electrons.collection = cms.InputTag("allLayer0ElectronsForTauAnalyses")
21     #allLayer0JetsForTauAnalyses.removeOverlaps.electrons.deltaR = cms.double(0.3) ##
22     #allLayer0JetsForTauAnalyses.removeOverlaps.electrons.cut = cms.string('pt > 10 & (trackIso()+ecalIso())/pt < 0.18')
23     #allLayer0JetsForTauAnalyses.removeOverlaps.electrons.flags = cms.vstring('Isolation/Tracker')
24    
25     allLayer0JetsForTauAnalyses.removeOverlaps.taus = cms.PSet() # disable default to remove overlap with electrons
26     #allLayer0JetsForTauAnalyses.removeOverlaps.taus.collection = cms.InputTag("allLayer0PFTausForTauAnalyses")
27     #allLayer0JetsForTauAnalyses.removeOverlaps.taus.deltaR = cms.double(0.3) ##
28     #allLayer0JetsForTauAnalyses.removeOverlaps.taus.cut = cms.string('')
29     #allLayer0JetsForTauAnalyses.removeOverlaps.taus.flags = cms.vstring('Isolation/Tracker')
30    
31    
32     ######################################################################################
33     #### Flavour ID
34     ######################################################################################
35     from PhysicsTools.PatAlgos.recoLayer0.jetFlavourId_cff import *
36    
37    
38     jetPartonAssociationForTauAnalyses = copy.deepcopy(jetPartonAssociation)
39     jetPartonAssociationForTauAnalyses.jets = cms.InputTag("allLayer0JetsForTauAnalyses")
40     jetPartonAssociationForTauAnalyses.partons = cms.InputTag("jetPartons")
41     jetPartonAssociationForTauAnalyses.coneSizeToAssociate = cms.double(0.3)
42    
43     jetFlavourAssociationForTauAnalyses = copy.deepcopy(jetFlavourAssociation)
44     jetFlavourAssociationForTauAnalyses.srcByReference = cms.InputTag("jetPartonAssociationForTauAnalyses")
45     jetFlavourAssociationForTauAnalyses.physicsDefinition = cms.bool(False)
46 amagnan 1.2 jetFlavourAssociationForTauAnalyses.definition = cms.int32(4)
47 amagnan 1.1
48     #PAT sequence for jet flavour identification
49     patJetFlavourIdForTauAnalyses = cms.Sequence(jetPartons
50     *jetPartonAssociationForTauAnalyses
51     *jetFlavourAssociationForTauAnalyses
52     )
53    
54    
55     ######################################################################################
56     #### b-tagging
57     ######################################################################################
58     from PhysicsTools.PatAlgos.recoLayer0.bTagging_cff import *
59    
60    
61     # Labels of Discriminators to use
62 amagnan 1.2 #patBTaggingDiscriminatorLabelsForTauAnalyses = cms.VInputTag(
63     # cms.InputTag("combinedSecondaryVertexBJetTags","RECO"),
64     # cms.InputTag("combinedSecondaryVertexMVABJetTags","RECO"),
65     # cms.InputTag("coneIsolationTauJetTags","RECO"),
66     # cms.InputTag("impactParameterMVABJetTags","RECO"),
67     # cms.InputTag("jetBProbabilityBJetTags","RECO"),
68     # cms.InputTag("jetProbabilityBJetTags","RECO"),
69     # cms.InputTag("simpleSecondaryVertexBJetTags","RECO"),
70     # cms.InputTag("softElectronBJetTags","RECO"),
71     # cms.InputTag("softMuonBJetTags","RECO"),
72     # cms.InputTag("softMuonNoIPBJetTags","RECO"),
73     # cms.InputTag("trackCountingHighEffBJetTags","RECO"),
74     # cms.InputTag("trackCountingHighPurBJetTags","RECO"),
75     #)
76 amagnan 1.1 # Labels of TagInfos to use
77 amagnan 1.2 #patBTaggingTagInfoLabelsForTauAnalyses = cms.VInputTag(
78     # cms.InputTag("secondaryVertexTagInfos","RECO"),
79     # cms.InputTag("softElectronTagInfos","RECO"),
80     # cms.InputTag("softMuonTagInfos","RECO"),
81     # cms.InputTag("impactParameterTagInfos","RECO"),
82     #)
83 amagnan 1.1
84    
85     # Need to convert all JetTags to ValueMap<double>
86     patAODBTagsForTauAnalyses = copy.deepcopy(patAODBTags)
87     patAODBTagsForTauAnalyses.collection = cms.InputTag("iterativeCone5CaloJets")
88 amagnan 1.2 patAODBTagsForTauAnalyses.associations = patBTaggingDiscriminatorLabels #ForTauAnalyses
89 amagnan 1.1 patAODBTagsForTauAnalyses.failSilently = cms.untracked.bool(True)
90    
91    
92     # Need to convert all JetTagInfoss to ValueMap<Ptr<BaseTagInfo>>
93     patAODTagInfosForTauAnalyses = copy.deepcopy(patAODTagInfos)
94     patAODTagInfosForTauAnalyses.collection = cms.InputTag("iterativeCone5CaloJets")
95 amagnan 1.2 patAODTagInfosForTauAnalyses.associations = patBTaggingTagInfoLabels #ForTauAnalyses
96 amagnan 1.1 patAODTagInfosForTauAnalyses.failSilently = cms.untracked.bool(True)
97    
98     layer0BTagsForTauAnalyses = copy.deepcopy(layer0BTags)
99     layer0BTagsForTauAnalyses.collection = cms.InputTag("allLayer0JetsForTauAnalyses")
100     layer0BTagsForTauAnalyses.backrefs = cms.InputTag("allLayer0JetsForTauAnalyses")
101     layer0BTagsForTauAnalyses.commonLabel = cms.InputTag("patAODBTagsForTauAnalyses")
102 amagnan 1.2 layer0BTagsForTauAnalyses.associations = patBTaggingDiscriminatorLabels #ForTauAnalyses
103 amagnan 1.1 layer0BTagsForTauAnalyses.failSilently = cms.untracked.bool(True)
104    
105     layer0TagInfosForTauAnalyses = copy.deepcopy(layer0TagInfos)
106     layer0TagInfosForTauAnalyses.collection = cms.InputTag("allLayer0JetsForTauAnalyses")
107     layer0TagInfosForTauAnalyses.backrefs = cms.InputTag("allLayer0JetsForTauAnalyses")
108     layer0TagInfosForTauAnalyses.commonLabel = cms.InputTag("patAODTagInfosForTauAnalyses")
109 amagnan 1.2 layer0TagInfosForTauAnalyses.associations = patBTaggingTagInfoLabels #ForTauAnalyses
110 amagnan 1.1 layer0TagInfosForTauAnalyses.failSilently = cms.untracked.bool(True)
111    
112     patAODBTaggingForTauAnalyses = cms.Sequence(patAODBTagsForTauAnalyses
113     *patAODTagInfosForTauAnalyses
114     )
115     patLayer0BTaggingForTauAnalyses = cms.Sequence(layer0BTagsForTauAnalyses
116     *layer0TagInfosForTauAnalyses
117     )
118    
119    
120     ######################################################################################
121     #### jet-tracks charge
122     ######################################################################################
123     from PhysicsTools.PatAlgos.recoLayer0.jetTracksCharge_cff import *
124    
125     patAODJetTracksAssociatorForTauAnalyses = copy.deepcopy(patAODJetTracksAssociator)
126     patAODJetTracksAssociatorForTauAnalyses.src = cms.InputTag("iterativeCone5CaloJets") ## the AOD jets given as input to the PAT jet cleaner
127     patAODJetTracksAssociatorForTauAnalyses.tracks = cms.InputTag("ic5JetTracksAssociatorAtVertex") ## any JetTracksAssociation
128     patAODJetTracksAssociatorForTauAnalyses.cut = cms.string('') # e.g. normalizedChi2 < 5
129    
130     ## Re-key from AOD jets to PAT Layer 0 jets
131     layer0JetTracksAssociatorForTauAnalyses = copy.deepcopy(layer0JetTracksAssociator)
132     layer0JetTracksAssociatorForTauAnalyses.collection = cms.InputTag("allLayer0JetsForTauAnalyses")
133     layer0JetTracksAssociatorForTauAnalyses.backrefs = cms.InputTag("allLayer0JetsForTauAnalyses")
134     layer0JetTracksAssociatorForTauAnalyses.association = cms.InputTag("patAODJetTracksAssociatorForTauAnalyses")
135    
136     ## Compute JET Charge
137     layer0JetChargeForTauAnalyses = copy.deepcopy(layer0JetCharge)
138     layer0JetChargeForTauAnalyses.src = cms.InputTag("allLayer0JetsForTauAnalyses") ## The Jets
139     layer0JetChargeForTauAnalyses.jetTracksAssociation = cms.InputTag("layer0JetTracksAssociatorForTauAnalyses")
140     # -- JetCharge parameters --
141     layer0JetChargeForTauAnalyses.var = cms.string('Pt')
142     layer0JetChargeForTauAnalyses.exp = cms.double(1.0)
143    
144     patLayer0JetTracksChargeForTauAnalyses = cms.Sequence(patAODJetTracksAssociatorForTauAnalyses
145     *layer0JetTracksAssociatorForTauAnalyses
146     *layer0JetChargeForTauAnalyses
147     )
148    
149    
150     ######################################################################################
151     #### jetMET corrections
152     ######################################################################################
153     from PhysicsTools.PatAlgos.recoLayer0.jetMETCorrections_cff import *
154     from JetMETCorrections.Configuration.L2L3Corrections_Summer08Redigi_cff import *
155    
156     #just to be sure it takes redigi values....
157     L2JetCorrectorIC5Calo.tagName = cms.string('Summer08Redigi_L2Relative_IC5Calo')
158     L3JetCorrectorIC5Calo.tagName = cms.string('Summer08Redigi_L3Absolute_IC5Calo')
159    
160     # change corrector to L2+L3
161     corMetType1Icone5ForTauAnalyses = copy.deepcopy(corMetType1Icone5)
162     corMetType1Icone5ForTauAnalyses.corrector = cms.string('L2L3JetCorrectorIC5Calo')
163    
164    
165     corMetType1Icone5MuonsForTauAnalyses = copy.deepcopy(corMetType1Icone5Muons)
166     corMetType1Icone5MuonsForTauAnalyses.uncorMETInputTag = cms.InputTag('corMetType1Icone5ForTauAnalyses')
167     corMetType1Icone5MuonsForTauAnalyses.muonsInputTag = cms.InputTag('goodMuonsforMETCorrection')
168    
169     # It would be better to get this config to JetMETCorrections/Type1MET/data/ at some point
170     corMetType1Icone5MuonsForTauAnalyses.TrackAssociatorParameters.useEcal = False ## RecoHits
171     corMetType1Icone5MuonsForTauAnalyses.TrackAssociatorParameters.useHcal = False ## RecoHits
172     corMetType1Icone5MuonsForTauAnalyses.TrackAssociatorParameters.useHO = False ## RecoHits
173     corMetType1Icone5MuonsForTauAnalyses.TrackAssociatorParameters.useCalo = True ## CaloTowers
174     corMetType1Icone5MuonsForTauAnalyses.TrackAssociatorParameters.useMuon = False ## RecoHits
175     corMetType1Icone5MuonsForTauAnalyses.TrackAssociatorParameters.truthMatch = False
176    
177     #define jetCorrFactors to associate
178     jetCorrFactorsForTauAnalyses = copy.deepcopy(jetCorrFactors)
179     jetCorrFactorsForTauAnalyses.useEMF = cms.bool(False) ## the use of emf in the JEC
180     ## choose sample type for flavor dependend corrections:
181     jetCorrFactorsForTauAnalyses.sampleType= cms.int32(0) ## 0: from dijet sample, 1: from ttbar sample
182     jetCorrFactorsForTauAnalyses.jetSource = cms.InputTag("iterativeCone5CaloJets") ## input collection of jets
183     ## tags for the jet correctors; when not available the string should be set to 'none'
184     jetCorrFactorsForTauAnalyses.L1Offset = cms.string('none')
185     jetCorrFactorsForTauAnalyses.L2Relative= cms.string('Summer08Redigi_L2Relative_IC5Calo')
186     jetCorrFactorsForTauAnalyses.L3Absolute= cms.string('Summer08Redigi_L3Absolute_IC5Calo')
187     jetCorrFactorsForTauAnalyses.L4EMF = cms.string('none')
188     jetCorrFactorsForTauAnalyses.L5Flavor = cms.string('none') #L5Flavor_IC5
189     jetCorrFactorsForTauAnalyses.L6UE = cms.string('none')
190     jetCorrFactorsForTauAnalyses.L7Parton = cms.string('none') #L7Parton_IC5
191    
192     # re-key jet energy corrections to layer 0 output
193     layer0JetCorrFactorsForTauAnalyses = copy.deepcopy(layer0JetCorrFactors)
194     layer0JetCorrFactorsForTauAnalyses.collection = cms.InputTag("allLayer0JetsForTauAnalyses")
195     layer0JetCorrFactorsForTauAnalyses.backrefs = cms.InputTag("allLayer0JetsForTauAnalyses")
196     layer0JetCorrFactorsForTauAnalyses.association = cms.InputTag("jetCorrFactorsForTauAnalyses")
197    
198    
199     # default PAT sequence for JetMET corrections before cleaners
200     patAODJetMETCorrectionsForTauAnalyses = cms.Sequence(jetCorrFactorsForTauAnalyses
201     +goodMuonsforMETCorrection
202     *corMetType1Icone5ForTauAnalyses
203     *corMetType1Icone5MuonsForTauAnalyses
204     )
205    
206     # default PAT sequence for JetMET corrections after cleaners
207     patLayer0JetMETCorrectionsForTauAnalyses = cms.Sequence(layer0JetCorrFactorsForTauAnalyses)
208    
209    
210    
211     ######################################################################################
212     #### trigger matching
213     ######################################################################################
214    
215    
216     ######################################################################################
217     #### MC matching : genjet and parton
218     ######################################################################################
219     from PhysicsTools.PatAlgos.mcMatchLayer0.jetMatch_cfi import *
220    
221     jetPartonMatchForTauAnalyses = copy.deepcopy(jetPartonMatch)
222     jetPartonMatchForTauAnalyses.src = cms.InputTag("allLayer0JetsForTauAnalyses") # RECO objects to match
223     jetPartonMatchForTauAnalyses.matched = cms.InputTag("genParticles") # mc-truth particle collection
224     jetPartonMatchForTauAnalyses.mcPdgId = cms.vint32(1, 2, 3, 4, 5, 21) # one or more PDG ID (quarks except top; gluons)
225     jetPartonMatchForTauAnalyses.mcStatus = cms.vint32(3) # PYTHIA status code (3 = hard scattering)
226     jetPartonMatchForTauAnalyses.checkCharge = cms.bool(False) # False = any value of the charge of MC and RECO is ok
227     jetPartonMatchForTauAnalyses.maxDeltaR = cms.double(0.4) # Minimum deltaR for the match
228     jetPartonMatchForTauAnalyses.maxDPtRel = cms.double(3.0) # Minimum deltaPt/Pt for the match
229     jetPartonMatchForTauAnalyses.resolveAmbiguities = cms.bool(True) # Forbid two RECO objects to match to the same GEN object
230     jetPartonMatchForTauAnalyses.resolveByMatchQuality = cms.bool(False) # False = just match input in order; True = pick lowest deltaR pair first
231    
232     jetGenJetMatchForTauAnalyses = copy.deepcopy(jetGenJetMatch)
233     jetGenJetMatchForTauAnalyses.src = cms.InputTag("allLayer0JetsForTauAnalyses") ## RECO jets (any View<Jet> is ok)
234     jetGenJetMatchForTauAnalyses.matched = cms.InputTag("iterativeCone5GenJetsNoNuBSM") ## GEN jets (must be GenJetCollection)
235     jetGenJetMatchForTauAnalyses.maxDeltaR = cms.double(0.4) # Minimum deltaR for the match
236     jetGenJetMatchForTauAnalyses.maxDPtRel = cms.double(3.0) # Minimum deltaPt/Pt for the match
237     jetGenJetMatchForTauAnalyses.resolveAmbiguities = cms.bool(True) # Forbid two RECO objects to match to the same GEN object
238     jetGenJetMatchForTauAnalyses.resolveByMatchQuality = cms.bool(False) # False = just match input in order; True = pick lowest deltaR pair first
239    
240    
241     ######################################################################################
242     #### PAT layer 1
243     ######################################################################################
244    
245     from PhysicsTools.PatAlgos.producersLayer1.jetProducer_cff import *
246    
247     allLayer1JetsForTauAnalyses = copy.deepcopy(allLayer1Jets)
248    
249     allLayer1JetsForTauAnalyses.jetSource = cms.InputTag("allLayer0JetsForTauAnalyses")
250    
251     allLayer1JetsForTauAnalyses.embedCaloTowers = cms.bool(True)
252    
253     allLayer1JetsForTauAnalyses.addJetCorrFactors = cms.bool(True)
254     allLayer1JetsForTauAnalyses.jetCorrFactorsSource = cms.InputTag("layer0JetCorrFactorsForTauAnalyses")
255    
256     allLayer1JetsForTauAnalyses.addBTagInfo = cms.bool(True)
257     allLayer1JetsForTauAnalyses.addDiscriminators = cms.bool(True)
258     allLayer1JetsForTauAnalyses.discriminatorModule = cms.InputTag("layer0BTagsForTauAnalyses")
259     allLayer1JetsForTauAnalyses.discriminatorNames = cms.vstring('*')
260     allLayer1JetsForTauAnalyses.addTagInfoRefs = cms.bool(True)
261     allLayer1JetsForTauAnalyses.tagInfoModule = cms.InputTag("layer0TagInfosForTauAnalyses")
262     allLayer1JetsForTauAnalyses.tagInfoNames = cms.vstring('secondaryVertexTagInfos','softElectronTagInfos','softMuonTagInfos','impactParameterTagInfos')
263    
264     allLayer1JetsForTauAnalyses.addAssociatedTracks = cms.bool(True)
265     allLayer1JetsForTauAnalyses.trackAssociationSource = cms.InputTag("layer0JetTracksAssociatorForTauAnalyses")
266    
267     allLayer1JetsForTauAnalyses.addJetCharge = cms.bool(True)
268     allLayer1JetsForTauAnalyses.jetChargeSource = cms.InputTag("layer0JetChargeForTauAnalyses")
269    
270     allLayer1JetsForTauAnalyses.addTrigMatch = cms.bool(False)
271     #allLayer1JetsForTauAnalyses.trigPrimMatch = cms.VInputTag(cms.InputTag("jetTrigMatchHLT1ElectronRelaxed"), cms.InputTag("jetTrigMatchHLT2jet"))
272    
273     allLayer1JetsForTauAnalyses.addGenPartonMatch = cms.bool(True)
274     allLayer1JetsForTauAnalyses.embedGenPartonMatch = cms.bool(True)
275     allLayer1JetsForTauAnalyses.genPartonMatch = cms.InputTag("jetPartonMatchForTauAnalyses")
276     allLayer1JetsForTauAnalyses.addGenJetMatch = cms.bool(True)
277     allLayer1JetsForTauAnalyses.genJetMatch = cms.InputTag("jetGenJetMatchForTauAnalyses")
278     allLayer1JetsForTauAnalyses.addPartonJetMatch = cms.bool(False)
279     allLayer1JetsForTauAnalyses.partonJetSource = cms.InputTag("NOT_IMPLEMENTED")
280    
281     allLayer1JetsForTauAnalyses.getJetMCFlavour = cms.bool(True)
282     allLayer1JetsForTauAnalyses.JetPartonMapSource = cms.InputTag("jetFlavourAssociationForTauAnalyses")
283    
284    
285    
286     ######################################################################################
287     #### Final PAT sequences
288     ######################################################################################
289    
290     patBeforeLayer0JetsForTauAnalyses = cms.Sequence( patAODBTaggingForTauAnalyses
291     #*patAODJetMETCorrectionsForTauAnalyses
292     #*patAODJetMETCorrections --> no need to do this here, as it is already done in met
293     )
294    
295     patLayer0JetsForTauAnalyses = cms.Sequence( patBeforeLayer0JetsForTauAnalyses
296     *allLayer0JetsForTauAnalyses
297     *patJetFlavourIdForTauAnalyses
298     *patLayer0BTaggingForTauAnalyses
299     *patLayer0JetMETCorrectionsForTauAnalyses
300     *patLayer0JetTracksChargeForTauAnalyses
301     *(jetPartonMatchForTauAnalyses+jetGenJetMatchForTauAnalyses)
302     )
303    
304     patLayer1JetsForTauAnalyses = cms.Sequence(allLayer1JetsForTauAnalyses)