TransientTrackingRecHit.h

Go to the documentation of this file.

#ifndef TransientTrackingRecHit_H
#define TransientTrackingRecHit_H

#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include <TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h>
#include "DataFormats/GeometrySurface/interface/Surface.h" 
//RC #include "DataFormats/Common/interface/OwnVector.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "DataFormats/GeometrySurface/interface/BlockWipedAllocator.h"

class GeomDetUnit;

class TransientTrackingRecHit : public TrackingRecHit, 
                public ReferenceCountedInEvent {
public:

  //RC typedef edm::OwnVector<const TransientTrackingRecHit>        RecHitContainer;

  typedef ReferenceCountingPointer< TransientTrackingRecHit>        RecHitPointer;
  typedef ConstReferenceCountingPointer< TransientTrackingRecHit>   ConstRecHitPointer;
  typedef std::vector<ConstRecHitPointer>                           RecHitContainer;
  typedef std::vector<ConstRecHitPointer>                           ConstRecHitContainer;

  explicit TransientTrackingRecHit(const GeomDet * geom=0, float weight=1., float annealing=1.) : 
    TrackingRecHit(geom ? geom->geographicalId().rawId() : 0), 
    geom_(geom), weight_(weight), annealing_(annealing),
    globalPosition_(0,0,0),
    globalError_(GlobalError::NullMatrix()),errorR_(0),errorZ_(0),errorRPhi_(0),
    hasGlobalPosition_(false), hasGlobalError_(false){}

  explicit TransientTrackingRecHit(const GeomDet * geom, DetId id, Type type=valid, float weight=1., float annealing=1. ) : 
    TrackingRecHit(id, type), 
    geom_(geom), weight_(weight), annealing_(annealing),
    globalPosition_(0,0,0),
    globalError_(GlobalError::NullMatrix()),errorR_(0),errorZ_(0),errorRPhi_(0),
    hasGlobalPosition_(false),hasGlobalError_(false){}

  explicit TransientTrackingRecHit(const GeomDet * geom, TrackingRecHit::id_type id, Type type=valid, float weight=1., float annealing=1. ) : 
    TrackingRecHit(id, type),
    geom_(geom),  weight_(weight), annealing_(annealing),
    globalPosition_(0,0,0),
    globalError_(GlobalError::NullMatrix()),errorR_(0),errorZ_(0),errorRPhi_(0),
    hasGlobalPosition_(false),hasGlobalError_(false){}
  
  explicit TransientTrackingRecHit(const GeomDet * geom, TrackingRecHit const & rh, float weight=1., float annealing=1. ) : 
    TrackingRecHit(rh.geographicalId(), rh.type()),
    geom_(geom), weight_(weight), annealing_(annealing),
    globalPosition_(0,0,0),
    globalError_(GlobalError::NullMatrix()),errorR_(0),errorZ_(0),errorRPhi_(0),
    hasGlobalPosition_(false),hasGlobalError_(false){}



  // Extension of the TrackingRecHit interface

  const GeomDet * det() const {return geom_;}
  virtual const Surface * surface() const {return &(geom_->surface());}

  virtual const GeomDetUnit * detUnit() const;

  virtual GlobalPoint globalPosition() const ;
  virtual GlobalError globalPositionError() const ;

  float errorGlobalR() const;
  float errorGlobalZ() const;
  float errorGlobalRPhi() const;

  virtual RecHitPointer clone (const TrajectoryStateOnSurface& ts) const;

  virtual bool canImproveWithTrack() const {return false;}

  virtual const TrackingRecHit * hit() const = 0;
  
  virtual ConstRecHitContainer transientHits() const;

  void setWeight(float weight){weight_ = weight;}

  float weight() const {return weight_;}
  

  void setAnnealingFactor(float annealing) {annealing_ = annealing;} 

  float getAnnealingFactor() const {return annealing_;} 

  virtual float clusterProbability() const { return 1; }

private:
  void setPositionErrors() const;
  
  const GeomDet * geom_ ;

  float weight_;
  float annealing_;

  // caching of some variable for fast access
  mutable GlobalPoint globalPosition_;  
  mutable GlobalError globalError_;
  mutable float errorR_,errorZ_,errorRPhi_;
  mutable bool hasGlobalPosition_;
  mutable bool hasGlobalError_;
 
  // hide the clone method for ReferenceCounted. Warning: this method is still 
  // accessible via the bas class TrackingRecHit interface!
  virtual TransientTrackingRecHit * clone() const = 0;

};

#endif