#include <L1TPUM.h>

Inheritance diagram for L1TPUM:

Public Member Functions
	L1TPUM (const edm::ParameterSet &ps)

	~L1TPUM () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

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

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

Private Attributes
std::string	histFolder_

MonitorElement *	nonZeroRegionsBx0_

MonitorElement *	nonZeroRegionsBxM2_

MonitorElement *	nonZeroRegionsBxP2_

MonitorElement *	regionBxEtSum_

MonitorElement *	regionBxPopulation_

MonitorElement *	regionsAvgEtBx0_

MonitorElement *	regionsAvgEtBxM2_

MonitorElement *	regionsAvgEtBxP2_

MonitorElement *	regionsAvgNonZeroEtBx0_

MonitorElement *	regionsAvgNonZeroEtBxM2_

MonitorElement *	regionsAvgNonZeroEtBxP2_

edm::EDGetTokenT< L1CaloRegionCollection >	regionSource_

std::vector< MonitorElement * >	regionsPUMEtaBx0_

std::vector< MonitorElement * >	regionsPUMEtaBxM2_

std::vector< MonitorElement * >	regionsPUMEtaBxP2_

MonitorElement *	regionsTotalEtBx0_

MonitorElement *	regionsTotalEtBxM2_

MonitorElement *	regionsTotalEtBxP2_

Detailed Description

Definition at line 25 of file L1TPUM.h.

Constructor & Destructor Documentation

L1TPUM::L1TPUM ( const edm::ParameterSet & ps )

Definition at line 30 of file L1TPUM.cc.

                                          :
    regionSource_(consumes<L1CaloRegionCollection>(ps.getParameter<edm::InputTag>("regionSource"))),
    histFolder_(ps.getParameter<std::string>("histFolder"))
 {
 }

L1TPUM::~L1TPUM ( )

override

Definition at line 36 of file L1TPUM.cc.

37 {

38 }

Member Function Documentation

void L1TPUM::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotected

Definition at line 44 of file L1TPUM.cc.

References conversionPostprocessing_cfi::etaBin, MonitorElement::Fill(), HcalObjRepresent::Fill(), nonZeroRegionsBx0_, nonZeroRegionsBxM2_, nonZeroRegionsBxP2_, regionBxEtSum_, regionBxPopulation_, regionsAvgEtBx0_, regionsAvgEtBxM2_, regionsAvgEtBxP2_, regionsAvgNonZeroEtBx0_, regionsAvgNonZeroEtBxM2_, regionsAvgNonZeroEtBxP2_, regionSource_, regionsPUMEtaBx0_, regionsPUMEtaBxM2_, regionsPUMEtaBxP2_, regionsTotalEtBx0_, regionsTotalEtBxM2_, and regionsTotalEtBxP2_.

 {
   edm::Handle<L1CaloRegionCollection> regionCollection;
   event.getByToken(regionSource_, regionCollection);
 
   int nonzeroRegionsBxP2{0};
   int nonzeroRegionsBx0{0};
   int nonzeroRegionsBxM2{0};
 
   float regionsTotalEtBxP2{0.};
   float regionsTotalEtBx0{0.};
   float regionsTotalEtBxM2{0.};
 
   for (const auto& region : *regionCollection) {
     if ( region.et() > 0 ) {
       if ( region.bx() == 0 ) {
         nonzeroRegionsBx0++;
         regionsTotalEtBx0 += region.et();
       }
       else if ( region.bx() == 2 ) {
         nonzeroRegionsBxP2++;
         regionsTotalEtBxP2 += region.et();
       }
       else if ( region.bx() == -2 ) {
         nonzeroRegionsBxM2++;
         regionsTotalEtBxM2 += region.et();
       }
     }
   }
 
   regionsTotalEtBxP2_->Fill(regionsTotalEtBxP2);
   regionsTotalEtBx0_->Fill(regionsTotalEtBx0);
   regionsTotalEtBxM2_->Fill(regionsTotalEtBxM2);
 
   regionsAvgEtBxP2_->Fill(regionsTotalEtBxP2/396.);
   regionsAvgEtBx0_->Fill(regionsTotalEtBx0/396.);
   regionsAvgEtBxM2_->Fill(regionsTotalEtBxM2/396.);
 
   regionsAvgNonZeroEtBxP2_->Fill(regionsTotalEtBxP2/nonzeroRegionsBxP2);
   regionsAvgNonZeroEtBx0_->Fill(regionsTotalEtBx0/nonzeroRegionsBx0);
   regionsAvgNonZeroEtBxM2_->Fill(regionsTotalEtBxM2/nonzeroRegionsBxM2);
 
   nonZeroRegionsBxP2_->Fill(nonzeroRegionsBxP2);
   nonZeroRegionsBx0_->Fill(nonzeroRegionsBx0);
   nonZeroRegionsBxM2_->Fill(nonzeroRegionsBxM2);
 
   for (const auto& region : *regionCollection) {
     size_t etaBin = region.gctEta();
     regionBxPopulation_->Fill(etaBin*18+region.gctPhi(), region.bx());
     regionBxEtSum_->Fill(etaBin*18+region.gctPhi(), region.bx(), region.et());
     if ( region.bx() == 0 )
       regionsPUMEtaBx0_[etaBin]->Fill(nonzeroRegionsBx0/PUMNORMALIZE, region.et());
     else if ( region.bx() == 2 )
       regionsPUMEtaBxP2_[etaBin]->Fill(nonzeroRegionsBxP2/PUMNORMALIZE, region.et());
     else if ( region.bx() == -2 )
       regionsPUMEtaBxM2_[etaBin]->Fill(nonzeroRegionsBxM2/PUMNORMALIZE, region.et());
   }
 }

void L1TPUM::bookHistograms	(	DQMStore::IBooker &	ibooker,
		const edm::Run &	run,
		const edm::EventSetup &	es
	)

overrideprotected

Definition at line 103 of file L1TPUM.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), histFolder_, nonZeroRegionsBx0_, nonZeroRegionsBxM2_, nonZeroRegionsBxP2_, R10BINS, R10MAX, R10MIN, regionBxEtSum_, regionBxPopulation_, regionsAvgEtBx0_, regionsAvgEtBxM2_, regionsAvgEtBxP2_, regionsAvgNonZeroEtBx0_, regionsAvgNonZeroEtBxM2_, regionsAvgNonZeroEtBxP2_, regionsPUMEtaBx0_, regionsPUMEtaBxM2_, regionsPUMEtaBxP2_, regionsTotalEtBx0_, regionsTotalEtBxM2_, regionsTotalEtBxP2_, and DQMStore::IBooker::setCurrentFolder().

 {
   ibooker.setCurrentFolder(histFolder_+"/BX0");
   regionsPUMEtaBx0_.resize(PUMETABINS);
   for (size_t ieta=0; ieta<PUMETABINS; ++ieta) {
     regionsTotalEtBx0_ = ibooker.book1D("regionsTotalEt", "Total ET distribution;Sum Rank;Counts", 200, 0, 2000);
     regionsAvgEtBx0_ = ibooker.book1D("regionsAvgEt", "Average Rank;Average Rank;Counts", R10BINS, R10MIN, R10MAX);
     regionsAvgNonZeroEtBx0_ = ibooker.book1D("regionsAvgNonZeroEt", "Average Rank >0;Average Rank Regions>0;Counts", R10BINS, R10MIN, R10MAX);
     nonZeroRegionsBx0_ = ibooker.book1D("nonZeroRegions", "Nonzero regions;Number Regions >0;Counts", 397, -0.5, 396.5);
     regionsPUMEtaBx0_[ieta] = ibooker.book2D("regionsPUMEta"+std::to_string(ieta), "PUM Bin rank distribution;PU bin;Rank", PUMBINS, PUMMIN, PUMMAX, R10BINS, R10MIN, R10MAX);
   }
 
   ibooker.setCurrentFolder(histFolder_+"/BXP2");
   regionsPUMEtaBxP2_.resize(PUMETABINS);
   for (size_t ieta=0; ieta<PUMETABINS; ++ieta) {
     regionsTotalEtBxP2_ = ibooker.book1D("regionsTotalEt", "Total ET distribution;Sum Rank;Counts", 200, 0, 2000);
     regionsAvgEtBxP2_ = ibooker.book1D("regionsAvgEt", "Average Rank;Average Rank;Counts", R10BINS, R10MIN, R10MAX);
     regionsAvgNonZeroEtBxP2_ = ibooker.book1D("regionsAvgNonZeroEt", "Average Rank >0;Average Rank Regions>0;Counts", R10BINS, R10MIN, R10MAX);
     nonZeroRegionsBxP2_ = ibooker.book1D("nonZeroRegions", "Nonzero regions;Number Regions >0;Counts", 397, -0.5, 396.5);
     regionsPUMEtaBxP2_[ieta] = ibooker.book2D("regionsPUMEta"+std::to_string(ieta), "PUM Bin rank distribution;PU bin;Rank", PUMBINS, PUMMIN, PUMMAX, R10BINS, R10MIN, R10MAX);
   }
 
   ibooker.setCurrentFolder(histFolder_+"/BXM2");
   regionsPUMEtaBxM2_.resize(PUMETABINS);
   for (size_t ieta=0; ieta<PUMETABINS; ++ieta) {
     regionsTotalEtBxM2_ = ibooker.book1D("regionsTotalEt", "Total ET distribution;Sum Rank;Counts", 200, 0, 2000);
     regionsAvgEtBxM2_ = ibooker.book1D("regionsAvgEt", "Average Rank;Average Rank;Counts", R10BINS, R10MIN, R10MAX);
     regionsAvgNonZeroEtBxM2_ = ibooker.book1D("regionsAvgNonZeroEt", "Average Rank >0;Average Rank Regions>0;Counts", R10BINS, R10MIN, R10MAX);
     nonZeroRegionsBxM2_ = ibooker.book1D("nonZeroRegions", "Nonzero regions;Number Regions >0;Counts", 397, -0.5, 396.5);
     regionsPUMEtaBxM2_[ieta] = ibooker.book2D("regionsPUMEta"+std::to_string(ieta), "PUM Bin rank distribution;PU bin;Rank", PUMBINS, PUMMIN, PUMMAX, R10BINS, R10MIN, R10MAX);
   }
 
   ibooker.setCurrentFolder(histFolder_+"/RegionBxInfo");
   regionBxPopulation_ = ibooker.book2D("regionBxPopulation", "Event counts per region per bunch crossing;Region index (18*eta+phi);BX index;Counts", 396, -0.5, 395.5, 5, -2.5, 2.5);
   regionBxEtSum_ = ibooker.book2D("regionBxEtSum", "Et per region per bunch crossing;Region index (18*eta+phi);BX index;Counts*et", 396, -0.5, 395.5, 5, -2.5, 2.5);
 }