ME0Segment.h

Go to the documentation of this file.

#ifndef GEMRecHit_ME0Segment_h
#define GEMRecHit_ME0Segment_h

#include "DataFormats/TrackingRecHit/interface/RecSegment.h"
#include "DataFormats/GEMRecHit/interface/ME0RecHitCollection.h"

#include <iosfwd>

class ME0DetId;

class ME0Segment final : public RecSegment {

 public:

  ME0Segment() : theChi2(0.){}
    
  ME0Segment(const std::vector<const ME0RecHit*>& proto_segment, const LocalPoint& origin, 
         const LocalVector& direction, const AlgebraicSymMatrix& errors, double chi2);

  ME0Segment(const std::vector<const ME0RecHit*>& proto_segment, const LocalPoint& origin, 
         const LocalVector& direction, const AlgebraicSymMatrix& errors, double chi2, double time, double timeErr);
  
  virtual ~ME0Segment();

  //--- Base class interface
  ME0Segment* clone() const { return new ME0Segment(*this); }

  LocalPoint localPosition() const { return theOrigin; }
  LocalError localPositionError() const ;
    
  LocalVector localDirection() const { return theLocalDirection; }
  LocalError localDirectionError() const ;

  AlgebraicVector parameters() const;

  AlgebraicSymMatrix parametersError() const { return theCovMatrix; }

  virtual AlgebraicMatrix projectionMatrix() const;

  virtual std::vector<const TrackingRecHit*> recHits() const;

  virtual std::vector<TrackingRecHit*> recHits();

  double chi2() const { return theChi2; };

  virtual int dimension() const { return 4; }

  virtual int degreesOfFreedom() const { return 2*nRecHits() - 4;}   

  //--- Extension of the interface
        
  const std::vector<ME0RecHit>& specificRecHits() const { return theME0RecHits; }

  int nRecHits() const { return theME0RecHits.size(); }        

  ME0DetId me0DetId() const { return  geographicalId(); }

  float time() const    { return theTimeValue; }
  float timeErr() const { return theTimeUncrt; }
    
  void print() const;       
    
 private:
    
  std::vector<ME0RecHit> theME0RecHits;
  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

};

std::ostream& operator<<(std::ostream& os, const ME0Segment& seg);

#endif