#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

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_

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
	)

override

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