#include <DQMRivetClient.h>

Inheritance diagram for DQMRivetClient:

Classes
struct	LumiOption

struct	NormOption

struct	ScaleFactorOption

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override

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

	DQMRivetClient (const edm::ParameterSet &pset)

void	endJob () override

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

void	normalizeToIntegral (const std::string &startDir, const std::string &histName, const std::string &normHistName)

void	normalizeToLumi (const std::string &startDir, const std::string &histName, const std::string &normHistName, double xsection)

void	scaleByFactor (const std::string &startDir, const std::string &histName, double factor)

	~DQMRivetClient () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources, edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
std::vector< LumiOption >	lumiOptions_

std::vector< NormOption >	normOptions_

std::string	outputFileName_

std::vector< ScaleFactorOption >	scaleOptions_

std::vector< std::string >	subDirs_

DQMStore *	theDQM

unsigned int	verbose_

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
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)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 12 of file DQMRivetClient.h.

Member Typedef Documentation

◆ DQMStore

typedef dqm::legacy::DQMStore DQMRivetClient::DQMStore

Definition at line 14 of file DQMRivetClient.h.

◆ MonitorElement

typedef dqm::legacy::MonitorElement DQMRivetClient::MonitorElement

Definition at line 15 of file DQMRivetClient.h.

Constructor & Destructor Documentation

◆ DQMRivetClient()

DQMRivetClient::DQMRivetClient ( const edm::ParameterSet & pset )

Definition at line 29 of file DQMRivetClient.cc.

References writedatasetfile::args, das-up-to-nevents::opt, or, muonDTDigis_cfi::pset, and HistogramManager_cfi::VPSet().

                                                        {
   typedef std::vector<edm::ParameterSet> VPSet;
   typedef std::vector<std::string> vstring;
   typedef boost::escaped_list_separator<char> elsc;
 
   usesResource("DQMStore");
 
   elsc commonEscapes("\\", " \t", "\'");
 
   // Parse Normalization commands
   vstring normCmds = pset.getUntrackedParameter<vstring>("normalizationToIntegral", vstring());
   for (vstring::const_iterator normCmd = normCmds.begin(); normCmd != normCmds.end(); ++normCmd) {
     if (normCmd->empty())
       continue;
     boost::tokenizer<elsc> tokens(*normCmd, commonEscapes);
 
     vector<string> args;
     for (boost::tokenizer<elsc>::const_iterator iToken = tokens.begin(); iToken != tokens.end(); ++iToken) {
       if (iToken->empty())
         continue;
       args.push_back(*iToken);
     }
 
     if (args.empty() or args.size() > 2) {
       LogInfo("DQMRivetClient") << "Wrong input to normCmds\n";
       continue;
     }
 
     NormOption opt;
     opt.name = args[0];
     opt.normHistName = args.size() == 2 ? args[1] : args[0];
 
     normOptions_.push_back(opt);
   }
 
   VPSet normSets = pset.getUntrackedParameter<VPSet>("normalizationToIntegralSets", VPSet());
   for (VPSet::const_iterator normSet = normSets.begin(); normSet != normSets.end(); ++normSet) {
     NormOption opt;
     opt.name = normSet->getUntrackedParameter<string>("name");
     opt.normHistName = normSet->getUntrackedParameter<string>("normalizedTo", opt.name);
 
     normOptions_.push_back(opt);
   }
 
   //normalize to lumi
   vstring lumiCmds = pset.getUntrackedParameter<vstring>("normalizationToLumi", vstring());
   for (vstring::const_iterator lumiCmd = lumiCmds.begin(); lumiCmd != lumiCmds.end(); ++lumiCmd) {
     if (lumiCmd->empty())
       continue;
     boost::tokenizer<elsc> tokens(*lumiCmd, commonEscapes);
 
     vector<string> args;
     for (boost::tokenizer<elsc>::const_iterator iToken = tokens.begin(); iToken != tokens.end(); ++iToken) {
       if (iToken->empty())
         continue;
       args.push_back(*iToken);
     }
 
     if (args.size() != 2) {
       LogInfo("DQMRivetClient") << "Wrong input to lumiCmds\n";
       continue;
     }
 
     DQMRivetClient::LumiOption opt;
     opt.name = args[0];
     opt.normHistName = args[1];
     opt.xsection = pset.getUntrackedParameter<double>("xsection", -1.);
     //opt.xsection = atof(args[2].c_str());
 
     //std::cout << opt.name << " " << opt.normHistName << " " << opt.xsection << std::endl;
     lumiOptions_.push_back(opt);
   }
 
   //multiply by a number
   vstring scaleCmds = pset.getUntrackedParameter<vstring>("scaleBy", vstring());
   for (vstring::const_iterator scaleCmd = scaleCmds.begin(); scaleCmd != scaleCmds.end(); ++scaleCmd) {
     if (scaleCmd->empty())
       continue;
     boost::tokenizer<elsc> tokens(*scaleCmd, commonEscapes);
 
     vector<string> args;
     for (boost::tokenizer<elsc>::const_iterator iToken = tokens.begin(); iToken != tokens.end(); ++iToken) {
       if (iToken->empty())
         continue;
       args.push_back(*iToken);
     }
 
     if (args.empty() or args.size() > 2) {
       LogInfo("DQMRivetClient") << "Wrong input to normCmds\n";
       continue;
     }
 
     ScaleFactorOption opt;
     opt.name = args[0];
     opt.scale = atof(args[1].c_str());
     scaleOptions_.push_back(opt);
   }
 
   outputFileName_ = pset.getUntrackedParameter<string>("outputFileName", "");
   subDirs_ = pset.getUntrackedParameter<vstring>("subDirs");
 }

◆ ~DQMRivetClient()

DQMRivetClient::~DQMRivetClient ( )

inlineoverride

Definition at line 21 of file DQMRivetClient.h.

21 {};

Member Function Documentation

◆ analyze()

void DQMRivetClient::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

inlineoverridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 23 of file DQMRivetClient.h.

23 {};

◆ beginRun()

void DQMRivetClient::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

inlineoverride

Definition at line 27 of file DQMRivetClient.h.

27 {}

◆ endJob()

void DQMRivetClient::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 191 of file DQMRivetClient.cc.

References LogTrace.

                             {
   // Update 2009-09-23
   // Migrated all code from here to endRun
 
   LogTrace("DQMRivetClient") << "inside of ::endJob()" << endl;
 }

◆ endRun()

void DQMRivetClient::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

override

Definition at line 131 of file DQMRivetClient.cc.

References TrackerOfflineValidation_Dqm_cff::dirName, PostProcessor_cff::normalizeToIntegral, PostProcessor_cff::normalizeToLumi, Utilities::operator, and PostProcessorHcalIsoTrack_cfi::subDir.

                                                                    {
   typedef vector<string> vstring;
 
   // Update 2009-09-23
   // Migrated all code from endJob to this function
   // endJob is not necessarily called in the proper sequence
   // and does not necessarily book histograms produced in
   // that step.
   // It more robust to do the histogram manipulation in
   // this endRun function
 
   theDQM = nullptr;
   theDQM = Service<DQMStore>().operator->();
 
   if (!theDQM) {
     LogInfo("DQMRivetClient") << "Cannot create DQMStore instance\n";
     return;
   }
 
   // Process wildcard in the sub-directory
   set<string> subDirSet;
 
   for (vstring::const_iterator iSubDir = subDirs_.begin(); iSubDir != subDirs_.end(); ++iSubDir) {
     string subDir = *iSubDir;
 
     if (subDir[subDir.size() - 1] == '/')
       subDir.erase(subDir.size() - 1);
 
     subDirSet.insert(subDir);
   }
 
   for (set<string>::const_iterator iSubDir = subDirSet.begin(); iSubDir != subDirSet.end(); ++iSubDir) {
     const string& dirName = *iSubDir;
     for (vector<NormOption>::const_iterator normOption = normOptions_.begin(); normOption != normOptions_.end();
          ++normOption) {
       normalizeToIntegral(dirName, normOption->name, normOption->normHistName);
     }
   }
 
   for (set<string>::const_iterator iSubDir = subDirSet.begin(); iSubDir != subDirSet.end(); ++iSubDir) {
     const string& dirName = *iSubDir;
     for (vector<LumiOption>::const_iterator lumiOption = lumiOptions_.begin(); lumiOption != lumiOptions_.end();
          ++lumiOption) {
       normalizeToLumi(dirName, lumiOption->name, lumiOption->normHistName, lumiOption->xsection);
     }
   }
 
   for (set<string>::const_iterator iSubDir = subDirSet.begin(); iSubDir != subDirSet.end(); ++iSubDir) {
     const string& dirName = *iSubDir;
     for (vector<ScaleFactorOption>::const_iterator scaleOption = scaleOptions_.begin();
          scaleOption != scaleOptions_.end();
          ++scaleOption) {
       scaleByFactor(dirName, scaleOption->name, scaleOption->scale);
     }
   }
 
   if (!outputFileName_.empty())
     theDQM->save(outputFileName_);
 }

◆ normalizeToIntegral()

void DQMRivetClient::normalizeToIntegral	(	const std::string &	startDir,
		const std::string &	histName,
		const std::string &	normHistName
	)

Definition at line 198 of file DQMRivetClient.cc.

References compareTotals::hist.

                                                                         {
   if (!theDQM->dirExists(startDir)) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Cannot find sub-directory " << startDir << endl;
     return;
   }
 
   theDQM->cd();
 
   ME* element = theDQM->get(startDir + "/" + histName);
   ME* normME = theDQM->get(startDir + "/" + normHistName);
 
   if (!element) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "No such element '" << histName << "' found\n";
     return;
   }
 
   if (!normME) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "No such element '" << normHistName << "' found\n";
     return;
   }
 
   TH1F* hist = element->getTH1F();
   if (!hist) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Cannot create TH1F from ME\n";
     return;
   }
 
   TH1F* normHist = normME->getTH1F();
   if (!normHist) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Cannot create TH1F from ME\n";
     return;
   }
 
   const double entries = normHist->Integral();
   if (entries != 0) {
     hist->Scale(1. / entries, "width");
   } else {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Zero entries in histogram\n";
   }
 
   return;
 }

◆ normalizeToLumi()

void DQMRivetClient::normalizeToLumi	(	const std::string &	startDir,
		const std::string &	histName,
		const std::string &	normHistName,
		double	xsection
	)

Definition at line 249 of file DQMRivetClient.cc.

References compareTotals::hist, PostProcessor_cff::normalizeToIntegral, and PostProcessor_cff::xsection.

                                                       {
   normalizeToIntegral(startDir, histName, normHistName);
   theDQM->cd();
   ME* element = theDQM->get(startDir + "/" + histName);
   TH1F* hist = element->getTH1F();
   if (!hist) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Cannot create TH1F from ME\n";
     return;
   }
   hist->Scale(xsection);
   return;
 }

◆ scaleByFactor()

void DQMRivetClient::scaleByFactor	(	const std::string &	startDir,
		const std::string &	histName,
		double	factor
	)

Definition at line 266 of file DQMRivetClient.cc.

References compareTotals::hist.

                                                                                                       {
   if (!theDQM->dirExists(startDir)) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Cannot find sub-directory " << startDir << endl;
     return;
   }
 
   theDQM->cd();
 
   ME* element = theDQM->get(startDir + "/" + histName);
 
   if (!element) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "No such element '" << histName << "' found\n";
     return;
   }
 
   TH1F* hist = element->getTH1F();
   if (!hist) {
     LogInfo("DQMRivetClient") << "normalizeToEntries() : "
                               << "Cannot create TH1F from ME\n";
     return;
   }
   hist->Scale(factor);
 }

Member Data Documentation

◆ lumiOptions_

std::vector<LumiOption> DQMRivetClient::lumiOptions_

private

Definition at line 55 of file DQMRivetClient.h.

◆ normOptions_

std::vector<NormOption> DQMRivetClient::normOptions_

private

Definition at line 54 of file DQMRivetClient.h.

◆ outputFileName_

std::string DQMRivetClient::outputFileName_

private

Definition at line 52 of file DQMRivetClient.h.

◆ scaleOptions_

std::vector<ScaleFactorOption> DQMRivetClient::scaleOptions_

private

Definition at line 56 of file DQMRivetClient.h.

◆ subDirs_

std::vector<std::string> DQMRivetClient::subDirs_

private

Definition at line 51 of file DQMRivetClient.h.

◆ theDQM

DQMStore* DQMRivetClient::theDQM

private

Definition at line 50 of file DQMRivetClient.h.

◆ verbose_

unsigned int DQMRivetClient::verbose_

private

Definition at line 48 of file DQMRivetClient.h.

Classes

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DQMStore

◆ MonitorElement

Constructor & Destructor Documentation

◆ DQMRivetClient()

◆ ~DQMRivetClient()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ normalizeToIntegral()

◆ normalizeToLumi()

◆ scaleByFactor()

Member Data Documentation

◆ lumiOptions_

◆ normOptions_

◆ outputFileName_

◆ scaleOptions_

◆ subDirs_

◆ theDQM

◆ verbose_