#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 &event, edm::EventSetup const &setup) final

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

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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 DQMEDAnalyzer
uint64_t	meId () const

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 DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 47 of file RPCTTUMonitor.h.

Constructor & Destructor Documentation

◆ RPCTTUMonitor()

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

explicit

Definition at line 6 of file RPCTTUMonitor.cc.

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

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

◆ ~RPCTTUMonitor()

RPCTTUMonitor::~RPCTTUMonitor ( )

override

Definition at line 18 of file RPCTTUMonitor.cc.

18 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 21 of file RPCTTUMonitor.cc.

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

References discriminateDecision(), discriminateGMT(), dqm::impl::MonitorElement::Fill(), L1GtTechnicalTriggerRecord::gtTechnicalTrigger(), iEvent, edm::HandleBase::isValid(), dqmdumpme::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().

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 214 of file RPCTTUMonitor.cc.

                                                                                                          {
   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("");
   }
 }

References dqm::implementation::IBooker::book1D(), dqmdumpme::k, m_bxDistDiffDt, m_bxDistDiffPac, m_dataVsemulator, m_maxttBits, m_ttBits, m_ttBitsDecisionData, m_ttBitsDecisionEmulator, dqm::implementation::NavigatorBase::setCurrentFolder(), and ttuFolder.

◆ discriminateDecision()

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

Definition at line 202 of file RPCTTUMonitor.cc.

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

References data, dqm::impl::MonitorElement::Fill(), and m_dataVsemulator.

Referenced by analyze().

◆ discriminateGMT()

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

Definition at line 124 of file RPCTTUMonitor.cc.

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