#include <L1TEfficiencyMuons_Offline.h>

Inheritance diagram for L1TEfficiencyMuons_Offline:

Public Member Functions
	L1TEfficiencyMuons_Offline (const edm::ParameterSet &ps)

virtual	~L1TEfficiencyMuons_Offline ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c) override

virtual void	bookControlHistos (DQMStore::IBooker &)

virtual void	bookEfficiencyHistos (DQMStore::IBooker &ibooker, int ptCut)

virtual void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &iSetup) override

virtual void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c)

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 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 Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	getMuonGmtPairs (edm::Handle< l1t::MuonBxCollection > &gmtCands)

const reco::Vertex	getPrimaryVertex (edm::Handle< reco::VertexCollection > &vertex, edm::Handle< reco::BeamSpot > &beamSpot)

void	getProbeMuons (edm::Handle< edm::TriggerResults > &trigResults, edm::Handle< trigger::TriggerEvent > &trigEvent)

void	getTightMuons (edm::Handle< reco::MuonCollection > &muons, const reco::Vertex &vertex)

bool	matchHlt (edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon *mu)

Private Attributes
edm::ESHandle< MagneticField >	m_BField

edm::EDGetTokenT< reco::BeamSpot >	m_BsInputTag

std::map< std::string, MonitorElement * >	m_ControlHistos

std::map< int, std::map< std::string, MonitorElement * > >	m_EfficiencyHistos

edm::EDGetTokenT< l1t::MuonBxCollection >	m_GmtInputTag

std::vector< int >	m_GmtPtCuts

HLTConfigProvider	m_hltConfig

BXVector< l1t::Muon >	m_L1tL1tMuons

std::vector< l1t::MuonBxCollection >	m_L1tMuons

float	m_MaxGmtMuonDR

float	m_MaxHltMuonDR

float	m_MaxMuonEta

std::vector< MuonGmtPair >	m_MuonGmtPairs

edm::EDGetTokenT< reco::MuonCollection >	m_MuonInputTag

std::vector< const reco::Muon * >	m_ProbeMuons

edm::ESHandle< Propagator >	m_propagatorAlong

edm::ESHandle< Propagator >	m_propagatorOpposite

std::vector< reco::MuonCollection >	m_RecoMuons

std::vector< reco::Muon >	m_RecoRecoMuons

std::vector< const reco::Muon * >	m_TightMuons

std::vector< int >	m_trigIndices

edm::EDGetTokenT< trigger::TriggerEvent >	m_trigInputTag

std::vector< std::string >	m_trigNames

std::string	m_trigProcess

edm::EDGetTokenT< edm::TriggerResults >	m_trigProcess_token

bool	m_verbose

edm::EDGetTokenT< reco::VertexCollection >	m_VtxInputTag

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 113 of file L1TEfficiencyMuons_Offline.h.

Constructor & Destructor Documentation

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

Definition at line 107 of file L1TEfficiencyMuons_Offline.cc.

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

                                                                              {
     m_verbose = ps.getUntrackedParameter<bool>("verbose");
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] ____________ Storage initialization ____________ " << endl;
 
     // Initializing config params
     m_GmtPtCuts = ps.getUntrackedParameter< vector<int> >("gmtPtCuts");
     m_MuonInputTag =  consumes<reco::MuonCollection>(ps.getUntrackedParameter<InputTag>("muonInputTag"));
     m_GmtInputTag  =  consumes<l1t::MuonBxCollection>(ps.getUntrackedParameter<InputTag>("gmtInputTag"));
     m_VtxInputTag =  consumes<VertexCollection>(ps.getUntrackedParameter<InputTag>("vtxInputTag"));
     m_BsInputTag  =  consumes<BeamSpot>(ps.getUntrackedParameter<InputTag>("bsInputTag"));
     m_trigInputTag = consumes<trigger::TriggerEvent>(ps.getUntrackedParameter<InputTag>("trigInputTag"));
     m_trigProcess  = ps.getUntrackedParameter<string>("trigProcess");
     m_trigProcess_token  = consumes<edm::TriggerResults>(ps.getUntrackedParameter<InputTag>("trigProcess_token"));
     m_trigNames    = ps.getUntrackedParameter<vector<string> >("triggerNames");
 
     // CB do we need them from cfi?
     m_MaxMuonEta   = 2.4;
     m_MaxGmtMuonDR = 0.7;
     m_MaxHltMuonDR = 0.1;
     // CB ignored at present
     //m_MinMuonDR    = 1.2;
 }

L1TEfficiencyMuons_Offline::~L1TEfficiencyMuons_Offline ( )

virtual

Definition at line 131 of file L1TEfficiencyMuons_Offline.cc.

131 { }

Member Function Documentation

void L1TEfficiencyMuons_Offline::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprivate

Definition at line 180 of file L1TEfficiencyMuons_Offline.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, BXVector< T >::begin(), gather_cfg::cout, BXVector< T >::end(), stringResolutionProvider_cfi::eta, edm::EventSetup::get(), edm::Event::getByToken(), RPCpg::mu, electronCleaner_cfi::muons, impactParameterTagInfos_cfi::primaryVertex, and EnergyCorrector::pt.

                                                                                            {
 
     Handle<reco::MuonCollection> muons;
     iEvent.getByToken(m_MuonInputTag, muons);
 
     Handle<BeamSpot> beamSpot;
     iEvent.getByToken(m_BsInputTag, beamSpot);
 
     Handle<VertexCollection> vertex;
     iEvent.getByToken(m_VtxInputTag, vertex);
 
     Handle<l1t::MuonBxCollection> gmtCands;
     iEvent.getByToken(m_GmtInputTag,gmtCands);
 
     Handle<edm::TriggerResults> trigResults;
     iEvent.getByToken(m_trigProcess_token,trigResults);
 
     edm::Handle<trigger::TriggerEvent> trigEvent;
     iEvent.getByToken(m_trigInputTag,trigEvent);
 
     eventSetup.get<IdealMagneticFieldRecord>().get(m_BField);
     eventSetup.get<TrackingComponentsRecord>().get("PropagatorWithMaterial",m_propagatorAlong);
     eventSetup.get<TrackingComponentsRecord>().get("PropagatorWithMaterialOpposite",m_propagatorOpposite);
 
     const Vertex primaryVertex = getPrimaryVertex(vertex,beamSpot);
 
     getTightMuons(muons,primaryVertex);
     getProbeMuons(trigResults,trigEvent); // CB add flag to run on orthogonal datasets (no T&P)
     getMuonGmtPairs(gmtCands);
 
 //    MuonCollection::const_iterator muonIt  = muons->begin();
 //    MuonCollection::const_iterator muonEnd = muons->end();
 
     vector<l1t::Muon> gmtContainer;// = gmtCands->getRecord(0).getGMTCands();                       
 
     for (auto mu = gmtCands->begin(0); mu != gmtCands->end(0); ++mu) {
         gmtContainer.push_back(*mu);
     }
 
 //    vector<l1t::Muon>::const_iterator gmtIt = gmtContainer.begin();
 //    vector<l1t::Muon>::const_iterator gmtEnd = gmtContainer.end();
 
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Computing efficiencies" << endl;
 
     vector<MuonGmtPair>::const_iterator muonGmtPairsIt  = m_MuonGmtPairs.begin();
     vector<MuonGmtPair>::const_iterator muonGmtPairsEnd = m_MuonGmtPairs.end();
 
     for(; muonGmtPairsIt!=muonGmtPairsEnd; ++muonGmtPairsIt) {
         float eta = muonGmtPairsIt->eta();
         float phi = muonGmtPairsIt->phi();
         float pt  = muonGmtPairsIt->pt();
 
         float gmtPt  = muonGmtPairsIt->gmtPt();
         int qual = muonGmtPairsIt->gmtQual();
 
         vector<int>::const_iterator gmtPtCutsIt  = m_GmtPtCuts.begin();
         vector<int>::const_iterator gmtPtCutsEnd = m_GmtPtCuts.end();
 
         for (; gmtPtCutsIt!=gmtPtCutsEnd; ++ gmtPtCutsIt) {
             int gmtPtCut = (*gmtPtCutsIt);
             bool gmtAboveCut = (gmtPt > gmtPtCut);
 
             stringstream ptCutToTag; ptCutToTag << gmtPtCut;
             string ptTag = ptCutToTag.str();
 
             if (fabs(eta) < m_MaxMuonEta) {
                 m_EfficiencyHistos[gmtPtCut]["EffvsPt" + ptTag + "Den"]->Fill(pt);
                 m_EfficiencyHistos[gmtPtCut]["EffvsPt_OPEN_" + ptTag + "Den"]->Fill(pt);
                 m_EfficiencyHistos[gmtPtCut]["EffvsPt_DOUBLE_" + ptTag + "Den"]->Fill(pt);
                 m_EfficiencyHistos[gmtPtCut]["EffvsPt_SINGLE_" + ptTag + "Den"]->Fill(pt);
 
                 if (gmtAboveCut) {
                     m_EfficiencyHistos[gmtPtCut]["EffvsPt" + ptTag + "Num"]->Fill(pt);
 
                     if (qual >= 4) m_EfficiencyHistos[gmtPtCut]["EffvsPt_OPEN_" + ptTag + "Num"]->Fill(pt);
                     if (qual >= 8) m_EfficiencyHistos[gmtPtCut]["EffvsPt_DOUBLE_" + ptTag + "Num"]->Fill(pt);
                     if (qual >= 12) m_EfficiencyHistos[gmtPtCut]["EffvsPt_SINGLE_" + ptTag + "Num"]->Fill(pt);
                 }
 
                 // efficiency in eta/phi at plateau
                 if (pt > 1.25*gmtPtCut) {       // efficiency in eta/phi at plateau
 
                     m_EfficiencyHistos[gmtPtCut]["EffvsPhi" + ptTag + "Den"]->Fill(phi);
                     m_EfficiencyHistos[gmtPtCut]["EffvsEta" + ptTag + "Den"]->Fill(eta);
 
                     m_EfficiencyHistos[gmtPtCut]["EffvsPhi_OPEN_" + ptTag + "Den"]->Fill(phi);
                     m_EfficiencyHistos[gmtPtCut]["EffvsEta_OPEN_" + ptTag + "Den"]->Fill(eta);
 
                     m_EfficiencyHistos[gmtPtCut]["EffvsPhi_DOUBLE_" + ptTag + "Den"]->Fill(phi);
                     m_EfficiencyHistos[gmtPtCut]["EffvsEta_DOUBLE_" + ptTag + "Den"]->Fill(eta);
 
                     m_EfficiencyHistos[gmtPtCut]["EffvsPhi_SINGLE_" + ptTag + "Den"]->Fill(phi);
                     m_EfficiencyHistos[gmtPtCut]["EffvsEta_SINGLE_" + ptTag + "Den"]->Fill(eta);
 
                     if (gmtAboveCut) {
                         m_EfficiencyHistos[gmtPtCut]["EffvsPhi" + ptTag + "Num"]->Fill(phi);
                         m_EfficiencyHistos[gmtPtCut]["EffvsEta" + ptTag + "Num"]->Fill(eta);
 
                         if (qual >= 4) {
                             m_EfficiencyHistos[gmtPtCut]["EffvsPhi_OPEN_" + ptTag + "Num"]->Fill(phi);
                             m_EfficiencyHistos[gmtPtCut]["EffvsEta_OPEN_" + ptTag + "Num"]->Fill(eta);
                         }
                         if (qual >= 8) {
                             m_EfficiencyHistos[gmtPtCut]["EffvsPhi_DOUBLE_" + ptTag + "Num"]->Fill(phi);
                             m_EfficiencyHistos[gmtPtCut]["EffvsEta_DOUBLE_" + ptTag + "Num"]->Fill(eta);
                         }
                         if (qual >= 12) {
                             m_EfficiencyHistos[gmtPtCut]["EffvsPhi_SINGLE_" + ptTag + "Num"]->Fill(phi);
                             m_EfficiencyHistos[gmtPtCut]["EffvsEta_SINGLE_" + ptTag + "Num"]->Fill(eta);
                         }
                     }
                 }
             }
         }
     }
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Computation finished" << endl;
 }

void L1TEfficiencyMuons_Offline::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	c
	)

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 170 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, edm::LuminosityBlockBase::id(), and edm::LuminosityBlockID::luminosityBlock().

                                                                                                            {
     if(m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Called beginLuminosityBlock at LS=" << lumiBlock.id().luminosityBlock() << endl;
 }

void L1TEfficiencyMuons_Offline::bookControlHistos ( DQMStore::IBooker & ibooker )

protectedvirtual

Definition at line 299 of file L1TEfficiencyMuons_Offline.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, dataset::name, Pi, and DQMStore::IBooker::setCurrentFolder().

                                                                            {
     if(m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Booking Control Plot Histos" << endl;
 
     ibooker.setCurrentFolder("L1T/Efficiency/Muons/Control");
 
     string name = "MuonGmtDeltaR";
     m_ControlHistos[name] = ibooker.book1D(name.c_str(),name.c_str(),25.,0.,2.5);
 
     name = "NTightVsAll";
     m_ControlHistos[name] = ibooker.book2D(name.c_str(),name.c_str(),5,-0.5,4.5,5,-0.5,4.5);
 
     name = "NProbesVsTight";
     m_ControlHistos[name] = ibooker.book2D(name.c_str(),name.c_str(),5,-0.5,4.5,5,-0.5,4.5);
 
     string name1 = "TagMuonEta_Histo";
     m_EfficiencyHistos[0][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
     string name2 = "TagMuonPhi_Histo";
     m_EfficiencyHistos[0][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
     string name3 = "TagMuonPt_Histo";
     m_EfficiencyHistos[0][name3] = ibooker.book1D(name3.c_str(),name3.c_str(),50,0.,100.);
 //*****
     name1 = "ProbeMuonEta_Histo";
     m_EfficiencyHistos[0][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
     name2 = "ProbeMuonPhi_Histo";
     m_EfficiencyHistos[0][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
     name3 = "ProbeMuonPt_Histo";
     m_EfficiencyHistos[0][name3] = ibooker.book1D(name3.c_str(),name3.c_str(),50,0.,100.);
 }

void L1TEfficiencyMuons_Offline::bookEfficiencyHistos	(	DQMStore::IBooker &	ibooker,
		int	ptCut
	)

protectedvirtual

Definition at line 330 of file L1TEfficiencyMuons_Offline.cc.

References DQMStore::IBooker::book1D(), gather_cfg::cout, Pi, HiGenCleaner_cff::ptCut, DQMStore::IBooker::setCurrentFolder(), and fw3dlego::xbins.

                                                                                          {
     if(m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Booking Efficiency Plot Histos for pt cut = " << ptCut << endl;
 
     stringstream ptCutToTag; ptCutToTag << ptCut;
     string ptTag = ptCutToTag.str();
 
     ibooker.setCurrentFolder("L1T/Efficiency/Muons/");
     float xbins[33] = {0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 55, 60, 65, 70, 80, 90, 100};
 
     string name1 = "EffvsPt" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),32,xbins);
     string name2 = "EffvsPt" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),32,xbins);
 
     name1 = "EffvsPt_OPEN_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),32,xbins);
     name2 = "EffvsPt_OPEN_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),32,xbins);
 
     name1 = "EffvsPt_DOUBLE_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),32, xbins);
     name2 = "EffvsPt_DOUBLE_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),32, xbins);
 
     name1 = "EffvsPt_SINGLE_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),32,xbins);
     name2 = "EffvsPt_SINGLE_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),32,xbins);
 
 
     name1 = "EffvsPhi" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),24,-TMath::Pi(),TMath::Pi());
     name2 = "EffvsPhi" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
 
     name1 = "EffvsEta" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
     name2 = "EffvsEta" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),50,-2.5,2.5);
 
 
     name1 = "EffvsPhi_OPEN_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),24,-TMath::Pi(),TMath::Pi());
     name2 = "EffvsPhi_OPEN_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
 
     name1 = "EffvsEta_OPEN_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
     name2 = "EffvsEta_OPEN_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),50,-2.5,2.5);
 
 
     name1 = "EffvsPhi_DOUBLE_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),24,-TMath::Pi(),TMath::Pi());
     name2 = "EffvsPhi_DOUBLE_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
 
     name1 = "EffvsEta_DOUBLE_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
     name2 = "EffvsEta_DOUBLE_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),50,-2.5,2.5);
 
 
     name1 = "EffvsPhi_SINGLE_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),24,-TMath::Pi(),TMath::Pi());
     name2 = "EffvsPhi_SINGLE_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
 
 
     name1 = "EffvsEta_SINGLE_" + ptTag + "Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
     name2 = "EffvsEta_SINGLE_" + ptTag + "Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),50,-2.5,2.5);
 }

void L1TEfficiencyMuons_Offline::bookHistograms	(	DQMStore::IBooker &	ibooker,
		const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Definition at line 140 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, and createfilelist::int.

                                                                                                                        {
     //book histos
     bookControlHistos(ibooker);
     vector<int>::const_iterator gmtPtCutsIt  = m_GmtPtCuts.begin();
     vector<int>::const_iterator gmtPtCutsEnd = m_GmtPtCuts.end();
 
     for (; gmtPtCutsIt!=gmtPtCutsEnd; ++ gmtPtCutsIt) {
         bookEfficiencyHistos(ibooker, (*gmtPtCutsIt));
     }
 
     vector<string>::const_iterator trigNamesIt  = m_trigNames.begin();
     vector<string>::const_iterator trigNamesEnd = m_trigNames.end();
 
     for (; trigNamesIt!=trigNamesEnd; ++trigNamesIt) {
         TString tNameTmp = TString(*trigNamesIt); // use TString as it handles regex
         TRegexp tNamePattern = TRegexp(tNameTmp,true);
         int tIndex = -1;
 
         for (unsigned ipath = 0; ipath < m_hltConfig.size(); ++ipath) {
             TString tmpName = TString(m_hltConfig.triggerName(ipath));
             if (tmpName.Contains(tNamePattern)) {
                 tIndex = int(ipath);
                 m_trigIndices.push_back(tIndex);
             }
         }
         if (tIndex < 0 && m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Warning: Could not find trigger " << (*trigNamesIt) << endl;   
     }
 }

void L1TEfficiencyMuons_Offline::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Definition at line 133 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout.

                                                                                           {
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Called beginRun." << endl;
     bool changed = true;
     m_hltConfig.init(run,iSetup,m_trigProcess,changed);
 }

void L1TEfficiencyMuons_Offline::dqmEndLuminosityBlock	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	c
	)

protectedvirtual

Definition at line 175 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, edm::LuminosityBlockBase::id(), and edm::LuminosityBlockID::luminosityBlock().

                                                                                                             {
     if(m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Called endLuminosityBlock at LS=" << lumiBlock.id().luminosityBlock() << endl;
 }

void L1TEfficiencyMuons_Offline::getMuonGmtPairs ( edm::Handle< l1t::MuonBxCollection > & gmtCands )

private

Definition at line 496 of file L1TEfficiencyMuons_Offline.cc.

References BXVector< T >::begin(), gather_cfg::cout, MuonGmtPair::dR(), BXVector< T >::end(), stringResolutionProvider_cfi::eta, RPCpg::mu, and EnergyCorrector::pt.

                                                                                           {
 
     m_MuonGmtPairs.clear();
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Getting muon GMT pairs" << endl;
 
     vector<const reco::Muon*>::const_iterator probeMuIt  = m_ProbeMuons.begin();
     vector<const reco::Muon*>::const_iterator probeMuEnd = m_ProbeMuons.end();
     vector<l1t::Muon> gmtContainer;
 
     for (auto mu = gmtCands->begin(0); mu != gmtCands->end(0); ++mu) {
         gmtContainer.push_back(*mu);
     }
 
     vector<l1t::Muon>::const_iterator gmtIt;
     vector<l1t::Muon>::const_iterator gmtEnd = gmtContainer.end();
 
     for (; probeMuIt!=probeMuEnd; ++probeMuIt) {
         float eta = (*probeMuIt)->eta();
         float phi = (*probeMuIt)->phi();
         float pt  = (*probeMuIt)->pt();
 
         m_EfficiencyHistos[0]["ProbeMuonEta_Histo"]->Fill(eta);
         m_EfficiencyHistos[0]["ProbeMuonPhi_Histo"]->Fill(phi);
         m_EfficiencyHistos[0]["ProbeMuonPt_Histo"]->Fill(pt);
 
         MuonGmtPair pairBestCand((*probeMuIt),0);
 //      pairBestCand.propagate(m_BField,m_propagatorAlong,m_propagatorOpposite);
         gmtIt = gmtContainer.begin();
 
         for(; gmtIt!=gmtEnd; ++gmtIt) {
             MuonGmtPair pairTmpCand((*probeMuIt),&(*gmtIt));
 //          pairTmpCand.propagate(m_BField,m_propagatorAlong,m_propagatorOpposite);     
 
             if ( (pairTmpCand.dR() < m_MaxGmtMuonDR) && (pairTmpCand.dR() < pairBestCand.dR() ) ) {
                 pairBestCand = pairTmpCand;
             }
 
         }
         m_MuonGmtPairs.push_back(pairBestCand);
         m_ControlHistos["MuonGmtDeltaR"]->Fill(pairBestCand.dR());
     }
 }

const reco::Vertex L1TEfficiencyMuons_Offline::getPrimaryVertex	(	edm::Handle< reco::VertexCollection > &	vertex,
		edm::Handle< reco::BeamSpot > &	beamSpot
	)

private

Definition at line 410 of file L1TEfficiencyMuons_Offline.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), edm::HandleBase::isValid(), reco::BeamSpot::position(), impactParameterTagInfos_cfi::primaryVertex, and reco::BeamSpot::sigmaZ().

                                                                {
     Vertex::Point posVtx;
     Vertex::Error errVtx;
 
     bool hasPrimaryVertex = false;
 
     if (vertex.isValid()) {
         vector<Vertex>::const_iterator vertexIt  = vertex->begin();
         vector<Vertex>::const_iterator vertexEnd = vertex->end();
 
         for (;vertexIt!=vertexEnd;++vertexIt) {
             if (vertexIt->isValid() && !vertexIt->isFake()) {
                 posVtx = vertexIt->position();
                 errVtx = vertexIt->error();
                 hasPrimaryVertex = true;
                 break;
             }
         }
     }
 
     if ( !hasPrimaryVertex ) {
         posVtx = beamSpot->position();
         errVtx(0,0) = beamSpot->BeamWidthX();
         errVtx(1,1) = beamSpot->BeamWidthY();
         errVtx(2,2) = beamSpot->sigmaZ();
     }
     const Vertex primaryVertex(posVtx,errVtx);
     return primaryVertex;
 }

void L1TEfficiencyMuons_Offline::getProbeMuons	(	edm::Handle< edm::TriggerResults > &	trigResults,
		edm::Handle< trigger::TriggerEvent > &	trigEvent
	)

private

Definition at line 458 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, stringResolutionProvider_cfi::eta, EnergyCorrector::pt, and mathSSE::sqrt().

                                                                        {
 
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] getting probe muons" << endl;
     m_ProbeMuons.clear();
 
     vector<const reco::Muon*>::const_iterator probeCandIt   = m_TightMuons.begin();
 //    vector<const reco::Muon*>::const_iterator probeMuIt  = m_ProbeMuons.begin();
     vector<const reco::Muon*>::const_iterator tightMuonsEnd = m_TightMuons.end();
 
     for (; probeCandIt!=tightMuonsEnd; ++probeCandIt) {
         bool tagHasTrig = false;
         vector<const reco::Muon*>::const_iterator tagCandIt  = m_TightMuons.begin();
         float deltar = 0.;
 
         for (; tagCandIt!=tightMuonsEnd; ++tagCandIt) {
             float eta = (*tagCandIt)->eta();
             float phi = (*tagCandIt)->phi();
             float pt  = (*tagCandIt)->pt();
             float dEta = eta - (*probeCandIt)->eta();
             float dPhi = phi - (*probeCandIt)->phi();
 
             deltar = sqrt(dEta*dEta + dPhi*dPhi);
 
             if ( (*tagCandIt) == (*probeCandIt) || (deltar<0.7) ) continue; // CB has a little bias for closed-by muons     
                 tagHasTrig |= matchHlt(trigEvent,(*tagCandIt));
             if (tagHasTrig) {
                 m_EfficiencyHistos[0]["TagMuonEta_Histo"]->Fill(eta);
                 m_EfficiencyHistos[0]["TagMuonPhi_Histo"]->Fill(phi);
                 m_EfficiencyHistos[0]["TagMuonPt_Histo"]->Fill(pt);
              }
         }
         if (tagHasTrig) m_ProbeMuons.push_back((*probeCandIt));
     }
     m_ControlHistos["NProbesVsTight"]->Fill(m_TightMuons.size(),m_ProbeMuons.size());
 }

void L1TEfficiencyMuons_Offline::getTightMuons	(	edm::Handle< reco::MuonCollection > &	muons,
		const reco::Vertex &	vertex
	)

private

Definition at line 442 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, and muon::isTightMuon().

                                                                                                             {
 
     if (m_verbose) cout << "[L1TEfficiencyMuons_Offline:] Getting tight muons" << endl;
     m_TightMuons.clear();
     MuonCollection::const_iterator muonIt  = muons->begin();
     MuonCollection::const_iterator muonEnd = muons->end();
 
     for(; muonIt!=muonEnd; ++muonIt) {
         if (muon::isTightMuon((*muonIt), vertex)) {
             m_TightMuons.push_back(&(*muonIt));
         }
     }
     m_ControlHistos["NTightVsAll"]->Fill(muons->size(),m_TightMuons.size());
 }

bool L1TEfficiencyMuons_Offline::matchHlt	(	edm::Handle< trigger::TriggerEvent > &	triggerEvent,
		const reco::Muon *	mu
	)

private

Definition at line 540 of file L1TEfficiencyMuons_Offline.cc.

References DEFINE_FWK_MODULE, reco::deltaR(), trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), and TriggerAnalyzer::trigObjs.

                                                                                                       {
 
     double matchDeltaR = 9999;
 
     TriggerObjectCollection trigObjs = triggerEvent->getObjects();
 
     vector<int>::const_iterator trigIndexIt  = m_trigIndices.begin();
     vector<int>::const_iterator trigIndexEnd = m_trigIndices.end();
 
     for(; trigIndexIt!=trigIndexEnd; ++trigIndexIt) {
         const vector<string> moduleLabels(m_hltConfig.moduleLabels(*trigIndexIt));
         const unsigned moduleIndex = m_hltConfig.size((*trigIndexIt))-2;
         const unsigned hltFilterIndex = triggerEvent->filterIndex(InputTag(moduleLabels[moduleIndex],"",m_trigProcess));
 
         if (hltFilterIndex < triggerEvent->sizeFilters()) {
             const Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
             const Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
             const unsigned nTriggers = triggerVids.size();
             for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
                 const TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
                 double dRtmp = deltaR((*mu),trigObject);
                 if (dRtmp < matchDeltaR) matchDeltaR = dRtmp;
             }
         }
     }
     return (matchDeltaR < m_MaxHltMuonDR);
 }