#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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

override

Definition at line 13 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 78 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()));
       }
     }
   }
 }