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

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

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

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

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

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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

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)

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 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::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)

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 74 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 12 of file RCTMonitor.cc.

12 {}

Member Function Documentation

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

override

Definition at line 83 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 14 of file RCTMonitor.cc.

References DQMStore::IBooker::book1D(), 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 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 103 of file RCTMonitor.cc.

References DEFINE_FWK_MODULE, DynamicScale(), 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()));
       }
     }
   }
 }