CMSSW_4_4_3_patch1/src/RecoTracker/TransientTrackingRecHit/interface/TSiStripMatchedRecHit.h

Go to the documentation of this file.

#ifndef RECOTRACKER_TRANSIENTRACKINGRECHIT_TSiStripMatchedRecHit_H
#define RECOTRACKER_TRANSIENTRACKINGRECHIT_TSiStripMatchedRecHit_H

#include "TrackingTools/TransientTrackingRecHit/interface/GenericTransientTrackingRecHit.h"
#include "TrackingTools/TransientTrackingRecHit/interface/HelpertRecHit2DLocalPos.h"
#include "RecoLocalTracker/SiStripRecHitConverter/interface/SiStripRecHitMatcher.h"
#include "RecoLocalTracker/ClusterParameterEstimator/interface/StripClusterParameterEstimator.h"
#include<memory>

#include "Geometry/TrackerGeometryBuilder/interface/GluedGeomDet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

class TSiStripMatchedRecHit : public GenericTransientTrackingRecHit{
public:

  virtual void getKfComponents( KfComponentsHolder & holder ) const {
      HelpertRecHit2DLocalPos().getKfComponents(holder, *hit(), *det()); 
  }

  virtual AlgebraicSymMatrix parametersError() const {
    return HelpertRecHit2DLocalPos().parError( localPositionError(), *det()); 
  }

  const GeomDetUnit* detUnit() const {return 0;}

  static RecHitPointer build( const GeomDet * geom, const TrackingRecHit * rh, 
                              const SiStripRecHitMatcher *matcher,
                              const StripClusterParameterEstimator* cpe=0,
                              float weight=1., float annealing=1.,
                              bool computeCoarseLocalPosition=false) {
    return RecHitPointer( new TSiStripMatchedRecHit( geom, rh, matcher,cpe, weight, annealing, computeCoarseLocalPosition));
  }

  static RecHitPointer build( const GeomDet * geom, std::auto_ptr<TrackingRecHit> rh, 
                              const SiStripRecHitMatcher *matcher,
                              const StripClusterParameterEstimator* cpe=0,
                              float weight=1., float annealing=1.,
                              bool computeCoarseLocalPosition=false) {
    return RecHitPointer( new TSiStripMatchedRecHit( geom, rh, matcher,cpe,weight, annealing, computeCoarseLocalPosition));
  }

  virtual RecHitPointer clone( const TrajectoryStateOnSurface& ts) const;
  virtual bool canImproveWithTrack() const {return (theMatcher != 0);}
  virtual ConstRecHitContainer  transientHits () const;

  struct DontCloneRecHit {};

  static void buildInPlace(std::auto_ptr<TSiStripMatchedRecHit> &memory,
                              const GeomDet * geom, const TrackingRecHit * rh,
                              const SiStripRecHitMatcher *matcher,
                              const StripClusterParameterEstimator* cpe=0,
                              float weight=1., float annealing=1.,
                              bool computeCoarseLocalPosition=false) {
        if (memory.get()) {
            memory->~TSiStripMatchedRecHit(); // call destructor
            new (memory.get()) TSiStripMatchedRecHit( geom, rh, matcher,cpe,weight, annealing, computeCoarseLocalPosition, DontCloneRecHit());
        } else {
            memory.reset(new TSiStripMatchedRecHit( geom, rh, matcher,cpe,weight, annealing, computeCoarseLocalPosition, DontCloneRecHit()));
        }
  }

  void clonePersistentHit()  { trackingRecHit_ = trackingRecHit_->clone(); }
  void clearPersistentHit() { trackingRecHit_ = 0; }

private:
  const SiStripRecHitMatcher* theMatcher; 
  const StripClusterParameterEstimator* theCPE;

private:
  TSiStripMatchedRecHit (const GeomDet * geom, const TrackingRecHit * rh, 
                         const SiStripRecHitMatcher *matcher,
                         const StripClusterParameterEstimator* cpe,
                         float weight, float annealing,
                         bool computeCoarseLocalPosition) : 
    GenericTransientTrackingRecHit(geom, *rh, weight, annealing), theMatcher(matcher),theCPE(cpe) {
    if (computeCoarseLocalPosition) ComputeCoarseLocalPosition();
  }

  TSiStripMatchedRecHit (const GeomDet * geom, std::auto_ptr<TrackingRecHit> rh,
                         const SiStripRecHitMatcher *matcher,
                         const StripClusterParameterEstimator* cpe,
                         float weight, float annealing,
                         bool computeCoarseLocalPosition) : 
    GenericTransientTrackingRecHit(geom, rh.release(), weight, annealing), theMatcher(matcher),theCPE(cpe) {
    if (computeCoarseLocalPosition) ComputeCoarseLocalPosition();
  }
  TSiStripMatchedRecHit (const GeomDet * geom, const TrackingRecHit * rh, 
                         const SiStripRecHitMatcher *matcher,
                         const StripClusterParameterEstimator* cpe,
                         float weight, float annealing,
                         bool computeCoarseLocalPosition,
                         const DontCloneRecHit &) : 
    GenericTransientTrackingRecHit(geom, const_cast<TrackingRecHit *>(rh), weight, annealing), theMatcher(matcher),theCPE(cpe) {
    if (computeCoarseLocalPosition) ComputeCoarseLocalPosition();
  }

private:
    void ComputeCoarseLocalPosition(){
      if (!theCPE || !theMatcher) return;
      const SiStripMatchedRecHit2D *orig = static_cast<const SiStripMatchedRecHit2D *> (trackingRecHit_);
      if (orig && !orig->hasPositionAndError()){
        LogDebug("TSiStripMatchedRecHit")<<"calculating coarse position/error.";
        const GeomDet *det = this->det();
        const GluedGeomDet *gdet = static_cast<const GluedGeomDet *> (det);
        LocalVector tkDir = det->surface().toLocal( det->position()-GlobalPoint(0,0,0));
        
        const SiStripMatchedRecHit2D* better=0;
        
        if(!orig->monoHit()->cluster().isNull()){
          const SiStripCluster& monoclust   = *orig->monoHit()->cluster();  
          const SiStripCluster& stereoclust = *orig->stereoHit()->cluster();
          
          StripClusterParameterEstimator::LocalValues lvMono = 
            theCPE->localParameters( monoclust, *gdet->monoDet());
          StripClusterParameterEstimator::LocalValues lvStereo = 
            theCPE->localParameters( stereoclust, *gdet->stereoDet());
          
          SiStripRecHit2D monoHit = SiStripRecHit2D( lvMono.first, lvMono.second,
                                                     gdet->monoDet()->geographicalId(),
                                                     orig->monoHit()->cluster());
          
          SiStripRecHit2D stereoHit = SiStripRecHit2D( lvStereo.first, lvStereo.second,
                                                       gdet->stereoDet()->geographicalId(),
                                                       orig->stereoHit()->cluster());
          better =  theMatcher->match(&monoHit,&stereoHit,gdet,tkDir);
        }else{
          const SiStripCluster& monoclust   = *orig->monoHit()->cluster_regional();  
          const SiStripCluster& stereoclust = *orig->stereoHit()->cluster_regional();
          StripClusterParameterEstimator::LocalValues lvMono = 
            theCPE->localParameters( monoclust, *gdet->monoDet());
          StripClusterParameterEstimator::LocalValues lvStereo = 
            theCPE->localParameters( stereoclust, *gdet->stereoDet());
          
          SiStripRecHit2D monoHit = SiStripRecHit2D( lvMono.first, lvMono.second,
                                                     gdet->monoDet()->geographicalId(),
                                                     orig->monoHit()->cluster_regional());
          
          SiStripRecHit2D stereoHit = SiStripRecHit2D( lvStereo.first, lvStereo.second,
                                                       gdet->stereoDet()->geographicalId(),
                                                       orig->stereoHit()->cluster_regional());
          better =  theMatcher->match(&monoHit,&stereoHit,gdet,tkDir);
          
        }
        if (!better) {
          edm::LogWarning("TSiStripMatchedRecHit")<<"could not get a matching rechit.";
        }else{
          trackingRecHit_ = better->clone();
        }
      }
    }

  virtual TSiStripMatchedRecHit* clone() const {
    return new TSiStripMatchedRecHit(*this);
  }

};



#endif