CSCSegment.h

Go to the documentation of this file.

#ifndef CSCRecHit_CSCSegment_h
#define CSCRecHit_CSCSegment_h

#include <DataFormats/TrackingRecHit/interface/RecSegment.h>
#include <DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h>

#include <iosfwd>

class CSCDetId;

class CSCSegment final : public RecSegment {

public:

    CSCSegment() : theChi2(0.), aME11a_duplicate(false) {}
    
    CSCSegment(const std::vector<const CSCRecHit2D*>& proto_segment, LocalPoint origin, 
            LocalVector direction, const AlgebraicSymMatrix& errors, double chi2);
  
    ~CSCSegment() override;

    //--- Base class interface
    CSCSegment* clone() const override { return new CSCSegment(*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<CSCRecHit2D>& specificRecHits() const { return theCSCRecHits; }

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

    CSCDetId cscDetId() const { return  geographicalId(); }

    void setDuplicateSegments(std::vector<CSCSegment*>& duplicates);

    bool isME11a_duplicate() const { return (!theDuplicateSegments.empty() ? true : false); }
    // a copy of the duplicated segments (ME1/1a only) 
    const std::vector< CSCSegment> & duplicateSegments() const { return theDuplicateSegments; } 
    
    bool testSharesAllInSpecificRecHits( const std::vector<CSCRecHit2D>& specificRecHits_1,
                     const std::vector<CSCRecHit2D>& specificRecHits_2,
                     CSCRecHit2D::SharedInputType) const;
    
    //bool sharesRecHits(CSCSegment  & anotherSegment, CSCRecHit2D::SharedInputType);
    // checks if ALL the rechits share the specific input (allWires, allStrips or all)
    bool sharesRecHits(const CSCSegment  & anotherSegment, CSCRecHit2D::SharedInputType sharesInput) const;
    // checks if ALL the rechits share SOME wire AND SOME strip input
    bool sharesRecHits(const CSCSegment  & anotherSegment) const;

    float time() const;
    
    void print() const;     
    
 private:
    
    std::vector<CSCRecHit2D> theCSCRecHits;
    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;
    bool aME11a_duplicate;
    std::vector<CSCSegment> theDuplicateSegments;// ME1/1a only
};

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

#endif // CSCRecHit_CSCSegment_h