L1TMicroGMTInputProducerFromGen Class Reference

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

Inheritance diagram for L1TMicroGMTInputProducerFromGen:

Public Member Functions
	L1TMicroGMTInputProducerFromGen (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

Static Private Member Functions
static bool	compareMuons (const RegionalMuonCand &, const RegionalMuonCand &)

Private Attributes
edm::EDGetTokenT< reco::GenParticleCollection >	genParticlesToken
int	m_currEvt
TRandom3	m_rnd

Static Private Attributes
static const int	m_maxMuons = 108

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 50 of file L1TMicroGMTInputProducerFromGen.cc.

Constructor & Destructor Documentation

L1TMicroGMTInputProducerFromGen()

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

explicit

Definition at line 79 of file L1TMicroGMTInputProducerFromGen.cc.

References genParticlesToken, and AlCaHLTBitMon_QueryRunRegistry::string.

    : m_currEvt(0), m_rnd(0) {
  //register your inputs:
  genParticlesToken = consumes<reco::GenParticleCollection>(std::string("genParticles"));
  //register your products
  produces<RegionalMuonCandBxCollection>("BarrelTFMuons");
  produces<RegionalMuonCandBxCollection>("OverlapTFMuons");
  produces<RegionalMuonCandBxCollection>("ForwardTFMuons");
  produces<MuonCaloSumBxCollection>("TriggerTowerSums");
}

Member Function Documentation

compareMuons()

bool L1TMicroGMTInputProducerFromGen::compareMuons ( const RegionalMuonCand &  mu1, const RegionalMuonCand &  mu2  )

staticprivate

Definition at line 94 of file L1TMicroGMTInputProducerFromGen.cc.

References l1t::RegionalMuonCand::processor().

Referenced by produce().

                                                                                                           {
  return mu1.processor() < mu2.processor();
}

fillDescriptions()

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

static

Definition at line 235 of file L1TMicroGMTInputProducerFromGen.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 L1TMicroGMTInputProducerFromGen::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 99 of file L1TMicroGMTInputProducerFromGen.cc.

References funct::abs(), l1t::bmtf, reco::Candidate::charge(), compareMuons(), l1t::emtf_neg, l1t::emtf_pos, HBHEDarkening_cff::energy, PVValHelper::eta, reco::Candidate::eta(), AJJGenJetFilter_cfi::genParticles, genParticlesToken, l1trig_cff::hwCharge, l1trig_cff::hwChargeValid, l1trig_cff::hwEta, l1trig_cff::hwPhi, l1trig_cff::hwPt, l1trig_cff::hwQual, mps_fire::i, iEvent, createfilelist::int, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogTrace, m_currEvt, m_maxMuons, m_rnd, PV_cfg::maxPt, eostools::move(), amptDefaultParameters_cff::mu, l1t::omtf_neg, l1t::omtf_pos, reco::Candidate::pdgId(), phi, reco::Candidate::phi(), Pi, reco::Candidate::pt(), jetUpdater_cfi::sort, reco::Candidate::status(), and Geom::twoPi().

                                                                                           {
  using namespace edm;

  std::unique_ptr<RegionalMuonCandBxCollection> barrelMuons(new RegionalMuonCandBxCollection());
  std::unique_ptr<RegionalMuonCandBxCollection> overlapMuons(new RegionalMuonCandBxCollection());
  std::unique_ptr<RegionalMuonCandBxCollection> endcapMuons(new RegionalMuonCandBxCollection());
  std::unique_ptr<MuonCaloSumBxCollection> towerSums(new MuonCaloSumBxCollection());

  std::vector<RegionalMuonCand> bmMuons;
  std::vector<RegionalMuonCand> omMuons;
  std::vector<RegionalMuonCand> emMuons;

  std::vector<int> muIndices;
  edm::Handle<reco::GenParticleCollection> genParticles;
  // Make sure that you can get genParticles
  if (iEvent.getByToken(genParticlesToken, genParticles)) {
    int cntr = 0;
    for (auto it = genParticles->cbegin(); it != genParticles->cend(); ++it) {
      const reco::Candidate& mcParticle = *it;
      if (abs(mcParticle.pdgId()) == 13 && mcParticle.status() == 1)
        muIndices.push_back(cntr);
      cntr++;
    }
  } else {
    LogTrace("GlobalMuon") << " GenParticleCollection not found." << std::endl;
  }

  RegionalMuonCand mu;
  MuonCaloSum tSum;
  // alternative scale (using full phi bit-width): 163.4521265553765f;
  const float phiToInt = 91.67324722093171f;
  // alternative scale: 100.0f;
  const float etaToInt = 90.9090909090f;
  const int maxPt = (1 << 9) - 1;
  int muCntr = 0;

  double twoPi = TMath::Pi() * 2.;

  for (auto it = muIndices.begin(); it != muIndices.end(); ++it) {
    // don't really care which muons are taken...
    // guess there ain't 108 generated anyways
    if (muCntr == m_maxMuons)
      break;
    int gen_idx = *it;
    const reco::Candidate& mcMuon = genParticles->at(gen_idx);
    double eta = mcMuon.eta();
    if (fabs(eta) > 2.45)
      continue;  // out of acceptance
    int hwPt = int(mcMuon.pt() * 2);
    hwPt = (hwPt < maxPt ? hwPt : maxPt);
    int hwEta = int(eta * etaToInt);
    double phi = mcMuon.phi();
    if (phi < 0)
      phi += twoPi;  // add 2*pi
    int hwPhi = (int(phi * phiToInt)) % 576;
    int hwQual = 8;
    int hwCharge = (mcMuon.charge() > 0) ? 0 : 1;
    int hwChargeValid = 1;

    mu.setHwPt(hwPt);

    tftype tf(tftype::bmtf);
    int globalWedgePhi = (hwPhi + 24) % 576;  // this sets CMS phi = 0 to -15 deg
    int localPhi = globalWedgePhi % 48;
    int processor = globalWedgePhi / 48 + 1;
    int globalSectorPhi = (hwPhi - 24);  // this sets CMS phi = 0 to +15 deg
    if (globalSectorPhi < 0) {
      globalSectorPhi += 576;
    }

    if (fabs(eta) > 0.8) {
      if (fabs(eta) < 1.2) {
        tf = (eta > 0 ? tftype::omtf_pos : tftype::omtf_neg);
        processor = globalSectorPhi / 96 + 1;
        localPhi = globalSectorPhi % 96;
      } else {
        tf = (eta > 0 ? tftype::emtf_pos : tftype::emtf_neg);
        processor = globalSectorPhi / 96 + 1;
        localPhi = globalSectorPhi % 96;
      }
    }
    mu.setHwPhi(localPhi);
    mu.setTFIdentifiers(processor, tf);

    mu.setHwEta(hwEta);
    mu.setHwSign(hwCharge);
    mu.setHwSignValid(hwChargeValid);
    mu.setHwQual(hwQual);

    if (fabs(eta) < 0.8 && bmMuons.size() < 36) {
      bmMuons.push_back(mu);
      muCntr++;
    } else if (fabs(eta) < 1.2 && omMuons.size() < 36) {
      omMuons.push_back(mu);
      muCntr++;
    } else if (emMuons.size() < 36) {
      emMuons.push_back(mu);
      muCntr++;
    }
  }

  std::sort(bmMuons.begin(), bmMuons.end(), L1TMicroGMTInputProducerFromGen::compareMuons);
  std::sort(omMuons.begin(), omMuons.end(), L1TMicroGMTInputProducerFromGen::compareMuons);
  std::sort(emMuons.begin(), emMuons.end(), L1TMicroGMTInputProducerFromGen::compareMuons);

  for (const auto& mu : bmMuons) {
    barrelMuons->push_back(0, mu);
  }

  for (const auto& mu : omMuons) {
    overlapMuons->push_back(0, mu);
  }

  for (const auto& mu : emMuons) {
    endcapMuons->push_back(0, mu);
  }

  for (int i = 0; i < 1008; ++i) {
    // from where could I take the tower energies?
    int energy = int(m_rnd.Gaus(12, 6));
    if (energy < 0)
      energy = 0;
    if (energy > 31)
      energy = 31;
    MuonCaloSum sum(energy, i / 28, i % 28, i);
    towerSums->push_back(0, sum);
  }

  iEvent.put(std::move(barrelMuons), "BarrelTFMuons");
  iEvent.put(std::move(overlapMuons), "OverlapTFMuons");
  iEvent.put(std::move(endcapMuons), "ForwardTFMuons");
  iEvent.put(std::move(towerSums), "TriggerTowerSums");
  m_currEvt++;
}

Member Data Documentation

genParticlesToken

edm::EDGetTokenT<reco::GenParticleCollection> L1TMicroGMTInputProducerFromGen::genParticlesToken

private

Definition at line 62 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by L1TMicroGMTInputProducerFromGen(), and produce().

m_currEvt

int L1TMicroGMTInputProducerFromGen::m_currEvt

private

Definition at line 63 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

m_maxMuons

const int L1TMicroGMTInputProducerFromGen::m_maxMuons = 108

staticprivate

Definition at line 64 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

m_rnd

TRandom3 L1TMicroGMTInputProducerFromGen::m_rnd

private

Definition at line 65 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().