#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

	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

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

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	hcalTopologyToken_

const edm::InputTag	preRecHitHBHE_

const edm::InputTag	recHitHBHE_

const HcalTopology *	theHBHETopology = nullptr

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

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 44 of file HcalCollapseAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalCollapseAnalyzer()

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

explicit

Definition at line 70 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)),
       hcalTopologyToken_{esConsumes<HcalTopology, HcalRecNumberingRecord>()} {
   // 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::~HcalCollapseAnalyzer ( )

inlineoverride

Definition at line 47 of file HcalCollapseAnalyzer.cc.

47 {}

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.

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

                                                                                       {
   if (verbosity_ > 0)
     edm::LogVerbatim("Collapse") << "Run " << iEvent.id().run() << " Event " << iEvent.id().event()
                                  << " starts ==============";
 
   theHBHETopology = &iSetup.getData(hcalTopologyToken_);
 
   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;
       }
     }
   }
 }

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 155 of file HcalCollapseAnalyzer.cc.

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

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

◆ fillDescriptions()

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

static

Definition at line 86 of file HcalCollapseAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_2022v12_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

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

Member Data Documentation

◆ doHB_

const bool HcalCollapseAnalyzer::doHB_

private

Definition at line 59 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

◆ doHE_

const bool HcalCollapseAnalyzer::doHE_

private

Definition at line 59 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

◆ h_balance

MonitorElement * HcalCollapseAnalyzer::h_balance

private

Definition at line 67 of file HcalCollapseAnalyzer.cc.

Referenced by analyze(), and bookHistograms().

◆ h_depth

MonitorElement * HcalCollapseAnalyzer::h_depth

private

Definition at line 66 of file HcalCollapseAnalyzer.cc.

Referenced by analyze(), and bookHistograms().

◆ h_frac

MonitorElement * HcalCollapseAnalyzer::h_frac

private

Definition at line 67 of file HcalCollapseAnalyzer.cc.

Referenced by analyze(), and bookHistograms().

◆ h_merge

MonitorElement* HcalCollapseAnalyzer::h_merge

private

Definition at line 66 of file HcalCollapseAnalyzer.cc.

Referenced by analyze(), and bookHistograms().

◆ h_sfrac

MonitorElement* HcalCollapseAnalyzer::h_sfrac

private

Definition at line 67 of file HcalCollapseAnalyzer.cc.

Referenced by analyze(), and bookHistograms().

◆ h_size

MonitorElement * HcalCollapseAnalyzer::h_size

private

Definition at line 66 of file HcalCollapseAnalyzer.cc.

Referenced by analyze(), and bookHistograms().

◆ hcalTopologyToken_

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> HcalCollapseAnalyzer::hcalTopologyToken_

private

Definition at line 64 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

◆ preRecHitHBHE_

const edm::InputTag HcalCollapseAnalyzer::preRecHitHBHE_

private

Definition at line 58 of file HcalCollapseAnalyzer.cc.

◆ recHitHBHE_

const edm::InputTag HcalCollapseAnalyzer::recHitHBHE_

private

Definition at line 58 of file HcalCollapseAnalyzer.cc.

◆ theHBHETopology

const HcalTopology* HcalCollapseAnalyzer::theHBHETopology = nullptr

private

Definition at line 60 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalCollapseAnalyzer::tok_hbhe_

private

Definition at line 62 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

◆ tok_prehbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalCollapseAnalyzer::tok_prehbhe_

private

Definition at line 63 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

◆ topFolderName_

const std::string HcalCollapseAnalyzer::topFolderName_

private

Definition at line 56 of file HcalCollapseAnalyzer.cc.

Referenced by bookHistograms().

◆ verbosity_

const int HcalCollapseAnalyzer::verbosity_

private

Definition at line 57 of file HcalCollapseAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalCollapseAnalyzer()

◆ ~HcalCollapseAnalyzer()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ fillDescriptions()

Member Data Documentation

◆ doHB_

◆ doHE_

◆ h_balance

◆ h_depth

◆ h_frac

◆ h_merge

◆ h_sfrac

◆ h_size

◆ hcalTopologyToken_

◆ preRecHitHBHE_

◆ recHitHBHE_

◆ theHBHETopology

◆ tok_hbhe_

◆ tok_prehbhe_

◆ topFolderName_

◆ verbosity_