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

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

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

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)

	~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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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

const edm::ESGetToken < L1RPCBxOrConfig, L1RPCBxOrConfigRcd >	m_BxOrConfigToken

unsigned long long	m_cacheID

const edm::ESGetToken < L1RPCConeBuilder, L1RPCConeBuilderRcd >	m_coneBuilderToken

const edm::ESGetToken < L1RPCConeDefinition, L1RPCConeDefinitionRcd >	m_coneDefinitionToken

const edm::ESGetToken < L1RPCConfig, L1RPCConfigRcd >	m_configToken

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

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

const edm::ESGetToken < L1RPCHsbConfig, L1RPCHsbConfigRcd >	m_hsbConfigToken

const edm::ESGetToken < L1RPCHwConfig, L1RPCHwConfigRcd >	m_hwConfigToken

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

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Implements RPC trigger emulation.

Author: Tomasz Fruboes

Definition at line 59 of file RPCTrigger.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file RPCTrigger.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_brlCandPutToken, m_brlLinksPutToken, m_BxOrConfigToken, m_cacheID, m_coneBuilderToken, m_coneDefinitionToken, m_configToken, m_fwdCandPutToken, m_fwdLinksPutToken, m_hsbConfigToken, m_hwConfigToken, m_label, m_rpcDigiToken, and AlCaHLTBitMon_QueryRunRegistry::string.

     : 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")},
       m_configToken(esConsumes<L1RPCConfig, L1RPCConfigRcd>()),
       m_coneBuilderToken(esConsumes<L1RPCConeBuilder, L1RPCConeBuilderRcd>()),
       m_coneDefinitionToken(esConsumes<L1RPCConeDefinition, L1RPCConeDefinitionRcd>()),
       m_hwConfigToken(esConsumes<L1RPCHwConfig, L1RPCHwConfigRcd>()),
       m_BxOrConfigToken(esConsumes<L1RPCBxOrConfig, L1RPCBxOrConfigRcd>()),
       m_hsbConfigToken(esConsumes<L1RPCHsbConfig, L1RPCHsbConfigRcd>()) {
   //MuonsGrabber is a singleton
   usesResource("MuonsGrabber");
 }

Member Function Documentation

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 141 of file RPCTrigger.cc.

Referenced by produce().

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

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

overridevirtual

Implements edm::one::EDProducerBase.

Definition at line 31 of file RPCTrigger.cc.

References edm::EventBase::bunchCrossing(), edm::eventsetup::EventSetupRecord::cacheIdentifier(), edm::Event::emplace(), edm::EventID::event(), muons_cff::finalMuons, edm::EventSetup::get(), edm::Event::getByToken(), RPCConeBuilderFromES::getConesFromES(), edm::EventSetup::getHandle(), giveFinallCandindates(), edm::EventBase::id(), RPCPacManager< TPacType >::init(), MuonsGrabber::Instance(), edm::HandleBase::isValid(), LogDebug, m_brlCandPutToken, m_brlLinksPutToken, m_BxOrConfigToken, m_cacheID, m_coneBuilderToken, m_coneDefinitionToken, m_configToken, m_fwdCandPutToken, m_fwdLinksPutToken, m_hsbConfigToken, m_hwConfigToken, m_label, m_pacManager, m_pacTrigger, m_rpcDigiToken, m_theLinksystemFromES, m_trigConfig, m_triggerDebug, eostools::move(), edm::ESHandle< class >::product(), L1RPCConfig_cff::rpcconf, MuonsGrabber::setRPCBasicTrigConfig(), MuonsGrabber::startNewEvent(), and MuonsGrabber::writeDataForRelativeBX().

                                                                       {
   //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.getHandle(m_configToken);
     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.getHandle(m_coneBuilderToken);
 
   edm::ESHandle<L1RPCConeDefinition> l1RPCConeDefinition = iSetup.getHandle(m_coneDefinitionToken);
 
   edm::ESHandle<L1RPCHwConfig> hwConfig = iSetup.getHandle(m_hwConfigToken);
 
   edm::ESHandle<L1RPCBxOrConfig> bxOrConfig = iSetup.getHandle(m_BxOrConfigToken);
 
   edm::ESHandle<L1RPCHsbConfig> hsbConfig = iSetup.getHandle(m_hsbConfigToken);
 
   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));
 }