Skip to content

Commit

Permalink
EMTF DQM fix to remove duplicate x-axis bin labels
Browse files Browse the repository at this point in the history
  • Loading branch information
@nickh2000 @mmusich
nickh2000 authored and mmusich committed Nov 2, 2022
1 parent fc52d5f commit 30b35ce
Showing 1 changed file with 27 additions and 26 deletions.
53 changes: 27 additions & 26 deletions DQM/L1TMonitor/src/L1TStage2EMTF.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -58,14 +58,14 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,

//cscOccupancy designed to match the cscDQM plot
cscDQMOccupancy = ibooker.book2D("cscDQMOccupancy", "CSC Chamber Occupancy", 42, 1, 43, 20, 0, 20);
cscDQMOccupancy->setAxisTitle("10#circ Chamber (N=neighbor)", 1);
cscDQMOccupancy->setAxisTitle("10#circ Chamber (Ni = Neighbor of Sector i)", 1);
int count = 0;
for (int xbin = 1; xbin < 43; ++xbin) {
cscDQMOccupancy->setBinLabel(xbin, std::to_string(xbin - count), 1);
if (xbin == 2 || xbin == 9 || xbin == 16 || xbin == 23 || xbin == 30 || xbin == 37) {
++xbin;
++count;
cscDQMOccupancy->setBinLabel(xbin, "N", 1);
cscDQMOccupancy->setBinLabel(xbin, "N" + std::to_string(count), 1);
}
}
for (int ybin = 1; ybin <= 10; ++ybin) {
Expand Down Expand Up @@ -107,10 +107,10 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
}

rpcHitOccupancy = ibooker.book2D("rpcHitOccupancy", "RPC Chamber Occupancy", 42, 1, 43, 12, 0, 12);
rpcHitOccupancy->setAxisTitle("Sector (N=neighbor)", 1);
rpcHitOccupancy->setAxisTitle("Sector (Ni = Neighbor of Sector i)", 1);
for (int bin = 1; bin <= 6; ++bin) {
rpcHitOccupancy->setBinLabel(bin * 7 - 6, std::to_string(bin), 1);
rpcHitOccupancy->setBinLabel(bin * 7, "N", 1);
rpcHitOccupancy->setBinLabel(bin * 7, "N" + std::to_string(bin), 1);
rpcHitOccupancy->setBinLabel(bin, "RE-" + rpc_label[bin - 1], 2);
rpcHitOccupancy->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1], 2);
}
Expand All @@ -135,14 +135,14 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
gemHitBX->setBinLabel(2, "GE+1/1", 2);

gemHitOccupancy = ibooker.book2D("gemHitOccupancy", "GEM Chamber Occupancy", 42, 1, 43, 2, 0, 2);
gemHitOccupancy->setAxisTitle("10#circ Chambers (N=neighbor)", 1);
gemHitOccupancy->setAxisTitle("10#circ Chambers (Ni = Neighbor of Sector i)", 1);
count = 0;
for (int xbin = 1; xbin < 43; ++xbin) {
gemHitOccupancy->setBinLabel(xbin, std::to_string(xbin - count), 1);
if (xbin == 2 || xbin == 9 || xbin == 16 || xbin == 23 || xbin == 30 || xbin == 37) {
++xbin;
++count;
gemHitOccupancy->setBinLabel(xbin, "N", 1);
gemHitOccupancy->setBinLabel(xbin, "N" + std::to_string(count), 1);
}
}

Expand Down Expand Up @@ -409,7 +409,7 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
24,
0,
24); // 8* (0-2) phi part + (0-7) eta part
gemChamberVFATBX[hist][bx - 1]->setAxisTitle("Chamber, " + label, 1);
gemChamberVFATBX[hist][bx - 1]->setAxisTitle("Chamber, (Ni = Neighbor of Sector i), " + label, 1);
gemChamberVFATBX[hist][bx - 1]->setAxisTitle("VFAT #", 2);

for (int bin = 1; bin <= 24; bin++)
Expand All @@ -421,31 +421,32 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
if (bin == 2 || bin == 9 || bin == 16 || bin == 23 || bin == 30 || bin == 37) {
++bin;
++count;
gemChamberVFATBX[hist][bx - 1]->setBinLabel(bin, "N", 1);
gemChamberVFATBX[hist][bx - 1]->setBinLabel(bin, "N" + std::to_string(count), 1);
}
}
gemChamberVFATBX[hist][bx - 1]->getTH2F()->GetXaxis()->SetCanExtend(false);
}
}
// CSC LCT and RPC Hit Timing
ibooker.setCurrentFolder(monitorDir + "/Timing");

cscTimingTot = ibooker.book2D("cscTimingTotal", "CSC Total BX ", 42, 1, 43, 20, 0, 20);
cscTimingTot->setAxisTitle("10#circ Chamber (N=neighbor)", 1);
cscTimingTot->setAxisTitle("10#circ Chamber, (Ni = Neighbor of Sector i)", 1);

rpcHitTimingTot = ibooker.book2D("rpcHitTimingTot", "RPC Chamber Occupancy ", 42, 1, 43, 12, 0, 12);
rpcHitTimingTot->setAxisTitle("Sector (N=neighbor)", 1);
rpcHitTimingTot->setAxisTitle("Sector (Ni = Neighbor of Sector i)", 1);

gemHitTimingTot =
ibooker.book2D("gemHitTimingTot", "GEM Chamber Occupancy ", 42, 1, 43, 2, 0, 2); // Add GEM Timing Oct 27 2020
gemHitTimingTot->setAxisTitle("10#circ Chamber (N=neighbor)", 1);
gemHitTimingTot->setAxisTitle("10#circ Chamber (Ni = Neighbor of Sector i)", 1);
const std::array<std::string, 5> nameBX{{"BXNeg1", "BXPos1", "BXNeg2", "BXPos2", "BX0"}};
const std::array<std::string, 5> labelBX{{"BX -1", "BX +1", "BX -2", "BX +2", "BX 0"}};

for (int hist = 0; hist < 5; ++hist) {
count = 0;
cscLCTTiming[hist] =
ibooker.book2D("cscLCTTiming" + nameBX[hist], "CSC Chamber Occupancy " + labelBX[hist], 42, 1, 43, 20, 0, 20);
cscLCTTiming[hist]->setAxisTitle("10#circ Chamber", 1);
cscLCTTiming[hist]->setAxisTitle("10#circ Chamber, (Ni = Neighbor of Sector i)", 1);

for (int xbin = 1; xbin < 43; ++xbin) {
cscLCTTiming[hist]->setBinLabel(xbin, std::to_string(xbin - count), 1);
Expand All @@ -454,9 +455,9 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
if (xbin == 2 || xbin == 9 || xbin == 16 || xbin == 23 || xbin == 30 || xbin == 37) {
++xbin;
++count;
cscLCTTiming[hist]->setBinLabel(xbin, "N", 1);
cscLCTTiming[hist]->setBinLabel(xbin, "N" + std::to_string(count), 1);
if (hist == 0)
cscTimingTot->setBinLabel(xbin, "N", 1);
cscTimingTot->setBinLabel(xbin, "N" + std::to_string(count), 1);
}
}

Expand All @@ -474,18 +475,18 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,

rpcHitTiming[hist] =
ibooker.book2D("rpcHitTiming" + nameBX[hist], "RPC Chamber Occupancy " + labelBX[hist], 42, 1, 43, 12, 0, 12);
rpcHitTiming[hist]->setAxisTitle("Sector (N=neighbor)", 1);
rpcHitTiming[hist]->setAxisTitle("Sector, (Ni=Neighbor of Sector i )", 1);
for (int bin = 1; bin < 7; ++bin) {
rpcHitTiming[hist]->setBinLabel(bin * 7 - 6, std::to_string(bin), 1);
rpcHitTiming[hist]->setBinLabel(bin * 7, "N", 1);
rpcHitTiming[hist]->setBinLabel(bin * 7, "N" + std::to_string(bin), 1);
rpcHitTiming[hist]->setBinLabel(bin, "RE-" + rpc_label[bin - 1], 2);
rpcHitTiming[hist]->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1], 2);
}
rpcHitTiming[hist]->getTH2F()->GetXaxis()->SetCanExtend(false); // Needed to stop multi-thread summing
if (hist == 0) {
for (int bin = 1; bin < 7; ++bin) {
rpcHitTimingTot->setBinLabel(bin * 7 - 6, std::to_string(bin), 1);
rpcHitTimingTot->setBinLabel(bin * 7, "N", 1);
rpcHitTimingTot->setBinLabel(bin * 7, "N" + std::to_string(bin), 1);
rpcHitTimingTot->setBinLabel(bin, "RE-" + rpc_label[bin - 1], 2);
rpcHitTimingTot->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1], 2);
}
Expand All @@ -496,7 +497,7 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
// Add GEM Timing Oct 27 2020
gemHitTiming[hist] =
ibooker.book2D("gemHitTiming" + nameBX[hist], "GEM Chamber Occupancy " + labelBX[hist], 42, 1, 43, 2, 0, 2);
gemHitTiming[hist]->setAxisTitle("10#circ Chamber", 1);
gemHitTiming[hist]->setAxisTitle("10#circ Chamber, (Ni = Neighbor of Sector i)", 1);
count = 0;
for (int xbin = 1; xbin < 43; ++xbin) {
gemHitTiming[hist]->setBinLabel(xbin, std::to_string(xbin - count), 1);
Expand All @@ -505,9 +506,9 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
if (xbin == 2 || xbin == 9 || xbin == 16 || xbin == 23 || xbin == 30 || xbin == 37) {
++xbin;
++count;
gemHitTiming[hist]->setBinLabel(xbin, "N", 1);
gemHitTiming[hist]->setBinLabel(xbin, "N" + std::to_string(count), 1);
if (hist == 0)
gemHitTimingTot->setBinLabel(xbin, "N", 1);
gemHitTimingTot->setBinLabel(xbin, "N" + std::to_string(count), 1);
}
}

Expand All @@ -523,13 +524,13 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
count = 0;
cscLCTTimingFrac[hist] = ibooker.book2D(
"cscLCTTimingFrac" + nameBX[hist], "CSC Chamber Occupancy " + labelBX[hist], 42, 1, 43, 20, 0, 20);
cscLCTTimingFrac[hist]->setAxisTitle("10#circ Chambers", 1);
cscLCTTimingFrac[hist]->setAxisTitle("10#circ Chambers, (Ni = Neighbor of Sector i)", 1);
for (int xbin = 1; xbin < 43; ++xbin) {
cscLCTTimingFrac[hist]->setBinLabel(xbin, std::to_string(xbin - count), 1);
if (xbin == 2 || xbin == 9 || xbin == 16 || xbin == 23 || xbin == 30 || xbin == 37) {
++xbin;
++count;
cscLCTTimingFrac[hist]->setBinLabel(xbin, "N", 1);
cscLCTTimingFrac[hist]->setBinLabel(xbin, "N" + std::to_string(count), 1);
}
}
for (int ybin = 1; ybin <= 10; ++ybin) {
Expand All @@ -540,25 +541,25 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,

rpcHitTimingFrac[hist] = ibooker.book2D(
"rpcHitTimingFrac" + nameBX[hist], "RPC Chamber Fraction in " + labelBX[hist], 42, 1, 43, 12, 0, 12);
rpcHitTimingFrac[hist]->setAxisTitle("Sector (N=neighbor)", 1);
rpcHitTimingFrac[hist]->setAxisTitle("Sector, (Ni = Neighbor of Sector i)", 1);
for (int bin = 1; bin < 7; ++bin) {
rpcHitTimingFrac[hist]->setBinLabel(bin * 7 - 6, std::to_string(bin), 1);
rpcHitTimingFrac[hist]->setBinLabel(bin * 7, "N", 1);
rpcHitTimingFrac[hist]->setBinLabel(bin * 7, "N" + std::to_string(bin), 1);
rpcHitTimingFrac[hist]->setBinLabel(bin, "RE-" + rpc_label[bin - 1], 2);
rpcHitTimingFrac[hist]->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1], 2);
}

// Add GEM Timing Oct 27 2020
gemHitTimingFrac[hist] =
ibooker.book2D("gemHitTimingFrac" + nameBX[hist], "GEM Chamber Occupancy " + labelBX[hist], 42, 1, 43, 2, 0, 2);
gemHitTimingFrac[hist]->setAxisTitle("10#circ Chambers", 1);
gemHitTimingFrac[hist]->setAxisTitle("10#circ Chambers, (Ni = Neighbor of Sector i)", 1);
count = 0;
for (int xbin = 1; xbin < 43; ++xbin) {
gemHitTimingFrac[hist]->setBinLabel(xbin, std::to_string(xbin - count), 1);
if (xbin == 2 || xbin == 9 || xbin == 16 || xbin == 23 || xbin == 30 || xbin == 37) {
++xbin;
++count;
gemHitTimingFrac[hist]->setBinLabel(xbin, "N", 1);
gemHitTimingFrac[hist]->setBinLabel(xbin, "N" + std::to_string(count), 1);
}
}
gemHitTimingFrac[hist]->setBinLabel(1, "GE-1/1", 2);
Expand Down

0 comments on commit 30b35ce

Please sign in to comment.