#include <HLTTauDQML1Plotter.h>

Inheritance diagram for HLTTauDQML1Plotter:

Public Member Functions
void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup, const HLTTauDQMOfflineObjects &refC)

void	bookHistograms (DQMStore::IBooker &iBooker)

	HLTTauDQML1Plotter (const edm::ParameterSet &, edm::ConsumesCollector &&cc, int phibins, double maxpt, double maxhighpt, bool ref, double dr, const std::string &dqmBaseFolder)

	~HLTTauDQML1Plotter ()

Private Attributes
const int	binsEt_

const int	binsEta_

const int	binsPhi_

const bool	doRefAnalysis_

MonitorElement *	firstIsoTauEt_

MonitorElement *	firstIsoTauEta_

MonitorElement *	firstIsoTauPhi_

MonitorElement *	firstTauEt_

MonitorElement *	firstTauEta_

MonitorElement *	firstTauPhi_

MonitorElement *	l1etmEt_

MonitorElement *	l1etmEtEffDenom_

MonitorElement *	l1etmEtEffNum_

double	l1ETMMin_

MonitorElement *	l1etmPhi_

MonitorElement *	l1isotauEt_

MonitorElement *	l1isotauEta_

MonitorElement *	l1isotauEtaEffDenom_

MonitorElement *	l1isotauEtaEffNum_

MonitorElement *	l1isotauEtEffDenom_

MonitorElement *	l1isotauEtEffNum_

MonitorElement *	l1isotauEtRes_

MonitorElement *	l1isotauHighEtEffDenom_

MonitorElement *	l1isotauHighEtEffNum_

MonitorElement *	l1isotauPhi_

MonitorElement *	l1isotauPhiEffDenom_

MonitorElement *	l1isotauPhiEffNum_

MonitorElement *	l1jetEt_

MonitorElement *	l1jetEta_

MonitorElement *	l1jetEtaEffDenom_

MonitorElement *	l1jetEtaEffNum_

MonitorElement *	l1jetEtEffDenom_

MonitorElement *	l1jetEtEffNum_

MonitorElement *	l1jetEtRes_

MonitorElement *	l1jetHighEtEffDenom_

MonitorElement *	l1jetHighEtEffNum_

double	l1JetMinEt_

MonitorElement *	l1jetPhi_

MonitorElement *	l1jetPhiEffDenom_

MonitorElement *	l1jetPhiEffNum_

edm::InputTag	l1stage2Sums_

edm::EDGetTokenT< l1t::EtSumBxCollection >	l1stage2SumsToken_

edm::InputTag	l1stage2Taus_

edm::EDGetTokenT< l1t::TauBxCollection >	l1stage2TausToken_

MonitorElement *	l1tauEt_

MonitorElement *	l1tauEta_

MonitorElement *	l1tauEtaEffDenom_

MonitorElement *	l1tauEtaEffNum_

MonitorElement *	l1tauEtEffDenom_

MonitorElement *	l1tauEtEffNum_

MonitorElement *	l1tauEtRes_

MonitorElement *	l1tauHighEtEffDenom_

MonitorElement *	l1tauHighEtEffNum_

MonitorElement *	l1tauPhi_

MonitorElement *	l1tauPhiEffDenom_

MonitorElement *	l1tauPhiEffNum_

const double	matchDeltaR_

const double	maxEta_

const double	maxHighPt_

const double	maxPt_

MonitorElement *	secondIsoTauEt_

MonitorElement *	secondIsoTauEta_

MonitorElement *	secondIsoTauPhi_

MonitorElement *	secondTauEt_

MonitorElement *	secondTauEta_

MonitorElement *	secondTauPhi_

Private Attributes inherited from HLTTauDQMPlotter
bool	configValid_

std::string	dqmFolder_

std::string	dqmFullFolder_

Additional Inherited Members
Private Member Functions inherited from HLTTauDQMPlotter
	HLTTauDQMPlotter (const edm::ParameterSet &pset, std::string dqmBaseFolder)

	HLTTauDQMPlotter (const std::string &dqmFolder, const std::string &dqmBaseFolder)

bool	isValid () const

	~HLTTauDQMPlotter ()

std::pair< bool, LV >	match (const LV &, const LVColl &, double)

const std::string &	triggerTag () const

Detailed Description

Definition at line 14 of file HLTTauDQML1Plotter.h.

Constructor & Destructor Documentation

HLTTauDQML1Plotter::HLTTauDQML1Plotter	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector &&	cc,
		int	phibins,
		double	maxpt,
		double	maxhighpt,
		bool	ref,
		double	dr,
		const std::string &	dqmBaseFolder
	)

Definition at line 17 of file HLTTauDQML1Plotter.cc.

References HLTTauDQMPlotter::configValid_, edm::ParameterSet::getUntrackedParameter(), l1ETMMin_, l1stage2Sums_, l1stage2SumsToken_, l1stage2Taus_, and l1stage2TausToken_.

                                                                                                                                                                                             :
   HLTTauDQMPlotter(ps, dqmBaseFolder),
   doRefAnalysis_(ref),
   matchDeltaR_(dr),
   maxPt_(maxpt),
   maxHighPt_(maxhighpt),
   binsEt_(ps.getUntrackedParameter<int>("EtHistoBins", 25)),
   binsEta_(ps.getUntrackedParameter<int>("EtaHistoBins", 14)),
   binsPhi_(phibins),
   maxEta_(getMaxEta(binsEta_, ps.getUntrackedParameter<double>("EtaHistoBinWidth", 0.348)))
 {
   if(!configValid_)
     return;
 
   //Process PSet
   l1stage2Taus_     = ps.getUntrackedParameter<edm::InputTag>("L1Taus");
   l1stage2TausToken_= cc.consumes<l1t::TauBxCollection>(l1stage2Taus_);
  
   l1stage2Sums_     = ps.getUntrackedParameter<edm::InputTag>("L1ETM");
   l1stage2SumsToken_= cc.consumes<l1t::EtSumBxCollection>(l1stage2Sums_);
   l1ETMMin_         = ps.getUntrackedParameter<double>("L1ETMMin");
 
   configValid_ = true;
 }

HLTTauDQML1Plotter::~HLTTauDQML1Plotter ( )

default

Referenced by bookHistograms().

Member Function Documentation

void HLTTauDQML1Plotter::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const HLTTauDQMOfflineObjects &	refC
	)

Definition at line 136 of file HLTTauDQML1Plotter.cc.

References a, b, BXVector< T >::begin(), doRefAnalysis_, edm::InputTag::encode(), BXVector< T >::end(), PVValHelper::eta, MonitorElement::Fill(), firstIsoTauEt_, firstIsoTauEta_, firstIsoTauPhi_, firstTauEt_, firstTauEta_, firstTauPhi_, edm::Event::getByToken(), BXVector< T >::getFirstBX(), BXVector< T >::getLastBX(), mps_fire::i, edm::HandleBase::isValid(), L1Analysis::kMissingEt, l1etmEt_, l1etmEtEffDenom_, l1etmEtEffNum_, l1ETMMin_, l1etmPhi_, l1isotauEt_, l1isotauEta_, l1isotauEtaEffDenom_, l1isotauEtaEffNum_, l1isotauEtEffDenom_, l1isotauEtEffNum_, l1isotauEtRes_, l1isotauPhi_, l1isotauPhiEffDenom_, l1isotauPhiEffNum_, l1stage2Sums_, l1stage2SumsToken_, l1stage2Taus_, l1stage2TausToken_, l1tauEt_, l1tauEta_, l1tauEtaEffDenom_, l1tauEtaEffNum_, l1tauEtEffDenom_, l1tauEtEffNum_, l1tauEtRes_, l1tauHighEtEffDenom_, l1tauHighEtEffNum_, l1tauPhi_, l1tauPhiEffDenom_, l1tauPhiEffNum_, funct::m, HLTTauDQMPlotter::match(), matchDeltaR_, HLTTauDQMOfflineObjects::met, phi, edm::Handle< T >::product(), EnergyCorrector::pt, secondIsoTauEt_, secondIsoTauEta_, secondIsoTauPhi_, secondTauEt_, secondTauEta_, secondTauPhi_, BXVector< T >::size(), jetUpdater_cfi::sort, metsig::tau, HLTTauDQMOfflineObjects::taus, and nano_cff::taus.

                                                                                                                            {
     if ( doRefAnalysis_ ) {
         //Tau reference
         for (auto const & tau : refC.taus) {
             l1tauEtEffDenom_->Fill(tau.pt());
             l1tauHighEtEffDenom_->Fill(tau.pt());
             
             l1tauEtaEffDenom_->Fill(tau.eta());
             
             l1tauPhiEffDenom_->Fill(tau.phi());
 
             l1isotauEtEffDenom_->Fill(tau.pt());
             l1isotauEtaEffDenom_->Fill(tau.eta());
             l1isotauPhiEffDenom_->Fill(tau.phi());
         }
     if(!refC.met.empty()) l1etmEtEffDenom_->Fill(refC.met[0].pt());
     }
 
     //Analyze L1 Objects (Tau+Jets)
     edm::Handle<l1t::TauBxCollection> taus;
     iEvent.getByToken(l1stage2TausToken_, taus);  
             
     edm::Handle<l1t::EtSumBxCollection> sums;
     iEvent.getByToken(l1stage2SumsToken_, sums);
 
     LVColl pathTaus;
     LVColl pathIsoTaus;
     
     //Set Variables for the threshold plot
     LVColl l1taus;
     LVColl l1isotaus;
     LVColl l1met;
 
     if(taus.isValid()) {
       for(auto const & i : *taus) {
         l1taus.push_back(i.p4());
         if(i.hwIso() > 0) l1isotaus.push_back(i.p4());
         if(!doRefAnalysis_) {
           l1tauEt_->Fill(i.et());
           l1tauEta_->Fill(i.eta());
           l1tauPhi_->Fill(i.phi());
           pathTaus.push_back(i.p4());
 
           l1isotauEt_->Fill(i.et());
           l1isotauEta_->Fill(i.eta());
           l1isotauPhi_->Fill(i.phi());
           if(i.hwIso() > 0) pathIsoTaus.push_back(i.p4());
         }
       }
     }
     else {
       edm::LogWarning("HLTTauDQMOffline") << "HLTTauDQML1Plotter::analyze: unable to read L1 tau collection " << l1stage2Taus_.encode();
     }
 
     if(sums.isValid() && sums.product()->size() > 0) {
       if(!doRefAnalysis_) {
         for (int ibx = sums->getFirstBX(); ibx <= sums->getLastBX(); ++ibx) {
           for (auto it=sums->begin(ibx); it!=sums->end(ibx); it++) {
             auto type = static_cast<int>( it->getType() );
             if(type == l1t::EtSum::EtSumType::kMissingEt) l1etmEt_->Fill(it->et());
           }
         }
       }
     }
     else {
       edm::LogWarning("HLTTauDQMOffline") << "HLTTauDQML1Plotter::analyze: unable to read L1 met collection " << l1stage2Sums_.encode();
     }
 
     //Now do the efficiency matching
     if ( doRefAnalysis_ ) {
         for (auto const & tau : refC.taus) {
             std::pair<bool,LV> m = match(tau,l1taus,matchDeltaR_);
             if ( m.first ) {
                 l1tauEt_->Fill(m.second.pt());
                 l1tauEta_->Fill(m.second.eta());
                 l1tauPhi_->Fill(m.second.phi());
 
                 l1tauEtEffNum_->Fill(tau.pt());
                 l1tauHighEtEffNum_->Fill(tau.pt());
                 l1tauEtaEffNum_->Fill(tau.eta());
                 l1tauPhiEffNum_->Fill(tau.phi());
 
                 l1tauEtRes_->Fill((m.second.pt()-tau.pt())/tau.pt());
 
                 pathTaus.push_back(m.second);
             }
             m = match(tau,l1isotaus,matchDeltaR_);
             if ( m.first ) {
                 l1isotauEt_->Fill(m.second.pt());
                 l1isotauEta_->Fill(m.second.eta());
                 l1isotauPhi_->Fill(m.second.phi());
 
                 l1isotauEtEffNum_->Fill(tau.pt());
                 l1isotauEtaEffNum_->Fill(tau.eta());
                 l1isotauPhiEffNum_->Fill(tau.phi());
 
                 l1isotauEtRes_->Fill((m.second.pt()-tau.pt())/tau.pt());
 
                 pathIsoTaus.push_back(m.second);
             }
 
         }
 
         if(sums.isValid() && sums.product()->size() > 0) {
           for (int ibx = sums->getFirstBX(); ibx <= sums->getLastBX(); ++ibx) {
             for (auto it=sums->begin(ibx); it!=sums->end(ibx); it++) {
               auto type = static_cast<int>( it->getType() );
               if(type == l1t::EtSum::EtSumType::kMissingEt) {
                 l1etmEt_->Fill(it->et());
                 l1etmPhi_->Fill(it->phi());
         
                 if( it->et() > l1ETMMin_){
                   l1etmEtEffNum_->Fill(it->et());
                 }
               }
             }
           }
         }
     }
     
     //Fill the Threshold Monitoring
     if(pathTaus.size() > 1) std::sort(pathTaus.begin(), pathTaus.end(), [](const LV& a, const LV& b) { return a.pt() > b.pt(); });
     if ( !pathTaus.empty() ) {
         firstTauEt_->Fill(pathTaus[0].pt());
         firstTauEta_->Fill(pathTaus[0].eta());
         firstTauPhi_->Fill(pathTaus[0].phi());
     }
     if ( pathTaus.size() > 1 ) {
         secondTauEt_->Fill(pathTaus[1].pt());
         secondTauEta_->Fill(pathTaus[1].eta());
         secondTauPhi_->Fill(pathTaus[1].phi());
     }
     if(pathIsoTaus.size() > 1) std::sort(pathIsoTaus.begin(), pathIsoTaus.end(), [](const LV& a, const LV& b) { return a.pt() > b.pt(); });
     if ( !pathIsoTaus.empty() ) {
         firstIsoTauEt_->Fill(pathIsoTaus[0].pt());
         firstIsoTauEta_->Fill(pathIsoTaus[0].eta());
         firstIsoTauPhi_->Fill(pathIsoTaus[0].phi());   
     }
     if ( pathIsoTaus.size() > 1 ) {
         secondIsoTauEt_->Fill(pathIsoTaus[1].pt());
         secondIsoTauEta_->Fill(pathIsoTaus[1].eta());  
         secondIsoTauPhi_->Fill(pathIsoTaus[1].phi());
     }
 }

void HLTTauDQML1Plotter::bookHistograms ( DQMStore::IBooker & iBooker )

Definition at line 42 of file HLTTauDQML1Plotter.cc.

References binsEt_, binsEta_, binsPhi_, DQMStore::IBooker::book1D(), edmScanValgrind::buffer, HLTTauDQMPlotter::configValid_, constexpr, doRefAnalysis_, firstIsoTauEt_, firstIsoTauEta_, firstIsoTauPhi_, firstTauEt_, firstTauEta_, firstTauPhi_, l1etmEt_, l1etmEtEffDenom_, l1etmEtEffNum_, l1etmPhi_, l1isotauEt_, l1isotauEta_, l1isotauEtaEffDenom_, l1isotauEtaEffNum_, l1isotauEtEffDenom_, l1isotauEtEffNum_, l1isotauEtRes_, l1isotauPhi_, l1isotauPhiEffDenom_, l1isotauPhiEffNum_, l1tauEt_, l1tauEta_, l1tauEtaEffDenom_, l1tauEtaEffNum_, l1tauEtEffDenom_, l1tauEtEffNum_, l1tauEtRes_, l1tauHighEtEffDenom_, l1tauHighEtEffNum_, l1tauPhi_, l1tauPhiEffDenom_, l1tauPhiEffNum_, maxEta_, maxHighPt_, trackingParticleSelector_cfi::maxPhi, maxPt_, trackingParticleSelector_cfi::minPhi, pi, secondIsoTauEt_, secondIsoTauEta_, secondIsoTauPhi_, secondTauEt_, secondTauEta_, secondTauPhi_, DQMStore::IBooker::setCurrentFolder(), HLTTauDQMPlotter::triggerTag(), and ~HLTTauDQML1Plotter().

                                                                 {
   if(!configValid_)
     return;
 
   // The L1 phi plot is asymmetric around 0 because of the discrete nature of L1 phi
   constexpr float pi = 3.1416f;
   constexpr float phiShift = pi/18; // half of 2pi/18 bin
   constexpr float minPhi = -pi+phiShift;
   constexpr float maxPhi = pi+phiShift;
 
   constexpr int BUFMAX = 256;
   char buffer[BUFMAX] = "";
 
   //Create the histograms
   iBooker.setCurrentFolder(triggerTag());
         
   l1tauEt_ = iBooker.book1D("L1TauEt","L1 #tau E_{T};L1 #tau E_{T};entries",binsEt_,0,maxPt_);
   l1tauEta_ = iBooker.book1D("L1TauEta","L1 #tau #eta;L1 #tau #eta;entries",binsEta_,-maxEta_,maxEta_);
   l1tauPhi_ = iBooker.book1D("L1TauPhi","L1 #tau #phi;L1 #tau #phi;entries",binsPhi_,minPhi,maxPhi);
 
   l1isotauEt_ = iBooker.book1D("L1IsoTauEt","L1 isolated #tau E_{T};L1 #tau E_{T};entries",binsEt_,0,maxPt_);
   l1isotauEta_ = iBooker.book1D("L1IsoTauEta","L1 isolated #tau #eta;L1 #tau #eta;entries",binsEta_,-maxEta_,maxEta_);
   l1isotauPhi_ = iBooker.book1D("L1IsoTauPhi","L1 isolated #tau #phi;L1 #tau #phi;entries",binsPhi_,minPhi,maxPhi);
 
   l1etmEt_  = iBooker.book1D("L1ETM","L1 ETM E_{T};L1 ETM E_{T};entries",binsEt_,0,maxPt_);
   l1etmPhi_ = iBooker.book1D("L1ETMPhi","L1 ETM #phi;L1 ETM #phi;entries",binsPhi_,minPhi,maxPhi);
         
   snprintf(buffer, BUFMAX, "L1 leading #tau E_{T};L1 #tau E_{T};entries");
   firstTauEt_ = iBooker.book1D("L1LeadTauEt", buffer, binsEt_, 0, maxPt_);
   snprintf(buffer, BUFMAX, "L1 leading #tau #eta;L1 #tau #eta;entries");
   firstTauEta_ = iBooker.book1D("L1LeadTauEta", buffer, binsEta_, -maxEta_, maxEta_);
   snprintf(buffer, BUFMAX, "L1 leading #tau #phi;L1 #tau #phi;entries");
   firstTauPhi_ = iBooker.book1D("L1LeadTauPhi", buffer, binsPhi_, minPhi, maxPhi);
         
   snprintf(buffer, BUFMAX, "L1 second-leading #tau E_{T};L1 #tau E_{T};entries");
   secondTauEt_ = iBooker.book1D("L1SecondTauEt", buffer, binsEt_, 0, maxPt_);
   snprintf(buffer, BUFMAX, "L1 second-leading #tau #eta;L1 #tau #eta;entries");
   secondTauEta_ = iBooker.book1D("L1SecondTauEta", buffer, binsEta_, -maxEta_, maxEta_);
   snprintf(buffer, BUFMAX, "L1 second-leading #tau #phi;L1 #tau #phi;entries");
   secondTauPhi_ = iBooker.book1D("L1SecondTauPhi", buffer, binsPhi_, minPhi, maxPhi);
 
   snprintf(buffer, BUFMAX, "L1 leading isolated #tau E_{T};L1 #tau E_{T};entries");
   firstIsoTauEt_ = iBooker.book1D("L1LeadIsoTauEt", buffer, binsEt_, 0, maxPt_);
   snprintf(buffer, BUFMAX, "L1 leading isolated #tau #eta;L1 #tau #eta;entries");
   firstIsoTauEta_ = iBooker.book1D("L1LeadIsoTauEta", buffer, binsEta_, -maxEta_, maxEta_);
   snprintf(buffer, BUFMAX, "L1 leading isolated #tau #phi;L1 #tau #phi;entries");
   firstIsoTauPhi_ = iBooker.book1D("L1LeadIsoTauPhi", buffer, binsPhi_, minPhi, maxPhi);
 
   snprintf(buffer, BUFMAX, "L1 second-leading isolated #tau E_{T};L1 #tau E_{T};entries");
   secondIsoTauEt_ = iBooker.book1D("L1SecondIsoTauEt", buffer, binsEt_, 0, maxPt_);
   snprintf(buffer, BUFMAX, "L1 second-leading isolated #tau #eta;L1 #tau #eta;entries");
   secondIsoTauEta_ = iBooker.book1D("L1SecondIsoTauEta", buffer, binsEta_, -maxEta_, maxEta_);
   snprintf(buffer, BUFMAX, "L1 second-leading isolated #tau #phi;L1 #tau #phi;entries");
   secondIsoTauPhi_ = iBooker.book1D("L1SecondIsoTauPhi", buffer, binsPhi_, minPhi, maxPhi);
         
   if (doRefAnalysis_) {
     l1tauEtRes_ = iBooker.book1D("L1TauEtResol","L1 #tau E_{T} resolution;[L1 #tau E_{T}-Ref #tau E_{T}]/Ref #tau E_{T};entries",60,-1,4);
     l1isotauEtRes_ = iBooker.book1D("L1IsoTauEtResol","L1 isolated #tau E_{T} resolution;[L1 #tau E_{T}-Ref #tau E_{T}]/Ref #tau E_{T};entries",60,-1,4);            
 
     iBooker.setCurrentFolder(triggerTag()+"/helpers");
             
     l1tauEtEffNum_ = iBooker.book1D("L1TauEtEffNum","L1 #tau E_{T} Efficiency;Ref #tau E_{T};entries",binsEt_,0,maxPt_);
     l1tauHighEtEffNum_ = iBooker.book1D("L1TauHighEtEffNum","L1 #tau E_{T} Efficiency (high E_{T});Ref #tau E_{T};entries",binsEt_,0,maxHighPt_);
             
     l1tauEtEffDenom_ = iBooker.book1D("L1TauEtEffDenom","L1 #tau E_{T} Denominator;Ref #tau E_{T};entries",binsEt_,0,maxPt_);
     l1tauHighEtEffDenom_ = iBooker.book1D("L1TauHighEtEffDenom","L1 #tau E_{T} Denominator (high E_{T});Ref #tau E_{T};Efficiency",binsEt_,0,maxHighPt_);
             
     l1tauEtaEffNum_ = iBooker.book1D("L1TauEtaEffNum","L1 #tau #eta Efficiency;Ref #tau #eta;entries",binsEta_,-maxEta_,maxEta_);
     l1tauEtaEffDenom_ = iBooker.book1D("L1TauEtaEffDenom","L1 #tau #eta Denominator;Ref #tau #eta;entries",binsEta_,-maxEta_,maxEta_);
             
     l1tauPhiEffNum_ = iBooker.book1D("L1TauPhiEffNum","L1 #tau #phi Efficiency;Ref #tau #phi;entries",binsPhi_,minPhi,maxPhi);
     l1tauPhiEffDenom_ = iBooker.book1D("L1TauPhiEffDenom","L1 #tau #phi Denominator;Ref #tau #phi;Efficiency",binsPhi_,minPhi,maxPhi);
 
     l1isotauEtEffNum_ = iBooker.book1D("L1IsoTauEtEffNum","L1 isolated #tau E_{T} Efficiency;Ref #tau E_{T};entries",binsEt_,0,maxPt_);
     l1isotauEtEffDenom_ = iBooker.book1D("L1IsoTauEtEffDenom","L1 isolated #tau E_{T} Denominator;Ref #tau E_{T};entries",binsEt_,0,maxPt_);
 
     l1isotauEtaEffNum_ = iBooker.book1D("L1IsoTauEtaEffNum","L1 isolated #tau #eta Efficiency;Ref #tau #eta;entries",binsEta_,-maxEta_,maxEta_);
     l1isotauEtaEffDenom_ = iBooker.book1D("L1IsoTauEtaEffDenom","L1 isolated #tau #eta Denominator;Ref #tau #eta;entries",binsEta_,-maxEta_,maxEta_);
     
     l1isotauPhiEffNum_ = iBooker.book1D("L1IsoTauPhiEffNum","L1 isolated #tau #phi Efficiency;Ref #tau #phi;entries",binsPhi_,minPhi,maxPhi);
     l1isotauPhiEffDenom_ = iBooker.book1D("L1IsoTauPhiEffDenom","L1 isolated #tau #phi Denominator;Ref #tau #phi;Efficiency",binsPhi_,minPhi,maxPhi);
 
     l1etmEtEffNum_ = iBooker.book1D("L1ETMEtEffNum", "L1 ETM Efficiency;Ref MET;entries",binsEt_, 0, maxPt_);
     l1etmEtEffDenom_ = iBooker.book1D("L1ETMEtEffDenom","L1 ETM Denominator;Ref MET;entries",binsEt_,0,maxPt_);
   }
 }