Inheritance diagram for HGCalSimHitsClient:

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

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

	HGCalSimHitsClient (const edm::ParameterSet &)

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

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

	~HGCalSimHitsClient () override

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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

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

	DQMEDHarvester ()

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

void	endJob () final

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &, edm::EventSetup const &) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) 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::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, ESInputTag const &tag)

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)

Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 31 of file HGCalSimHitsClient.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 50 of file HGCalSimHitsClient.cc.

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

HGCalSimHitsClient::~HGCalSimHitsClient ( )

inlineoverride

Definition at line 42 of file HGCalSimHitsClient.cc.

References beginRun(), EnergyCorrector::c, dqmEndJob(), cuy::ib, findQualityFiles::run, runClient_(), and simHitsEndjob().

42 {}

Member Function Documentation

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

override

Definition at line 55 of file HGCalSimHitsClient.cc.

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

Referenced by ~HGCalSimHitsClient().

                                                                                 { 
 
   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::LogVerbatim("HGCalValidation") << "Initialize HGCalSimHitsClient for "
                     << nameDetector_ << " : " << layers_;
 }

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

overridevirtual

Implements DQMEDHarvester.

Definition at line 73 of file HGCalSimHitsClient.cc.

References runClient_().

Referenced by ~HGCalSimHitsClient().

                                                                            {
   runClient_(ib, ig);
 }

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

virtual

Definition at line 77 of file HGCalSimHitsClient.cc.

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

Referenced by dqmEndJob(), and ~HGCalSimHitsClient().

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

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

loop over timeslice

Definition at line 108 of file HGCalSimHitsClient.cc.

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

Referenced by runClient_(), and ~HGCalSimHitsClient().

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