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

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob ()

virtual void	beginRun (const edm::EventSetup &)

virtual void	endJob ()

Private Attributes
DQMStore *	dbe

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::InputTag	m_gmtReadoutLabel

edm::InputTag	m_gtReadoutLabel

int	m_maxttBits

std::vector< int >	m_RPCcandidatesBx

edm::InputTag	m_rpcDigiLabel

edm::InputTag	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::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

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)

Detailed Description

Definition at line 71 of file RPCTTUMonitor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 7 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     = iConfig.getParameter<edm::InputTag>("GTReadoutRcd");
   m_gmtReadoutLabel    = iConfig.getParameter<edm::InputTag>("GMTReadoutRcd");
   m_rpcTechTrigEmu     = iConfig.getParameter<edm::InputTag>("L1TTEmuBitsLabel");
   m_ttBits             = iConfig.getParameter< std::vector<unsigned> >("BitNumbers");
   m_maxttBits          = m_ttBits.size();
   
 }

RPCTTUMonitor::~RPCTTUMonitor ( )

Definition at line 19 of file RPCTTUMonitor.cc.

 {
   
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 26 of file RPCTTUMonitor.cc.

References edm::EventBase::bunchCrossing(), discriminateDecision(), discriminateGMT(), MonitorElement::Fill(), edm::Event::getByLabel(), 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, and m_ttBitsDecisionEmulator.

 {
   
   //..............................................................................................
   // Data .
   edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
   iEvent.getByLabel( m_gtReadoutLabel, gtRecord);
   
   if ( !gtRecord.isValid() ) {
     edm::LogError("RPCTTUMonitor") << "can't find L1GlobalTriggerRecord with label: " 
                                     << m_gtReadoutLabel << '\n';
     return;
   }
   
   // Emulator .
   edm::Handle< L1GtTechnicalTriggerRecord > emuTTRecord;
   iEvent.getByLabel( m_rpcTechTrigEmu , emuTTRecord);
   
   if ( !emuTTRecord.isValid() ) {
     edm::LogError("RPCTTUMonitor") << "can't find L1GtTechnicalTriggerRecord (emulator) with label: " 
                                     << m_rpcTechTrigEmu << '\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.size() <= 0 ) 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::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 242 of file RPCTTUMonitor.cc.

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

                              {
 
   dbe = edm::Service<DQMStore>().operator->();
   dbe->showDirStructure();
 
   
   dbe->setCurrentFolder(ttuFolder);
 
   
   m_ttBitsDecisionData = dbe->book1D("TechTrigger.Bits.Data",
                                              "Technical Trigger bits : Summary",
                                               10, 23, 33 );
 
   m_ttBitsDecisionEmulator =  dbe->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] =  dbe->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] =  dbe->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] =  dbe->book1D(hname.str().c_str(),
                                               "Comparison between emulator and TT decisions", 
                                               10, 0, 10 );
    
    hname.str("");
    
  }
 }

void RPCTTUMonitor::beginRun ( const edm::EventSetup & iSetup )

privatevirtual

Definition at line 289 of file RPCTTUMonitor.cc.

                                                        {
   
 }

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

Definition at line 221 of file RPCTTUMonitor.cc.

References MonitorElement::Fill(), and m_dataVsemulator.

Referenced by analyze().

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

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

Definition at line 138 of file RPCTTUMonitor.cc.

References edm::Event::getByLabel(), 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.getByLabel(m_gmtReadoutLabel,pCollection);
   
   if ( ! pCollection.isValid() ) {
     edm::LogError("discriminateGMT") << "can't find L1MuGMTReadoutCollection with label "
                                      << m_gmtReadoutLabel ;
     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;
   
 }

void RPCTTUMonitor::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 295 of file RPCTTUMonitor.cc.

References dbe, outputFile, and DQMStore::save().

                       {
 
   if(outputFile != "")
     dbe->save(outputFile);
 
   dbe=0;  
 }

Member Data Documentation

DQMStore* RPCTTUMonitor::dbe

private

Definition at line 93 of file RPCTTUMonitor.h.

Referenced by beginJob(), and endJob().

MonitorElement* RPCTTUMonitor::m_bxDistDiffDt[8]

private

Definition at line 90 of file RPCTTUMonitor.h.

Referenced by analyze(), and beginJob().

MonitorElement* RPCTTUMonitor::m_bxDistDiffPac[8]

private

Definition at line 89 of file RPCTTUMonitor.h.

Referenced by analyze(), and beginJob().

MonitorElement* RPCTTUMonitor::m_dataVsemulator[8]

private

Definition at line 91 of file RPCTTUMonitor.h.

Referenced by beginJob(), and discriminateDecision().

std::vector<int> RPCTTUMonitor::m_DTcandidatesBx

private

Definition at line 104 of file RPCTTUMonitor.h.

Referenced by analyze(), and discriminateGMT().

bool RPCTTUMonitor::m_dtTrigger

private

Definition at line 100 of file RPCTTUMonitor.h.

Referenced by analyze(), and discriminateGMT().

std::vector<int> RPCTTUMonitor::m_GMTcandidatesBx

private

Definition at line 103 of file RPCTTUMonitor.h.

Referenced by analyze(), and discriminateGMT().

edm::InputTag RPCTTUMonitor::m_gmtReadoutLabel

private

Definition at line 109 of file RPCTTUMonitor.h.

Referenced by discriminateGMT(), and RPCTTUMonitor().

edm::InputTag RPCTTUMonitor::m_gtReadoutLabel

private

Definition at line 108 of file RPCTTUMonitor.h.

Referenced by analyze(), and RPCTTUMonitor().

int RPCTTUMonitor::m_maxttBits

private

Definition at line 97 of file RPCTTUMonitor.h.

Referenced by analyze(), beginJob(), and RPCTTUMonitor().

std::vector<int> RPCTTUMonitor::m_RPCcandidatesBx

private

Definition at line 105 of file RPCTTUMonitor.h.

edm::InputTag RPCTTUMonitor::m_rpcDigiLabel

private

Definition at line 107 of file RPCTTUMonitor.h.

edm::InputTag RPCTTUMonitor::m_rpcTechTrigEmu

private

Definition at line 110 of file RPCTTUMonitor.h.

Referenced by analyze(), and RPCTTUMonitor().

bool RPCTTUMonitor::m_rpcTrigger

private

Definition at line 101 of file RPCTTUMonitor.h.

Referenced by analyze(), and discriminateGMT().

std::vector<unsigned> RPCTTUMonitor::m_ttBits

private

Definition at line 98 of file RPCTTUMonitor.h.

Referenced by analyze(), beginJob(), and RPCTTUMonitor().

MonitorElement* RPCTTUMonitor::m_ttBitsDecisionData

private

Definition at line 87 of file RPCTTUMonitor.h.

Referenced by analyze(), and beginJob().

MonitorElement* RPCTTUMonitor::m_ttBitsDecisionEmulator

private

Definition at line 88 of file RPCTTUMonitor.h.

Referenced by analyze(), and beginJob().

std::string RPCTTUMonitor::outputFile

private

Definition at line 95 of file RPCTTUMonitor.h.

Referenced by endJob(), and RPCTTUMonitor().

std::string RPCTTUMonitor::ttuFolder

private

Definition at line 94 of file RPCTTUMonitor.h.

Referenced by beginJob(), and RPCTTUMonitor().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation