TValidTrackingRecHit.h

Go to the documentation of this file.

#ifndef TValidTrackingRecHit_H
#define TValidTrackingRecHit_H

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

class GeomDetUnit;


class TValidTrackingRecHit : public TransientTrackingRecHit {
public:
  
  typedef ReferenceCountingPointer<TransientTrackingRecHit>        RecHitPointer;
  typedef ConstReferenceCountingPointer<TransientTrackingRecHit>   ConstRecHitPointer;
  typedef std::vector<ConstRecHitPointer>                           RecHitContainer;
  typedef std::vector<ConstRecHitPointer>                           ConstRecHitContainer;

  TValidTrackingRecHit(const GeomDet * geom) : 
    TransientTrackingRecHit(geom->geographicalId()), geom_(geom) {}


  template<typename... Args>
  TValidTrackingRecHit(const GeomDet * geom, Args && ...args) : 
    TransientTrackingRecHit(std::forward<Args>(args)...), geom_(geom) {}

  // to be moved in children
  TrackingRecHit * cloneHit() const { return hit()->clone();}

  // Extension of the TrackingRecHit interface
  virtual const GeomDet * det() const GCC11_FINAL {return geom_;}
  virtual const Surface * surface() const GCC11_FINAL {return &(det()->surface());}


  virtual GlobalPoint globalPosition() const GCC11_FINAL {
      return surface()->toGlobal(localPosition());
  }
  
  GlobalError globalPositionError() const GCC11_FINAL { return ErrorFrameTransformer().transform( localPositionError(), *surface() );}
  float errorGlobalR() const GCC11_FINAL { return std::sqrt(globalPositionError().rerr(globalPosition()));}
  float errorGlobalZ() const GCC11_FINAL { return std::sqrt(globalPositionError().czz()); }
  float errorGlobalRPhi() const GCC11_FINAL { return globalPosition().perp()*sqrt(globalPositionError().phierr(globalPosition())); }

  // once cache removed will obsolete the above
  TrackingRecHitGlobalState globalState() const {
    GlobalError  
      globalError = ErrorFrameTransformer::transform( localPositionError(), *surface() );
    auto gp = globalPosition();
    float r = gp.perp();
    float errorRPhi = r*std::sqrt(float(globalError.phierr(gp))); 
    float errorR = std::sqrt(float(globalError.rerr(gp)));
    float errorZ = std::sqrt(float(globalError.czz()));
    return (TrackingRecHitGlobalState){
      gp.basicVector(), r, gp.barePhi(),
    errorR,errorZ,errorRPhi
    };
  }


  virtual bool canImproveWithTrack() const {return false;}

 
  virtual float clusterProbability() const { return 1.f; }

private:
  
  const GeomDet * geom_ ;

 
  // hide the clone method for ReferenceCounted. Warning: this method is still 
  // accessible via the bas class TrackingRecHit interface!
  virtual TValidTrackingRecHit * clone() const = 0;

};

#endif