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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () 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

(

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

override

destructor

Definition at line 31 of file HLTMuonL2SelectorForL3IO.cc.

{}

Member Function Documentation

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

static

default values

Definition at line 100 of file HLTMuonL2SelectorForL3IO.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_2018_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);
}

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, reco::TrackBase::eta(), edm::Event::getByToken(), reco::MuonTrackLinks::globalTrack(), cutBasedMuonId_MuonPOG_V0_cff::globalTrack, reco::TrackBase::hitPattern(), cmsLHEtoEOSManager::l, l2Src_, l3linkToken_, l3OISrc_, MainPageGenerator::link, electronStore::links, max_NormalizedChi2_, max_PtDifference_, metname, min_Nhits_, min_NmuonHits_, eostools::move(), reco::TrackBase::normalizedChi2(), reco::TrackBase::numberOfValidHits(), reco::HitPattern::numberOfValidMuonHits(), reco::TrackBase::phi(), reco::TrackBase::pt(), edm::Event::put(), mps_fire::result, reco::MuonTrackLinks::standAloneTrack(), 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

const bool HLTMuonL2SelectorForL3IO::applyL3Filters_

private

Definition at line 45 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

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

private

Definition at line 42 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

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

private

Definition at line 44 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

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

private

Definition at line 43 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const double HLTMuonL2SelectorForL3IO::max_NormalizedChi2_

private

Definition at line 46 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const double HLTMuonL2SelectorForL3IO::max_PtDifference_

private

Definition at line 46 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const int HLTMuonL2SelectorForL3IO::min_Nhits_

private

Definition at line 47 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const int HLTMuonL2SelectorForL3IO::min_NmuonHits_

private

Definition at line 47 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().