#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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
static const 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 51 of file L1TMicroGMTInputProducerFromGen.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 87 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 100 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 269 of file L1TMicroGMTInputProducerFromGen.cc.

270 {

271 }

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

overrideprivate

Definition at line 257 of file L1TMicroGMTInputProducerFromGen.cc.

258 {

259 }

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

staticprivate

Definition at line 112 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 275 of file L1TMicroGMTInputProducerFromGen.cc.

276 {

277 }

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

overrideprivate

Definition at line 263 of file L1TMicroGMTInputProducerFromGen.cc.

264 {

265 }

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

static

Definition at line 281 of file L1TMicroGMTInputProducerFromGen.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, 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 119 of file L1TMicroGMTInputProducerFromGen.cc.

References funct::abs(), l1t::bmtf, reco::Candidate::charge(), compareMuons(), l1t::emtf_neg, l1t::emtf_pos, PVValHelper::eta, reco::Candidate::eta(), GenHFHadronMatcher_cfi::genParticles, genParticlesToken, edm::Event::getByToken(), mps_fire::i, createfilelist::int, LogTrace, m_currEvt, m_maxMuons, m_rnd, MuonErrorMatrixAnalyzer_cfi::maxPt, eostools::move(), RPCpg::mu, l1t::omtf_neg, l1t::omtf_pos, reco::Candidate::pdgId(), phi, reco::Candidate::phi(), Pi, reco::Candidate::pt(), edm::Event::put(), l1t::RegionalMuonCand::setHwEta(), l1t::RegionalMuonCand::setHwPhi(), l1t::RegionalMuonCand::setHwPt(), l1t::RegionalMuonCand::setHwQual(), l1t::RegionalMuonCand::setHwSign(), l1t::RegionalMuonCand::setHwSignValid(), l1t::RegionalMuonCand::setTFIdentifiers(), 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++;
 
 }