#include <DQMOffline/Hcal/src/HcalCollapseAnalyzer.cc>

Inheritance diagram for HcalCollapseAnalyzer:

Public Member Functions
	HcalCollapseAnalyzer (const edm::ParameterSet &)

	~HcalCollapseAnalyzer () override

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

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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

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

Private Attributes
const bool	doHB_

const bool	doHE_

MonitorElement *	h_balance

MonitorElement *	h_depth

MonitorElement *	h_frac

MonitorElement *	h_merge

MonitorElement *	h_sfrac

MonitorElement *	h_size

const edm::InputTag	preRecHitHBHE_

const edm::InputTag	recHitHBHE_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< HBHERecHitCollection >	tok_prehbhe_

const std::string	topFolderName_

const int	verbosity_

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

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

Description: Studies the collapsing of HcalRecHits

Implementation: <Notes on="" implementation>="">

Definition at line 47 of file HcalCollapseAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 72 of file HcalCollapseAnalyzer.cc.

References doHB_, doHE_, preRecHitHBHE_, recHitHBHE_, tok_hbhe_, tok_prehbhe_, and verbosity_.

                                                                          : 
   topFolderName_(iConfig.getParameter<std::string>("topFolderName")),
   verbosity_(iConfig.getUntrackedParameter<int>("verbosity",0)),
   recHitHBHE_(iConfig.getUntrackedParameter<edm::InputTag>("recHitHBHE",edm::InputTag("hbhereco"))),
   preRecHitHBHE_(iConfig.getUntrackedParameter<edm::InputTag>("preRecHitHBHE",edm::InputTag("hbheprereco"))),
   doHE_(iConfig.getUntrackedParameter<bool>("doHE", true)),
   doHB_(iConfig.getUntrackedParameter<bool>("doHB", false)) {
 
   // define tokens for access
   tok_hbhe_    = consumes<HBHERecHitCollection>(recHitHBHE_);
   tok_prehbhe_ = consumes<HBHERecHitCollection>(preRecHitHBHE_);
   
   edm::LogVerbatim("Collapse") << "Verbosity " << verbosity_ << " with tags " 
                    << recHitHBHE_  << " and " << preRecHitHBHE_
                    << " and Do " << doHB_ << ":" << doHE_;
 }

HcalCollapseAnalyzer::~HcalCollapseAnalyzer ( )

inlineoverride

Definition at line 51 of file HcalCollapseAnalyzer.cc.

References analyze(), bookHistograms(), and fillDescriptions().

51 { }

Member Function Documentation

void HcalCollapseAnalyzer::analyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 100 of file HcalCollapseAnalyzer.cc.

References egammaForCoreTracking_cff::depth, doHB_, doHE_, edm::EventID::event(), MonitorElement::Fill(), spr::find(), cropTnPTrees::frac, edm::EventSetup::get(), edm::Event::getByToken(), h_balance, h_depth, h_frac, h_merge, h_sfrac, h_size, hcalLocalReco_cff::hbheprereco, HiIsolationCommonParameters_cff::hbhereco, HcalBarrel, HcalEndcap, edm::EventBase::id(), hit::id, photons_cff::ids, edm::HandleBase::isValid(), gen::k, edm::Handle< T >::product(), edm::ESHandle< T >::product(), edm::EventID::run(), edm::SortedCollection< T, SORT >::size(), tok_hbhe_, tok_prehbhe_, HcalTopology::unmergeDepthDetId(), and verbosity_.

Referenced by ~HcalCollapseAnalyzer().

                                                                                       {
   if (verbosity_ > 0) 
     edm::LogVerbatim("Collapse") << "Run " << iEvent.id().run() << " Event "
                  << iEvent.id().event() << " starts ==============";
   
   edm::ESHandle<HcalTopology> htopo;
   iSetup.get<HcalRecNumberingRecord>().get(htopo);
   const HcalTopology* theHBHETopology = htopo.product();
 
   edm::Handle<HBHERecHitCollection> hbhereco;
   iEvent.getByToken(tok_hbhe_,      hbhereco);
   edm::Handle<HBHERecHitCollection> hbheprereco;
   iEvent.getByToken(tok_prehbhe_,   hbheprereco);
   if (verbosity_ > 0) {
     edm::LogVerbatim("Collapse") << "Handle Reco "    << hbhereco << " Size " 
                  << hbhereco->size();
     edm::LogVerbatim("Collapse") << "Handle PreReco " << hbheprereco <<" Size "
                  << hbheprereco->size();
   }
   if (hbhereco.isValid() && hbheprereco.isValid()) {
     const HBHERecHitCollection* recohbhe    = hbhereco.product();
     const HBHERecHitCollection* prerecohbhe = hbheprereco.product();
     if (verbosity_ > 0) 
       edm::LogVerbatim("Collapse") << "Size of hbhereco: " << recohbhe->size()
                    << " and hbheprereco: " << prerecohbhe->size();
     double sfrac = (prerecohbhe->size() == 0) ? 1 :
       ((double)(recohbhe->size()))/((double)(prerecohbhe->size()));
     h_sfrac->Fill(sfrac);
     h_size->Fill(recohbhe->size());
     for (const auto & hit : (*recohbhe)) {
       HcalDetId id = hit.id();
       if (((id.subdet() == HcalEndcap) && doHE_) ||
       ((id.subdet() == HcalBarrel) && doHB_)) {
     h_depth->Fill(id.depth());
     std::vector<HcalDetId> ids;
     theHBHETopology->unmergeDepthDetId(id, ids);
     if (verbosity_ > 0) {
       edm::LogVerbatim("Collapse") << id << " is derived from " << ids.size() 
                        << " components";
       for (unsigned int k=0; k<ids.size(); ++k) 
         edm::LogVerbatim("Collapse") << "[" << k << "] " << ids[k];
     }
     h_merge->Fill(ids.size());
     double energy = hit.energy();
     double etot(0);
     unsigned int k(0);
     for (const auto phit : (*prerecohbhe)) {
       if (std::find(ids.begin(),ids.end(),phit.id()) != ids.end()) {
         etot += phit.energy();
         double frac = (energy > 0) ? phit.energy()/energy : 1.0; 
         h_frac->Fill(frac);
         if (verbosity_ > 0) 
           edm::LogVerbatim("Collapse") << "Frac[" << k << "] " << frac 
                        << " for " <<phit.id();
         ++k;
       }
     }
     double frac = (energy > 0) ? etot/energy : 1.0; 
         h_balance->Fill(frac);
     if (verbosity_ > 0) 
       edm::LogVerbatim("Collapse") << "Total energy " << energy << ":" << etot
                        << ":" << frac;
       }
     }
   }
 }

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 167 of file HcalCollapseAnalyzer.cc.

References DQMStore::IBooker::book1D(), DEFINE_FWK_MODULE, h_balance, h_depth, h_frac, h_merge, h_sfrac, h_size, DQMStore::IBooker::setCurrentFolder(), and topFolderName_.

Referenced by ~HcalCollapseAnalyzer().

                                              {
   // Book histograms
   ibooker.setCurrentFolder(topFolderName_);
   h_merge  = ibooker.book1D("h_merge","Number of hits merged",10,0.0,10.0);
   h_size   = ibooker.book1D("h_size","Size of the RecHit collection",100,500.0,1500.0);
   h_depth  = ibooker.book1D("h_depth","Depth of the Id's used",10,0.0,10.0);
   h_sfrac  = ibooker.book1D("h_sfrac","Ratio of sizes of preRecHit and RecHit collections",150,0.0,1.5);
   h_frac   = ibooker.book1D("h_frac", "Fraction of energy before collapse",150,0.0,1.5);
   h_balance= ibooker.book1D("h_balance", "Balance of energy between pre- and post-collapse",100,0.5,1.5);
 }

void HcalCollapseAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 89 of file HcalCollapseAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~HcalCollapseAnalyzer().

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("topFolderName",           "HcalCollapse");
   desc.addUntracked<int>("verbosity",              0);
   desc.addUntracked<edm::InputTag>("recHitHBHE",   edm::InputTag("hbhereco"));
   desc.addUntracked<edm::InputTag>("preRecHitHBHE",edm::InputTag("hbheprereco"));
   desc.addUntracked<bool>("doHE",                  true);
   desc.addUntracked<bool>("doHB",                  false);
   descriptions.add("hcalCollapseAnalyzer",desc);
 }