#include <L1TPUM.h>

Inheritance diagram for L1TPUM:

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

	~L1TPUM () 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	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

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)

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

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)

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_

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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);
 }