/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/RecoTracker/TransientTrackingRecHit/interface/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