ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/ExtraJetReco_cff.py
Revision: 1.16
Committed: Wed Feb 13 18:19:26 2013 UTC (12 years, 2 months ago) by yilmaz
Content type: text/x-python
Branch: MAIN
Changes since 1.15: +1 -0 lines
Log Message:
TracksJets and rho calculation

File Contents

# User Rev Content
1 frankma 1.1 import FWCore.ParameterSet.Config as cms
2    
3     # reco jet with russian pileup subtraction
4     from RecoHI.HiJetAlgos.HiRecoJets_cff import *
5     from RecoHI.HiJetAlgos.HiRecoPFJets_cff import *
6 yjlee 1.10 from RecoJets.JetProducers.JetIDParams_cfi import *
7    
8 yilmaz 1.2 akPu3CaloJets = cms.EDProducer(
9     "FastjetJetProducer",
10     HiCaloJetParameters,
11     AnomalousCellParameters,
12     MultipleAlgoIteratorBlock,
13     jetAlgorithm = cms.string("AntiKt"),
14 yilmaz 1.12 rParam = cms.double(0.3),
15 yilmaz 1.2 )
16 yilmaz 1.12
17 yilmaz 1.2 akPu5CaloJets.radiusPU = 0.5
18    
19 frankma 1.1 iterativeConePu5CaloJets.doPVCorrection = cms.bool(True)
20     iterativeConePu5CaloJets.srcPVs = 'hiSelectedVertex'
21    
22 yjlee 1.15 iterativeConePu5CaloJets.jetPtMin = 15.0
23     ak5PFJets.jetPtMin = 15.0
24 yilmaz 1.7
25 frankma 1.1 akPu5PFJets = ak5PFJets.clone()
26     akPu5PFJets.src = 'PFTowers'
27     akPu5PFJets.jetType = 'BasicJet'
28     akPu5PFJets.doPUOffsetCorr = True
29     akPu5PFJets.sumRecHits = False
30    
31     akPu3PFJets = akPu5PFJets.clone()
32     akPu3PFJets.rParam = cms.double(0.3)
33    
34     # pileup subtraction jet exclusion pt min
35     iterativeConePu5CaloJets.puPtMin = cms.double(10.0)
36     akPu5PFJets.puPtMin = cms.double(25.0)
37     akPu3PFJets.puPtMin = cms.double(15.0)
38 yilmaz 1.2 akPu5CaloJets.puPtMin = cms.double(10.0)
39     akPu3CaloJets.puPtMin = cms.double(10.0)
40    
41 yilmaz 1.3 akPu5PFJets.doRhoFastjet = False
42     akPu5PFJets.doAreaFastjet = False
43    
44     akPu3PFJets.doRhoFastjet = False
45     akPu3PFJets.doAreaFastjet = False
46    
47 yilmaz 1.2 akPu5CaloJets.doRhoFastjet = False
48     akPu5CaloJets.doAreaFastjet = False
49     akPu5CaloJets.doPUOffsetCorr = True
50    
51     akPu3CaloJets.doRhoFastjet = False
52     akPu3CaloJets.doAreaFastjet = False
53     akPu3CaloJets.doPUOffsetCorr = True
54    
55    
56 yilmaz 1.4 ### Extra extended algos & sequence
57 yjlee 1.15 akPu6PFJets = akPu3PFJets.clone(rParam = 0.6, jetPtMin = 15)
58     akPu6CaloJets = akPu3CaloJets.clone(rParam = 0.6, jetPtMin = 15)
59     akPu5PFJets = akPu3PFJets.clone(rParam = 0.5, jetPtMin = 15)
60     akPu5CaloJets = akPu3CaloJets.clone(rParam = 0.5, jetPtMin = 15)
61     akPu4PFJets = akPu3PFJets.clone(rParam = 0.4, jetPtMin = 15)
62     akPu4CaloJets = akPu3CaloJets.clone(rParam = 0.4, jetPtMin = 15)
63     akPu2PFJets = akPu3PFJets.clone(rParam = 0.2, jetPtMin = 15)
64     akPu2CaloJets = akPu3CaloJets.clone(rParam = 0.2, jetPtMin = 15)
65     akPu1PFJets = akPu3PFJets.clone(rParam = 0.1, jetPtMin = 15)
66     akPu1CaloJets = akPu3CaloJets.clone(rParam = 0.1, jetPtMin = 15)
67    
68     iterativeCone5CaloJets = iterativeConePu5CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
69     ak6PFJets = akPu6PFJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
70     ak6CaloJets = akPu6CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
71     ak5PFJets = akPu5PFJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
72     ak5CaloJets = akPu5CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
73     ak4PFJets = akPu4PFJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
74     ak4CaloJets = akPu4CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
75     ak3PFJets = akPu3PFJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
76     ak3CaloJets = akPu3CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
77     ak2PFJets = akPu2PFJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
78     ak2CaloJets = akPu2CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
79     ak1PFJets = akPu1PFJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
80     ak1CaloJets = akPu1CaloJets.clone(doPUOffsetCorr = False, jetPtMin = 15)
81 yilmaz 1.12
82    
83     akPu1PFJets.puPtMin = 5
84     akPu2PFJets.puPtMin = 10
85     akPu3PFJets.puPtMin = 15
86     akPu4PFJets.puPtMin = 20
87     akPu5PFJets.puPtMin = 25
88     akPu6PFJets.puPtMin = 30
89    
90     akPu1CaloJets.puPtMin = 2
91     akPu2CaloJets.puPtMin = 4
92     akPu3CaloJets.puPtMin = 6
93     akPu4CaloJets.puPtMin = 8
94     akPu5CaloJets.puPtMin = 10
95     akPu6CaloJets.puPtMin = 12
96 yilmaz 1.5
97 yilmaz 1.16
98 yilmaz 1.9 #akPu1PFJets.puPtMin = cms.double(5.0)
99 yilmaz 1.12 #akPu2PFJets.puPtMin = cms.double(5.0)
100     #akPu3PFJets.puPtMin = cms.double(5.0)
101     #akPu4PFJets.puPtMin = cms.double(5.0)
102     #akPu5PFJets.puPtMin = cms.double(5.0)
103     #akPu6PFJets.puPtMin = cms.double(5.0)
104    
105     #akPu1CaloJets.puPtMin = cms.double(5.0)
106     #akPu2CaloJets.puPtMin = cms.double(5.0)
107     #akPu3CaloJets.puPtMin = cms.double(5.0)
108     #akPu4CaloJets.puPtMin = cms.double(5.0)
109     #akPu5CaloJets.puPtMin = cms.double(5.0)
110     #akPu6CaloJets.puPtMin = cms.double(5.0)
111 yilmaz 1.5
112 yjlee 1.10 # jet ID producer
113     ak1CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak1CaloJets'))
114     ak2CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak2CaloJets'))
115     ak3CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak3CaloJets'))
116     ak4CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak4CaloJets'))
117     ak5CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak5CaloJets'))
118     ak6CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak6CaloJets'))
119     akPu1CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu1CaloJets'))
120     akPu2CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu2CaloJets'))
121     akPu3CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu3CaloJets'))
122     akPu4CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu4CaloJets'))
123     akPu5CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu5CaloJets'))
124     akPu6CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu6CaloJets'))
125    
126     ak1PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak1PFJets'))
127     ak2PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak2PFJets'))
128     ak3PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak3PFJets'))
129     ak4PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak4PFJets'))
130     ak5PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak5PFJets'))
131     ak6PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak6PFJets'))
132     akPu1PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu1PFJets'))
133     akPu2PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu2PFJets'))
134     akPu3PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu3PFJets'))
135     akPu4PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu4PFJets'))
136     akPu5PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu5PFJets'))
137     akPu6PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu6PFJets'))
138    
139    
140 yilmaz 1.4 recoAk1to6 = cms.Sequence( akPu1PFJets * akPu2PFJets *akPu3PFJets * akPu4PFJets * akPu5PFJets * akPu6PFJets *
141     ak1PFJets * ak2PFJets *ak3PFJets * ak4PFJets * ak5PFJets * ak6PFJets *
142     akPu1CaloJets * akPu2CaloJets *akPu3CaloJets * akPu4CaloJets * akPu5CaloJets * akPu6CaloJets *
143     ak1CaloJets * ak2CaloJets *ak3CaloJets * ak4CaloJets * ak5CaloJets * ak6CaloJets
144     )
145    
146 yilmaz 1.13 recoAk1to6ID = cms.Sequence(
147     #akPu1PFJetID * akPu2PFJetID *akPu3PFJetID * akPu4PFJetID * akPu5PFJetID * akPu6PFJetID *
148     # ak1PFJetID * ak2PFJetID *ak3PFJetID * ak4PFJetID * ak5PFJetID * ak6PFJetID *
149 yjlee 1.10 akPu1CaloJetID * akPu2CaloJetID *akPu3CaloJetID * akPu4CaloJetID * akPu5CaloJetID * akPu6CaloJetID *
150     ak1CaloJetID * ak2CaloJetID *ak3CaloJetID * ak4CaloJetID * ak5CaloJetID * ak6CaloJetID
151     )
152    
153 yjlee 1.14 recoAk2to5 = cms.Sequence( akPu2PFJets *akPu3PFJets * akPu4PFJets * akPu5PFJets *
154     ak2PFJets *ak3PFJets * ak4PFJets * ak5PFJets *
155     akPu2CaloJets *akPu3CaloJets * akPu4CaloJets * akPu5CaloJets *
156     ak2CaloJets *ak3CaloJets * ak4CaloJets * ak5CaloJets
157     )
158    
159     recoAk2to5ID = cms.Sequence(
160     #akPu1PFJetID * akPu2PFJetID *akPu3PFJetID * akPu4PFJetID * akPu5PFJetID * akPu6PFJetID *
161     # ak1PFJetID * ak2PFJetID *ak3PFJetID * ak4PFJetID * ak5PFJetID * ak6PFJetID *
162     akPu2CaloJetID *akPu3CaloJetID * akPu4CaloJetID * akPu5CaloJetID *
163     ak2CaloJetID *ak3CaloJetID * ak4CaloJetID * ak5CaloJetID
164     )
165    
166 yilmaz 1.13 recoJetsWithID = cms.Sequence(recoAk1to6*recoAk1to6ID)
167 yjlee 1.14 recoJetsWithID2to5 = cms.Sequence(recoAk2to5*recoAk2to5ID)
168 yilmaz 1.13
169 yilmaz 1.4
170    
171    
172    
173    
174    
175