#include <L1TMuonDQMOffline.h>

Inheritance diagram for L1TMuonDQMOffline:

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

	~L1TMuonDQMOffline () override

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

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)

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

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

	~EDAnalyzerBase () 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

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

virtual	~EDConsumerBase () noexcept(false)

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

virtual void	bookControlHistos (DQMStore::IBooker &)

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

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

virtual void	bookResolutionHistos (DQMStore::IBooker &ibooker)

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< Control, MonitorElement * >	m_ControlHistos

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

std::vector< double >	m_effVsEtaBins

std::vector< double >	m_effVsPhiBins

std::vector< double >	m_effVsPtBins

edm::EDGetTokenT< l1t::MuonBxCollection >	m_GmtInputTag

std::vector< int >	m_GmtPtCuts

std::string	m_HistFolder

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::map< Resol, MonitorElement * >	m_ResolutionHistos

double	m_TagPtCut

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 115 of file L1TMuonDQMOffline.h.

Member Enumeration Documentation

enum L1TMuonDQMOffline::Control

private

Enumerator
CONTROL_MuonGmtDeltaR
CONTROL_NTightVsAll
CONTROL_NProbesVsTight

Definition at line 177 of file L1TMuonDQMOffline.h.

         {
             CONTROL_MuonGmtDeltaR,
             CONTROL_NTightVsAll,
             CONTROL_NProbesVsTight,
         };

enum L1TMuonDQMOffline::Resol

private

Enumerator
RESOL_Pt
RESOL_1overPt
RESOL_Eta
RESOL_Phi
RESOL_Charge
RESOL_Pt_OPEN
RESOL_1overPt_OPEN
RESOL_Eta_OPEN
RESOL_Phi_OPEN
RESOL_Charge_OPEN
RESOL_Pt_DOUBLE
RESOL_1overPt_DOUBLE
RESOL_Eta_DOUBLE
RESOL_Phi_DOUBLE
RESOL_Charge_DOUBLE
RESOL_Pt_SINGLE
RESOL_1overPt_SINGLE
RESOL_Eta_SINGLE
RESOL_Phi_SINGLE
RESOL_Charge_SINGLE

Definition at line 150 of file L1TMuonDQMOffline.h.

         {
             RESOL_Pt,
             RESOL_1overPt,
             RESOL_Eta,
             RESOL_Phi,
             RESOL_Charge,
 
             RESOL_Pt_OPEN,
             RESOL_1overPt_OPEN,
             RESOL_Eta_OPEN,
             RESOL_Phi_OPEN,
             RESOL_Charge_OPEN,
 
             RESOL_Pt_DOUBLE,
             RESOL_1overPt_DOUBLE,
             RESOL_Eta_DOUBLE,
             RESOL_Phi_DOUBLE,
             RESOL_Charge_DOUBLE,
 
             RESOL_Pt_SINGLE,
             RESOL_1overPt_SINGLE,
             RESOL_Eta_SINGLE,
             RESOL_Phi_SINGLE,
             RESOL_Charge_SINGLE,
         };

Constructor & Destructor Documentation

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

Definition at line 106 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, m_MaxGmtMuonDR, m_MaxHltMuonDR, m_MaxMuonEta, and m_verbose.

                                                             :
     m_verbose(ps.getUntrackedParameter<bool>("verbose")),
     m_HistFolder(ps.getUntrackedParameter<string>("histFolder")),
     m_GmtPtCuts(ps.getUntrackedParameter< vector<int> >("gmtPtCuts")),
     m_TagPtCut(ps.getUntrackedParameter<double>("tagPtCut")),
     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")),
     m_effVsPtBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPtBins")),
     m_effVsPhiBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPhiBins")),
     m_effVsEtaBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsEtaBins"))
 {
     if (m_verbose) cout << "[L1TMuonDQMOffline:] ____________ Storage initialization ____________ " << endl;
     // 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;
 }

L1TMuonDQMOffline::~L1TMuonDQMOffline ( )

override

Definition at line 133 of file L1TMuonDQMOffline.cc.

133 { }

Member Function Documentation

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

overrideprivate

Definition at line 183 of file L1TMuonDQMOffline.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, BXVector< T >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, gather_cfg::cout, BXVector< T >::end(), PVValHelper::eta, edm::EventSetup::get(), edm::Event::getByToken(), getMuonGmtPairs(), getPrimaryVertex(), getProbeMuons(), getTightMuons(), m_BField, m_BsInputTag, m_EfficiencyHistos, m_GmtInputTag, m_GmtPtCuts, m_MaxMuonEta, m_MuonGmtPairs, m_MuonInputTag, m_propagatorAlong, m_propagatorOpposite, m_ResolutionHistos, m_trigInputTag, m_trigProcess_token, m_verbose, m_VtxInputTag, RPCpg::mu, nano_cff::muons, phi, impactParameterTagInfos_cfi::primaryVertex, EnergyCorrector::pt, RESOL_1overPt, RESOL_1overPt_DOUBLE, RESOL_1overPt_OPEN, RESOL_1overPt_SINGLE, RESOL_Charge, RESOL_Charge_DOUBLE, RESOL_Charge_OPEN, RESOL_Charge_SINGLE, RESOL_Eta, RESOL_Eta_DOUBLE, RESOL_Eta_OPEN, RESOL_Eta_SINGLE, RESOL_Phi, RESOL_Phi_DOUBLE, RESOL_Phi_OPEN, RESOL_Phi_SINGLE, RESOL_Pt, RESOL_Pt_DOUBLE, RESOL_Pt_OPEN, and RESOL_Pt_SINGLE.

                                                                                   {
 
     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 << "[L1TMuonDQMOffline:] Computing efficiencies" << endl;
 
     vector<MuonGmtPair>::const_iterator muonGmtPairsIt  = m_MuonGmtPairs.begin();
     vector<MuonGmtPair>::const_iterator muonGmtPairsEnd = m_MuonGmtPairs.end();
 
     for(; muonGmtPairsIt!=muonGmtPairsEnd; ++muonGmtPairsIt) {
         float pt  = muonGmtPairsIt->pt();
         float eta = muonGmtPairsIt->eta();
         float phi = muonGmtPairsIt->phi();
         int charge = muonGmtPairsIt->charge();
 
         float gmtPt  = muonGmtPairsIt->gmtPt();
         float gmtEta  = muonGmtPairsIt->gmtEta();
         float gmtPhi  = muonGmtPairsIt->gmtPhi();
         int gmtCharge  = muonGmtPairsIt->gmtCharge();
         int qual = muonGmtPairsIt->gmtQual();
 
         vector<int>::const_iterator gmtPtCutsIt  = m_GmtPtCuts.begin();
         vector<int>::const_iterator gmtPtCutsEnd = m_GmtPtCuts.end();
 
         if ( (fabs(eta) < m_MaxMuonEta) && (gmtPt > 0) ) {
             m_ResolutionHistos[RESOL_Pt]->Fill(pt - gmtPt);
             m_ResolutionHistos[RESOL_1overPt]->Fill(1/pt - 1/gmtPt);
             m_ResolutionHistos[RESOL_Eta]->Fill(eta - gmtEta);
             m_ResolutionHistos[RESOL_Phi]->Fill(phi - gmtPhi);
             m_ResolutionHistos[RESOL_Charge]->Fill(charge - gmtCharge);
             if (qual >= 4) {
                 m_ResolutionHistos[RESOL_Pt_OPEN]->Fill(pt - gmtPt);
                 m_ResolutionHistos[RESOL_1overPt_OPEN]->Fill(1/pt - 1/gmtPt);
                 m_ResolutionHistos[RESOL_Eta_OPEN]->Fill(eta - gmtEta);
                 m_ResolutionHistos[ RESOL_Phi_OPEN]->Fill(phi - gmtPhi);
                 m_ResolutionHistos[RESOL_Charge_OPEN]->Fill(charge - gmtCharge);
                 if (qual >= 8) {
                     m_ResolutionHistos[RESOL_Pt_DOUBLE]->Fill(pt - gmtPt);
                     m_ResolutionHistos[RESOL_1overPt_DOUBLE]->Fill(1/pt - 1/gmtPt);
                     m_ResolutionHistos[RESOL_Eta_DOUBLE]->Fill(eta - gmtEta);
                     m_ResolutionHistos[RESOL_Phi_DOUBLE]->Fill(phi - gmtPhi);
                     m_ResolutionHistos[RESOL_Charge_DOUBLE]->Fill(charge - gmtCharge);
                     if (qual >= 12) {
                         m_ResolutionHistos[RESOL_Pt_SINGLE]->Fill(pt - gmtPt);
                         m_ResolutionHistos[RESOL_1overPt_SINGLE]->Fill(1/pt - 1/gmtPt);
                         m_ResolutionHistos[RESOL_Eta_SINGLE]->Fill(eta - gmtEta);
                         m_ResolutionHistos[RESOL_Phi_SINGLE]->Fill(phi - gmtPhi);
                         m_ResolutionHistos[RESOL_Charge_SINGLE]->Fill(charge - gmtCharge);
                     }
                 }
             }
         }
 
         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 << "[L1TMuonDQMOffline:] Computation finished" << endl;
 }

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

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 173 of file L1TMuonDQMOffline.cc.

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

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

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

protectedvirtual

Definition at line 335 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), CONTROL_MuonGmtDeltaR, CONTROL_NProbesVsTight, CONTROL_NTightVsAll, gather_cfg::cout, m_ControlHistos, m_EfficiencyHistos, m_HistFolder, m_verbose, dataset::name, Pi, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

                                                                   {
     if(m_verbose) cout << "[L1TMuonDQMOffline:] Booking Control Plot Histos" << endl;
 
     ibooker.setCurrentFolder(m_HistFolder+"/control_variables");
 
     string name = "MuonGmtDeltaR";
     m_ControlHistos[CONTROL_MuonGmtDeltaR] = ibooker.book1D(name.c_str(),name.c_str(),25.,0.,2.5);
 
     name = "NTightVsAll";
     m_ControlHistos[CONTROL_NTightVsAll] = ibooker.book2D(name.c_str(),name.c_str(),16,-0.5,15.5,16,-0.5,15.5);
 
     name = "NProbesVsTight";
     m_ControlHistos[CONTROL_NProbesVsTight] = ibooker.book2D(name.c_str(),name.c_str(),8,-0.5,7.5,8,-0.5,7.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 L1TMuonDQMOffline::bookEfficiencyHistos	(	DQMStore::IBooker &	ibooker,
		int	ptCut
	)

protectedvirtual

Definition at line 366 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), gather_cfg::cout, m_EfficiencyHistos, m_effVsEtaBins, m_effVsPhiBins, m_effVsPtBins, m_HistFolder, m_verbose, HiGenCleaner_cff::ptCut, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

                                                                                 {
     if(m_verbose) cout << "[L1TMuonDQMOffline:] Booking Efficiency Plot Histos for pt cut = " << ptCut << endl;
 
     stringstream ptCutToTag; ptCutToTag << ptCut;
     string ptTag = ptCutToTag.str();
 
     ibooker.setCurrentFolder(m_HistFolder+"/numerators_and_denominators");
 
     std::vector<float> effVsPtBins(m_effVsPtBins.begin(), m_effVsPtBins.end());
     int nEffVsPtBins = effVsPtBins.size() - 1;
     float* ptBinsArray = &(effVsPtBins[0]);
 
     string name1 = "EffvsPt_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
     string name2 = "EffvsPt_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
 
     name1 = "EffvsPt_OPEN_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
     name2 = "EffvsPt_OPEN_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
 
     name1 = "EffvsPt_DOUBLE_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
     name2 = "EffvsPt_DOUBLE_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
 
     name1 = "EffvsPt_SINGLE_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
     name2 = "EffvsPt_SINGLE_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
 
 
     std::vector<float> effVsPhiBins(m_effVsPhiBins.begin(), m_effVsPhiBins.end());
     int nEffVsPhiBins = effVsPhiBins.size() - 1;
     float* phiBinsArray = &(effVsPhiBins[0]);
 
     std::vector<float> effVsEtaBins(m_effVsEtaBins.begin(), m_effVsEtaBins.end());
     int nEffVsEtaBins = effVsEtaBins.size() - 1;
     float* etaBinsArray = &(effVsEtaBins[0]);
 
     name1 = "EffvsPhi_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
     name2 = "EffvsPhi_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
 
     name1 = "EffvsEta_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
     name2 = "EffvsEta_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
 
 
     name1 = "EffvsPhi_OPEN_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
     name2 = "EffvsPhi_OPEN_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
 
     name1 = "EffvsEta_OPEN_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
     name2 = "EffvsEta_OPEN_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
 
 
     name1 = "EffvsPhi_DOUBLE_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
     name2 = "EffvsPhi_DOUBLE_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
 
     name1 = "EffvsEta_DOUBLE_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
     name2 = "EffvsEta_DOUBLE_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
 
 
     name1 = "EffvsPhi_SINGLE_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
     name2 = "EffvsPhi_SINGLE_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
 
     name1 = "EffvsEta_SINGLE_" + ptTag + "_Den";
     m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
     name2 = "EffvsEta_SINGLE_" + ptTag + "_Num";
     m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
 }

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

overrideprotected

Definition at line 142 of file L1TMuonDQMOffline.cc.

References bookControlHistos(), bookEfficiencyHistos(), bookResolutionHistos(), gather_cfg::cout, createfilelist::int, m_GmtPtCuts, m_hltConfig, m_trigIndices, m_trigNames, m_verbose, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().

                                                                                                               {
     //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));
     }
     bookResolutionHistos(ibooker);
 
     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 << "[L1TMuonDQMOffline:] Warning: Could not find trigger " << (*trigNamesIt) << endl;
     }
 }

void L1TMuonDQMOffline::bookResolutionHistos ( DQMStore::IBooker & ibooker )

protectedvirtual

Definition at line 455 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), gather_cfg::cout, m_HistFolder, m_ResolutionHistos, m_verbose, dataset::name, RESOL_1overPt, RESOL_1overPt_DOUBLE, RESOL_1overPt_OPEN, RESOL_1overPt_SINGLE, RESOL_Charge, RESOL_Charge_DOUBLE, RESOL_Charge_OPEN, RESOL_Charge_SINGLE, RESOL_Eta, RESOL_Eta_DOUBLE, RESOL_Eta_OPEN, RESOL_Eta_SINGLE, RESOL_Phi, RESOL_Phi_DOUBLE, RESOL_Phi_OPEN, RESOL_Phi_SINGLE, RESOL_Pt, RESOL_Pt_DOUBLE, RESOL_Pt_OPEN, RESOL_Pt_SINGLE, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

                                                                      {
     if(m_verbose) cout << "[L1TMuonOffline:] Booking Resolution Plot Histos" << endl;
 
     ibooker.setCurrentFolder(m_HistFolder+"/resolution");
 
     string name = "Pt_Res";
     m_ResolutionHistos[RESOL_Pt] = ibooker.book1D(name.c_str(),name.c_str(),100,-50.,50.);
     name = "Pt_Res_OPEN";
     m_ResolutionHistos[RESOL_Pt_OPEN] = ibooker.book1D(name.c_str(),name.c_str(),100,-50.,50.);
     name = "Pt_Res_DOUBLE";
     m_ResolutionHistos[RESOL_Pt_DOUBLE] = ibooker.book1D(name.c_str(),name.c_str(),100,-50.,50.);
     name = "Pt_Res_SINGLE";
     m_ResolutionHistos[RESOL_Pt_SINGLE] = ibooker.book1D(name.c_str(),name.c_str(),100,-50.,50.);
     name = "1overPt_Res";
     m_ResolutionHistos[RESOL_1overPt] = ibooker.book1D(name.c_str(),name.c_str(),50,-0.05,0.05);
     name = "1overPt_Res_OPEN";
     m_ResolutionHistos[RESOL_1overPt_OPEN] = ibooker.book1D(name.c_str(),name.c_str(),50,-0.05,0.05);
     name = "1overPt_Res_DOUBLE";
     m_ResolutionHistos[RESOL_1overPt_DOUBLE] = ibooker.book1D(name.c_str(),name.c_str(),50,-0.05,0.05);
     name = "1overPt_Res_SINGLE";
     m_ResolutionHistos[RESOL_1overPt_SINGLE] = ibooker.book1D(name.c_str(),name.c_str(),50,-0.05,0.05);
     name = "Eta_Res";
     m_ResolutionHistos[RESOL_Eta] = ibooker.book1D(name.c_str(),name.c_str(),100,-0.1,0.1);
     name = "Eta_Res_OPEN";
     m_ResolutionHistos[RESOL_Eta_OPEN] = ibooker.book1D(name.c_str(),name.c_str(),100,-0.1,0.1);
     name = "Eta_Res_DOUBLE";
     m_ResolutionHistos[RESOL_Eta_DOUBLE] = ibooker.book1D(name.c_str(),name.c_str(),100,-0.1,0.1);
     name = "Eta_Res_SINGLE";
     m_ResolutionHistos[RESOL_Eta_SINGLE] = ibooker.book1D(name.c_str(),name.c_str(),100,-0.1,0.1);
     name = "Phi_Res";
     m_ResolutionHistos[RESOL_Phi] = ibooker.book1D(name.c_str(),name.c_str(),96,-0.2,0.2);
     name = "Phi_Res_OPEN";
     m_ResolutionHistos[RESOL_Phi_OPEN] = ibooker.book1D(name.c_str(),name.c_str(),96,-0.2,0.2);
     name = "Phi_Res_DOUBLE";
     m_ResolutionHistos[RESOL_Phi_DOUBLE] = ibooker.book1D(name.c_str(),name.c_str(),96,-0.2,0.2);
     name = "Phi_Res_SINGLE";
     m_ResolutionHistos[RESOL_Phi_SINGLE] = ibooker.book1D(name.c_str(),name.c_str(),96,-0.2,0.2);
     name = "Charge_Res";
     m_ResolutionHistos[RESOL_Charge] = ibooker.book1D(name.c_str(),name.c_str(),5,-2,3);
     name = "Charge_Res_OPEN";
     m_ResolutionHistos[RESOL_Charge_OPEN] = ibooker.book1D(name.c_str(),name.c_str(),5,-2,3);
     name = "Charge_Res_DOUBLE";
     m_ResolutionHistos[RESOL_Charge_DOUBLE] = ibooker.book1D(name.c_str(),name.c_str(),5,-2,3);
     name = "Charge_Res_SINGLE";
     m_ResolutionHistos[RESOL_Charge_SINGLE] = ibooker.book1D(name.c_str(),name.c_str(),5,-2,3);
 }

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

overrideprotected

Definition at line 135 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, HLTConfigProvider::init(), m_hltConfig, m_trigProcess, and m_verbose.

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

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

protectedvirtual

Definition at line 178 of file L1TMuonDQMOffline.cc.

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

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

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

private

Definition at line 610 of file L1TMuonDQMOffline.cc.

References BXVector< T >::begin(), CONTROL_MuonGmtDeltaR, gather_cfg::cout, MuonGmtPair::dR(), BXVector< T >::end(), PVValHelper::eta, m_ControlHistos, m_EfficiencyHistos, m_MaxGmtMuonDR, m_MuonGmtPairs, m_ProbeMuons, m_verbose, RPCpg::mu, phi, and EnergyCorrector::pt.

Referenced by analyze().

                                                                                  {
 
     m_MuonGmtPairs.clear();
     if (m_verbose) cout << "[L1TMuonDQMOffline:] 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),nullptr);
 //      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[CONTROL_MuonGmtDeltaR]->Fill(pairBestCand.dR());
     }
 }

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

private

Definition at line 507 of file L1TMuonDQMOffline.cc.

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

Referenced by analyze().

                                                                {
     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 L1TMuonDQMOffline::getProbeMuons	(	edm::Handle< edm::TriggerResults > &	trigResults,
		edm::Handle< trigger::TriggerEvent > &	trigEvent
	)

private

Definition at line 555 of file L1TMuonDQMOffline.cc.

References CONTROL_NProbesVsTight, gather_cfg::cout, SoftLeptonByDistance_cfi::distance, PVValHelper::eta, m_ControlHistos, m_EfficiencyHistos, m_ProbeMuons, m_TagPtCut, m_TightMuons, m_verbose, matchHlt(), phi, EnergyCorrector::pt, and mathSSE::sqrt().

Referenced by analyze().

                                                                        {
 
     if (m_verbose) cout << "[L1TMuonDQMOffline:] getting probe muons" << endl;
     m_ProbeMuons.clear();
     std::vector<const reco::Muon*>  tagMuonsInHist;
 
     tagMuonsInHist.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 isProbe = false;
         vector<const reco::Muon*>::const_iterator tagCandIt  = m_TightMuons.begin();
         float deltar = 0.;
 
         for (; tagCandIt!=tightMuonsEnd; ++tagCandIt) {
             bool tagMuonAlreadyInHist = false;
             bool tagHasTrig = false;
             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)) && (pt > m_TagPtCut);
             isProbe |= tagHasTrig;
             if (tagHasTrig) {
                 if (std::distance(m_TightMuons.begin(), m_TightMuons.end()) > 2 ) {
                     for (vector<const reco::Muon*>::const_iterator tagMuonsInHistIt = tagMuonsInHist.begin(); tagMuonsInHistIt!=tagMuonsInHist.end(); ++tagMuonsInHistIt) {
                         if ( (*tagCandIt) == (*tagMuonsInHistIt) )  {
                             tagMuonAlreadyInHist = true;
                            break;
                         }
                     }
                     if (tagMuonAlreadyInHist == false) tagMuonsInHist.push_back((*tagCandIt));
                 }
                 if (tagMuonAlreadyInHist == false) {
                     m_EfficiencyHistos[0]["TagMuonEta_Histo"]->Fill(eta);
                     m_EfficiencyHistos[0]["TagMuonPhi_Histo"]->Fill(phi);
                     m_EfficiencyHistos[0]["TagMuonPt_Histo"]->Fill(pt);
                 }
             }
         }
         if (isProbe) m_ProbeMuons.push_back((*probeCandIt));
     }
     m_ControlHistos[CONTROL_NProbesVsTight]->Fill(m_TightMuons.size(),m_ProbeMuons.size());
 }

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

private

Definition at line 539 of file L1TMuonDQMOffline.cc.

References CONTROL_NTightVsAll, gather_cfg::cout, muon::isTightMuon(), m_ControlHistos, m_TightMuons, and m_verbose.

Referenced by analyze().

                                                                                                    {
 
     if (m_verbose) cout << "[L1TMuonDQMOffline:] 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[CONTROL_NTightVsAll]->Fill(muons->size(),m_TightMuons.size());
 }

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

private

Definition at line 654 of file L1TMuonDQMOffline.cc.

References DEFINE_FWK_MODULE, deltaR(), trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), m_hltConfig, m_MaxHltMuonDR, m_trigIndices, m_trigProcess, HLTConfigProvider::moduleLabels(), HLTConfigProvider::size(), and TriggerAnalyzer::trigObjs.

Referenced by getProbeMuons().

                                                                                              {
 
     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);
 }

Member Data Documentation

edm::ESHandle<MagneticField> L1TMuonDQMOffline::m_BField

private

Definition at line 146 of file L1TMuonDQMOffline.h.

Referenced by analyze().

edm::EDGetTokenT<reco::BeamSpot> L1TMuonDQMOffline::m_BsInputTag

private

Definition at line 207 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::map<Control, MonitorElement*> L1TMuonDQMOffline::m_ControlHistos

private

Definition at line 187 of file L1TMuonDQMOffline.h.

Referenced by bookControlHistos(), getMuonGmtPairs(), getProbeMuons(), and getTightMuons().

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

private

Definition at line 185 of file L1TMuonDQMOffline.h.

Referenced by analyze(), bookControlHistos(), bookEfficiencyHistos(), getMuonGmtPairs(), and getProbeMuons().

std::vector<double> L1TMuonDQMOffline::m_effVsEtaBins

private

Definition at line 214 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos().

std::vector<double> L1TMuonDQMOffline::m_effVsPhiBins

private

Definition at line 213 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos().

std::vector<double> L1TMuonDQMOffline::m_effVsPtBins

private

Definition at line 212 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos().

edm::EDGetTokenT<l1t::MuonBxCollection> L1TMuonDQMOffline::m_GmtInputTag

private

Definition at line 205 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<int> L1TMuonDQMOffline::m_GmtPtCuts

private

Definition at line 202 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookHistograms().

std::string L1TMuonDQMOffline::m_HistFolder

private

Definition at line 201 of file L1TMuonDQMOffline.h.

Referenced by bookControlHistos(), bookEfficiencyHistos(), and bookResolutionHistos().

HLTConfigProvider L1TMuonDQMOffline::m_hltConfig

private

Definition at line 144 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), dqmBeginRun(), and matchHlt().

BXVector<l1t::Muon> L1TMuonDQMOffline::m_L1tL1tMuons

private

Definition at line 197 of file L1TMuonDQMOffline.h.

std::vector<l1t::MuonBxCollection> L1TMuonDQMOffline::m_L1tMuons

private

Definition at line 195 of file L1TMuonDQMOffline.h.

float L1TMuonDQMOffline::m_MaxGmtMuonDR

private

Definition at line 219 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and L1TMuonDQMOffline().

float L1TMuonDQMOffline::m_MaxHltMuonDR

private

Definition at line 220 of file L1TMuonDQMOffline.h.

Referenced by L1TMuonDQMOffline(), and matchHlt().

float L1TMuonDQMOffline::m_MaxMuonEta

private

Definition at line 218 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::vector<MuonGmtPair> L1TMuonDQMOffline::m_MuonGmtPairs

private

Definition at line 192 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and getMuonGmtPairs().

edm::EDGetTokenT<reco::MuonCollection> L1TMuonDQMOffline::m_MuonInputTag

private

Definition at line 204 of file L1TMuonDQMOffline.h.

Referenced by analyze().

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

private

Definition at line 191 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and getProbeMuons().

edm::ESHandle<Propagator> L1TMuonDQMOffline::m_propagatorAlong

private

Definition at line 147 of file L1TMuonDQMOffline.h.

Referenced by analyze().

edm::ESHandle<Propagator> L1TMuonDQMOffline::m_propagatorOpposite

private

Definition at line 148 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<reco::MuonCollection> L1TMuonDQMOffline::m_RecoMuons

private

Definition at line 194 of file L1TMuonDQMOffline.h.

std::vector<reco::Muon> L1TMuonDQMOffline::m_RecoRecoMuons

private

Definition at line 196 of file L1TMuonDQMOffline.h.

std::map<Resol, MonitorElement*> L1TMuonDQMOffline::m_ResolutionHistos

private

Definition at line 186 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookResolutionHistos().

double L1TMuonDQMOffline::m_TagPtCut

private

Definition at line 203 of file L1TMuonDQMOffline.h.

Referenced by getProbeMuons().

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

private

Definition at line 190 of file L1TMuonDQMOffline.h.

Referenced by getProbeMuons(), and getTightMuons().

std::vector<int> L1TMuonDQMOffline::m_trigIndices

private

Definition at line 216 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), and matchHlt().

edm::EDGetTokenT<trigger::TriggerEvent> L1TMuonDQMOffline::m_trigInputTag

private

Definition at line 208 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<std::string> L1TMuonDQMOffline::m_trigNames

private

Definition at line 211 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms().

std::string L1TMuonDQMOffline::m_trigProcess

private

Definition at line 209 of file L1TMuonDQMOffline.h.

Referenced by dqmBeginRun(), and matchHlt().

edm::EDGetTokenT<edm::TriggerResults> L1TMuonDQMOffline::m_trigProcess_token

private

Definition at line 210 of file L1TMuonDQMOffline.h.

Referenced by analyze().

bool L1TMuonDQMOffline::m_verbose

private

Definition at line 200 of file L1TMuonDQMOffline.h.

Referenced by analyze(), beginLuminosityBlock(), bookControlHistos(), bookEfficiencyHistos(), bookHistograms(), bookResolutionHistos(), dqmBeginRun(), dqmEndLuminosityBlock(), getMuonGmtPairs(), getProbeMuons(), getTightMuons(), and L1TMuonDQMOffline().

edm::EDGetTokenT<reco::VertexCollection> L1TMuonDQMOffline::m_VtxInputTag

private

Definition at line 206 of file L1TMuonDQMOffline.h.

Referenced by analyze().

Private Types
enum	Control { CONTROL_MuonGmtDeltaR, CONTROL_NTightVsAll, CONTROL_NProbesVsTight }

enum	Resol { RESOL_Pt, RESOL_1overPt, RESOL_Eta, RESOL_Phi, RESOL_Charge, RESOL_Pt_OPEN, RESOL_1overPt_OPEN, RESOL_Eta_OPEN, RESOL_Phi_OPEN, RESOL_Charge_OPEN, RESOL_Pt_DOUBLE, RESOL_1overPt_DOUBLE, RESOL_Eta_DOUBLE, RESOL_Phi_DOUBLE, RESOL_Charge_DOUBLE, RESOL_Pt_SINGLE, RESOL_1overPt_SINGLE, RESOL_Eta_SINGLE, RESOL_Phi_SINGLE, RESOL_Charge_SINGLE }

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation