#include <HGCalRecHitsClient.h>

Inheritance diagram for HGCalRecHitsClient:

Public Member Functions
virtual void	beginRun (const edm::Run &run, const edm::EventSetup &c)

void	dqmEndJob (DQMStore::IBooker &ib, DQMStore::IGetter &ig)

	HGCalRecHitsClient (const edm::ParameterSet &)

int	recHitsEndjob (const std::vector< MonitorElement * > &hcalMEs)

virtual void	runClient_ (DQMStore::IBooker &ib, DQMStore::IGetter &ig)

virtual	~HGCalRecHitsClient ()

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

virtual void	beginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDHarvester (void)

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endJob () final

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

virtual void	endRun (edm::Run const &, edm::EventSetup const &)

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
unsigned int	layers_

std::string	nameDetector_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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)

Detailed Description

Definition at line 33 of file HGCalRecHitsClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file HGCalRecHitsClient.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), nameDetector_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

                                                                      {
   
   verbosity_     = iConfig.getUntrackedParameter<int>("Verbosity",0);
   nameDetector_  = iConfig.getParameter<std::string>("DetectorName");
 }

HGCalRecHitsClient::~HGCalRecHitsClient ( )

virtual

Definition at line 23 of file HGCalRecHitsClient.cc.

23 { }

Member Function Documentation

void HGCalRecHitsClient::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	c
	)

virtual

Definition at line 25 of file HGCalRecHitsClient.cc.

References edm::EventSetup::get(), HcalDDDRecConstants::getMaxDepth(), HGCalDDDConstants::layers(), layers_, and nameDetector_.

                                                                                 { 
   if (nameDetector_ == "HCal") {
     edm::ESHandle<HcalDDDRecConstants> pHRNDC;
     iSetup.get<HcalRecNumberingRecord>().get( pHRNDC );
     const HcalDDDRecConstants *hcons  = &(*pHRNDC);
     layers_ = hcons->getMaxDepth(1);
   } else {
     edm::ESHandle<HGCalDDDConstants>  pHGDC;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, pHGDC);
     const HGCalDDDConstants & hgcons_ = (*pHGDC);
     layers_ = hgcons_.layers(true);
   }
 }

void HGCalRecHitsClient::dqmEndJob	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

virtual

Implements DQMEDHarvester.

Definition at line 39 of file HGCalRecHitsClient.cc.

References runClient_().

                                                                            {
   runClient_(ib, ig);
 }

int HGCalRecHitsClient::recHitsEndjob ( const std::vector< MonitorElement * > & hcalMEs )

Definition at line 72 of file HGCalRecHitsClient.cc.

References DEFINE_FWK_MODULE, mps_fire::i, createfilelist::int, layers_, dataset::name, and nevent.

Referenced by runClient_().

                                                                                 {
   std::vector<MonitorElement*> energy_;
   std::vector<MonitorElement*> EtaPhi_Plus_;
   std::vector<MonitorElement*> EtaPhi_Minus_;
   std::vector<MonitorElement*> HitOccupancy_Plus_;
   std::vector<MonitorElement*> HitOccupancy_Minus_;
   std::vector<MonitorElement*> MeanHitOccupancy_Plus_;
   std::vector<MonitorElement*> MeanHitOccupancy_Minus_;
 
   std::ostringstream name;
   double nevent;
   int nbinsx, nbinsy;
 
   for (unsigned int ilayer = 0; ilayer < layers_; ilayer++ ){ 
     name.str(""); name << "energy_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
     energy_.push_back(hgcalMEs[ih]);
       }
     }
     
     //normalization
     nevent = energy_.at(ilayer)->getEntries();
     nbinsx = energy_.at(ilayer)->getNbinsX();
     for(int i=1; i <= nbinsx; i++) {
       double binValue = energy_.at(ilayer)->getBinContent(i)/nevent;
       energy_.at(ilayer)->setBinContent(i,binValue);
     }
     
     //EtaPhi 2d plots
     name.str(""); name << "EtaPhi_Plus_" << "layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
     EtaPhi_Plus_.push_back(hgcalMEs[ih]);
       }
     }
 
     name.str(""); name << "EtaPhi_Minus_" << "layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
         EtaPhi_Minus_.push_back(hgcalMEs[ih]);
       }
     }
 
     //normalization EtaPhi
     nevent = EtaPhi_Plus_.at(ilayer)->getEntries();
     nbinsx = EtaPhi_Plus_.at(ilayer)->getNbinsX();
     nbinsy = EtaPhi_Plus_.at(ilayer)->getNbinsY();
     for(int i=1; i<= nbinsx; ++i) {
       for(int j=1; j<= nbinsy; ++j) {
     double binValue = EtaPhi_Plus_.at(ilayer)->getBinContent(i, j)/nevent;
     EtaPhi_Plus_.at(ilayer)->setBinContent(i, j, binValue);
       }
     }
     
     nevent =  EtaPhi_Minus_.at(ilayer)->getEntries();
     nbinsx =  EtaPhi_Minus_.at(ilayer)->getNbinsX();
     nbinsy = EtaPhi_Plus_.at(ilayer)->getNbinsY();
     for(int i=1; i<= nbinsx; ++i) {
       for(int j=1; j<= nbinsy; ++j) {
         double binValue =  EtaPhi_Minus_.at(ilayer)->getBinContent(i, j)/nevent;
     EtaPhi_Minus_.at(ilayer)->setBinContent(i, j, binValue);
       }
     }
     
     //HitOccupancy 
     name.str(""); name << "HitOccupancy_Plus_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
     HitOccupancy_Plus_.push_back(hgcalMEs[ih]);
       }
     }
     
     name.str(""); name << "HitOccupancy_Minus_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
     HitOccupancy_Minus_.push_back(hgcalMEs[ih]);
       }
     }
     
     //normalization of hit occupancy histos
     nevent = HitOccupancy_Plus_.at(ilayer)->getEntries();
     nbinsx = HitOccupancy_Plus_.at(ilayer)->getNbinsX();
     for(int i=1; i<= nbinsx; ++i) {
       double binValue = HitOccupancy_Plus_.at(ilayer)->getBinContent(i)/nevent;
       HitOccupancy_Plus_.at(ilayer)->setBinContent(i, binValue);
     }
     
     nevent = HitOccupancy_Minus_.at(ilayer)->getEntries();
     nbinsx = HitOccupancy_Minus_.at(ilayer)->getNbinsX();
     for(int i=1; i<= nbinsx; ++i) {
       double binValue = HitOccupancy_Minus_.at(ilayer)->getBinContent(i)/nevent;
       HitOccupancy_Minus_.at(ilayer)->setBinContent(i, binValue);
     }
     
   }//loop over layers
   
   name.str(""); name << "SUMOfRecHitOccupancy_Plus";
   for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
       MeanHitOccupancy_Plus_.push_back(hgcalMEs[ih]);
       unsigned int indx = MeanHitOccupancy_Plus_.size()-1;
       for (int ilayer=0; ilayer < (int)layers_; ++ilayer) {
     double meanVal = HitOccupancy_Plus_.at(ilayer)->getMean();
     MeanHitOccupancy_Plus_[indx]->setBinContent(ilayer+1, meanVal);
       }
       break;
     }
   }
   
   name.str(""); name << "SUMOfRecHitOccupancy_Plus";
   for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
       MeanHitOccupancy_Minus_.push_back(hgcalMEs[ih]);
       unsigned indx = MeanHitOccupancy_Minus_.size()-1;
       for (int ilayer=0; ilayer < (int)layers_; ++ilayer) {
     double meanVal = HitOccupancy_Minus_.at(ilayer)->getMean();
     MeanHitOccupancy_Minus_[indx]->setBinContent(ilayer+1, meanVal);
       }
       break;
     }
   }
   
   return 1;
 }

void HGCalRecHitsClient::runClient_	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

virtual

Definition at line 43 of file HGCalRecHitsClient.cc.

References DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), mps_fire::i, nameDetector_, recHitsEndjob(), DQMStore::IGetter::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

Referenced by dqmEndJob().

                                                                             {
 
   ig.setCurrentFolder("/"); 
   if (verbosity_>0) edm::LogInfo("HGCalValidation") << "\nrunClient";
   std::vector<MonitorElement*> hgcalMEs;
   std::vector<std::string> fullDirPath = ig.getSubdirs();
 
   for (unsigned int i=0; i<fullDirPath.size(); i++) {
     if (verbosity_>0) 
       edm::LogInfo("HGCalValidation") << "\nfullPath: " << fullDirPath.at(i);
     ig.setCurrentFolder(fullDirPath.at(i));
     std::vector<std::string> fullSubDirPath = ig.getSubdirs();
 
     for (unsigned int j=0; j<fullSubDirPath.size(); j++) {
       if (verbosity_>1) 
     edm::LogInfo("HGCalValidation") << "fullSubPath: " << fullSubDirPath.at(j);
       std::string nameDirectory = "HGCalRecHitsV/"+nameDetector_;
       if (strcmp(fullSubDirPath.at(j).c_str(), nameDirectory.c_str()) == 0) {
         hgcalMEs = ig.getContents(fullSubDirPath.at(j));
         if (verbosity_>1) 
       edm::LogInfo("HGCalValidation") << "hgcalMES size : " << hgcalMEs.size();
         if (!recHitsEndjob(hgcalMEs)) 
       edm::LogWarning("HGCalValidation") << "\nError in RecHitsEndjob!";
       }
     }
 
   }
 }