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

  virtual ~DTRecSegment2D();

  /* Operations */ 

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


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

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

  virtual AlgebraicMatrix projectionMatrix() const {
    if ( !isInitialized) initialize();
    return theProjectionMatrix;
  }
    
  virtual int dimension() const { return 2;}
    
  virtual LocalPoint localPosition() const {return thePosition; }
  
  virtual LocalError localPositionError() const ;
  
  virtual LocalVector localDirection() const { return theDirection; }

  virtual LocalError localDirectionError() const;

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

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

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

  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 bool isInitialized;
  static AlgebraicMatrix theProjectionMatrix;
  
  void initialize() const {
    isInitialized=true;
    theProjectionMatrix = AlgebraicMatrix( 2, 5, 0);
    theProjectionMatrix[0][1]=1;
    theProjectionMatrix[1][3]=1;
  }
  
  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