#include <Validation/HcalDigisClient/src/HcalDigisClient.cc>

Inheritance diagram for HcalDigisClient:

Classes
struct	HistLim

Public Member Functions
	HcalDigisClient (const edm::ParameterSet &)

	~HcalDigisClient () 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
void	book1D (DQMStore::IBooker &ib, std::string name, int n, double min, double max)

void	book1D (DQMStore::IBooker &ib, std::string name, const HistLim &limX)

void	book2D (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)

void	bookPf (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)

void	bookPf (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY, const char *option)

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

void	fill1D (std::string name, double X, double weight=1)

void	fill2D (std::string name, double X, double Y, double weight=1)

void	fillPf (std::string name, double X, double Y)

int	HcalDigisEndjob (const std::vector< MonitorElement * > &hcalMEs, std::string subdet_, DQMStore::IBooker &ib)

double	integralMETH2D (MonitorElement *ME, int i0, int i1, int j0, int j1)

MonitorElement *	monitor (std::string name)

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

void	scaleMETH2D (MonitorElement *ME, double s)

std::string	str (int x)

Private Attributes
std::string	dirName_

std::map< std::string, MonitorElement * > *	msm_

std::string	outputFile_

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 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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 27 of file HcalDigisClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 23 of file HcalDigisClient.cc.

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

                                                                {
     outputFile_ = iConfig.getUntrackedParameter<std::string > ("outputFile", "HcalDigisClient.root");
     dirName_ = iConfig.getParameter<std::string > ("DQMDirName");
     msm_ = new std::map<std::string, MonitorElement*>();
 
 
 }

HcalDigisClient::~HcalDigisClient ( )

override

Definition at line 32 of file HcalDigisClient.cc.

References msm_.

                                   {
     delete msm_;
 }

Member Function Documentation

void HcalDigisClient::book1D	(	DQMStore::IBooker &	ib,
		std::string	name,
		int	n,
		double	min,
		double	max
	)

inlineprivate

Definition at line 55 of file HcalDigisClient.h.

References DQMStore::IBooker::book1D(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and dataset::name.

Referenced by HcalDigisEndjob().

                                                                                     {
         if (!msm_->count(name)) (*msm_)[name] = ib.book1D(name.c_str(), name.c_str(), n, min, max);
     }

void HcalDigisClient::book1D	(	DQMStore::IBooker &	ib,
		std::string	name,
		const HistLim &	limX
	)

inlineprivate

Definition at line 59 of file HcalDigisClient.h.

References DQMStore::IBooker::book1D(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and dataset::name.

                                                                           {
         if (!msm_->count(name)) (*msm_)[name] = ib.book1D(name.c_str(), name.c_str(), limX.n, limX.min, limX.max);
     }

void HcalDigisClient::book2D	(	DQMStore::IBooker &	ib,
		std::string	name,
		const HistLim &	limX,
		const HistLim &	limY
	)

inlineprivate

Definition at line 67 of file HcalDigisClient.h.

References DQMStore::IBooker::book2D(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and dataset::name.

                                                                                                {
         if (!msm_->count(name)) (*msm_)[name] = ib.book2D(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max);
     }

void HcalDigisClient::bookPf	(	DQMStore::IBooker &	ib,
		std::string	name,
		const HistLim &	limX,
		const HistLim &	limY
	)

inlineprivate

Definition at line 75 of file HcalDigisClient.h.

References DQMStore::IBooker::bookProfile(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and dataset::name.

                                                                                                {
         if (!msm_->count(name)) (*msm_)[name] = ib.bookProfile(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max);
     }

void HcalDigisClient::bookPf	(	DQMStore::IBooker &	ib,
		std::string	name,
		const HistLim &	limX,
		const HistLim &	limY,
		const char *	option
	)

inlineprivate

Definition at line 79 of file HcalDigisClient.h.

References DQMStore::IBooker::bookProfile(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, dataset::name, and TSGForRoadSearch_cfi::option.

                                                                                                                    {
         if (!msm_->count(name)) (*msm_)[name] = ib.bookProfile(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max, option);
     }

void HcalDigisClient::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

inlineoverrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 35 of file HcalDigisClient.h.

References dirName_, runClient(), and DQMStore::IGetter::setCurrentFolder().

                                                                                   {
       igetter.setCurrentFolder(dirName_); // This sets the DQMStore (should apply to ibooker as well
       runClient(ibooker, igetter);
     }

void HcalDigisClient::fill1D	(	std::string	name,
		double	X,
		double	weight = `1`
	)

inlineprivate

Definition at line 63 of file HcalDigisClient.h.

References msm_.

                                                              {
         msm_->find(name)->second->Fill(X, weight);
     }

void HcalDigisClient::fill2D	(	std::string	name,
		double	X,
		double	Y,
		double	weight = `1`
	)

inlineprivate

Definition at line 71 of file HcalDigisClient.h.

References msm_.

                                                                        {
         msm_->find(name)->second->Fill(X, Y, weight);
     }

void HcalDigisClient::fillPf	(	std::string	name,
		double	X,
		double	Y
	)

inlineprivate

Definition at line 83 of file HcalDigisClient.h.

References integralMETH2D(), msm_, alignCSCRings::s, scaleMETH2D(), str(), AlCaHLTBitMon_QueryRunRegistry::string, and x.

                                                     {
         msm_->find(name)->second->Fill(X, Y);
     }

int HcalDigisClient::HcalDigisEndjob	(	const std::vector< MonitorElement * > &	hcalMEs,
		std::string	subdet_,
		DQMStore::IBooker &	ib
	)

private

Definition at line 63 of file HcalDigisClient.cc.

References book1D(), ws_sso_content_reader::content, egammaForCoreTracking_cff::depth, particleFlowClusterHBHE_cfi::depths, MillePedeFileConverter_cfg::e, end, ev, spr::find(), objects.autophobj::float, MonitorElement::getEntries(), MonitorElement::getTH1F(), mps_fire::i, msm_, position, MonitorElement::setBinContent(), MonitorElement::setBinError(), mathSSE::sqrt(), str(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by runClient().

                                                                                                                       {
 
     using namespace std;
     string strtmp;
 
 
     MonitorElement * nevtot(nullptr);
 
     std::vector<MonitorElement*> ieta_iphi_occupancy_maps;
     std::vector<std::string> depthID;
 
     // std::cout << " Number of histos " <<     hcalMEs.size() << std::endl;
 
     for (unsigned int ih = 0; ih < hcalMEs.size(); ih++) {
       
          if (hcalMEs[ih]->getName() == "nevtot"){
               nevtot = hcalMEs[ih];
               continue;
          }
 
          //We search the occupancy maps corresponding to this subdetector
          if ( (hcalMEs[ih]->getName().find("HcalDigiTask_ieta_iphi_occupancy_map_depth") != std::string::npos)
             &&(hcalMEs[ih]->getName().find(subdet_) != std::string::npos) ){
 
                     ieta_iphi_occupancy_maps.push_back(hcalMEs[ih]);
 
                     std::string start = "depth";
                     std::string end = "_H";
 
                     int position = hcalMEs[ih]->getName().find(start) + start.length();
                     int length = hcalMEs[ih]->getName().find(end) - position;
 
                     depthID.push_back(hcalMEs[ih]->getName().substr(position, length));
 
                     continue;
          }
 
 
     }
 
 
     if ( hcalMEs.empty()){
       edm::LogError("HcalDigisClient") << "No nevtot or maps histo found..."; 
       return 0;
     }
     if( !nevtot){
       edm::LogError("HcalDigisClient") << "No nevtot histoo found...";
       return 0;
     }
     if( ieta_iphi_occupancy_maps.empty()){
       edm::LogError("HcalDigisClient") << "No maps histos found...";
       return 0;
     }
 
 
     int ev = nevtot->getEntries();
 
 
     if(ev <= 0) {
       edm::LogError("HcalDigisClient") << "normalization factor <= 0!"; 
       return 0;
     }
 
     float fev = (float) nevtot->getEntries();
 
     int depths = ieta_iphi_occupancy_maps.size();
 
     HistLim ietaLim(85, -42.5, 42.5);
 
     for (int depth = 1; depth <= depths; depth++) {
         strtmp = "HcalDigiTask_occupancy_vs_ieta_depth" + str(depth) + "_" + subdet_;
         book1D(ib,strtmp, ietaLim);
     }
 
     std::vector<float> sumphi(depths,0);
     std::vector<float> sumphie(depths,0);
 
 
     float phi_factor;
     float cnorm;
     float enorm;
 
 
     for(int depth = 1; depth <= depths; depth++){
        int nx = ieta_iphi_occupancy_maps[depth-1]->getNbinsX();
        int ny = ieta_iphi_occupancy_maps[depth-1]->getNbinsY();
 
     
        for (int i = 1; i <= nx; i++) {
            for (int j = 1; j <= ny; j++) {
 
                // occupancies
                cnorm = ieta_iphi_occupancy_maps[depth-1]->getBinContent(i, j) / fev;
                enorm = ieta_iphi_occupancy_maps[depth-1]->getBinError(i, j) / fev;
                ieta_iphi_occupancy_maps[depth-1]->setBinContent(i, j, cnorm);
                ieta_iphi_occupancy_maps[depth-1]->setBinError(i, j, enorm);
 
            } //for loop over NbinsYU
        } //for loop over NbinsX     
     } //for loop over the occupancy maps
 
     for (int i = 1; i <= 82; i++) {
 
         int ieta = i - 42; // -41 -1, 0 40
         if (ieta >= 0) ieta += 1; // -41 -1, 1 41  - to make it detector-like
 
         if (ieta >= -20 && ieta <= 20) {
             phi_factor = 72.; 
         } else {
             if (ieta >= 40 || ieta <= -40) 
                 phi_factor = 18.; 
             else 
                 phi_factor = 36.; 
         }
     
         //zero the sumphi and sumphie vector at the start of each ieta ring
         sumphi.assign(depths,0);
         sumphie.assign(depths,0);
 
         for (int iphi = 1; iphi <= 72; iphi++) {
             for(int depth = 1; depth <= depths; depth++){
                 int binIeta = ieta_iphi_occupancy_maps[depth-1]->getTH2F()->GetXaxis()->FindBin(ieta);
                 int binIphi = ieta_iphi_occupancy_maps[depth-1]->getTH2F()->GetYaxis()->FindBin(iphi);
                 
                 float content = ieta_iphi_occupancy_maps[depth-1]->getBinContent(binIeta,binIphi);
                 float econtent = ieta_iphi_occupancy_maps[depth-1]->getBinError(binIeta,binIphi);
 
                 sumphi[depth-1] += content;
                 sumphie[depth-1] += econtent*econtent;
 
             }//for loop over depths
         }//for loop over phi
  
         //double deta = double(ieta);
 
         // occupancies vs ieta
         for(int depth = 1; depth <= depths; depth++){
             strtmp = "HcalDigiTask_occupancy_vs_ieta_depth" + depthID[depth-1] + "_" + subdet_;
             MonitorElement* ME = msm_->find(strtmp)->second;
             int ietabin = ME->getTH1F()->GetXaxis()->FindBin(float(ieta));
 
             if (sumphi[depth-1]>1.e-30){
                 cnorm = sumphi[depth-1] / phi_factor;
                 enorm = sqrt(sumphie[depth-1]) / phi_factor;
                 ME->setBinContent(ietabin,cnorm);
                 ME->setBinError(ietabin,enorm);
             }
         }
     } // end of i-loop
 
   return 1;
 }

double HcalDigisClient::integralMETH2D	(	MonitorElement *	ME,
		int	i0,
		int	i1,
		int	j0,
		int	j1
	)

private

Definition at line 227 of file HcalDigisClient.cc.

References MonitorElement::getBinContent(), and mps_fire::i.

Referenced by fillPf().

                                                                                          {
     double sum(0);
     for (int i = i0; i <= i1; i++) {
         for (int j = j0; j <= j1; j++) {
             sum += ME->getBinContent(i, j);
         }
     }
 
     return sum;
 }

MonitorElement * HcalDigisClient::monitor ( std::string name )

private

Definition at line 216 of file HcalDigisClient.cc.

References msm_.

                                                        {
     if (!msm_->count(name)) return nullptr;
     else return msm_->find(name)->second;
 }

void HcalDigisClient::runClient	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

privatevirtual

Definition at line 37 of file HcalDigisClient.cc.

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

Referenced by dqmEndJob().

                                                                         {
     ig.setCurrentFolder(dirName_);
     std::vector<MonitorElement*> hcalMEs;
     // Since out folders are fixed to three, we can just go over these three folders
     // i.e., CaloTowersV/CaloTowersTask, HcalRecHitsV/HcalRecHitTask, NoiseRatesV/NoiseRatesTask.
     std::vector<std::string> fullPathHLTFolders = ig.getSubdirs();
     for (unsigned int i = 0; i < fullPathHLTFolders.size(); i++) {
         ig.setCurrentFolder(fullPathHLTFolders[i]);
         std::vector<std::string> fullSubPathHLTFolders = ig.getSubdirs();
         for (unsigned int j = 0; j < fullSubPathHLTFolders.size(); j++) {
             if (strcmp(fullSubPathHLTFolders[j].c_str(), "HcalDigisV/HcalDigiTask") == 0) {
                 hcalMEs = ig.getContents(fullSubPathHLTFolders[j]);
                 ig.setCurrentFolder("HcalDigisV/HcalDigiTask");
                 if (!HcalDigisEndjob(hcalMEs, "HB", ib)) 
                     edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HB"; 
                 if (!HcalDigisEndjob(hcalMEs, "HE", ib)) 
                     edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HE"; 
                 if (!HcalDigisEndjob(hcalMEs, "HO", ib)) 
                     edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HO"; 
                 if (!HcalDigisEndjob(hcalMEs, "HF", ib)) 
                     edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HF";             
             }
         }
     }
 }

void HcalDigisClient::scaleMETH2D	(	MonitorElement *	ME,
		double	s
	)

private

Definition at line 238 of file HcalDigisClient.cc.

References ws_sso_content_reader::content, DEFINE_FWK_MODULE, relativeConstraints::error, MonitorElement::getBinContent(), MonitorElement::getBinError(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), mps_fire::i, alignCSCRings::s, MonitorElement::setBinContent(), and MonitorElement::setBinError().

Referenced by fillPf().

                                                               {
     int nx = ME->getNbinsX();
     int ny = ME->getNbinsY();
 
     double content(0);
     double error(0);
     for (int i = 1; i <= nx; i++) {
         for (int j = 1; j <= ny; j++) {
             content = ME->getBinContent(i, j);
             error = ME->getBinError(i, j);
             content *= s;
             error *= s;
             ME->setBinContent(i, j, content);
             ME->setBinError(i, j, error);
         }
     }
 }

std::string HcalDigisClient::str ( int x )

private

Definition at line 221 of file HcalDigisClient.cc.

References MillePedeFileConverter_cfg::out, and x.

Referenced by fillPf(), and HcalDigisEndjob().

                                     {
     std::stringstream out;
     out << x;
     return out.str();
 }

Member Data Documentation

std::string HcalDigisClient::dirName_

private

Definition at line 93 of file HcalDigisClient.h.

Referenced by dqmEndJob(), HcalDigisClient(), and runClient().

std::map<std::string, MonitorElement*>* HcalDigisClient::msm_

private

Definition at line 91 of file HcalDigisClient.h.

Referenced by book1D(), book2D(), bookPf(), fill1D(), fill2D(), fillPf(), HcalDigisClient(), HcalDigisEndjob(), monitor(), and ~HcalDigisClient().

std::string HcalDigisClient::outputFile_

private

Definition at line 92 of file HcalDigisClient.h.

Referenced by HcalDigisClient().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation