#include <L1TEfficiencyHarvesting.h>

Public Member Functions
void	book (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter)

void	computeEfficiency (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter)

	L1TEfficiencyPlotHandler (const edm::ParameterSet &ps, std::string plotName)

	L1TEfficiencyPlotHandler (const L1TEfficiencyPlotHandler &handler)

	~L1TEfficiencyPlotHandler ()

Private Attributes
std::string	denominatorDir_

std::string	denominatorSuffix_

MonitorElement *	h_efficiency_

std::string	numeratorDir_

std::string	numeratorSuffix_

std::string	outputDir_

std::string	plotName_

Detailed Description

Definition at line 36 of file L1TEfficiencyHarvesting.h.

Constructor & Destructor Documentation

dqmoffline::l1t::L1TEfficiencyPlotHandler::L1TEfficiencyPlotHandler	(	const edm::ParameterSet &	ps,
		std::string	plotName
	)

Definition at line 17 of file L1TEfficiencyHarvesting.cc.

                                                                                               :
         numeratorDir_(ps.getUntrackedParameter < std::string > ("numeratorDir")),
         denominatorDir_(ps.getUntrackedParameter < std::string > ("denominatorDir", numeratorDir_)),
         outputDir_(ps.getUntrackedParameter < std::string > ("outputDir", numeratorDir_)),
         plotName_(plotName),
         numeratorSuffix_(ps.getUntrackedParameter < std::string > ("numeratorSuffix", "Num")),
         denominatorSuffix_(ps.getUntrackedParameter < std::string > ("denominatorSuffix", "Den")),
         h_efficiency_()
 {
 
 }

dqmoffline::l1t::L1TEfficiencyPlotHandler::L1TEfficiencyPlotHandler ( const L1TEfficiencyPlotHandler & handler )

Definition at line 29 of file L1TEfficiencyHarvesting.cc.

                                                                                           :
         numeratorDir_(handler.numeratorDir_),
         denominatorDir_(handler.denominatorDir_),
         outputDir_(handler.outputDir_),
         plotName_(handler.plotName_),
         numeratorSuffix_(handler.numeratorSuffix_),
         denominatorSuffix_(handler.denominatorSuffix_),
         h_efficiency_(handler.h_efficiency_)
 {
 
 }

dqmoffline::l1t::L1TEfficiencyPlotHandler::~L1TEfficiencyPlotHandler ( )

inline

Definition at line 44 of file L1TEfficiencyHarvesting.h.

References book(), and computeEfficiency().

45 {

46 }

Member Function Documentation

void dqmoffline::l1t::L1TEfficiencyPlotHandler::book	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

Definition at line 41 of file L1TEfficiencyHarvesting.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), denominatorDir_, denominatorSuffix_, MonitorElement::DQM_KIND_TH1D, MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, DQMStore::IGetter::get(), MonitorElement::getTH1(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), h_efficiency_, fastTrackerRecHitType::is2D(), callgraph::kind, MonitorElement::kind(), pileupDistInMC::num, numeratorDir_, numeratorSuffix_, outputDir_, plotName_, DQMStore::IBooker::setCurrentFolder(), MonitorElement::setEfficiencyFlag(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~L1TEfficiencyPlotHandler().

 {
 
   edm::LogInfo("L1TEfficiencyPlotHandler") << "Booking efficiency histogram for " << outputDir_ << " and " << plotName_
       << endl;
 
   std::string numeratorName = numeratorDir_ + "/" + plotName_ + numeratorSuffix_;
   std::string denominatorName = denominatorDir_ + "/" + plotName_ + denominatorSuffix_;
   MonitorElement *num = igetter.get(numeratorName);
   MonitorElement *den = igetter.get(denominatorName);
 
   if (!num || !den) {
 
     edm::LogWarning("L1TEfficiencyPlotHandler")
         << (!num && !den ? numeratorName + " && " + denominatorName : !num ? numeratorName : denominatorName)
         << " not gettable. Quitting booking" << endl;
     return;
   }
 
   TH1 *numH = num->getTH1();
   TH1 *denH = den->getTH1();
 
   if (!numH || !denH) {
 
     edm::LogWarning("L1TEfficiencyPlotHandler")
         << (!numH && !denH ? numeratorName + " && " + denominatorName : !num ? numeratorName : denominatorName)
         << " is not TH1F. Quitting booking" << endl;
 
     return;
 
   }
 
   if (numH->GetNbinsX() != denH->GetNbinsX()) {
     edm::LogWarning(
       "L1TEfficiencyPlotHandler") << " # X bins in " <<
       numeratorName << " and " << denominatorName <<
       " are different. Quitting booking" << endl;
     return;
   }
 
   MonitorElement::Kind kind = num->kind();
   bool is1D = kind == MonitorElement::DQM_KIND_TH1F ||
               kind == MonitorElement::DQM_KIND_TH1D;
   bool is2D = kind == MonitorElement::DQM_KIND_TH2F ||
               kind == MonitorElement::DQM_KIND_TH2D;
 
   if (is2D) {
     if (numH->GetNbinsY() != denH->GetNbinsY()) {
       edm::LogWarning(
         "L1TEfficiencyPlotHandler") << " # Y bins in " << numeratorName <<
         " and " << denominatorName << " are different. Quitting booking" <<
         endl;
       return;
     }
   }
 
   ibooker.setCurrentFolder(outputDir_);
   if (is1D) {
     h_efficiency_ = ibooker.book1D(plotName_, den->getTH1F());
   }
   else if (is2D) {
     h_efficiency_ = ibooker.book2D(plotName_, den->getTH2F());
   }
   h_efficiency_->setEfficiencyFlag();
 }

void dqmoffline::l1t::L1TEfficiencyPlotHandler::computeEfficiency	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

Definition at line 107 of file L1TEfficiencyHarvesting.cc.

References denominatorDir_, denominatorSuffix_, DQMStore::IGetter::get(), MonitorElement::getTH1(), h_efficiency_, pileupDistInMC::num, numeratorDir_, numeratorSuffix_, and plotName_.

Referenced by ~L1TEfficiencyPlotHandler().

 {
   if (!h_efficiency_)
     return;
 
   edm::LogInfo("L1TEfficiencyPlotHandler") << " Computing efficiency for " << plotName_ << endl;
 
   MonitorElement *num = igetter.get(numeratorDir_ + "/" + plotName_ + numeratorSuffix_);
   MonitorElement *den = igetter.get(denominatorDir_ + "/" + plotName_ + denominatorSuffix_);
 
   TH1 *numH = num->getTH1();
   TH1 *denH = den->getTH1();
   TH1 *effH = h_efficiency_->getTH1();
 
   effH->Divide(numH, denH, 1.0, 1.0, "B");
 }