DTRecSegment2D.h

Go to the documentation of this file.

#ifndef TrackingRecHit_DTRecSegment2D_h
#define TrackingRecHit_DTRecSegment2D_h

/* Base Class Headers */
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"

/* Collaborating Class Declarations */
#include "DataFormats/GeometrySurface/interface/LocalError.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"

#include "DataFormats/DTRecHit/interface/DTRecHit1D.h"
/* C++ Headers */
#include <iosfwd>

/* Fwd declaration */
class DTSegmentUpdator;

/* ====================================================================== */

/* Class DTRecSegment2D Interface */

class DTRecSegment2D : public RecSegment{

 public:

  DTRecSegment2D(): theChi2(0.0), theT0(0.), theVdrift(0.) {}
  
  DTRecSegment2D(DetId id, const std::vector<DTRecHit1D>& hits) ;
  
  DTRecSegment2D(DetId id, 
         LocalPoint &position, LocalVector &direction,
         AlgebraicSymMatrix & covMatrix, double chi2, 
         std::vector<DTRecHit1D> &hits1D);

  ~DTRecSegment2D() override;

  /* Operations */ 

  DTRecSegment2D* clone() const override { return new DTRecSegment2D(*this);}


  AlgebraicVector parameters() const override {
    return param( localPosition(), localDirection());
  }

  // The parameter error matrix 
  AlgebraicSymMatrix parametersError() const override;

  AlgebraicMatrix projectionMatrix() const override {
    return theProjectionMatrix;
  }
    
  int dimension() const override { return 2;}
    
  LocalPoint localPosition() const override {return thePosition; }
  
  LocalError localPositionError() const override ;
  
  LocalVector localDirection() const override { return theDirection; }

  LocalError localDirectionError() const override;

  double chi2() const override { return theChi2; }
    
  int degreesOfFreedom() const override ;

  // Access to component RecHits (if any)
  std::vector<const TrackingRecHit*> recHits() const override ;

  // Non-const access to component RecHits (if any)
  std::vector<TrackingRecHit*> recHits() override ;

  std::vector<DTRecHit1D> specificRecHits() const ;
  
  AlgebraicSymMatrix covMatrix() const {return theCovMatrix;}

  double t0() const {return theT0;}
  bool ist0Valid() const {return (theT0 > -998.) ? true : false;}

  double vDrift() const {return theVdrift;}

 protected:
  friend class DTSegmentUpdator;
  void setPosition(const LocalPoint& pos);
  void setDirection(const LocalVector& dir);
  void setCovMatrix(const AlgebraicSymMatrix& cov);
  void setChi2(const double& chi2);
  void update(std::vector<DTRecHit1D> & updatedRecHits);
  void setT0(const double& t0);
  void setVdrift(const double& vdrift);

  LocalPoint  thePosition;  // in SL frame
  LocalVector theDirection; // in SL frame
  
  AlgebraicSymMatrix theCovMatrix; // the covariance matrix

  double theChi2;           // chi2 of the fit
  double theT0;             // T0 as coming from the fit
  double theVdrift;             // vDrift as coming from the fit

  std::vector<DTRecHit1D> theHits; // the hits with defined R/L
  

 private:

  static const AlgebraicMatrix theProjectionMatrix;
  
  AlgebraicVector param( const LocalPoint& lp, const LocalVector& lv) const {
    AlgebraicVector result(2);
    result[1]=lp.x();
    result[0]=lv.x()/lv.z();
    return result;
  }

};
std::ostream& operator<<(std::ostream& os, const DTRecSegment2D& seg);
#endif // TrackingRecHit_DTRecSegment2D_h