DTSLRecCluster.h

Go to the documentation of this file.

#ifndef DTSLRECCLUSTER_H
#define DTSLRECCLUSTER_H

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

/* Collaborating Class Declarations */
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"
#include "DataFormats/DTRecHit/interface/DTRecHit1DPair.h"

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

/* ====================================================================== */

/* Class DTSLRecCluster Interface */

class DTSLRecCluster : public RecHit1D {

  public:

    /* Constructor */ 
    DTSLRecCluster() {}

    DTSLRecCluster(const DTSuperLayerId id, const std::vector<DTRecHit1DPair>& pair) ;

    DTSLRecCluster(const DTSuperLayerId id,
                   const LocalPoint&,
                   const LocalError&,
                   const std::vector<DTRecHit1DPair>& pair) ;

    /* Destructor */ 
    virtual ~DTSLRecCluster() {}

    /* Operations */ 
    virtual DTSLRecCluster* clone() const { return new DTSLRecCluster(*this); }

    DTSuperLayerId superLayerId() const { return theSlid; }

    virtual AlgebraicVector parameters() const {
      return param( localPosition());
    }

    // The parameter error matrix 
    virtual AlgebraicSymMatrix parametersError() const {
      return parError( localPositionError());
    }

    virtual AlgebraicMatrix projectionMatrix() const {
      if ( !isInitialized) initialize();
      return theProjectionMatrix;
    }

    virtual int dimension() const { return 2;}
    virtual LocalPoint localPosition() const { return thePos; }
    virtual LocalError localPositionError() const { return thePosError; }

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

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

    std::vector<DTRecHit1DPair> specificRecHits() const { return thePairs; }

    int nHits() const { return thePairs.size() ; }

  private:
    DTSuperLayerId theSlid;

    LocalPoint thePos;
    LocalError thePosError;

    std::vector<DTRecHit1DPair> thePairs;

  private:
    static bool isInitialized;
    static AlgebraicMatrix theProjectionMatrix;

    void initialize() const {
      isInitialized=true;
      theProjectionMatrix = AlgebraicMatrix( 2, 5, 0);
      theProjectionMatrix[0][1]=1;
    }

    AlgebraicVector param( const LocalPoint& lp) const {
      AlgebraicVector result(1);
      result[1]=lp.x();
      return result;
    }

    AlgebraicSymMatrix parError( const LocalError& le) const {
      AlgebraicSymMatrix m(1);
      m[0][0]=le.xx();
      return m;
    }

  protected:

};
std::ostream& operator<<(std::ostream& os, const DTSLRecCluster& seg);
#endif // DTSLRECCLUSTER_H