ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CMGWPrimeGroup/python/selectors_cff.py
Revision: 1.21
Committed: Fri Jan 25 17:57:24 2013 UTC (12 years, 3 months ago) by fantasia
Content type: text/x-python
Branch: MAIN
Changes since 1.20: +16 -7 lines
Log Message:
print outs

File Contents

# User Rev Content
1 fantasia 1.1 import FWCore.ParameterSet.Config as cms
2 fantasia 1.5
3     #####################
4     ###### Muons ######
5     #####################
6    
7 fantasia 1.1 #Note: These cuts are <= and >=
8     muonSelectors = cms.PSet(
9 fantasia 1.5 VBTF = cms.PSet(
10 fantasia 1.1 minPt = cms.untracked.double(10.),
11     maxEta = cms.untracked.double(2.4),
12     minIsGlobal = cms.untracked.int32(1),
13     maxDxy = cms.untracked.double(0.2),
14     maxNormalizedChi2 = cms.untracked.double(10.0),
15     minNTrackerHits = cms.untracked.int32(11),
16     minNPixelHits = cms.untracked.int32(1),
17     minNMuonHits = cms.untracked.int32(1),
18     minNMatches = cms.untracked.int32(2),
19     ),
20 fantasia 1.17 PFIso = cms.PSet(
21     minPt = cms.untracked.double(10.),
22 fantasia 1.21 maxEta = cms.untracked.double(2.2),
23     # maxEta = cms.untracked.double(2.4),
24 fantasia 1.20 minIsGblOrTrk = cms.untracked.int32(1),
25 fantasia 1.21 minIsGlobal = cms.untracked.int32(1), #Cory: Remove for highpt
26     ## minIsPF = cms.untracked.int32(1),
27 fantasia 1.17 maxDxy = cms.untracked.double(0.2),
28 fantasia 1.19 maxDz = cms.untracked.double(0.5),
29 fantasia 1.21 maxNormalizedChi2 = cms.untracked.double(10.0), #Cory: Remove for highpt
30     maxRelPtErr = cms.untracked.double(0.3),
31 fantasia 1.17 minNPixelHits = cms.untracked.int32(1),
32 fantasia 1.21 minNMuonHits = cms.untracked.int32(1), #Cory: This causes problems at high pt
33 fantasia 1.19 minNTrackerLayers = cms.untracked.int32(6),
34     #minNTrackerLayers = cms.untracked.int32(9),
35     #minNTrackerHits = cms.untracked.int32(11),
36 fantasia 1.17 minNMatches = cms.untracked.int32(2),
37     #maxPFIso = cms.untracked.double(???),#Added below in clones
38     ),
39 fantasia 1.1 WZLoose = cms.PSet(),
40     WZRelaxed = cms.PSet(),
41     WZTight = cms.PSet(),
42 fantasia 1.5 HadVZLoose = cms.PSet(),
43     HadVZTight = cms.PSet(),
44 fantasia 1.1 exotica = cms.PSet(
45 fantasia 1.9 maxEta = cms.untracked.double(2.4),
46     minIsGlobal = cms.untracked.int32(1),
47     minIsTracker = cms.untracked.int32(1),
48     maxDxy = cms.untracked.double(0.2),
49 cleonido 1.10 # maxNormalizedChi2 = cms.untracked.double(10.), not needed for TeV muons
50     # minNTrackerHits = cms.untracked.int32(11), deprecated
51 fantasia 1.9 minNPixelHits = cms.untracked.int32(1),
52     minNMuonHits = cms.untracked.int32(1),
53     minNMatches = cms.untracked.int32(2),
54 cleonido 1.10 minNTrackerLayers = cms.untracked.int32(9),
55 fantasia 1.9 # minTrackerValidFrac = cms.untracked.double(),
56 fantasia 1.1 ),
57     )
58 fantasia 1.17 ###WprimeWZ
59 fantasia 1.9 muonSelectors.WZRelaxed = muonSelectors.VBTF.clone(
60 fantasia 1.18 )
61 fantasia 1.19 muonSelectors.WZRelaxedPt20 = muonSelectors.WZRelaxed.clone(
62 fantasia 1.9 minPt = 20.
63     )
64 fantasia 1.19 muonSelectors.WZLoose = muonSelectors.WZRelaxed.clone(
65     maxIso03 = cms.untracked.double(0.15)
66     )
67     muonSelectors.WZTight = muonSelectors.WZRelaxed.clone(
68 fantasia 1.18 maxIso03 = cms.untracked.double(0.1)
69     )
70 fantasia 1.19 muonSelectors.WZTightPt20 = muonSelectors.WZTight.clone(
71 fantasia 1.9 minPt = 20.,
72     )
73 fantasia 1.17 ###EWKWZ
74     muonSelectors.EWKWZRelaxed = muonSelectors.PFIso.clone(
75 fantasia 1.18 )
76 fantasia 1.19 muonSelectors.EWKWZRelaxedPt20 = muonSelectors.EWKWZRelaxed.clone(
77 fantasia 1.17 minPt = 20.
78     )
79 fantasia 1.19 muonSelectors.EWKWZLoose = muonSelectors.EWKWZRelaxed.clone(
80     maxPFIso = cms.untracked.double(0.2)
81 fantasia 1.20 )
82     muonSelectors.EWKWZLoosePt20 = muonSelectors.EWKWZLoose.clone(
83     minPt = 20.,
84 fantasia 1.19 )
85     muonSelectors.EWKWZTight = muonSelectors.EWKWZRelaxed.clone(
86 fantasia 1.18 maxPFIso = cms.untracked.double(0.12)
87     )
88 fantasia 1.19 muonSelectors.EWKWZTightPt20 = muonSelectors.EWKWZTight.clone(
89 fantasia 1.17 minPt = 20.,
90     )
91 fantasia 1.18
92 fantasia 1.17 ###HadVZ
93 fantasia 1.11 muonSelectors.HadVZLoose = muonSelectors.exotica.clone(
94 fantasia 1.12 minPt = cms.untracked.double(20.),
95 fladias 1.16 # maxIso03 = cms.untracked.double(0.1)
96 fantasia 1.9 )
97 fantasia 1.11 muonSelectors.HadVZTight = muonSelectors.exotica.clone(
98 fantasia 1.12 minPt = cms.untracked.double(35.),
99 fladias 1.16 # maxIso03 = cms.untracked.double(0.1)
100 fantasia 1.9 )
101 fantasia 1.5
102     #####################
103     #### Electrons ####
104     #####################
105    
106 fantasia 1.1 cutsMissingHits = [0, 0, 0, 0, 0, 0]
107     cutsConvDist = [0., 0., 0.02, 0.02, 0.02, 0.02]
108     cutsConvDcot = [0., 0., 0.02, 0.02, 0.02, 0.02]
109 fantasia 1.2 cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.070, 0.030, 0.016]#2010
110     #cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.040, 0.030, 0.016]#2011
111 fantasia 1.1 cutsEBSigmaIEtaIEta = [0.012, 0.01, 0.01, 0.01, 0.01, 0.01]
112     cutsEBDeltaPhi = [0.800, 0.071, 0.039, 0.027, 0.020, 0.020]
113     cutsEBDeltaEta = [0.007, 0.007, 0.005, 0.005, 0.004, 0.004]
114     #cutsEECombRelIso = [0.100, 0.051, 0.042, 0.033, 0.016, 0.008]#2011
115     cutsEECombRelIso = [0.100, 0.051, 0.042, 0.060, 0.016, 0.008]#2010
116     cutsEESigmaIEtaIEta = [0.031, 0.031, 0.031, 0.031, 0.031, 0.031]
117     cutsEEDeltaPhi = [0.7, 0.047, 0.028, 0.021, 0.021, 0.021]
118     cutsEEDeltaEta = [0.011, 0.011, 0.007, 0.006, 0.005, 0.004]
119    
120     electronSelectors = cms.PSet(
121     WZLoose = cms.PSet(),
122 fantasia 1.18 WZRelaxed95 = cms.PSet(),
123     WZRelaxed80 = cms.PSet(),
124 fantasia 1.1 WZTight = cms.PSet(),
125 fantasia 1.19 exotica = cms.PSet(
126     joint = cms.PSet(
127     minPt = cms.untracked.double(10.),
128     minIsEcalDriven = cms.untracked.int32(1),
129     maxDeltaPhi = cms.untracked.double(0.06),
130     maxHoverE = cms.untracked.double(0.05),
131     #maxTrackIso = cms.untracked.double(5.),
132     maxPFIso = cms.untracked.double(0.15),
133     maxMissingHits = cms.untracked.int32(1),
134     ),
135     barrel = cms.PSet(
136     maxDeltaEta = cms.untracked.double(0.005),
137     #maxSigmaIEtaIEta = cms.untracked.double(),
138     minPassEX5overE55 = cms.untracked.int32(1),
139     #minPassEMHadDepth1Iso = cms.untracked.int32(1),
140     maxd0 = cms.untracked.double(0.02),
141     ),
142     endcap = cms.PSet(
143     maxDeltaEta = cms.untracked.double(0.007),
144     maxSigmaIEtaIEta = cms.untracked.double(0.03),
145     minPassEX5overE55 = cms.untracked.int32(1),
146     #minPassEMHadDepth1Iso = cms.untracked.int32(1),
147     maxd0 = cms.untracked.double(0.05),
148     ),
149     ),
150 fantasia 1.18 mva = cms.PSet(
151     joint = cms.PSet(
152     minPt = cms.untracked.double(10.),
153     minpassMVAPresel = cms.untracked.int32(1),
154     minpassMVATrig = cms.untracked.int32(1),
155     #maxPFIso = cms.untracked.double(???),#Added below in clones
156     ),
157     barrel = cms.PSet(
158     ),
159     endcap = cms.PSet(
160     ),
161     ),
162     EWKWZLoose = cms.PSet(),
163     EWKWZTight = cms.PSet(),
164 fantasia 1.19
165     CiC2012Loose = cms.PSet(
166     joint = cms.PSet(
167     minPt = cms.untracked.double(10.),
168     ),
169     barrel = cms.PSet(
170     maxMissingHits = cms.untracked.int32(1),
171     minVtxFitProb = cms.untracked.double(1e-6),
172     maxSigmaIEtaIEta = cms.untracked.double(0.01),
173     maxDeltaPhi = cms.untracked.double(0.15),
174     maxDeltaEta = cms.untracked.double(0.007),
175     maxHoverE = cms.untracked.double(0.12),
176     maxd0 = cms.untracked.double(0.02),
177 fantasia 1.20 #maxdz = cms.untracked.double(200.),
178     #maxdz = cms.untracked.double(0.2),#proper cut
179 fantasia 1.19 maxfabsdiffEp = cms.untracked.double(0.05),
180     maxPFIso = cms.untracked.double(0.15),
181     ),
182     endcap = cms.PSet(
183     maxMissingHits = cms.untracked.int32(1),
184     minVtxFitProb = cms.untracked.double(1e-6),
185     maxSigmaIEtaIEta = cms.untracked.double(0.03),
186     maxDeltaPhi = cms.untracked.double(0.10),
187     maxDeltaEta = cms.untracked.double(0.009),
188     maxHoverE = cms.untracked.double(0.10),
189     maxd0 = cms.untracked.double(0.02),
190 fantasia 1.20 #maxdz = cms.untracked.double(200.),
191     #maxdz = cms.untracked.double(0.2),#proper cut
192 fantasia 1.19 maxfabsdiffEp = cms.untracked.double(0.05),
193     maxPFIso = cms.untracked.double(0.15),
194     ),
195     ),
196 fantasia 1.21 CiC2012LoosePt20 = cms.PSet(),
197 fantasia 1.19 CiC2012Medium = cms.PSet(
198     joint = cms.PSet(
199 fantasia 1.21 minPt = cms.untracked.double(10.),
200 fantasia 1.19 ),
201     barrel = cms.PSet(
202     maxMissingHits = cms.untracked.int32(1),
203     minVtxFitProb = cms.untracked.double(1e-6),
204     maxSigmaIEtaIEta = cms.untracked.double(0.01),
205     maxDeltaPhi = cms.untracked.double(0.06),
206     maxDeltaEta = cms.untracked.double(0.004),
207     maxHoverE = cms.untracked.double(0.12),
208     maxd0 = cms.untracked.double(0.02),
209     #maxdz = cms.untracked.double(0.1),
210     maxfabsdiffEp = cms.untracked.double(0.05),
211     maxPFIso = cms.untracked.double(0.15),
212     ),
213     endcap = cms.PSet(
214     maxMissingHits = cms.untracked.int32(1),
215     minVtxFitProb = cms.untracked.double(1e-6),
216     maxSigmaIEtaIEta = cms.untracked.double(0.03),
217     maxDeltaPhi = cms.untracked.double(0.03),
218     maxDeltaEta = cms.untracked.double(0.007),
219     maxHoverE = cms.untracked.double(0.10),
220     maxd0 = cms.untracked.double(0.02),
221     #maxdz = cms.untracked.double(0.1),
222     maxfabsdiffEp = cms.untracked.double(0.05),
223     maxPFIso = cms.untracked.double(0.15),
224     ),
225     ),
226 fantasia 1.21 CiC2012MediumPt20 = cms.PSet(),
227 fantasia 1.1 )
228    
229     for i, s in enumerate(["wp95", "wp90", "wp85", "wp80", "wp70", "wp60"]):
230     pset = cms.PSet(
231 fantasia 1.18 joint = cms.PSet(
232     minPt = cms.untracked.double(10.),
233     ),
234 fantasia 1.1 barrel = cms.PSet(
235     maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
236     minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an OR of these two cuts below
237     maxCombRelIso = cms.untracked.double(cutsEBCombRelIso[i]),
238     maxSigmaIEtaIEta = cms.untracked.double(cutsEBSigmaIEtaIEta[i]),
239     maxDeltaPhi = cms.untracked.double(cutsEBDeltaPhi[i]),
240     maxDeltaEta = cms.untracked.double(cutsEBDeltaEta[i]),
241     ),
242     endcap = cms.PSet(
243     maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
244     minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an or of these two cuts below
245     maxCombRelIso = cms.untracked.double(cutsEECombRelIso[i]),
246     maxSigmaIEtaIEta = cms.untracked.double(cutsEESigmaIEtaIEta[i]),
247     maxDeltaPhi = cms.untracked.double(cutsEEDeltaPhi[i]),
248     maxDeltaEta = cms.untracked.double(cutsEEDeltaEta[i]),
249     ),
250     )
251     setattr(electronSelectors, s, pset)
252    
253     electronSelectors.WZLoose = electronSelectors.wp95.clone()
254 fantasia 1.9
255 fantasia 1.18 electronSelectors.WZRelaxed95 = electronSelectors.wp95.clone()
256     electronSelectors.WZRelaxed95.barrel.maxCombRelIso = cms.untracked.double(999999999.)
257     electronSelectors.WZRelaxed95.endcap.maxCombRelIso = cms.untracked.double(999999999.)
258    
259     electronSelectors.WZRelaxed80 = electronSelectors.wp80.clone()
260     electronSelectors.WZRelaxed80.barrel.maxCombRelIso = cms.untracked.double(999999999.)
261     electronSelectors.WZRelaxed80.endcap.maxCombRelIso = cms.untracked.double(999999999.)
262    
263     electronSelectors.WZRelaxed80Pt20 = electronSelectors.wp80.clone()
264     electronSelectors.WZRelaxed80Pt20.joint.minPt = cms.untracked.double(20.)
265     electronSelectors.WZRelaxed80Pt20.barrel.maxCombRelIso = cms.untracked.double(999999999.)
266     electronSelectors.WZRelaxed80Pt20.endcap.maxCombRelIso = cms.untracked.double(999999999.)
267 fantasia 1.1
268     electronSelectors.WZTight = electronSelectors.wp80.clone()
269    
270 fantasia 1.18 electronSelectors.WZTightPt20 = electronSelectors.wp80.clone()
271     electronSelectors.WZTightPt20.joint.minPt = cms.untracked.double(20.)
272    
273     ##EWKWZ
274     electronSelectors.EWKWZLoose = electronSelectors.mva.clone()
275     electronSelectors.EWKWZLoose.joint.maxPFIso = cms.untracked.double(0.1)
276    
277     electronSelectors.EWKWZRelaxed = electronSelectors.mva.clone()
278    
279     electronSelectors.EWKWZRelaxedPt20 = electronSelectors.mva.clone()
280     electronSelectors.EWKWZRelaxedPt20.joint.minPt = 20.
281    
282     electronSelectors.EWKWZTight = electronSelectors.mva.clone()
283     electronSelectors.EWKWZTight.joint.maxPFIso = cms.untracked.double(0.1)
284    
285     electronSelectors.EWKWZTightPt20 = electronSelectors.mva.clone()
286     electronSelectors.EWKWZTightPt20.joint.minPt = 20.
287     electronSelectors.EWKWZTightPt20.joint.maxPFIso = cms.untracked.double(0.1)
288 fantasia 1.9
289 fantasia 1.21 #CiC2012
290     electronSelectors.CiC2012LoosePt20 = electronSelectors.CiC2012Loose.clone()
291     electronSelectors.CiC2012LoosePt20.joint.minPt = 20.
292    
293     electronSelectors.CiC2012MediumPt20 = electronSelectors.CiC2012Medium.clone()
294     electronSelectors.CiC2012MediumPt20.joint.minPt = 20.
295 fantasia 1.19
296     ##Exotica
297     electronSelectors.exoticaPt20 = electronSelectors.exotica.clone()
298     electronSelectors.exoticaPt20.joint.minPt = cms.untracked.double(20.)
299    
300    
301    
302 fantasia 1.4 ####################
303     ##### Jets #######
304     ####################
305     jetSelectors = cms.PSet(
306 fantasia 1.8 Base = cms.PSet(
307 fladias 1.14 minPt = cms.untracked.double(30.),
308 fantasia 1.4 maxEta = cms.untracked.double(2.4),
309     maxNHF = cms.untracked.double(0.99),
310     maxNEF = cms.untracked.double(0.99),
311 fantasia 1.6 minNDaughters = cms.untracked.int32(1),
312     minCHF = cms.untracked.double(0.0),
313 fantasia 1.4 maxCEF = cms.untracked.double(0.99),
314 fantasia 1.6 minCMult = cms.untracked.int32(0),
315 fantasia 1.4 ),
316 fantasia 1.17 Pat = cms.PSet(
317     minPt = cms.untracked.double(30.),
318     maxEta = cms.untracked.double(3.0),
319     ),
320 fantasia 1.4 )
321 fantasia 1.8
322     ####################
323     #### Photons #######
324     ####################
325     photonSelectors = cms.PSet(
326     Base = cms.PSet(#Made up values
327     minPt = cms.untracked.double(10.),
328     maxEta = cms.untracked.double(5),
329     maxECalIso = cms.untracked.double(0.1),
330     maxHCalIso = cms.untracked.double(0.1),
331     maxTrkIso = cms.untracked.double(0.1),
332     maxHoE = cms.untracked.double(0.05),
333     maxSigmaee = cms.untracked.double(9e9),
334     minHasSeed = cms.untracked.bool(True),
335     ),
336     )
337