GEMCSCSegment.h

Go to the documentation of this file.

#ifndef GEMRecHit_GEMCSCSegment_h
#define GEMRecHit_GEMCSCSegment_h

#include "DataFormats/TrackingRecHit/interface/RecSegment.h"
     
#include "DataFormats/GEMRecHit/interface/GEMRecHitCollection.h"
#include "DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h"
#include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
     
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iosfwd>

class CSCDetId;

class GEMCSCSegment final : public RecSegment {

public:

    GEMCSCSegment() : theChi2(0.) {}
    
    GEMCSCSegment(const CSCSegment* csc_segment, const std::vector<const GEMRecHit*> gem_rhs, LocalPoint origin, LocalVector direction, AlgebraicSymMatrix errors, double chi2);
      
    ~GEMCSCSegment() override;

    //--- Base class interface
    GEMCSCSegment* clone() const override { return new GEMCSCSegment(*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;

    double chi2() const override { return theChi2; };

    int dimension() const override { return 4; }

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

    int nRecHits() const 
    { 
      return (theGEMRecHits.size() + theCSCSegment.specificRecHits().size()); 
      edm::LogVerbatim("GEMCSCSegment")<< "[GEMCSCSegment :: nRecHits] CSC RecHits: " <<theCSCSegment.specificRecHits().size()
                       << " + GEM RecHits: " << theGEMRecHits.size() << " = " <<(theGEMRecHits.size() + theCSCSegment.specificRecHits().size());
    }
    

    //--- Return the constituents in different ways
    const CSCSegment                        cscSegment() const { return theCSCSegment; }
    const std::vector<GEMRecHit>&           gemRecHits() const { return theGEMRecHits; }
    const std::vector<CSCRecHit2D>&         cscRecHits() const { return theCSCSegment.specificRecHits(); }
    std::vector<const TrackingRecHit*> recHits() const override;
    std::vector<TrackingRecHit*>       recHits() override;

    CSCDetId cscDetId() const { return  geographicalId(); }
    
    void print() const;     
    
 private:
    std::vector<GEMRecHit> theGEMRecHits;      // store GEM Rechits
    CSCSegment theCSCSegment;                  // store CSC RecHits and store CSC Segment 
                                               // eventually we have to disentangle if later on we decide 
                                               // not to have a one-to-one relationship anymore
                                               // i.e. if we allow the GEMCSC segment to modify the
                                               // (selection of the rechits of the) CSC segment
    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;

};

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

#endif