DTRecSegment4D.h

Go to the documentation of this file.

#ifndef DTRecHit_DTRecSegment4D_h
#define DTRecHit_DTRecSegment4D_h

/* Base Class Headers */
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"

/* Collaborating Class Declarations */
#include "DataFormats/DTRecHit/interface/DTSLRecSegment2D.h"
#include "DataFormats/DTRecHit/interface/DTChamberRecSegment2D.h"

/* C++ Headers */
#include <iosfwd>

class DTRecSegment4D : public RecSegment {

 public:
  friend class DTSegmentUpdator;
  DTRecSegment4D() : theProjection(none), theDimension(0) {}
  
  DTRecSegment4D(const DTChamberRecSegment2D& phiSeg, const DTSLRecSegment2D& zedSeg, const LocalPoint& posZInCh, const LocalVector& dirZInCh);

  DTRecSegment4D(const DTChamberRecSegment2D& phiSeg);

  DTRecSegment4D(const DTSLRecSegment2D& zedSeg, const LocalPoint& posZInCh, const LocalVector& dirZInCh);

  ~DTRecSegment4D() override ;

  //--- Base class interface

  DTRecSegment4D* clone() const override { return new DTRecSegment4D(*this);}

  AlgebraicVector parameters() const override ;

  AlgebraicSymMatrix parametersError() const override ;

  AlgebraicMatrix projectionMatrix() const override;

  LocalPoint localPosition() const override { return thePosition;}

  LocalError localPositionError() const override ;

  LocalVector localDirection() const override { return theDirection; }

  LocalError localDirectionError() const override ;

  // Chi2 of the segment fit
  double chi2() const override ;
  
  // Degrees of freedom of the segment fit
  int degreesOfFreedom() const override ;

  // Dimension (in parameter space)
  int dimension() const override { return theDimension; }

  // Access to component RecHits (if any)
  std::vector<const TrackingRecHit*> recHits() const override ;

  // Non-const access to component RecHits (if any)
  std::vector<TrackingRecHit*> recHits() override ;


  //--- Extension of the interface

  
  bool hasPhi() const {return (theProjection==full || theProjection==phi);}
  
  bool hasZed() const {return (theProjection==full || theProjection==Z);}
  
  const DTChamberRecSegment2D *phiSegment() const {
    return hasPhi()? &thePhiSeg: nullptr;
  }
    
  const DTSLRecSegment2D *zSegment() const {
    return hasZed()? &theZedSeg : nullptr;
  }
    
  void setPosition(LocalPoint pos) { thePosition = pos; }

  void setDirection(LocalVector dir) { theDirection = dir; }

  void setCovMatrix(const AlgebraicSymMatrix& mat) { theCovMatrix = mat; }

  virtual DTChamberId chamberId() const;
    
 private:
  enum Projection {full, phi, Z, none};
  Projection theProjection;

  DTChamberRecSegment2D *phiSegment() {return &thePhiSeg;}
    
  DTSLRecSegment2D *zSegment() {return &theZedSeg;}

  LocalPoint thePosition;   // in chamber frame
  LocalVector theDirection; // in chamber frame

  void setCovMatrixForZed(const LocalPoint& posZInCh);
    
  // the covariance matrix, has the following meaning
  // mat[0][0]=sigma (dx/dz)
  // mat[1][1]=sigma (dy/dz)
  // mat[2][2]=sigma (x)
  // mat[3][3]=sigma (y)
  // mat[0][2]=cov(dx/dz,x)
  // mat[1][3]=cov(dy/dz,y)
  AlgebraicSymMatrix theCovMatrix; 

  DTChamberRecSegment2D thePhiSeg;
  DTSLRecSegment2D theZedSeg;

  int theDimension; // the dimension of this rechit

};

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

#endif // DTRecHit_DTRecSegment4D_h