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 */
  ~DTSLRecCluster() override {}

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

  DTSuperLayerId superLayerId() const { return theSlid; }

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

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

  AlgebraicMatrix projectionMatrix() const override { return theProjectionMatrix; }

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

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

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

  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 const AlgebraicMatrix theProjectionMatrix;

  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