#include <HLTBTagHarvestingAnalyzer.h>

Inheritance diagram for HLTBTagHarvestingAnalyzer:

Public Member Functions
TH1F	calculateEfficiency1D (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, TH1 &num, TH1 &den, std::string name)

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

bool	GetNumDenumerators (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, std::string num, std::string den, TH1 &ptrnum, TH1 &ptrden, int type)

	HLTBTagHarvestingAnalyzer (const edm::ParameterSet &)

void	mistagrate (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, TH1F num, TH1F den, std::string effName)

void	modulesrate (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, TH1F num, TH1F den, std::string effName)

	~HLTBTagHarvestingAnalyzer () override

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 Types
typedef unsigned int	flavour_t

typedef std::vector< flavour_t >	flavours_t

enum	HCALSpecials { HEP17, HEP18, HEM17 }

Private Attributes
std::map< std::string, MonitorElement * >	H1_

std::map< HLTBTagHarvestingAnalyzer::HCALSpecials, std::string >	HCALSpecialsNames

std::vector< std::string >	hltPathNames_

std::vector< std::string >	m_histoName

std::vector< flavours_t >	m_mcFlavours

std::vector< std::string >	m_mcLabels

bool	m_mcMatching

double	m_minTag

std::string	mainFolder_

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

Code used to produce DQM validation plots for b-tag at HLT. This class read the plots producted by HLTBTagPerformanceAnalyzer and make plots of: b-tag efficiency vs discr, b-tag efficiency vs jet pt, b-tag efficiency vs mistag rate

Definition at line 27 of file HLTBTagHarvestingAnalyzer.h.

Member Typedef Documentation

typedef unsigned int HLTBTagHarvestingAnalyzer::flavour_t

private

Definition at line 49 of file HLTBTagHarvestingAnalyzer.h.

typedef std::vector<flavour_t> HLTBTagHarvestingAnalyzer::flavours_t

private

Definition at line 50 of file HLTBTagHarvestingAnalyzer.h.

Member Enumeration Documentation

enum HLTBTagHarvestingAnalyzer::HCALSpecials

private

Enumerator
HEP17
HEP18
HEM17

Definition at line 60 of file HLTBTagHarvestingAnalyzer.h.

60 { HEP17, HEP18, HEM17 };

HLTBTagHarvestingAnalyzer::HEM17

Definition: HLTBTagHarvestingAnalyzer.h:60

HLTBTagHarvestingAnalyzer::HEP18

Definition: HLTBTagHarvestingAnalyzer.h:60

HLTBTagHarvestingAnalyzer::HEP17

Definition: HLTBTagHarvestingAnalyzer.h:60

Constructor & Destructor Documentation

HLTBTagHarvestingAnalyzer::HLTBTagHarvestingAnalyzer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 3 of file HLTBTagHarvestingAnalyzer.cc.

References edm::ParameterSet::getParameter(), HCALSpecialsNames, HEM17, HEP17, HEP18, hltPathNames_, m_histoName, m_mcLabels, m_minTag, mainFolder_, CaloTowersParam_cfi::mc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                    {
   // getParameter
   mainFolder_ = iConfig.getParameter<std::string>("mainFolder");
   hltPathNames_ = iConfig.getParameter<std::vector<std::string>>("HLTPathNames");
   edm::ParameterSet mc = iConfig.getParameter<edm::ParameterSet>("mcFlavours");
   m_mcLabels = mc.getParameterNamesForType<std::vector<unsigned int>>();
   m_histoName = iConfig.getParameter<std::vector<std::string>>("histoName");
   m_minTag = iConfig.getParameter<double>("minTag");
 
   HCALSpecialsNames[HEP17] = "HEP17";
   HCALSpecialsNames[HEP18] = "HEP18";
   HCALSpecialsNames[HEM17] = "HEM17";
 }

HLTBTagHarvestingAnalyzer::~HLTBTagHarvestingAnalyzer ( )

override

Definition at line 17 of file HLTBTagHarvestingAnalyzer.cc.

17 {}

Member Function Documentation

TH1F HLTBTagHarvestingAnalyzer::calculateEfficiency1D	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		TH1 &	num,
		TH1 &	den,
		std::string	name
	)

Definition at line 340 of file HLTBTagHarvestingAnalyzer.cc.

References DQMStore::IBooker::book1D(), edmIntegrityCheck::d, DEFINE_FWK_MODULE, MillePedeFileConverter_cfg::e, mps_fire::i, SiStripPI::max, gen::n, and MonitorElement::setEfficiencyFlag().

Referenced by dqmEndJob().

                                                                                                {
   // calculate the efficiency as num/den ratio
   TH1F eff;
   if (num.GetXaxis()->GetXbins()->GetSize() == 0) {
     eff = TH1F(effName.c_str(),
                effName.c_str(),
                num.GetXaxis()->GetNbins(),
                num.GetXaxis()->GetXmin(),
                num.GetXaxis()->GetXmax());
   } else {
     eff = TH1F(effName.c_str(), effName.c_str(), num.GetXaxis()->GetNbins(), num.GetXaxis()->GetXbins()->GetArray());
   }
   eff.SetTitle(effName.c_str());
   eff.SetXTitle(num.GetXaxis()->GetTitle());
   eff.SetYTitle("Efficiency");
   eff.SetOption("PE");
   eff.SetLineColor(2);
   eff.SetLineWidth(2);
   eff.SetMarkerStyle(20);
   eff.SetMarkerSize(0.8);
   eff.GetYaxis()->SetRangeUser(-0.001, 1.001);
   for (int i = 1; i <= num.GetNbinsX(); i++) {
     double d, n, err;
     d = den.GetBinContent(i);
     n = num.GetBinContent(i);
     double e;
     if (d != 0) {
       e = n / d;
       err = std::max(e - TEfficiency::ClopperPearson(d, n, 0.683, false),
                      TEfficiency::ClopperPearson(d, n, 0.683, true) - e);
       // err =  sqrt(e*(1-e)/d); //from binomial standard deviation
     } else {
       e = 0;
       err = 0;
     }
     eff.SetBinContent(i, e);
     eff.SetBinError(i, err);
   }
 
   MonitorElement *me;
   me = ibooker.book1D(effName, &eff);
   me->setEfficiencyFlag();
 
   return eff;
 }

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

overridevirtual

save efficiency_vs_disc_HEP17 / efficiency_vs_disc_HEM17 plots

for mc labels

save mistagrate vs b-eff plots

save mistagrate vs b-eff plots for modules (HEP17 etc.)

save mistagrate_HEP17 / mistagrate_HEM17 plots

for triggers

Implements DQMEDHarvester.

Definition at line 21 of file HLTBTagHarvestingAnalyzer.cc.

References calculateEfficiency1D(), GenHFHadronMatcher_cff::flavour, GetNumDenumerators(), HCALSpecialsNames, HEM17, HEP17, hltPathNames_, mps_fire::i, label, edm::errors::LogicError, m_histoName, m_mcLabels, mainFolder_, mistagrate(), modulesrate(), pileupDistInMC::num, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                             {
   using namespace edm;
   Exception excp(errors::LogicError);
   std::string dqmFolder_hist;
 
   // for each hltPath and for each flavour, do the "b-tag efficiency vs jet pt"
   // and "b-tag efficiency vs mistag rate" plots
   for (unsigned int ind = 0; ind < hltPathNames_.size(); ind++) {
     dqmFolder_hist = Form("%s/Discriminator/%s", mainFolder_.c_str(), hltPathNames_[ind].c_str());
     std::string effDir = Form("%s/Discriminator/%s/efficiency", mainFolder_.c_str(), hltPathNames_[ind].c_str());
     std::string relationsDir = Form("%s/Discriminator/%s/HEP17_HEM17", mainFolder_.c_str(), hltPathNames_[ind].c_str());
     ibooker.setCurrentFolder(effDir);
     TH1 *den = nullptr;
     TH1 *num = nullptr;
     std::map<std::string, TH1F> effics;
     std::map<std::string, bool> efficsOK;
     std::map<std::string, std::map<HCALSpecials, TH1F>> efficsmod;
     std::map<std::string, std::map<HCALSpecials, bool>> efficsmodOK;
     for (unsigned int i = 0; i < m_mcLabels.size(); ++i) {
       bool isOK = false;
       std::string label = m_histoName.at(ind) + "__";  //"JetTag__";
       std::string flavour = m_mcLabels[i];
       label += flavour;
       isOK =
           GetNumDenumerators(ibooker, igetter, dqmFolder_hist + "/" + label, dqmFolder_hist + "/" + label, num, den, 0);
       if (isOK) {
         // do the 'b-tag efficiency vs discr' plot
         effics[flavour] = calculateEfficiency1D(ibooker, igetter, *num, *den, label + "_efficiency_vs_disc");
         efficsOK[flavour] = isOK;
       }
       // for modules (HEP17 etc.)
       for (auto j : HCALSpecialsNames) {
         ibooker.setCurrentFolder(dqmFolder_hist + "/" + j.second + "/efficiency");
         isOK = GetNumDenumerators(ibooker,
                                   igetter,
                                   dqmFolder_hist + "/" + j.second + "/" + label,
                                   dqmFolder_hist + "/" + j.second + "/" + label,
                                   num,
                                   den,
                                   0);
         if (isOK) {
           // do the 'b-tag efficiency vs discr' plot
           efficsmod[flavour][j.first] =
               calculateEfficiency1D(ibooker, igetter, *num, *den, label + "_efficiency_vs_disc");
           efficsmodOK[flavour][j.first] = isOK;
         }
       }
       ibooker.setCurrentFolder(effDir);
       label = m_histoName.at(ind) + "___";
       std::string labelEta = label;
       std::string labelPhi = label;
       label += flavour + "_disc_pT";
       labelEta += flavour + "_disc_eta";
       labelPhi += flavour + "_disc_phi";
       isOK =
           GetNumDenumerators(ibooker, igetter, dqmFolder_hist + "/" + label, dqmFolder_hist + "/" + label, num, den, 1);
       if (isOK) {
         // do the 'b-tag efficiency vs pT' plot
         TH1F eff = calculateEfficiency1D(ibooker, igetter, *num, *den, label + "_efficiency_vs_pT");
       }
       isOK = GetNumDenumerators(
           ibooker, igetter, dqmFolder_hist + "/" + labelEta, dqmFolder_hist + "/" + labelEta, num, den, 2);
       if (isOK) {
         // do the 'b-tag efficiency vs Eta' plot
         TH1F eff = calculateEfficiency1D(ibooker, igetter, *num, *den, labelEta + "_efficiency_vs_eta");
       }
       isOK = GetNumDenumerators(
           ibooker, igetter, dqmFolder_hist + "/" + labelPhi, dqmFolder_hist + "/" + labelPhi, num, den, 2);
       if (isOK) {
         // do the 'b-tag efficiency vs Phi' plot
         TH1F eff = calculateEfficiency1D(ibooker, igetter, *num, *den, labelPhi + "_efficiency_vs_phi");
       }
 
       ibooker.setCurrentFolder(relationsDir);
       if (efficsmodOK[flavour][HEP17] && efficsmodOK[flavour][HEM17])
         modulesrate(ibooker,
                     igetter,
                     &efficsmod[flavour][HEP17],
                     &efficsmod[flavour][HEM17],
                     m_histoName.at(ind) + "_" + flavour + "_HEP17_HEM17_effs_vs_disc_rate");
       ibooker.setCurrentFolder(effDir);
 
     }  
 
     if (efficsOK["b"] && efficsOK["c"])
       mistagrate(ibooker, igetter, &effics["b"], &effics["c"], m_histoName.at(ind) + "_b_c_mistagrate");
     if (efficsOK["b"] && efficsOK["light"])
       mistagrate(ibooker, igetter, &effics["b"], &effics["light"], m_histoName.at(ind) + "_b_light_mistagrate");
     if (efficsOK["b"] && efficsOK["g"])
       mistagrate(ibooker, igetter, &effics["b"], &effics["g"], m_histoName.at(ind) + "_b_g_mistagrate");
 
     for (auto j : HCALSpecialsNames) {
       ibooker.setCurrentFolder(dqmFolder_hist + "/" + j.second + "/efficiency");
       if (efficsmodOK["b"][j.first] && efficsmodOK["c"][j.first])
         mistagrate(ibooker,
                    igetter,
                    &efficsmod["b"][j.first],
                    &efficsmod["c"][j.first],
                    m_histoName.at(ind) + "_b_c_mistagrate");
       if (efficsmodOK["b"][j.first] && efficsmodOK["light"][j.first])
         mistagrate(ibooker,
                    igetter,
                    &efficsmod["b"][j.first],
                    &efficsmod["light"][j.first],
                    m_histoName.at(ind) + "_b_light_mistagrate");
       if (efficsmodOK["b"][j.first] && efficsmodOK["g"][j.first])
         mistagrate(ibooker,
                    igetter,
                    &efficsmod["b"][j.first],
                    &efficsmod["g"][j.first],
                    m_histoName.at(ind) + "_b_g_mistagrate");
     }
 
     ibooker.setCurrentFolder(relationsDir);
     bool isOK = false;
     isOK = GetNumDenumerators(ibooker,
                               igetter,
                               dqmFolder_hist + "/HEP17/efficiency/" + m_histoName.at(ind) + "_b_c_mistagrate",
                               dqmFolder_hist + "/HEM17/efficiency/" + m_histoName.at(ind) + "_b_c_mistagrate",
                               num,
                               den,
                               3);
     if (isOK)
       modulesrate(ibooker, igetter, (TH1F *)num, (TH1F *)den, m_histoName.at(ind) + "_HEP17_HEM17_b_c_mistagrate");
     isOK = GetNumDenumerators(ibooker,
                               igetter,
                               dqmFolder_hist + "/HEP17/efficiency/" + m_histoName.at(ind) + "_b_light_mistagrate",
                               dqmFolder_hist + "/HEM17/efficiency/" + m_histoName.at(ind) + "_b_light_mistagrate",
                               num,
                               den,
                               3);
     if (isOK)
       modulesrate(ibooker, igetter, (TH1F *)num, (TH1F *)den, m_histoName.at(ind) + "_HEP17_HEM17_b_light_mistagrate");
     isOK = GetNumDenumerators(ibooker,
                               igetter,
                               dqmFolder_hist + "/HEP17/efficiency/" + m_histoName.at(ind) + "_b_g_mistagrate",
                               dqmFolder_hist + "/HEM17/efficiency/" + m_histoName.at(ind) + "_b_g_mistagrate",
                               num,
                               den,
                               3);
     if (isOK)
       modulesrate(ibooker, igetter, (TH1F *)num, (TH1F *)den, m_histoName.at(ind) + "_HEP17_HEM17_b_g_mistagrate");
   }  
 }

bool HLTBTagHarvestingAnalyzer::GetNumDenumerators	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		std::string	num,
		std::string	den,
		TH1 *&	ptrnum,
		TH1 *&	ptrden,
		int	type
	)

numerator preparing

denominator preparing

numerator preparing

denominator preparing

Definition at line 170 of file HLTBTagHarvestingAnalyzer.cc.

References DQMStore::IGetter::get(), MonitorElement::getTH1(), MonitorElement::getTH2F(), edm::errors::LogicError, m_minTag, and cms::Exception::raise().

Referenced by dqmEndJob().

                                                              {
   using namespace edm;
   /*
      possible types:
      type =0 for eff_vs_discriminator
      type =1 for eff_vs_pT
      type =2 for eff_vs_eta or eff_vs_phi
      type =3 for HEP17 / HEM17 mistagrate relation
    */
   MonitorElement *denME = nullptr;
   MonitorElement *numME = nullptr;
   denME = igetter.get(den);
   numME = igetter.get(num);
   Exception excp(errors::LogicError);
 
   if (denME == nullptr || numME == nullptr) {
     excp << "Plots not found:\n";
     if (denME == nullptr)
       excp << den << "\n";
     if (numME == nullptr)
       excp << num << "\n";
     excp.raise();
   }
 
   if (type == 0)  // efficiency_vs_discr: fill "ptrnum" with the cumulative function of
                   // the DQM plots contained in "num" and "ptrden" with a flat function
   {
     TH1 *numH1 = numME->getTH1();
     TH1 *denH1 = denME->getTH1();
     ptrden = (TH1 *)denH1->Clone("denominator");
     ptrnum = (TH1 *)numH1->Clone("numerator");
 
     ptrnum->SetBinContent(1, numH1->Integral());
     ptrden->SetBinContent(1, numH1->Integral());
     for (int j = 2; j <= numH1->GetNbinsX(); j++) {
       ptrnum->SetBinContent(j, numH1->Integral() - numH1->Integral(1, j - 1));
       ptrden->SetBinContent(j, numH1->Integral());
     }
   }
 
   if (type == 1)  // efficiency_vs_pT: fill "ptrden" with projection of the plots
                   // contained in "den" and fill "ptrnum" with projection of the
                   // plots contained in "num", having btag>m_minTag
   {
     TH2F *numH2 = numME->getTH2F();
     TH2F *denH2 = denME->getTH2F();
 
     TCutG *cutg_num = new TCutG("cutg_num", 4);
     cutg_num->SetPoint(0, m_minTag, 0);
     cutg_num->SetPoint(1, m_minTag, 9999);
     cutg_num->SetPoint(2, 1.1, 9999);
     cutg_num->SetPoint(3, 1.1, 0);
     ptrnum = numH2->ProjectionY("numerator", 0, -1, "[cutg_num]");
 
     TCutG *cutg_den = new TCutG("cutg_den", 4);
     cutg_den->SetPoint(0, -10.1, 0);
     cutg_den->SetPoint(1, -10.1, 9999);
     cutg_den->SetPoint(2, 1.1, 9999);
     cutg_den->SetPoint(3, 1.1, 0);
     ptrden = denH2->ProjectionY("denumerator", 0, -1, "[cutg_den]");
     delete cutg_num;
     delete cutg_den;
   }
 
   if (type == 2)  // efficiency_vs_eta: fill "ptrden" with projection of the
                   // plots contained in "den" and fill "ptrnum" with projection
                   // of the plots contained in "num", having btag>m_minTag
   {
     TH2F *numH2 = numME->getTH2F();
     TH2F *denH2 = denME->getTH2F();
 
     TCutG *cutg_num = new TCutG("cutg_num", 4);
     cutg_num->SetPoint(0, m_minTag, -10);
     cutg_num->SetPoint(1, m_minTag, 10);
     cutg_num->SetPoint(2, 1.1, 10);
     cutg_num->SetPoint(3, 1.1, -10);
     ptrnum = numH2->ProjectionY("numerator", 0, -1, "[cutg_num]");
 
     TCutG *cutg_den = new TCutG("cutg_den", 4);
     cutg_den->SetPoint(0, -10.1, -10);
     cutg_den->SetPoint(1, -10.1, 10);
     cutg_den->SetPoint(2, 1.1, 10);
     cutg_den->SetPoint(3, 1.1, -10);
     ptrden = denH2->ProjectionY("denumerator", 0, -1, "[cutg_den]");
     delete cutg_num;
     delete cutg_den;
   }
 
   if (type == 3)  // mistagrate HEP17 / HEM17 relation: fill "ptrnum" with HEP17
                   // mistagrate and "ptrden" with HEM17 mistagrate
   {
     ptrden = denME->getTH1();
     ptrnum = numME->getTH1();
   }
   return true;
 }

void HLTBTagHarvestingAnalyzer::mistagrate	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		TH1F *	num,
		TH1F *	den,
		std::string	effName
	)

Definition at line 277 of file HLTBTagHarvestingAnalyzer.cc.

References DQMStore::IBooker::book1D(), mps_fire::i, and MonitorElement::setEfficiencyFlag().

Referenced by dqmEndJob().

                                                                                                  {
   // do the efficiency_vs_mistag_rate plot
   TH1F *eff;
   eff = new TH1F(effName.c_str(), effName.c_str(), 100, 0, 1);
   eff->SetTitle(effName.c_str());
   eff->SetXTitle("b-effficiency");
   eff->SetYTitle("mistag rate");
   eff->SetOption("E");
   eff->SetLineColor(2);
   eff->SetLineWidth(2);
   eff->SetMarkerStyle(20);
   eff->SetMarkerSize(0.8);
   eff->GetYaxis()->SetRangeUser(0.001, 1.001);
   eff->GetXaxis()->SetRangeUser(-0.001, 1.001);
   eff->SetStats(kFALSE);
 
   // for each bin in the discr -> find efficiency and mistagrate -> put them in
   // a plot
   for (int i = 1; i <= num->GetNbinsX(); i++) {
     double beff = num->GetBinContent(i);
     double miseff = den->GetBinContent(i);
     double miseffErr = den->GetBinError(i);
     int binX = eff->GetXaxis()->FindBin(beff);
     if (eff->GetBinContent(binX) != 0)
       continue;
     eff->SetBinContent(binX, miseff);
     eff->SetBinError(binX, miseffErr);
   }
   MonitorElement *me;
   me = ibooker.book1D(effName, eff);
   me->setEfficiencyFlag();
 
   delete eff;
   return;
 }

void HLTBTagHarvestingAnalyzer::modulesrate	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		TH1F *	num,
		TH1F *	den,
		std::string	effName
	)

Definition at line 314 of file HLTBTagHarvestingAnalyzer.cc.

References DQMStore::IBooker::book1D(), and MonitorElement::setEfficiencyFlag().

Referenced by dqmEndJob().

                                                                                                  {
   // do the eff_vs_disc_HEP17 / eff_vs_disc_HEM17 plot
   TH1F *eff = new TH1F(*num);
   // eff = new TH1F(effName.c_str(),effName.c_str(),100,0,1);
   eff->Divide(den);
   eff->SetTitle(effName.c_str());
   eff->SetXTitle(num->GetXaxis()->GetTitle());
   eff->SetYTitle("");
   eff->SetOption("E");
   eff->SetLineColor(2);
   eff->SetLineWidth(2);
   eff->SetMarkerStyle(20);
   eff->SetMarkerSize(0.8);
   eff->GetYaxis()->SetRangeUser(0.001, 2.001);
   // eff->GetXaxis()->SetRangeUser(-0.001,1.001);
   eff->SetStats(kFALSE);
 
   MonitorElement *me;
   me = ibooker.book1D(effName, eff);
   me->setEfficiencyFlag();
 
   delete eff;
   return;
 }