HcalNoiseRates Class Reference

#include <HcalNoiseRates.h>

Inheritance diagram for HcalNoiseRates:

Public Member Functions
	HcalNoiseRates (const edm::ParameterSet &)
	~HcalNoiseRates () 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 () 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)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	endJob () 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::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 45 of file HcalNoiseRates.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file HcalNoiseRates.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");
  tok_rbx_ = consumes<reco::HcalNoiseRBXCollection>(rbxCollName_);
  minRBXEnergy_  = iConfig.getUntrackedParameter<double>("minRBXEnergy");
  minHitEnergy_  = iConfig.getUntrackedParameter<double>("minHitEnergy");

  useAllHistos_  = iConfig.getUntrackedParameter<bool>("useAllHistos", false);

  //Hcal Noise Summary
  noisetoken_ = consumes<HcalNoiseSummary>(iConfig.getParameter<edm::InputTag>("noiselabel"));
}

HcalNoiseRates::~HcalNoiseRates ( )

override

Definition at line 100 of file HcalNoiseRates.cc.

{
}

Member Function Documentation

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

overrideprivate

Definition at line 111 of file HcalNoiseRates.cc.

References 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 HcalNoiseRates::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 178 of file HcalNoiseRates.cc.

{}

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 35 of file HcalNoiseRates.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_.

{ 
 
  
  ibooker.setCurrentFolder("HcalNoiseRatesD/HcalNoiseRatesTask");
  
  Char_t histo[100];

  // book histograms

  //Lumi block is not drawn; the rest are
  if (useAllHistos_){
    sprintf  (histo, "hLumiBlockCount" );
    hLumiBlockCount_ = ibooker.book1D(histo, histo, 1, -0.5, 0.5);
  }
  
  sprintf  (histo, "hRBXEnergy" );
  hRBXEnergy_ = ibooker.book1D(histo, histo, 300, 0, 3000);

  sprintf  (histo, "hRBXEnergyType1" );
  hRBXEnergyType1_ = ibooker.book1D(histo, histo, 300, 0, 3000);

  sprintf  (histo, "hRBXEnergyType2" );
  hRBXEnergyType2_ = ibooker.book1D(histo, histo, 300, 0, 3000);

  sprintf  (histo, "hRBXEnergyType3" );
  hRBXEnergyType3_ = ibooker.book1D(histo, histo, 300, 0, 3000);

  sprintf  (histo, "hRBXNHits" );
  hRBXNHits_ = ibooker.book1D(histo, histo, 73,-0.5,72.5);

  //HcalNoiseSummary

  sprintf (histo, "nNNumChannels");
  nNNumChannels_ = ibooker.book1D(histo, histo, 100, 0, 100);
  sprintf (histo, "nNSumE");
  nNSumE_ = ibooker.book1D(histo, histo , 100, 0, 5000);
  sprintf (histo, "nNSumEt");
  nNSumEt_ = ibooker.book1D(histo, histo , 100, 0, 2000);

  sprintf (histo, "sNNumChannels");
  sNNumChannels_ = ibooker.book1D(histo, histo, 100, 0, 100);
  sprintf (histo, "sNSumE");
  sNSumE_ = ibooker.book1D(histo, histo , 100, 0, 5000);
  sprintf (histo, "sNSumEt");
  sNSumEt_ = ibooker.book1D(histo, histo , 100, 0, 2000);

  sprintf (histo, "iNNumChannels");
  iNNumChannels_ = ibooker.book1D(histo, histo, 100, 0, 100);
  sprintf (histo, "iNSumE");
  iNSumE_ = ibooker.book1D(histo, histo , 100, 0, 5000);
  sprintf (histo, "iNSumEt");
  iNSumEt_ = ibooker.book1D(histo, histo , 100, 0, 2000);

  sprintf (histo, "hNoise_maxZeros");
  hNoise_maxZeros_ = ibooker.book1D(histo, histo, 80, 0, 80);
  sprintf (histo, "hNoise_maxHPDHits");
  hNoise_maxHPDHits_ = ibooker.book1D(histo, histo, 20, 0, 20);
  sprintf (histo, "hNoise_maxHPDNoOtherHits");
  hNoise_maxHPDNoOtherHits_ = ibooker.book1D(histo, histo, 20, 0, 20);

}

void HcalNoiseRates::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 182 of file HcalNoiseRates.cc.

References DEFINE_FWK_MODULE, MonitorElement::Fill(), hLumiBlockCount_, lumiCountMap_, and useAllHistos_.

                       {

  if (useAllHistos_) hLumiBlockCount_->Fill(0.0, lumiCountMap_.size()); 

}

Member Data Documentation

MonitorElement* HcalNoiseRates::hLumiBlockCount_

private

Definition at line 69 of file HcalNoiseRates.h.

Referenced by bookHistograms(), and endJob().

MonitorElement* HcalNoiseRates::hNoise_maxHPDHits_

private

Definition at line 94 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hNoise_maxHPDNoOtherHits_

private

Definition at line 95 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hNoise_maxZeros_

private

Definition at line 93 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hRBXEnergy_

private

Definition at line 70 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hRBXEnergyType1_

private

Definition at line 71 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hRBXEnergyType2_

private

Definition at line 72 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hRBXEnergyType3_

private

Definition at line 73 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::hRBXNHits_

private

Definition at line 74 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::iNNumChannels_

private

Definition at line 89 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::iNSumE_

private

Definition at line 90 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::iNSumEt_

private

Definition at line 91 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

std::map<int, int> HcalNoiseRates::lumiCountMap_

private

Definition at line 77 of file HcalNoiseRates.h.

Referenced by analyze(), and endJob().

double HcalNoiseRates::minHitEnergy_

private

Definition at line 63 of file HcalNoiseRates.h.

Referenced by analyze(), and HcalNoiseRates().

double HcalNoiseRates::minRBXEnergy_

private

Definition at line 62 of file HcalNoiseRates.h.

Referenced by analyze(), and HcalNoiseRates().

MonitorElement* HcalNoiseRates::nNNumChannels_

private

Definition at line 81 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::nNSumE_

private

Definition at line 82 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::nNSumEt_

private

Definition at line 83 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<HcalNoiseSummary> HcalNoiseRates::noisetoken_

private

Definition at line 67 of file HcalNoiseRates.h.

Referenced by analyze(), and HcalNoiseRates().

std::string HcalNoiseRates::outputFile_

private

Definition at line 57 of file HcalNoiseRates.h.

Referenced by HcalNoiseRates().

edm::InputTag HcalNoiseRates::rbxCollName_

private

Definition at line 60 of file HcalNoiseRates.h.

Referenced by analyze(), and HcalNoiseRates().

MonitorElement* HcalNoiseRates::sNNumChannels_

private

Definition at line 85 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::sNSumE_

private

Definition at line 86 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalNoiseRates::sNSumEt_

private

Definition at line 87 of file HcalNoiseRates.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 61 of file HcalNoiseRates.h.

Referenced by analyze(), and HcalNoiseRates().

bool HcalNoiseRates::useAllHistos_

private

Definition at line 64 of file HcalNoiseRates.h.

Referenced by bookHistograms(), endJob(), and HcalNoiseRates().