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/PatternTools/interface/MeasurementEstimator.h>
#include <TrackingTools/PatternTools/interface/TrajectoryMeasurement.h>


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


class TkGluedMeasurementDet : public MeasurementDet {
public:

  TkGluedMeasurementDet( const GluedGeomDet* gdet,const SiStripRecHitMatcher* matcher, const StripClusterParameterEstimator* cpe);
  void init(const MeasurementDet* monoDet,
        const MeasurementDet* stereoDet);

  virtual RecHitContainer recHits( const TrajectoryStateOnSurface&) const;

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

  virtual std::vector<TrajectoryMeasurement> 
  fastMeasurements( const TrajectoryStateOnSurface& stateOnThisDet, 
            const TrajectoryStateOnSurface& startingState, 
            const Propagator&, 
            const MeasurementEstimator&) const;

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

  bool isActive() const {return monoDet()->isActive() && stereoDet()->isActive(); }
         
  bool hasBadComponents( const TrajectoryStateOnSurface &tsos ) const {
    return (monoDet()->hasBadComponents(tsos) || stereoDet()->hasBadComponents(tsos));}

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


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

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

  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_.push_back(
            TSiStripMatchedRecHit::build( geomDet_, std::auto_ptr<TrackingRecHit>(hit.clone()), 
                              matcher_,cpe_)
            );
      hasNewHits_ = true; 
    }
    void addProjected(const TransientTrackingRecHit& hit,
              const GlobalVector & gdir) ;
    SiStripRecHitMatcher::Collector & collector() { return collector_; }
    bool hasNewMatchedHits() const { return hasNewHits_;  }
    void clearNewMatchedHitsFlag() { hasNewHits_ = false; }
  private: 
    const GeomDet              * geomDet_;
    const SiStripRecHitMatcher * matcher_;
    const StripClusterParameterEstimator* cpe_;
    RecHitContainer       & 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,
                    std::vector<TrajectoryMeasurement> & target) ;
    void add(SiStripMatchedRecHit2D const& hit) ;
    void addProjected(const TransientTrackingRecHit& hit,
              const GlobalVector & gdir) ;
    
    SiStripRecHitMatcher::Collector & collector() { return collector_; }
    bool hasNewMatchedHits() const { return hasNewHits_;  }
    void clearNewMatchedHitsFlag() { hasNewHits_ = false; }
  private: 
    const GeomDet              * geomDet_;
    const SiStripRecHitMatcher * matcher_;
    const StripClusterParameterEstimator* cpe_;
    const TrajectoryStateOnSurface & stateOnThisDet_;
    const MeasurementEstimator     & est_;
    std::vector<TrajectoryMeasurement> & target_;
    SiStripRecHitMatcher::Collector collector_;       
    bool hasNewHits_;
  };
  
  
  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 TransientTrackingRecHit& 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