#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

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

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

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

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

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

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

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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 () noexcept(false) 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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::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 Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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)

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

edm::EDPutTokenT< DQMToken >	runToken_

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
	)

override

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::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::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::dqmstoreimpl::DQMStore::IBooker::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 DEFINE_FWK_MODULE, 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()));
       }
     }
   }
 }