#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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

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

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
	)

overrideprotected

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