#include <CaloTowersDQMClient.h>

Inheritance diagram for CaloTowersDQMClient:

Public Member Functions
void	beginJob (void) override

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

	CaloTowersDQMClient (const edm::ParameterSet &)

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

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

	~CaloTowersDQMClient () 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 Attributes
edm::ParameterSet	conf_

bool	debug_

std::string	dirName_

std::string	dirNameJet_

std::string	dirNameMET_

std::string	outputFile_

bool	verbose_

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 38 of file CaloTowersDQMClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 12 of file CaloTowersDQMClient.cc.

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

                                                                        : conf_(iConfig) {
   outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile", "myfile.root");
   debug_ = false;
   verbose_ = false;
   dirName_ = iConfig.getParameter<std::string>("DQMDirName");
 }

CaloTowersDQMClient::~CaloTowersDQMClient ( )

override

Definition at line 19 of file CaloTowersDQMClient.cc.

19 {}

Member Function Documentation

void CaloTowersDQMClient::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 21 of file CaloTowersDQMClient.cc.

21 {}

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

override

Definition at line 23 of file CaloTowersDQMClient.cc.

23 {}

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

Definition at line 27 of file CaloTowersDQMClient.cc.

References gather_cfg::cout, MillePedeFileConverter_cfg::e, objects.autophobj::float, MonitorElement::getBinContent(), MonitorElement::getBinError(), MonitorElement::getEntries(), getName(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), mps_fire::i, nevent, funct::pow(), MonitorElement::setBinContent(), MonitorElement::setBinError(), and verbose_.

Referenced by dqmEndJob().

                                                                                     {
   int useAllHistos = 0;
   MonitorElement *Ntowers_vs_ieta = nullptr;
   MonitorElement *mapEnergy_N = nullptr, *mapEnergy_E = nullptr, *mapEnergy_H = nullptr, *mapEnergy_EH = nullptr;
   MonitorElement *occupancy_map = nullptr, *occupancy_vs_ieta = nullptr;
   for (unsigned int ih = 0; ih < hcalMEs.size(); ih++) {
     if (strcmp(hcalMEs[ih]->getName().c_str(), "Ntowers_per_event_vs_ieta") == 0) {
       Ntowers_vs_ieta = hcalMEs[ih];
     }
     if (strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_Nentries") == 0) {
       mapEnergy_N = hcalMEs[ih];
     }
     if (strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_H") == 0) {
       useAllHistos++;
       mapEnergy_H = hcalMEs[ih];
     }
     if (strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_E") == 0) {
       useAllHistos++;
       mapEnergy_E = hcalMEs[ih];
     }
     if (strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_EH") == 0) {
       useAllHistos++;
       mapEnergy_EH = hcalMEs[ih];
     }
     if (strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_occupancy") == 0) {
       occupancy_map = hcalMEs[ih];
     }
     if (strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_occupancy_vs_ieta") == 0) {
       occupancy_vs_ieta = hcalMEs[ih];
     }
   }
   if (useAllHistos != 0 && useAllHistos != 3)
     return 0;
 
   double nevent = mapEnergy_N->getEntries();
   if (verbose_)
     std::cout << "nevent : " << nevent << std::endl;
 
   // mean number of towers per ieta
   int nx = Ntowers_vs_ieta->getNbinsX();
   float cont;
   float conte;
   float fev = float(nevent);
 
   for (int i = 1; i <= nx; i++) {
     cont = Ntowers_vs_ieta->getBinContent(i) / fev;
     conte = pow(Ntowers_vs_ieta->getBinContent(i), 0.5) / fev;
     Ntowers_vs_ieta->setBinContent(i, cont);
     Ntowers_vs_ieta->setBinError(i, conte);
   }
 
   // mean energies & occupancies evaluation
 
   nx = mapEnergy_N->getNbinsX();
   int ny = mapEnergy_N->getNbinsY();
   float cnorm;
   float cnorme;
   float phi_factor;
 
   for (int i = 1; i <= nx; i++) {
     float sumphi = 0.;
 
     for (int j = 1; j <= ny; j++) {
       // Emean
       cnorm = mapEnergy_N->getBinContent(i, j);
       // Phi histos are not used in the macros
       if (cnorm > 0.000001 && useAllHistos) {
         cont = mapEnergy_E->getBinContent(i, j) / cnorm;
         conte = mapEnergy_E->getBinError(i, j) / cnorm;
         mapEnergy_E->setBinContent(i, j, cont);
         mapEnergy_E->setBinError(i, j, conte);
 
         cont = mapEnergy_H->getBinContent(i, j) / cnorm;
         conte = mapEnergy_H->getBinError(i, j) / cnorm;
         mapEnergy_H->setBinContent(i, j, cont);
         mapEnergy_H->setBinError(i, j, conte);
 
         cont = mapEnergy_EH->getBinContent(i, j) / cnorm;
         conte = mapEnergy_EH->getBinError(i, j) / cnorm;
         mapEnergy_EH->setBinContent(i, j, cont);
         mapEnergy_EH->setBinError(i, j, conte);
       }
 
       // Occupancy (needed for occupancy vs ieta)
       cont = occupancy_map->getBinContent(i, j);
       conte = occupancy_map->getBinError(i, j);
       if (fev > 0. && cnorm > 1.e-30) {
         occupancy_map->setBinContent(i, j, cont / fev);
         occupancy_map->setBinError(i, j, conte / fev);
       }
 
       sumphi += cont;
 
     }  // end of iphy cycle (j)
 
     // Occupancy vs ieta histo is drawn
     // phi-factor evaluation for occupancy_vs_ieta calculation
     int ieta = i - 43;  // should be the same as int ieta =
                         // int(occupancy_vs_ieta->getBinCenter(i));
 
     if (ieta >= -20 && ieta <= 20) {
       phi_factor = 72.;
     } else {
       if (ieta >= 40 || ieta <= -40) {
         phi_factor = 18.;
       } else
         phi_factor = 36.;
     }
 
     cnorm = sumphi / phi_factor;
     cnorme = pow(sumphi, 0.5) / phi_factor;
     if (fev > 0. && cnorm > 1.e-30) {
       occupancy_vs_ieta->setBinContent(i, cnorm / fev);
       occupancy_vs_ieta->setBinError(i, cnorme / fev);
     }
 
   }  // end of ieta cycle (i)
 
   return 1;
 }

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

overridevirtual

Implements DQMEDHarvester.

Definition at line 148 of file CaloTowersDQMClient.cc.

References CaloTowersEndjob(), gather_cfg::cout, DEFINE_FWK_MODULE, dirName_, DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), mps_fire::i, DQMStore::IGetter::setCurrentFolder(), and verbose_.

                                                                                       {
   igetter.setCurrentFolder(dirName_);
   if (verbose_)
     std::cout << "\nrunClient" << std::endl;
 
   std::vector<MonitorElement *> hcalMEs;
 
   // Since out folders are fixed to three, we can just go over these three
   // folders i.e., CaloTowersD/CaloTowersTask, HcalRecHitsD/HcalRecHitTask,
   // NoiseRatesV/NoiseRatesTask.
   std::vector<std::string> fullPathHLTFolders = igetter.getSubdirs();
   for (unsigned int i = 0; i < fullPathHLTFolders.size(); i++) {
     if (verbose_)
       std::cout << "\nfullPath: " << fullPathHLTFolders[i] << std::endl;
     igetter.setCurrentFolder(fullPathHLTFolders[i]);
 
     std::vector<std::string> fullSubPathHLTFolders = igetter.getSubdirs();
     for (unsigned int j = 0; j < fullSubPathHLTFolders.size(); j++) {
       if (verbose_)
         std::cout << "fullSub: " << fullSubPathHLTFolders[j] << std::endl;
 
       if (strcmp(fullSubPathHLTFolders[j].c_str(), "CaloTowersD/CaloTowersTask") == 0) {
         hcalMEs = igetter.getContents(fullSubPathHLTFolders[j]);
         if (verbose_)
           std::cout << "hltMES size : " << hcalMEs.size() << std::endl;
         if (!CaloTowersEndjob(hcalMEs))
           std::cout << "\nError in CaloTowersEndjob!" << std::endl << std::endl;
       }
     }
   }
 }

Member Data Documentation

edm::ParameterSet CaloTowersDQMClient::conf_

private

Definition at line 41 of file CaloTowersDQMClient.h.

bool CaloTowersDQMClient::debug_

private

Definition at line 44 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient().

std::string CaloTowersDQMClient::dirName_

private

Definition at line 46 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient(), and dqmEndJob().

std::string CaloTowersDQMClient::dirNameJet_

private

Definition at line 47 of file CaloTowersDQMClient.h.

std::string CaloTowersDQMClient::dirNameMET_

private

Definition at line 48 of file CaloTowersDQMClient.h.

std::string CaloTowersDQMClient::outputFile_

private

Definition at line 40 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient().

bool CaloTowersDQMClient::verbose_

private

Definition at line 43 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient(), CaloTowersEndjob(), and dqmEndJob().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation