CSCSegAlgoDF.h

Go to the documentation of this file.

#ifndef CSCSegment_CSCSegAlgoDF_h
#define CSCSegment_CSCSegAlgoDF_h
 
#include <RecoLocalMuon/CSCSegment/src/CSCSegmentAlgorithm.h>
#include <DataFormats/CSCRecHit/interface/CSCRecHit2D.h>
 
#include <deque>
#include <vector>
 
class CSCSegAlgoPreClustering;
class CSCSegAlgoShowering;
class CSCSegFit;
 
class CSCSegAlgoDF : public CSCSegmentAlgorithm {
public:
 
  typedef std::vector<int> LayerIndex;
  typedef std::vector<const CSCRecHit2D*> ChamberHitContainer;
  typedef std::vector<const CSCRecHit2D*>::const_iterator ChamberHitContainerCIt;
  typedef std::deque<bool> BoolContainer;
 
  explicit CSCSegAlgoDF(const edm::ParameterSet& ps);
 
  ~CSCSegAlgoDF() override;
 
  std::vector<CSCSegment> buildSegments(const ChamberHitContainer& rechits);
 
  std::vector<CSCSegment> run(const CSCChamber* aChamber, const ChamberHitContainer& rechits) override;
 
private:
 
  void tryAddingHitsToSegment(const ChamberHitContainer& rechitsInChamber,
                              const ChamberHitContainerCIt i1,
                              const ChamberHitContainerCIt i2,
                              const LayerIndex& layerIndex);
 
  void flagHitsAsUsed(const ChamberHitContainer& rechitsInChamber);
 
  void pruneFromResidual(void);
 
  bool isHitNearSegment(const CSCRecHit2D* h) const;
  bool addHit(const CSCRecHit2D* hit, int layer);
  void updateParameters(void);
  bool hasHitOnLayer(int layer) const;
  void compareProtoSegment(const CSCRecHit2D* h, int layer);
  void dumpSegment(const CSCSegment& seg) const;
 
  // Member variables
  const std::string myName;
  const CSCChamber* theChamber;
  BoolContainer usedHits;
 
  ChamberHitContainer closeHits;
 
  ChamberHitContainer protoSegment;
  ChamberHitContainer secondSeedHits;
 
  // input from .cfi file
  bool debug;
  bool preClustering;
  int minHitsForPreClustering;
  //  bool   testSeg;
  bool Pruning;
  int minLayersApart;
  int nHitsPerClusterIsShower;
  //  float  nSigmaFromSegment;
  int minHitsPerSegment;
  //  int    muonsPerChamberMax;
  double dRPhiFineMax;
  double dPhiFineMax;
  float tanPhiMax;
  float tanThetaMax;
  float chi2Max;
  float maxRatioResidual;
 
  CSCSegAlgoPreClustering* preCluster_;
  CSCSegAlgoShowering* showering_;
  CSCSegFit* sfit_;
};
 
#endif