#include <EmDQMPostProcessor.h>

Inheritance diagram for EmDQMPostProcessor:

Public Member Functions
TProfile *	dividehistos (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string &num, const std::string &denom, const std::string &out, const std::string &label, const std::string &titel="")

TH2F *	dividehistos2D (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string &num, const std::string &denom, const std::string &out, const std::string &label, const std::string &titel="")

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

	EmDQMPostProcessor (const edm::ParameterSet &pset)

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

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

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

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

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 Member Functions
TH2F *	get2DHistogram (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string &histoPath)

TH1F *	getHistogram (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string &histoPath)

Static Private Member Functions
static void	Efficiency (int passing, int total, double level, double &mode, double &lowerBound, double &upperBound)

Private Attributes
std::string	dataSet_

bool	ignoreEmpty

bool	noPhiPlots

bool	normalizeToReco

std::string	subDir_

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

Constructor & Destructor Documentation

◆ EmDQMPostProcessor()

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

Definition at line 17 of file EmDQMPostProcessor.cc.

References dataSet_, ignoreEmpty, noPhiPlots, normalizeToReco, muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, and subDir_.

                                                                   {
   subDir_ = pset.getUntrackedParameter<std::string>("subDir");
 
   dataSet_ = pset.getUntrackedParameter<std::string>("dataSet", "unknown");
 
   noPhiPlots = pset.getUntrackedParameter<bool>("noPhiPlots", true);
 
   normalizeToReco = pset.getUntrackedParameter<bool>("normalizeToReco", false);
 
   ignoreEmpty = pset.getUntrackedParameter<bool>("ignoreEmpty", true);
 }

◆ ~EmDQMPostProcessor()

EmDQMPostProcessor::~EmDQMPostProcessor ( )

inlineoverride

Definition at line 10 of file EmDQMPostProcessor.h.

10 {};

Member Function Documentation

◆ dividehistos()

TProfile * EmDQMPostProcessor::dividehistos	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const std::string &	num,
		const std::string &	denom,
		const std::string &	out,
		const std::string &	label,
		const std::string &	titel = `""`
	)

Definition at line 396 of file EmDQMPostProcessor.cc.

References dqm::implementation::IBooker::bookProfile(), makePileupJSON::denom, MillePedeFileConverter_cfg::e, Efficiency(), submitPVResolutionJobs::err, getHistogram(), dqm::legacy::MonitorElement::getTProfile(), LaserClient_cfi::high, mps_fire::i, createfilelist::int, label, LaserClient_cfi::low, EgammaValidation_cff::num, MillePedeFileConverter_cfg::out, dataset::outName, dqm::impl::MonitorElement::setEfficiencyFlag(), and mathSSE::sqrt().

Referenced by dqmEndJob().

                                                                    {
   // std::cout << numName <<std::endl;
 
   TH1F *num = getHistogram(ibooker, igetter, numName);
 
   // std::cout << denomName << std::endl;
   TH1F *denom = getHistogram(ibooker, igetter, denomName);
 
   if (num == nullptr)
     edm::LogWarning("EmDQMPostProcessor") << "numerator histogram " << numName << " does not exist";
 
   if (denom == nullptr)
     edm::LogWarning("EmDQMPostProcessor") << "denominator histogram " << denomName << " does not exist";
 
   // Check if histograms actually exist
 
   if (num == nullptr || denom == nullptr)
     return nullptr;
 
   MonitorElement *meOut = ibooker.bookProfile(
       outName, titel, num->GetXaxis()->GetNbins(), num->GetXaxis()->GetXmin(), num->GetXaxis()->GetXmax(), 0., 1.2);
   meOut->setEfficiencyFlag();
   TProfile *out = meOut->getTProfile();
   out->GetXaxis()->SetTitle(label.c_str());
   out->SetYTitle("Efficiency");
   out->SetOption("PE");
   out->SetLineColor(2);
   out->SetLineWidth(2);
   out->SetMarkerStyle(20);
   out->SetMarkerSize(0.8);
   out->SetStats(kFALSE);
   for (int i = 1; i <= num->GetNbinsX(); i++) {
     double e, low, high;
     Efficiency((int)num->GetBinContent(i), (int)denom->GetBinContent(i), 0.683, e, low, high);
     double err = e - low > high - e ? e - low : high - e;
     // here is the trick to store info in TProfile:
     out->SetBinContent(i, e);
     out->SetBinEntries(i, 1);
     out->SetBinError(i, sqrt(e * e + err * err));
   }
 
   return out;
 }

◆ dividehistos2D()

TH2F * EmDQMPostProcessor::dividehistos2D	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const std::string &	num,
		const std::string &	denom,
		const std::string &	out,
		const std::string &	label,
		const std::string &	titel = `""`
	)

Definition at line 446 of file EmDQMPostProcessor.cc.

References dqm::implementation::IBooker::book2D(), makePileupJSON::denom, get2DHistogram(), dqm::legacy::MonitorElement::getTH2F(), label, EgammaValidation_cff::num, MillePedeFileConverter_cfg::out, and dataset::outName.

Referenced by dqmEndJob().

                                                                  {
   // std::cout << numName <<std::endl;
   TH2F *num = get2DHistogram(ibooker, igetter, numName);
   // std::cout << denomName << std::endl;
   TH2F *denom = get2DHistogram(ibooker, igetter, denomName);
   if (num == nullptr)
     edm::LogWarning("EmDQMPostProcessor") << "2D numerator histogram " << numName << " does not exist";
 
   if (denom == nullptr)
     edm::LogWarning("EmDQMPostProcessor") << "2D denominator histogram " << denomName << " does not exist";
 
   // Check if histograms actually exist
   if (num == nullptr || denom == nullptr)
     return nullptr;
 
   MonitorElement *meOut = ibooker.book2D(outName,
                                          titel,
                                          num->GetXaxis()->GetNbins(),
                                          num->GetXaxis()->GetXmin(),
                                          num->GetXaxis()->GetXmax(),
                                          num->GetYaxis()->GetNbins(),
                                          num->GetYaxis()->GetXmin(),
                                          num->GetYaxis()->GetXmax());
   TH2F *out = meOut->getTH2F();
   out->Add(num);
   out->Divide(denom);
   out->GetXaxis()->SetTitle(label.c_str());
   out->SetYTitle("#phi");
   out->SetXTitle("#eta");
   out->SetOption("COLZ");
   out->SetStats(kFALSE);
 
   return out;
 }

◆ dqmEndJob()

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

overridevirtual

Implements DQMEDHarvester.

Definition at line 31 of file EmDQMPostProcessor.cc.

References newFWLiteAna::bin, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), dqm::implementation::NavigatorBase::cd(), dataSet_, DeadROC_duringRun::dir, dqm::implementation::IGetter::dirExists(), dividehistos(), dividehistos2D(), Efficiency(), relativeConstraints::error, ALCARECOTkAlBeamHalo_cff::filter, pdwgDoubleElectron_cfi::filterName, getHistogram(), dqm::implementation::IGetter::getSubdirs(), dqm::legacy::MonitorElement::getTProfile(), dqm::implementation::NavigatorBase::goUp(), HltBtagPostValidation_cff::histoName, ignoreEmpty, createfilelist::int, noPhiPlots, normalizeToReco, l1ctLayer1_cff::nPhoton, nano_postproc::postfix, dqm::implementation::NavigatorBase::pwd(), dqm::impl::MonitorElement::setEfficiencyFlag(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, subDir_, personalPlayback::subdirectories, groupFilesInBlocks::temp, dqmMemoryStats::total, cscTnPEfficiencyTask_cfi::trigName, relativeConstraints::value, trigObjTnPSource_cfi::var, and varNames.

                                                                                      {
   // go to the directory to be processed
   if (igetter.dirExists(subDir_))
     ibooker.cd(subDir_);
   else {
     edm::LogWarning("EmDQMPostProcessor") << "cannot find directory: " << subDir_ << " , skipping";
     return;
   }
 
   //--------------------
   // with respect to what are (some) efficiencies calculated ?
   std::string shortReferenceName;
   if (normalizeToReco)
     shortReferenceName = "RECO";
   else
     shortReferenceName = "gen";
   //--------------------
 
   // loop over all triggers/samples//
 
   // store dataset name (if defined) in output file
   // DQMStore:bookString seems to put a key in the file which is
   // of the form <dataSet>s=....</dataSet> which is not very convenient
   // (it points to a null pointer, one would have to loop over
   // all keys of the corresponding directory in the ROOT file
   // and check whether it is of the desired form and then parse
   // it from this string...).
   //
   // So we store the name of the dataset as the title of a histogram,
   // which is much easier to access...
   // TH1D *dataSetNameHisto =
   ibooker.book1D("DataSetNameHistogram", dataSet_, 1, 0, 1);
 
   std::vector<std::string> subdirectories = igetter.getSubdirs();
   // Do everything twice: once for mc-matched histos,   //
   // once for unmatched histos                          //
 
   std::vector<std::string> postfixes;
   postfixes.push_back("");               // unmatched histograms
   postfixes.push_back("_RECO_matched");  // for data
   // we put this on the list even when we're running on
   // data (where there is no generator information).
   // The first test in the loop will then fail and
   // the iteration is skipped.
   postfixes.push_back("_MC_matched");
 
   std::vector<TProfile *> allElePaths;
   int nEle = 0;
   int nPhoton = 0;
 
   // find the number of electron and photon paths
   for (std::vector<std::string>::iterator dir = subdirectories.begin(); dir != subdirectories.end(); ++dir) {
     if (dir->find("Ele") != std::string::npos || dir->find("_SC") != std::string::npos)
       ++nEle;
     else if (dir->find("Photon") != std::string::npos)
       ++nPhoton;
   }
 
   std::vector<TProfile *> allPhotonPaths;
   for (std::vector<std::string>::iterator postfix = postfixes.begin(); postfix != postfixes.end(); postfix++) {
     bool pop = false;
     int elePos = 1;
     int photonPos = 1;
 
     // computer per-event efficiencies //
 
     std::string histoName = "efficiency_by_step" + *postfix;
     std::string baseName = "total_eff" + *postfix;
 
     std::string allEleHistoName = "EfficiencyByPath_Ele" + *postfix;
     std::string allEleHistoLabel = "Efficiency_for_each_validated_electron_path" + *postfix;
     allElePaths.push_back(
         new TProfile(allEleHistoName.c_str(), allEleHistoLabel.c_str(), nEle, 0., (double)nEle, 0., 1.2));
     std::string allPhotonHistoName = "EfficiencyByPath_Photon" + *postfix;
     std::string allPhotonHistoLabel = "Efficiency_for_each_validated_photon_path" + *postfix;
     allPhotonPaths.push_back(
         new TProfile(allPhotonHistoName.c_str(), allPhotonHistoLabel.c_str(), nPhoton, 0., (double)nPhoton, 0., 1.2));
 
     for (std::vector<std::string>::iterator dir = subdirectories.begin(); dir != subdirectories.end(); dir++) {
       ibooker.cd(*dir);
 
       // get the current trigger name
       std::string trigName = dir->substr(dir->rfind("/") + 1);
       trigName = trigName.replace(trigName.rfind("_DQM"), 4, "");
 
       // Get the gen-level (or reco, for data) plots
       std::string genName;
 
       // only generate efficiency plots if there are generated/recoed events
       // selected but label the bin in the overview plot, even if the bin is
       // empty
       if (ignoreEmpty) {
         if (normalizeToReco)
           genName = ibooker.pwd() + "/reco_et";
         else
           genName = ibooker.pwd() + "/gen_et";
         TH1F *genHist = getHistogram(ibooker, igetter, genName);
         if (genHist->GetEntries() == 0) {
           edm::LogInfo("EmDQMPostProcessor")
               << "Zero events were selected for path '" << trigName << "'. No efficiency plots will be generated.";
           if (trigName.find("Ele") != std::string::npos || trigName.find("_SC") != std::string::npos) {
             allElePaths.back()->GetXaxis()->SetBinLabel(elePos, trigName.c_str());
             ++elePos;
           } else if (trigName.find("Photon") != std::string::npos) {
             allPhotonPaths.back()->GetXaxis()->SetBinLabel(photonPos, trigName.c_str());
             ++photonPos;
           }
           ibooker.goUp();
           continue;
         }
       }
 
       TH1F *basehist = getHistogram(ibooker, igetter, ibooker.pwd() + "/" + baseName);
       if (basehist == nullptr) {
         // edm::LogWarning("EmDQMPostProcessor") << "histogram " <<
         // (ibooker.pwd() + "/" + baseName) << " does not exist, skipping
         // postfix
         // '" << *postfix << "'";
         pop = true;
         ibooker.goUp();
         continue;
       }
       // at least one histogram with postfix was found
       pop = false;
 
       ibooker.goUp();
       MonitorElement *meTotal = ibooker.bookProfile(trigName + "__" + histoName,
                                                     trigName + "__" + histoName,
                                                     basehist->GetXaxis()->GetNbins(),
                                                     basehist->GetXaxis()->GetXmin(),
                                                     basehist->GetXaxis()->GetXmax(),
                                                     0.,
                                                     1.2);
       meTotal->setEfficiencyFlag();
       TProfile *total = meTotal->getTProfile();
       ibooker.cd(*dir);
       total->GetXaxis()->SetBinLabel(1, basehist->GetXaxis()->GetBinLabel(1));
 
       //       std::vector<std::string> mes = igetter.getMEs();
       //       for(std::vector<std::string>::iterator me = mes.begin() ;me!=
       //       mes.end(); me++ )
       //    std::cout <<*me <<std::endl;
       //       std::cout <<std::endl;
 
       double value = 0;
       double errorh = 0, errorl = 0, error = 0;
       // compute stepwise total efficiencies
       for (int bin = total->GetNbinsX() - 2; bin > 1; bin--) {
         value = 0;
         errorl = 0;
         errorh = 0;
         error = 0;
         if (basehist->GetBinContent(bin - 1) != 0) {
           Efficiency(
               (int)basehist->GetBinContent(bin), (int)basehist->GetBinContent(bin - 1), 0.683, value, errorl, errorh);
           error = value - errorl > errorh - value ? value - errorl : errorh - value;
         }
         total->SetBinContent(bin, value);
         total->SetBinEntries(bin, 1);
         total->SetBinError(bin, sqrt(value * value + error * error));
         total->GetXaxis()->SetBinLabel(bin, basehist->GetXaxis()->GetBinLabel(bin));
       }
 
       // set first bin to L1 efficiency
       if (basehist->GetBinContent(basehist->GetNbinsX()) != 0) {
         Efficiency((int)basehist->GetBinContent(1),
                    (int)basehist->GetBinContent(basehist->GetNbinsX()),
                    0.683,
                    value,
                    errorl,
                    errorh);
         error = value - errorl > errorh - value ? value - errorl : errorh - value;
       } else {
         value = 0;
         error = 0;
       }
       total->SetBinContent(1, value);
       total->SetBinEntries(1, 1);
       total->SetBinError(1, sqrt(value * value + error * error));
 
       // total efficiency relative to gen or reco
       if (basehist->GetBinContent(basehist->GetNbinsX()) != 0) {
         Efficiency((int)basehist->GetBinContent(basehist->GetNbinsX() - 2),
                    (int)basehist->GetBinContent(basehist->GetNbinsX()),
                    0.683,
                    value,
                    errorl,
                    errorh);
         error = value - errorl > errorh - value ? value - errorl : errorh - value;
       } else {
         value = 0;
         error = 0;
       }
       total->SetBinContent(total->GetNbinsX(), value);
       total->SetBinEntries(total->GetNbinsX(), 1);
       total->SetBinError(total->GetNbinsX(), sqrt(value * value + error * error));
       total->GetXaxis()->SetBinLabel(total->GetNbinsX(), ("total efficiency rel. " + shortReferenceName).c_str());
 
       // total efficiency relative to L1
       if (basehist->GetBinContent(1) != 0) {
         Efficiency((int)basehist->GetBinContent(basehist->GetNbinsX() - 2),
                    (int)basehist->GetBinContent(1),
                    0.683,
                    value,
                    errorl,
                    errorh);
         error = value - errorl > errorh - value ? value - errorl : errorh - value;
       } else {
         value = 0;
         error = 0;
       }
       total->SetBinContent(total->GetNbinsX() - 1, value);
       total->SetBinError(total->GetNbinsX() - 1, sqrt(value * value + error * error));
       total->SetBinEntries(total->GetNbinsX() - 1, 1);
       total->GetXaxis()->SetBinLabel(total->GetNbinsX() - 1, "total efficiency rel. L1");
 
       //----------------------------------------
 
       // compute per-object efficiencies       //
       // MonitorElement *eff, *num, *denom, *genPlot, *effVsGen, *effL1VsGen;
       std::vector<std::string> varNames;
       varNames.push_back("et");
       varNames.push_back("eta");
 
       if (!noPhiPlots) {
         varNames.push_back("phi");
       }
       varNames.push_back("etaphi");
 
       std::string filterName;
       std::string filterName2;
       std::string denomName;
       std::string numName;
 
       // Get the L1 over gen filter first
       filterName2 = total->GetXaxis()->GetBinLabel(1);
       // loop over variables (eta/phi/et)
       for (std::vector<std::string>::iterator var = varNames.begin(); var != varNames.end(); var++) {
         numName = ibooker.pwd() + "/" + filterName2 + *var + *postfix;
 
         if (normalizeToReco)
           genName = ibooker.pwd() + "/reco_" + *var;
         else
           genName = ibooker.pwd() + "/gen_" + *var;
 
         if ((*var).find("etaphi") != std::string::npos) {
           if (!dividehistos2D(ibooker,
                               igetter,
                               numName,
                               genName,
                               "efficiency_" + filterName2 + "_vs_" + *var + *postfix,
                               *var,
                               "eff. of" + filterName2 + " vs " + *var + *postfix))
             break;
         } else if (!dividehistos(ibooker,
                                  igetter,
                                  numName,
                                  genName,
                                  "efficiency_" + filterName2 + "_vs_" + *var + *postfix,
                                  *var,
                                  "eff. of" + filterName2 + " vs " + *var + *postfix))
           break;
       }  // loop over variables
 
       // get the filter names from the bin-labels of the master-histogram
       for (int filter = 1; filter < total->GetNbinsX() - 2; filter++) {
         filterName = total->GetXaxis()->GetBinLabel(filter);
         filterName2 = total->GetXaxis()->GetBinLabel(filter + 1);
 
         // loop over variables (eta/et/phi)
         for (std::vector<std::string>::iterator var = varNames.begin(); var != varNames.end(); var++) {
           numName = ibooker.pwd() + "/" + filterName2 + *var + *postfix;
           denomName = ibooker.pwd() + "/" + filterName + *var + *postfix;
 
           // Is this the last filter? Book efficiency vs gen (or reco, for data)
           // level
           std::string temp = *postfix;
           if (filter == total->GetNbinsX() - 3 && temp.find("matched") != std::string::npos) {
             if (normalizeToReco)
               genName = ibooker.pwd() + "/reco_" + *var;
             else
               genName = ibooker.pwd() + "/gen_" + *var;
 
             if ((*var).find("etaphi") != std::string::npos) {
               if (!dividehistos2D(ibooker,
                                   igetter,
                                   numName,
                                   genName,
                                   "final_eff_vs_" + *var,
                                   *var,
                                   "Efficiency Compared to " + shortReferenceName + " vs " + *var))
                 break;
             } else if (!dividehistos(ibooker,
                                      igetter,
                                      numName,
                                      genName,
                                      "final_eff_vs_" + *var,
                                      *var,
                                      "Efficiency Compared to " + shortReferenceName + " vs " + *var))
               break;
           }
 
           if ((*var).find("etaphi") != std::string::npos) {
             if (!dividehistos2D(ibooker,
                                 igetter,
                                 numName,
                                 denomName,
                                 "efficiency_" + filterName2 + "_vs_" + *var + *postfix,
                                 *var,
                                 "efficiency_" + filterName2 + "_vs_" + *var + *postfix))
               break;
           } else if (!dividehistos(ibooker,
                                    igetter,
                                    numName,
                                    denomName,
                                    "efficiency_" + filterName2 + "_vs_" + *var + *postfix,
                                    *var,
                                    "efficiency_" + filterName2 + "_vs_" + *var + *postfix))
             break;
 
         }  // loop over variables
       }    // loop over monitoring modules within path
 
       ibooker.goUp();
 
       // fill overall efficiency histograms
       double totCont = total->GetBinContent(total->GetNbinsX());
       double totErr = total->GetBinError(total->GetNbinsX());
       if (trigName.find("Ele") != std::string::npos || trigName.find("_SC") != std::string::npos) {
         allElePaths.back()->SetBinContent(elePos, totCont);
         allElePaths.back()->SetBinEntries(elePos, 1);
         allElePaths.back()->SetBinError(elePos, sqrt(totCont * totCont + totErr * totErr));
         allElePaths.back()->GetXaxis()->SetBinLabel(elePos, trigName.c_str());
         ++elePos;
       } else if (trigName.find("Photon") != std::string::npos) {
         allPhotonPaths.back()->SetBinContent(photonPos, totCont);
         allPhotonPaths.back()->SetBinEntries(photonPos, 1);
         allPhotonPaths.back()->SetBinError(photonPos, sqrt(totCont * totCont + totErr * totErr));
         allPhotonPaths.back()->GetXaxis()->SetBinLabel(photonPos, trigName.c_str());
         ++photonPos;
       }
 
     }  // loop over dirs
     if (pop) {
       allElePaths.pop_back();
       allPhotonPaths.pop_back();
     } else {
       allElePaths.back()->GetXaxis()->SetLabelSize(0.03);
       allPhotonPaths.back()->GetXaxis()->SetLabelSize(0.03);
       ibooker.bookProfile(allEleHistoName, allElePaths.back())->setEfficiencyFlag();
       ibooker.bookProfile(allPhotonHistoName, allPhotonPaths.back())->setEfficiencyFlag();
     }
   }  // loop over postfixes
 }

◆ Efficiency()

void EmDQMPostProcessor::Efficiency	(	int	passing,
		int	total,
		double	level,
		double &	mode,
		double &	lowerBound,
		double &	upperBound
	)

staticprivate

a replacement for the function TGraphAsymmErrors::Efficiency(..) used with earlier versions of ROOT (this functionality has been moved to a separate class TEfficiency)

Definition at line 509 of file EmDQMPostProcessor.cc.

References personalPlayback::level, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, TagProbeFitTreeAnalyzer_cfi::passing, and dqmMemoryStats::total.

Referenced by dividehistos(), and dqmEndJob().

                                                                                                 {
   // protection (see also TGraphAsymmErrors::Efficiency(..), mimick the old
   // behaviour )
   if (total == 0) {
     mode = 0.5;
     lowerBound = 0;
     upperBound = 1;
     return;
   }
 
   mode = passing / ((double)total);
 
   // note that the order of the two first arguments ('total' and 'passed') is
   // the opposited with respect to the earlier TGraphAsymmErrors::Efficiency(..)
   // method
   //
   // see http://root.cern.ch/root/html/TEfficiency.html#compare for
   // an illustration of the coverage of the different methods provided by
   // TEfficiency
 
   lowerBound = TEfficiency::Wilson(total, passing, level, false);
   upperBound = TEfficiency::Wilson(total, passing, level, true);
 }

◆ get2DHistogram()

TH2F * EmDQMPostProcessor::get2DHistogram	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const std::string &	histoPath
	)

private

Definition at line 497 of file EmDQMPostProcessor.cc.

References dqm::implementation::IGetter::get(), and dqm::legacy::MonitorElement::getTH2F().

Referenced by dividehistos2D().

                                                                      {
   MonitorElement *monElement = igetter.get(histoPath);
   if (monElement != nullptr)
     return monElement->getTH2F();
   else
     return nullptr;
 }

◆ getHistogram()

TH1F * EmDQMPostProcessor::getHistogram	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const std::string &	histoPath
	)

private

convenience method to get a histogram but checks first whether the corresponding MonitorElement is non-null.

Returns: null if the MonitorElement is null

Definition at line 487 of file EmDQMPostProcessor.cc.

References dqm::implementation::IGetter::get(), and dqm::legacy::MonitorElement::getTH1F().

Referenced by dividehistos(), and dqmEndJob().

                                                                    {
   MonitorElement *monElement = igetter.get(histoPath);
   if (monElement != nullptr)
     return monElement->getTH1F();
   else
     return nullptr;
 }

Member Data Documentation

◆ dataSet_

std::string EmDQMPostProcessor::dataSet_

private

dataset with which these histograms were produced. This is set by a user parameter in the configuration file.

It is just used for writing it to the DQM output file. Useful to remember with which dataset a histogram file was produced. This code does not do much with this information (apart from copying it to the output file) but it can be used when generating reports.

Definition at line 61 of file EmDQMPostProcessor.h.

Referenced by dqmEndJob(), and EmDQMPostProcessor().

◆ ignoreEmpty

bool EmDQMPostProcessor::ignoreEmpty

private

Definition at line 42 of file EmDQMPostProcessor.h.

Referenced by dqmEndJob(), and EmDQMPostProcessor().

◆ noPhiPlots

bool EmDQMPostProcessor::noPhiPlots

private

read from the configuration: if set to true, efficiencies are calculated with respect to reconstructed objects (instead of generated objects). This is e.g. a useful option when running on data.

Definition at line 40 of file EmDQMPostProcessor.h.

Referenced by dqmEndJob(), and EmDQMPostProcessor().

◆ normalizeToReco

bool EmDQMPostProcessor::normalizeToReco

private

Definition at line 41 of file EmDQMPostProcessor.h.

Referenced by dqmEndJob(), and EmDQMPostProcessor().

◆ subDir_

std::string EmDQMPostProcessor::subDir_

private

Definition at line 49 of file EmDQMPostProcessor.h.

Referenced by dqmEndJob(), and EmDQMPostProcessor().

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EmDQMPostProcessor()

◆ ~EmDQMPostProcessor()

Member Function Documentation

◆ dividehistos()

◆ dividehistos2D()

◆ dqmEndJob()

◆ Efficiency()

◆ get2DHistogram()

◆ getHistogram()

Member Data Documentation

◆ dataSet_

◆ ignoreEmpty

◆ noPhiPlots

◆ normalizeToReco

◆ subDir_