#include <RPCTTUMonitor.h>

Inheritance diagram for RPCTTUMonitor:

Public Member Functions
void	discriminateDecision (bool, bool, int)

int	discriminateGMT (const edm::Event &iEvent, const edm::EventSetup &iSetup)

	RPCTTUMonitor (const edm::ParameterSet &)

	~RPCTTUMonitor () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	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 () 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

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)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

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)

Private Attributes
MonitorElement *	m_bxDistDiffDt [8]

MonitorElement *	m_bxDistDiffPac [8]

MonitorElement *	m_dataVsemulator [8]

std::vector< int >	m_DTcandidatesBx

bool	m_dtTrigger

std::vector< int >	m_GMTcandidatesBx

edm::EDGetTokenT< L1MuGMTReadoutCollection >	m_gmtReadoutLabel

edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	m_gtReadoutLabel

int	m_maxttBits

std::vector< int >	m_RPCcandidatesBx

edm::EDGetTokenT< L1GtTechnicalTriggerRecord >	m_rpcTechTrigEmu

bool	m_rpcTrigger

std::vector< unsigned >	m_ttBits

MonitorElement *	m_ttBitsDecisionData

MonitorElement *	m_ttBitsDecisionEmulator

std::string	outputFile

std::string	ttuFolder

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 48 of file RPCTTUMonitor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file RPCTTUMonitor.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_gmtReadoutLabel, m_gtReadoutLabel, m_maxttBits, m_rpcTechTrigEmu, m_ttBits, outputFile, AlCaHLTBitMon_QueryRunRegistry::string, and ttuFolder.

                                                           {
 
   ttuFolder    = iConfig.getUntrackedParameter<std::string>("TTUFolder", "RPC/TTU");
   outputFile    = iConfig.getUntrackedParameter<std::string>("OutPutFile", ""); 
 
   m_gtReadoutLabel     = consumes<L1GlobalTriggerReadoutRecord>(iConfig.getParameter<edm::InputTag>("GTReadoutRcd"));
   m_gmtReadoutLabel    = consumes<L1MuGMTReadoutCollection>(iConfig.getParameter<edm::InputTag>("GMTReadoutRcd"));
   m_rpcTechTrigEmu     = consumes<L1GtTechnicalTriggerRecord>(iConfig.getParameter<edm::InputTag>("L1TTEmuBitsLabel"));
   
 
   m_ttBits             = iConfig.getParameter< std::vector<unsigned> >("BitNumbers");
   m_maxttBits          = m_ttBits.size();
   
 }

RPCTTUMonitor::~RPCTTUMonitor ( )

override

Definition at line 21 of file RPCTTUMonitor.cc.

21 {}

Member Function Documentation

void RPCTTUMonitor::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprotected

Definition at line 25 of file RPCTTUMonitor.cc.

References edm::EventBase::bunchCrossing(), discriminateDecision(), discriminateGMT(), MonitorElement::Fill(), edm::Event::getByToken(), L1GtTechnicalTriggerRecord::gtTechnicalTrigger(), edm::HandleBase::isValid(), gen::k, m_bxDistDiffDt, m_bxDistDiffPac, m_DTcandidatesBx, m_dtTrigger, m_GMTcandidatesBx, m_gtReadoutLabel, m_maxttBits, m_rpcTechTrigEmu, m_rpcTrigger, m_ttBits, m_ttBitsDecisionData, m_ttBitsDecisionEmulator, and L1GlobalTriggerReadoutRecord::technicalTriggerWord().

 {
   
   //..............................................................................................
   // Data .
   edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
   iEvent.getByToken( m_gtReadoutLabel, gtRecord);
   
   if ( !gtRecord.isValid() ) {
     edm::LogError("RPCTTUMonitor") << "can nout find L1GlobalTriggerRecord \n";
     return;
   }
   
   // Emulator .
   edm::Handle< L1GtTechnicalTriggerRecord > emuTTRecord;
   iEvent.getByToken( m_rpcTechTrigEmu , emuTTRecord);
   
   if ( !emuTTRecord.isValid() ) {
     edm::LogError("RPCTTUMonitor") << "can not find L1GtTechnicalTriggerRecord (emulator) \n";
     return;
   }
   
   //..............................................................................................
   //
   //Timing difference between RPC-PAT and DT
   
   int dGMT(0);
   dGMT = discriminateGMT( iEvent , iSetup );
   if ( dGMT < 0 ) return;
   
   std::map<int,bool> ttuDec;
   std::map<int,bool>::iterator decItr;
   
   int bxX = iEvent.bunchCrossing(); // ... 1 to 3564
   
   for( int k=0; k < m_maxttBits; ++k) {
     for( int iebx=0; iebx<=2; iebx++) {
       const TechnicalTriggerWord gtTTWord = gtRecord->technicalTriggerWord(iebx-1);
       ttuDec[iebx-1] = gtTTWord[ 24+k ];
     }
   
     //. RPC
     if ( m_rpcTrigger ) {
       
       int ndec(0);
       int bx1 = (bxX - m_GMTcandidatesBx[0]);
       for( decItr = ttuDec.begin(); decItr != ttuDec.end(); ++decItr ){
         if ( (*decItr).second ) {
           int bx2 = (*decItr).first;
           float bxdiffPacTT = 1.0*( bx1 - bx2);
           m_bxDistDiffPac[k]->Fill( bxdiffPacTT );
           ++ndec;
         }
       }
     }
     
     //.. DT
     if ( m_dtTrigger ) {
       
       int ndec(0);
       int bx1 = (bxX - m_DTcandidatesBx[0]);
       for( decItr = ttuDec.begin(); decItr != ttuDec.end(); ++decItr ){
         if ( (*decItr).second ) {
           int bx2 = (*decItr).first;
           float bxdiffDtTT = 1.0*( bx1 - bx2);
           m_bxDistDiffDt[k]->Fill( bxdiffDtTT );
           ++ndec;
         }
       }
     }
     ttuDec.clear();
   
   }
 
   m_GMTcandidatesBx.clear();
   m_DTcandidatesBx.clear();
     
   //..............................................................................................
   //
   //... For Data Emulator comparison
 
   const TechnicalTriggerWord gtTTWord = gtRecord->technicalTriggerWord();
 
   std::vector<L1GtTechnicalTrigger> ttVec = emuTTRecord->gtTechnicalTrigger();
 
   std::vector<unsigned>::iterator bitsItr;
   int k = 0;
   //int m_BxWindow = 0;
   bool hasDataTrigger = false;
   bool hasEmulatorTrigger = false;
 
   if ( ttVec.empty() ) return;
   
   for ( bitsItr = m_ttBits.begin(); bitsItr != m_ttBits.end(); ++bitsItr ) {
     
     hasDataTrigger = gtTTWord.at( (*bitsItr) );
     m_ttBitsDecisionData->Fill( (*bitsItr), (int)hasDataTrigger  );
     
     hasEmulatorTrigger = ttVec[k].gtTechnicalTriggerResult();
     m_ttBitsDecisionEmulator->Fill( ttVec[k].gtTechnicalTriggerBitNumber(), (int)hasEmulatorTrigger );
     
     discriminateDecision(hasDataTrigger ,hasEmulatorTrigger  , k );
         
     ++k;
     
   }
   
 
 }

void RPCTTUMonitor::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	r,
		edm::EventSetup const &	iSetup
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 235 of file RPCTTUMonitor.cc.

References DQMStore::IBooker::book1D(), gen::k, m_bxDistDiffDt, m_bxDistDiffPac, m_dataVsemulator, m_maxttBits, m_ttBits, m_ttBitsDecisionData, m_ttBitsDecisionEmulator, DQMStore::IBooker::setCurrentFolder(), and ttuFolder.

                                                                                                              {
 
   ibooker.setCurrentFolder(ttuFolder);
 
   
   m_ttBitsDecisionData = ibooker.book1D("TechTrigger.Bits.Data",
                                              "Technical Trigger bits : Summary",
                                               10, 23, 33 );
 
   m_ttBitsDecisionEmulator =  ibooker.book1D("TechTrigger.Bits.Emulator",
                                                  "Technical Trigger bits : Summary",
                                                  10, 23, 33 );
  for( int k=0; k < m_maxttBits; ++k) {
     
    std::ostringstream hname;
    
    hname << "BX.diff.PAC-TTU.bit." << m_ttBits[k];
    
    m_bxDistDiffPac[k] =  ibooker.book1D(hname.str().c_str(),
                                              "Timing difference between PAC and TTU",
                                              7, -3, 3);
    
    hname.str("");
    
    hname << "BX.diff.DT-TTU.bit." << m_ttBits[k];
    
    m_bxDistDiffDt[k] =  ibooker.book1D(hname.str().c_str(),
                                             "Timing difference between DT and TTU",
                                             7, -3, 3);
    
    hname.str("");
    
    hname << "Emu.Ttu.Compare.bit." << m_ttBits[k];
    
    m_dataVsemulator[k] =  ibooker.book1D(hname.str().c_str(),
                                               "Comparison between emulator and TT decisions", 
                                               10, 0, 10 );
    
    hname.str("");
    
  }
 
  
 }

void RPCTTUMonitor::discriminateDecision	(	bool	data,
		bool	emu,
		int	indx
	)

Definition at line 218 of file RPCTTUMonitor.cc.

References MonitorElement::Fill(), and m_dataVsemulator.

Referenced by analyze().

                                                                          {
  
   if ( data == 1 && emu == 1 ){
     m_dataVsemulator[indx]->Fill( 1 );
   }else if ( data == 1 && emu == 0 ){
     m_dataVsemulator[indx]->Fill( 3 );
   }else if ( data == 0 && emu == 1 ){
     m_dataVsemulator[indx]->Fill( 5 );
   }else if ( data == 0 && emu == 0 ){
     m_dataVsemulator[indx]->Fill( 7 );
   }
   
 }

int RPCTTUMonitor::discriminateGMT	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 135 of file RPCTTUMonitor.cc.

References edm::Event::getByToken(), L1MuGMTReadoutCollection::getRecords(), edm::HandleBase::isValid(), m_DTcandidatesBx, m_dtTrigger, m_GMTcandidatesBx, m_gmtReadoutLabel, m_rpcTrigger, and edm::Handle< T >::product().

Referenced by analyze().

 {
   
   //.............................................................................................
   
   edm::Handle<L1MuGMTReadoutCollection> pCollection;
   iEvent.getByToken(m_gmtReadoutLabel,pCollection);
   
   if ( ! pCollection.isValid() ) {
     edm::LogError("discriminateGMT") << "can't find L1MuGMTReadoutCollection with label \n";
 
     return -1; 
   }
   
   //.............................................................................................
   
   int gmtDec(0);
   
   bool rpcBar_l1a  = false;
   bool dtBar_l1a   = false;
   
   m_dtTrigger  = false;
   m_rpcTrigger = false;
   
   // get GMT readout collection
   const L1MuGMTReadoutCollection * gmtRC = pCollection.product();
   
   // get record vector
   std::vector<L1MuGMTReadoutRecord>::const_iterator RRItr;
   std::vector<L1MuGMTReadoutRecord> gmt_records = gmtRC->getRecords();
   
   edm::LogInfo("DiscriminateGMT") << "nRecords: " << gmt_records.size() << '\n';
   
   for( RRItr = gmt_records.begin(); RRItr != gmt_records.end(); ++RRItr ) {
     
     int BxInEvent = RRItr->getBxInEvent();
     int BxInEventNew = RRItr->getBxNr();
     
     // RPC barrel muon candidates
     int nrpcB = 0;
     int ndtB  = 0;
     
     std::vector<L1MuRegionalCand> BrlRpcCands = RRItr->getBrlRPCCands();
     std::vector<L1MuRegionalCand> BrlDtCands  = RRItr->getDTBXCands ();
     
     std::vector<L1MuRegionalCand>::const_iterator RCItr;
     
     for( RCItr = BrlRpcCands.begin(); RCItr !=BrlRpcCands.end(); ++RCItr) {
       if ( !(*RCItr).empty() ) {
         
         m_GMTcandidatesBx.push_back( BxInEventNew );
         
         nrpcB++;
       }
     }
     
     for( RCItr = BrlDtCands.begin(); RCItr !=BrlDtCands.end(); ++RCItr) {
       if ( !(*RCItr).empty() ) {
         m_DTcandidatesBx.push_back( BxInEventNew );
         ndtB++;
       }
     }
     
     if( BxInEvent == 0 && nrpcB > 0) rpcBar_l1a = true;
     if( BxInEvent == 0 && ndtB > 0) dtBar_l1a = true;
     
   }
   
   if( rpcBar_l1a ) { 
     gmtDec = 1;
     m_rpcTrigger = true;
     
   }
   
   if( dtBar_l1a ) { 
     gmtDec = 2;
     m_dtTrigger = true;
   }
   
   return gmtDec;
   
 }