AlignmentMuonSelector.h

Go to the documentation of this file.

#ifndef Alignment_CommonAlignmentAlgorithm_AlignmentMuonSelector_h
#define Alignment_CommonAlignmentAlgorithm_AlignmentMuonSelector_h

#include "DataFormats/MuonReco/interface/Muon.h"
#include "CommonTools/RecoAlgos/interface/MuonSelector.h"
#include <vector>

namespace edm { class Event; }

class AlignmentMuonSelector 
{

 public:

  typedef std::vector<const reco::Muon*> Muons; 

  AlignmentMuonSelector(const edm::ParameterSet & cfg);

  ~AlignmentMuonSelector();

  Muons select(const Muons& muons, const edm::Event& evt) const;

 private:

  Muons basicCuts(const Muons& muons) const;

  Muons theNHighestPtMuons(const Muons& muons) const;
  
  Muons theBestMassPairCombinationMuons(const Muons& muons) const;
  
  struct ComparePt {
    bool operator()( const reco::Muon* t1, const reco::Muon* t2 ) const {
      return t1->pt()> t2->pt();
    }
  };
  ComparePt ptComparator;

  bool applyBasicCuts,applyNHighestPt,applyMultiplicityFilter,applyMassPairFilter;
  int nHighestPt,minMultiplicity;
  double pMin,pMax,ptMin,ptMax,etaMin,etaMax,phiMin,phiMax;
  double nHitMinSA,nHitMaxSA,chi2nMaxSA;
  double nHitMinGB,nHitMaxGB,chi2nMaxGB;
  double nHitMinTO,nHitMaxTO,chi2nMaxTO;
  double minMassPair,maxMassPair;

};

#endif