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