Skip to content

Commit 2f29cd9

Browse files
mmusichmmusich
committed
restrict PCA z-component in the tracks matched to PV (dz<0.1)
1 parent 2e956e4 commit 2f29cd9

File tree

4 files changed

+37
-3
lines changed

4 files changed

+37
-3
lines changed

‎DQM/TrackingMonitor/python/TrackingMonitor_cfi.py‎

Lines changed: 5 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -315,7 +315,11 @@
315315
VZBin = cms.int32(100),
316316
VZMax = cms.double(30.0),
317317
VZMin = cms.double(-30.0),
318-
318+
319+
# PCA z position (to PV)
320+
VZ_PVMax = cms.double(30.0),
321+
VZ_PVMin = cms.double(-30.0),
322+
319323
# PCA z position for profile
320324
VZBinProf = cms.int32(100),
321325
VZMaxProf = cms.double(0.2),

‎DQM/TrackingMonitor/src/TrackAnalyzer.cc‎

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -362,6 +362,9 @@ void TrackAnalyzer::bookHistosForHitProperties(DQMStore::IBooker& ibooker) {
362362
double VZMin = conf_->getParameter<double>("VZMin");
363363
double VZMax = conf_->getParameter<double>("VZMax");
364364

365+
double VZ_PVMin = conf_->getParameter<double>("VZ_PVMin");
366+
double VZ_PVMax = conf_->getParameter<double>("VZ_PVMax");
367+
365368
ibooker.setCurrentFolder(TopFolder_);
366369

367370
// book the Hit Property histograms
@@ -682,7 +685,7 @@ void TrackAnalyzer::bookHistosForHitProperties(DQMStore::IBooker& ibooker) {
682685

683686
histname = "zPointOfClosestApproachToPV_";
684687
zPointOfClosestApproachToPV =
685-
ibooker.book1D(histname + CategoryName, histname + CategoryName, VZBin, VZMin, VZMax);
688+
ibooker.book1D(histname + CategoryName, histname + CategoryName, VZBin, VZ_PVMin, VZ_PVMax);
686689
zPointOfClosestApproachToPV->setAxisTitle("z component of Track PCA to pv line (cm)", 1);
687690
zPointOfClosestApproachToPV->setAxisTitle("Number of Tracks", 2);
688691
}

‎DQM/TrackingMonitorSource/python/TrackCollections2monitor_cff.py‎

Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -7,6 +7,8 @@
77
trackPtN = {}
88
trackPtMin = {}
99
trackPtMax = {}
10+
pcaZtoPVMax = {}
11+
pcaZtoPVMin = {}
1012
doPlotsPCA = {}
1113
numCutString = {}
1214
denCutString = {}
@@ -37,6 +39,8 @@
3739
trackPtN ['generalTracks'] = cms.int32(100)
3840
trackPtMin ['generalTracks'] = cms.double(0.)
3941
trackPtMax ['generalTracks'] = cms.double(100.)
42+
pcaZtoPVMax ['generalTracks'] = cms.double(30.)
43+
pcaZtoPVMin ['generalTracks'] = cms.double(-30.)
4044
doPlotsPCA ['generalTracks'] = cms.bool(False)
4145
numCutString ['generalTracks'] = cms.string("quality('highPurity')") # num := den + quality('highPurity')
4246
denCutString ['generalTracks'] = cms.string("") # den := kinematics cuts
@@ -76,6 +80,8 @@
7680
trackPtN ['highPurityPtRange0to1'] = cms.int32(10)
7781
trackPtMin ['highPurityPtRange0to1'] = cms.double(0.)
7882
trackPtMax ['highPurityPtRange0to1'] = cms.double(1.)
83+
pcaZtoPVMax ['highPurityPtRange0to1'] = cms.double(30.)
84+
pcaZtoPVMin ['highPurityPtRange0to1'] = cms.double(-30.)
7985
numCutString ['highPurityPtRange0to1'] = cms.string(" pt >= 0 & pt < 1 & quality('highPurity')") # num := den + quality('highPurity') [it is the same as the main selection, but just to be sure ...]
8086
denCutString ['highPurityPtRange0to1'] = cms.string(" pt >= 0 & pt < 1 ") # den := kinematics cut
8187
doPlotsPCA ['highPurityPtRange0to1'] = cms.bool(False)
@@ -109,6 +115,8 @@
109115
trackPtN ['highPurityPtRange1to10'] = cms.int32(10)
110116
trackPtMin ['highPurityPtRange1to10'] = cms.double(1.)
111117
trackPtMax ['highPurityPtRange1to10'] = cms.double(10.)
118+
pcaZtoPVMax ['highPurityPtRange1to10'] = cms.double(30.)
119+
pcaZtoPVMin ['highPurityPtRange1to10'] = cms.double(-30.)
112120
numCutString ['highPurityPtRange1to10'] = cms.string(" pt >= 1 & pt < 10 & quality('highPurity')") # num := den + quality('highPurity') [it is the same as the main selection, but just to be sure ...]
113121
denCutString ['highPurityPtRange1to10'] = cms.string(" pt >= 1 & pt < 10 ") # den := kinematics cuts
114122
doGoodTracksPlots ['highPurityPtRange1to10'] = cms.bool(True)
@@ -141,6 +149,8 @@
141149
trackPtN ['highPurityPt10'] = cms.int32(100)
142150
trackPtMin ['highPurityPt10'] = cms.double(10.)
143151
trackPtMax ['highPurityPt10'] = cms.double(110.)
152+
pcaZtoPVMax ['highPurityPt10'] = cms.double(30.)
153+
pcaZtoPVMin ['highPurityPt10'] = cms.double(-30.)
144154
numCutString ['highPurityPt10'] = cms.string(" pt >= 10 & quality('highPurity')") # num := den + quality('highPurity') [it is the same as the main selection, but just to be sure ...]
145155
denCutString ['highPurityPt10'] = cms.string(" pt >= 10 ") # den := kinematics cuts
146156
doGoodTracksPlots ['highPurityPt10'] = cms.bool(True)
@@ -174,6 +184,8 @@
174184
trackPtN ['highPurityPt1'] = cms.int32(100)
175185
trackPtMin ['highPurityPt1'] = cms.double(0.)
176186
trackPtMax ['highPurityPt1'] = cms.double(100.)
187+
pcaZtoPVMax ['highPurityPt1'] = cms.double(30.)
188+
pcaZtoPVMin ['highPurityPt1'] = cms.double(-30.)
177189
doPlotsPCA ['highPurityPt1'] = cms.bool(True)
178190
numCutString ['highPurityPt1'] = cms.string(" pt >= 1 & quality('highPurity')") # num := den + quality('highPurity') [it is the same as the main selection, but just to be sure ...]
179191
denCutString ['highPurityPt1'] = cms.string(" pt >= 1 ") # den := kinematics cut
@@ -208,6 +220,8 @@
208220
trackPtN ['highPurityPt1Eta2p5to3p0'] = cms.int32(100)
209221
trackPtMin ['highPurityPt1Eta2p5to3p0'] = cms.double(0.)
210222
trackPtMax ['highPurityPt1Eta2p5to3p0'] = cms.double(100.)
223+
pcaZtoPVMax ['highPurityPt1Eta2p5to3p0'] = cms.double(30.)
224+
pcaZtoPVMin ['highPurityPt1Eta2p5to3p0'] = cms.double(-30.)
211225
doPlotsPCA ['highPurityPt1Eta2p5to3p0'] = cms.bool(True)
212226
numCutString ['highPurityPt1Eta2p5to3p0'] = cms.string(" pt >= 1 & abs(eta) > 2.5 & quality('highPurity')") # num := den + quality('highPurity') [it is the same as the main selection, but just to be sure ...]
213227
denCutString ['highPurityPt1Eta2p5to3p0'] = cms.string(" pt >= 1 & abs(eta) > 2.5") # den := kinematics cut
@@ -289,6 +303,8 @@
289303
trackPtN ['highPurityPV0p1'] = cms.int32(100)
290304
trackPtMin ['highPurityPV0p1'] = cms.double(0.)
291305
trackPtMax ['highPurityPV0p1'] = cms.double(100.)
306+
pcaZtoPVMax ['highPurityPV0p1'] = cms.double(0.2)
307+
pcaZtoPVMin ['highPurityPV0p1'] = cms.double(-0.2)
292308
doPlotsPCA ['highPurityPV0p1'] = cms.bool(True)
293309
numCutString ['highPurityPV0p1'] = cms.string("quality('highPurity')") # num := den + quality('highPurity')
294310
denCutString ['highPurityPV0p1'] = cms.string("") # den := kinematic cuts
@@ -324,6 +340,8 @@
324340
trackPtN ['hiConformalPixelTracksQA'] = cms.int32(100)
325341
trackPtMin ['hiConformalPixelTracksQA'] = cms.double(0.)
326342
trackPtMax ['hiConformalPixelTracksQA'] = cms.double(10.)
343+
pcaZtoPVMax ['hiConformalPixelTracksQA'] = cms.double(30.)
344+
pcaZtoPVMin ['hiConformalPixelTracksQA'] = cms.double(-30.)
327345
numCutString ['hiConformalPixelTracksQA'] = cms.string(" pt >= 0 ")
328346
denCutString ['hiConformalPixelTracksQA'] = cms.string(" pt >= 0 ")
329347
doPlotsPCA ['hiConformalPixelTracksQA'] = cms.bool(False)

‎DQM/TrackingMonitorSource/python/TrackingSourceConfig_Tier0_cff.py‎

Lines changed: 10 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -57,6 +57,8 @@
5757
locals()[label].TrackPBin = trackPtN[tracks]
5858
locals()[label].TrackPMin = trackPtMin[tracks]
5959
locals()[label].TrackPMax = trackPtMax[tracks]
60+
locals()[label].VZ_PVMax = pcaZtoPVMax[tracks]
61+
locals()[label].VZ_PVMin = pcaZtoPVMin[tracks]
6062
locals()[label].doDCAPlots = doPlotsPCA[tracks]
6163
locals()[label].doDCAwrtPVPlots = doPlotsPCA[tracks]
6264
locals()[label].doDCAwrt000Plots = doPlotsPCA[tracks]
@@ -98,6 +100,8 @@
98100
locals()[label].TrackPBin = trackPtN[tracks]
99101
locals()[label].TrackPMin = trackPtMin[tracks]
100102
locals()[label].TrackPMax = trackPtMax[tracks]
103+
locals()[label].VZ_PVMax = pcaZtoPVMax[tracks]
104+
locals()[label].VZ_PVMin = pcaZtoPVMin[tracks]
101105
locals()[label].doDCAPlots = doPlotsPCA[tracks]
102106
locals()[label].doDCAwrtPVPlots = doPlotsPCA[tracks]
103107
locals()[label].doDCAwrt000Plots = doPlotsPCA[tracks]
@@ -138,6 +142,8 @@
138142
locals()[label].TrackPBin = trackPtN[tracks]
139143
locals()[label].TrackPMin = trackPtMin[tracks]
140144
locals()[label].TrackPMax = trackPtMax[tracks]
145+
locals()[label].VZ_PVMax = pcaZtoPVMax[tracks]
146+
locals()[label].VZ_PVMin = pcaZtoPVMin[tracks]
141147
locals()[label].doDCAPlots = doPlotsPCA[tracks]
142148
locals()[label].doDCAwrtPVPlots = doPlotsPCA[tracks]
143149
locals()[label].doDCAwrt000Plots = doPlotsPCA[tracks]
@@ -178,6 +184,8 @@
178184
locals()[label].TrackPBin = trackPtN[tracks]
179185
locals()[label].TrackPMin = trackPtMin[tracks]
180186
locals()[label].TrackPMax = trackPtMax[tracks]
187+
locals()[label].VZ_PVMax = pcaZtoPVMax[tracks]
188+
locals()[label].VZ_PVMin = pcaZtoPVMin[tracks]
181189
locals()[label].doDCAPlots = doPlotsPCA[tracks]
182190
locals()[label].doDCAwrtPVPlots = doPlotsPCA[tracks]
183191
locals()[label].doDCAwrt000Plots = doPlotsPCA[tracks]
@@ -219,6 +227,8 @@
219227
locals()[label].TrackPBin = trackPtN[tracks]
220228
locals()[label].TrackPMin = trackPtMin[tracks]
221229
locals()[label].TrackPMax = trackPtMax[tracks]
230+
locals()[label].VZ_PVMax = pcaZtoPVMax[tracks]
231+
locals()[label].VZ_PVMin = pcaZtoPVMin[tracks]
222232
locals()[label].doDCAPlots = doPlotsPCA[tracks]
223233
locals()[label].doDCAwrtPVPlots = doPlotsPCA[tracks]
224234
locals()[label].doDCAwrt000Plots = doPlotsPCA[tracks]
@@ -450,4 +460,3 @@ def _copyIfExists(mod, pset, name):
450460
TrackingDQMSourceTier0MinBias += primaryVertexResolution
451461

452462
TrackingDQMSourceTier0MinBias += dqmInfoTracking
453-

0 commit comments

Comments
 (0)