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

Inheritance diagram for L1TMicroGMTInputProducerFromGen:

Public Member Functions
	L1TMicroGMTInputProducerFromGen (const edm::ParameterSet &)

	~L1TMicroGMTInputProducerFromGen ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	beginJob () override

virtual void	beginLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override

virtual void	beginRun (const edm::Run &, edm::EventSetup const &) override

virtual void	endJob () override

virtual void	endLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override

virtual void	endRun (const edm::Run &, edm::EventSetup const &) override

virtual 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::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

Definition at line 102 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::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 259 of file L1TMicroGMTInputProducerFromGen.cc.

260 {

261 }

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

overrideprivatevirtual

Definition at line 282 of file L1TMicroGMTInputProducerFromGen.cc.

283 {

284 }

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

overrideprivatevirtual

Definition at line 270 of file L1TMicroGMTInputProducerFromGen.cc.

271 {

272 }

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

staticprivate

Definition at line 114 of file L1TMicroGMTInputProducerFromGen.cc.

References l1t::RegionalMuonCand::processor().

Referenced by produce().

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

void L1TMicroGMTInputProducerFromGen::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 265 of file L1TMicroGMTInputProducerFromGen.cc.

265 {

266 }

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

overrideprivatevirtual

Definition at line 288 of file L1TMicroGMTInputProducerFromGen.cc.

289 {

290 }

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

overrideprivatevirtual

Definition at line 276 of file L1TMicroGMTInputProducerFromGen.cc.

277 {

278 }

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

static

Definition at line 294 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
	)

overrideprivatevirtual

Definition at line 121 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++;
 
 }