#include <ESOccupancyTask.h>

Inheritance diagram for ESOccupancyTask:

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

	~ESOccupancyTask () 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 () 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 Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

Private Attributes
int	eCount_

MonitorElement *	hE1E2_ [2]

MonitorElement *	hEnDensity_ [2][2]

MonitorElement *	hEng_ [2][2]

MonitorElement *	hEvEng_ [2][2]

MonitorElement *	hGoodRecNHit_ [2][2]

MonitorElement *	hRecNHit_ [2][2]

MonitorElement *	hRecOCC_ [2][2]

MonitorElement *	hSelEnDensity_ [2][2]

MonitorElement *	hSelEng_ [2][2]

MonitorElement *	hSelOCC_ [2][2]

std::string	prefixME_

edm::EDGetTokenT< ESRecHitCollection >	rechittoken_

int	runNum_

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

Definition at line 14 of file ESOccupancyTask.h.

Constructor & Destructor Documentation

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

Definition at line 25 of file ESOccupancyTask.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and mps_fire::i.

                                                           {
 
   rechittoken_ = consumes<ESRecHitCollection>(ps.getParameter<InputTag>("RecHitLabel"));
   prefixME_ = ps.getUntrackedParameter<string>("prefixME", "EcalPreshower"); 
   
   eCount_ = 0;
   
   //Histogram init  
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       hRecOCC_[i][j] = nullptr;
       hRecNHit_[i][j] = nullptr;
       hEng_[i][j] = nullptr;
       hEvEng_[i][j] = nullptr;
       hEnDensity_[i][i] = nullptr;
       hGoodRecNHit_[i][j] = nullptr;
 
       hSelEng_[i][j] = nullptr;
       hSelOCC_[i][j] = nullptr;
       hSelEnDensity_[i][j] = nullptr;
     }
   
   for (int i = 0; i<2; ++i) 
     hE1E2_[i]=nullptr;
 }

ESOccupancyTask::~ESOccupancyTask ( )

inlineoverride

Definition at line 19 of file ESOccupancyTask.h.

References analyze(), and bookHistograms().

19 {}

Member Function Documentation

void ESOccupancyTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 114 of file ESOccupancyTask.cc.

References edm::SortedCollection< T, SORT >::begin(), DEFINE_FWK_MODULE, edm::SortedCollection< T, SORT >::end(), edm::Event::getByToken(), mps_fire::i, edm::EventBase::id(), edm::EventID::run(), ESDetId::zside(), and ecaldqm::zside().

Referenced by ~ESOccupancyTask().

                                                                             {
 
    runNum_ = e.id().run();
    eCount_++;
 
    // RecHits
    int zside, plane, ix, iy;
    int sum_RecHits[2][2], sum_GoodRecHits[2][2];
    float sum_Energy[2][2];
 
    for (int i = 0; i < 2; ++i) 
      for (int j = 0; j < 2; ++j) {
        sum_RecHits[i][j] = 0;
        sum_GoodRecHits[i][j] = 0;
        sum_Energy[i][j] = 0;
      }
    
    Handle<ESRecHitCollection> ESRecHit;
    if ( e.getByToken(rechittoken_, ESRecHit) ) {
      
      for (ESRecHitCollection::const_iterator hitItr = ESRecHit->begin(); hitItr != ESRecHit->end(); ++hitItr) {
        
        ESDetId id = ESDetId(hitItr->id());
        
        zside = id.zside();
        plane = id.plane();
        ix    = id.six();
        iy    = id.siy();
  
        int i = (zside==1)? 0:1;
        int j = plane-1;
        
        sum_RecHits[i][j]++;
        sum_Energy[i][j] += hitItr->energy();
        hRecOCC_[i][j]->Fill(ix, iy);
        if (hitItr->energy() != 0) {
      hEng_[i][j]->Fill(hitItr->energy());
      hEnDensity_[i][j]->Fill(ix, iy, hitItr->energy());
 
      if (hitItr->recoFlag()==14 || hitItr->recoFlag()==1 || (hitItr->recoFlag()<=10 && hitItr->recoFlag()>=5)) continue;
      sum_GoodRecHits[i][j]++;
      hSelEng_[i][j]->Fill(hitItr->energy());
      hSelEnDensity_[i][j]->Fill(ix, iy, hitItr->energy());
      hSelOCC_[i][j]->Fill(ix, iy);
        }
        
      }
    } else {
      LogWarning("ESOccupancyTask") << "RecHitCollection not available";
    }
 
    //Fill histograms after a event
    for (int i = 0; i < 2; ++i) 
      for (int j = 0; j < 2; ++j) {
 
        hRecNHit_[i][j]->Fill(sum_RecHits[i][j]);
        hGoodRecNHit_[i][j]->Fill(sum_GoodRecHits[i][j]);
        hEvEng_[i][j]->Fill(sum_Energy[i][j]);
        
        //Save eCount_ for Scaling
        hRecOCC_[i][j]->setBinContent(40,40,eCount_);
        hEnDensity_[i][j]->setBinContent(40,40,eCount_);
 
        hSelOCC_[i][j]->setBinContent(40,40,eCount_);
        hSelEnDensity_[i][j]->setBinContent(40,40,eCount_);
      }
 
    hE1E2_[0]->Fill(sum_Energy[0][0], sum_Energy[0][1]);
    hE1E2_[1]->Fill(sum_Energy[1][0], sum_Energy[1][1]);
 
 }

void ESOccupancyTask::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 52 of file ESOccupancyTask.cc.

References DQMStore::IBooker::book1DD(), DQMStore::IBooker::book2D(), trackerHits::histo, mps_fire::i, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

Referenced by ~ESOccupancyTask().

 {
   iBooker.setCurrentFolder(prefixME_ + "/ESOccupancyTask");
   
   //Booking Histograms
   //Notice: Change ESRenderPlugin under DQM/RenderPlugins/src if you change this histogram name.
   char histo[200];
   for (int i=0 ; i<2; ++i) 
     for (int j=0 ; j<2; ++j) {
       int iz = (i==0)? 1:-1;
       sprintf(histo, "ES RecHit 2D Occupancy Z %d P %d", iz, j+1);
       hRecOCC_[i][j] = iBooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       hRecOCC_[i][j]->setAxisTitle("Si X", 1);
       hRecOCC_[i][j]->setAxisTitle("Si Y", 2);
       
       //Bin 40,40 is used to save eumber of event for scaling.
       sprintf(histo, "ES Energy Density Z %d P %d", iz, j+1);
       hEnDensity_[i][j] = iBooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       hEnDensity_[i][j]->setAxisTitle("Si X", 1);
       hEnDensity_[i][j]->setAxisTitle("Si Y", 2);
       
       sprintf(histo, "ES Num of RecHits Z %d P %d", iz, j+1);
       hRecNHit_[i][j] = iBooker.book1DD(histo, histo, 60, 0, 1920);
       hRecNHit_[i][j]->setAxisTitle("# of RecHits", 1);
       hRecNHit_[i][j]->setAxisTitle("Num of Events", 2);
       
       sprintf(histo, "ES Num of Good RecHits Z %d P %d", iz, j+1);
       hGoodRecNHit_[i][j] = iBooker.book1DD(histo, histo, 60, 0, 1920);
       hGoodRecNHit_[i][j]->setAxisTitle("# of good RecHits", 1);
       hGoodRecNHit_[i][j]->setAxisTitle("Num of Events", 2);
       
       sprintf(histo, "ES RecHit Energy Z %d P %d", iz, j+1);
       hEng_[i][j] = iBooker.book1DD(histo, histo, 50, 0, 0.001);
       hEng_[i][j]->setAxisTitle("RecHit Energy", 1);
       hEng_[i][j]->setAxisTitle("Num of ReHits", 2);
       
       sprintf(histo, "ES Event Energy Z %d P %d", iz, j+1);
       hEvEng_[i][j] = iBooker.book1DD(histo, histo, 50, 0, 0.1);
       hEvEng_[i][j]->setAxisTitle("Event Energy", 1);
       hEvEng_[i][j]->setAxisTitle("Num of Events", 2);
 
       // histograms with selected hits
       sprintf(histo, "ES RecHit Energy with selected hits Z %d P %d", iz, j+1);
       hSelEng_[i][j] = iBooker.book1DD(histo, histo, 50, 0, 0.001);
       hSelEng_[i][j]->setAxisTitle("RecHit Energy", 1);
       hSelEng_[i][j]->setAxisTitle("Num of ReHits", 2);
 
       sprintf(histo, "ES Occupancy with selected hits Z %d P %d", iz, j+1);
       hSelOCC_[i][j] = iBooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       hSelOCC_[i][j]->setAxisTitle("Si X", 1);
       hSelOCC_[i][j]->setAxisTitle("Si Y", 2);
 
       sprintf(histo, "ES Energy Density with selected hits Z %d P %d", iz, j+1);
       hSelEnDensity_[i][j] = iBooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       hSelEnDensity_[i][j]->setAxisTitle("Si X", 1);
       hSelEnDensity_[i][j]->setAxisTitle("Si Y", 2);
     }
 
    hE1E2_[0] = iBooker.book2D("ES+ EP1 vs EP2", "ES+ EP1 vs EP2", 50, 0, 0.1, 50, 0, 0.1);
    hE1E2_[1] = iBooker.book2D("ES- EP1 vs EP2", "ES- EP1 vs EP2", 50, 0, 0.1, 50, 0, 0.1);
 }

Member Data Documentation

int ESOccupancyTask::eCount_

private

Definition at line 41 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hE1E2_[2]

private

Definition at line 39 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hEnDensity_[2][2]

private

Definition at line 33 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hEng_[2][2]

private

Definition at line 37 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hEvEng_[2][2]

private

Definition at line 38 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hGoodRecNHit_[2][2]

private

Definition at line 35 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hRecNHit_[2][2]

private

Definition at line 32 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hRecOCC_[2][2]

private

Definition at line 30 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hSelEnDensity_[2][2]

private

Definition at line 34 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hSelEng_[2][2]

private

Definition at line 36 of file ESOccupancyTask.h.

MonitorElement* ESOccupancyTask::hSelOCC_[2][2]

private

Definition at line 31 of file ESOccupancyTask.h.

std::string ESOccupancyTask::prefixME_

private

Definition at line 28 of file ESOccupancyTask.h.

edm::EDGetTokenT<ESRecHitCollection> ESOccupancyTask::rechittoken_

private

Definition at line 27 of file ESOccupancyTask.h.

int ESOccupancyTask::runNum_

private

Definition at line 41 of file ESOccupancyTask.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation