GEMSegment.h

Go to the documentation of this file.

#ifndef GEMRecHit_GEMSegment_h
#define GEMRecHit_GEMSegment_h
 
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"
#include "DataFormats/GEMRecHit/interface/GEMRecHitCollection.h"
 
#include <iosfwd>
 
class GEMDetId;
 
class GEMSegment final : public RecSegment {
public:
  GEMSegment() : theChi2(0.) {}
 
  GEMSegment(const std::vector<const GEMRecHit*>& proto_segment,
             const LocalPoint& origin,
             const LocalVector& direction,
             const AlgebraicSymMatrix& errors,
             double chi2);
 
  GEMSegment(const std::vector<const GEMRecHit*>& proto_segment,
             const LocalPoint& origin,
             const LocalVector& direction,
             const AlgebraicSymMatrix& errors,
             double chi2,
             float bx);
 
  GEMSegment(const std::vector<const GEMRecHit*>& proto_segment,
             const LocalPoint& origin,
             const LocalVector& direction,
             const AlgebraicSymMatrix& errors,
             double chi2,
             double time,
             double timeErr);
 
  ~GEMSegment() override;
 
  //--- Base class interface
  GEMSegment* clone() const override { return new GEMSegment(*this); }
 
  LocalPoint localPosition() const override { return theOrigin; }
  LocalError localPositionError() const override;
 
  LocalVector localDirection() const override { return theLocalDirection; }
  LocalError localDirectionError() const override;
 
  AlgebraicVector parameters() const override;
 
  AlgebraicSymMatrix parametersError() const override { return theCovMatrix; }
 
  AlgebraicMatrix projectionMatrix() const override;
 
  std::vector<const TrackingRecHit*> recHits() const override;
 
  std::vector<TrackingRecHit*> recHits() override;
 
  double chi2() const override { return theChi2; };
 
  int dimension() const override { return 4; }
 
  int degreesOfFreedom() const override { return 2 * nRecHits() - 4; }
 
  //--- Extension of the interface
 
  const std::vector<GEMRecHit>& specificRecHits() const { return theGEMRecHits; }
 
  int nRecHits() const { return theGEMRecHits.size(); }
 
  GEMDetId gemDetId() const { return geographicalId(); }
 
  float time() const { return theTimeValue; }
  float timeErr() const { return theTimeUncrt; }
  float bunchX() const { return theBX; }
  void print() const;
 
private:
  std::vector<GEMRecHit> theGEMRecHits;
  LocalPoint theOrigin;             // in chamber frame - the GeomDet local coordinate system
  LocalVector theLocalDirection;    // in chamber frame - the GeomDet local coordinate system
  AlgebraicSymMatrix theCovMatrix;  // the covariance matrix
  double theChi2;                   // the Chi squared of the segment fit
  double theTimeValue;              // the best time estimate of the segment
  double theTimeUncrt;              // the uncertainty on the time estimation
  float theBX;                      // the bunch crossing
};
 
std::ostream& operator<<(std::ostream& os, const GEMSegment& seg);
 
#endif