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

EMTF DQM plot updates #21459

Merged
merged 3 commits into from Dec 12, 2017
Merged

EMTF DQM plot updates #21459

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
11c7990
EMTF DQM modified CSC LCT and RPC hit occupancy plots
abrinke1 Nov 24, 2017
f8bd914
EMTF DQM fix typos
abrinke1 Nov 24, 2017
35fd4a9
EMTF DQM indent fixes, simplify chamber_bin
abrinke1 Nov 27, 2017
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
10 changes: 8 additions & 2 deletions DQM/Integration/python/clients/l1tstage2_dqm_sourceclient-live_cfg.py
View file
Open in desktop
Expand Up @@ -9,12 +9,18 @@
# Live Online DQM in P5
process.load("DQM.Integration.config.inputsource_cfi")

# Testing in lxplus
#process.load("DQM.Integration.config.fileinputsource_cfi")
# # Testing in lxplus
# process.load("DQM.Integration.config.fileinputsource_cfi")
# process.load("FWCore.MessageLogger.MessageLogger_cfi")
# process.MessageLogger.cerr.FwkReport.reportEvery = 1

# Required to load Global Tag
process.load("DQM.Integration.config.FrontierCondition_GT_cfi")

# # Condition for lxplus: change and possibly customise the GT
# from Configuration.AlCa.GlobalTag import GlobalTag as gtCustomise
# process.GlobalTag = gtCustomise(process.GlobalTag, 'auto:run2_data', '')

# Required to load EcalMappingRecord
process.load("Configuration.StandardSequences.GeometryRecoDB_cff")

Expand Down
9 changes: 6 additions & 3 deletions DQM/L1TMonitor/interface/L1TStage2EMTF.h
View file
Open in desktop
Expand Up @@ -47,9 +47,11 @@ class L1TStage2EMTF : public DQMEDAnalyzer {
MonitorElement* cscChamberStrip[20];
MonitorElement* cscChamberWire[20];
MonitorElement* cscLCTOccupancy;
MonitorElement* cscDQMOccupancy;
MonitorElement* cscLCTTiming[5];
MonitorElement* cscLCTTimingFrac[4];

MonitorElement* cscLCTTimingFrac[5];
MonitorElement* cscTimingTot;

MonitorElement* emtfnTracks;
MonitorElement* emtfTracknHits;
MonitorElement* emtfTrackBX;
Expand All @@ -73,7 +75,8 @@ class L1TStage2EMTF : public DQMEDAnalyzer {
MonitorElement* rpcHitBX;
MonitorElement* rpcHitOccupancy;
MonitorElement* rpcHitTiming[5];
MonitorElement* rpcHitTimingFrac[4];
MonitorElement* rpcHitTimingFrac[5];
MonitorElement* rpcHitTimingTot;
MonitorElement* rpcHitPhi[12];
MonitorElement* rpcHitTheta[12];
MonitorElement* rpcChamberPhi[12];
Expand Down
195 changes: 139 additions & 56 deletions DQM/L1TMonitor/src/L1TStage2EMTF.cc
View file
Open in desktop
Expand Up @@ -59,7 +59,24 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
for (unsigned int ybin = 0; ybin < binNames.size(); ++ybin) {
cscLCTOccupancy->setBinLabel(ybin+1, binNames[ybin], 2);
}


//cscOccupancy designed to match the cscDQM plot
cscDQMOccupancy = ibooker.book2D("cscDQMOccupancy", "CSC Chamber Occupancy", 42, 1, 43, 20, 0, 20);
cscDQMOccupancy->setAxisTitle("10 degree Chamber", 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);
}
}
for (int ybin = 1; ybin <= 10; ++ybin) {
cscDQMOccupancy->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
cscDQMOccupancy->setBinLabel(21 - ybin, "ME+" + suffix_label[ybin - 1], 2);
}

mpcLinkErrors = ibooker.book2D("mpcLinkErrors", "MPC Link Errors", 54, 1, 55, 12, -6, 6);
mpcLinkErrors->setAxisTitle("Sector (CSCID 1-9 Unlabelled)", 1);
for (int xbin = 1; xbin < 7; ++xbin) {
Expand Down Expand Up @@ -92,10 +109,11 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
rpcHitBX->setBinLabel(13 - ybin, "RE+" + rpc_label[ybin - 1], 2);
}

rpcHitOccupancy = ibooker.book2D("rpcHitOccupancy", "RPC Chamber Occupancy", 36, 1, 37, 12, 0, 12);
rpcHitOccupancy = ibooker.book2D("rpcHitOccupancy", "RPC Chamber Occupancy", 42, 1, 43, 12, 0, 12);
rpcHitOccupancy->setAxisTitle("Sector", 1);
for (int bin = 1; bin < 7; ++bin) {
rpcHitOccupancy->setBinLabel(bin*6 - 5, std::to_string(bin), 1);
rpcHitOccupancy->setBinLabel(bin*7 - 6, std::to_string(bin), 1);
rpcHitOccupancy->setBinLabel(bin*7, "N", 1);
rpcHitOccupancy->setBinLabel(bin, "RE-" + rpc_label[bin - 1], 2);
rpcHitOccupancy->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1],2);
}
Expand Down Expand Up @@ -171,9 +189,9 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,

if (hist < 6) {
nChambs = (i % 2) ? 18 : 36;
} else if (hist > 13) {
nChambs = (i % 2) ? 36 : 18;
} else {
} else if (hist >13) {
nChambs = (i % 2) ? 36 : 18;
} else {
nChambs = 36;
}

Expand Down Expand Up @@ -237,52 +255,85 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,

// 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 degree Chambers", 1);

rpcHitTimingTot = ibooker.book2D("rpcHitTimingTot", "RPC Chamber Occupancy ", 42, 1, 43, 12, 0, 12);
rpcHitTimingTot->setAxisTitle("Sector", 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) {

cscLCTTiming[hist] = ibooker.book2D("cscLCTTiming" + nameBX[hist], "CSC Chamber Occupancy " + labelBX[hist], 54, 1, 55, 12, -6, 6);
cscLCTTiming[hist]->setAxisTitle("Sector (CSCID 1-9 Unlabelled)", 1);
for (int xbin = 1; xbin < 7; ++xbin) {
cscLCTTiming[hist]->setBinLabel(xbin * 9 - 8, std::to_string(xbin), 1);
count = 0;
cscLCTTiming[hist] = ibooker.book2D("cscLCTTiming" + nameBX[hist], "CSC Chamber Occupancy " + labelBX[hist], 42, 1, 43, 20, 0, 20);
cscLCTTiming[hist]->setAxisTitle("10 degree Chambers", 1);

for (int xbin=1; xbin < 43; ++xbin) {
cscLCTTiming[hist]->setBinLabel(xbin, std::to_string(xbin-count), 1);
if (hist==0) cscTimingTot->setBinLabel(xbin, std::to_string(xbin-count), 1);//only fill once in the loop
if (xbin==2 || xbin==9 || xbin==16 || xbin==23 || xbin==30 ||xbin==37 ) {
++xbin;
++count;
cscLCTTiming[hist]->setBinLabel(xbin, "N", 1);
if (hist==0) cscTimingTot->setBinLabel(xbin, "N", 1);
}
}
for (unsigned int ybin = 0; ybin < binNames.size(); ++ybin) {
cscLCTTiming[hist]->setBinLabel(ybin+1, binNames[ybin], 2);

for (int ybin = 1; ybin <= 10; ++ybin) {
cscLCTTiming[hist]->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
cscLCTTiming[hist]->setBinLabel(21 - ybin, "ME+" + suffix_label[ybin - 1], 2);
if (hist==0) cscTimingTot->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
if (hist==0) cscTimingTot->setBinLabel(21 - ybin, "ME+" + suffix_label[ybin - 1], 2);
}

rpcHitTiming[hist] = ibooker.book2D("rpcHitTiming" + nameBX[hist], "RPC Chamber Occupancy " + labelBX[hist], 36, 1, 37, 12, 0, 12);
rpcHitTiming[hist] = ibooker.book2D("rpcHitTiming" + nameBX[hist], "RPC Chamber Occupancy " + labelBX[hist], 42, 1, 43, 12, 0, 12);
rpcHitTiming[hist]->setAxisTitle("Sector", 1);
for (int bin = 1; bin < 7; ++bin) {
rpcHitTiming[hist]->setBinLabel(bin*6 - 5, std::to_string(bin), 1);
rpcHitTiming[hist]->setBinLabel(bin*7 - 6, std::to_string(bin), 1);
rpcHitTiming[hist]->setBinLabel(bin*7, "N", 1);
rpcHitTiming[hist]->setBinLabel(bin, "RE-" + rpc_label[bin - 1], 2);
rpcHitTiming[hist]->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1],2);
}

if (hist == 4) continue; // Don't book for BX = 0

cscLCTTimingFrac[hist] = ibooker.book2D("cscLCTTimingFrac" + nameBX[hist], "CSC Chamber Fraction in " + labelBX[hist], 54, 1, 55, 12, -6, 6);
cscLCTTimingFrac[hist]->setEfficiencyFlag();
cscLCTTimingFrac[hist]->setAxisTitle("Sector (CSCID 1-9 Unlabelled)", 1);
for (int xbin = 1; xbin < 7; ++xbin) {
cscLCTTimingFrac[hist]->setBinLabel(xbin * 9 - 8, std::to_string(xbin), 1);
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, "RE-" + rpc_label[bin - 1], 2);
rpcHitTimingTot->setBinLabel(13 - bin, "RE+" + rpc_label[bin - 1],2);
}
}
for (unsigned int ybin = 0; ybin < binNames.size(); ++ybin) {
cscLCTTimingFrac[hist]->setBinLabel(ybin+1, binNames[ybin], 2);
//if (hist == 4) continue; // Don't book for BX = 0

count = 0;
cscLCTTimingFrac[hist] = ibooker.book2D("cscLCTTimingFrac" + nameBX[hist], "CSC Chamber Occupancy " + labelBX[hist], 42, 1, 43, 20, 0, 20);
cscLCTTimingFrac[hist]->setAxisTitle("10 degree Chambers", 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);
}
}

rpcHitTimingFrac[hist] = ibooker.book2D("rpcHitTimingFrac" + nameBX[hist], "RPC Chamber Fraction in " + labelBX[hist], 36, 1, 37, 12, 0, 12);
rpcHitTimingFrac[hist]->setEfficiencyFlag();
for (int ybin = 1; ybin <= 10; ++ybin) {
cscLCTTimingFrac[hist]->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
cscLCTTimingFrac[hist]->setBinLabel(21 - ybin, "ME+" + suffix_label[ybin - 1], 2);
}

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

}

rpcHitTimingInTrack = ibooker.book2D("rpcHitTimingInTrack", "RPC Hit Timing (matched to BX 0 track)", 7, -3, 4, 12, 0, 12);
rpcHitTimingInTrack = ibooker.book2D("rpcHitTimingInTrack", "RPC Hit Timing (matched to track in BX 0)", 7, -3, 4, 12, 0, 12);
rpcHitTimingInTrack->setAxisTitle("BX", 1);
for (int xbin = 1, xbin_label = -3; xbin <= 7; ++xbin, ++xbin_label) {
rpcHitTimingInTrack->setBinLabel(xbin, std::to_string(xbin_label), 1);
Expand Down Expand Up @@ -340,6 +391,18 @@ void L1TStage2EMTF::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
}
}

// CSCOccupancy chamber mapping for neighbor inclusive plots
int chamber_bin (int station, int ring, int chamber) {
int chamber_bin_index = 0;
if (station > 1 && (ring % 2) == 1) {
chamber_bin_index = (chamber * 2) + ((chamber + 1) / 3);
} else {
chamber_bin_index = chamber + ((chamber + 3) / 6);
}
return chamber_bin_index;
};


void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {

if (verbose) edm::LogInfo("L1TStage2EMTF") << "L1TStage2EMTF: analyze..." << std::endl;
Expand Down Expand Up @@ -428,7 +491,13 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
hist_index = histIndexCSC.at( {station, ring} );
if (endcap > 0) hist_index = 19 - hist_index;
cscLCTBX->Fill(Hit->BX(), hist_index);
float evt_wgt = (Hit->Station() > 1 && Hit->Ring() == 1) ? 0.5 : 1.0;
if (Hit->Neighbor() == false) {
//Map for cscDQMOccupancy plot
cscDQMOccupancy->Fill(chamber_bin(station,ring,chamber),hist_index,evt_wgt);
if (station>1 && (ring % 2)==1) {
cscDQMOccupancy->Fill(chamber_bin(station,ring,chamber)-1,hist_index,evt_wgt);
}
cscLCTStrip[hist_index]->Fill(strip);
cscLCTWire[hist_index]->Fill(wire);
cscChamberStrip[hist_index]->Fill(chamber, strip);
Expand All @@ -442,7 +511,10 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
// Map neighbor chambers to "fake" CSC IDs: 1/3 --> 1, 1/6 --> 2, 1/9 --> 3, 2/3 --> 4, 2/9 --> 5, etc.
int cscid_n = (station == 1 ? (cscid / 3) : (station * 2) + ((cscid - 3) / 6) );
cscLCTOccupancy->Fill(cscid_n + cscid_offset, endcap * 5.5);
}
}
if (Hit->Neighbor() == true) {
cscDQMOccupancy->Fill(sector*7-4,hist_index,evt_wgt);
}
}

if (Hit->Is_RPC() == true) {
Expand All @@ -456,8 +528,10 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
rpcHitTheta[hist_index]->Fill(Hit->Theta_fp() / 4);
rpcChamberPhi[hist_index]->Fill(chamber, Hit->Phi_fp() / 4);
rpcChamberTheta[hist_index]->Fill(chamber, Hit->Theta_fp() / 4);
rpcHitOccupancy->Fill((Hit->Sector_RPC() - 1) * 6 + Hit->Subsector(), hist_index + 0.5);
}
rpcHitOccupancy->Fill((Hit->Sector_RPC() - 1) * 7 + Hit->Subsector(), hist_index + 0.5);
} else if (Hit->Neighbor() == true) {
rpcHitOccupancy->Fill((Hit->Sector_RPC() - 1) * 7 + 7, hist_index + 0.5);
}
}
}

Expand Down Expand Up @@ -505,14 +579,14 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
if (numHits < 2 || numHits > 4) continue;
l1t::EMTFHitCollection tmp_hits = Track->Hits();
int numHitsInTrack_BX0 = 0;
unsigned int hist_index = 4 - numHits;
unsigned int hist_index2 = 4 - numHits;

for (const auto & iTrkHit: Track->Hits()) {
if (iTrkHit.Is_CSC() == true) {
emtfTrackBXVsCSCLCT[hist_index]->Fill(iTrkHit.BX(), Track->BX());
emtfTrackBXVsCSCLCT[hist_index2]->Fill(iTrkHit.BX(), Track->BX());
}
else if (iTrkHit.Is_RPC() == true) {
emtfTrackBXVsRPCHit[hist_index]->Fill(iTrkHit.BX(), Track->BX());
emtfTrackBXVsRPCHit[hist_index2]->Fill(iTrkHit.BX(), Track->BX());
}
}

Expand All @@ -527,34 +601,38 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
for (const auto & TrkHit: Track->Hits()) {

int trackHitBX = TrkHit.BX();
int cscid = TrkHit.CSC_ID();
//int cscid = TrkHit.CSC_ID();
int ring = TrkHit.Ring();
int station = TrkHit.Station();
int sector = TrkHit.Sector();
int subsector = TrkHit.Subsector();
int cscid_offset = (sector - 1) * 9;
//int cscid_offset = (sector - 1) * 9;//no longer needed after new time plots (maybe useful for future plots)
int neighbor = TrkHit.Neighbor();
int endcap = TrkHit.Endcap();
int chamber = TrkHit.Chamber();

int hist_index = 0;
int hist_index = 0;
float evt_wgt = (TrkHit.Station() > 1 && TrkHit.Ring() == 1) ? 0.5 : 1.0;
// Maps CSC BX from -2 to 2 to monitor element cscLCTTIming
const std::map<int, int> histIndexBX = {{0, 4}, {-1, 0}, {1, 1}, {-2, 2}, {2, 3}};
if (std::abs(trackHitBX) > 2) continue; // Should never happen, but just to be safe ...

if (TrkHit.Is_CSC() == true) {
hist_index = histIndexCSC.at( {station, ring} );
if (endcap > 0) hist_index = 19 - hist_index;
if (neighbor == false) {
if (subsector == 1) {
cscLCTTiming[histIndexBX.at(trackHitBX)]->Fill(cscid + cscid_offset, endcap * (station - 0.5));
}
else {
cscLCTTiming[histIndexBX.at(trackHitBX)]->Fill(cscid + cscid_offset, endcap * (station + 0.5));
}
cscLCTTiming[histIndexBX.at(trackHitBX)]->Fill(chamber_bin(station,ring,chamber),hist_index,evt_wgt);
cscTimingTot->Fill(chamber_bin(station,ring,chamber),hist_index,evt_wgt);
if (station>1 && (ring % 2)==1) {
cscLCTTiming[histIndexBX.at(trackHitBX)]->Fill(chamber_bin(station,ring,chamber)-1,hist_index,evt_wgt);
cscTimingTot->Fill(chamber_bin(station,ring,chamber)-1,hist_index,evt_wgt);
}
}
else {
// Map neighbor chambers to "fake" CSC IDs: 1/3 --> 1, 1/6 --> 2, 1/9 --> 3, 2/3 --> 4, 2/9 --> 5, etc.
int cscid_n = (station == 1 ? (cscid / 3) : (station * 2) + ((cscid - 3) / 6) );
cscLCTTiming[histIndexBX.at(trackHitBX)]->Fill(cscid_n + cscid_offset, endcap * 5.5);
//int cscid_n = (station == 1 ? (cscid / 3) : (station * 2) + ((cscid - 3) / 6) );
cscLCTTiming[histIndexBX.at(trackHitBX)]->Fill(sector*7-4,hist_index,evt_wgt);
cscTimingTot->Fill(sector*7-4,hist_index,evt_wgt);
}

// Fill RPC timing with matched CSC LCTs
Expand All @@ -569,9 +647,9 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {

hist_index = histIndexRPC.at( {Hit->Station(), Hit->Ring()} );
if (Hit->Endcap() > 0) hist_index = 11 - hist_index;
rpcHitTimingInTrack->Fill(Hit->BX(), hist_index + 0.5);
rpcHitTiming[histIndexBX.at(Hit->BX())]->Fill((Hit->Sector_RPC() - 1) * 6 + Hit->Subsector(), hist_index + 0.5);
rpcHitTimingInTrack->Fill(Hit->BX(), hist_index + 0.5);
rpcHitTiming[histIndexBX.at(Hit->BX())]->Fill((Hit->Sector_RPC() - 1) * 7 + Hit->Subsector(), hist_index + 0.5);
rpcHitTimingTot->Fill((Hit->Sector_RPC() - 1) * 7 + Hit->Subsector(), hist_index + 0.5);
} // End loop: for (auto Hit = HitCollection->begin(); Hit != HitCollection->end(); ++Hit)
} // End conditional: if (trackHitBX == 0 && ring == 2)
} // End conditional: if (TrkHit.Is_CSC() == true)
Expand All @@ -584,9 +662,14 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
if (endcap > 0) hist_index = 11 - hist_index;

rpcHitTimingInTrack->Fill(trackHitBX, hist_index + 0.5);
rpcHitTiming[histIndexBX.at(trackHitBX)]->Fill((TrkHit.Sector_RPC() - 1) * 6 + subsector, hist_index + 0.5);
rpcHitTiming[histIndexBX.at(trackHitBX)]->Fill((TrkHit.Sector_RPC() - 1) * 7 + subsector, hist_index + 0.5);
rpcHitTimingTot->Fill((TrkHit.Sector_RPC() - 1) * 7 + subsector, hist_index + 0.5);
} // End conditional: if (TrkHit.Is_RPC() == true && neighbor == false)

if (TrkHit.Is_RPC() == true && neighbor == true) {
hist_index = histIndexRPC.at( {station, ring} );
if (endcap > 0) hist_index = 11 - hist_index;
rpcHitTiming[histIndexBX.at(trackHitBX)]->Fill((TrkHit.Sector_RPC() - 1) * 7, hist_index + 0.5);
}
} // End loop: for (int iHit = 0; iHit < numHits; ++iHit)
}
//////////////////////////////////////////////////
Expand All @@ -596,9 +679,9 @@ void L1TStage2EMTF::analyze(const edm::Event& e, const edm::EventSetup& c) {
} // End loop: for (auto Track = TrackCollection->begin(); Track != TrackCollection->end(); ++Track)

// CSC LCT and RPC Hit Timing Efficieny
for (int hist_index = 0; hist_index < 4; ++hist_index) {
cscLCTTimingFrac[hist_index]->getTH2F()->Divide(cscLCTTiming[hist_index]->getTH2F(), cscLCTTiming[4]->getTH2F());
rpcHitTimingFrac[hist_index]->getTH2F()->Divide(rpcHitTiming[hist_index]->getTH2F(), rpcHitTiming[4]->getTH2F());
for (int hist_index = 0; hist_index < 5; ++hist_index) {
cscLCTTimingFrac[hist_index]->getTH2F()->Divide(cscLCTTiming[hist_index]->getTH2F(), cscTimingTot->getTH2F());
rpcHitTimingFrac[hist_index]->getTH2F()->Divide(rpcHitTiming[hist_index]->getTH2F(), rpcHitTimingTot->getTH2F());
}

// Regional Muon Candidates
Expand Down