TkPixelMeasurementDet.h

Go to the documentation of this file.

#ifndef TkPixelMeasurementDet_H
#define TkPixelMeasurementDet_H

#include "TrackingTools/MeasurementDet/interface/MeasurementDet.h"
#include "RecoLocalTracker/ClusterParameterEstimator/interface/PixelClusterParameterEstimator.h"
//#include "DataFormats/SiPixelCluster/interface/SiPixelClusterFwd.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"

class TransientTrackingRecHit;
class LocalTrajectoryParameters;

class TkPixelMeasurementDet : public MeasurementDet {
public:

  typedef edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> SiPixelClusterRef;
  
  typedef edmNew::DetSet<SiPixelCluster> detset;
  typedef detset::const_iterator const_iterator;
  typedef PixelClusterParameterEstimator::LocalValues    LocalValues;

  TkPixelMeasurementDet( const GeomDet* gdet,
             const PixelClusterParameterEstimator* cpe);

  void update( const detset & detSet, 
           const edm::Handle<edmNew::DetSetVector<SiPixelCluster> > & h,
           unsigned int id ) { 
    detSet_ = detSet; 
    handle_ = h;
    id_ = id;
    empty = false;
    activeThisEvent_ = true;
  }

  void setEmpty(){empty = true; activeThisEvent_ = true; }

  virtual ~TkPixelMeasurementDet() { }

  virtual RecHitContainer recHits( const TrajectoryStateOnSurface& ) const;

  virtual std::vector<TrajectoryMeasurement> 
  fastMeasurements( const TrajectoryStateOnSurface& stateOnThisDet, 
            const TrajectoryStateOnSurface& startingState, 
            const Propagator&, 
            const MeasurementEstimator&) const;

  const PixelGeomDetUnit& specificGeomDet() const {return static_cast<PixelGeomDetUnit const &>(fastGeomDet());}

  TransientTrackingRecHit::RecHitPointer 
  buildRecHit( const SiPixelClusterRef & cluster,
           const LocalTrajectoryParameters & ltp) const;

  void setActive(bool active) { activeThisPeriod_ = active; activeThisEvent_ = true; if (!active) empty = true; }
  void setActiveThisEvent(bool active) { activeThisEvent_ = active;  if (!active) empty = true; }
  bool isActive() const { return activeThisEvent_ && activeThisPeriod_; }

  bool hasBadComponents( const TrajectoryStateOnSurface &tsos ) const ; 

  void setBadRocPositions(std::vector< LocalPoint > & positions) { badRocPositions_.swap(positions); }
  void clearBadRocPositions() { badRocPositions_.clear(); }
private:

  const PixelClusterParameterEstimator* theCPE;
  edm::Handle<edmNew::DetSetVector<SiPixelCluster> > handle_;
  detset detSet_;
  std::vector< LocalPoint > badRocPositions_;
  const std::vector<bool>* skipClusters_;
  unsigned int id_;
  bool empty;
  bool activeThisEvent_, activeThisPeriod_;
  
 public:

  inline bool accept(SiPixelClusterRefNew & r) const {
    
    if(0==skipClusters_ || skipClusters_->empty()) return true;
    if (r.key()>=skipClusters_->size()){
      edm::LogError("IndexMisMatch")<<r.key()<<" is larger than: "<<skipClusters_->size()<<" no skipping done";
      return true;
    }
    return not (*skipClusters_)[r.key()];
  }

  void setClusterToSkip(const std::vector<bool>*skip ){
     skipClusters_ = skip;
  }
};

#endif