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, LocalPoint origin, 
           LocalVector direction, AlgebraicSymMatrix errors, double chi2);

   ME0Segment(const std::vector<const ME0RecHit*>& proto_segment, LocalPoint origin, 
          LocalVector direction, 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