Inheritance diagram for HLTGenValClient:

Classes
struct	EfficOption

Public Member Functions
void	computeEfficiency (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string &dirName, const std::string &efficMEName, const std::string &efficMETitle, const std::string &numeratorMEName, const std::string &denominatorMEName)

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

void	dqmEndLuminosityBlock (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override

void	dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &) override

	HLTGenValClient (const edm::ParameterSet &pset)

	~HLTGenValClient () override

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginJob () override

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &) final

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

virtual void	dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)

	DQMEDHarvester (edm::ParameterSet const &iConfig)

	DQMEDHarvester ()

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final

void	endProcessBlockProduce (edm::ProcessBlock &) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &es) final

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	findAllSubdirectories (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, std::string dir, std::set< std::string > *myList, const TString &pattern)

void	genericEff (TH1 denom, TH1 numer, MonitorElement *efficiencyHist)

void	makeAllPlots (DQMStore::IBooker &, DQMStore::IGetter &)

Private Attributes
std::vector< EfficOption >	efficOptions_

bool	isWildcardUsed_

bool	makeGlobalEffPlot_

TPRegexp	metacharacters_

TPRegexp	nonPerlWildcard_

std::string	outputFileName_

bool	runOnEndJob_

bool	runOnEndLumi_

const std::string	separator_ = "__"

std::vector< std::string >	subDirs_

DQMStore *	theDQM

unsigned int	verbose_

Additional Inherited Members
Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

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

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

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

Protected Attributes inherited from DQMEDHarvester
DQMStore *	dqmstore_

edm::GetterOfProducts< DQMToken >	jobmegetter_

edm::EDPutTokenT< DQMToken >	jobToken_

edm::GetterOfProducts< DQMToken >	lumimegetter_

edm::EDPutTokenT< DQMToken >	lumiToken_

edm::GetterOfProducts< DQMToken >	runmegetter_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 36 of file HLTGenValClient.cc.

Constructor & Destructor Documentation

◆ HLTGenValClient()

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

Definition at line 89 of file HLTGenValClient.cc.

References isWildcardUsed_, makeGlobalEffPlot_, outputFileName_, muonDTDigis_cfi::pset, runOnEndJob_, runOnEndLumi_, AlCaHLTBitMon_QueryRunRegistry::string, subDirs_, and verbose_.

     : metacharacters_("[\\^\\$\\.\\*\\+\\?\\|\\(\\)\\{\\}\\[\\]]"), nonPerlWildcard_("\\w\\*|^\\*") {
   boost::escaped_list_separator<char> commonEscapes("\\", " \t", "\'");
 
   verbose_ = pset.getUntrackedParameter<unsigned int>("verbose", 0);
   runOnEndLumi_ = pset.getUntrackedParameter<bool>("runOnEndLumi", false);
   runOnEndJob_ = pset.getUntrackedParameter<bool>("runOnEndJob", true);
   makeGlobalEffPlot_ = pset.getUntrackedParameter<bool>("makeGlobalEffienciesPlot", true);
 
   outputFileName_ = pset.getUntrackedParameter<std::string>("outputFileName", "");
   subDirs_ = pset.getUntrackedParameter<std::vector<std::string>>("subDirs");
 
   isWildcardUsed_ = false;
 }

◆ ~HLTGenValClient()

HLTGenValClient::~HLTGenValClient ( )

inlineoverride

Definition at line 39 of file HLTGenValClient.cc.

39 {};

Member Function Documentation

◆ computeEfficiency()

void HLTGenValClient::computeEfficiency	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const std::string &	dirName,
		const std::string &	efficMEName,
		const std::string &	efficMETitle,
		const std::string &	numeratorMEName,
		const std::string &	denominatorMEName
	)

Definition at line 244 of file HLTGenValClient.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::book3D(), dqm::implementation::NavigatorBase::cd(), dqm::implementation::IGetter::dirExists(), TrackerOfflineValidation_Dqm_cff::dirName, genericEff(), dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getEntries(), dqm::legacy::MonitorElement::getTH1(), dqm::legacy::MonitorElement::getTH1F(), isWildcardUsed_, makeGlobalEffPlot_, dqm::implementation::NavigatorBase::setCurrentFolder(), dqm::impl::MonitorElement::setEfficiencyFlag(), dqm::impl::MonitorElement::setEntries(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

Referenced by makeAllPlots().

                                                                             {
   // checking if directory exists
   if (!igetter.dirExists(dirName)) {
     if (verbose_ >= 2 || (verbose_ == 1 && !isWildcardUsed_)) {
       edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                        << "Cannot find sub-directory " << dirName << std::endl;
     }
     return;
   }
 
   ibooker.cd();
 
   // getting input MEs
   HLTGenValClient::MonitorElement* denominatorME = igetter.get(dirName + "/" + denominatorMEName);
   HLTGenValClient::MonitorElement* numeratorME = igetter.get(dirName + "/" + numeratorMEName);
 
   // checking of input MEs exist
   if (!denominatorME) {
     if (verbose_ >= 2 || (verbose_ == 1 && !isWildcardUsed_)) {
       edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                        << "No denominator-ME '" << denominatorMEName << "' found\n";
     }
     return;
   }
   if (!numeratorME) {
     if (verbose_ >= 2 || (verbose_ == 1 && !isWildcardUsed_)) {
       edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                        << "No numerator-ME '" << numeratorMEName << "' found\n";
     }
     return;
   }
 
   // Treat everything as the base class, TH1
   TH1* hDenominator = denominatorME->getTH1();
   TH1* hNumerator = numeratorME->getTH1();
 
   // check if TH1 extraction has succeeded
   if (!hDenominator || !hNumerator) {
     if (verbose_ >= 2 || (verbose_ == 1 && !isWildcardUsed_)) {
       edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                        << "Cannot create TH1 from ME\n";
     }
     return;
   }
 
   // preparing efficiency output path and name
   std::string efficDir = dirName;
   std::string newEfficMEName = efficMEName;
   std::string::size_type shiftPos;
   if (std::string::npos != (shiftPos = efficMEName.rfind('/'))) {
     efficDir += "/" + efficMEName.substr(0, shiftPos);
     newEfficMEName.erase(0, shiftPos + 1);
   }
   ibooker.setCurrentFolder(efficDir);
 
   // creating the efficiency MonitorElement
   HLTGenValClient::MonitorElement* efficME = nullptr;
 
   // We need to know what kind of TH1 we have
   // That information is obtained from the class name of the hDenominator
   // Then we use the appropriate booking function
   TH1* efficHist = static_cast<TH1*>(hDenominator->Clone(newEfficMEName.c_str()));
   efficHist->SetDirectory(nullptr);
   efficHist->SetTitle(efficMETitle.c_str());
   TClass* myHistClass = efficHist->IsA();
   std::string histClassName = myHistClass->GetName();
   if (histClassName == "TH1F") {
     efficME = ibooker.book1D(newEfficMEName, (TH1F*)efficHist);
   } else if (histClassName == "TH2F") {
     efficME = ibooker.book2D(newEfficMEName, (TH2F*)efficHist);
   } else if (histClassName == "TH3F") {
     efficME = ibooker.book3D(newEfficMEName, (TH3F*)efficHist);
   }
   delete efficHist;
 
   // checking whether efficME was succesfully created
   if (!efficME) {
     edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                      << "Cannot book effic-ME from the DQM\n";
     return;
   }
 
   // actually calculating the efficiency and filling the ME
   genericEff(hDenominator, hNumerator, efficME);
   efficME->setEntries(denominatorME->getEntries());
 
   // Putting total efficiency in "GLobal efficiencies" histogram
   if (makeGlobalEffPlot_) {
     // getting global efficiency ME
     HLTGenValClient::MonitorElement* globalEfficME = igetter.get(efficDir + "/globalEfficiencies");
     if (!globalEfficME)  // in case it does not exist yet, we create it
       globalEfficME = ibooker.book1D("globalEfficiencies", "Global efficiencies", 1, 0, 1);
     if (!globalEfficME) {  // error handling in case creation failed
       edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                        << "Cannot book globalEffic-ME from the DQM\n";
       return;
     }
     globalEfficME->setEfficiencyFlag();
 
     // extracting histogram
     TH1F* hGlobalEffic = globalEfficME->getTH1F();
     if (!hGlobalEffic) {
       edm::LogError("HLTGenValClient") << "computeEfficiency() : "
                                        << "Cannot create TH1F from ME, globalEfficME\n";
       return;
     }
 
     // getting total counts
     const float nDenominatorAll = hDenominator->GetEntries();
     const float nNumeratorAll = hNumerator->GetEntries();
 
     // calculating total efficiency
     float efficAll = 0;
     float errorAll = 0;
     efficAll = nDenominatorAll ? nNumeratorAll / nDenominatorAll : 0;
     errorAll = nDenominatorAll && efficAll < 1 ? sqrt(efficAll * (1 - efficAll) / nDenominatorAll) : 0;
 
     // Filling the histogram bin
     const int iBin = hGlobalEffic->Fill(newEfficMEName.c_str(), 0);
     hGlobalEffic->SetBinContent(iBin, efficAll);
     hGlobalEffic->SetBinError(iBin, errorAll);
   }
 }

◆ dqmEndJob()

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

inlineoverridevirtual

Implements DQMEDHarvester.

Definition at line 46 of file HLTGenValClient.cc.

46 {};

◆ dqmEndLuminosityBlock()

void HLTGenValClient::dqmEndLuminosityBlock	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	c
	)

overridevirtual

Reimplemented from DQMEDHarvester.

Definition at line 104 of file HLTGenValClient.cc.

References makeAllPlots(), and runOnEndLumi_.

                                                                     {
   if (runOnEndLumi_) {
     makeAllPlots(ibooker, igetter);
   }
 }

◆ dqmEndRun()

void HLTGenValClient::dqmEndRun	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from DQMEDHarvester.

Definition at line 113 of file HLTGenValClient.cc.

References makeAllPlots(), Utilities::operator, outputFileName_, runOnEndJob_, dqm::implementation::DQMStore::save(), and theDQM.

                                                       {
   // Create new MEs in endRun, even though we are requested to do it in endJob.
   // This gives the QTests a chance to run, before summaries are created in
   // endJob. The negative side effect is that we cannot run the GenericClient
   // for plots produced in Harvesting, but that seems rather rare.
   //
   // It is important that this is still save in the presence of multiple runs,
   // first because in multi-run harvesting, we accumulate statistics over all
   // runs and have full statistics at the endRun of the last run, and second,
   // because we set the efficiencyFlag so any further aggregation should produce
   // correct results. Also, all operations should be idempotent; running them
   // more than once does no harm.
 
   theDQM = edm::Service<DQMStore>().operator->();
 
   if (runOnEndJob_) {
     makeAllPlots(ibooker, igetter);
   }
 
   if (!outputFileName_.empty())
     theDQM->save(outputFileName_);
 }

◆ findAllSubdirectories()

void HLTGenValClient::findAllSubdirectories	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		std::string	dir,
		std::set< std::string > *	myList,
		const TString &	pattern = `TString("")`
	)

private

Definition at line 376 of file HLTGenValClient.cc.

References dqm::implementation::NavigatorBase::cd(), DeadROC_duringRun::dir, dqm::implementation::IGetter::dirExists(), TrackerOfflineValidation_Dqm_cff::dirName, dqm::implementation::IGetter::getSubdirs(), PatBasicFWLiteJetAnalyzer_Selector_cfg::myList, nonPerlWildcard_, and fetchall_from_DQM_v2::regexp.

Referenced by makeAllPlots().

                                                                                    {
   TString patternTmp = _pattern;
 
   // checking if directory exists
   if (!igetter.dirExists(dir)) {
     edm::LogError("HLTGenValClient") << " HLTGenValClient::findAllSubdirectories ==> Missing folder " << dir << " !!!";
     return;
   }
 
   // replacing wildcards
   if (patternTmp != "") {
     if (patternTmp.Contains(nonPerlWildcard_))
       patternTmp.ReplaceAll("*", ".*");
     TPRegexp regexp(patternTmp);
     ibooker.cd(dir);
     std::vector<std::string> foundDirs = igetter.getSubdirs();
     for (const auto& iDir : foundDirs) {
       TString dirName = iDir.substr(iDir.rfind('/') + 1, iDir.length());
       if (dirName.Contains(regexp))
         findAllSubdirectories(ibooker, igetter, iDir, myList);
     }
   } else if (igetter.dirExists(dir)) {
     // we have found a subdirectory - adding it to the list
     myList->insert(dir);
 
     // moving into the found subdirectory and recursively continue
     ibooker.cd(dir);
     findAllSubdirectories(ibooker, igetter, dir, myList, "*");
 
   } else {
     // error handling in case found directory does not exist
     edm::LogError("HLTGenValClient") << "Trying to find sub-directories of " << dir << " failed because " << dir
                                      << " does not exist";
   }
   return;
 }

◆ genericEff()

void HLTGenValClient::genericEff	(	TH1 *	denom,
		TH1 *	numer,
		MonitorElement *	efficiencyHist
	)

private

Definition at line 418 of file HLTGenValClient.cc.

References makePileupJSON::denom, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinError(), dqm::impl::MonitorElement::setEfficiencyFlag(), and mathSSE::sqrt().

Referenced by computeEfficiency().

                                                                                        {
   // looping over all bins. Up to three dimentions can be handled
   // in case of less dimensions, the inner for loops are excecuted only once
   for (int iBinX = 1; iBinX < denom->GetNbinsX() + 1; iBinX++) {
     for (int iBinY = 1; iBinY < denom->GetNbinsY() + 1; iBinY++) {
       for (int iBinZ = 1; iBinZ < denom->GetNbinsZ() + 1; iBinZ++) {
         int globalBinNum = denom->GetBin(iBinX, iBinY, iBinZ);
 
         // getting numerator and denominator values
         float numerVal = numer->GetBinContent(globalBinNum);
         float denomVal = denom->GetBinContent(globalBinNum);
 
         // calculating effiency
         float effVal = 0;
         effVal = denomVal ? numerVal / denomVal : 0;
 
         // calculating error
         float errVal = 0;
         errVal = (denomVal && (effVal <= 1)) ? sqrt(effVal * (1 - effVal) / denomVal) : 0;
 
         // inserting value into the efficiency histogram
         efficiencyHist->setBinContent(globalBinNum, effVal);
         efficiencyHist->setBinError(globalBinNum, errVal);
         efficiencyHist->setEfficiencyFlag();
       }
     }
   }
 }

◆ makeAllPlots()

void HLTGenValClient::makeAllPlots	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

private

Definition at line 140 of file HLTGenValClient.cc.

References computeEfficiency(), Skims_PA_cff::content, relmon_rootfiles_spy::contents, TrackerOfflineValidation_Dqm_cff::dirName, efficOptions_, findAllSubdirectories(), dqm::implementation::IGetter::getAllContents(), isWildcardUsed_, metacharacters_, Skims_PA_cff::name, runTheMatrix::opt, topSingleLeptonDQM_PU_cfi::pattern, separator_, AlCaHLTBitMon_QueryRunRegistry::string, PostProcessorHcalIsoTrack_cfi::subDir, subDirs_, makeGlobalPositionRcd_cfg::tag, and unpackBuffers-CaloStage2::token.

Referenced by dqmEndLuminosityBlock(), and dqmEndRun().

                                                                                      {
   // Process wildcard in the sub-directory
   std::set<std::string> subDirSet;
   for (auto& subDir : subDirs_) {
     if (subDir[subDir.size() - 1] == '/')
       subDir.erase(subDir.size() - 1);
 
     if (TString(subDir).Contains(metacharacters_)) {
       isWildcardUsed_ = true;
 
       const std::string::size_type shiftPos = subDir.rfind('/');
       const std::string searchPath = subDir.substr(0, shiftPos);
       const std::string pattern = subDir.substr(shiftPos + 1, subDir.length());
 
       findAllSubdirectories(ibooker, igetter, searchPath, &subDirSet, pattern);
 
     } else {
       subDirSet.insert(subDir);
     }
   }
 
   // loop through all sub-directories
   // from the current implementation of the HLTGenValSource, we expect all histograms in a single directory
   // however, this module is also capable of handling sub-directories, if needed
   for (std::set<std::string>::const_iterator iSubDir = subDirSet.begin(); iSubDir != subDirSet.end(); ++iSubDir) {
     const std::string& dirName = *iSubDir;
 
     // construct efficiency options automatically from systematically names histograms^
     const auto contents = igetter.getAllContents(dirName);
     for (const auto& content : contents) {
       // splitting the input string
       std::string name = content->getName();
       std::vector<std::string> seglist;
       size_t pos = 0;
       std::string token;
       while ((pos = name.find(separator_)) != std::string::npos) {
         token = name.substr(0, pos);
         seglist.push_back(token);
         name.erase(0, pos + separator_.length());
       }
       seglist.push_back(name);
 
       if (seglist.size() == 4 ||
           seglist.size() ==
               5) {  // this should be the only "proper" files we want to look at. 5 means that a custom tag was set!
         if (seglist.at(2) == "GEN")
           continue;  // this is the "before" hist, we won't create an effiency from this alone
 
         // if a fifth entry exists, it is expected to be the custom tag
         std::string tag = "";
         if (seglist.size() == 5)
           tag = seglist.at(4);
 
         // first we determing whether we have the 1D or 2D case
         if (seglist.at(3).rfind("2D", 0) == 0) {
           // 2D case
           EfficOption opt;
           opt.name = seglist.at(0) + separator_ + seglist.at(1) + separator_ + seglist.at(2) + separator_ +
                      seglist.at(3) + separator_ + "eff";  // efficiency histogram name
           opt.title = seglist.at(0) + " " + seglist.at(1) + " " + seglist.at(2) + " " + seglist.at(3) +
                       " efficiency";           // efficiency histogram title
           opt.numerator = content->getName();  // numerator histogram (after a filter)
           opt.denominator = seglist.at(0) + separator_ + seglist.at(1) + separator_ + "GEN" + separator_ +
                             seglist.at(3);  // denominator histogram (before all filters)
 
           efficOptions_.push_back(opt);
 
         } else {
           // 1D case
           EfficOption opt;
           opt.name = seglist.at(0) + separator_ + seglist.at(1) + separator_ + seglist.at(2) + separator_ +
                      seglist.at(3) + separator_ + "eff";  // efficiency histogram name
           opt.title = seglist.at(0) + " " + seglist.at(1) + " " + seglist.at(2) + " " + seglist.at(3) +
                       " efficiency";           // efficiency histogram title
           opt.numerator = content->getName();  // numerator histogram (after a filter)
           opt.denominator = seglist.at(0) + separator_ + seglist.at(1) + separator_ + "GEN" + separator_ +
                             seglist.at(3);  // denominator histogram (before all filters)
 
           // propagating the custom tag to the efficiency
           if (!tag.empty()) {
             opt.name += separator_ + tag;
             opt.title += " " + tag;
             opt.denominator += separator_ + tag;
           }
 
           efficOptions_.push_back(opt);
         }
       }
     }
 
     // now that we have all EfficOptions, we create the histograms
     for (const auto& efficOption : efficOptions_) {
       computeEfficiency(ibooker,
                         igetter,
                         dirName,
                         efficOption.name,
                         efficOption.title,
                         efficOption.numerator,
                         efficOption.denominator);
     }
   }
 }

Member Data Documentation

◆ efficOptions_

std::vector<EfficOption> HLTGenValClient::efficOptions_

private

Definition at line 74 of file HLTGenValClient.cc.

Referenced by makeAllPlots().

◆ isWildcardUsed_

bool HLTGenValClient::isWildcardUsed_

private

Definition at line 68 of file HLTGenValClient.cc.

Referenced by computeEfficiency(), HLTGenValClient(), and makeAllPlots().

◆ makeGlobalEffPlot_

bool HLTGenValClient::makeGlobalEffPlot_

private

Definition at line 67 of file HLTGenValClient.cc.

Referenced by computeEfficiency(), and HLTGenValClient().

◆ metacharacters_

TPRegexp HLTGenValClient::metacharacters_

private

Definition at line 62 of file HLTGenValClient.cc.

Referenced by makeAllPlots().

◆ nonPerlWildcard_

TPRegexp HLTGenValClient::nonPerlWildcard_

private

Definition at line 63 of file HLTGenValClient.cc.

Referenced by findAllSubdirectories().

◆ outputFileName_

std::string HLTGenValClient::outputFileName_

private

Definition at line 72 of file HLTGenValClient.cc.

Referenced by dqmEndRun(), and HLTGenValClient().

◆ runOnEndJob_

bool HLTGenValClient::runOnEndJob_

private

Definition at line 66 of file HLTGenValClient.cc.

Referenced by dqmEndRun(), and HLTGenValClient().

◆ runOnEndLumi_

bool HLTGenValClient::runOnEndLumi_

private

Definition at line 65 of file HLTGenValClient.cc.

Referenced by dqmEndLuminosityBlock(), and HLTGenValClient().

◆ separator_

const std::string HLTGenValClient::separator_ = "__"

private

Definition at line 76 of file HLTGenValClient.cc.

Referenced by makeAllPlots().

◆ subDirs_

std::vector<std::string> HLTGenValClient::subDirs_

private

Definition at line 71 of file HLTGenValClient.cc.

Referenced by HLTGenValClient(), and makeAllPlots().

◆ theDQM

DQMStore* HLTGenValClient::theDQM

private

Definition at line 70 of file HLTGenValClient.cc.

Referenced by dqmEndRun().

◆ verbose_

unsigned int HLTGenValClient::verbose_

private

Definition at line 64 of file HLTGenValClient.cc.

Referenced by computeEfficiency(), and HLTGenValClient().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HLTGenValClient()

◆ ~HLTGenValClient()

Member Function Documentation

◆ computeEfficiency()

◆ dqmEndJob()

◆ dqmEndLuminosityBlock()

◆ dqmEndRun()

◆ findAllSubdirectories()

◆ genericEff()

◆ makeAllPlots()

Member Data Documentation

◆ efficOptions_

◆ isWildcardUsed_

◆ makeGlobalEffPlot_

◆ metacharacters_

◆ nonPerlWildcard_

◆ outputFileName_

◆ runOnEndJob_

◆ runOnEndLumi_

◆ separator_

◆ subDirs_

◆ theDQM

◆ verbose_