#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

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 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 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 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 21 of file HcalSimHitsClient.cc.

21 { }

Member Function Documentation

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

override

Definition at line 23 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 37 of file HcalSimHitsClient.cc.

References runClient_().

                                                                          {
   runClient_(ib,ig);
 }

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

private

Definition at line 177 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 41 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 73 of file HcalSimHitsClient.cc.

References objects.autophobj::float, MonitorElement::getBinContent(), MonitorElement::getEntries(), getHistogramTypes(), 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;
 }