#include <RCTMonitor.h>

Inheritance diagram for RCTMonitor:

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

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

void	FillRCT (const edm::Event &, const edm::EventSetup &)

	RCTMonitor (const edm::ParameterSet &)

	~RCTMonitor () 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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
MonitorElement *	m_rctIsoEmOccEtaPhi1

MonitorElement *	m_rctIsoEmOccEtaPhi10

MonitorElement *	m_rctIsoEmRank1

MonitorElement *	m_rctIsoEmRank10

MonitorElement *	m_rctIsoEmRankEtaPhi1

MonitorElement *	m_rctIsoEmRankEtaPhi10

MonitorElement *	m_rctNonIsoEmOccEtaPhi1

MonitorElement *	m_rctNonIsoEmOccEtaPhi10

MonitorElement *	m_rctNonIsoEmRank1

MonitorElement *	m_rctNonIsoEmRank10

MonitorElement *	m_rctNonIsoEmRankEtaPhi1

MonitorElement *	m_rctNonIsoEmRankEtaPhi10

MonitorElement *	m_rctRegionEt

MonitorElement *	m_rctRegionsEtEtaPhi

MonitorElement *	m_rctRegionsOccEtaPhi

MonitorElement *	m_rctRelaxedEmOccEtaPhi1

MonitorElement *	m_rctRelaxedEmOccEtaPhi10

MonitorElement *	m_rctRelaxedEmRank1

MonitorElement *	m_rctRelaxedEmRank10

MonitorElement *	m_rctRelaxedEmRankEtaPhi1

MonitorElement *	m_rctRelaxedEmRankEtaPhi10

edm::EDGetTokenT < L1CaloEmCollection >	m_rctSourceToken_

MonitorElement *	m_rctTauVetoEtaPhi

Static Private Attributes
static const unsigned int	ETABINS = 22

static const float	ETAMAX = 21.5

static const float	ETAMIN = -0.5

static const unsigned int	L1EETABINS = 22

static const float	L1EETAMAX = 5

static const float	L1EETAMIN = -5

static const unsigned int	L1EPHIBINS = 18

static const float	L1EPHIMAX = M_PI

static const float	L1EPHIMIN = -M_PI

static const unsigned int	METPHIBINS = 72

static const float	METPHIMAX = 71.5

static const float	METPHIMIN = -0.5

static const unsigned int	PHIBINS = 18

static const float	PHIMAX = 17.5

static const float	PHIMIN = -0.5

static const unsigned int	R10BINS = 1024

static const float	R10MAX = 1023.5

static const float	R10MIN = -0.5

static const unsigned int	R12BINS = 4096

static const float	R12MAX = 4095.5

static const float	R12MIN = -0.5

static const unsigned int	R6BINS = 64

static const float	R6MAX = 63.5

static const float	R6MIN = -0.5

static const unsigned int	RTPBINS = 101

static const float	RTPMAX = 100.5

static const float	RTPMIN = -0.5

static const unsigned int	TEVBINS = 1001

static const float	TEVMAX = 1000.5

static const float	TEVMIN = -0.5

static const unsigned int	TPETABINS = 65

static const float	TPETAMAX = 32.5

static const float	TPETAMIN = -32.5

static const unsigned int	TPPHIBINS = 72

static const float	TPPHIMAX = 72.5

static const float	TPPHIMIN = 0.5

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 >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Description: DQM monitor for the Regional Calorimeter Trigger

Definition at line 73 of file RCTMonitor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file RCTMonitor.cc.

References edm::ParameterSet::getUntrackedParameter(), and m_rctSourceToken_.

                                                      {
   // set Token(-s)
   m_rctSourceToken_ = consumes<L1CaloEmCollection>(iConfig.getUntrackedParameter<edm::InputTag>("rctSource"));
 }

RCTMonitor::~RCTMonitor ( )

override

Definition at line 11 of file RCTMonitor.cc.

11 {}

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 59 of file RCTMonitor.cc.

References FillRCT().

                                                                             {
   // Fill histograms
   FillRCT(iEvent, iSetup);
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 13 of file RCTMonitor.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), ETABINS, ETAMAX, ETAMIN, m_rctIsoEmOccEtaPhi1, m_rctIsoEmOccEtaPhi10, m_rctIsoEmRank1, m_rctIsoEmRank10, m_rctIsoEmRankEtaPhi1, m_rctIsoEmRankEtaPhi10, m_rctNonIsoEmOccEtaPhi1, m_rctNonIsoEmOccEtaPhi10, m_rctNonIsoEmRank1, m_rctNonIsoEmRank10, m_rctNonIsoEmRankEtaPhi1, m_rctNonIsoEmRankEtaPhi10, m_rctRegionEt, m_rctRegionsEtEtaPhi, m_rctRegionsOccEtaPhi, m_rctRelaxedEmOccEtaPhi1, m_rctRelaxedEmOccEtaPhi10, m_rctRelaxedEmRank1, m_rctRelaxedEmRank10, m_rctRelaxedEmRankEtaPhi1, m_rctRelaxedEmRankEtaPhi10, m_rctTauVetoEtaPhi, PHIBINS, PHIMAX, PHIMIN, R10BINS, R10MAX, R10MIN, R6BINS, R6MAX, R6MIN, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                {
   // Book RCT histograms
   iBooker.setCurrentFolder("RCT");
 
   m_rctIsoEmRankEtaPhi1 =
       iBooker.book2D("RctIsoEmRankEtaPhi", "ISO EM RANK", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctIsoEmOccEtaPhi1 =
       iBooker.book2D("RctIsoEmOccEtaPhi", "ISO EM OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctIsoEmRank1 = iBooker.book1D("RctIsoEmRank", "ISO EM RANK", R6BINS, R6MIN, R6MAX);
   m_rctIsoEmRankEtaPhi10 =
       iBooker.book2D("RctIsoEmRankEtaPhi10", "ISO EM RANK", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctIsoEmOccEtaPhi10 =
       iBooker.book2D("RctIsoEmOccEtaPhi10", "ISO EM OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctIsoEmRank10 = iBooker.book1D("RctIsoEmRank10", "ISO EM RANK", R6BINS, R6MIN, R6MAX);
 
   m_rctNonIsoEmRankEtaPhi1 =
       iBooker.book2D("RctNonIsoEmRankEtaPhi", "NON-ISO EM RANK", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctNonIsoEmOccEtaPhi1 =
       iBooker.book2D("RctNonIsoEmOccEtaPhi", "NON-ISO EM OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctNonIsoEmRank1 = iBooker.book1D("RctNonIsoEmRank", "NON-ISO EM RANK", R6BINS, R6MIN, R6MAX);
   m_rctNonIsoEmRankEtaPhi10 =
       iBooker.book2D("RctNonIsoEmRankEtaPhi10", "NON-ISO EM RANK", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctNonIsoEmOccEtaPhi10 = iBooker.book2D(
       "RctNonIsoEmOccEtaPhi10", "NON-ISO EM OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctNonIsoEmRank10 = iBooker.book1D("RctNonIsoEmRank10", "NON-ISO EM RANK", R6BINS, R6MIN, R6MAX);
 
   m_rctRelaxedEmRankEtaPhi1 =
       iBooker.book2D("RctRelaxedEmRankEtaPhi", "RELAXED EM RANK", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctRelaxedEmOccEtaPhi1 =
       iBooker.book2D("RctRelaxedEmOccEtaPhi", "RELAXED EM OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctRelaxedEmRank1 = iBooker.book1D("RctRelaxedEmRank", "RELAXED EM RANK", R6BINS, R6MIN, R6MAX);
   m_rctRelaxedEmRankEtaPhi10 =
       iBooker.book2D("RctRelaxedEmRankEtaPhi", "RELAXED EM RANK", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctRelaxedEmOccEtaPhi10 = iBooker.book2D(
       "RctRelaxedEmOccEtaPhi10", "RELAXED EM OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctRelaxedEmRank10 = iBooker.book1D("RctRelaxedEmRank", "RELAXED EM RANK", R6BINS, R6MIN, R6MAX);
 
   m_rctRegionsEtEtaPhi =
       iBooker.book2D("RctRegionsEtEtaPhi", "REGION E_{T}", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctRegionsOccEtaPhi =
       iBooker.book2D("RctRegionsOccEtaPhi", "REGION OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctTauVetoEtaPhi =
       iBooker.book2D("RctTauVetoEtaPhi", "TAU VETO OCCUPANCY", PHIBINS, PHIMIN, PHIMAX, ETABINS, ETAMIN, ETAMAX);
   m_rctRegionEt = iBooker.book1D("RctRegionEt", "REGION E_{T}", R10BINS, R10MIN, R10MAX);
 }

void RCTMonitor::FillRCT	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 78 of file RCTMonitor.cc.

References DynamicScale(), dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), m_rctIsoEmOccEtaPhi1, m_rctIsoEmOccEtaPhi10, m_rctIsoEmRank1, m_rctIsoEmRankEtaPhi1, m_rctNonIsoEmOccEtaPhi1, m_rctNonIsoEmOccEtaPhi10, m_rctNonIsoEmRank1, m_rctNonIsoEmRankEtaPhi1, m_rctRelaxedEmOccEtaPhi1, m_rctRelaxedEmOccEtaPhi10, m_rctRelaxedEmRank1, m_rctRelaxedEmRankEtaPhi1, and m_rctSourceToken_.

Referenced by analyze().

                                                                             {
   // Get the RCT digis
   edm::Handle<L1CaloEmCollection> em;
 
   iEvent.getByToken(m_rctSourceToken_, em);
 
   // Isolated and non-isolated EM with cut at >1 GeV
   for (L1CaloEmCollection::const_iterator iem = em->begin(); iem != em->end(); iem++) {
     if (iem->rank() > 1.) {   // applies the 1 GeV cut
       if (iem->isolated()) {  // looks for isolated EM candidates only
         m_rctIsoEmRank1->Fill(iem->rank());
         // std::cout << "Just to show what is there " << iem->rank() <<
         // std::endl ;
         m_rctIsoEmRankEtaPhi1->Fill(iem->regionId().iphi(), iem->regionId().ieta(), iem->rank());
         m_rctIsoEmOccEtaPhi1->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
         m_rctRelaxedEmRankEtaPhi1->Fill(iem->regionId().iphi(), iem->regionId().ieta(), iem->rank());
         m_rctRelaxedEmOccEtaPhi1->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
         m_rctRelaxedEmRank1->Fill(iem->rank());
       } else {  // instructions for Non-isolated EM candidates
         m_rctNonIsoEmRank1->Fill(iem->rank());
         m_rctNonIsoEmRankEtaPhi1->Fill(iem->regionId().iphi(), iem->regionId().ieta(), iem->rank());
         m_rctNonIsoEmOccEtaPhi1->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
         m_rctRelaxedEmRankEtaPhi1->Fill(iem->regionId().iphi(), iem->regionId().ieta(), iem->rank());
         m_rctRelaxedEmOccEtaPhi1->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
         m_rctRelaxedEmRank1->Fill(iem->rank());
       }
     }
     if (iem->rank() > 10.) {  // applies the 10 GeV cut
       if (iem->isolated()) {  // looks for isolated EM candidates only
         m_rctIsoEmOccEtaPhi10->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
         m_rctRelaxedEmOccEtaPhi10->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
       } else {  // instructions for Non-isolated EM candidates
         m_rctNonIsoEmOccEtaPhi10->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
         m_rctRelaxedEmOccEtaPhi10->Fill(
             iem->regionId().iphi(), iem->regionId().ieta(), DynamicScale(iem->regionId().ieta()));
       }
     }
   }
 }