L1TMuonCaloSumProducer Class Reference

#include <L1Trigger/L1TGlobalMuon/plugins/L1TMuonCaloSumProducer.cc>

Inheritance diagram for L1TMuonCaloSumProducer:

Public Member Functions
	L1TMuonCaloSumProducer (const edm::ParameterSet &)
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::InputTag	m_caloLabel
edm::EDGetTokenT< CaloTowerBxCollection >	m_caloTowerToken

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: takes generated muons and fills them in the expected collections for the MicroGMT

Implementation: [Notes on implementation]

Definition at line 47 of file L1TMuonCaloSumProducer.cc.

Constructor & Destructor Documentation

L1TMuonCaloSumProducer()

L1TMuonCaloSumProducer::L1TMuonCaloSumProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 71 of file L1TMuonCaloSumProducer.cc.

References edm::ParameterSet::getParameter().

                                                                             {
  //register your inputs:
  m_caloLabel = iConfig.getParameter<edm::InputTag>("caloStage2Layer2Label");
  m_caloTowerToken = consumes<CaloTowerBxCollection>(m_caloLabel);
  //register your products
  produces<MuonCaloSumBxCollection>("TriggerTowerSums");
  produces<MuonCaloSumBxCollection>("TriggerTower2x2s");
}

Member Function Documentation

fillDescriptions()

void L1TMuonCaloSumProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 184 of file L1TMuonCaloSumProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                        {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

produce()

void L1TMuonCaloSumProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 85 of file L1TMuonCaloSumProducer.cc.

References funct::abs(), nano_mu_digi_cff::bx, eleHcalExtractorBlocks_cff::caloTowers, l1t::MuonCaloSum::etBits(), l1trig_cff::hwEta, l1t::L1Candidate::hwEta(), l1trig_cff::hwPhi, l1t::L1Candidate::hwPhi(), l1trig_cff::hwPt, l1t::L1Candidate::hwPt(), hcalRecHitTable_cff::ieta, iEvent, hcalRecHitTable_cff::iphi, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, eostools::move(), and l1t::MuonCaloSum::setEtBits().

                                                                                  {
  using namespace edm;
  std::unique_ptr<MuonCaloSumBxCollection> towerSums(std::make_unique<MuonCaloSumBxCollection>());
  std::unique_ptr<MuonCaloSumBxCollection> tower2x2s(std::make_unique<MuonCaloSumBxCollection>());

  edm::Handle<CaloTowerBxCollection> caloTowers;

  if (iEvent.getByToken(m_caloTowerToken, caloTowers)) {
    int detamax = 4;
    int dphimax = 4;

    const int iFirstBx = caloTowers->getFirstBX();
    const int iLastBx = caloTowers->getLastBX();

    // set BX range for sums
    towerSums->setBXRange(iFirstBx, iLastBx);
    tower2x2s->setBXRange(iFirstBx, iLastBx);

    for (int bx = iFirstBx; bx <= iLastBx; ++bx) {
      std::map<int, MuonCaloSum> sums;
      std::map<int, MuonCaloSum> regs;

      for (auto it = caloTowers->begin(bx); it != caloTowers->end(bx); ++it) {
        const CaloTower& twr = *it;
        int hwEta = twr.hwEta();
        if (std::abs(hwEta) > 27) {
          continue;
        }
        int hwPt = twr.hwPt();
        if (hwPt < 1) {
          continue;
        }
        int hwPhi = twr.hwPhi();

        // calculating tower2x2s
        int ieta2x2 = (hwEta + 27) / 2;
        int iphi2x2 = hwPhi / 2;
        int muon_idx = iphi2x2 * 28 + ieta2x2;
        if (regs.count(muon_idx) == 0) {
          regs[muon_idx] = MuonCaloSum(hwPt, iphi2x2, ieta2x2, muon_idx);
        } else {
          regs.at(muon_idx).setEtBits(regs.at(muon_idx).etBits() + hwPt);
        }

        // std::cout << "iphi; phi " << hwPhi << "; " << phi << " .. ieta; eta" << hwEta << "; " << twr.eta() << std::endl;

        // calculating towerSums
        int ietamax = hwEta + detamax + 1;
        for (int ieta = hwEta - detamax; ieta < ietamax; ++ieta) {
          if (std::abs(ieta) > 27) {
            continue;
          }
          int ietamu = (ieta + 27) / 2;
          int iphimax = hwPhi + dphimax + 1;
          for (int iphi = hwPhi - dphimax; iphi < iphimax; ++iphi) {
            int iphiwrapped = iphi;
            if (iphiwrapped < 0) {
              iphiwrapped += 72;
            } else if (iphiwrapped > 71) {
              iphiwrapped -= 72;
            }
            int iphimu = iphiwrapped / 2;
            int idxmu = iphimu * 28 + ietamu;
            if (sums.count(idxmu) == 0) {
              sums[idxmu] = MuonCaloSum(hwPt, iphimu, ietamu, idxmu);
            } else {
              sums.at(idxmu).setEtBits(sums.at(idxmu).etBits() + hwPt);
            }
          }
        }
      }

      // fill towerSums output collection for this BX
      for (auto it = sums.begin(); it != sums.end(); ++it) {
        if (it->second.etBits() > 0) {
          MuonCaloSum sum = MuonCaloSum(it->second);
          // convert Et to correct scale:
          if (sum.etBits() > 31) {
            sum.setEtBits(31);
          }
          towerSums->push_back(bx, sum);
        }
      }
      // fill tower2x2s output collection for this BX
      for (auto it = regs.begin(); it != regs.end(); ++it) {
        if (it->second.etBits() > 0) {
          tower2x2s->push_back(bx, it->second);
        }
      }
    }
  } else {
    LogWarning("GlobalMuon") << "CaloTowers not found. Producing empty collections." << std::endl;
  }

  iEvent.put(std::move(towerSums), "TriggerTowerSums");
  iEvent.put(std::move(tower2x2s), "TriggerTower2x2s");
}

Member Data Documentation

m_caloLabel

edm::InputTag L1TMuonCaloSumProducer::m_caloLabel

private

Definition at line 57 of file L1TMuonCaloSumProducer.cc.

m_caloTowerToken

edm::EDGetTokenT<CaloTowerBxCollection> L1TMuonCaloSumProducer::m_caloTowerToken

private

Definition at line 56 of file L1TMuonCaloSumProducer.cc.