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

Backport muon shower (un)packing and DQM to 12_2_X #37225

Merged
merged 18 commits into from
Mar 14, 2022
Merged

Backport muon shower (un)packing and DQM to 12_2_X #37225

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
18 commits
Select commit Hold shift + click to select a range
0e7d663
Use the same interface to set TF identifiers for RegionalMuonShower a…
dinyar Sep 23, 2021
892b3db
Add (un)packing functions for RegionalMuonShower
dinyar Sep 23, 2021
2f3ae1c
Add tokens and collections for unpacked muon showers
dinyar Sep 23, 2021
200ae97
(Un)pack regional muon showers
dinyar Sep 23, 2021
37281c1
(Un)packing of hadronic showers at the uGMT output
dinyar Sep 27, 2021
925ab4f
Add muon shower emulators to Phase2 sequence
dinyar Sep 29, 2021
1133d2e
Explicitly name constant indicating number of output copies
dinyar Dec 7, 2021
d9783db
Update shower (un)packing function with new specification
dinyar Dec 12, 2021
41b46f6
Fix error inhibiting muon shower packing
dinyar Jan 19, 2022
9a78fe0
Select correct word in payload for muon showers
dinyar Jan 26, 2022
a8c6d4f
Don't write Run-2 EMTF debug data for Run-3
dinyar Jan 26, 2022
dafa165
Only pack muon showers if they're valid
dinyar Jan 26, 2022
f046d30
Fill old fields with something to avoid L1TNtuple crashes
dinyar Jan 26, 2022
69399a9
Make MuonShower data format more descriptive
dinyar Feb 1, 2022
841964f
Enable unpacking of uGMT output muon showers
dinyar Jan 31, 2022
bca02af
Enable shower unpacking for the uGT inputs
dinyar Jan 31, 2022
a30d4d8
Only fill RegionalMuonShower if it's valid
dinyar Feb 1, 2022
26727cc
Add muon shower plots to uGMT online DQM
dinyar Jan 18, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
29 changes: 22 additions & 7 deletions DQM/L1TMonitor/interface/L1TStage2uGMT.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,6 +3,8 @@

#include "DataFormats/L1Trigger/interface/Muon.h"
#include "DataFormats/L1TMuon/interface/RegionalMuonCand.h"
#include "DataFormats/L1Trigger/interface/MuonShower.h"
#include "DataFormats/L1TMuon/interface/RegionalMuonShower.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTConfiguration.h"

#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
Expand All @@ -27,14 +29,17 @@ class L1TStage2uGMT : public DQMEDAnalyzer {
private:
l1t::tftype getTfOrigin(const int tfMuonIndex);

edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> ugmtBMTFToken;
edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> ugmtOMTFToken;
edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> ugmtEMTFToken;
edm::EDGetTokenT<l1t::MuonBxCollection> ugmtMuonToken;
std::string monitorDir;
bool emul;
bool verbose;
edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> ugmtBMTFToken_;
edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> ugmtOMTFToken_;
edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> ugmtEMTFToken_;
edm::EDGetTokenT<l1t::MuonBxCollection> ugmtMuonToken_;
edm::EDGetTokenT<l1t::RegionalMuonShowerBxCollection> ugmtEMTFShowerToken_;
edm::EDGetTokenT<l1t::MuonShowerBxCollection> ugmtMuonShowerToken_;
std::string monitorDir_;
bool emul_;
bool verbose_;
bool displacedQuantities_;
bool hadronicShowers_;

const float etaScale_;
const float phiScale_;
Expand Down Expand Up @@ -94,6 +99,10 @@ class L1TStage2uGMT : public DQMEDAnalyzer {
MonitorElement* ugmtEMTFMuMuDPhi;
MonitorElement* ugmtEMTFMuMuDR;

MonitorElement* ugmtEMTFShowerTypeOccupancyPerSector;
MonitorElement* ugmtEMTFShowerTypeOccupancyPerBx;
MonitorElement* ugmtEMTFShowerSectorOccupancyPerBx;

MonitorElement* ugmtBOMTFposMuMuDEta;
MonitorElement* ugmtBOMTFposMuMuDPhi;
MonitorElement* ugmtBOMTFposMuMuDR;
Expand Down Expand Up @@ -161,6 +170,9 @@ class L1TStage2uGMT : public DQMEDAnalyzer {
MonitorElement* ugmtMuonBXvshwIso;
MonitorElement* ugmtMuonChargevsLink;

// Output shower plots
MonitorElement* ugmtMuonShowerTypeOccupancyPerBx;

// muon correlations
MonitorElement* ugmtMuMuInvMass;
MonitorElement* ugmtMuMuInvMassAtVtx;
Expand Down Expand Up @@ -200,6 +212,9 @@ class L1TStage2uGMT : public DQMEDAnalyzer {
MonitorElement* ugmtMuMuDEtaEneg;
MonitorElement* ugmtMuMuDPhiEneg;
MonitorElement* ugmtMuMuDREneg;

static constexpr unsigned IDX_TIGHT_SHOWER{2};
static constexpr unsigned IDX_NOMINAL_SHOWER{1};
};

#endif
7 changes: 5 additions & 2 deletions DQM/L1TMonitor/python/L1TStage2uGMT_cfi.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,13 +7,16 @@
bmtfProducer = cms.InputTag("gmtStage2Digis", "BMTF"),
omtfProducer = cms.InputTag("gmtStage2Digis", "OMTF"),
emtfProducer = cms.InputTag("gmtStage2Digis", "EMTF"),
emtfShowerProducer = cms.InputTag("gmtStage2Digis", "EMTF"),
muonProducer = cms.InputTag("gmtStage2Digis", "Muon"),
muonShowerProducer = cms.InputTag("gmtStage2Digis", "MuonShower"),
monitorDir = cms.untracked.string("L1T/L1TStage2uGMT"),
emulator = cms.untracked.bool(False),
verbose = cms.untracked.bool(False),
displacedQuantities = cms.untracked.bool(False)
displacedQuantities = cms.untracked.bool(False),
hadronicShowers = cms.untracked.bool(False)
)

## Era: Run3_2021; Displaced muons from BMTF used in uGMT from Run-3
from Configuration.Eras.Modifier_stage2L1Trigger_2021_cff import stage2L1Trigger_2021
stage2L1Trigger_2021.toModify(l1tStage2uGMT, displacedQuantities = cms.untracked.bool(True))
stage2L1Trigger_2021.toModify(l1tStage2uGMT, displacedQuantities = cms.untracked.bool(True), hadronicShowers = cms.untracked.bool(True))
4 changes: 2 additions & 2 deletions DQM/L1TMonitor/src/L1TStage2Shower.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -166,8 +166,8 @@ void L1TStage2Shower::analyze(const edm::Event& e, const edm::EventSetup& c) {
if (not Shower.isValid())
continue;
if (Shower.isOneNominalInTime() or Shower.isTwoLooseInTime() or Shower.isOneTightInTime()) {
int endcap = Shower.endcap();
int sector = Shower.sector();
int endcap = Shower.trackFinderType() == l1t::tftype::emtf_pos ? 1 : -1;
int sector = Shower.processor() + 1;
if (Shower.isOneTightInTime())
emtfShowerTypeOccupancy->Fill(sector, (endcap == 1) ? 7.5 : 0.5);
if (Shower.isTwoLooseInTime())
Expand Down
164 changes: 141 additions & 23 deletions DQM/L1TMonitor/src/L1TStage2uGMT.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,18 +1,22 @@
#include "DQM/L1TMonitor/interface/L1TStage2uGMT.h"

L1TStage2uGMT::L1TStage2uGMT(const edm::ParameterSet& ps)
: ugmtMuonToken(consumes<l1t::MuonBxCollection>(ps.getParameter<edm::InputTag>("muonProducer"))),
monitorDir(ps.getUntrackedParameter<std::string>("monitorDir")),
emul(ps.getUntrackedParameter<bool>("emulator")),
verbose(ps.getUntrackedParameter<bool>("verbose")),
: ugmtMuonToken_(consumes<l1t::MuonBxCollection>(ps.getParameter<edm::InputTag>("muonProducer"))),
ugmtMuonShowerToken_(consumes<l1t::MuonShowerBxCollection>(ps.getParameter<edm::InputTag>("muonShowerProducer"))),
monitorDir_(ps.getUntrackedParameter<std::string>("monitorDir")),
emul_(ps.getUntrackedParameter<bool>("emulator")),
verbose_(ps.getUntrackedParameter<bool>("verbose")),
displacedQuantities_(ps.getUntrackedParameter<bool>("displacedQuantities")),
hadronicShowers_(ps.getUntrackedParameter<bool>("hadronicShowers")),
etaScale_(0.010875), // eta scale (CMS DN-2015/017)
phiScale_(0.010908) // phi scale (2*pi/576 HW values)
{
if (!emul) {
ugmtBMTFToken = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("bmtfProducer"));
ugmtOMTFToken = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("omtfProducer"));
ugmtEMTFToken = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("emtfProducer"));
if (!emul_) {
ugmtBMTFToken_ = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("bmtfProducer"));
ugmtOMTFToken_ = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("omtfProducer"));
ugmtEMTFToken_ = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("emtfProducer"));
ugmtEMTFShowerToken_ =
consumes<l1t::RegionalMuonShowerBxCollection>(ps.getParameter<edm::InputTag>("emtfShowerProducer"));
}
}

Expand All @@ -21,23 +25,26 @@ L1TStage2uGMT::~L1TStage2uGMT() {}
void L1TStage2uGMT::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
desc.add<edm::InputTag>("muonProducer")->setComment("uGMT output muons.");
;

desc.add<edm::InputTag>("bmtfProducer")->setComment("RegionalMuonCands from BMTF.");
desc.add<edm::InputTag>("omtfProducer")->setComment("RegionalMuonCands from OMTF.");
desc.add<edm::InputTag>("emtfProducer")->setComment("RegionalMuonCands from EMTF.");
desc.add<edm::InputTag>("muonShowerProducer")->setComment("uGMT output showers.");
desc.add<edm::InputTag>("emtfShowerProducer")->setComment("RegionalMuonShowers from EMTF.");
desc.addUntracked<std::string>("monitorDir", "")
->setComment("Target directory in the DQM file. Will be created if not existing.");
desc.addUntracked<bool>("emulator", false)
->setComment("Create histograms for muonProducer input only. xmtfProducer inputs are ignored.");
desc.addUntracked<bool>("verbose", false);
desc.addUntracked<bool>("displacedQuantities", false);
desc.addUntracked<bool>("hadronicShowers", false);
descriptions.add("l1tStage2uGMT", desc);
}

void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&, const edm::EventSetup&) {
if (!emul) {
if (!emul_) {
// BMTF Input
ibooker.setCurrentFolder(monitorDir + "/BMTFInput");
ibooker.setCurrentFolder(monitorDir_ + "/BMTFInput");

ugmtBMTFBX = ibooker.book1D("ugmtBMTFBX", "uGMT BMTF Input BX", 7, -3.5, 3.5);
ugmtBMTFBX->setAxisTitle("BX", 1);
Expand Down Expand Up @@ -98,7 +105,7 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
ugmtBMTFMuMuDR->setAxisTitle("#DeltaR", 1);

// OMTF Input
ibooker.setCurrentFolder(monitorDir + "/OMTFInput");
ibooker.setCurrentFolder(monitorDir_ + "/OMTFInput");

ugmtOMTFBX = ibooker.book1D("ugmtOMTFBX", "uGMT OMTF Input BX", 7, -3.5, 3.5);
ugmtOMTFBX->setAxisTitle("BX", 1);
Expand Down Expand Up @@ -165,7 +172,7 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
ugmtOMTFMuMuDR->setAxisTitle("#DeltaR", 1);

// EMTF Input
ibooker.setCurrentFolder(monitorDir + "/EMTFInput");
ibooker.setCurrentFolder(monitorDir_ + "/EMTFInput");

ugmtEMTFBX = ibooker.book1D("ugmtEMTFBX", "uGMT EMTF Input BX", 7, -3.5, 3.5);
ugmtEMTFBX->setAxisTitle("BX", 1);
Expand Down Expand Up @@ -240,8 +247,56 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
ugmtEMTFMuMuDR = ibooker.book1D("ugmtEMTFMuMuDR", "uGMT EMTF input muons #DeltaR between sectors", 50, 0., 0.5);
ugmtEMTFMuMuDR->setAxisTitle("#DeltaR", 1);

// EMTF muon showers
if (hadronicShowers_) {
ibooker.setCurrentFolder(monitorDir_ + "/EMTFInput/Muon showers");

ugmtEMTFShowerTypeOccupancyPerSector = ibooker.book2D(
"ugmtEMTFShowerTypeOccupancyPerSector", "Shower type occupancy per sector", 12, 1, 13, 2, 1, 3);
ugmtEMTFShowerTypeOccupancyPerSector->setAxisTitle("Processor", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(12, "+6", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(11, "+5", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(10, "+4", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(9, "+3", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(8, "+2", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(7, "+1", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(6, "-6", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(5, "-5", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(4, "-4", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(3, "-3", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(2, "-2", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(1, "-1", 1);
ugmtEMTFShowerTypeOccupancyPerSector->setAxisTitle("Shower type", 2);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(IDX_TIGHT_SHOWER, "Tight", 2);
ugmtEMTFShowerTypeOccupancyPerSector->setBinLabel(IDX_NOMINAL_SHOWER, "Nominal", 2);

ugmtEMTFShowerTypeOccupancyPerBx =
ibooker.book2D("ugmtEMTFShowerTypeOccupancyPerBx", "Shower type occupancy per BX", 7, -3.5, 3.5, 2, 1, 3);
ugmtEMTFShowerTypeOccupancyPerBx->setAxisTitle("BX", 1);
ugmtEMTFShowerTypeOccupancyPerBx->setAxisTitle("Shower type", 2);
ugmtEMTFShowerTypeOccupancyPerBx->setBinLabel(IDX_TIGHT_SHOWER, "Tight", 2);
ugmtEMTFShowerTypeOccupancyPerBx->setBinLabel(IDX_NOMINAL_SHOWER, "Nominal", 2);

ugmtEMTFShowerSectorOccupancyPerBx = ibooker.book2D(
"ugmtEMTFShowerSectorOccupancyPerBx", "Shower BX occupancy per sector", 7, -3.5, 3.5, 12, 1, 13);
ugmtEMTFShowerSectorOccupancyPerBx->setAxisTitle("BX", 1);
ugmtEMTFShowerSectorOccupancyPerBx->setAxisTitle("Processor", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(12, "+6", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(11, "+5", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(10, "+4", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(9, "+3", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(8, "+2", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(7, "+1", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(6, "-6", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(5, "-5", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(4, "-4", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(3, "-3", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(2, "-2", 2);
ugmtEMTFShowerSectorOccupancyPerBx->setBinLabel(1, "-1", 2);
}

// inter-TF muon correlations
ibooker.setCurrentFolder(monitorDir + "/muon_correlations");
ibooker.setCurrentFolder(monitorDir_ + "/muon_correlations");

ugmtBOMTFposMuMuDEta =
ibooker.book1D("ugmtBOMTFposMuMuDEta", "uGMT input muons #Delta#eta between BMTF and OMTF+", 100, -0.5, 0.5);
Expand Down Expand Up @@ -293,9 +348,9 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
}

// Subsystem Monitoring and Muon Output
ibooker.setCurrentFolder(monitorDir);
ibooker.setCurrentFolder(monitorDir_);

if (!emul) {
if (!emul_) {
ugmtBMTFBXvsProcessor =
ibooker.book2D("ugmtBXvsProcessorBMTF", "uGMT BMTF Input BX vs Processor", 12, -0.5, 11.5, 5, -2.5, 2.5);
ugmtBMTFBXvsProcessor->setAxisTitle("Wedge", 1);
Expand Down Expand Up @@ -523,6 +578,15 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
}
ugmtMuonChargevsLink->setAxisTitle("Charge", 2);

if (hadronicShowers_) {
ugmtMuonShowerTypeOccupancyPerBx =
ibooker.book2D("ugmtMuonShowerTypeOccupancyPerBx", "Shower type occupancy per BX", 7, -3.5, 3.5, 2, 1, 3);
ugmtMuonShowerTypeOccupancyPerBx->setAxisTitle("BX", 1);
ugmtMuonShowerTypeOccupancyPerBx->setAxisTitle("Shower type", 2);
ugmtMuonShowerTypeOccupancyPerBx->setBinLabel(IDX_TIGHT_SHOWER, "Tight", 2);
ugmtMuonShowerTypeOccupancyPerBx->setBinLabel(IDX_NOMINAL_SHOWER, "Nominal", 2);
}

ugmtMuonBXvshwPt = ibooker.book2D("ugmtMuonBXvshwPt", "uGMT Muon BX vs HW p_{T}", 128, -0.5, 511.5, 5, -2.5, 2.5);
ugmtMuonBXvshwPt->setAxisTitle("Hardware p_{T}", 1);
ugmtMuonBXvshwPt->setAxisTitle("BX", 2);
Expand Down Expand Up @@ -554,7 +618,7 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
ugmtMuonBXvshwIso->setAxisTitle("BX", 2);

// muon correlations
ibooker.setCurrentFolder(monitorDir + "/muon_correlations");
ibooker.setCurrentFolder(monitorDir_ + "/muon_correlations");

ugmtMuMuInvMass = ibooker.book1D("ugmtMuMuInvMass", "uGMT dimuon invariant mass", 200, 0., 200.);
ugmtMuMuInvMass->setAxisTitle("m(#mu#mu) [GeV]", 1);
Expand Down Expand Up @@ -680,12 +744,12 @@ void L1TStage2uGMT::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&,
}

void L1TStage2uGMT::analyze(const edm::Event& e, const edm::EventSetup& c) {
if (verbose)
if (verbose_)
edm::LogInfo("L1TStage2uGMT") << "L1TStage2uGMT: analyze..." << std::endl;

if (!emul) {
if (!emul_) {
edm::Handle<l1t::RegionalMuonCandBxCollection> BMTFBxCollection;
e.getByToken(ugmtBMTFToken, BMTFBxCollection);
e.getByToken(ugmtBMTFToken_, BMTFBxCollection);

ugmtBMTFnMuons->Fill(BMTFBxCollection->size(0));

Expand Down Expand Up @@ -734,7 +798,7 @@ void L1TStage2uGMT::analyze(const edm::Event& e, const edm::EventSetup& c) {
}

edm::Handle<l1t::RegionalMuonCandBxCollection> OMTFBxCollection;
e.getByToken(ugmtOMTFToken, OMTFBxCollection);
e.getByToken(ugmtOMTFToken_, OMTFBxCollection);

ugmtOMTFnMuons->Fill(OMTFBxCollection->size(0));

Expand Down Expand Up @@ -789,7 +853,7 @@ void L1TStage2uGMT::analyze(const edm::Event& e, const edm::EventSetup& c) {
}

edm::Handle<l1t::RegionalMuonCandBxCollection> EMTFBxCollection;
e.getByToken(ugmtEMTFToken, EMTFBxCollection);
e.getByToken(ugmtEMTFToken_, EMTFBxCollection);

ugmtEMTFnMuons->Fill(EMTFBxCollection->size(0));

Expand Down Expand Up @@ -847,6 +911,38 @@ void L1TStage2uGMT::analyze(const edm::Event& e, const edm::EventSetup& c) {
}
}

// Fill shower plots
if (hadronicShowers_) {
edm::Handle<l1t::RegionalMuonShowerBxCollection> EMTFShowersBxCollection;
e.getByToken(ugmtEMTFShowerToken_, EMTFShowersBxCollection);

for (int itBX = EMTFShowersBxCollection->getFirstBX(); itBX <= EMTFShowersBxCollection->getLastBX(); ++itBX) {
for (l1t::RegionalMuonShowerBxCollection::const_iterator shower = EMTFShowersBxCollection->begin(itBX);
shower != EMTFShowersBxCollection->end(itBX);
++shower) {
if (not shower->isValid()) {
continue;
}
if (shower->isOneNominalInTime()) {
ugmtEMTFShowerSectorOccupancyPerBx->Fill(
itBX, shower->processor() + 1 + (shower->trackFinderType() == l1t::tftype::emtf_pos ? 6 : 0));
ugmtEMTFShowerTypeOccupancyPerSector->Fill(
shower->processor() + 1 + (shower->trackFinderType() == l1t::tftype::emtf_pos ? 6 : 0),
IDX_NOMINAL_SHOWER);
ugmtEMTFShowerTypeOccupancyPerBx->Fill(itBX, IDX_NOMINAL_SHOWER);
}
if (shower->isOneTightInTime()) {
ugmtEMTFShowerSectorOccupancyPerBx->Fill(
itBX, shower->processor() + 1 + (shower->trackFinderType() == l1t::tftype::emtf_pos ? 6 : 0));
ugmtEMTFShowerTypeOccupancyPerSector->Fill(
shower->processor() + 1 + (shower->trackFinderType() == l1t::tftype::emtf_pos ? 6 : 0),
IDX_TIGHT_SHOWER);
ugmtEMTFShowerTypeOccupancyPerBx->Fill(itBX, IDX_TIGHT_SHOWER);
}
}
}
}

// barrel-overlap muon correlations
int firstBxBO = (BMTFBxCollection->getFirstBX() < OMTFBxCollection->getFirstBX()) ? OMTFBxCollection->getFirstBX()
: BMTFBxCollection->getFirstBX();
Expand Down Expand Up @@ -921,7 +1017,7 @@ void L1TStage2uGMT::analyze(const edm::Event& e, const edm::EventSetup& c) {
}

edm::Handle<l1t::MuonBxCollection> MuonBxCollection;
e.getByToken(ugmtMuonToken, MuonBxCollection);
e.getByToken(ugmtMuonToken_, MuonBxCollection);

ugmtnMuons->Fill(MuonBxCollection->size(0));

Expand Down Expand Up @@ -1062,6 +1158,28 @@ void L1TStage2uGMT::analyze(const edm::Event& e, const edm::EventSetup& c) {
}
}
}

// Fill shower plots
if (hadronicShowers_) {
edm::Handle<l1t::MuonShowerBxCollection> showersBxCollection;
e.getByToken(ugmtMuonShowerToken_, showersBxCollection);

for (int itBX = showersBxCollection->getFirstBX(); itBX <= showersBxCollection->getLastBX(); ++itBX) {
for (l1t::MuonShowerBxCollection::const_iterator shower = showersBxCollection->begin(itBX);
shower != showersBxCollection->end(itBX);
++shower) {
if (not shower->isValid()) {
continue;
}
if (shower->isOneNominalInTime()) {
ugmtMuonShowerTypeOccupancyPerBx->Fill(itBX, IDX_NOMINAL_SHOWER);
}
if (shower->isOneTightInTime()) {
ugmtMuonShowerTypeOccupancyPerBx->Fill(itBX, IDX_TIGHT_SHOWER);
}
}
}
}
}

l1t::tftype L1TStage2uGMT::getTfOrigin(const int tfMuonIndex) {
Expand Down
Loading