#include <PFClient_JetRes.h>

Inheritance diagram for PFClient_JetRes:

Public Member Functions
	PFClient_JetRes (const edm::ParameterSet &parameterSet)

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 Member Functions
void	createEfficiencyPlots (DQMStore::IBooker &, DQMStore::IGetter &, std::string &folder, std::string &name)

void	createResolutionPlots (DQMStore::IBooker &, DQMStore::IGetter &, std::string &folder, std::string &name)

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

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

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

void	getHistogramParameters (MonitorElement *me_slice, double &avarage, double &rms, double &mean, double &sigma)

Private Attributes
std::vector< std::string >	effHistogramNames_

bool	efficiencyFlag_

std::vector< std::string >	folderNames_

std::vector< std::string >	histogramNames_

std::vector< int >	PtBins_

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

Constructor & Destructor Documentation

PFClient_JetRes::PFClient_JetRes ( const edm::ParameterSet & parameterSet )

Definition at line 18 of file PFClient_JetRes.cc.

References effHistogramNames_, efficiencyFlag_, folderNames_, edm::ParameterSet::getParameter(), histogramNames_, and PtBins_.

                                                                     {
   folderNames_ = parameterSet.getParameter<std::vector<std::string>>("FolderNames");
   histogramNames_ = parameterSet.getParameter<std::vector<std::string>>("HistogramNames");
   efficiencyFlag_ = parameterSet.getParameter<bool>("CreateEfficiencyPlots");
   effHistogramNames_ = parameterSet.getParameter<std::vector<std::string>>("HistogramNamesForEfficiencyPlots");
   PtBins_ = parameterSet.getParameter<std::vector<int>>("VariablePtBins");
 }

Member Function Documentation

void PFClient_JetRes::createEfficiencyPlots	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		std::string &	folder,
		std::string &	name
	)

private

Definition at line 205 of file PFClient_JetRes.cc.

Referenced by doEfficiency().

                                                              {
   MonitorElement *me1 = igetter.get(folder + "/" + name);
   MonitorElement *me2 = igetter.get(folder + "/" + name + "ref_");
   if (!me1 || !me2)
     return;
   MonitorElement *me_eff;
   if ((me1->kind() == MonitorElement::DQM_KIND_TH1F) && (me1->kind() == MonitorElement::DQM_KIND_TH1F)) {
     TH1 *th1 = me1->getTH1F();
     size_t nbinx = me1->getNbinsX();
 
     float xmin = th1->GetXaxis()->GetXmin();
     float xmax = th1->GetXaxis()->GetXmax();
     std::string xtit = me1->getAxisTitle(1);
     std::string tit_new;
     tit_new = ";" + xtit + ";Efficiency";
 
     ibooker.setCurrentFolder(folder);
     me_eff = ibooker.book1D("efficiency_" + name, tit_new, nbinx, xmin, xmax);
 
     double efficiency;
     me_eff->Reset();
     me_eff->setEfficiencyFlag();
     for (size_t ix = 1; ix < nbinx + 1; ++ix) {
       float val1 = me1->getBinContent(ix);
       float val2 = me2->getBinContent(ix);
       if (val2 > 0.0)
         efficiency = val1 / val2;
       else
         efficiency = 0;
       me_eff->setBinContent(ix, efficiency);
     }
   }
 }

void PFClient_JetRes::createResolutionPlots	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		std::string &	folder,
		std::string &	name
	)

private

Definition at line 70 of file PFClient_JetRes.cc.

References DQMStore::IBooker::book1D(), MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TH2S, DQMStore::IGetter::get(), getHistogramParameters(), MonitorElement::getNbinsY(), MonitorElement::getTH2F(), MonitorElement::kind(), SiStripPI::mean, PVValHelper::pT, PtBins_, SiStripPI::rms, MonitorElement::setBinContent(), MonitorElement::setBinError(), DQMStore::IBooker::setCurrentFolder(), MonitorElement::setEfficiencyFlag(), MonitorElement::setEntries(), AlCaHLTBitMon_QueryRunRegistry::string, fw3dlego::xbins, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::ymin.

Referenced by doSummaries().

                                                              {
   MonitorElement *me = igetter.get(folder + "/" + name);
   if (!me)
     return;
 
   MonitorElement *pT[PtBins_.size() - 1];
   std::vector<double> pTEntries(PtBins_.size() - 1, 0);
 
   // std::vector<std::string> pTRange (PtBins_.size() -1) ;
   // char* pTRange[PtBins_.size() -1] ;
   std::vector<TString> pTRange(PtBins_.size() - 1);
   // float pTCenter[PtBins_.size() -1] ;
 
   MonitorElement *me_average;
   MonitorElement *me_rms;
   MonitorElement *me_mean;
   MonitorElement *me_sigma;
 
   if ((me->kind() == MonitorElement::DQM_KIND_TH2F) || (me->kind() == MonitorElement::DQM_KIND_TH2S) ||
       (me->kind() == MonitorElement::DQM_KIND_TH2D)) {
     TH2 *th = me->getTH2F();
     // size_t nbinx = me->getNbinsX();
     size_t nbinx = PtBins_.size() - 1;
     size_t nbiny = me->getNbinsY();
 
     float ymin = th->GetYaxis()->GetXmin();
     float ymax = th->GetYaxis()->GetXmax();
     std::string xtit = th->GetXaxis()->GetTitle();
     // std::string ytit = th->GetYaxis()->GetTitle();
     std::string ytit = "#Deltap_{T}/p_{T}";
 
     float *xbins = new float[nbinx + 1];
     for (size_t ix = 1; ix < nbinx + 1; ++ix) {
       // xbins[ix-1] = th->GetBinLowEdge(ix);
       xbins[ix - 1] = PtBins_[ix - 1];
       // snprintf(pTRange[ix-1].data(), 15, "Pt%d_%d", PtBins_[ix-1],
       // PtBins_[ix]);
       pTRange[ix - 1] = TString::Format("Pt%d_%d", PtBins_[ix - 1], PtBins_[ix]);
       if (name == "BRdelta_et_Over_et_VS_et_")
         pTRange[ix - 1] = TString::Format("BRPt%d_%d", PtBins_[ix - 1], PtBins_[ix]);
       else if (name == "ERdelta_et_Over_et_VS_et_")
         pTRange[ix - 1] = TString::Format("ERPt%d_%d", PtBins_[ix - 1], PtBins_[ix]);
 
       // pTCenter[ix-1] = (PtBins_[ix] - PtBins_[ix-1]) / 2. ;
       if (ix == nbinx) {
         // xbins[ix] = th->GetXaxis()->GetBinUpEdge(ix);
         xbins[ix] = PtBins_[ix];
       }
     }
 
     std::string tit_new;
     ibooker.setCurrentFolder(folder);
     // MonitorElement* me_slice = ibooker.book1D("PFlowSlice", "PFlowSlice",
     // nbiny, ymin, ymax);
 
     tit_new = "Average " + ytit + ";" + xtit + ";Average_" + ytit;
     me_average = ibooker.book1D("average_" + name, tit_new, nbinx, xbins);
     me_average->setEfficiencyFlag();
     tit_new = "RMS " + ytit + ";" + xtit + ";RMS_" + ytit;
     me_rms = ibooker.book1D("rms_" + name, tit_new, nbinx, xbins);
     me_rms->setEfficiencyFlag();
     tit_new = ";" + xtit + ";Mean_" + ytit;
     me_mean = ibooker.book1D("mean_" + name, tit_new, nbinx, xbins);
     me_mean->setEfficiencyFlag();
     tit_new = ";" + xtit + ";Sigma_" + ytit;
     me_sigma = ibooker.book1D("sigma_" + name, tit_new, nbinx, xbins);
     me_sigma->setEfficiencyFlag();
 
     double average, rms, mean, sigma;
     for (size_t ix = 1; ix < nbinx + 1; ++ix) {
       // me_slice->Reset();
       if (name == "delta_et_Over_et_VS_et_")
         pT[ix - 1] = ibooker.book1D(
             pTRange[ix - 1], TString::Format("Total %s;%s;Events", ytit.data(), ytit.data()), nbiny, ymin, ymax);
       if (name == "BRdelta_et_Over_et_VS_et_")
         pT[ix - 1] = ibooker.book1D(
             pTRange[ix - 1], TString::Format("Barrel %s;%s;Events", ytit.data(), ytit.data()), nbiny, ymin, ymax);
       else if (name == "ERdelta_et_Over_et_VS_et_")
         pT[ix - 1] = ibooker.book1D(
             pTRange[ix - 1], TString::Format("Endcap %s;%s;Events", ytit.data(), ytit.data()), nbiny, ymin, ymax);
 
       for (size_t iy = 0; iy <= nbiny + 1; ++iy)  // add under and overflow
         if (th->GetBinContent(ix, iy)) {
           // me_slice->setBinContent(iy,th->GetBinContent(ix,iy));
           pT[ix - 1]->setBinContent(iy, th->GetBinContent(ix, iy));
           pT[ix - 1]->setBinError(iy, th->GetBinError(ix, iy));
           pTEntries[ix - 1] += th->GetBinContent(ix, iy);
         }
 
       pT[ix - 1]->setEntries(pTEntries[ix - 1]);
 
       // getHistogramParameters(me_slice, average, rms, mean, sigma);
       getHistogramParameters(pT[ix - 1], average, rms, mean, sigma);
       me_average->setBinContent(ix, average);
       me_rms->setBinContent(ix, rms);
       me_mean->setBinContent(ix, mean);
       me_sigma->setBinContent(ix, sigma);
     }
     // if (me_slice) igetter.removeElement(me_slice->getName());
     delete[] xbins;
   }
 }

void PFClient_JetRes::doEfficiency	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

private

Definition at line 54 of file PFClient_JetRes.cc.

References createEfficiencyPlots(), effHistogramNames_, folderNames_, callgraph::path, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dqmEndJob().

                                                                                      {
   for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin(); ifolder != folderNames_.end();
        ifolder++) {
     std::string path = "ParticleFlow/" + (*ifolder);
 
     for (std::vector<std::string>::const_iterator ihist = effHistogramNames_.begin(); ihist != effHistogramNames_.end();
          ihist++) {
       std::string hname = (*ihist);
       createEfficiencyPlots(ibooker, igetter, path, hname);
     }
   }
 }

void PFClient_JetRes::doSummaries	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

private

Definition at line 38 of file PFClient_JetRes.cc.

References createResolutionPlots(), folderNames_, histogramNames_, callgraph::path, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dqmEndJob().

                                                                                     {
   for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin(); ifolder != folderNames_.end();
        ifolder++) {
     std::string path = "ParticleFlow/" + (*ifolder);
 
     for (std::vector<std::string>::const_iterator ihist = histogramNames_.begin(); ihist != histogramNames_.end();
          ihist++) {
       std::string hname = (*ihist);
       createResolutionPlots(ibooker, igetter, path, hname);
     }
   }
 }

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

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 29 of file PFClient_JetRes.cc.

References doEfficiency(), doSummaries(), and efficiencyFlag_.

                                                                                   {
   doSummaries(ibooker, igetter);
   if (efficiencyFlag_)
     doEfficiency(ibooker, igetter);
 }

void PFClient_JetRes::getHistogramParameters	(	MonitorElement *	me_slice,
		double &	avarage,
		double &	rms,
		double &	mean,
		double &	sigma
	)

private

Definition at line 179 of file PFClient_JetRes.cc.

References MonitorElement::DQM_KIND_TH1F, MonitorElement::getMean(), MonitorElement::getRMS(), MonitorElement::getTH1F(), and MonitorElement::kind().

Referenced by createResolutionPlots().

                                                                                          {
   average = 0.0;
   rms = 0.0;
   mean = 0.0;
   sigma = 0.0;
 
   if (!me_slice)
     return;
   if (me_slice->kind() == MonitorElement::DQM_KIND_TH1F) {
     average = me_slice->getMean();
     rms = me_slice->getRMS();
     TH1F *th_slice = me_slice->getTH1F();
     if (th_slice && th_slice->GetEntries() > 0) {
       // need our own copy for thread safety
       TF1 gaus("mygaus", "gaus");
       th_slice->Fit(&gaus, "Q0 SERIAL");
       sigma = gaus.GetParameter(2);
       mean = gaus.GetParameter(1);
     }
   }
 }