L1MuonMatcherAlgo Class Reference

Matcher of reconstructed objects to L1 Muons. More...

#include "MuonAnalysis/MuonAssociators/interface/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...
	L1MuonMatcherAlgo (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...
	~L1MuonMatcherAlgo ()

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 (const T &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_
SortBy	sortBy_
bool	useStage2L1_

Detailed Description

Matcher of reconstructed objects to L1 Muons.

"HLTriggerOffline/Muon/interface/L1MuonMatcherAlgo.h"

Author

Giovanni Petrucciani

Giovanni Petrucciani

Version

Id:

L1MuonMatcherAlgo.h,v 1.8 2010/07/01 07:40:52 gpetrucc Exp

Definition at line 29 of file L1MuonMatcherAlgo.h.

Member Typedef Documentation

typedef StringCutObjectSelector<reco::Candidate, true> L1MuonMatcherAlgo::L1Selector

private

Definition at line 245 of file L1MuonMatcherAlgo.h.

Member Enumeration Documentation

enum L1MuonMatcherAlgo::SortBy

private

Sort by deltaPhi or deltaEta instead of deltaR.

Enumerator
SortByDeltaR
SortByDeltaPhi
SortByDeltaEta
SortByPt
SortByQual

Definition at line 253 of file L1MuonMatcherAlgo.h.

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

Constructor & Destructor Documentation

L1MuonMatcherAlgo::L1MuonMatcherAlgo ( const edm::ParameterSet &  iConfig, edm::ConsumesCollector  iCollector  )

explicit

Definition at line 5 of file L1MuonMatcherAlgo.cc.

References Exception, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), sortBy_, SortByDeltaEta, SortByDeltaPhi, SortByDeltaR, SortByPt, SortByQual, and AlCaHLTBitMon_QueryRunRegistry::string.

    : prop_(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")
        << "L1MuonMatcherAlgo: Can't set more than one 'sortBy<XXX>' parameter to True.\n";
  if (sortBy == "deltaPhi") {
    if (reqEta || reqPt || reqQual)
      throw cms::Exception("Configuration")
          << "L1MuonMatcherAlgo: 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")
          << "L1MuonMatcherAlgo: 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")
          << "L1MuonMatcherAlgo: 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")
          << "L1MuonMatcherAlgo: 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")
          << "L1MuonMatcherAlgo: 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;
}

L1MuonMatcherAlgo::~L1MuonMatcherAlgo

(

)

Definition at line 68 of file L1MuonMatcherAlgo.cc.

{}

Member Function Documentation

TrajectoryStateOnSurface L1MuonMatcherAlgo::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 38 of file L1MuonMatcherAlgo.h.

References PropagateToMuon::extrapolate(), and prop_.

Referenced by match(), and matchGeneric().

{ return prop_.extrapolate(tk); }

TrajectoryStateOnSurface L1MuonMatcherAlgo::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 41 of file L1MuonMatcherAlgo.h.

References PropagateToMuon::extrapolate(), and prop_.

{ return prop_.extrapolate(tk); }

TrajectoryStateOnSurface L1MuonMatcherAlgo::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 44 of file L1MuonMatcherAlgo.h.

References PropagateToMuon::extrapolate(), and prop_.

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

TrajectoryStateOnSurface L1MuonMatcherAlgo::extrapolate ( const FreeTrajectoryState &  state ) const

inline

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

Definition at line 49 of file L1MuonMatcherAlgo.h.

References PropagateToMuon::extrapolate(), and prop_.

{ return prop_.extrapolate(state); }

template<class T >

int L1MuonMatcherAlgo::genericQuality ( const T &  cand ) const

inlineprivate

Definition at line 271 of file L1MuonMatcherAlgo.h.

                                                                {
  return 0;
}

void L1MuonMatcherAlgo::init

(

const edm::EventSetup &

iSetup

)

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

Definition at line 70 of file L1MuonMatcherAlgo.cc.

References PropagateToMuon::init(), and prop_.

Referenced by pat::HLTL1MuonMatcher::produce().

{ prop_.init(iSetup); }

bool L1MuonMatcherAlgo::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 58 of file L1MuonMatcherAlgo.h.

References extrapolate(), and TrajectoryStateOnSurface::isValid().

Referenced by match(), and matchGeneric().

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

bool L1MuonMatcherAlgo::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 69 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

bool L1MuonMatcherAlgo::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 80 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

bool L1MuonMatcherAlgo::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 72 of file L1MuonMatcherAlgo.cc.

References deltaEta_, srCondWrite_cfg::deltaPhi, deltaPhi_, reco::deltaR2(), deltaR2_, PV3DBase< T, PVType, FrameType >::eta(), reco::LeafCandidate::eta(), TrajectoryStateOnSurface::globalPosition(), l1PhiOffset_, PV3DBase< T, PVType, FrameType >::phi(), reco::LeafCandidate::phi(), preselectionCut_, and mathSSE::sqrt().

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

int L1MuonMatcherAlgo::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 102 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

int L1MuonMatcherAlgo::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 114 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

int L1MuonMatcherAlgo::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 126 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

int L1MuonMatcherAlgo::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 90 of file L1MuonMatcherAlgo.cc.

References matchGeneric(), and preselectionCut_.

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

int L1MuonMatcherAlgo::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 149 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

int L1MuonMatcherAlgo::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 161 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

int L1MuonMatcherAlgo::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 173 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and match().

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

int L1MuonMatcherAlgo::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 97 of file L1MuonMatcherAlgo.cc.

References matchGeneric(), and preselectionCut_.

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

template<typename Collection , typename Selector >

int L1MuonMatcherAlgo::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 198 of file L1MuonMatcherAlgo.h.

References extrapolate(), and TrajectoryStateOnSurface::isValid().

Referenced by match(), matchGeneric(), and pat::HLTL1MuonMatcher::produce().

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

template<typename Collection , typename Selector >

int L1MuonMatcherAlgo::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 213 of file L1MuonMatcherAlgo.h.

References extrapolate(), TrajectoryStateOnSurface::isValid(), and matchGeneric().

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

template<typename Collection , typename Selector >

int L1MuonMatcherAlgo::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 276 of file L1MuonMatcherAlgo.h.

References deltaEta_, srCondWrite_cfg::deltaPhi, deltaPhi_, reco::deltaR2(), deltaR2_, PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, l1PhiOffset_, match(), HLT_FULL_cff::minDeltaEta, dqmiodumpmetadata::n, getGTfromDQMFile::obj, PV3DBase< T, PVType, FrameType >::phi(), EgammaValidation_Wenu_cff::sel, sortBy_, SortByDeltaEta, SortByDeltaPhi, SortByDeltaR, SortByQual, and mathSSE::sqrt().

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

PropagateToMuon& L1MuonMatcherAlgo::propagatorToMuon ( )

inline

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

Definition at line 52 of file L1MuonMatcherAlgo.h.

References prop_.

{ return prop_; }

const PropagateToMuon& L1MuonMatcherAlgo::propagatorToMuon ( ) const

inline

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

Definition at line 54 of file L1MuonMatcherAlgo.h.

References prop_.

{ return prop_; }

void L1MuonMatcherAlgo::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 235 of file L1MuonMatcherAlgo.h.

References l1PhiOffset_.

{ l1PhiOffset_ = l1PhiOffset; }

Member Data Documentation

double L1MuonMatcherAlgo::deltaEta_

private

Definition at line 250 of file L1MuonMatcherAlgo.h.

Referenced by match(), and matchGeneric().

double L1MuonMatcherAlgo::deltaPhi_

private

Definition at line 250 of file L1MuonMatcherAlgo.h.

Referenced by match(), and matchGeneric().

double L1MuonMatcherAlgo::deltaR2_

private

Matching cuts.

Definition at line 250 of file L1MuonMatcherAlgo.h.

Referenced by match(), and matchGeneric().

double L1MuonMatcherAlgo::l1PhiOffset_

private

offset to be added to the L1 phi before the match

Definition at line 257 of file L1MuonMatcherAlgo.h.

Referenced by match(), matchGeneric(), and setL1PhiOffset().

L1Selector L1MuonMatcherAlgo::preselectionCut_

private

Preselection cut to apply to L1 candidates before matching.

Definition at line 247 of file L1MuonMatcherAlgo.h.

Referenced by match().

PropagateToMuon L1MuonMatcherAlgo::prop_

private

Definition at line 241 of file L1MuonMatcherAlgo.h.

Referenced by extrapolate(), init(), and propagatorToMuon().

SortBy L1MuonMatcherAlgo::sortBy_

private

Definition at line 254 of file L1MuonMatcherAlgo.h.

Referenced by L1MuonMatcherAlgo(), and matchGeneric().

bool L1MuonMatcherAlgo::useStage2L1_

private

Definition at line 243 of file L1MuonMatcherAlgo.h.