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