CMS 3D CMS Logo
singleTopDQM_miniAOD_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 looseMuonCut = " (isGlobalMuon || isTrackerMuon) && isPFMuon"
4 looseIsoCut = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.25"
5 
6 tightMuonCut = " isGlobalMuon && isPFMuon && globalTrack.normalizedChi2 < 10. && globalTrack.hitPattern.numberOfValidMuonHits > 0 && " + \
7  "numberOfMatchedStations > 1 && innerTrack.hitPattern.numberOfValidPixelHits > 0 && innerTrack.hitPattern.trackerLayersWithMeasurement > 5"
8 # CB PV cut!
9 tightIsoCut = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.15"
10 
11 
12 EletightIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.1"
13 ElelooseIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.15"
14 
15 
16 looseElecCut = "((full5x5_sigmaIetaIeta < 0.011 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00477 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.222 && hadronicOverEm < 0.298 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.241 && gsfTrack.hitPattern.numberOfHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) < 1.479) || (full5x5_sigmaIetaIeta() < 0.0314 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00868 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.213 && hadronicOverEm < 0.101 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.14 && gsfTrack.hitPattern.numberOfHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) > 1.479))"
17 
18 elecIPcut = "(abs(gsfTrack.d0)<0.05 & abs(gsfTrack.dz)<0.1 & abs(superCluster.eta) < 1.479)||(abs(gsfTrack.d0)<0.1 && abs(gsfTrack.dz)<0.2 && abs(superCluster.eta) > 1.479)"
19 
20 
21 tightElecCut = "((full5x5_sigmaIetaIeta < 0.00998 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00308 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.0816 && hadronicOverEm < 0.0414 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.0129 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) < 1.479) || (full5x5_sigmaIetaIeta() < 0.0292 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00605 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.0394 && hadronicOverEm < 0.0641 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.0129 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) > 1.479))"
22 
23 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
24 
25 singleTopMuonMediumDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
26  setup = cms.PSet(
27 
28  directory = cms.string("Physics/Top/SingleTopMuonMediumDQM_miniAOD/"),
29  sources = cms.PSet(
30  muons = cms.InputTag("slimmedMuons"),
31  elecs = cms.InputTag("slimmedElectrons"),
32  jets = cms.InputTag("slimmedJets"),
33  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
34  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
35  ),
36 
37  monitoring = cms.PSet(
38  verbosity = cms.string("DEBUG")
39  ),
40 
41  pvExtras = cms.PSet(
42  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
43  ),
44  elecExtras = cms.PSet(
45  select = cms.string(tightElecCut + "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
46  rho = cms.InputTag("fixedGridRhoFastjetAll"),
47  #isolation = cms.string(ElelooseIsoCut),
48  ),
49 
50  muonExtras = cms.PSet(
51  select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
52  isolation = cms.string(looseIsoCut)
53  ),
54  jetExtras = cms.PSet(
55 
56  select = cms.string("pt>30 & abs(eta)<2.4"), # & neutralEmEnergyFraction >0.01 & chargedEmEnergyFraction>0.01"),
57  jetBTaggers = cms.PSet(
58  trackCountingEff = cms.PSet(
59  label = cms.InputTag("trackCountingHighEffBJetTags" ),
60  workingPoint = cms.double(1.25)
61  ),
62  trackCountingPur = cms.PSet(
63  label = cms.InputTag("trackCountingHighPurBJetTags" ),
64  workingPoint = cms.double(3.41)
65  ),
66  secondaryVertex = cms.PSet(
67  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
68  workingPoint = cms.double(2.05)
69  ),
70  combinedSecondaryVertex = cms.PSet(
71  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
72  workingPoint = cms.double(0.898)
73  )
74  ),
75  ),
76  massExtras = cms.PSet(
77  lowerEdge = cms.double( 70.),
78  upperEdge = cms.double(110.)
79  ),
80  ),
81  preselection = cms.PSet(
82  vertex = cms.PSet(
83  src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
84  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
85  )
86  ),
87 
88  selection = cms.VPSet(
89  cms.PSet(
90  label = cms.string("muons:step0"),
91  src = cms.InputTag("slimmedMuons"),
92  select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut), #tightMuonCut +"&&"+ tightIsoCut + " && pt>20 & abs(eta)<2.1"), # CB what about iso? CD Added tightIso
93  min = cms.int32(1),
94  #max = cms.int32(1),
95  ),
96  cms.PSet(
97  label = cms.string("jets:step1"),
98  src = cms.InputTag("slimmedJets"),
99  select = cms.string("pt>30 & abs(eta)<2.4 "),
100  min = cms.int32(2),
101  ),
102  cms.PSet(
103  label = cms.string("met:step2"),
104  src = cms.InputTag("slimmedMETs"),
105  select = cms.string("pt>30"),
106  #min = cms.int32(2),
107  ),
108  )
109 )
110 
111 singleTopElectronMediumDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
112 
113  setup = cms.PSet(
114 
115  directory = cms.string("Physics/Top/SingleTopElectronMediumDQM_miniAOD/"),
116 
117  sources = cms.PSet(
118  muons = cms.InputTag("slimmedMuons"),
119  elecs = cms.InputTag("slimmedElectrons"),
120  jets = cms.InputTag("slimmedJets"),
121  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
122  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
123 
124  ),
125 
126  monitoring = cms.PSet(
127  verbosity = cms.string("DEBUG")
128  ),
129 
130  pvExtras = cms.PSet(
131  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
132  ),
133  elecExtras = cms.PSet(
134  select = cms.string(tightElecCut + " && pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
135  #select = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
136  rho = cms.InputTag("fixedGridRhoFastjetAll"),
137  #isolation = cms.string(ElelooseIsoCut),
138  ),
139  muonExtras = cms.PSet(
140  select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
141  isolation = cms.string(looseIsoCut)
142  ),
143 
144  jetExtras = cms.PSet(
145 
146  select = cms.string("pt>30 & abs(eta)<2.4"),
147  jetBTaggers = cms.PSet(
148  trackCountingEff = cms.PSet(
149  label = cms.InputTag("trackCountingHighEffBJetTags" ),
150  workingPoint = cms.double(1.25)
151  ),
152  trackCountingPur = cms.PSet(
153  label = cms.InputTag("trackCountingHighPurBJetTags" ),
154  workingPoint = cms.double(3.41)
155  ),
156  secondaryVertex = cms.PSet(
157  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
158  workingPoint = cms.double(2.05)
159  ),
160  combinedSecondaryVertex = cms.PSet(
161  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
162  workingPoint = cms.double(0.898)
163  ),
164  ),
165  ),
166  massExtras = cms.PSet(
167  lowerEdge = cms.double( 70.),
168  upperEdge = cms.double(110.)
169  ),
170 
171  ),
172 
173  preselection = cms.PSet(
174  vertex = cms.PSet(
175  src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
176  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
177  )
178  ),
179 
180  selection = cms.VPSet(
181  cms.PSet(
182  label = cms.string("elecs:step0"),
183  src = cms.InputTag("slimmedElectrons"),
184  select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) &&" + tightElecCut),
185  min = cms.int32(1),
186  max = cms.int32(1),
187  ),
188  cms.PSet(
189  label = cms.string("jets:step1"),
190  src = cms.InputTag("slimmedJets"),
191  select = cms.string("pt>30 & abs(eta)<2.4 "),
192  min = cms.int32(2),
193  ),
194  cms.PSet(
195  label = cms.string("met:step2"),
196  src = cms.InputTag("slimmedMETs"),
197  select = cms.string("pt>30"),
198  #min = cms.int32(2),
199  ),
200  )
201 )
DQMEDAnalyzer
Definition: DQMEDAnalyzer.py:1
Generated for CMSSW Reference Manual by doxygen 1.8.16