#include <HGCalDigiClient.h>

Inheritance diagram for HGCalDigiClient:

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

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

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

	HGCalDigiClient (const edm::ParameterSet &)

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

virtual	~HGCalDigiClient ()

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 &	itemsToGetFromEvent () 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
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 31 of file HGCalDigiClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 15 of file HGCalDigiClient.cc.

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

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

HGCalDigiClient::~HGCalDigiClient ( )

virtual

Definition at line 21 of file HGCalDigiClient.cc.

21 { }

Member Function Documentation

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

virtual

Definition at line 23 of file HGCalDigiClient.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);
   }
 }

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

Definition at line 71 of file HGCalDigiClient.cc.

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

Referenced by runClient_().

                                                                            {
 
   std::vector<MonitorElement*> charge_;
   std::vector<MonitorElement*> DigiOccupancy_XY_;
   std::vector<MonitorElement*> ADC_;
   std::vector<MonitorElement*> DigiOccupancy_Plus_;
   std::vector<MonitorElement*> DigiOccupancy_Minus_;
   std::vector<MonitorElement*> MeanDigiOccupancy_Plus_;
   std::vector<MonitorElement*> MeanDigiOccupancy_Minus_;
   std::ostringstream name;
   double nevent;
   int nbinsx, nbinsy;
   
   for (int ilayer = 0; ilayer < layers_; ilayer++ ){ 
     //charge
     name.str(""); name << "charge_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0)
     charge_.push_back(hgcalMEs[ih]);
     }
     //normalization
     nevent = charge_.at(ilayer)->getEntries();
     nbinsx = charge_.at(ilayer)->getNbinsX();
     for(int i=1; i <= nbinsx; i++) {
       double binValue = charge_.at(ilayer)->getBinContent(i)/nevent;
       charge_.at(ilayer)->setBinContent(i,binValue);
     }
     
     //XY 2d plots
     name.str(""); name << "DigiOccupancy_XY_" << "layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0)
     DigiOccupancy_XY_.push_back(hgcalMEs[ih]);
     }
     
     //normalization of XY 2d
     nevent = DigiOccupancy_XY_.at(ilayer)->getEntries();
     nbinsx = DigiOccupancy_XY_.at(ilayer)->getNbinsX();
     nbinsy = DigiOccupancy_XY_.at(ilayer)->getNbinsY();
     for(int i=1; i<= nbinsx; ++i) {
       for(int j=1; j<= nbinsy; ++j) {
     double binValue = DigiOccupancy_XY_.at(ilayer)->getBinContent(i, j)/nevent;
     DigiOccupancy_XY_.at(ilayer)->setBinContent(i, j, binValue);
       }
     }
     
     //ADC
     name.str(""); name << "ADC_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0)
     ADC_.push_back(hgcalMEs[ih]);
     }
     
     //normalization of ADC Histos
     nevent = ADC_.at(ilayer)->getEntries();
     nbinsx = ADC_.at(ilayer)->getNbinsX();
     for(int i=1; i<= nbinsx; ++i) {
       double binValue = ADC_.at(ilayer)->getBinContent(i)/nevent;
       ADC_.at(ilayer)->setBinContent(i, binValue);
     }
     
     //Digi Occupancy
     name.str(""); name << "DigiOccupancy_Plus_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
         DigiOccupancy_Plus_.push_back(hgcalMEs[ih]);
       }
     }
     
     name.str(""); name << "DigiOccupancy_Minus_layer_" << ilayer;
     for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
       if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
         DigiOccupancy_Minus_.push_back(hgcalMEs[ih]);
       }
     }
     
     //normalization of Digi Occupancy histos
     nevent = DigiOccupancy_Plus_.at(ilayer)->getEntries();
     nbinsx = DigiOccupancy_Plus_.at(ilayer)->getNbinsX();
     for(int i=1; i<= nbinsx; ++i) {
       double binValue = DigiOccupancy_Plus_.at(ilayer)->getBinContent(i)/nevent;
       DigiOccupancy_Plus_.at(ilayer)->setBinContent(i, binValue);
     }
     
     nevent = DigiOccupancy_Minus_.at(ilayer)->getEntries();
     nbinsx = DigiOccupancy_Minus_.at(ilayer)->getNbinsX();
     for(int i=1; i<= nbinsx; ++i) {
       double binValue = DigiOccupancy_Minus_.at(ilayer)->getBinContent(i)/nevent;
       DigiOccupancy_Minus_.at(ilayer)->setBinContent(i, binValue);
     }
 
   }//loop over layers
   
   name.str(""); name << "SUMOfDigiOccupancy_Plus";
   for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0) {
       MeanDigiOccupancy_Plus_.push_back(hgcalMEs[ih]);
       unsigned indx = MeanDigiOccupancy_Plus_.size()-1;
       for(int ilayer=0; ilayer < (int)layers_; ++ilayer) {
     double meanVal = DigiOccupancy_Plus_.at(ilayer)->getMean();
     MeanDigiOccupancy_Plus_[indx]->setBinContent(ilayer+1, meanVal);
       }
       break;
     }
   }
 
   name.str(""); name << "SUMOfDigiOccupancy_Plus";
   for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
       MeanDigiOccupancy_Minus_.push_back(hgcalMEs[ih]);
       unsigned indx = MeanDigiOccupancy_Minus_.size()-1;
       for(int ilayer=0; ilayer < (int)layers_; ++ilayer) {
     double meanVal = DigiOccupancy_Minus_.at(ilayer)->getMean();
     MeanDigiOccupancy_Minus_[indx]->setBinContent(ilayer+1, meanVal);
       }
       break;
     }
   }
   return 1;
 }

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

virtual

Implements DQMEDHarvester.

Definition at line 38 of file HGCalDigiClient.cc.

References runClient_().

                                                                         {
   runClient_(ib, ig);
 }

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

Definition at line 42 of file HGCalDigiClient.cc.

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

Referenced by dqmEndJob().

                                                                          {
 
   ig.setCurrentFolder("/"); 
   if (verbosity_) 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_) 
       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_) 
     edm::LogInfo("HGCalValidation") << "fullSubPath: " << fullSubDirPath.at(j);
       std::string nameDirectory = "HGCalDigiV/"+nameDetector_;
       if (strcmp(fullSubDirPath.at(j).c_str(), nameDirectory.c_str()) == 0) {
         hgcalMEs = ig.getContents(fullSubDirPath.at(j));
         if (verbosity_) 
       edm::LogInfo("HGCalValidation") << "hgcalMES size : " << hgcalMEs.size();
         if (!digisEndjob(hgcalMEs)) 
       edm::LogInfo("HGCalValidation") << "\nError in DigisEndjob!";
       }
     }
 
   }
 }