L1MuonMatcherAlgoT< Tr > Class Template Reference

#include <L1MuonMatcherAlgo.h>

Public Member Functions
TrajectoryStateOnSurface	extrapolate (const reco::Track &tk) const
	Extrapolate reco::Track to the muon station 2, return an invalid TSOS if it fails. More...
TrajectoryStateOnSurface	extrapolate (const reco::Candidate &tk) const
	Extrapolate reco::Candidate to the muon station 2, return an invalid TSOS if it fails. More...
TrajectoryStateOnSurface	extrapolate (const SimTrack &tk, const edm::SimVertexContainer &vtx) const
	Extrapolate a SimTrack to the muon station 2, return an invalid TSOS if it fails. Requires SimVertices to know where to start from. More...
TrajectoryStateOnSurface	extrapolate (const FreeTrajectoryState &state) const
	Extrapolate a FreeTrajectoryState to the muon station 2, return an invalid TSOS if it fails. More...
void	init (const edm::EventSetup &iSetup)
	Call this method at the beginning of each run, to initialize geometry, magnetic field and propagators. More...
	L1MuonMatcherAlgoT (const edm::ParameterSet &iConfig, edm::ConsumesCollector)
bool	match (const reco::Track &tk, const l1extra::L1MuonParticle &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
bool	match (const reco::Candidate &c, const l1extra::L1MuonParticle &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
bool	match (const SimTrack &tk, const edm::SimVertexContainer &vtxs, const l1extra::L1MuonParticle &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
bool	match (TrajectoryStateOnSurface &propagated, const l1extra::L1MuonParticle &l1, float &deltaR, float &deltaPhi) const
int	match (const reco::Track &tk, const std::vector< l1extra::L1MuonParticle > &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
int	match (const reco::Candidate &c, const std::vector< l1extra::L1MuonParticle > &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
int	match (const SimTrack &tk, const edm::SimVertexContainer &vtxs, const std::vector< l1extra::L1MuonParticle > &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
int	match (TrajectoryStateOnSurface &propagated, const std::vector< l1extra::L1MuonParticle > &l1, float &deltaR, float &deltaPhi) const
int	match (const reco::Track &tk, const std::vector< l1t::Muon > &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
int	match (const reco::Candidate &c, const std::vector< l1t::Muon > &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
int	match (const SimTrack &tk, const edm::SimVertexContainer &vtxs, const std::vector< l1t::Muon > &l1, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
int	match (TrajectoryStateOnSurface &propagated, const std::vector< l1t::Muon > &l1, float &deltaR, float &deltaPhi) const
template<typename Collection , typename Selector >
int	matchGeneric (const reco::Track &tk, const Collection &l1, const Selector &sel, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
template<typename Collection , typename Selector >
int	matchGeneric (const reco::Candidate &c, const Collection &l1, const Selector &sel, float &deltaR, float &deltaPhi, TrajectoryStateOnSurface &propagated) const
template<typename Collection , typename Selector >
int	matchGeneric (TrajectoryStateOnSurface &propagated, const Collection &l1, const Selector &sel, float &deltaR, float &deltaPhi) const
PropagateToMuon &	propagatorToMuon ()
	Return the propagator to second muon station (in case it's needed) More...
const PropagateToMuon &	propagatorToMuon () const
	Return the propagator to second muon station (in case it's needed) More...
void	setL1PhiOffset (double l1PhiOffset)
	Add this offset to the L1 phi before doing the match, to correct for different scales in L1 vs offline. More...
	~L1MuonMatcherAlgoT ()=default

Private Types
typedef StringCutObjectSelector< reco::Candidate, true >	L1Selector
enum	SortBy {   SortByDeltaR = 0, SortByDeltaPhi, SortByDeltaEta, SortByPt,   SortByQual }
	Sort by deltaPhi or deltaEta instead of deltaR. More...

Private Member Functions
template<class T >
int	genericQuality (T const &cand) const
int	genericQuality (l1extra::L1MuonParticle const &cand) const
int	genericQuality (l1t::Muon const &cand) const

Private Attributes
double	deltaEta_
double	deltaPhi_
double	deltaR2_
	Matching cuts. More...
double	l1PhiOffset_
	offset to be added to the L1 phi before the match More...
L1Selector	preselectionCut_
	Preselection cut to apply to L1 candidates before matching. More...
PropagateToMuon	prop_
PropagateToMuonSetupT< Tr >	propSetup_
SortBy	sortBy_
bool	useStage2L1_

Detailed Description

template<edm::Transition Tr>
class L1MuonMatcherAlgoT< Tr >

Definition at line 27 of file L1MuonMatcherAlgo.h.

Member Typedef Documentation

L1Selector

template<edm::Transition Tr>

typedef StringCutObjectSelector<reco::Candidate, true> L1MuonMatcherAlgoT< Tr >::L1Selector

private

Definition at line 249 of file L1MuonMatcherAlgo.h.

Member Enumeration Documentation

SortBy

template<edm::Transition Tr>

enum L1MuonMatcherAlgoT::SortBy

private

Sort by deltaPhi or deltaEta instead of deltaR.

Enumerator
SortByDeltaR
SortByDeltaPhi
SortByDeltaEta
SortByPt
SortByQual

Definition at line 257 of file L1MuonMatcherAlgo.h.

{ SortByDeltaR = 0, SortByDeltaPhi, SortByDeltaEta, SortByPt, SortByQual };

Constructor & Destructor Documentation

L1MuonMatcherAlgoT()

template<edm::Transition Tr>

L1MuonMatcherAlgoT< Tr >::L1MuonMatcherAlgoT ( const edm::ParameterSet &  iConfig, edm::ConsumesCollector  iCollector  )

explicit

Definition at line 326 of file L1MuonMatcherAlgo.h.

    : propSetup_(iConfig, iCollector),
      useStage2L1_(iConfig.existsAs<bool>("useStage2L1") ? iConfig.getParameter<bool>("useStage2L1") : false),
      preselectionCut_(
          (iConfig.existsAs<std::string>("preselection")) ? iConfig.getParameter<std::string>("preselection") : ""),
      deltaR2_(std::pow(iConfig.getParameter<double>("maxDeltaR"), 2)),
      deltaPhi_(iConfig.existsAs<double>("maxDeltaPhi") ? iConfig.getParameter<double>("maxDeltaPhi") : 10),
      deltaEta_(iConfig.existsAs<double>("maxDeltaEta") ? iConfig.getParameter<double>("maxDeltaEta") : 10),
      l1PhiOffset_(iConfig.existsAs<double>("l1PhiOffset") ? iConfig.getParameter<double>("l1PhiOffset") : 0) {
  bool reqQual = iConfig.existsAs<bool>("sortByQual") && iConfig.getParameter<bool>("sortByQual");
  bool reqPhi = iConfig.existsAs<bool>("sortByDeltaPhi") && iConfig.getParameter<bool>("sortByDeltaPhi");
  bool reqEta = iConfig.existsAs<bool>("sortByDeltaEta") && iConfig.getParameter<bool>("sortByDeltaEta");
  bool reqPt = iConfig.existsAs<bool>("sortByPt") && iConfig.getParameter<bool>("sortByPt");
  std::string sortBy = iConfig.existsAs<std::string>("sortBy") ? iConfig.getParameter<std::string>("sortBy") : "";
  if (reqPhi + reqEta + reqPt + reqQual > 1)
    throw cms::Exception("Configuration")
        << "L1MuonMatcherAlgoT: Can't set more than one 'sortBy<XXX>' parameter to True.\n";
  if (sortBy == "deltaPhi") {
    if (reqEta || reqPt || reqQual)
      throw cms::Exception("Configuration") << "L1MuonMatcherAlgoT: Can't set sortBy = 'deltaPhi' and set also another "
                                               "'sortBy<XXX>' parameter to True.\n";
    else
      reqPhi = true;
  }
  if (sortBy == "deltaEta") {
    if (reqPhi || reqPt || reqQual)
      throw cms::Exception("Configuration") << "L1MuonMatcherAlgoT: Can't set sortBy = 'deltaEta' and set also another "
                                               "'sortBy<XXX>' parameter to True.\n";
    else
      reqEta = true;
  }
  if (sortBy == "pt") {
    if (reqPhi || reqEta || reqQual)
      throw cms::Exception("Configuration")
          << "L1MuonMatcherAlgoT: Can't set sortBy = 'pt' and set also another 'sortBy<XXX>' parameter to True.\n";
    else
      reqPt = true;
  }
  if (sortBy == "quality") {
    if (reqPhi || reqEta || reqPt)
      throw cms::Exception("Configuration")
          << "L1MuonMatcherAlgoT: Can't set sortBy = 'quality' and set also another 'sortBy<XXX>' parameter to True.\n";
    else
      reqQual = true;
  }
  if (sortBy == "deltaR") {
    if (reqPhi || reqEta || reqPt || reqQual)
      throw cms::Exception("Configuration")
          << "L1MuonMatcherAlgoT: Can't set sortBy = 'deltaR' and set also another 'sortBy<XXX>' parameter to True.\n";
  }
  // so, if we're here there's no ambiguity in what the user may want. either everything is false, or exactly one req is true.
  if (reqEta)
    sortBy_ = SortByDeltaEta;
  else if (reqPhi)
    sortBy_ = SortByDeltaPhi;
  else if (reqQual)
    sortBy_ = SortByQual;
  else if (reqPt)
    sortBy_ = SortByPt;
  else
    sortBy_ = SortByDeltaR;
}

~L1MuonMatcherAlgoT()

template<edm::Transition Tr>

L1MuonMatcherAlgoT< Tr >::~L1MuonMatcherAlgoT ( )

default

Member Function Documentation

extrapolate() [1/4]

template<edm::Transition Tr>

TrajectoryStateOnSurface L1MuonMatcherAlgoT< Tr >::extrapolate ( const reco::Track &  tk ) const

inline

Extrapolate reco::Track to the muon station 2, return an invalid TSOS if it fails.

Definition at line 36 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgoT< edm::Transition::BeginRun >::match(), and L1MuonMatcherAlgoT< edm::Transition::BeginRun >::matchGeneric().

{ return prop_.extrapolate(tk); }

extrapolate() [2/4]

template<edm::Transition Tr>

TrajectoryStateOnSurface L1MuonMatcherAlgoT< Tr >::extrapolate ( const reco::Candidate &  tk ) const

inline

Extrapolate reco::Candidate to the muon station 2, return an invalid TSOS if it fails.

Definition at line 39 of file L1MuonMatcherAlgo.h.

{ return prop_.extrapolate(tk); }

extrapolate() [3/4]

template<edm::Transition Tr>

TrajectoryStateOnSurface L1MuonMatcherAlgoT< Tr >::extrapolate ( const SimTrack &  tk, const edm::SimVertexContainer &  vtx  ) const

inline

Extrapolate a SimTrack to the muon station 2, return an invalid TSOS if it fails. Requires SimVertices to know where to start from.

Definition at line 42 of file L1MuonMatcherAlgo.h.

                                                                                                   {
    return prop_.extrapolate(tk, vtx);
  }

extrapolate() [4/4]

template<edm::Transition Tr>

TrajectoryStateOnSurface L1MuonMatcherAlgoT< Tr >::extrapolate ( const FreeTrajectoryState &  state ) const

inline

Extrapolate a FreeTrajectoryState to the muon station 2, return an invalid TSOS if it fails.

Definition at line 47 of file L1MuonMatcherAlgo.h.

{ return prop_.extrapolate(state); }

genericQuality() [1/3]

template<edm::Transition Tr>

template<class T >

int L1MuonMatcherAlgoT< Tr >::genericQuality ( T const &  cand ) const

inlineprivate

Definition at line 237 of file L1MuonMatcherAlgo.h.

                                          {
    return 0;
  }

genericQuality() [2/3]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::genericQuality ( l1extra::L1MuonParticle const &  cand ) const

inlineprivate

Definition at line 241 of file L1MuonMatcherAlgo.h.

{ return cand.gmtMuonCand().rank(); }

genericQuality() [3/3]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::genericQuality ( l1t::Muon const &  cand ) const

inlineprivate

Definition at line 242 of file L1MuonMatcherAlgo.h.

{ return cand.hwQual(); }

init()

template<edm::Transition Tr>

void L1MuonMatcherAlgoT< Tr >::init

(

const edm::EventSetup &

iSetup

)

Call this method at the beginning of each run, to initialize geometry, magnetic field and propagators.

Definition at line 390 of file L1MuonMatcherAlgo.h.

Referenced by HLTMuonPlotter::beginRun().

                                                             {
  prop_ = propSetup_.init(iSetup);
}

match() [1/12]

template<edm::Transition Tr>

bool L1MuonMatcherAlgoT< Tr >::match ( const reco::Track &  tk, const l1extra::L1MuonParticle &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Try to match one track to one L1. Return true if succeeded (and update deltaR, deltaPhi and propagated TSOS accordingly) The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 56 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgoT< edm::Transition::BeginRun >::match().

                                                         {
    propagated = extrapolate(tk);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : false;
  }

match() [2/12]

template<edm::Transition Tr>

bool L1MuonMatcherAlgoT< Tr >::match ( const reco::Candidate &  c, const l1extra::L1MuonParticle &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Try to match one track to one L1. Return true if succeeded (and update deltaR, deltaPhi and propagated TSOS accordingly) The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 67 of file L1MuonMatcherAlgo.h.

                                                         {
    propagated = extrapolate(c);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : false;
  }

match() [3/12]

template<edm::Transition Tr>

bool L1MuonMatcherAlgoT< Tr >::match ( const SimTrack &  tk, const edm::SimVertexContainer &  vtxs, const l1extra::L1MuonParticle &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Try to match one simtrack to one L1. Return true if succeeded (and update deltaR, deltaPhi and propagated TSOS accordingly) The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 78 of file L1MuonMatcherAlgo.h.

                                                         {
    propagated = extrapolate(tk, vtxs);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : false;
  }

match() [4/12]

template<edm::Transition Tr>

bool L1MuonMatcherAlgoT< Tr >::match	(	TrajectoryStateOnSurface &	propagated,
		const l1extra::L1MuonParticle &	l1,
		float &	deltaR,
		float &	deltaPhi
	)		const

Try to match one track to one L1. Return true if succeeded (and update deltaR, deltaPhi accordingly) The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 395 of file L1MuonMatcherAlgo.h.

                                                          {
  if (preselectionCut_(l1)) {
    GlobalPoint pos = propagated.globalPosition();
    double thisDeltaPhi = ::deltaPhi(double(pos.phi()), l1.phi() + l1PhiOffset_);
    double thisDeltaEta = pos.eta() - l1.eta();
    double thisDeltaR2 = ::deltaR2(double(pos.eta()), double(pos.phi()), l1.eta(), l1.phi() + l1PhiOffset_);
    if ((fabs(thisDeltaPhi) < deltaPhi_) && (fabs(thisDeltaEta) < deltaEta_) && (thisDeltaR2 < deltaR2_)) {
      deltaR = std::sqrt(thisDeltaR2);
      deltaPhi = thisDeltaPhi;
      return true;
    }
  }
  return false;
}

match() [5/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match ( const reco::Track &  tk, const std::vector< l1extra::L1MuonParticle > &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 100 of file L1MuonMatcherAlgo.h.

                                                        {
    propagated = extrapolate(tk);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : -1;
  }

match() [6/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match ( const reco::Candidate &  c, const std::vector< l1extra::L1MuonParticle > &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 112 of file L1MuonMatcherAlgo.h.

                                                        {
    propagated = extrapolate(c);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : -1;
  }

match() [7/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match ( const SimTrack &  tk, const edm::SimVertexContainer &  vtxs, const std::vector< l1extra::L1MuonParticle > &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 124 of file L1MuonMatcherAlgo.h.

                                                        {
    propagated = extrapolate(tk, vtxs);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : -1;
  }

match() [8/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match	(	TrajectoryStateOnSurface &	propagated,
		const std::vector< l1extra::L1MuonParticle > &	l1,
		float &	deltaR,
		float &	deltaPhi
	)		const

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi accordingly) Returns -1 if the match fails The preselection cut on L1, if specified in the config, is applied before the match

Definition at line 414 of file L1MuonMatcherAlgo.h.

                                                         {
  return matchGeneric(propagated, l1s, preselectionCut_, deltaR, deltaPhi);
}

match() [9/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match ( const reco::Track &  tk, const std::vector< l1t::Muon > &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to stage2 L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails The preselection cut on stage2 L1, if specified in the config, is applied before the match

Definition at line 147 of file L1MuonMatcherAlgo.h.

                                                        {
    propagated = extrapolate(tk);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : -1;
  }

match() [10/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match ( const reco::Candidate &  c, const std::vector< l1t::Muon > &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to stage2 L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails The preselection cut on stage2 L1, if specified in the config, is applied before the match

Definition at line 159 of file L1MuonMatcherAlgo.h.

                                                        {
    propagated = extrapolate(c);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : -1;
  }

match() [11/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match ( const SimTrack &  tk, const edm::SimVertexContainer &  vtxs, const std::vector< l1t::Muon > &  l1, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to stage2 L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails The preselection cut on stage 2 L1, if specified in the config, is applied before the match

Definition at line 171 of file L1MuonMatcherAlgo.h.

                                                        {
    propagated = extrapolate(tk, vtxs);
    return propagated.isValid() ? match(propagated, l1, deltaR, deltaPhi) : -1;
  }

match() [12/12]

template<edm::Transition Tr>

int L1MuonMatcherAlgoT< Tr >::match	(	TrajectoryStateOnSurface &	propagated,
		const std::vector< l1t::Muon > &	l1,
		float &	deltaR,
		float &	deltaPhi
	)		const

Find the best match to stage 2 L1, and return its index in the vector (and update deltaR, deltaPhi accordingly) Returns -1 if the match fails The preselection cut on stage 2 L1, if specified in the config, is applied before the match

Definition at line 422 of file L1MuonMatcherAlgo.h.

                                                         {
  return matchGeneric(propagated, l1s, preselectionCut_, deltaR, deltaPhi);
}

matchGeneric() [1/3]

template<edm::Transition Tr>

template<typename Collection , typename Selector >

int L1MuonMatcherAlgoT< Tr >::matchGeneric ( const reco::Track &  tk, const Collection &  l1, const Selector &  sel, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails Only the objects passing the selector will be allowed for the match. If you don't need a selector, just use an AnySelector (CommonTools/Utils) which accepts everything

Definition at line 196 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgoT< edm::Transition::BeginRun >::matchGeneric().

                                                               {
    propagated = extrapolate(tk);
    return propagated.isValid() ? matchGeneric(propagated, l1, sel, deltaR, deltaPhi) : -1;
  }

matchGeneric() [2/3]

template<edm::Transition Tr>

template<typename Collection , typename Selector >

int L1MuonMatcherAlgoT< Tr >::matchGeneric ( const reco::Candidate &  c, const Collection &  l1, const Selector &  sel, float &  deltaR, float &  deltaPhi, TrajectoryStateOnSurface &  propagated  ) const

inline

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi and propagated TSOS accordingly) Returns -1 if the match fails Only the objects passing the selector will be allowed for the match. If you don't need a selector, just use an AnySelector (CommonTools/Utils) which accepts everything

Definition at line 211 of file L1MuonMatcherAlgo.h.

                                                               {
    propagated = extrapolate(c);
    return propagated.isValid() ? matchGeneric(propagated, l1, sel, deltaR, deltaPhi) : -1;
  }

matchGeneric() [3/3]

template<edm::Transition Tr>

template<typename Collection , typename Selector >

int L1MuonMatcherAlgoT< Tr >::matchGeneric	(	TrajectoryStateOnSurface &	propagated,
		const Collection &	l1,
		const Selector &	sel,
		float &	deltaR,
		float &	deltaPhi
	)		const

Find the best match to L1, and return its index in the vector (and update deltaR, deltaPhi accordingly) Returns -1 if the match fails Only the objects passing the selector will be allowed for the match. The selector defaults to an AnySelector (CommonTools/Utils) which just accepts everything

Definition at line 266 of file L1MuonMatcherAlgo.h.

                                                                {
  typedef typename Collection::value_type obj;

  int match = -1;
  double minDeltaPhi = deltaPhi_;
  double minDeltaEta = deltaEta_;
  double minDeltaR2 = deltaR2_;
  double maxPt = -1.0;
  int maxQual = 0;
  GlobalPoint pos = propagated.globalPosition();
  for (int i = 0, n = l1s.size(); i < n; ++i) {
    const obj &l1 = l1s[i];
    if (sel(l1)) {
      double thisDeltaPhi = ::deltaPhi(double(pos.phi()), l1.phi() + l1PhiOffset_);
      double thisDeltaEta = pos.eta() - l1.eta();
      double thisDeltaR2 = ::deltaR2(double(pos.eta()), double(pos.phi()), l1.eta(), l1.phi() + l1PhiOffset_);
      int thisQual = genericQuality<obj>(l1);
      double thisPt = l1.pt();
      if ((fabs(thisDeltaPhi) < deltaPhi_) && (fabs(thisDeltaEta) < deltaEta_) &&
          (thisDeltaR2 < deltaR2_)) {  // check all
        bool betterMatch = (match == -1);
        switch (sortBy_) {
          case SortByDeltaR:
            betterMatch = (thisDeltaR2 < minDeltaR2);
            break;
          case SortByDeltaEta:
            betterMatch = (fabs(thisDeltaEta) < fabs(minDeltaEta));
            break;
          case SortByDeltaPhi:
            betterMatch = (fabs(thisDeltaPhi) < fabs(minDeltaPhi));
            break;
          case SortByPt:
            betterMatch = (thisPt > maxPt);
            break;
          // Quality is an int, adding sorting by pT in case of identical qualities
          case SortByQual:
            betterMatch = (thisQual > maxQual || ((thisQual == maxQual) && (thisPt > maxPt)));
            break;
        }
        if (betterMatch) {  // sort on one
          match = i;
          deltaR = std::sqrt(thisDeltaR2);
          deltaPhi = thisDeltaPhi;
          minDeltaR2 = thisDeltaR2;
          minDeltaEta = thisDeltaEta;
          minDeltaPhi = thisDeltaPhi;
          maxQual = thisQual;
          maxPt = thisPt;
        }
      }
    }
  }
  return match;
}

propagatorToMuon() [1/2]

template<edm::Transition Tr>

PropagateToMuon& L1MuonMatcherAlgoT< Tr >::propagatorToMuon ( )

inline

Return the propagator to second muon station (in case it's needed)

Definition at line 50 of file L1MuonMatcherAlgo.h.

{ return prop_; }

propagatorToMuon() [2/2]

template<edm::Transition Tr>

const PropagateToMuon& L1MuonMatcherAlgoT< Tr >::propagatorToMuon ( ) const

inline

Return the propagator to second muon station (in case it's needed)

Definition at line 52 of file L1MuonMatcherAlgo.h.

{ return prop_; }

setL1PhiOffset()

template<edm::Transition Tr>

void L1MuonMatcherAlgoT< Tr >::setL1PhiOffset ( double  l1PhiOffset )

inline

Add this offset to the L1 phi before doing the match, to correct for different scales in L1 vs offline.

Definition at line 233 of file L1MuonMatcherAlgo.h.

{ l1PhiOffset_ = l1PhiOffset; }

Member Data Documentation

deltaEta_

template<edm::Transition Tr>

double L1MuonMatcherAlgoT< Tr >::deltaEta_

private

Definition at line 254 of file L1MuonMatcherAlgo.h.

deltaPhi_

template<edm::Transition Tr>

double L1MuonMatcherAlgoT< Tr >::deltaPhi_

private

Definition at line 254 of file L1MuonMatcherAlgo.h.

deltaR2_

template<edm::Transition Tr>

double L1MuonMatcherAlgoT< Tr >::deltaR2_

private

Matching cuts.

Definition at line 254 of file L1MuonMatcherAlgo.h.

l1PhiOffset_

template<edm::Transition Tr>

double L1MuonMatcherAlgoT< Tr >::l1PhiOffset_

private

offset to be added to the L1 phi before the match

Definition at line 261 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgoT< edm::Transition::BeginRun >::setL1PhiOffset().

preselectionCut_

template<edm::Transition Tr>

L1Selector L1MuonMatcherAlgoT< Tr >::preselectionCut_

private

Preselection cut to apply to L1 candidates before matching.

Definition at line 251 of file L1MuonMatcherAlgo.h.

prop_

template<edm::Transition Tr>

PropagateToMuon L1MuonMatcherAlgoT< Tr >::prop_

private

Definition at line 245 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgoT< edm::Transition::BeginRun >::extrapolate(), and L1MuonMatcherAlgoT< edm::Transition::BeginRun >::propagatorToMuon().

propSetup_

template<edm::Transition Tr>

PropagateToMuonSetupT<Tr> L1MuonMatcherAlgoT< Tr >::propSetup_

private

Definition at line 244 of file L1MuonMatcherAlgo.h.

sortBy_

template<edm::Transition Tr>

SortBy L1MuonMatcherAlgoT< Tr >::sortBy_

private

Definition at line 258 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgoT< edm::Transition::BeginRun >::L1MuonMatcherAlgoT().

useStage2L1_

template<edm::Transition Tr>

bool L1MuonMatcherAlgoT< Tr >::useStage2L1_

private

Definition at line 247 of file L1MuonMatcherAlgo.h.