#include <DQMOffline/Hcal/src/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

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
	)

overrideprivate

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

overrideprivate

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