Skip to content

Commit

Permalink
add SiStripHitEfficiency Summary plots In DQM file
Browse files Browse the repository at this point in the history
  • Loading branch information
@mmusich
mmusich committed Dec 2, 2022
1 parent 82bcefc commit e0879e0
Showing 1 changed file with 156 additions and 119 deletions.
275 changes: 156 additions & 119 deletions CalibTracker/SiStripHitEfficiency/plugins/SiStripHitEfficiencyHarvester.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +88,7 @@ class SiStripHitEfficiencyHarvester : public DQMEDHarvester {
void printAndWriteBadModules(const SiStripQuality& quality, const SiStripDetInfo& detInfo) const;
bool checkMapsValidity(const std::vector<MonitorElement*>& maps, const std::string& type) const;
unsigned int countTotalHits(const std::vector<MonitorElement*>& maps); /* to check if TK was ON */
void makeSummary(DQMStore::IGetter& getter, TFileService& fs) const;
void makeSummary(DQMStore::IGetter& getter, DQMStore::IBooker& booker) const;
void makeSummaryVsBX(DQMStore::IGetter& getter, TFileService& fs) const;
void makeSummaryVsLumi(DQMStore::IGetter& getter) const;
void makeSummaryVsCM(DQMStore::IGetter& getter, TFileService& fs) const;
Expand Down Expand Up @@ -182,6 +182,10 @@ unsigned int SiStripHitEfficiencyHarvester::countTotalHits(const std::vector<Mon
void SiStripHitEfficiencyHarvester::dqmEndJob(DQMStore::IBooker& booker, DQMStore::IGetter& getter) {
if (!isAtPCL_) {
edm::Service<TFileService> fs;
if (!fs.isAvailable()) {
throw cms::Exception("BadConfig") << "TFileService unavailable: "
<< "please add it to config file";
}

if (doStoreOnTree_) {
// store information per DetId in the output tree
Expand Down Expand Up @@ -480,14 +484,16 @@ void SiStripHitEfficiencyHarvester::dqmEndJob(DQMStore::IBooker& booker, DQMStor
// << " lumiBlock " << e.luminosityBlock() << " time " << e.time().value() << "\n-----------------\n";
printAndWriteBadModules(pQuality, detInfo); // TODO

if (!isAtPCL_) {
edm::Service<TFileService> fs;
makeSummary(getter, *fs); // TODO
// make summary plots
makeSummary(getter, booker);
makeSummaryVsLumi(getter); // TODO

/*
if (!isAtPCL_) {
makeSummaryVsBX(getter, *fs); // TODO
makeSummaryVsCM(getter, *fs); // TODO
}

makeSummaryVsLumi(getter); // TODO
}
*/
}

void SiStripHitEfficiencyHarvester::printTotalStatistics(
Expand Down Expand Up @@ -556,7 +562,7 @@ void SiStripHitEfficiencyHarvester::writeBadStripPayload(const SiStripQuality& q
}
}

void SiStripHitEfficiencyHarvester::makeSummary(DQMStore::IGetter& getter, TFileService& fs) const {
void SiStripHitEfficiencyHarvester::makeSummary(DQMStore::IGetter& getter, DQMStore::IBooker& booker) const {
// use goodlayer_total/found and alllayer_total/found, collapse side and/or ring if needed

unsigned int nLayers = 34;
Expand All @@ -569,20 +575,23 @@ void SiStripHitEfficiencyHarvester::makeSummary(DQMStore::IGetter& getter, TFile
nLayers = 20;
}

TH1F* found = fs.make<TH1F>("found", "found", nLayers + 1, 0, nLayers + 1);
TH1F* all = fs.make<TH1F>("all", "all", nLayers + 1, 0, nLayers + 1);
TH1F* found2 = fs.make<TH1F>("found2", "found2", nLayers + 1, 0, nLayers + 1);
TH1F* all2 = fs.make<TH1F>("all2", "all2", nLayers + 1, 0, nLayers + 1);
// come back to the main folder and create a final efficiency folder
booker.setCurrentFolder(fmt::format("{}/FinalEfficiency", inputFolder_));
MonitorElement* found = booker.book1D("found", "found", nLayers + 1, 0, nLayers + 1);
MonitorElement* all = booker.book1D("all", "all", nLayers + 1, 0, nLayers + 1);
MonitorElement* found2 = booker.book1D("found2", "found", nLayers + 1, 0, nLayers + 1);
MonitorElement* all2 = booker.book1D("all2", "all2", nLayers + 1, 0, nLayers + 1);

// first bin only to keep real data off the y axis so set to -1
found->SetBinContent(0, -1);
all->SetBinContent(0, 1);
found->setBinContent(0, -1);
all->setBinContent(0, 1);

// new ROOT version: TGraph::Divide don't handle null or negative values
for (unsigned int i = 1; i < nLayers + 2; ++i) {
found->SetBinContent(i, 1e-6);
all->SetBinContent(i, 1);
found2->SetBinContent(i, 1e-6);
all2->SetBinContent(i, 1);
found->setBinContent(i, 1e-6);
all->setBinContent(i, 1);
found2->setBinContent(i, 1e-6);
all2->setBinContent(i, 1);
}

TCanvas* c7 = new TCanvas("c7", " test ", 10, 10, 800, 600);
Expand All @@ -596,14 +605,14 @@ void SiStripHitEfficiencyHarvester::makeSummary(DQMStore::IGetter& getter, TFile
LOGPRINT << "Fill only good modules layer " << i << ": S = " << goodlayerfound[i]
<< " B = " << goodlayertotal[i];
if (goodlayertotal[i] > 5) {
found->SetBinContent(i, goodlayerfound[i]);
all->SetBinContent(i, goodlayertotal[i]);
found->setBinContent(i, goodlayerfound[i]);
all->setBinContent(i, goodlayertotal[i]);
}

LOGPRINT << "Filling all modules layer " << i << ": S = " << alllayerfound[i] << " B = " << alllayertotal[i];
if (alllayertotal[i] > 5) {
found2->SetBinContent(i, alllayerfound[i]);
all2->SetBinContent(i, alllayertotal[i]);
found2->setBinContent(i, alllayerfound[i]);
all2->setBinContent(i, alllayertotal[i]);
}
}

Expand All @@ -613,139 +622,169 @@ void SiStripHitEfficiencyHarvester::makeSummary(DQMStore::IGetter& getter, TFile
LOGPRINT << "Fill only good modules layer " << i << ": S = " << goodlayerfound[i] + goodlayerfound[i + 3]
<< " B = " << goodlayertotal[i] + goodlayertotal[i + 3];
if (goodlayertotal[i] + goodlayertotal[i + 3] > 5) {
found->SetBinContent(i, goodlayerfound[i] + goodlayerfound[i + 3]);
all->SetBinContent(i, goodlayertotal[i] + goodlayertotal[i + 3]);
found->setBinContent(i, goodlayerfound[i] + goodlayerfound[i + 3]);
all->setBinContent(i, goodlayertotal[i] + goodlayertotal[i + 3]);
}
LOGPRINT << "Filling all modules layer " << i << ": S = " << alllayerfound[i] + alllayerfound[i + 3]
<< " B = " << alllayertotal[i] + alllayertotal[i + 3];
if (alllayertotal[i] + alllayertotal[i + 3] > 5) {
found2->SetBinContent(i, alllayerfound[i] + alllayerfound[i + 3]);
all2->SetBinContent(i, alllayertotal[i] + alllayertotal[i + 3]);
found2->setBinContent(i, alllayerfound[i] + alllayerfound[i + 3]);
all2->setBinContent(i, alllayertotal[i] + alllayertotal[i + 3]);
}
}
for (unsigned int i = 17; i < 17 + nTEClayers_; ++i) { // TEC disks
LOGPRINT << "Fill only good modules layer " << i - 3
<< ": S = " << goodlayerfound[i] + goodlayerfound[i + nTEClayers_]
<< " B = " << goodlayertotal[i] + goodlayertotal[i + nTEClayers_];
if (goodlayertotal[i] + goodlayertotal[i + nTEClayers_] > 5) {
found->SetBinContent(i - 3, goodlayerfound[i] + goodlayerfound[i + nTEClayers_]);
all->SetBinContent(i - 3, goodlayertotal[i] + goodlayertotal[i + nTEClayers_]);
found->setBinContent(i - 3, goodlayerfound[i] + goodlayerfound[i + nTEClayers_]);
all->setBinContent(i - 3, goodlayertotal[i] + goodlayertotal[i + nTEClayers_]);
}
LOGPRINT << "Filling all modules layer " << i - 3
<< ": S = " << alllayerfound[i] + alllayerfound[i + nTEClayers_]
<< " B = " << alllayertotal[i] + alllayertotal[i + nTEClayers_];
if (alllayertotal[i] + alllayertotal[i + nTEClayers_] > 5) {
found2->SetBinContent(i - 3, alllayerfound[i] + alllayerfound[i + nTEClayers_]);
all2->SetBinContent(i - 3, alllayertotal[i] + alllayertotal[i + nTEClayers_]);
found2->setBinContent(i - 3, alllayerfound[i] + alllayerfound[i + nTEClayers_]);
all2->setBinContent(i - 3, alllayertotal[i] + alllayertotal[i + nTEClayers_]);
}
}
}

found->Sumw2();
all->Sumw2();
found->getTH1F()->Sumw2();
all->getTH1F()->Sumw2();

found2->Sumw2();
all2->Sumw2();
found2->getTH1F()->Sumw2();
all2->getTH1F()->Sumw2();

TGraphAsymmErrors* gr = fs.make<TGraphAsymmErrors>(nLayers + 1);
gr->SetName("eff_good");
gr->BayesDivide(found, all);
MonitorElement* eff_all = booker.book1D("eff_all", "efficiency for all modules", nLayers + 1, 0, nLayers + 1);
MonitorElement* eff_good = booker.book1D("eff_good", "efficiency for good modules", nLayers + 1, 0, nLayers + 1);

TGraphAsymmErrors* gr2 = fs.make<TGraphAsymmErrors>(nLayers + 1);
gr2->SetName("eff_all");
gr2->BayesDivide(found2, all2);
for (int i = 1; i < found->getNbinsX(); i++) {
const auto& den_all = all2->getBinContent(i);
const auto& num_all = found2->getBinContent(i);
const auto& den_good = all->getBinContent(i);
const auto& num_good = found->getBinContent(i);

for (unsigned int j = 0; j < nLayers + 1; j++) {
gr->SetPointError(j, 0., 0., gr->GetErrorYlow(j), gr->GetErrorYhigh(j));
gr2->SetPointError(j, 0., 0., gr2->GetErrorYlow(j), gr2->GetErrorYhigh(j));
if (den_all > 0.) {
eff_all->setBinContent(i, num_all / den_all);
}
if (den_good > 0.) {
eff_good->setBinContent(i, num_good / den_good);
}
}

gr->GetXaxis()->SetLimits(0, nLayers);
gr->SetMarkerColor(2);
gr->SetMarkerSize(1.2);
gr->SetLineColor(2);
gr->SetLineWidth(4);
gr->SetMarkerStyle(20);
gr->SetMinimum(effPlotMin_);
gr->SetMaximum(1.001);
gr->GetYaxis()->SetTitle("Efficiency");
gStyle->SetTitleFillColor(0);
gStyle->SetTitleBorderSize(0);
gr->SetTitle(title_.c_str());

gr2->GetXaxis()->SetLimits(0, nLayers);
gr2->SetMarkerColor(1);
gr2->SetMarkerSize(1.2);
gr2->SetLineColor(1);
gr2->SetLineWidth(4);
gr2->SetMarkerStyle(21);
gr2->SetMinimum(effPlotMin_);
gr2->SetMaximum(1.001);
gr2->GetYaxis()->SetTitle("Efficiency");
gr2->SetTitle(title_.c_str());

for (Long_t k = 1; k < nLayers + 1; k++) {
TString label;
if (showEndcapSides_)
label = ::layerSideName(k, showRings_, nTEClayers_);
else
label = ::layerName(k, showRings_, nTEClayers_);
if (!showTOB6TEC9_) {
if (k == 10)
label = "";
if (!showRings_ && k == nLayers)
label = "";
if (!showRings_ && showEndcapSides_ && k == 25)
label = "";
if (!isAtPCL_) {
// if TFileService is not avaible, just go on
edm::Service<TFileService> fs;
if (!fs.isAvailable()) {
throw cms::Exception("BadConfig") << "TFileService unavailable: "
<< "please add it to config file";
}
if (!showRings_) {
if (showEndcapSides_) {
gr->GetXaxis()->SetBinLabel(((k + 1) * 100 + 2) / (nLayers)-4, label);
gr2->GetXaxis()->SetBinLabel(((k + 1) * 100 + 2) / (nLayers)-4, label);
} else {
gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-6, label);
gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-6, label);

TGraphAsymmErrors* gr = (*fs).make<TGraphAsymmErrors>(nLayers + 1);
gr->SetName("eff_good");
gr->BayesDivide(found->getTH1F(), all->getTH1F());

TGraphAsymmErrors* gr2 = (*fs).make<TGraphAsymmErrors>(nLayers + 1);
gr2->SetName("eff_all");
gr2->BayesDivide(found2->getTH1F(), all2->getTH1F());

for (unsigned int j = 0; j < nLayers + 1; j++) {
gr->SetPointError(j, 0., 0., gr->GetErrorYlow(j), gr->GetErrorYhigh(j));
gr2->SetPointError(j, 0., 0., gr2->GetErrorYlow(j), gr2->GetErrorYhigh(j));
}

gr->GetXaxis()->SetLimits(0, nLayers);
gr->SetMarkerColor(2);
gr->SetMarkerSize(1.2);
gr->SetLineColor(2);
gr->SetLineWidth(4);
gr->SetMarkerStyle(20);
gr->SetMinimum(effPlotMin_);
gr->SetMaximum(1.001);
gr->GetYaxis()->SetTitle("Efficiency");
gStyle->SetTitleFillColor(0);
gStyle->SetTitleBorderSize(0);
gr->SetTitle("SiStripHitEfficiency by Layer");

gr2->GetXaxis()->SetLimits(0, nLayers);
gr2->SetMarkerColor(1);
gr2->SetMarkerSize(1.2);
gr2->SetLineColor(1);
gr2->SetLineWidth(4);
gr2->SetMarkerStyle(21);
gr2->SetMinimum(effPlotMin_);
gr2->SetMaximum(1.001);
gr2->GetYaxis()->SetTitle("Efficiency");
gr2->SetTitle("SiStripHitEfficiency by Layer");

for (Long_t k = 1; k < nLayers + 1; k++) {
TString label;
if (showEndcapSides_)
label = ::layerSideName(k, showRings_, nTEClayers_);
else
label = ::layerName(k, showRings_, nTEClayers_);
if (!showTOB6TEC9_) {
if (k == 10)
label = "";
if (!showRings_ && k == nLayers)
label = "";
if (!showRings_ && showEndcapSides_ && k == 25)
label = "";
}
} else {
if (showEndcapSides_) {
gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-4, label);
gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-4, label);
if (!showRings_) {
if (showEndcapSides_) {
gr->GetXaxis()->SetBinLabel(((k + 1) * 100 + 2) / (nLayers)-4, label);
gr2->GetXaxis()->SetBinLabel(((k + 1) * 100 + 2) / (nLayers)-4, label);
} else {
gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-6, label);
gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-6, label);
}
} else {
gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-7, label);
gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-7, label);
if (showEndcapSides_) {
gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-4, label);
gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-4, label);
} else {
gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-7, label);
gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-7, label);
}
}
}
}

gr->Draw("AP");
gr->GetXaxis()->SetNdivisions(36);

c7->cd();
TPad* overlay = new TPad("overlay", "", 0, 0, 1, 1);
overlay->SetFillStyle(4000);
overlay->SetFillColor(0);
overlay->SetFrameFillStyle(4000);
overlay->Draw("same");
overlay->cd();
if (!showOnlyGoodModules_)
gr2->Draw("AP");

TLegend* leg = new TLegend(0.70, 0.27, 0.88, 0.40);
leg->AddEntry(gr, "Good Modules", "p");
if (!showOnlyGoodModules_)
leg->AddEntry(gr2, "All Modules", "p");
leg->SetTextSize(0.020);
leg->SetFillColor(0);
leg->Draw("same");

c7->SaveAs("Summary.png");
c7->SaveAs("Summary.root");
gr->Draw("AP");
gr->GetXaxis()->SetNdivisions(36);

c7->cd();
TPad* overlay = new TPad("overlay", "", 0, 0, 1, 1);
overlay->SetFillStyle(4000);
overlay->SetFillColor(0);
overlay->SetFrameFillStyle(4000);
overlay->Draw("same");
overlay->cd();
if (!showOnlyGoodModules_)
gr2->Draw("AP");

TLegend* leg = new TLegend(0.70, 0.27, 0.88, 0.40);
leg->AddEntry(gr, "Good Modules", "p");
if (!showOnlyGoodModules_)
leg->AddEntry(gr2, "All Modules", "p");
leg->SetTextSize(0.020);
leg->SetFillColor(0);
leg->Draw("same");

c7->SaveAs("Summary.png");
c7->SaveAs("Summary.root");
} // if it's not run at PCL
}

// not yet implemented
void SiStripHitEfficiencyHarvester::makeSummaryVsBX(DQMStore::IGetter& getter, TFileService& fs) const {
// use found/totalVsBx_layer%i [0,23)
}

// not yet impelement
void SiStripHitEfficiencyHarvester::makeSummaryVsCM(DQMStore::IGetter& getter, TFileService& fs) const {}

void SiStripHitEfficiencyHarvester::makeSummaryVsLumi(DQMStore::IGetter& getter) const {
for (unsigned int iLayer = 1; iLayer != (showRings_ ? 20 : 22); ++iLayer) {
auto hfound = getter.get(fmt::format("{}/VsLumi/layerfound_vsLumi_layer_{}", inputFolder_, iLayer))->getTH1F();
Expand Down Expand Up @@ -779,8 +818,6 @@ void SiStripHitEfficiencyHarvester::makeSummaryVsLumi(DQMStore::IGetter& getter)
// continue
}

void SiStripHitEfficiencyHarvester::makeSummaryVsCM(DQMStore::IGetter& getter, TFileService& fs) const {}

namespace {
void setBadComponents(int i,
int comp,
Expand Down

0 comments on commit e0879e0

Please sign in to comment.