CMGWPrimeGroup/python/selectors_cff.py

import FWCore.ParameterSet.Config as cms

#####################
######  Muons  ######
#####################

#Note: These cuts are <= and >= 
muonSelectors = cms.PSet(
    VBTF = cms.PSet(
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(2.4),
       minIsGlobal = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
       maxNormalizedChi2 = cms.untracked.double(10.0),
       minNTrackerHits = cms.untracked.int32(11),
       minNPixelHits = cms.untracked.int32(1),
       minNMuonHits = cms.untracked.int32(1),
       minNMatches = cms.untracked.int32(2),
    ),
    WZLoose = cms.PSet(),
    WZRelaxed = cms.PSet(),
    WZTight = cms.PSet(),
    HadVZLoose = cms.PSet(),
    HadVZTight = cms.PSet(),
    exotica = cms.PSet(
       maxEta = cms.untracked.double(2.4),
       minIsGlobal = cms.untracked.int32(1),
       minIsTracker = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
#       maxNormalizedChi2 = cms.untracked.double(10.), not needed for TeV muons
#       minNTrackerHits = cms.untracked.int32(11), deprecated
       minNPixelHits = cms.untracked.int32(1),
       minNMuonHits = cms.untracked.int32(1),
       minNMatches = cms.untracked.int32(2),
       minNTrackerLayers = cms.untracked.int32(9),
#       minTrackerValidFrac = cms.untracked.double(),
    ),
    )
muonSelectors.WZLoose = muonSelectors.VBTF.clone(
    minPt = 10.,
    maxIso03 = cms.untracked.double(0.15)
    )
muonSelectors.WZRelaxed = muonSelectors.VBTF.clone(
    minPt = 20.
    )
muonSelectors.WZTight = muonSelectors.VBTF.clone(
    minPt = 20.,
    maxIso03 = cms.untracked.double(0.1)
    )
muonSelectors.HadVZLoose = muonSelectors.exotica.clone(
    minPt = cms.untracked.double(20.),
    )
muonSelectors.HadVZTight = muonSelectors.exotica.clone(
    minPt = cms.untracked.double(35.),
    )
#muonSelectors.HadVZTight.remove(maxNormalizedChi2)

#####################
####  Electrons  ####
#####################

cutsMissingHits = [0, 0, 0, 0, 0, 0]
cutsConvDist = [0., 0., 0.02, 0.02, 0.02, 0.02]
cutsConvDcot = [0., 0., 0.02, 0.02, 0.02, 0.02]
cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.070, 0.030, 0.016]#2010
#cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.040, 0.030, 0.016]#2011
cutsEBSigmaIEtaIEta = [0.012, 0.01, 0.01, 0.01, 0.01, 0.01]
cutsEBDeltaPhi = [0.800, 0.071, 0.039, 0.027, 0.020, 0.020]
cutsEBDeltaEta = [0.007, 0.007, 0.005, 0.005, 0.004, 0.004]
#cutsEECombRelIso = [0.100, 0.051, 0.042, 0.033, 0.016, 0.008]#2011
cutsEECombRelIso = [0.100, 0.051, 0.042, 0.060, 0.016, 0.008]#2010
cutsEESigmaIEtaIEta = [0.031, 0.031, 0.031, 0.031, 0.031, 0.031]
cutsEEDeltaPhi = [0.7, 0.047, 0.028, 0.021, 0.021, 0.021]
cutsEEDeltaEta = [0.011, 0.011, 0.007, 0.006, 0.005, 0.004]

electronSelectors = cms.PSet(
    WZLoose = cms.PSet(),
    WZRelaxed = cms.PSet(),
    WZTight = cms.PSet(),
    exotica = cms.PSet(),
    )

for i, s in enumerate(["wp95", "wp90", "wp85", "wp80", "wp70", "wp60"]):
    pset = cms.PSet(
        barrel = cms.PSet(
           maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
           minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an OR of these two cuts below
           maxCombRelIso = cms.untracked.double(cutsEBCombRelIso[i]),
           maxSigmaIEtaIEta = cms.untracked.double(cutsEBSigmaIEtaIEta[i]),
           maxDeltaPhi = cms.untracked.double(cutsEBDeltaPhi[i]),
           maxDeltaEta = cms.untracked.double(cutsEBDeltaEta[i]),
           ),
        endcap = cms.PSet(
           maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
           minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an or of these two cuts below
           maxCombRelIso = cms.untracked.double(cutsEECombRelIso[i]),
           maxSigmaIEtaIEta = cms.untracked.double(cutsEESigmaIEtaIEta[i]),
           maxDeltaPhi = cms.untracked.double(cutsEEDeltaPhi[i]),
           maxDeltaEta = cms.untracked.double(cutsEEDeltaEta[i]),
           ),
        )
    setattr(electronSelectors, s, pset)
    
electronSelectors.WZLoose = electronSelectors.wp95.clone()
electronSelectors.WZLoose.barrel.minPt = cms.untracked.double(10.)
electronSelectors.WZLoose.endcap.minPt = cms.untracked.double(10.)

electronSelectors.WZRelaxed = electronSelectors.wp80.clone()
electronSelectors.WZRelaxed.barrel.minPt = cms.untracked.double(20.)
electronSelectors.WZRelaxed.endcap.minPt = cms.untracked.double(20.)
electronSelectors.WZRelaxed.barrel.maxCombRelIso = cms.untracked.double(999999999.)
electronSelectors.WZRelaxed.endcap.maxCombRelIso = cms.untracked.double(999999999.)
#electronSelectors.WZRelaxed.barrel.remove(maxCombRelIso)
#electronSelectors.WZRelaxed.endcap.remove(maxCombRelIso)

electronSelectors.WZTight = electronSelectors.wp80.clone()
electronSelectors.WZTight.barrel.minPt = cms.untracked.double(20.)
electronSelectors.WZTight.endcap.minPt = cms.untracked.double(20.)


####################
#####  Jets  #######
####################
jetSelectors = cms.PSet(
    Base = cms.PSet(
       minPt = cms.untracked.double(30.),
       maxEta = cms.untracked.double(2.4),
       maxNHF = cms.untracked.double(0.99),
       maxNEF = cms.untracked.double(0.99),
       minNDaughters = cms.untracked.int32(1),
       minCHF = cms.untracked.double(0.0),
       maxCEF = cms.untracked.double(0.99),
       minCMult = cms.untracked.int32(0),
       ),
    )

####################
#### Photons #######
####################
photonSelectors = cms.PSet(
    Base = cms.PSet(#Made up values
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(5),
       maxECalIso = cms.untracked.double(0.1),
       maxHCalIso = cms.untracked.double(0.1),
       maxTrkIso = cms.untracked.double(0.1),
       maxHoE = cms.untracked.double(0.05),
       maxSigmaee = cms.untracked.double(9e9),
       minHasSeed = cms.untracked.bool(True),       
       ),
    )

Revision:	1.12
Committed:	Tue Nov 8 12:29:29 2011 UTC (13 years, 5 months ago) by fantasia
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	V00-03-60
Changes since 1.11:	+2 -4 lines
Log Message:	moved VZ to TeV exotica muon selectors
#	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			WZLoose = cms.PSet(),
21			WZRelaxed = cms.PSet(),
22			WZTight = cms.PSet(),
23	fantasia	1.5	HadVZLoose = cms.PSet(),
24			HadVZTight = cms.PSet(),
25	fantasia	1.1	exotica = cms.PSet(
26	fantasia	1.9	maxEta = cms.untracked.double(2.4),
27			minIsGlobal = cms.untracked.int32(1),
28			minIsTracker = cms.untracked.int32(1),
29			maxDxy = cms.untracked.double(0.2),
30	cleonido	1.10	# maxNormalizedChi2 = cms.untracked.double(10.), not needed for TeV muons
31			# minNTrackerHits = cms.untracked.int32(11), deprecated
32	fantasia	1.9	minNPixelHits = cms.untracked.int32(1),
33			minNMuonHits = cms.untracked.int32(1),
34			minNMatches = cms.untracked.int32(2),
35	cleonido	1.10	minNTrackerLayers = cms.untracked.int32(9),
36	fantasia	1.9	# minTrackerValidFrac = cms.untracked.double(),
37	fantasia	1.1	),
38			)
39	fantasia	1.9	muonSelectors.WZLoose = muonSelectors.VBTF.clone(
40			minPt = 10.,
41			maxIso03 = cms.untracked.double(0.15)
42			)
43			muonSelectors.WZRelaxed = muonSelectors.VBTF.clone(
44			minPt = 20.
45			)
46			muonSelectors.WZTight = muonSelectors.VBTF.clone(
47			minPt = 20.,
48			maxIso03 = cms.untracked.double(0.1)
49			)
50	fantasia	1.11	muonSelectors.HadVZLoose = muonSelectors.exotica.clone(
51	fantasia	1.12	minPt = cms.untracked.double(20.),
52	fantasia	1.9	)
53	fantasia	1.11	muonSelectors.HadVZTight = muonSelectors.exotica.clone(
54	fantasia	1.12	minPt = cms.untracked.double(35.),
55	fantasia	1.9	)
56			#muonSelectors.HadVZTight.remove(maxNormalizedChi2)
57	fantasia	1.5
58			#####################
59			#### Electrons ####
60			#####################
61
62	fantasia	1.1	cutsMissingHits = [0, 0, 0, 0, 0, 0]
63			cutsConvDist = [0., 0., 0.02, 0.02, 0.02, 0.02]
64			cutsConvDcot = [0., 0., 0.02, 0.02, 0.02, 0.02]
65	fantasia	1.2	cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.070, 0.030, 0.016]#2010
66			#cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.040, 0.030, 0.016]#2011
67	fantasia	1.1	cutsEBSigmaIEtaIEta = [0.012, 0.01, 0.01, 0.01, 0.01, 0.01]
68			cutsEBDeltaPhi = [0.800, 0.071, 0.039, 0.027, 0.020, 0.020]
69			cutsEBDeltaEta = [0.007, 0.007, 0.005, 0.005, 0.004, 0.004]
70			#cutsEECombRelIso = [0.100, 0.051, 0.042, 0.033, 0.016, 0.008]#2011
71			cutsEECombRelIso = [0.100, 0.051, 0.042, 0.060, 0.016, 0.008]#2010
72			cutsEESigmaIEtaIEta = [0.031, 0.031, 0.031, 0.031, 0.031, 0.031]
73			cutsEEDeltaPhi = [0.7, 0.047, 0.028, 0.021, 0.021, 0.021]
74			cutsEEDeltaEta = [0.011, 0.011, 0.007, 0.006, 0.005, 0.004]
75
76			electronSelectors = cms.PSet(
77			WZLoose = cms.PSet(),
78			WZRelaxed = cms.PSet(),
79			WZTight = cms.PSet(),
80			exotica = cms.PSet(),
81			)
82
83			for i, s in enumerate(["wp95", "wp90", "wp85", "wp80", "wp70", "wp60"]):
84			pset = cms.PSet(
85			barrel = cms.PSet(
86			maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
87			minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an OR of these two cuts below
88			maxCombRelIso = cms.untracked.double(cutsEBCombRelIso[i]),
89			maxSigmaIEtaIEta = cms.untracked.double(cutsEBSigmaIEtaIEta[i]),
90			maxDeltaPhi = cms.untracked.double(cutsEBDeltaPhi[i]),
91			maxDeltaEta = cms.untracked.double(cutsEBDeltaEta[i]),
92			),
93			endcap = cms.PSet(
94			maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
95			minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an or of these two cuts below
96			maxCombRelIso = cms.untracked.double(cutsEECombRelIso[i]),
97			maxSigmaIEtaIEta = cms.untracked.double(cutsEESigmaIEtaIEta[i]),
98			maxDeltaPhi = cms.untracked.double(cutsEEDeltaPhi[i]),
99			maxDeltaEta = cms.untracked.double(cutsEEDeltaEta[i]),
100			),
101			)
102			setattr(electronSelectors, s, pset)
103
104			electronSelectors.WZLoose = electronSelectors.wp95.clone()
105	fantasia	1.3	electronSelectors.WZLoose.barrel.minPt = cms.untracked.double(10.)
106			electronSelectors.WZLoose.endcap.minPt = cms.untracked.double(10.)
107	fantasia	1.9
108	fantasia	1.1	electronSelectors.WZRelaxed = electronSelectors.wp80.clone()
109	fantasia	1.3	electronSelectors.WZRelaxed.barrel.minPt = cms.untracked.double(20.)
110			electronSelectors.WZRelaxed.endcap.minPt = cms.untracked.double(20.)
111	fantasia	1.1	electronSelectors.WZRelaxed.barrel.maxCombRelIso = cms.untracked.double(999999999.)
112			electronSelectors.WZRelaxed.endcap.maxCombRelIso = cms.untracked.double(999999999.)
113	fantasia	1.9	#electronSelectors.WZRelaxed.barrel.remove(maxCombRelIso)
114			#electronSelectors.WZRelaxed.endcap.remove(maxCombRelIso)
115	fantasia	1.1
116			electronSelectors.WZTight = electronSelectors.wp80.clone()
117	fantasia	1.3	electronSelectors.WZTight.barrel.minPt = cms.untracked.double(20.)
118			electronSelectors.WZTight.endcap.minPt = cms.untracked.double(20.)
119	fantasia	1.1
120	fantasia	1.9
121	fantasia	1.4	####################
122			##### Jets #######
123			####################
124			jetSelectors = cms.PSet(
125	fantasia	1.8	Base = cms.PSet(
126	fantasia	1.6	minPt = cms.untracked.double(30.),
127	fantasia	1.4	maxEta = cms.untracked.double(2.4),
128			maxNHF = cms.untracked.double(0.99),
129			maxNEF = cms.untracked.double(0.99),
130	fantasia	1.6	minNDaughters = cms.untracked.int32(1),
131			minCHF = cms.untracked.double(0.0),
132	fantasia	1.4	maxCEF = cms.untracked.double(0.99),
133	fantasia	1.6	minCMult = cms.untracked.int32(0),
134	fantasia	1.4	),
135			)
136	fantasia	1.8
137			####################
138			#### Photons #######
139			####################
140			photonSelectors = cms.PSet(
141			Base = cms.PSet(#Made up values
142			minPt = cms.untracked.double(10.),
143			maxEta = cms.untracked.double(5),
144			maxECalIso = cms.untracked.double(0.1),
145			maxHCalIso = cms.untracked.double(0.1),
146			maxTrkIso = cms.untracked.double(0.1),
147			maxHoE = cms.untracked.double(0.05),
148			maxSigmaee = cms.untracked.double(9e9),
149			minHasSeed = cms.untracked.bool(True),
150			),
151			)
152