HLTMuonL2SelectorForL3IO Class Reference

#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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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.

References LogTrace.

    : 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>();
}

~HLTMuonL2SelectorForL3IO()

HLTMuonL2SelectorForL3IO::~HLTMuonL2SelectorForL3IO ( )

override

destructor

Definition at line 31 of file HLTMuonL2SelectorForL3IO.cc.

{}

Member Function Documentation

fillDescriptions()

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

static

default values

Definition at line 100 of file HLTMuonL2SelectorForL3IO.cc.

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

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

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.

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.

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

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().