MuonArbitrationMethods.h

Go to the documentation of this file.

#ifndef MuonIdentification_MuonArbitrationMethods_h
#define MuonIdentification_MuonArbitrationMethods_h
 
#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "DataFormats/MuonReco/interface/MuonChamberMatch.h"
 
#include <utility>
 
// Author: Jake Ribnik (UCSB)
 
struct SortMuonSegmentMatches {
  SortMuonSegmentMatches(unsigned int flag) { flag_ = flag; }
  bool operator()(std::pair<reco::MuonChamberMatch*, reco::MuonSegmentMatch*> p1,
                  std::pair<reco::MuonChamberMatch*, reco::MuonSegmentMatch*> p2) {
    reco::MuonChamberMatch* cm1 = p1.first;
    reco::MuonSegmentMatch* sm1 = p1.second;
    reco::MuonChamberMatch* cm2 = p2.first;
    reco::MuonSegmentMatch* sm2 = p2.second;
 
    if (flag_ == reco::MuonSegmentMatch::BestInChamberByDX || flag_ == reco::MuonSegmentMatch::BestInStationByDX ||
        flag_ == reco::MuonSegmentMatch::BelongsToTrackByDX)
      return fabs(sm1->x - cm1->x) < fabs(sm2->x - cm2->x);
    if (flag_ == reco::MuonSegmentMatch::BestInChamberByDR || flag_ == reco::MuonSegmentMatch::BestInStationByDR ||
        flag_ == reco::MuonSegmentMatch::BelongsToTrackByDR) {
      if ((!sm1->hasZed()) || (!sm2->hasZed()))  // no y information so return dx
        return fabs(sm1->x - cm1->x) < fabs(sm2->x - cm2->x);
      return sqrt(pow(sm1->x - cm1->x, 2) + pow(sm1->y - cm1->y, 2)) <
             sqrt(pow(sm2->x - cm2->x, 2) + pow(sm2->y - cm2->y, 2));
    }
    if (flag_ == reco::MuonSegmentMatch::BestInChamberByDXSlope ||
        flag_ == reco::MuonSegmentMatch::BestInStationByDXSlope ||
        flag_ == reco::MuonSegmentMatch::BelongsToTrackByDXSlope)
      return fabs(sm1->dXdZ - cm1->dXdZ) < fabs(sm2->dXdZ - cm2->dXdZ);
    if (flag_ == reco::MuonSegmentMatch::BestInChamberByDRSlope ||
        flag_ == reco::MuonSegmentMatch::BestInStationByDRSlope ||
        flag_ == reco::MuonSegmentMatch::BelongsToTrackByDRSlope) {
      if ((!sm1->hasZed()) || (!sm2->hasZed()))  // no y information so return dx
        return fabs(sm1->dXdZ - cm1->dXdZ) < fabs(sm2->dXdZ - cm2->dXdZ);
      return sqrt(pow(sm1->dXdZ - cm1->dXdZ, 2) + pow(sm1->dYdZ - cm1->dYdZ, 2)) <
             sqrt(pow(sm2->dXdZ - cm2->dXdZ, 2) + pow(sm2->dYdZ - cm2->dYdZ, 2));
    }
 
    return false;  // is this appropriate? fix this
  }
 
  unsigned int flag_;
};
 
#endif