#include <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 &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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

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 &)

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 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 >
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

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

DQMOffline/Hcal/src/HcalCollapseAnalyzer.cc

Description: Studies the collapsing of HcalRecHits

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

Definition at line 47 of file HcalCollapseAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalCollapseAnalyzer()

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

explicit

Definition at line 71 of file HcalCollapseAnalyzer.cc.

     : 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_;
 }

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

◆ ~HcalCollapseAnalyzer()

HcalCollapseAnalyzer::~HcalCollapseAnalyzer ( )

inlineoverride

Definition at line 50 of file HcalCollapseAnalyzer.cc.

50 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 97 of file HcalCollapseAnalyzer.cc.

                                                                                       {
   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->empty()) ? 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;
       }
     }
   }
 }

References LEDCalibrationChannels::depth, doHB_, doHE_, edm::SortedCollection< T, SORT >::empty(), HCALHighEnergyHPDFilter_cfi::energy, dqm::impl::MonitorElement::Fill(), spr::find(), cropTnPTrees::frac, edm::EventSetup::get(), get, h_balance, h_depth, h_frac, h_merge, h_sfrac, h_size, hcalLocalReco_cff::hbheprereco, hcalcalib_dqm_sourceclient-live_cfg::hbhereco, HcalBarrel, HcalEndcap, hit::id, iEvent, dqmdumpme::k, edm::ESHandle< T >::product(), edm::SortedCollection< T, SORT >::size(), tok_hbhe_, tok_prehbhe_, HcalTopology::unmergeDepthDetId(), and verbosity_.

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 157 of file HcalCollapseAnalyzer.cc.

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

References dqm::implementation::IBooker::book1D(), h_balance, h_depth, h_frac, h_merge, h_sfrac, h_size, dqm::implementation::NavigatorBase::setCurrentFolder(), and topFolderName_.

◆ fillDescriptions()

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

static

Definition at line 86 of file HcalCollapseAnalyzer.cc.

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