---

TransientTrackingRecHit.cc

Go to the documentation of this file.

#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHit.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/ErrorFrameTransformer.h"

#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"

const GeomDetUnit * TransientTrackingRecHit::detUnit() const
{
  return dynamic_cast<const GeomDetUnit*>(det());
}


GlobalPoint TransientTrackingRecHit::globalPosition() const {
  if(! hasGlobalPosition_){
    globalPosition_ = surface()->toGlobal(localPosition());
    hasGlobalPosition_ = true;
    return globalPosition_;
  }else{
    return globalPosition_;
  }
}

GlobalError TransientTrackingRecHit::globalPositionError() const {
  if(! hasGlobalError_){
    setPositionErrors();
    return globalError_;
  }else{
    return globalError_;
  }

}   

float 
TransientTrackingRecHit::errorGlobalR() const {
  if(!hasGlobalError_){
    setPositionErrors();
    return errorR_;
  }else{
    return errorR_;
  }
}

float 
TransientTrackingRecHit::errorGlobalZ() const {
  if(!hasGlobalError_){
    setPositionErrors();
    return errorZ_;
  }else{
    return errorZ_;
  }
}

float 
TransientTrackingRecHit::errorGlobalRPhi() const {
  if(!hasGlobalError_){
    setPositionErrors();
    return errorRPhi_;
  }else{
    return errorRPhi_;
  }
}

void
TransientTrackingRecHit::setPositionErrors() const {
  globalError_ = ErrorFrameTransformer().transform( localPositionError(), *surface() );
  errorRPhi_ = globalPosition().perp()*sqrt(globalError_.phierr(globalPosition())); 
  errorR_ = sqrt(globalError_.rerr(globalPosition()));
  errorZ_ = sqrt(globalError_.czz());
  hasGlobalError_ = true;
}

TransientTrackingRecHit::ConstRecHitContainer TransientTrackingRecHit::transientHits() const 
{
  // no components by default
  return ConstRecHitContainer();
}

TransientTrackingRecHit::RecHitPointer 
TransientTrackingRecHit::clone( const TrajectoryStateOnSurface& ts) const {
  return RecHitPointer(const_cast<TransientTrackingRecHit*>(this));
}

---