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