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

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

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

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

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

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

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

Implements RPC trigger emulation.

Author: Tomasz Fruboes

Definition at line 50 of file RPCTrigger.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file RPCTrigger.cc.

References edm::ParameterSet::getUntrackedParameter(), m_brlCandPutToken, m_brlLinksPutToken, m_cacheID, m_fwdCandPutToken, m_fwdLinksPutToken, 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")}
 {
   //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 172 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
	)

override

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

Member Data Documentation

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

private

Definition at line 76 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

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

private

Definition at line 79 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

unsigned long long RPCTrigger::m_cacheID

private

Definition at line 68 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

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

private

Definition at line 77 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

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

private

Definition at line 80 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

const std::string RPCTrigger::m_label

private

Definition at line 73 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

RPCPacManager<RPCPacData> RPCTrigger::m_pacManager

private

Definition at line 61 of file RPCTrigger.h.

Referenced by produce().

std::unique_ptr<RPCPacTrigger> RPCTrigger::m_pacTrigger

private

Definition at line 65 of file RPCTrigger.h.

Referenced by produce().

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

private

Definition at line 74 of file RPCTrigger.h.

Referenced by produce(), and RPCTrigger().

RPCConeBuilderFromES RPCTrigger::m_theLinksystemFromES

private

Definition at line 59 of file RPCTrigger.h.

Referenced by produce().

std::unique_ptr<RPCBasicTrigConfig> RPCTrigger::m_trigConfig

private

Definition at line 63 of file RPCTrigger.h.

Referenced by produce().

const int RPCTrigger::m_triggerDebug

private

Definition at line 67 of file RPCTrigger.h.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation