L1TMicroGMTInputProducerFromGen Class Reference

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

Inheritance diagram for L1TMicroGMTInputProducerFromGen:

Public Member Functions
	L1TMicroGMTInputProducerFromGen (const edm::ParameterSet &)
	~L1TMicroGMTInputProducerFromGen () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Private Member Functions
void	beginLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override
void	beginRun (const edm::Run &, edm::EventSetup const &) override
void	endLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override
void	endRun (const edm::Run &, edm::EventSetup const &) override
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
const static int	m_maxMuons = 108

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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 85 of file L1TMicroGMTInputProducerFromGen.cc.

    : 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");
}

References genParticlesToken, and AlCaHLTBitMon_QueryRunRegistry::string.

~L1TMicroGMTInputProducerFromGen()

L1TMicroGMTInputProducerFromGen::~L1TMicroGMTInputProducerFromGen ( )

override

Definition at line 96 of file L1TMicroGMTInputProducerFromGen.cc.

                                                                  {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

beginLuminosityBlock()

void L1TMicroGMTInputProducerFromGen::beginLuminosityBlock ( const edm::LuminosityBlock &  , edm::EventSetup const &    )

overrideprivate

Definition at line 252 of file L1TMicroGMTInputProducerFromGen.cc.

{}

beginRun()

void L1TMicroGMTInputProducerFromGen::beginRun ( const edm::Run &  , edm::EventSetup const &    )

overrideprivate

Definition at line 246 of file L1TMicroGMTInputProducerFromGen.cc.

{}

compareMuons()

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

staticprivate

Definition at line 105 of file L1TMicroGMTInputProducerFromGen.cc.

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

References l1t::RegionalMuonCand::processor().

Referenced by produce().

endLuminosityBlock()

void L1TMicroGMTInputProducerFromGen::endLuminosityBlock ( const edm::LuminosityBlock &  , edm::EventSetup const &    )

overrideprivate

Definition at line 255 of file L1TMicroGMTInputProducerFromGen.cc.

{}

endRun()

void L1TMicroGMTInputProducerFromGen::endRun ( const edm::Run &  , edm::EventSetup const &    )

overrideprivate

Definition at line 249 of file L1TMicroGMTInputProducerFromGen.cc.

{}

fillDescriptions()

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

static

Definition at line 258 of file L1TMicroGMTInputProducerFromGen.cc.

                                                                                                 {
  //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);
}

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

produce()

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

overrideprivate

Definition at line 110 of file L1TMicroGMTInputProducerFromGen.cc.

                                                                                           {
  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++;
}

References funct::abs(), l1t::bmtf, reco::Candidate::charge(), compareMuons(), l1t::emtf_neg, l1t::emtf_pos, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, reco::Candidate::eta(), genParticles2HepMC_cfi::genParticles, genParticlesToken, mps_fire::i, iEvent, createfilelist::int, LogTrace, m_currEvt, m_maxMuons, m_rnd, MuonErrorMatrixAnalyzer_cfi::maxPt, eostools::move(), amptDefaultParameters_cff::mu, l1t::omtf_neg, l1t::omtf_pos, reco::Candidate::pdgId(), phi, reco::Candidate::phi(), Pi, reco::Candidate::pt(), reco::Candidate::status(), and Geom::twoPi().

Member Data Documentation

genParticlesToken

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

private

Definition at line 68 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by L1TMicroGMTInputProducerFromGen(), and produce().

m_currEvt

int L1TMicroGMTInputProducerFromGen::m_currEvt

private

Definition at line 69 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

m_maxMuons

const static int L1TMicroGMTInputProducerFromGen::m_maxMuons = 108

staticprivate

Definition at line 70 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

m_rnd

TRandom3 L1TMicroGMTInputProducerFromGen::m_rnd

private

Definition at line 71 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().