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 GCC11_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