#include <HGCalSimHitsClient.h>

Inheritance diagram for HGCalSimHitsClient:

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

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

	HGCalSimHitsClient (const edm::ParameterSet &)

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

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

virtual	~HGCalSimHitsClient ()

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
unsigned int	layers_

std::string	nameDetector_

int	nTimes_

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 37 of file HGCalSimHitsClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file HGCalSimHitsClient.cc.

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

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

HGCalSimHitsClient::~HGCalSimHitsClient ( )

virtual

Definition at line 23 of file HGCalSimHitsClient.cc.

23 { }

Member Function Documentation

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

virtual

Definition at line 25 of file HGCalSimHitsClient.cc.

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

                                                                                 { 
 
   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(false);
   }
   if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << "Initialize HGCalSimHitsClient for " 
                     << nameDetector_ << " : " << layers_ <<"\n";
 }

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

virtual

Implements DQMEDHarvester.

Definition at line 43 of file HGCalSimHitsClient.cc.

References runClient_().

                                                                            {
   runClient_(ib, ig);
 }

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

virtual

Definition at line 47 of file HGCalSimHitsClient.cc.

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

Referenced by dqmEndJob().

                                                                             {
 
   ig.setCurrentFolder("/"); 
   if (verbosity_>0) edm::LogInfo("HGCalValidation") << " runClient\n";
   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") << "fullPath: " << fullDirPath.at(i) << "\n";
     ig.setCurrentFolder(fullDirPath.at(i));
     std::vector<std::string> fullSubDirPath = ig.getSubdirs();
 
     for (unsigned int j=0; j<fullSubDirPath.size(); j++) {
       if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << "fullSubPath: " << fullSubDirPath.at(j) << "\n";
       std::string nameDirectory = "HGCalSimHitsV/"+nameDetector_;
 
       if (strcmp(fullSubDirPath.at(j).c_str(), nameDirectory.c_str()) == 0) {
         hgcalMEs = ig.getContents(fullSubDirPath.at(j));
     if (verbosity_>0) 
       edm::LogInfo("HGCalValidation") << "hgcalMES size : " << hgcalMEs.size() << "\n";
         if (!simHitsEndjob(hgcalMEs)) 
       edm::LogWarning("HGCalValidation") << "\nError in SimhitsEndjob!";
       }
     }
   }
 }

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

loop over timeslice

Definition at line 76 of file HGCalSimHitsClient.cc.

References DEFINE_FWK_MODULE, MonitorElement::getMean(), i, createfilelist::int, j, layers_, dataset::name, nevent, nTimes_, and MonitorElement::setBinContent().

Referenced by runClient_().

                                                                                 {
 
   std::vector<MonitorElement*> energy_[6];
   std::vector<MonitorElement*> EtaPhi_Plus_,       EtaPhi_Minus_;
   std::vector<MonitorElement*> HitOccupancy_Plus_, HitOccupancy_Minus_;
   MonitorElement              *MeanHitOccupancy_Plus_, *MeanHitOccupancy_Minus_;
 
   std::ostringstream name;
   double nevent;
   int nbinsx, nbinsy;
   for (unsigned int ilayer = 0; ilayer < layers_; ilayer++ ){ 
     for (int itimeslice = 0; itimeslice < nTimes_; itimeslice++ ) {
       //Energy
       name.str(""); name << "energy_time_" << itimeslice << "_layer_" << ilayer;
       for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0){
       energy_[itimeslice].push_back(hgcalMEs[ih]);
     }
       }
       //normalization
       nevent = energy_[itimeslice].at(ilayer)->getEntries();
       nbinsx = energy_[itimeslice].at(ilayer)->getNbinsX();
       for(int i=1; i <= nbinsx; i++) {
     double binValue = energy_[itimeslice].at(ilayer)->getBinContent(i)/nevent;
     energy_[itimeslice].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]);
       }
     }
 
     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 << "MeanHitOccupancy_Plus";
   for (unsigned int ih=0; ih<hgcalMEs.size(); ih++) {
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0) {
       MeanHitOccupancy_Plus_ = hgcalMEs[ih];
       for (int ilayer=0; ilayer < (int)layers_; ++ilayer) {
     double meanVal = HitOccupancy_Plus_.at(ilayer)->getMean();
     MeanHitOccupancy_Plus_->setBinContent(ilayer+1, meanVal);
       }
       break;
     }
   }
 
   name.str(""); name << "MeanHitOccupancy_Minus";
   for(unsigned int ih=0; ih<hgcalMEs.size(); ih++){
     if (strcmp(hgcalMEs[ih]->getName().c_str(), name.str().c_str()) == 0) {
       MeanHitOccupancy_Minus_ = hgcalMEs[ih];
       for (int ilayer=0; ilayer < (int)layers_; ++ilayer) {
     double meanVal = HitOccupancy_Minus_.at(ilayer)->getMean();
     MeanHitOccupancy_Minus_->setBinContent(ilayer+1, meanVal);
       }
       break;
     }
   }
 
   return 1;
 }