#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=default

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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_

const double	minHitEnergy_

const double	minRBXEnergy_

MonitorElement *	nNNumChannels_

MonitorElement *	nNSumE_

MonitorElement *	nNSumEt_

const edm::EDGetTokenT < HcalNoiseSummary >	noisetoken_

const std::string	outputFile_

const edm::InputTag	rbxCollName_

MonitorElement *	sNNumChannels_

MonitorElement *	sNSumE_

MonitorElement *	sNSumEt_

const edm::EDGetTokenT < reco::HcalNoiseRBXCollection >	tok_rbx_

const bool	useAllHistos_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 41 of file NoiseRates.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file NoiseRates.cc.

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

NoiseRates::~NoiseRates ( )

overridedefault

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 93 of file NoiseRates.cc.

References relval_parameters_module::energy, Exception, dqm::impl::MonitorElement::Fill(), edm::Event::getHandle(), patZpeak::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(), edmLumisInFiles::summary, and tok_rbx_.

                                                                              {
   // get the lumi section
   int lumiSection = iEvent.luminosityBlock();
   lumiCountMap_[lumiSection]++;
 
   // get the RBX Noise collection
   const edm::Handle<reco::HcalNoiseRBXCollection> &handle = iEvent.getHandle(tok_rbx_);
   if (!handle.isValid()) {
     throw edm::Exception(edm::errors::ProductNotFound)
         << " could not find HcalNoiseRBXCollection named " << rbxCollName_ << ".\n";
     return;
   }
 
   // get the Noise summary object
   const edm::Handle<HcalNoiseSummary> &summary_h = iEvent.getHandle(noisetoken_);
   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 34 of file NoiseRates.cc.

References dqm::implementation::IBooker::book1D(), timingPdfMaker::histo, hLumiBlockCount_, hNoise_maxHPDHits_, hNoise_maxHPDNoOtherHits_, hNoise_maxZeros_, hRBXEnergy_, hRBXEnergyType1_, hRBXEnergyType2_, hRBXEnergyType3_, hRBXNHits_, iNNumChannels_, iNSumE_, iNSumEt_, nNNumChannels_, nNSumE_, nNSumEt_, dqm::implementation::NavigatorBase::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);
 }