#include <TriggerMatchEfficiencyPlotter.h>

Inheritance diagram for TriggerMatchEfficiencyPlotter:

Public Member Functions
	TriggerMatchEfficiencyPlotter (const edm::ParameterSet &ps)
	Constructor. More...

	~TriggerMatchEfficiencyPlotter () override
	Destructor. More...

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

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)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
	DQM Client Diagnostic. More...

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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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)

Private Attributes
MonitorElement *	h_eff_Path1_eta_tight

MonitorElement *	h_eff_Path1_phi_tight

MonitorElement *	h_eff_Path1_pt_tight

MonitorElement *	h_eff_Path2_eta_tight

MonitorElement *	h_eff_Path2_phi_tight

MonitorElement *	h_eff_Path2_pt_tight

edm::ParameterSet	parameters

DQMStore *	theDbe

std::string	theFolder

std::string	triggerhistName1_

std::string	triggerhistName2_

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 Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 31 of file TriggerMatchEfficiencyPlotter.h.

Constructor & Destructor Documentation

TriggerMatchEfficiencyPlotter::TriggerMatchEfficiencyPlotter ( const edm::ParameterSet & ps )

Constructor.

Definition at line 34 of file TriggerMatchEfficiencyPlotter.cc.

References gather_cfg::cout, and edm::ParameterSet::getParameter().

                                                                                      {
 #ifdef DEBUG
   cout << "TriggerMatchEfficiencyPlotter(): Constructor " << endl;
 #endif
   parameters = ps;
 
   triggerhistName1_ = parameters.getParameter<string>("triggerhistName1");
   triggerhistName2_ = parameters.getParameter<string>("triggerhistName2"); 
   theFolder = parameters.getParameter<string>("folder");
 }

TriggerMatchEfficiencyPlotter::~TriggerMatchEfficiencyPlotter ( )

override

Destructor.

Definition at line 44 of file TriggerMatchEfficiencyPlotter.cc.

44 {}

Member Function Documentation

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

overrideprotectedvirtual

DQM Client Diagnostic.

— Tight Muon trigger mathcing efficiency

Implements DQMEDHarvester.

Definition at line 46 of file TriggerMatchEfficiencyPlotter.cc.

References DQMStore::IBooker::book1D(), DQMStore::IGetter::get(), MonitorElement::getTH1F(), specificLumi::inputdir, DQMStore::IBooker::setCurrentFolder(), and MonitorElement::setEfficiencyFlag().

                                                                                                   {
   
   ibooker.setCurrentFolder(theFolder);
   
   // efficiency plot
   h_eff_Path1_eta_tight      = ibooker.book1D("matchingEff_eta_"+triggerhistName1_+"_tight",triggerhistName1_+":matching Eff. vs #eta",8, -2.5, 2.5);
   h_eff_Path1_pt_tight      = ibooker.book1D("matchingEff_pt_"+triggerhistName1_+"_tight",triggerhistName1_+":matching Eff. vs pt",10, 20, 220);
   h_eff_Path1_phi_tight      = ibooker.book1D("matchingEff_phi_"+triggerhistName1_+"_tight",triggerhistName1_+":matching Eff. vs #phi", 8, -3.0, 3.0);
   h_eff_Path2_eta_tight      = ibooker.book1D("matchingEff_eta_"+triggerhistName2_+"_tight",triggerhistName2_+":matching Eff. vs #eta",8, -2.5, 2.5);
   h_eff_Path2_pt_tight      = ibooker.book1D("matchingEff_pt_"+triggerhistName2_+"_tight",triggerhistName2_+":matching Eff. vs pt",10, 20, 220);
   h_eff_Path2_phi_tight      = ibooker.book1D("matchingEff_phi_"+triggerhistName2_+"_tight",triggerhistName2_+":matching Eff. vs #phi", 8, -3.0, 3.0);
 
 
 
   // This prevents this ME to be normalized when drawn into the GUI
   h_eff_Path1_eta_tight->setEfficiencyFlag();
   h_eff_Path1_pt_tight->setEfficiencyFlag();
   h_eff_Path1_phi_tight->setEfficiencyFlag();
   h_eff_Path2_eta_tight->setEfficiencyFlag();
   h_eff_Path2_pt_tight->setEfficiencyFlag();
   h_eff_Path2_phi_tight->setEfficiencyFlag();
 
   // AXIS TITLES....
   h_eff_Path1_eta_tight      ->setAxisTitle("#eta",         1);  
   h_eff_Path1_pt_tight       ->setAxisTitle("pt",         1);
   h_eff_Path1_phi_tight       ->setAxisTitle("#phi",         1);  
   h_eff_Path2_eta_tight      ->setAxisTitle("#eta",         1);
   h_eff_Path2_pt_tight       ->setAxisTitle("pt",         1);
   h_eff_Path2_phi_tight       ->setAxisTitle("#phi",         1);
 
 
 
   string inputdir = "Muons_miniAOD/TriggerMatchMonitor/EfficiencyInput";
   string numpath_eta_path1 = inputdir+"/passHLT"+triggerhistName1_+"_eta_Tight";
   string denpath_eta_path1 = inputdir+"/totalHLT"+triggerhistName1_+"_eta_Tight";
 
   string numpath_pt_path1 = inputdir+"/passHLT"+triggerhistName1_+"_pt_Tight";
   string denpath_pt_path1 = inputdir+"/totalHLT"+triggerhistName1_+"_pt_Tight";
   
   string numpath_phi_path1 = inputdir+"/passHLT"+triggerhistName1_+"_phi_Tight";
   string denpath_phi_path1 = inputdir+"/totalHLT"+triggerhistName1_+"_phi_Tight"; 
 
   string numpath_eta_path2 = inputdir+"/passHLT"+triggerhistName2_+"_eta_Tight";
   string denpath_eta_path2 = inputdir+"/totalHLT"+triggerhistName2_+"_eta_Tight";
 
   string numpath_pt_path2 = inputdir+"/passHLT"+triggerhistName2_+"_pt_Tight";
   string denpath_pt_path2 = inputdir+"/totalHLT"+triggerhistName2_+"_pt_Tight";
 
   string numpath_phi_path2 = inputdir+"/passHLT"+triggerhistName2_+"_phi_Tight";
   string denpath_phi_path2 = inputdir+"/totalHLT"+triggerhistName2_+"_phi_Tight";
   
   MonitorElement *Numerator_eta_path1   = igetter.get(numpath_eta_path1);
   MonitorElement *Denominator_eta_path1 = igetter.get(denpath_eta_path1);
 
   MonitorElement *Numerator_pt_path1   = igetter.get(numpath_pt_path1);
   MonitorElement *Denominator_pt_path1 = igetter.get(denpath_pt_path1);
 
   MonitorElement *Numerator_phi_path1   = igetter.get(numpath_phi_path1);
   MonitorElement *Denominator_phi_path1 = igetter.get(denpath_phi_path1);
 
   MonitorElement *Numerator_eta_path2   = igetter.get(numpath_eta_path2);
   MonitorElement *Denominator_eta_path2 = igetter.get(denpath_eta_path2);
 
   MonitorElement *Numerator_pt_path2   = igetter.get(numpath_pt_path2);
   MonitorElement *Denominator_pt_path2 = igetter.get(denpath_pt_path2);
 
   MonitorElement *Numerator_phi_path2   = igetter.get(numpath_phi_path2);
   MonitorElement *Denominator_phi_path2 = igetter.get(denpath_phi_path2); 
 
   if (Numerator_eta_path1 && Denominator_eta_path1){
     TH1F *h_numerator_eta_path1   = Numerator_eta_path1->getTH1F();
     TH1F *h_denominator_eta_path1 = Denominator_eta_path1->getTH1F();
     TH1F *h_eff_eta_path1         = h_eff_Path1_eta_tight->getTH1F();
     
     if (h_eff_eta_path1->GetSumw2N() == 0) h_eff_eta_path1->Sumw2();  
     h_eff_eta_path1->Divide(h_numerator_eta_path1, h_denominator_eta_path1, 1., 1., "B");
   }
 
   if (Numerator_pt_path1 && Denominator_pt_path1){
     TH1F *h_numerator_pt_path1   = Numerator_pt_path1->getTH1F();
     TH1F *h_denominator_pt_path1 = Denominator_pt_path1->getTH1F();
     TH1F *h_eff_pt_path1         = h_eff_Path1_pt_tight->getTH1F();
 
     if (h_eff_pt_path1->GetSumw2N() == 0) h_eff_pt_path1->Sumw2();
     h_eff_pt_path1->Divide(h_numerator_pt_path1, h_denominator_pt_path1, 1., 1., "B");
   }
 
   if (Numerator_phi_path1 && Denominator_phi_path1){
     TH1F *h_numerator_phi_path1   = Numerator_phi_path1->getTH1F();
     TH1F *h_denominator_phi_path1 = Denominator_phi_path1->getTH1F();
     TH1F *h_eff_phi_path1         = h_eff_Path1_phi_tight->getTH1F();
 
     if (h_eff_phi_path1->GetSumw2N() == 0) h_eff_phi_path1->Sumw2();
     h_eff_phi_path1->Divide(h_numerator_phi_path1, h_denominator_phi_path1, 1., 1., "B");
   }
  
   //trigger path2
   if (Numerator_eta_path2 && Denominator_eta_path2){
     TH1F *h_numerator_eta_path2   = Numerator_eta_path2->getTH1F();
     TH1F *h_denominator_eta_path2 = Denominator_eta_path2->getTH1F();
     TH1F *h_eff_eta_path2         = h_eff_Path2_eta_tight->getTH1F();
 
     if (h_eff_eta_path2->GetSumw2N() == 0) h_eff_eta_path2->Sumw2();
     h_eff_eta_path2->Divide(h_numerator_eta_path2, h_denominator_eta_path2, 1., 1., "B");
   }
 
   if (Numerator_pt_path2 && Denominator_pt_path2){
     TH1F *h_numerator_pt_path2   = Numerator_pt_path2->getTH1F();
     TH1F *h_denominator_pt_path2 = Denominator_pt_path2->getTH1F();
     TH1F *h_eff_pt_path2         = h_eff_Path2_pt_tight->getTH1F();
 
     if (h_eff_pt_path2->GetSumw2N() == 0) h_eff_pt_path2->Sumw2();
     h_eff_pt_path2->Divide(h_numerator_pt_path2, h_denominator_pt_path2, 1., 1., "B");
   }
 
   if (Numerator_phi_path2 && Denominator_phi_path2){
     TH1F *h_numerator_phi_path2   = Numerator_phi_path2->getTH1F();
     TH1F *h_denominator_phi_path2 = Denominator_phi_path2->getTH1F();
     TH1F *h_eff_phi_path2         = h_eff_Path2_phi_tight->getTH1F();
 
     if (h_eff_phi_path2->GetSumw2N() == 0) h_eff_phi_path2->Sumw2();
     h_eff_phi_path2->Divide(h_numerator_phi_path2, h_denominator_phi_path2, 1., 1., "B");
   }
 
 
 }