Implements RPC trigger emulation. More...

#include <RPCTrigger.h>

Inheritance diagram for RPCTrigger:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

	RPCTrigger (const edm::ParameterSet &)

Public Member Functions inherited from edm::one::EDProducer< edm::one::SharedResources >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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 &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::vector< L1MuRegionalCand >	giveFinallCandindates (const L1RpcTBMuonsVec &finalMuons, int type, int bx, edm::Handle< RPCDigiCollection > rpcDigis, std::vector< RPCDigiL1Link > &retRPCDigiLink)
	Returns vector of L1MuRegionalCand (input of L1GMT) More...

Private Attributes
const edm::EDPutTokenT< std::vector< L1MuRegionalCand > >	m_brlCandPutToken

const edm::EDPutTokenT< std::vector< RPCDigiL1Link > >	m_brlLinksPutToken

unsigned long long	m_cacheID

const edm::EDPutTokenT< std::vector< L1MuRegionalCand > >	m_fwdCandPutToken

const edm::EDPutTokenT< std::vector< RPCDigiL1Link > >	m_fwdLinksPutToken

const std::string	m_label

RPCPacManager< RPCPacData >	m_pacManager

std::unique_ptr< RPCPacTrigger >	m_pacTrigger

const edm::EDGetTokenT< RPCDigiCollection >	m_rpcDigiToken

RPCConeBuilderFromES	m_theLinksystemFromES

std::unique_ptr< RPCBasicTrigConfig >	m_trigConfig

const int	m_triggerDebug

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

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

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

Detailed Description

Implements RPC trigger emulation.

Author: Tomasz Fruboes

Definition at line 46 of file RPCTrigger.h.

Constructor & Destructor Documentation

◆ RPCTrigger()

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

explicit

Definition at line 20 of file RPCTrigger.cc.

     : m_trigConfig(),
       m_pacTrigger(),
       // 0 - no debug
       // 1 - dump to xml
       m_triggerDebug{iConfig.getUntrackedParameter<int>("RPCTriggerDebug", 0) == 1 ? 1 : 0},
       m_cacheID{0},
       m_label{iConfig.getParameter<std::string>("label")},
       m_rpcDigiToken{consumes<RPCDigiCollection>(m_label)},
  
       m_brlCandPutToken{produces<std::vector<L1MuRegionalCand> >("RPCb")},
       m_fwdCandPutToken{produces<std::vector<L1MuRegionalCand> >("RPCf")},
  
       m_brlLinksPutToken{produces<std::vector<RPCDigiL1Link> >("RPCb")},
       m_fwdLinksPutToken{produces<std::vector<RPCDigiL1Link> >("RPCf")} {
   //MuonsGrabber is a singleton
   usesResource("MuonsGrabber");
 }

References edm::ParameterSet::getUntrackedParameter().

Member Function Documentation

◆ giveFinallCandindates()

std::vector< L1MuRegionalCand > RPCTrigger::giveFinallCandindates	(	const L1RpcTBMuonsVec &	finalMuons,
		int	type,
		int	bx,
		edm::Handle< RPCDigiCollection >	rpcDigis,
		std::vector< RPCDigiL1Link > &	retRPCDigiLink
	)

private

Returns vector of L1MuRegionalCand (input of L1GMT)

Note: - type is defined in L1MuRegionalCand 1 - barell, 3 - forward

Definition at line 155 of file RPCTrigger.cc.

                                                                                                           {
   std::vector<L1MuRegionalCand> RPCCand;
  
   for (unsigned int iMu = 0; iMu < finalMuons.size(); iMu++) {
     if (finalMuons[iMu].getPtCode() == 0) {
       continue;
     }
  
     RPCDigiL1Link newDigiLink;
  
     //std::cout << "######################################## " << std::endl;
     //std::cout << finalMuons[iMu].getPhiAddr() << " " << finalMuons[iMu].getEtaAddr() << std::endl;
     RPCMuon::TDigiLinkVec digiIVec = finalMuons[iMu].getDigiIdxVec();
     // Here the iteration has to be the same as in
     short int digiIndex = 0;
     RPCDigiCollection::DigiRangeIterator detUnitIt;
     for (detUnitIt = rpcDigis->begin(); detUnitIt != rpcDigis->end(); ++detUnitIt) {
       const RPCDetId& id = (*detUnitIt).first;
       uint32_t rawId = id.rawId();
       const RPCDigiCollection::Range& range = (*detUnitIt).second;
  
       for (RPCDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt) {
         ++digiIndex;
  
         RPCMuon::TDigiLinkVec::iterator it = digiIVec.begin();
         for (; it != digiIVec.end(); ++it) {
           if (digiIndex == it->m_digiIdx) {
             newDigiLink.setLink(it->m_layer + 1, rawId, digiIt->strip(), digiIt->bx());
             //std::cout << type << " " << iMu << " layer: " << it->m_layer <<  " index " << it->m_digiIdx << std::endl;
             //std::cout << "         " << id << " " << " |bx " << digiIt->bx() << " strip " << digiIt->strip() << std::endl;
           }
         }
       }
     }
     retRPCDigiLink.push_back(newDigiLink);
  
     L1MuRegionalCand l1Cand;
  
     l1Cand.setBx(bx);
  
     l1Cand.setQualityPacked(finalMuons[iMu].getQuality());
     l1Cand.setPtPacked(finalMuons[iMu].getPtCode());
  
     l1Cand.setType(type);
  
     int charge = finalMuons[iMu].getSign();
  
     if (charge == 0)  // negative
       l1Cand.setChargePacked(1);
     else
       l1Cand.setChargePacked(0);
  
     //RPCConst::l1RpcConeCrdnts cone = finalMuons[iMu].getConeCrdnts();
  
     /*
     int pac = cone.m_LogSector*12+cone.m_LogSegment;
     const float pi = 3.14159265;
     const float offset = 5*(2*pi/360); // redefinition! Defined also in RPCRingFromRolls::phiMapCompare
     float phi = 2*pi*pac/144-offset;
     if (phi<0)
       phi+=2*pi;
     
     l1Cand.setPhiValue(phi);
     */
  
     //Note: pac numbering begins at 5 deg and goes from 1 to 144.
     // we want phi values from 0 to 2.5 deg to be phiPacked=0
     // max phiPacked value is 143 (see CMS IN 2004-022)
     //int phiPacked = (finalMuons[iMu].getPhiAddr()+2)%144;
     int phiPacked = finalMuons[iMu].getPhiAddr();
     l1Cand.setPhiPacked(phiPacked);
     /*
     float eta = RPCConst::etaFromTowerNum(cone.m_Tower);
     l1Cand.setEtaValue(eta);
 */
     //Note: etaAddr is packed in special way: see CMS IN 2004-022
     signed short etaAddr = finalMuons[iMu].getEtaAddr();  //
     //    signed short etaAddr = finalMuons[iMu].getEtaAddr()-16; // -16..16
     //    bool etaNegative = false;
     //    if (etaAddr < 0){
     //      etaNegative = true;
     //      etaAddr = ~(-etaAddr)+1; // convert to negative :)
     //    }
  
     //    etaAddr &= 63; // 6 bits only
  
     l1Cand.setEtaPacked(etaAddr);
     l1Cand.setChargeValid(true);
  
     /*    
     std::cout<< std::endl << "RBMuon::" << finalMuons[iMu].getEtaAddr() << " " 
              << finalMuons[iMu].getPhiAddr() << std::endl ;
     std::cout<< "cand " <<  l1Cand.eta_packed() << " " 
              << l1Cand.phi_packed() << std::endl ;
     */
  
     RPCCand.push_back(l1Cand);
  
     LogDebug("RPCTrigger") << "Found muonf of pt " << finalMuons[iMu].getPtCode() << " bx " << l1Cand.bx()
                            << " L1Charge " << l1Cand.charge_packed() << " ql " << l1Cand.quality() << " fp "
                            << finalMuons[iMu].getFiredPlanes() << " b/f " << l1Cand.type_idx() << " phi "
                            << l1Cand.phi_packed() << " eta "
                            << l1Cand.eta_packed()
                            //<< " eta l1 " << l1Cand.etaValue() // will drop out soon
                            << " killed " << finalMuons[iMu].wasKilled();
   }
  
   return RPCCand;
 }

Referenced by produce().

◆ produce()

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

overridevirtual

Implements edm::one::EDProducerBase.

Definition at line 39 of file RPCTrigger.cc.

                                                                       {
   //static int ev=1;
   //std::cout << "----------------------------------- "  << ev++ << std::endl;
   if (m_triggerDebug == 1)
     MuonsGrabber::Instance().startNewEvent(iEvent.id().event(), iEvent.bunchCrossing());
  
   if (m_cacheID != iSetup.get<L1RPCConfigRcd>().cacheIdentifier()) {
     //std::cout << " New pats: " << iSetup.get<L1RPCConfigRcd>().cacheIdentifier() << std::endl ;
     m_cacheID = iSetup.get<L1RPCConfigRcd>().cacheIdentifier();
  
     edm::ESHandle<L1RPCConfig> conf;
     iSetup.get<L1RPCConfigRcd>().get(conf);
     const L1RPCConfig* rpcconf = conf.product();
  
     m_pacManager.init(rpcconf);
     m_trigConfig = std::make_unique<RPCBasicTrigConfig>(&m_pacManager);
     m_trigConfig->setDebugLevel(m_triggerDebug);
  
     m_pacTrigger = std::make_unique<RPCPacTrigger>(m_trigConfig.get());
  
     if (m_triggerDebug == 1)
       MuonsGrabber::Instance().setRPCBasicTrigConfig(m_trigConfig.get());
   }
  
   edm::Handle<RPCDigiCollection> rpcDigis;
   iEvent.getByToken(m_rpcDigiToken, rpcDigis);
  
   std::vector<L1MuRegionalCand> candBarell;
   std::vector<L1MuRegionalCand> candForward;
   if (!rpcDigis.isValid()) {
     LogDebug("RPCTrigger")
         << "\nWarning: RPCDigiCollection with input tag " << m_label
         << "\nrequested in configuration, but not found in the event. Emulator will produce empty collection \n ";
  
     iEvent.emplace(m_brlCandPutToken, std::move(candBarell));
     iEvent.emplace(m_fwdCandPutToken, std::move(candForward));
  
     return;
   }
  
   if (rpcDigis->begin() == rpcDigis->end()) {
     LogDebug("RPCTrigger") << "\nWarning: RPCDigiCollection with input tag " << m_label
                            << "\n seems to be empty for this event. Emulator will run on empty collection ";
   }
  
   std::vector<RPCDigiL1Link> brlLinks;
   std::vector<RPCDigiL1Link> fwdLinks;
  
   edm::ESHandle<L1RPCConeBuilder> coneBuilder;
   iSetup.get<L1RPCConeBuilderRcd>().get(coneBuilder);
  
   edm::ESHandle<L1RPCConeDefinition> l1RPCConeDefinition;
   iSetup.get<L1RPCConeDefinitionRcd>().get(l1RPCConeDefinition);
  
   edm::ESHandle<L1RPCHwConfig> hwConfig;
   iSetup.get<L1RPCHwConfigRcd>().get(hwConfig);
  
   edm::ESHandle<L1RPCBxOrConfig> bxOrConfig;
   iSetup.get<L1RPCBxOrConfigRcd>().get(bxOrConfig);
  
   edm::ESHandle<L1RPCHsbConfig> hsbConfig;
   iSetup.get<L1RPCHsbConfigRcd>().get(hsbConfig);
  
   for (int iBx = -1; iBx < 2; ++iBx) {
     L1RpcLogConesVec ActiveCones =
         m_theLinksystemFromES.getConesFromES(rpcDigis, coneBuilder, l1RPCConeDefinition, bxOrConfig, hwConfig, iBx);
  
     L1RpcTBMuonsVec2 finalMuons = m_pacTrigger->runEvent(ActiveCones, hsbConfig);
  
     //int maxFiredPlanes = 0;
  
     /*
     for (unsigned int i=0;i<ActiveCones.size();i++){
         int fpCnt = ActiveCones[i].getFiredPlanesCnt();
         if (fpCnt > maxFiredPlanes)
           maxFiredPlanes = fpCnt;
     }
   
     // Fill out the products
     // finalMuons[0]=barell, finalMuons[1]=endcap
     LogDebug("RPCTrigger") << "---Filling candindates in new event--- " 
                           << maxFiredPlanes << std::endl;
 <<<<<<< RPCTrigger.cc
     */
  
     //    std::vector<L1MuRegionalCand> RPCb = giveFinallCandindates(finalMuons[0],1, iBx);
     //    std::vector<L1MuRegionalCand> RPCf = giveFinallCandindates(finalMuons[1],3, iBx);
     std::vector<RPCDigiL1Link> dlBrl;
     std::vector<RPCDigiL1Link> dlFwd;
     std::vector<L1MuRegionalCand> RPCb = giveFinallCandindates(finalMuons[0], 1, iBx, rpcDigis, dlBrl);
     std::vector<L1MuRegionalCand> RPCf = giveFinallCandindates(finalMuons[1], 3, iBx, rpcDigis, dlFwd);
  
     brlLinks.insert(brlLinks.end(), dlBrl.begin(), dlBrl.end());
     fwdLinks.insert(fwdLinks.end(), dlFwd.begin(), dlFwd.end());
  
     candBarell.insert(candBarell.end(), RPCb.begin(), RPCb.end());
     candForward.insert(candForward.end(), RPCf.begin(), RPCf.end());
  
     if (m_triggerDebug == 1)
       MuonsGrabber::Instance().writeDataForRelativeBX(iBx);
   }
  
   iEvent.emplace(m_fwdLinksPutToken, std::move(fwdLinks));
   iEvent.emplace(m_brlLinksPutToken, std::move(brlLinks));
   iEvent.emplace(m_brlCandPutToken, std::move(candBarell));
   iEvent.emplace(m_fwdCandPutToken, std::move(candForward));
 }

References edm::eventsetup::EventSetupRecord::cacheIdentifier(), muons_cff::finalMuons, edm::EventSetup::get(), get, RPCConeBuilderFromES::getConesFromES(), giveFinallCandindates(), iEvent, RPCPacManager< TPacType >::init(), MuonsGrabber::Instance(), edm::HandleBase::isValid(), LogDebug, m_brlCandPutToken, m_brlLinksPutToken, m_cacheID, m_fwdCandPutToken, m_fwdLinksPutToken, m_label, m_pacManager, m_pacTrigger, m_rpcDigiToken, m_theLinksystemFromES, m_trigConfig, m_triggerDebug, eostools::move(), edm::ESHandle< T >::product(), L1RPCConfig_cff::rpcconf, MuonsGrabber::setRPCBasicTrigConfig(), MuonsGrabber::startNewEvent(), and MuonsGrabber::writeDataForRelativeBX().

Member Data Documentation

◆ m_brlCandPutToken

const edm::EDPutTokenT<std::vector<L1MuRegionalCand> > RPCTrigger::m_brlCandPutToken

private

Definition at line 75 of file RPCTrigger.h.

Referenced by produce().

◆ m_brlLinksPutToken

const edm::EDPutTokenT<std::vector<RPCDigiL1Link> > RPCTrigger::m_brlLinksPutToken

private

Definition at line 78 of file RPCTrigger.h.

Referenced by produce().

◆ m_cacheID

unsigned long long RPCTrigger::m_cacheID

private

Definition at line 64 of file RPCTrigger.h.

Referenced by produce().

◆ m_fwdCandPutToken

const edm::EDPutTokenT<std::vector<L1MuRegionalCand> > RPCTrigger::m_fwdCandPutToken

private

Definition at line 76 of file RPCTrigger.h.

Referenced by produce().

◆ m_fwdLinksPutToken

const edm::EDPutTokenT<std::vector<RPCDigiL1Link> > RPCTrigger::m_fwdLinksPutToken

private

Definition at line 79 of file RPCTrigger.h.

Referenced by produce().

◆ m_label

const std::string RPCTrigger::m_label

private

Definition at line 72 of file RPCTrigger.h.

Referenced by produce().

◆ m_pacManager

RPCPacManager<RPCPacData> RPCTrigger::m_pacManager

private

Definition at line 57 of file RPCTrigger.h.

Referenced by produce().

◆ m_pacTrigger

std::unique_ptr<RPCPacTrigger> RPCTrigger::m_pacTrigger

private

Definition at line 61 of file RPCTrigger.h.

Referenced by produce().

◆ m_rpcDigiToken

const edm::EDGetTokenT<RPCDigiCollection> RPCTrigger::m_rpcDigiToken

private

Definition at line 73 of file RPCTrigger.h.

Referenced by produce().

◆ m_theLinksystemFromES

RPCConeBuilderFromES RPCTrigger::m_theLinksystemFromES

private

Definition at line 55 of file RPCTrigger.h.

Referenced by produce().

◆ m_trigConfig

std::unique_ptr<RPCBasicTrigConfig> RPCTrigger::m_trigConfig

private

Definition at line 59 of file RPCTrigger.h.

Referenced by produce().

◆ m_triggerDebug

const int RPCTrigger::m_triggerDebug

private

Definition at line 63 of file RPCTrigger.h.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RPCTrigger()

Member Function Documentation

◆ giveFinallCandindates()

◆ produce()

Member Data Documentation

◆ m_brlCandPutToken

◆ m_brlLinksPutToken

◆ m_cacheID

◆ m_fwdCandPutToken

◆ m_fwdLinksPutToken

◆ m_label

◆ m_pacManager

◆ m_pacTrigger

◆ m_rpcDigiToken

◆ m_theLinksystemFromES

◆ m_trigConfig

◆ m_triggerDebug