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 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);

  virtual ~CSCSegAlgoDF();

  std::vector<CSCSegment> buildSegments(const ChamberHitContainer& rechits);

  std::vector<CSCSegment> run(const CSCChamber* aChamber, const ChamberHitContainer& rechits); 

private:


  void tryAddingHitsToSegment( const ChamberHitContainer& rechitsInChamber,
                               const ChamberHitContainerCIt i1, 
                               const ChamberHitContainerCIt i2,
                               const LayerIndex& layerIndex);

  void flagHitsAsUsed(const ChamberHitContainer& rechitsInChamber);
 
  void pruneFromResidual();


  void orderSecondSeed( GlobalPoint gp1, 
                                  const ChamberHitContainerCIt i1,
                                  const ChamberHitContainerCIt i2,
                              const ChamberHitContainer& rechits,
                                  const LayerIndex& layerIndex );

    
  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);
  CLHEP::HepMatrix derivativeMatrix(void) const;
  AlgebraicSymMatrix weightMatrix(void) const;
  AlgebraicSymMatrix calculateError(void) const;
  void flipErrors(AlgebraicSymMatrix&) const;

  // Member variables
  const std::string myName; 
  const CSCChamber* theChamber;
  BoolContainer usedHits;

  ChamberHitContainer closeHits;

  ChamberHitContainer protoSegment;
  ChamberHitContainer secondSeedHits;
  float       protoSlope_u;
  float       protoSlope_v;
  LocalPoint  protoIntercept;       
  double      protoChi2;
  LocalVector protoDirection;

  // 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_;

};

#endif