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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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 submitPVResolutionJobs::desc.

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().