#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

	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	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
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 ( const edm::ParameterSet & iConfig )

explicit

Definition at line 85 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");
 }

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

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

overrideprivate

Definition at line 252 of file L1TMicroGMTInputProducerFromGen.cc.

252 {}

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

overrideprivate

Definition at line 246 of file L1TMicroGMTInputProducerFromGen.cc.

246 {}

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

staticprivate

Definition at line 105 of file L1TMicroGMTInputProducerFromGen.cc.

References l1t::RegionalMuonCand::processor().

Referenced by produce().

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

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

overrideprivate

Definition at line 255 of file L1TMicroGMTInputProducerFromGen.cc.

255 {}

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

overrideprivate

Definition at line 249 of file L1TMicroGMTInputProducerFromGen.cc.

249 {}

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

static

Definition at line 258 of file L1TMicroGMTInputProducerFromGen.cc.

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

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

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

overrideprivate

Definition at line 110 of file L1TMicroGMTInputProducerFromGen.cc.

References funct::abs(), l1t::bmtf, reco::Candidate::charge(), compareMuons(), l1t::emtf_neg, l1t::emtf_pos, relval_parameters_module::energy, PVValHelper::eta, reco::Candidate::eta(), genParticleCandidates2GenParticles_cfi::genParticles, genParticlesToken, edm::Event::getByToken(), mps_fire::i, LogTrace, m_currEvt, m_maxMuons, m_rnd, eostools::move(), RPCpg::mu, l1t::omtf_neg, l1t::omtf_pos, reco::Candidate::pdgId(), phi, reco::Candidate::phi(), Pi, reco::Candidate::pt(), edm::Event::put(), 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++;
 }