#include <HcalSimHitsClient.h>

Inheritance diagram for HcalSimHitsClient:

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

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

	HcalSimHitsClient (const edm::ParameterSet &)

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

	~HcalSimHitsClient () 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::vector< std::string >	getHistogramTypes ()

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

Private Attributes
std::string	dirName_

const HcalDDDRecConstants *	hcons

int	maxDepthHB_

int	maxDepthHE_

int	maxDepthHF_

int	maxDepthHO_

bool	verbose_

Static Private Attributes
static const int	nTime = 4

static const int	nType1 = 4

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 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 36 of file HcalSimHitsClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file HcalSimHitsClient.cc.

References dirName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

                                                                    {
   dirName_ = iConfig.getParameter<std::string>("DQMDirName");
   verbose_ = iConfig.getUntrackedParameter<bool>("Verbosity", false);
 }

HcalSimHitsClient::~HcalSimHitsClient ( )

override

Definition at line 19 of file HcalSimHitsClient.cc.

19 {}

Member Function Documentation

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

override

Definition at line 21 of file HcalSimHitsClient.cc.

References edm::EventSetup::get(), HcalDDDRecConstants::getMaxDepth(), hcons, maxDepthHB_, maxDepthHE_, maxDepthHF_, and maxDepthHO_.

                                                                           {
   edm::ESHandle<HcalDDDRecConstants> pHRNDC;
   c.get<HcalRecNumberingRecord>().get(pHRNDC);
   hcons = &(*pHRNDC);
   maxDepthHB_ = hcons->getMaxDepth(0);
   maxDepthHE_ = hcons->getMaxDepth(1);
   maxDepthHF_ = hcons->getMaxDepth(2);
   maxDepthHO_ = hcons->getMaxDepth(3);
 
   edm::LogInfo("HitsValidationHcal") << " Maximum Depths HB:" << maxDepthHB_ << " HE:" << maxDepthHE_
                                      << " HO:" << maxDepthHO_ << " HF:" << maxDepthHF_;
 }

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

overridevirtual

Implements DQMEDHarvester.

Definition at line 34 of file HcalSimHitsClient.cc.

References runClient_().

34 { runClient_(ib, ig); }

HcalSimHitsClient::runClient_

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

Definition: HcalSimHitsClient.cc:36

std::vector< std::string > HcalSimHitsClient::getHistogramTypes ( )

private

Definition at line 172 of file HcalSimHitsClient.cc.

References DEFINE_FWK_MODULE, egammaForCoreTracking_cff::depth, gen::k, SiStripPI::max, CMSBoostedTauSeedingParameters_cfi::maxDepth, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

Referenced by SimHitsEndjob().

                                                           {
   int maxDepth = std::max(maxDepthHB_, maxDepthHE_);
   if (verbose_)
     edm::LogInfo("HitsValidationHcal") << "Max depth 1st step:: " << maxDepth;
   maxDepth = std::max(maxDepth, maxDepthHF_);
   if (verbose_)
     edm::LogInfo("HitsValidationHcal") << "Max depth 2nd step:: " << maxDepth;
   maxDepth = std::max(maxDepth, maxDepthHO_);
   if (verbose_)
     edm::LogInfo("HitsValidationHcal") << "Max depth 3rd step:: " << maxDepth;
   std::vector<std::string> divisions;
   char name1[20];
 
   // first overall Hcal
   for (int depth = 0; depth < maxDepth; ++depth) {
     sprintf(name1, "HC%d", depth);
     divisions.push_back(std::string(name1));
   }
   // HB
   for (int depth = 0; depth < maxDepthHB_; ++depth) {
     sprintf(name1, "HB%d", depth);
     divisions.push_back(std::string(name1));
   }
   // HE
   for (int depth = 0; depth < maxDepthHE_; ++depth) {
     sprintf(name1, "HE%d+z", depth);
     divisions.push_back(std::string(name1));
     sprintf(name1, "HE%d-z", depth);
     divisions.push_back(std::string(name1));
   }
   // HO
   {
     int depth = maxDepthHO_;
     sprintf(name1, "HO%d", depth);
     divisions.push_back(std::string(name1));
   }
   // HF (first absorber, then different types of abnormal hits)
   std::string hfty1[4] = {"A", "W", "B", "J"};
   for (int k = 0; k < 4; ++k) {
     for (int depth = 0; depth < maxDepthHF_; ++depth) {
       sprintf(name1, "HF%s%d+z", hfty1[k].c_str(), depth);
       divisions.push_back(std::string(name1));
       sprintf(name1, "HF%s%d-z", hfty1[k].c_str(), depth);
       divisions.push_back(std::string(name1));
     }
   }
   return divisions;
 }

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

virtual

Definition at line 36 of file HcalSimHitsClient.cc.

References dirName_, DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), mps_fire::i, DQMStore::IGetter::setCurrentFolder(), SimHitsEndjob(), and verbose_.

Referenced by dqmEndJob().

                                                                            {
   ig.setCurrentFolder(dirName_);
 
   if (verbose_)
     edm::LogInfo("HitsValidationHcal") << "runClient";
 
   std::vector<MonitorElement *> hcalMEs;
 
   std::vector<std::string> fullPathHLTFolders = ig.getSubdirs();
   for (unsigned int i = 0; i < fullPathHLTFolders.size(); i++) {
     if (verbose_)
       edm::LogInfo("HitsValidationHcal") << "fullPath: " << fullPathHLTFolders[i];
     ig.setCurrentFolder(fullPathHLTFolders[i]);
 
     std::vector<std::string> fullSubPathHLTFolders = ig.getSubdirs();
     for (unsigned int j = 0; j < fullSubPathHLTFolders.size(); j++) {
       if (verbose_)
         edm::LogInfo("HitsValidationHcal") << "fullSub: " << fullSubPathHLTFolders[j];
 
       if (strcmp(fullSubPathHLTFolders[j].c_str(), "HcalHitsV/SimHitsValidationHcal") == 0) {
         hcalMEs = ig.getContents(fullSubPathHLTFolders[j]);
         if (verbose_)
           edm::LogInfo("HitsValidationHcal") << "hltMES size : " << hcalMEs.size();
         if (!SimHitsEndjob(hcalMEs))
           edm::LogWarning("HitsValidationHcal") << "Error in SimhitEndjob!";
       }
     }
   }
 }

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

private

Definition at line 68 of file HcalSimHitsClient.cc.

References objects.autophobj::float, MonitorElement::getBinContent(), MonitorElement::getEntries(), getHistogramTypes(), getName(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), mps_fire::i, gen::k, dataset::name, nevent, nTime, nType1, MonitorElement::setBinContent(), AlCaHLTBitMon_QueryRunRegistry::string, ntuplemaker::time, and verbose_.

Referenced by runClient_().

                                                                                {
   std::vector<std::string> divisions = getHistogramTypes();
   MonitorElement *Occupancy_map[nTime][divisions.size()];
   MonitorElement *Energy[nType1], *Time_weighteden[nType1];
   MonitorElement *HitEnergyvsieta[divisions.size()], *HitTimevsieta[divisions.size()];
 
   std::string time[nTime] = {"25", "50", "100", "250"};
   std::string detdivision[nType1] = {"HB", "HE", "HF", "HO"};
   char name[40], name1[40], name2[40], name3[40], name4[40];
 
   for (int k = 0; k < nType1; k++) {
     Energy[k] = nullptr;
     Time_weighteden[k] = nullptr;
     for (unsigned int ih = 0; ih < hcalMEs.size(); ih++) {
       sprintf(name1, "Energy_%s", detdivision[k].c_str());
       sprintf(name2, "Time_Enweighted_%s", detdivision[k].c_str());
       if (strcmp(hcalMEs[ih]->getName().c_str(), name1) == 0) {
         Energy[k] = hcalMEs[ih];
       }
       if (strcmp(hcalMEs[ih]->getName().c_str(), name2) == 0) {
         Time_weighteden[k] = hcalMEs[ih];
       }
     }
   }
 
   for (int i = 0; i < nTime; i++) {
     for (unsigned int j = 0; j < divisions.size(); j++) {
       for (unsigned int ih = 0; ih < hcalMEs.size(); ih++) {
         sprintf(name, "HcalHitE%s%s", time[i].c_str(), divisions[j].c_str());
         if (strcmp(hcalMEs[ih]->getName().c_str(), name) == 0) {
           Occupancy_map[i][j] = hcalMEs[ih];
         }
       }
     }
   }
 
   for (unsigned int k = 0; k < divisions.size(); k++) {
     HitEnergyvsieta[k] = nullptr;
     HitTimevsieta[k] = nullptr;
     for (unsigned int ih = 0; ih < hcalMEs.size(); ih++) {
       sprintf(name3, "HcalHitEta%s", divisions[k].c_str());
       sprintf(name4, "HcalHitTimeAEta%s", divisions[k].c_str());
       if (strcmp(hcalMEs[ih]->getName().c_str(), name3) == 0) {
         HitEnergyvsieta[k] = hcalMEs[ih];
       }
       if (strcmp(hcalMEs[ih]->getName().c_str(), name4) == 0) {
         HitTimevsieta[k] = hcalMEs[ih];
       }
     }
   }
 
   // mean energy
 
   double nevent = Energy[0]->getEntries();
   if (verbose_)
     edm::LogInfo("HitsValidationHcal") << "nevent : " << nevent;
 
   float cont[nTime][divisions.size()];
   float en[nType1], tme[nType1];
   float hitenergy[divisions.size()], hittime[divisions.size()];
   float fev = float(nevent);
 
   for (int dettype = 0; dettype < nType1; dettype++) {
     int nx1 = Energy[dettype]->getNbinsX();
     for (int i = 0; i <= nx1; i++) {
       en[dettype] = Energy[dettype]->getBinContent(i) / fev;
       Energy[dettype]->setBinContent(i, en[dettype]);
     }
     int nx2 = Time_weighteden[dettype]->getNbinsX();
     for (int i = 0; i <= nx2; i++) {
       tme[dettype] = Time_weighteden[dettype]->getBinContent(i) / fev;
       Time_weighteden[dettype]->setBinContent(i, tme[dettype]);
     }
   }
 
   for (unsigned int dettype = 0; dettype < divisions.size(); dettype++) {
     int nx1 = HitEnergyvsieta[dettype]->getNbinsX();
     for (int i = 0; i <= nx1; i++) {
       hitenergy[dettype] = HitEnergyvsieta[dettype]->getBinContent(i) / fev;
       HitEnergyvsieta[dettype]->setBinContent(i, hitenergy[dettype]);
     }
     int nx2 = HitTimevsieta[dettype]->getNbinsX();
     for (int i = 0; i <= nx2; i++) {
       hittime[dettype] = HitTimevsieta[dettype]->getBinContent(i) / fev;
       HitTimevsieta[dettype]->setBinContent(i, hittime[dettype]);
     }
   }
 
   for (int itime = 0; itime < nTime; itime++) {
     for (unsigned int det = 0; det < divisions.size(); det++) {
       int ny = Occupancy_map[itime][det]->getNbinsY();
       int nx = Occupancy_map[itime][det]->getNbinsX();
       for (int i = 1; i < nx + 1; i++) {
         for (int j = 1; j < ny + 1; j++) {
           cont[itime][det] = Occupancy_map[itime][det]->getBinContent(i, j) / fev;
           Occupancy_map[itime][det]->setBinContent(i, j, cont[itime][det]);
         }
       }
     }
   }
 
   return 1;
 }