Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

update track candidate DQM code #10126

Merged
merged 7 commits into from Jul 28, 2015
Merged

update track candidate DQM code #10126

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
19a8068
add track candidate phi vs eta plot
mtosi Jul 10, 2015
7a3cdf9
delete debugging line
mtosi Jul 10, 2015
2355ae7
add jetCoreRegionalStep and muonSeededStepOutInDisplaced
mtosi Jul 10, 2015
6ebf5d8
exclude muonSeededStepOutInDisplaced
mtosi Jul 10, 2015
56b8aa3
restore file in release
mtosi Jul 10, 2015
d9f41e4
restore file in release
mtosi Jul 10, 2015
6478127
fix bug
mtosi Jul 13, 2015
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
4 changes: 3 additions & 1 deletion DQM/TrackingMonitor/interface/TrackBuildingAnalyzer.h
View file
Open in desktop
Expand Up @@ -82,7 +82,8 @@ class TrackBuildingAnalyzer
MonitorElement* TrackCandPt;
MonitorElement* TrackCandEta;
MonitorElement* TrackCandPhi;
MonitorElement* TrackCandTheta;
MonitorElement* TrackCandPhiVsEta;
MonitorElement* TrackCandTheta;
MonitorElement* TrackCandQ;
MonitorElement* TrackCandDxy;
MonitorElement* TrackCandDz;
Expand All @@ -96,6 +97,7 @@ class TrackBuildingAnalyzer
bool doAllPlots;
bool doAllSeedPlots;
bool doTCPlots;
bool doAllTCPlots;
bool doPT;
bool doETA;
bool doPHI;
Expand Down
1 change: 1 addition & 0 deletions DQM/TrackingMonitor/python/TrackingMonitor_cfi.py
View file
Open in desktop
Expand Up @@ -54,6 +54,7 @@
doPrimaryVertexPlots = cms.bool(False),
doSeedParameterHistos = cms.bool(False),
doTrackCandHistos = cms.bool(False),
doAllTrackCandHistos = cms.bool(False),
doDCAPlots = cms.bool(False),
doDCAwrtPVPlots = cms.bool(False),
doDCAwrt000Plots = cms.bool(False),
Expand Down
89 changes: 50 additions & 39 deletions DQM/TrackingMonitor/src/TrackBuildingAnalyzer.cc
View file
Open in desktop
Expand Up @@ -101,6 +101,7 @@ void TrackBuildingAnalyzer::initHisto(DQMStore::IBooker & ibooker)
doAllPlots = conf_.getParameter<bool>("doAllPlots");
doAllSeedPlots = conf_.getParameter<bool>("doSeedParameterHistos");
doTCPlots = conf_.getParameter<bool>("doTrackCandHistos");
doAllTCPlots = conf_.getParameter<bool>("doAllTrackCandHistos");
doPT = conf_.getParameter<bool>("doSeedPTHisto");
doETA = conf_.getParameter<bool>("doSeedETAHisto");
doPHI = conf_.getParameter<bool>("doSeedPHIHisto");
Expand Down Expand Up @@ -221,40 +222,49 @@ void TrackBuildingAnalyzer::initHisto(DQMStore::IBooker & ibooker)
TrackCandPhi->setAxisTitle("Track Candidate #phi", 1);
TrackCandPhi->setAxisTitle("Number of Track Candidates", 2);

histname = "TrackCandTheta_"+tcProducer.label() + "_";
TrackCandTheta = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, ThetaBin, ThetaMin, ThetaMax);
TrackCandTheta->setAxisTitle("Track Candidate #theta", 1);
TrackCandTheta->setAxisTitle("Number of Track Candidates", 2);
if (doTheta) {
histname = "TrackCandTheta_"+tcProducer.label() + "_";
TrackCandTheta = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, ThetaBin, ThetaMin, ThetaMax);
TrackCandTheta->setAxisTitle("Track Candidate #theta", 1);
TrackCandTheta->setAxisTitle("Number of Track Candidates", 2);
}

histname = "TrackCandQ_"+tcProducer.label() + "_";
TrackCandQ = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TrackQBin, TrackQMin, TrackQMax);
TrackCandQ->setAxisTitle("Track Candidate Charge", 1);
TrackCandQ->setAxisTitle("Number of Track Candidates",2);

histname = "TrackCandDxy_"+tcProducer.label() + "_";
TrackCandDxy = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TCDxyBin, TCDxyMin, TCDxyMax);
TrackCandDxy->setAxisTitle("Track Candidate d_{xy} (cm)", 1);
TrackCandDxy->setAxisTitle("Number of Track Candidates",2);

histname = "TrackCandDz_"+tcProducer.label() + "_";
TrackCandDz = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TCDzBin, TCDzMin, TCDzMax);
TrackCandDz->setAxisTitle("Track Candidate d_{z} (cm)", 1);
TrackCandDz->setAxisTitle("Number of Track Candidates",2);

histname = "NumberOfRecHitsPerTrackCand_"+tcProducer.label() + "_";
NumberOfRecHitsPerTrackCand = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TCHitBin, TCHitMin, TCHitMax);
NumberOfRecHitsPerTrackCand->setAxisTitle("Number of RecHits per Track Candidate", 1);
NumberOfRecHitsPerTrackCand->setAxisTitle("Number of Track Candidates",2);

histname = "NumberOfRecHitsPerTrackCandVsPhiProfile_"+tcProducer.label() + "_";
NumberOfRecHitsPerTrackCandVsPhiProfile = ibooker.bookProfile(histname+CatagoryName, histname+CatagoryName, PhiBin, PhiMin, PhiMax, TCHitBin, TCHitMin, TCHitMax,"s");
NumberOfRecHitsPerTrackCandVsPhiProfile->setAxisTitle("Track Candidate #phi",1);
NumberOfRecHitsPerTrackCandVsPhiProfile->setAxisTitle("Number of RecHits of each Track Candidate",2);
if (doAllTCPlots) {
histname = "TrackCandQ_"+tcProducer.label() + "_";
TrackCandQ = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TrackQBin, TrackQMin, TrackQMax);
TrackCandQ->setAxisTitle("Track Candidate Charge", 1);
TrackCandQ->setAxisTitle("Number of Track Candidates",2);

histname = "TrackCandDxy_"+tcProducer.label() + "_";
TrackCandDxy = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TCDxyBin, TCDxyMin, TCDxyMax);
TrackCandDxy->setAxisTitle("Track Candidate d_{xy} (cm)", 1);
TrackCandDxy->setAxisTitle("Number of Track Candidates",2);

histname = "TrackCandDz_"+tcProducer.label() + "_";
TrackCandDz = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TCDzBin, TCDzMin, TCDzMax);
TrackCandDz->setAxisTitle("Track Candidate d_{z} (cm)", 1);
TrackCandDz->setAxisTitle("Number of Track Candidates",2);

histname = "NumberOfRecHitsPerTrackCand_"+tcProducer.label() + "_";
NumberOfRecHitsPerTrackCand = ibooker.book1D(histname+CatagoryName, histname+CatagoryName, TCHitBin, TCHitMin, TCHitMax);
NumberOfRecHitsPerTrackCand->setAxisTitle("Number of RecHits per Track Candidate", 1);
NumberOfRecHitsPerTrackCand->setAxisTitle("Number of Track Candidates",2);

histname = "NumberOfRecHitsPerTrackCandVsEtaProfile_"+tcProducer.label() + "_";
NumberOfRecHitsPerTrackCandVsEtaProfile = ibooker.bookProfile(histname+CatagoryName, histname+CatagoryName, EtaBin, EtaMin, EtaMax, TCHitBin, TCHitMin, TCHitMax,"s");
NumberOfRecHitsPerTrackCandVsEtaProfile->setAxisTitle("Track Candidate #eta",1);
NumberOfRecHitsPerTrackCandVsEtaProfile->setAxisTitle("Number of RecHits of each Track Candidate",2);
histname = "NumberOfRecHitsPerTrackCandVsPhiProfile_"+tcProducer.label() + "_";
NumberOfRecHitsPerTrackCandVsPhiProfile = ibooker.bookProfile(histname+CatagoryName, histname+CatagoryName, PhiBin, PhiMin, PhiMax, TCHitBin, TCHitMin, TCHitMax,"s");
NumberOfRecHitsPerTrackCandVsPhiProfile->setAxisTitle("Track Candidate #phi",1);
NumberOfRecHitsPerTrackCandVsPhiProfile->setAxisTitle("Number of RecHits of each Track Candidate",2);

histname = "NumberOfRecHitsPerTrackCandVsEtaProfile_"+tcProducer.label() + "_";
NumberOfRecHitsPerTrackCandVsEtaProfile = ibooker.bookProfile(histname+CatagoryName, histname+CatagoryName, EtaBin, EtaMin, EtaMax, TCHitBin, TCHitMin, TCHitMax,"s");
NumberOfRecHitsPerTrackCandVsEtaProfile->setAxisTitle("Track Candidate #eta",1);
NumberOfRecHitsPerTrackCandVsEtaProfile->setAxisTitle("Number of RecHits of each Track Candidate",2);
}

histname = "TrackCandPhiVsEta_"+tcProducer.label() + "_";
TrackCandPhiVsEta = ibooker.book2D(histname+CatagoryName, histname+CatagoryName, EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax);
TrackCandPhiVsEta->setAxisTitle("Track Candidate #eta", 1);
TrackCandPhiVsEta->setAxisTitle("Track Candidate #phi", 2);
}

}
Expand Down Expand Up @@ -355,15 +365,16 @@ void TrackBuildingAnalyzer::analyze

if (doTCPlots){
// fill the ME's
TrackCandQ->Fill( state.charge() );
if (doAllTCPlots) TrackCandQ->Fill( state.charge() );
TrackCandPt->Fill( pt );
TrackCandEta->Fill( eta );
TrackCandPhi->Fill( phi );
TrackCandTheta->Fill( theta );
TrackCandDxy->Fill( dxy );
TrackCandDz->Fill( dz );
NumberOfRecHitsPerTrackCand->Fill( numberOfHits );
NumberOfRecHitsPerTrackCandVsEtaProfile->Fill( eta, numberOfHits );
NumberOfRecHitsPerTrackCandVsPhiProfile->Fill( phi, numberOfHits );
TrackCandPhiVsEta->Fill( eta, phi );
if (doTheta) TrackCandTheta->Fill( theta );
if (doAllTCPlots) TrackCandDxy->Fill( dxy );
if (doAllTCPlots) TrackCandDz->Fill( dz );
if (doAllTCPlots) NumberOfRecHitsPerTrackCand->Fill( numberOfHits );
if (doAllTCPlots) NumberOfRecHitsPerTrackCandVsEtaProfile->Fill( eta, numberOfHits );
if (doAllTCPlots) NumberOfRecHitsPerTrackCandVsPhiProfile->Fill( phi, numberOfHits );
}
}
22 changes: 21 additions & 1 deletion DQM/TrackingMonitorSource/python/IterTrackingModules4seedMonitoring_cfi.py
View file
Open in desktop
Expand Up @@ -12,7 +12,7 @@
clusterMax = {}

seedInputTag ['initialStep'] = cms.InputTag("initialStepSeeds")
trackCandInputTag['initialStep'] = cms.InputTag("initialStepSeeds")
trackCandInputTag['initialStep'] = cms.InputTag("initialStepTrackCandidates")
trackSeedSizeBin ['initialStep'] = cms.int32(100) # could be 50 ?
trackSeedSizeMin ['initialStep'] = cms.double(0)
trackSeedSizeMax ['initialStep'] = cms.double(5000)
Expand Down Expand Up @@ -92,12 +92,32 @@
clusterBin ['muonSeededStepOutIn'] = cms.int32(500)
clusterMax ['muonSeededStepOutIn'] = cms.double(100000)

seedInputTag ['muonSeededStepOutInDisplaced'] = cms.InputTag("muonSeededSeedsOutInDisplaced")
trackCandInputTag['muonSeededStepOutInDisplaced'] = cms.InputTag("muonSeededTrackCandidatesOutInDisplacedg")
trackSeedSizeBin ['muonSeededStepOutInDisplaced'] = cms.int32(15)
trackSeedSizeMin ['muonSeededStepOutInDisplaced'] = cms.double(-0.5)
trackSeedSizeMax ['muonSeededStepOutInDisplaced'] = cms.double(14.5)
clusterLabel ['muonSeededStepOutInDisplaced'] = cms.vstring('Strip')
clusterBin ['muonSeededStepOutInDisplaced'] = cms.int32(500)
clusterMax ['muonSeededStepOutInDisplaced'] = cms.double(100000)

seedInputTag ['jetCoreRegionalStep'] = cms.InputTag("jetCoreRegionalStepSeeds")
trackCandInputTag['jetCoreRegionalStep'] = cms.InputTag("jetCoreRegionalStepTrackCandidates")
trackSeedSizeBin ['jetCoreRegionalStep'] = cms.int32(15)
trackSeedSizeMin ['jetCoreRegionalStep'] = cms.double(-0.5)
trackSeedSizeMax ['jetCoreRegionalStep'] = cms.double(14.5)
clusterLabel ['jetCoreRegionalStep'] = cms.vstring('Tot')
clusterBin ['jetCoreRegionalStep'] = cms.int32(500)
clusterMax ['jetCoreRegionalStep'] = cms.double(100000)

selectedIterTrackingStep.extend( ['initialStep'] )
selectedIterTrackingStep.extend( ['lowPtTripletStep'] )
selectedIterTrackingStep.extend( ['pixelPairStep'] )
selectedIterTrackingStep.extend( ['detachedTripletStep'] )
selectedIterTrackingStep.extend( ['mixedTripletStep'] )
selectedIterTrackingStep.extend( ['pixelLessStep'] )
selectedIterTrackingStep.extend( ['tobTecStep'] )
selectedIterTrackingStep.extend( ['jetCoreRegionalStep'] )
selectedIterTrackingStep.extend( ['muonSeededStepInOut'] )
selectedIterTrackingStep.extend( ['muonSeededStepOutIn'] )
#selectedIterTrackingStep.extend( ['muonSeededStepOutInDisplaced'] )
4 changes: 3 additions & 1 deletion DQM/TrackingMonitorSource/python/TrackingSourceConfig_Tier0_cff.py
View file
Open in desktop
Expand Up @@ -55,7 +55,9 @@
from DQM.TrackingMonitorSource.IterTrackingModules4seedMonitoring_cfi import *
for step in selectedIterTrackingStep :
label = 'TrackSeedMon'+str(step)
locals()[label] = DQM.TrackingMonitor.TrackingMonitorSeed_cfi.TrackMonSeed.clone()
locals()[label] = DQM.TrackingMonitor.TrackingMonitorSeed_cfi.TrackMonSeed.clone(
doTrackCandHistos = cms.bool(True)
)
locals()[label].TrackProducer = cms.InputTag("generalTracks")
locals()[label].FolderName = cms.string("Tracking/TrackParameters/generalTracks")
locals()[label].SeedProducer = seedInputTag[step]
Expand Down