ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/ExtraJetReco_cff.py
Revision: 1.18
Committed: Fri Feb 15 16:22:42 2013 UTC (12 years, 2 months ago) by yilmaz
Content type: text/x-python
Branch: MAIN
CVS Tags: HiForest_V02_85, HiForest_V02_84, HiForest_V02_83, HiForest_V02_82, HiForest_V02_81, HiForest_V02_80, HiForest_V02_79, HiForest_V02_78, HiForest_V02_77, HiForest_V02_76, HEAD
Changes since 1.17: +6 -6 lines
Log Message:
bug fix eta window

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 yilmaz 1.17 ak1PFJets0 = ak1PFJets.clone(src = 'particleFlow',
84     jetType = 'PFJet'
85     )
86    
87     ak1PFJets0.doAreaFastjet = cms.bool(True)
88     ak1PFJets0.doRhoFastjet = cms.bool(True)
89     ak1PFJets0.doFastJetNonUniform = cms.bool(True)
90     ak1PFJets0.puCenters = cms.vdouble(-5,-4.5,-4,-3.5,-3,-2.5,-2,-1.5,-1,-0.5,0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5)
91     ak1PFJets0.puWidth = cms.double(0.5)
92     ak1PFJets0.nExclude = cms.uint32(0)
93     ak1PFJets0.Rho_EtaMax = cms.double(5.5)
94     ak1PFJets0.Ghost_EtaMax = cms.double(5.5)
95    
96     kt1PFJets0 = ak1PFJets0.clone(jetAlgorithm = cms.string("Kt"))
97    
98     ak1PFJets05 = ak1PFJets0.clone(rParam = cms.double(0.1),puWidth = 0.5)
99     ak2PFJets05 = ak1PFJets0.clone(rParam = cms.double(0.2),puWidth = 0.5)
100     ak3PFJets05 = ak1PFJets0.clone(rParam = cms.double(0.3),puWidth = 0.5)
101     ak4PFJets05 = ak1PFJets0.clone(rParam = cms.double(0.4),puWidth = 0.5)
102     ak5PFJets05 = ak1PFJets0.clone(rParam = cms.double(0.5),puWidth = 0.5)
103     ak6PFJets05 = ak1PFJets0.clone(rParam = cms.double(0.6),puWidth = 0.5)
104    
105     ak1PFJets1 = ak1PFJets0.clone(rParam = cms.double(0.1),puWidth = 2)
106     ak2PFJets1 = ak1PFJets0.clone(rParam = cms.double(0.2),puWidth = 2)
107     ak3PFJets1 = ak1PFJets0.clone(rParam = cms.double(0.3),puWidth = 2)
108     ak4PFJets1 = ak1PFJets0.clone(rParam = cms.double(0.4),puWidth = 2)
109     ak5PFJets1 = ak1PFJets0.clone(rParam = cms.double(0.5),puWidth = 2)
110     ak6PFJets1 = ak1PFJets0.clone(rParam = cms.double(0.6),puWidth = 2)
111    
112     ak1PFJets2 = ak1PFJets0.clone(rParam = cms.double(0.1),puWidth = 2)
113     ak2PFJets2 = ak1PFJets0.clone(rParam = cms.double(0.2),puWidth = 2)
114     ak3PFJets2 = ak1PFJets0.clone(rParam = cms.double(0.3),puWidth = 2)
115     ak4PFJets2 = ak1PFJets0.clone(rParam = cms.double(0.4),puWidth = 2)
116     ak5PFJets2 = ak1PFJets0.clone(rParam = cms.double(0.5),puWidth = 2)
117     ak6PFJets2 = ak1PFJets0.clone(rParam = cms.double(0.6),puWidth = 2)
118    
119    
120     kt1PFJets05 = kt1PFJets0.clone(rParam = cms.double(0.1),puWidth = 0.5)
121     kt2PFJets05 = kt1PFJets0.clone(rParam = cms.double(0.2),puWidth = 0.5)
122     kt3PFJets05 = kt1PFJets0.clone(rParam = cms.double(0.3),puWidth = 0.5)
123     kt4PFJets05 = kt1PFJets0.clone(rParam = cms.double(0.4),puWidth = 0.5)
124     kt5PFJets05 = kt1PFJets0.clone(rParam = cms.double(0.5),puWidth = 0.5)
125     kt6PFJets05 = kt1PFJets0.clone(rParam = cms.double(0.6),puWidth = 0.5)
126    
127 yilmaz 1.18 kt1PFJets1 = kt1PFJets0.clone(rParam = cms.double(0.1),puWidth = 1)
128     kt2PFJets1 = kt1PFJets0.clone(rParam = cms.double(0.2),puWidth = 1)
129     kt3PFJets1 = kt1PFJets0.clone(rParam = cms.double(0.3),puWidth = 1)
130     kt4PFJets1 = kt1PFJets0.clone(rParam = cms.double(0.4),puWidth = 1)
131     kt5PFJets1 = kt1PFJets0.clone(rParam = cms.double(0.5),puWidth = 1)
132     kt6PFJets1 = kt1PFJets0.clone(rParam = cms.double(0.6),puWidth = 1)
133 yilmaz 1.17
134     kt1PFJets2 = kt1PFJets0.clone(rParam = cms.double(0.1),puWidth = 2)
135     kt2PFJets2 = kt1PFJets0.clone(rParam = cms.double(0.2),puWidth = 2)
136     kt3PFJets2 = kt1PFJets0.clone(rParam = cms.double(0.3),puWidth = 2)
137     kt4PFJets2 = kt1PFJets0.clone(rParam = cms.double(0.4),puWidth = 2)
138     kt5PFJets2 = kt1PFJets0.clone(rParam = cms.double(0.5),puWidth = 2)
139     kt6PFJets2 = kt1PFJets0.clone(rParam = cms.double(0.6),puWidth = 2)
140    
141    
142    
143 yilmaz 1.12 akPu1PFJets.puPtMin = 5
144     akPu2PFJets.puPtMin = 10
145     akPu3PFJets.puPtMin = 15
146     akPu4PFJets.puPtMin = 20
147     akPu5PFJets.puPtMin = 25
148     akPu6PFJets.puPtMin = 30
149    
150     akPu1CaloJets.puPtMin = 2
151     akPu2CaloJets.puPtMin = 4
152     akPu3CaloJets.puPtMin = 6
153     akPu4CaloJets.puPtMin = 8
154     akPu5CaloJets.puPtMin = 10
155     akPu6CaloJets.puPtMin = 12
156 yilmaz 1.5
157 yilmaz 1.16
158 yilmaz 1.9 #akPu1PFJets.puPtMin = cms.double(5.0)
159 yilmaz 1.12 #akPu2PFJets.puPtMin = cms.double(5.0)
160     #akPu3PFJets.puPtMin = cms.double(5.0)
161     #akPu4PFJets.puPtMin = cms.double(5.0)
162     #akPu5PFJets.puPtMin = cms.double(5.0)
163     #akPu6PFJets.puPtMin = cms.double(5.0)
164    
165     #akPu1CaloJets.puPtMin = cms.double(5.0)
166     #akPu2CaloJets.puPtMin = cms.double(5.0)
167     #akPu3CaloJets.puPtMin = cms.double(5.0)
168     #akPu4CaloJets.puPtMin = cms.double(5.0)
169     #akPu5CaloJets.puPtMin = cms.double(5.0)
170     #akPu6CaloJets.puPtMin = cms.double(5.0)
171 yilmaz 1.5
172 yjlee 1.10 # jet ID producer
173     ak1CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak1CaloJets'))
174     ak2CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak2CaloJets'))
175     ak3CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak3CaloJets'))
176     ak4CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak4CaloJets'))
177     ak5CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak5CaloJets'))
178     ak6CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak6CaloJets'))
179     akPu1CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu1CaloJets'))
180     akPu2CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu2CaloJets'))
181     akPu3CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu3CaloJets'))
182     akPu4CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu4CaloJets'))
183     akPu5CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu5CaloJets'))
184     akPu6CaloJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu6CaloJets'))
185    
186     ak1PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak1PFJets'))
187     ak2PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak2PFJets'))
188     ak3PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak3PFJets'))
189     ak4PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak4PFJets'))
190     ak5PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak5PFJets'))
191     ak6PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('ak6PFJets'))
192     akPu1PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu1PFJets'))
193     akPu2PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu2PFJets'))
194     akPu3PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu3PFJets'))
195     akPu4PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu4PFJets'))
196     akPu5PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu5PFJets'))
197     akPu6PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akPu6PFJets'))
198    
199    
200 yilmaz 1.4 recoAk1to6 = cms.Sequence( akPu1PFJets * akPu2PFJets *akPu3PFJets * akPu4PFJets * akPu5PFJets * akPu6PFJets *
201     ak1PFJets * ak2PFJets *ak3PFJets * ak4PFJets * ak5PFJets * ak6PFJets *
202     akPu1CaloJets * akPu2CaloJets *akPu3CaloJets * akPu4CaloJets * akPu5CaloJets * akPu6CaloJets *
203     ak1CaloJets * ak2CaloJets *ak3CaloJets * ak4CaloJets * ak5CaloJets * ak6CaloJets
204     )
205    
206 yilmaz 1.13 recoAk1to6ID = cms.Sequence(
207     #akPu1PFJetID * akPu2PFJetID *akPu3PFJetID * akPu4PFJetID * akPu5PFJetID * akPu6PFJetID *
208     # ak1PFJetID * ak2PFJetID *ak3PFJetID * ak4PFJetID * ak5PFJetID * ak6PFJetID *
209 yjlee 1.10 akPu1CaloJetID * akPu2CaloJetID *akPu3CaloJetID * akPu4CaloJetID * akPu5CaloJetID * akPu6CaloJetID *
210     ak1CaloJetID * ak2CaloJetID *ak3CaloJetID * ak4CaloJetID * ak5CaloJetID * ak6CaloJetID
211     )
212    
213 yjlee 1.14 recoAk2to5 = cms.Sequence( akPu2PFJets *akPu3PFJets * akPu4PFJets * akPu5PFJets *
214     ak2PFJets *ak3PFJets * ak4PFJets * ak5PFJets *
215     akPu2CaloJets *akPu3CaloJets * akPu4CaloJets * akPu5CaloJets *
216     ak2CaloJets *ak3CaloJets * ak4CaloJets * ak5CaloJets
217     )
218    
219     recoAk2to5ID = cms.Sequence(
220     #akPu1PFJetID * akPu2PFJetID *akPu3PFJetID * akPu4PFJetID * akPu5PFJetID * akPu6PFJetID *
221     # ak1PFJetID * ak2PFJetID *ak3PFJetID * ak4PFJetID * ak5PFJetID * ak6PFJetID *
222     akPu2CaloJetID *akPu3CaloJetID * akPu4CaloJetID * akPu5CaloJetID *
223     ak2CaloJetID *ak3CaloJetID * ak4CaloJetID * ak5CaloJetID
224     )
225    
226 yilmaz 1.17 recoFastJets = cms.Sequence(ak1PFJets05+ak2PFJets05+ak3PFJets05+ak4PFJets05+ak5PFJets05+ak6PFJets05
227     +
228     ak1PFJets1+ak2PFJets1+ak3PFJets1+ak4PFJets1+ak5PFJets1+ak6PFJets1
229     +
230     ak1PFJets2+ak2PFJets2+ak3PFJets2+ak4PFJets2+ak5PFJets2+ak6PFJets2
231     +
232     kt1PFJets05+kt2PFJets05+kt3PFJets05+kt4PFJets05+kt5PFJets05+kt6PFJets05
233     +
234     kt1PFJets1+kt2PFJets1+kt3PFJets1+kt4PFJets1+kt5PFJets1+kt6PFJets1
235     +
236     kt1PFJets2+kt2PFJets2+kt3PFJets2+kt4PFJets2+kt5PFJets2+kt6PFJets2
237     )
238    
239 yilmaz 1.13 recoJetsWithID = cms.Sequence(recoAk1to6*recoAk1to6ID)
240 yjlee 1.14 recoJetsWithID2to5 = cms.Sequence(recoAk2to5*recoAk2to5ID)
241 yilmaz 1.13
242 yilmaz 1.4
243    
244    
245    
246    
247    
248