#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

	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

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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)

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)

void	getMuonGmtPairs (edm::Handle< L1MuGMTReadoutCollection > &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< BeamSpot >	m_BsInputTag

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

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

edm::EDGetTokenT < L1MuGMTReadoutCollection >	m_GmtInputTag

std::vector< int >	m_GmtPtCuts

HLTConfigProvider	m_hltConfig

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< 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 < 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 126 of file L1TEfficiencyMuons_Offline.h.

Constructor & Destructor Documentation

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

Definition at line 133 of file L1TEfficiencyMuons_Offline.cc.

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

                                                                              {
 
   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<L1MuGMTReadoutCollection>(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<string>("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 164 of file L1TEfficiencyMuons_Offline.cc.

164 { }

Member Function Documentation

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

privatevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 238 of file L1TEfficiencyMuons_Offline.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, gather_cfg::cout, eta, edm::EventSetup::get(), edm::Event::getByToken(), patZpeak::muons, phi(), HLT_25ns14e33_v1_cff::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<L1MuGMTReadoutCollection> 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("SmartPropagatorAny",m_propagatorAlong);
   eventSetup.get<TrackingComponentsRecord>().get("SmartPropagatorAnyOpposite",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);
 
   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();
 
     // unmatched gmt cands have gmtPt = -1.
     float gmtPt  = muonGmtPairsIt->gmtPt();
 
     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);
     if (gmtAboveCut) m_EfficiencyHistos[gmtPtCut]["EffvsPt" + ptTag + "Num"]->Fill(pt);
 
     if (pt > gmtPtCut + 8.) { // efficiency in eta/phi at plateau
       
       m_EfficiencyHistos[gmtPtCut]["EffvsPhi" + ptTag + "Den"]->Fill(phi);
       m_EfficiencyHistos[gmtPtCut]["EffvsEta" + ptTag + "Den"]->Fill(eta);
       
       if (gmtAboveCut) { 
         m_EfficiencyHistos[gmtPtCut]["EffvsPhi" + ptTag + "Num"]->Fill(phi);
         m_EfficiencyHistos[gmtPtCut]["EffvsEta" + ptTag + "Num"]->Fill(eta);
       }
     }
       }
     }
   }
 
   cout << "[L1TEfficiencyMuons_Offline:] Computation finished" << endl;
 
   
 }

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

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 216 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 320 of file L1TEfficiencyMuons_Offline.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, mergeVDriftHistosByStation::name, 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);
   
 }

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

protectedvirtual

Definition at line 339 of file L1TEfficiencyMuons_Offline.cc.

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

                                                                                          { 
   
   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/");
 
   string effTag[2] = {"Den", "Num"};
   
   for(int iEffTag=0; iEffTag<2; ++ iEffTag) {
     string name = "EffvsPt" + ptTag + effTag[iEffTag];
     m_EfficiencyHistos[ptCut][name] = ibooker.book1D(name.c_str(),name.c_str(),16,0.,40.);
     
     name = "EffvsPhi" + ptTag + effTag[iEffTag];
     m_EfficiencyHistos[ptCut][name] = ibooker.book1D(name.c_str(),name.c_str(),12,0.,2*TMath::Pi());
     
     name = "EffvsEta" + ptTag + effTag[iEffTag];
     m_EfficiencyHistos[ptCut][name] = ibooker.book1D(name.c_str(),name.c_str(),12,-2.4,2.4);
   }
   
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 176 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout.

                                                                                                                        {
 
   //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
	)

protectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 166 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 227 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< L1MuGMTReadoutCollection > & gmtCands )

private

Definition at line 466 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, MuonGmtPair::dR(), MuonGmtPair::gmtPt(), and MuonGmtPair::propagate().

                                                                                                {
 
   m_MuonGmtPairs.clear();
   
   cout << "[L1TEfficiencyMuons_Offline:] Getting muon GMT pairs" << endl;  
 
   vector<const Muon*>::const_iterator probeMuIt  = m_ProbeMuons.begin();
   vector<const Muon*>::const_iterator probeMuEnd = m_ProbeMuons.end();
 
   vector<L1MuGMTExtendedCand> gmtContainer = gmtCands->getRecord(0).getGMTCands();
   
   vector<L1MuGMTExtendedCand>::const_iterator gmtIt;
   vector<L1MuGMTExtendedCand>::const_iterator gmtEnd = gmtContainer.end();
   
   for (; probeMuIt!=probeMuEnd; ++probeMuIt) {
     
     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.gmtPt() > pairBestCand.gmtPt())
     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 368 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout, edm::HandleBase::isValid(), and HLT_25ns14e33_v1_cff::primaryVertex.

                                                                {
   
   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 ) {
 
     if(m_verbose){
       cout << "[L1TEfficiencyMuons_Offline:] PrimaryVertex not found, use BeamSpot position instead" << endl;
     }
     
     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 437 of file L1TEfficiencyMuons_Offline.cc.

References gather_cfg::cout.

                                                                        {
 
   cout << "[L1TEfficiencyMuons_Offline:] getting probe muons" << endl;  
 
   m_ProbeMuons.clear();
   
   vector<const Muon*>::const_iterator probeCandIt   = m_TightMuons.begin();
   vector<const Muon*>::const_iterator tightMuonsEnd = m_TightMuons.end();
 
   for (; probeCandIt!=tightMuonsEnd; ++probeCandIt) {
     
     bool tagHasTrig = false;
     vector<const Muon*>::const_iterator tagCandIt  = m_TightMuons.begin();
     
     for (; tagCandIt!=tightMuonsEnd; ++tagCandIt) {
       if ((*tagCandIt) == (*probeCandIt)) continue; // CB has a little bias for closed-by muons
       tagHasTrig |= matchHlt(trigEvent,(*tagCandIt));
     }
     
     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 416 of file L1TEfficiencyMuons_Offline.cc.

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

                                                                                                             {
 
   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 506 of file L1TEfficiencyMuons_Offline.cc.

References reco::deltaR(), HLT_25ns14e33_v1_cff::InputTag, 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);
 
 }