TkGluedMeasurementDet.h

Go to the documentation of this file.

#ifndef TkGluedMeasurementDet_H
#define TkGluedMeasurementDet_H

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

class GluedGeomDet;
//class SiStripRecHitMatcher;
#include "RecoLocalTracker/SiStripRecHitConverter/interface/SiStripRecHitMatcher.h"
#include "RecoLocalTracker/ClusterParameterEstimator/interface/StripClusterParameterEstimator.h"
#include "RecoTracker/TransientTrackingRecHit/interface/TSiStripMatchedRecHit.h"
#include <TrackingTools/DetLayers/interface/MeasurementEstimator.h>
#include <TrackingTools/PatternTools/interface/TrajectoryMeasurement.h>

#include "FWCore/Utilities/interface/Visibility.h"

class dso_hidden TkGluedMeasurementDet final : public MeasurementDet {
public:
  TkGluedMeasurementDet(const GluedGeomDet* gdet,
                        const SiStripRecHitMatcher* matcher,
                        const StripClusterParameterEstimator* cpe);
  void init(const MeasurementDet* monoDet, const MeasurementDet* stereoDet, const TrackerTopology* tTopo);

  RecHitContainer recHits(const TrajectoryStateOnSurface&, const MeasurementTrackerEvent& data) const override;

  // simple hits
  bool recHits(SimpleHitContainer& result,
               const TrajectoryStateOnSurface& stateOnThisDet,
               const MeasurementEstimator&,
               const MeasurementTrackerEvent& data) const override;

  const GluedGeomDet& specificGeomDet() const { return static_cast<GluedGeomDet const&>(fastGeomDet()); }

  bool measurements(const TrajectoryStateOnSurface& stateOnThisDet,
                    const MeasurementEstimator& est,
                    const MeasurementTrackerEvent& data,
                    TempMeasurements& result) const override;

  const TkStripMeasurementDet* monoDet() const { return theMonoDet; }
  const TkStripMeasurementDet* stereoDet() const { return theStereoDet; }

  unsigned int rawId() const { return fastGeomDet().geographicalId(); }

  bool isActive(const MeasurementTrackerEvent& data) const override {
    return monoDet()->isActive(data) && stereoDet()->isActive(data);
  }

  bool hasBadComponents(const TrajectoryStateOnSurface& tsos, const MeasurementTrackerEvent& data) const override {
    return (monoDet()->hasBadComponents(tsos, data) || stereoDet()->hasBadComponents(tsos, data));
  }

private:
  const SiStripRecHitMatcher* theMatcher;
  const StripClusterParameterEstimator* theCPE;
  const TkStripMeasurementDet* theMonoDet;
  const TkStripMeasurementDet* theStereoDet;
  const TrackerTopology* theTopology;

  template <typename Collector>
  void doubleMatch(const TrajectoryStateOnSurface& ts,
                   const MeasurementTrackerEvent& data,
                   Collector& collector) const dso_internal;

  template <typename Collector>
  void collectRecHits(const TrajectoryStateOnSurface&,
                      const MeasurementTrackerEvent& data,
                      Collector& coll) const dso_internal;

  // for TTRH
  class dso_internal HitCollectorForRecHits {
  public:
    typedef SiStripRecHitMatcher::Collector Collector;
    HitCollectorForRecHits(const GeomDet* geomDet,
                           const SiStripRecHitMatcher* matcher,
                           const StripClusterParameterEstimator* cpe,
                           RecHitContainer& target);
    void add(SiStripMatchedRecHit2D const& hit) {
      target_.emplace_back(hit.cloneSH());
      hasNewHits_ = true;
    }
    void addProjected(const TrackingRecHit& hit, const GlobalVector& gdir);
    SiStripRecHitMatcher::Collector& collector() { return collector_; }
    bool hasNewMatchedHits() const { return hasNewHits_; }
    void clearNewMatchedHitsFlag() { hasNewHits_ = false; }
    static bool filter() { return false; }  
    size_t size() const { return target_.size(); }

    static const MeasurementEstimator& estimator() {
      static const MeasurementEstimator* dummy = nullptr;
      return *dummy;
    }

  private:
    const GeomDet* geomDet_;
    const SiStripRecHitMatcher* matcher_;
    const StripClusterParameterEstimator* cpe_;
    RecHitContainer& target_;
    SiStripRecHitMatcher::Collector collector_;
    bool hasNewHits_;
  };

  // for TRH
  class dso_internal HitCollectorForSimpleHits {
  public:
    typedef SiStripRecHitMatcher::Collector Collector;
    HitCollectorForSimpleHits(const GeomDet* geomDet,
                              const SiStripRecHitMatcher* matcher,
                              const StripClusterParameterEstimator* cpe,
                              const TrajectoryStateOnSurface& stateOnThisDet,
                              const MeasurementEstimator& est,
                              SimpleHitContainer& target);
    void add(SiStripMatchedRecHit2D const& hit);
    void addProjected(const TrackingRecHit& hit, const GlobalVector& gdir);
    SiStripRecHitMatcher::Collector& collector() { return collector_; }
    bool hasNewMatchedHits() const { return hasNewHits_; }
    void clearNewMatchedHitsFlag() { hasNewHits_ = false; }
    bool filter() const {
      return matcher_->preFilter();
    }  // if true mono-colection will been filter using the estimator before matching
    size_t size() const { return target_.size(); }
    const MeasurementEstimator& estimator() { return est_; }

  private:
    const GeomDet* geomDet_;
    const SiStripRecHitMatcher* matcher_;
    const StripClusterParameterEstimator* cpe_;
    const TrajectoryStateOnSurface& stateOnThisDet_;
    const MeasurementEstimator& est_;
    SimpleHitContainer& target_;
    SiStripRecHitMatcher::Collector collector_;
    bool hasNewHits_;
  };

  class dso_internal HitCollectorForFastMeasurements {
  public:
    typedef TransientTrackingRecHit::RecHitPointer RecHitPointer;
    typedef SiStripRecHitMatcher::Collector Collector;

    HitCollectorForFastMeasurements(const GeomDet* geomDet,
                                    const SiStripRecHitMatcher* matcher,
                                    const StripClusterParameterEstimator* cpe,
                                    const TrajectoryStateOnSurface& stateOnThisDet,
                                    const MeasurementEstimator& est,
                                    TempMeasurements& target);
    void add(SiStripMatchedRecHit2D const& hit);
    void addProjected(const TrackingRecHit& hit, const GlobalVector& gdir);

    SiStripRecHitMatcher::Collector& collector() { return collector_; }
    bool hasNewMatchedHits() const { return hasNewHits_; }
    void clearNewMatchedHitsFlag() { hasNewHits_ = false; }
    bool filter() const {
      return matcher_->preFilter();
    }  // if true mono-colection will been filter using the estimator before matching
    size_t size() const { return target_.size(); }
    const MeasurementEstimator& estimator() { return est_; }

  private:
    const GeomDet* geomDet_;
    const SiStripRecHitMatcher* matcher_;
    const StripClusterParameterEstimator* cpe_;
    const TrajectoryStateOnSurface& stateOnThisDet_;
    const MeasurementEstimator& est_;
    TempMeasurements& target_;
    SiStripRecHitMatcher::Collector collector_;
    bool hasNewHits_;
  };

  RecHitContainer projectOnGluedDet(const std::vector<SiStripRecHit2D>& hits,
                                    const TrajectoryStateOnSurface& ts) const dso_internal;
  template <typename HitCollector>
  void projectOnGluedDet(HitCollector& collector,
                         const std::vector<SiStripRecHit2D>& hits,
                         const GlobalVector& gdir) const dso_internal;

  RecHitContainer projectOnGluedDet(const RecHitContainer& hits, const TrajectoryStateOnSurface& ts) const dso_internal;
  template <typename HitCollector>
  void projectOnGluedDet(HitCollector& collector,
                         const RecHitContainer& hits,
                         const GlobalVector& gdir) const dso_internal;

  void checkProjection(const TrajectoryStateOnSurface& ts,
                       const RecHitContainer& monoHits,
                       const RecHitContainer& stereoHits) const;
  void checkHitProjection(const TrackingRecHit& hit,
                          const TrajectoryStateOnSurface& ts,
                          const GeomDet& det) const dso_internal;

  bool testStrips(const TrajectoryStateOnSurface& tsos,
                  const BoundPlane& gluedPlane,
                  const TkStripMeasurementDet& mdet) const dso_internal;
};

#endif