MeasurementDet.h

Go to the documentation of this file.

#ifndef MeasurementDet_H
#define MeasurementDet_H

#include "TrackingTools/MeasurementDet/interface/TempMeasurements.h"

#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "DataFormats/TrackerRecHit2D/interface/BaseTrackerRecHit.h"
#include "DataFormats/TrackingRecHit/interface/InvalidTrackingRecHit.h"

class TrajectoryStateOnSurface;
class Propagator;
class MeasurementEstimator;
class MeasurementTrackerEvent;

class MeasurementDet {
public:
  typedef tracking::TempMeasurements TempMeasurements;
  typedef TrackingRecHit::ConstRecHitContainer RecHitContainer;

  using SimpleHitContainer = std::vector<BaseTrackerRecHit*>;

  MeasurementDet(const GeomDet* gdet)
      : theGeomDet(gdet),
        theMissingHit(std::make_shared<InvalidTrackingRecHit>(fastGeomDet(), TrackingRecHit::missing)),
        theInactiveHit(std::make_shared<InvalidTrackingRecHit>(fastGeomDet(), TrackingRecHit::inactive)) {}

  virtual ~MeasurementDet() = default;
  virtual RecHitContainer recHits(const TrajectoryStateOnSurface&, const MeasurementTrackerEvent&) const = 0;

  // use a MeasurementEstimator to filter the hits (same algo as below..)
  // default as above
  virtual bool recHits(const TrajectoryStateOnSurface& stateOnThisDet,
                       const MeasurementEstimator&,
                       const MeasurementTrackerEvent& data,
                       RecHitContainer& result,
                       std::vector<float>&) const {
    result = recHits(stateOnThisDet, data);
    return !result.empty();
  }

  // default for non-tracker dets...
  virtual bool recHits(SimpleHitContainer& result,
                       const TrajectoryStateOnSurface& stateOnThisDet,
                       const MeasurementEstimator&,
                       const MeasurementTrackerEvent& data) const {
    return false;
  }

  std::vector<TrajectoryMeasurement> fastMeasurements(const TrajectoryStateOnSurface& stateOnThisDet,
                                                      const TrajectoryStateOnSurface&,
                                                      const Propagator&,
                                                      const MeasurementEstimator& est,
                                                      const MeasurementTrackerEvent& data) const {
    TempMeasurements tmps;
    measurements(stateOnThisDet, est, data, tmps);
    std::vector<TrajectoryMeasurement> result;
    result.reserve(tmps.size());
    int index[tmps.size()];
    tmps.sortIndex(index);
    for (std::size_t i = 0; i != tmps.size(); ++i) {
      auto j = index[i];
      result.emplace_back(stateOnThisDet, std::move(tmps.hits[j]), tmps.distances[j]);
    }
    return result;
  }

  // return false if missing ( if inactive is true and one hit)
  virtual bool measurements(const TrajectoryStateOnSurface& stateOnThisDet,
                            const MeasurementEstimator& est,
                            const MeasurementTrackerEvent& data,
                            TempMeasurements& result) const = 0;

  const GeomDet& fastGeomDet() const { return *theGeomDet; }
  virtual const GeomDet& geomDet() const { return *theGeomDet; }

  const Surface& surface() const { return geomDet().surface(); }

  const Surface::PositionType& position() const { return geomDet().position(); }

  virtual bool isActive(const MeasurementTrackerEvent& data) const = 0;
  virtual bool hasBadComponents(const TrajectoryStateOnSurface& tsos, const MeasurementTrackerEvent& data) const = 0;

private:
  const GeomDet* theGeomDet;

protected:
  TrackingRecHit::ConstRecHitPointer theMissingHit;
  TrackingRecHit::ConstRecHitPointer theInactiveHit;
};

#endif