#include <HLTMuonL2SelectorForL3IO.h>

Inheritance diagram for HLTMuonL2SelectorForL3IO:

Public Member Functions
	HLTMuonL2SelectorForL3IO (const edm::ParameterSet &)
	constructor with config More...

void	produce (edm::Event &, const edm::EventSetup &) override
	select muons More...

	~HLTMuonL2SelectorForL3IO () override
	destructor More...

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
	default values More...

Private Attributes
const bool	applyL3Filters_

const edm::EDGetTokenT< reco::TrackCollection >	l2Src_

const edm::EDGetTokenT< reco::MuonTrackLinksCollection >	l3linkToken_

const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	l3OISrc_

const double	max_NormalizedChi2_

const double	max_PtDifference_

const int	min_Nhits_

const int	min_NmuonHits_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
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

Detailed Description

L2 muon selector for L3 IO: finds L2 muons not previous converted into (good) L3 muons

Author: Benjamin Radburn-Smith, Santiago Folgueras - Purdue University

Definition at line 27 of file HLTMuonL2SelectorForL3IO.h.

Constructor & Destructor Documentation

◆ HLTMuonL2SelectorForL3IO()

HLTMuonL2SelectorForL3IO::HLTMuonL2SelectorForL3IO ( const edm::ParameterSet & iConfig )

constructor with config

Definition at line 17 of file HLTMuonL2SelectorForL3IO.cc.

     : l2Src_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("l2Src"))),
       l3OISrc_(consumes<reco::RecoChargedCandidateCollection>(iConfig.getParameter<edm::InputTag>("l3OISrc"))),
       l3linkToken_(consumes<reco::MuonTrackLinksCollection>(iConfig.getParameter<edm::InputTag>("InputLinks"))),
       applyL3Filters_(iConfig.getParameter<bool>("applyL3Filters")),
       max_NormalizedChi2_(iConfig.getParameter<double>("MaxNormalizedChi2")),
       max_PtDifference_(iConfig.getParameter<double>("MaxPtDifference")),
       min_Nhits_(iConfig.getParameter<int>("MinNhits")),
       min_NmuonHits_(iConfig.getParameter<int>("MinNmuonHits")) {
   LogTrace("Muon|RecoMuon|HLTMuonL2SelectorForL3IO") << "constructor called";
   produces<reco::TrackCollection>();
 }

References LogTrace.

◆ ~HLTMuonL2SelectorForL3IO()

HLTMuonL2SelectorForL3IO::~HLTMuonL2SelectorForL3IO ( )

override

destructor

Definition at line 31 of file HLTMuonL2SelectorForL3IO.cc.

31 {}

Member Function Documentation

◆ fillDescriptions()

void HLTMuonL2SelectorForL3IO::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

default values

Definition at line 100 of file HLTMuonL2SelectorForL3IO.cc.

                                                                                           {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("l2Src", edm::InputTag("hltL2Muons", "UpdatedAtVtx"));
   desc.add<edm::InputTag>("l3OISrc", edm::InputTag("hltNewOIL3MuonCandidates"));
   desc.add<edm::InputTag>("InputLinks", edm::InputTag("hltNewOIL3MuonsLinksCombination"));
   desc.add<bool>("applyL3Filters", true);
   desc.add<int>("MinNhits", 1);
   desc.add<double>("MaxNormalizedChi2", 20.0);
   desc.add<int>("MinNmuonHits", 0);
   desc.add<double>("MaxPtDifference", 999.0);  //relative difference
   descriptions.add("HLTMuonL2SelectorForL3IO", desc);
 }

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

◆ produce()

void HLTMuonL2SelectorForL3IO::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

select muons

create collection of L2 muons not already reconstructed as L3 muons

Definition at line 34 of file HLTMuonL2SelectorForL3IO.cc.

                                                                                     {
   const std::string metname = "Muon|RecoMuon|HLTMuonL2SelectorForL3IO";
  
   //  IN:
   edm::Handle<reco::TrackCollection> l2muonH;
   iEvent.getByToken(l2Src_, l2muonH);
  
   edm::Handle<reco::RecoChargedCandidateCollection> l3muonH;
   iEvent.getByToken(l3OISrc_, l3muonH);
  
   // Read Links collection:
   edm::Handle<reco::MuonTrackLinksCollection> links;
   iEvent.getByToken(l3linkToken_, links);
  
   //    OUT:
   std::unique_ptr<reco::TrackCollection> result(new reco::TrackCollection());
  
   for (unsigned int il2 = 0; il2 != l2muonH->size(); ++il2) {
     reco::TrackRef l2muRef(l2muonH, il2);
     bool re_do_this_L2 = true;
  
     for (unsigned int il3 = 0; il3 != l3muonH->size(); ++il3) {
       reco::RecoChargedCandidateRef cand(l3muonH, il3);
       reco::TrackRef tk = cand->track();
  
       bool useThisLink = false;
       for (unsigned int l(0); l < links->size() && !useThisLink; ++l) {
         const reco::MuonTrackLinks* link = &links->at(l);
  
         // Check if the L3 link matches the L3 candidate
         const reco::Track& globalTrack = *link->globalTrack();
         float dR2 = deltaR2(tk->eta(), tk->phi(), globalTrack.eta(), globalTrack.phi());
         if (dR2 < 0.02 * 0.02 and std::abs(tk->pt() - globalTrack.pt()) < 0.001 * tk->pt()) {
           useThisLink = true;
         }
  
         if (!useThisLink)
           continue;
  
         // Check whether the stand-alone track matches a L2, if not, we will re-use this L2
         const reco::TrackRef staTrack = link->standAloneTrack();
         if (l2muRef == staTrack)
           re_do_this_L2 = false;
  
         // Check the quality of the reconstructed L3, if poor quality, we will re-use this L2
         if (staTrack == l2muRef && applyL3Filters_) {
           re_do_this_L2 = true;
           const reco::Track& globalTrack = *link->globalTrack();
           if (globalTrack.numberOfValidHits() < min_Nhits_)
             continue;  // cut on number of hits
           if (globalTrack.normalizedChi2() > max_NormalizedChi2_)
             continue;  //normalizedChi2 cut
           if (globalTrack.hitPattern().numberOfValidMuonHits() < min_NmuonHits_)
             continue;  //min muon hits cut
           if (std::abs(globalTrack.pt() - l2muRef->pt()) > max_PtDifference_ * globalTrack.pt())
             continue;  // pt difference
           re_do_this_L2 = false;
         }
       }
     }
     if (re_do_this_L2)
       result->push_back(*l2muRef);  // used the L2 if no L3 if matched or if the matched L3 has poor quality cuts.
   }
   iEvent.put(std::move(result));
 }

References funct::abs(), applyL3Filters_, HLTMuonOfflineAnalyzer_cfi::deltaR2, cutBasedMuonId_MuonPOG_V0_cff::globalTrack, iEvent, cmsLHEtoEOSManager::l, l2Src_, l3linkToken_, l3OISrc_, MainPageGenerator::link, electronStore::links, max_NormalizedChi2_, max_PtDifference_, metname, min_Nhits_, min_NmuonHits_, eostools::move(), mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

◆ applyL3Filters_

const bool HLTMuonL2SelectorForL3IO::applyL3Filters_

private

Definition at line 45 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ l2Src_

const edm::EDGetTokenT<reco::TrackCollection> HLTMuonL2SelectorForL3IO::l2Src_

private

Definition at line 42 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ l3linkToken_

const edm::EDGetTokenT<reco::MuonTrackLinksCollection> HLTMuonL2SelectorForL3IO::l3linkToken_

private

Definition at line 44 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ l3OISrc_

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTMuonL2SelectorForL3IO::l3OISrc_

private

Definition at line 43 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ max_NormalizedChi2_

const double HLTMuonL2SelectorForL3IO::max_NormalizedChi2_

private

Definition at line 46 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ max_PtDifference_

const double HLTMuonL2SelectorForL3IO::max_PtDifference_

private

Definition at line 46 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ min_Nhits_

const int HLTMuonL2SelectorForL3IO::min_Nhits_

private

Definition at line 47 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

◆ min_NmuonHits_

const int HLTMuonL2SelectorForL3IO::min_NmuonHits_

private

Definition at line 47 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HLTMuonL2SelectorForL3IO()

◆ ~HLTMuonL2SelectorForL3IO()

Member Function Documentation

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ applyL3Filters_

◆ l2Src_

◆ l3linkToken_

◆ l3OISrc_

◆ max_NormalizedChi2_

◆ max_PtDifference_

◆ min_Nhits_

◆ min_NmuonHits_