MuonMesh.h

Go to the documentation of this file.

#ifndef _MuonIdentification_MuonMesh_h_
#define _MuonIdentification_MuonMesh_h_
//
// Creates a mesh of muons connected by overlapping segments
// Original author: Lindsey Gray
//

#include <vector>
#include <utility>
#include <map>
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"

class CSCGeometry;

class MuonMesh {

  typedef std::map<reco::Muon*,
                   std::vector<std::pair<reco::Muon*,
                                         std::pair<reco::MuonChamberMatch*,
                                                   reco::MuonSegmentMatch*
                                                  >
                                        >
                              >
                  > MeshType;

  typedef std::vector<std::pair<reco::Muon*,
                                std::pair<reco::MuonChamberMatch*,
                                          reco::MuonSegmentMatch*
                                         >
                               >
                     > AssociationType;

 public:
  
  MuonMesh(const edm::ParameterSet&);
  
  void runMesh(std::vector<reco::Muon>* p) {fillMesh(p); pruneMesh();}

  void clearMesh() { mesh_.clear(); }

  void setCSCGeometry(const CSCGeometry* pg) { geometry_ = pg; } 

  bool isDuplicateOf(const CSCSegmentRef& lhs, const CSCSegmentRef& rhs) const;
  bool isDuplicateOf(const std::pair<CSCDetId,CSCSegmentRef>& rhs, 
             const std::pair<CSCDetId,CSCSegmentRef>& lhs) const;
  bool isClusteredWith(const std::pair<CSCDetId,CSCSegmentRef>& lhs, 
               const std::pair<CSCDetId,CSCSegmentRef>& rhs) const;

 private:

  void fillMesh(std::vector<reco::Muon>*);

  void pruneMesh();

  
  // implement to remove cases where two segments in the same
  // chamber overlap within 2 sigma of ALL of their errors
  bool withinTwoSigma(const std::pair<CSCDetId,CSCSegmentRef>& rhs, 
              const std::pair<CSCDetId,CSCSegmentRef>& lhs) const { return false; }
  
  
  
  

  MeshType mesh_;

  // geometry cache for segment arbitration
   const CSCGeometry* geometry_;
   
   // do various cleanings?
   const bool doME1a, doOverlaps, doClustering;
   // overlap and clustering parameters
   const double OverlapDPhi, OverlapDTheta, ClusterDPhi, ClusterDTheta;
};

#endif