#include <NoiseRates.h>

Inheritance diagram for NoiseRates:

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

	NoiseRates (const edm::ParameterSet &)

	~NoiseRates () override

Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () 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

Private Attributes
MonitorElement *	hLumiBlockCount_

MonitorElement *	hNoise_maxHPDHits_

MonitorElement *	hNoise_maxHPDNoOtherHits_

MonitorElement *	hNoise_maxZeros_

MonitorElement *	hRBXEnergy_

MonitorElement *	hRBXEnergyType1_

MonitorElement *	hRBXEnergyType2_

MonitorElement *	hRBXEnergyType3_

MonitorElement *	hRBXNHits_

MonitorElement *	iNNumChannels_

MonitorElement *	iNSumE_

MonitorElement *	iNSumEt_

std::map< int, int >	lumiCountMap_

double	minHitEnergy_

double	minRBXEnergy_

MonitorElement *	nNNumChannels_

MonitorElement *	nNSumE_

MonitorElement *	nNSumEt_

edm::EDGetTokenT< HcalNoiseSummary >	noisetoken_

std::string	outputFile_

edm::InputTag	rbxCollName_

MonitorElement *	sNNumChannels_

MonitorElement *	sNSumE_

MonitorElement *	sNSumEt_

edm::EDGetTokenT< reco::HcalNoiseRBXCollection >	tok_rbx_

bool	useAllHistos_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
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)

Detailed Description

Definition at line 43 of file NoiseRates.h.

Constructor & Destructor Documentation

NoiseRates::NoiseRates ( const edm::ParameterSet & iConfig )

explicit

Definition at line 17 of file NoiseRates.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), minHitEnergy_, minRBXEnergy_, noisetoken_, outputFile_, rbxCollName_, AlCaHLTBitMon_QueryRunRegistry::string, tok_rbx_, and useAllHistos_.

 {
 
   // DQM ROOT output
   outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile","myfile.root");
 
   // set parameters
   rbxCollName_   = iConfig.getUntrackedParameter<edm::InputTag>("rbxCollName");
   minRBXEnergy_  = iConfig.getUntrackedParameter<double>("minRBXEnergy");
   minHitEnergy_  = iConfig.getUntrackedParameter<double>("minHitEnergy");
 
   tok_rbx_ = consumes<reco::HcalNoiseRBXCollection>(rbxCollName_);
 
   useAllHistos_  = iConfig.getUntrackedParameter<bool>("useAllHistos", false);
 
   //Hcal Noise Summary
   noisetoken_ = consumes<HcalNoiseSummary>(iConfig.getParameter<edm::InputTag>("noiselabel"));
 
 }

NoiseRates::~NoiseRates ( )

override

Definition at line 38 of file NoiseRates.cc.

39 {

40 }

Member Function Documentation

void NoiseRates::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	evSetup
	)

overrideprivate

Definition at line 112 of file NoiseRates.cc.

References DEFINE_FWK_MODULE, Exception, MonitorElement::Fill(), edm::Event::getByToken(), cmsBatch::handle, hNoise_maxHPDHits_, hNoise_maxHPDNoOtherHits_, hNoise_maxZeros_, hRBXEnergy_, hRBXEnergyType1_, hRBXEnergyType2_, hRBXEnergyType3_, hRBXNHits_, iNNumChannels_, iNSumE_, iNSumEt_, HcalNoiseSummary::isolatedNoiseSumE(), HcalNoiseSummary::isolatedNoiseSumEt(), edm::HandleBase::isValid(), lumiCountMap_, edm::EventBase::luminosityBlock(), HcalNoiseSummary::maxHPDHits(), HcalNoiseSummary::maxHPDNoOtherHits(), HcalNoiseSummary::maxZeros(), minHitEnergy_, minRBXEnergy_, HcalNoiseSummary::NegativeNoiseSumE(), HcalNoiseSummary::NegativeNoiseSumEt(), nhits, nNNumChannels_, nNSumE_, nNSumEt_, noisetoken_, HcalNoiseSummary::numIsolatedNoiseChannels(), HcalNoiseSummary::numNegativeNoiseChannels(), reco::HcalNoiseRBX::numRecHits(), HcalNoiseSummary::numSpikeNoiseChannels(), edm::errors::ProductNotFound, rbxCollName_, reco::HcalNoiseRBX::recHitEnergy(), sNNumChannels_, sNSumE_, sNSumEt_, HcalNoiseSummary::spikeNoiseSumE(), HcalNoiseSummary::spikeNoiseSumEt(), heppy_report::summary, and tok_rbx_.

 {
 
   // get the lumi section
   int lumiSection = iEvent.luminosityBlock();
   lumiCountMap_[lumiSection]++;
 
   // get the RBX Noise collection
   edm::Handle<reco::HcalNoiseRBXCollection> handle;
   iEvent.getByToken(tok_rbx_,handle);
   if(!handle.isValid()) {
     throw edm::Exception(edm::errors::ProductNotFound)
       << " could not find HcalNoiseRBXCollection named " << rbxCollName_ << ".\n";
     return;
   }
 
   // get the Noise summary object
   edm::Handle<HcalNoiseSummary> summary_h;
   iEvent.getByToken(noisetoken_, summary_h);
   if(!summary_h.isValid()) {
     throw edm::Exception(edm::errors::ProductNotFound) << " could not find HcalNoiseSummary.\n";
     return;
   }
   const HcalNoiseSummary summary = *summary_h;
 
   //Fill the Noise Summary histograms
   nNNumChannels_->Fill(summary.numNegativeNoiseChannels());
   nNSumE_->Fill(summary.NegativeNoiseSumE());
   nNSumEt_->Fill(summary.NegativeNoiseSumEt());
 
   sNNumChannels_->Fill(summary.numSpikeNoiseChannels());
   sNSumE_->Fill(summary.spikeNoiseSumE());
   sNSumEt_->Fill(summary.spikeNoiseSumEt());
 
   iNNumChannels_->Fill(summary.numIsolatedNoiseChannels());
   iNSumE_->Fill(summary.isolatedNoiseSumE());
   iNSumEt_->Fill(summary.isolatedNoiseSumEt());  
 
   hNoise_maxZeros_->Fill(summary.maxZeros());
   hNoise_maxHPDHits_->Fill(summary.maxHPDHits());
   hNoise_maxHPDNoOtherHits_->Fill(summary.maxHPDNoOtherHits());
 
   // loop over the RBXs and fill the histograms
   for(reco::HcalNoiseRBXCollection::const_iterator it=handle->begin(); it!=handle->end(); ++it) {
     const reco::HcalNoiseRBX &rbx=(*it);
 
     double energy = rbx.recHitEnergy(minHitEnergy_);
 
     int nhits = rbx.numRecHits(minHitEnergy_);
 
     if(energy < minRBXEnergy_) continue;
 
     hRBXEnergy_->Fill(energy);
     
     if      (nhits <= 9)  hRBXEnergyType1_->Fill(energy);
     else if (nhits <= 18) hRBXEnergyType2_->Fill(energy);
     else                  hRBXEnergyType3_->Fill(energy);
     
     hRBXNHits_->Fill(nhits);
     
   }   // done looping over RBXs
 
 }

void NoiseRates::bookHistograms	(	DQMStore::IBooker &	ib,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 47 of file NoiseRates.cc.

References DQMStore::IBooker::book1D(), trackerHits::histo, hLumiBlockCount_, hNoise_maxHPDHits_, hNoise_maxHPDNoOtherHits_, hNoise_maxZeros_, hRBXEnergy_, hRBXEnergyType1_, hRBXEnergyType2_, hRBXEnergyType3_, hRBXNHits_, iNNumChannels_, iNSumE_, iNSumEt_, nNNumChannels_, nNSumE_, nNSumEt_, DQMStore::IBooker::setCurrentFolder(), sNNumChannels_, sNSumE_, sNSumEt_, and useAllHistos_.

 {
 
   ib.setCurrentFolder("NoiseRatesV/NoiseRatesTask");
 
   // book histograms
   Char_t histo[100];
 
   //Lumi block is not drawn; the rest are
   if (useAllHistos_){
     sprintf  (histo, "hLumiBlockCount" );
     hLumiBlockCount_ = ib.book1D(histo, histo, 1, -0.5, 0.5);
   }
   
   sprintf  (histo, "hRBXEnergy" );
   hRBXEnergy_ = ib.book1D(histo, histo, 300, 0, 3000);
 
   sprintf  (histo, "hRBXEnergyType1" );
   hRBXEnergyType1_ = ib.book1D(histo, histo, 300, 0, 3000);
 
   sprintf  (histo, "hRBXEnergyType2" );
   hRBXEnergyType2_ = ib.book1D(histo, histo, 300, 0, 3000);
 
   sprintf  (histo, "hRBXEnergyType3" );
   hRBXEnergyType3_ = ib.book1D(histo, histo, 300, 0, 3000);
 
   sprintf  (histo, "hRBXNHits" );
   hRBXNHits_ = ib.book1D(histo, histo, 73,-0.5,72.5);
 
   //HcalNoiseSummary
 
   sprintf (histo, "nNNumChannels");
   nNNumChannels_ = ib.book1D(histo, histo, 100, 0, 100);
   sprintf (histo, "nNSumE");
   nNSumE_ = ib.book1D(histo, histo , 100, 0, 5000);
   sprintf (histo, "nNSumEt");
   nNSumEt_ = ib.book1D(histo, histo , 100, 0, 2000);
 
   sprintf (histo, "sNNumChannels");
   sNNumChannels_ = ib.book1D(histo, histo, 100, 0, 100);
   sprintf (histo, "sNSumE");
   sNSumE_ = ib.book1D(histo, histo , 100, 0, 5000);
   sprintf (histo, "sNSumEt");
   sNSumEt_ = ib.book1D(histo, histo , 100, 0, 2000);
 
   sprintf (histo, "iNNumChannels");
   iNNumChannels_ = ib.book1D(histo, histo, 100, 0, 100);
   sprintf (histo, "iNSumE");
   iNSumE_ = ib.book1D(histo, histo , 100, 0, 5000);
   sprintf (histo, "iNSumEt");
   iNSumEt_ = ib.book1D(histo, histo , 100, 0, 2000);
 
   sprintf (histo, "hNoise_maxZeros");
   hNoise_maxZeros_ = ib.book1D(histo, histo, 80, 0, 80);
   sprintf (histo, "hNoise_maxHPDHits");
   hNoise_maxHPDHits_ = ib.book1D(histo, histo, 20, 0, 20);
   sprintf (histo, "hNoise_maxHPDNoOtherHits");
   hNoise_maxHPDNoOtherHits_ = ib.book1D(histo, histo, 20, 0, 20);
 
 
 
 }