hltriggeroffline::L1MuonMatcherAlgo Class Reference

#include <L1MuonMatcherAlgo.h>

Public Types
using	ESTokens = PropagateToMuon::ESTokens

Public Member Functions
TrajectoryStateOnSurface	extrapolate (const reco::Track &tk) const
TrajectoryStateOnSurface	extrapolate (const reco::Candidate &tk) const
TrajectoryStateOnSurface	extrapolate (const FreeTrajectoryState &state) const
void	init (const edm::EventSetup &iSetup)
	L1MuonMatcherAlgo (const edm::ParameterSet &iConfig, const ESTokens &iTokens)
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 (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 (TrajectoryStateOnSurface &propagated, const std::vector< l1extra::L1MuonParticle > &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...
	~L1MuonMatcherAlgo ()

Private Types
typedef StringCutObjectSelector < l1extra::L1MuonParticle >	L1Selector

Private Attributes
double	deltaPhi_
double	deltaR2_
	Matching cuts. More...
L1Selector	preselectionCut_
	Preselection cut to apply to L1 candidates before matching. More...
PropagateToMuon	prop_
bool	sortByDeltaPhi_
	Sort by deltaPhi instead of deltaR. More...

Detailed Description

Definition at line 29 of file L1MuonMatcherAlgo.h.

Member Typedef Documentation

using hltriggeroffline::L1MuonMatcherAlgo::ESTokens = PropagateToMuon::ESTokens

Definition at line 31 of file L1MuonMatcherAlgo.h.

typedef StringCutObjectSelector<l1extra::L1MuonParticle> hltriggeroffline::L1MuonMatcherAlgo::L1Selector

private

Definition at line 171 of file L1MuonMatcherAlgo.h.

Constructor & Destructor Documentation

L1MuonMatcherAlgo::L1MuonMatcherAlgo ( const edm::ParameterSet &  iConfig, const ESTokens &  iTokens  )

explicit

Definition at line 7 of file L1MuonMatcherAlgo.cc.

    : prop_(iConfig, tokens),
      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),
      sortByDeltaPhi_(iConfig.existsAs<bool>("sortByDeltaPhi") ? iConfig.getParameter<bool>("sortByDeltaPhi") : false) {
}

L1MuonMatcherAlgo::~L1MuonMatcherAlgo

(

)

Definition at line 16 of file L1MuonMatcherAlgo.cc.

{}

Member Function Documentation

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

References hltriggeroffline::PropagateToMuon::extrapolate(), and prop_.

Referenced by match(), and matchGeneric().

{ return prop_.extrapolate(tk); }

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

References hltriggeroffline::PropagateToMuon::extrapolate(), and prop_.

{ return prop_.extrapolate(tk); }

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

References hltriggeroffline::PropagateToMuon::extrapolate(), and prop_.

{ return prop_.extrapolate(state); }

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 18 of file L1MuonMatcherAlgo.cc.

References hltriggeroffline::PropagateToMuon::init(), and prop_.

Referenced by HLTMuonPlotter::beginRun().

{ prop_.init(iSetup); }

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

References srCondWrite_cfg::deltaPhi, deltaPhi_, reco::deltaR2(), deltaR2_, PV3DBase< T, PVType, FrameType >::eta(), reco::LeafCandidate::eta(), TrajectoryStateOnSurface::globalPosition(), 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());
    double thisDeltaR2 = ::deltaR2(double(pos.eta()), double(pos.phi()), l1.eta(), l1.phi());
    if ((fabs(thisDeltaPhi) < deltaPhi_) && (thisDeltaR2 < deltaR2_)) {
      deltaR = std::sqrt(thisDeltaR2);
      deltaPhi = thisDeltaPhi;
      return true;
    }
  }
  return false;
}

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

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

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

int hltriggeroffline::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 106 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	(	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 37 of file L1MuonMatcherAlgo.cc.

References matchGeneric(), and preselectionCut_.

                                                    {
  return matchGeneric(propagated, l1s, preselectionCut_, deltaR, deltaPhi);
  /*
      int match = -1;
      double minDeltaPhi = deltaPhi_;
      double minDeltaR2  = deltaR2_;
      GlobalPoint pos = propagated.globalPosition();
      for (int i = 0, n = l1s.size(); i < n; ++i) {
          const l1extra::L1MuonParticle &l1 = l1s[i];
          if (preselectionCut_(l1)) {
              double thisDeltaPhi = ::deltaPhi(double(pos.phi()),  l1.phi());
              double thisDeltaR2  = ::deltaR2(double(pos.eta()),
     double(pos.phi()), l1.eta(), l1.phi()); if ((fabs(thisDeltaPhi) <
     deltaPhi_) && (thisDeltaR2 < deltaR2_)) { // check both if (sortByDeltaPhi_
     ? (fabs(thisDeltaPhi) < fabs(minDeltaPhi)) : (thisDeltaR2 < minDeltaR2)) {
     // sort on one match = i; deltaR   = std::sqrt(thisDeltaR2); deltaPhi =
     thisDeltaPhi; if (sortByDeltaPhi_) minDeltaPhi = thisDeltaPhi; else
     minDeltaR2 = thisDeltaR2;
                  }
              }
          }
      }
      return match;
  */
}

template<typename Collection , typename Selector >

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

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

Referenced by match(), and matchGeneric().

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

template<typename Collection , typename Selector >

int hltriggeroffline::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 146 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 183 of file L1MuonMatcherAlgo.h.

References srCondWrite_cfg::deltaPhi, deltaPhi_, reco::deltaR2(), deltaR2_, PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, match(), dqmiodumpmetadata::n, getGTfromDQMFile::obj, PV3DBase< T, PVType, FrameType >::phi(), EgammaValidation_Wenu_cff::sel, sortByDeltaPhi_, and mathSSE::sqrt().

                                                             {
    typedef typename Collection::value_type obj;
    int match = -1;
    double minDeltaPhi = deltaPhi_;
    double minDeltaR2 = deltaR2_;
    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());
        double thisDeltaR2 = ::deltaR2(double(pos.eta()), double(pos.phi()), l1.eta(), l1.phi());
        if ((fabs(thisDeltaPhi) < deltaPhi_) && (thisDeltaR2 < deltaR2_)) {                               // check both
          if (sortByDeltaPhi_ ? (fabs(thisDeltaPhi) < fabs(minDeltaPhi)) : (thisDeltaR2 < minDeltaR2)) {  // sort on one
            match = i;
            deltaR = std::sqrt(thisDeltaR2);
            deltaPhi = thisDeltaPhi;
            if (sortByDeltaPhi_)
              minDeltaPhi = thisDeltaPhi;
            else
              minDeltaR2 = thisDeltaR2;
          }
        }
      }
    }
    return match;
  }

PropagateToMuon& hltriggeroffline::L1MuonMatcherAlgo::propagatorToMuon ( )

inline

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

Definition at line 53 of file L1MuonMatcherAlgo.h.

References prop_.

{ return prop_; }

const PropagateToMuon& hltriggeroffline::L1MuonMatcherAlgo::propagatorToMuon ( ) const

inline

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

Definition at line 55 of file L1MuonMatcherAlgo.h.

References prop_.

{ return prop_; }

Member Data Documentation

double hltriggeroffline::L1MuonMatcherAlgo::deltaPhi_

private

Definition at line 176 of file L1MuonMatcherAlgo.h.

Referenced by match(), and matchGeneric().

double hltriggeroffline::L1MuonMatcherAlgo::deltaR2_

private

Matching cuts.

Definition at line 176 of file L1MuonMatcherAlgo.h.

Referenced by match(), and matchGeneric().

L1Selector hltriggeroffline::L1MuonMatcherAlgo::preselectionCut_

private

Preselection cut to apply to L1 candidates before matching.

Definition at line 173 of file L1MuonMatcherAlgo.h.

Referenced by match().

PropagateToMuon hltriggeroffline::L1MuonMatcherAlgo::prop_

private

Definition at line 169 of file L1MuonMatcherAlgo.h.

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

bool hltriggeroffline::L1MuonMatcherAlgo::sortByDeltaPhi_

private

Sort by deltaPhi instead of deltaR.

Definition at line 179 of file L1MuonMatcherAlgo.h.

Referenced by matchGeneric().