TSiStripRecHit2DLocalPos.h

Go to the documentation of this file.

#ifndef RECOTRACKER_TRANSIENTRACKINGRECHIT_TSiStripRecHit2DLocalPos_H
#define RECOTRACKER_TRANSIENTRACKINGRECHIT_TSiStripRecHit2DLocalPos_H

#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHit.h"
#include "RecoLocalTracker/ClusterParameterEstimator/interface/StripClusterParameterEstimator.h"
#include "TrackingTools/TransientTrackingRecHit/interface/HelpertRecHit2DLocalPos.h"
#include "DataFormats/Common/interface/RefGetter.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

class TSiStripRecHit2DLocalPos : public TransientTrackingRecHit{
public:

  typedef SiStripRecHit2D::ClusterRef SiStripClusterRef;
  
  typedef edm::LazyGetter<SiStripCluster>::value_ref  SiStripRegionalClusterRef;

  virtual ~TSiStripRecHit2DLocalPos() {}


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

  virtual AlgebraicVector parameters() const {return theHitData.parameters();}
  virtual AlgebraicSymMatrix parametersError() const {
    return HelpertRecHit2DLocalPos().parError( theHitData.localPositionError(), *det()); 
    //    return theHitData->parametersError();
  }
  
  virtual AlgebraicMatrix projectionMatrix() const {return theHitData.projectionMatrix();}
  virtual int dimension() const {return theHitData.dimension();}

  virtual LocalPoint localPosition() const {return theHitData.localPosition();}
  virtual LocalError localPositionError() const {return theHitData.localPositionError();}

  virtual const TrackingRecHit * hit() const {return &theHitData;};
  
  virtual std::vector<const TrackingRecHit*> recHits() const {
    return hit()->recHits();
  }
  virtual std::vector<TrackingRecHit*> recHits() {
    return theHitData.recHits();
  }

  virtual const GeomDetUnit* detUnit() const;

  virtual bool canImproveWithTrack() const {return true;}

  //RC virtual TSiStripRecHit2DLocalPos* clone(const TrajectoryStateOnSurface& ts) const;
  virtual RecHitPointer clone(const TrajectoryStateOnSurface& ts) const;

  // Extension of the TransientTrackingRecHit interface

  const SiStripRecHit2D* specificHit() const {return &theHitData;};
  const StripClusterParameterEstimator* cpe() const {return theCPE;}

  static RecHitPointer build( const GeomDet * geom, const SiStripRecHit2D* rh,
                  const StripClusterParameterEstimator* cpe,
                  float weight=1., float annealing=1.,
                  bool computeCoarseLocalPosition=false) {
    return RecHitPointer( new TSiStripRecHit2DLocalPos( geom, rh, cpe, weight, annealing,computeCoarseLocalPosition));
  }

  static RecHitPointer build( const LocalPoint& pos, const LocalError& err,
                  const GeomDet* det,
                  const SiStripClusterRef clust,
                  const StripClusterParameterEstimator* cpe,
                  float weight=1., float annealing=1.) {
    return RecHitPointer( new TSiStripRecHit2DLocalPos( pos, err, det, clust, cpe, weight, annealing));
  }

  static RecHitPointer build( const LocalPoint& pos, const LocalError& err,
                  const GeomDet* det,
                  const SiStripRegionalClusterRef clust,
                  const StripClusterParameterEstimator* cpe,
                  float weight=1., float annealing=1.) {
    return RecHitPointer( new TSiStripRecHit2DLocalPos( pos, err, det, clust, cpe, weight, annealing));
  }



private:

  SiStripRecHit2D              theHitData;
  const StripClusterParameterEstimator* theCPE;

  TSiStripRecHit2DLocalPos (const GeomDet * geom, const SiStripRecHit2D* rh,
                const StripClusterParameterEstimator* cpe,
                float weight, float annealing,
                bool computeCoarseLocalPosition) : 
    TransientTrackingRecHit(geom, weight, annealing), theCPE(cpe) 
    {
      if (rh->hasPositionAndError() || !computeCoarseLocalPosition)
    theHitData = SiStripRecHit2D(*rh);
      else{
    if (computeCoarseLocalPosition && !cpe){
      edm::LogError("TSiStripRecHit2DLocalPos")<<" trying to compute coarse local position but CPE is not provided. Not computing local position from disk for the transient tracking rechit.";
      theHitData = SiStripRecHit2D(*rh);
    }
    else{
      const GeomDetUnit* gdu = dynamic_cast<const GeomDetUnit*>(geom);
      LogDebug("TSiStripRecHit2DLocalPos")<<"calculating coarse position/error.";
      if (gdu){
        if (rh->cluster().isNonnull()){
          StripClusterParameterEstimator::LocalValues lval= theCPE->localParameters(*rh->cluster(), *gdu);
          theHitData = SiStripRecHit2D(lval.first, lval.second, geom->geographicalId(),rh->cluster());
        }else{
          StripClusterParameterEstimator::LocalValues lval= theCPE->localParameters(*rh->cluster_regional(), *gdu);
          theHitData = SiStripRecHit2D(lval.first, lval.second, geom->geographicalId(),rh->cluster_regional());
        }
      }else{
        edm::LogError("TSiStripRecHit2DLocalPos")<<" geomdet does not cast into geomdet unit. cannot create strip local parameters.";
        theHitData = SiStripRecHit2D(*rh);
      }
    }
      }
    }

  TSiStripRecHit2DLocalPos( const LocalPoint& pos, const LocalError& err,
                const GeomDet* det,
                const SiStripClusterRef clust,
                const StripClusterParameterEstimator* cpe,
                float weight, float annealing) :
    TransientTrackingRecHit(det, weight, annealing), theHitData(pos, err, det->geographicalId(), clust), 
    theCPE(cpe){} 

  //  TSiStripRecHit2DLocalPos( const TSiStripRecHit2DLocalPos& other ) :
  //  TransientTrackingRecHit( other.det()), 
  //  theHitData( other.specificHit()->clone()),
  //  theCPE( other.cpe()) {}

  TSiStripRecHit2DLocalPos( const LocalPoint& pos, const LocalError& err,
                const GeomDet* det,
                const SiStripRegionalClusterRef clust,              
                const StripClusterParameterEstimator* cpe,
                float weight, float annealing) :
    TransientTrackingRecHit(det, weight, annealing), theHitData(pos, err, det->geographicalId(), clust), 
    theCPE(cpe){} 
  
  

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

  virtual ConstRecHitContainer transientHits() const;

};

#endif