#include <HcalCollapseAnalyzer.cc>

Inheritance diagram for HcalCollapseAnalyzer:

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

	~HcalCollapseAnalyzer () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::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_

Detailed Description

DQMOffline/Hcal/src/HcalCollapseAnalyzer.cc

Description: Studies the collapsing of HcalRecHits

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

Definition at line 48 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
	)

overrideprivate

Definition at line 98 of file HcalCollapseAnalyzer.cc.

References egammaForCoreTracking_cff::depth, doHB_, doHE_, edm::SortedCollection< T, SORT >::empty(), edm::EventID::event(), MonitorElement::Fill(), spr::find(), DivergingColor::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, 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->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;
       }
     }
   }
 }

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

overrideprivate

Definition at line 158 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 87 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);
 }