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"
#include "RecoTracker/MeasurementDet/interface/MeasurementTrackerEvent.h"
#include "RecoTracker/MeasurementDet/src/TkMeasurementDetSet.h"

class TrackingRecHit;
class LocalTrajectoryParameters;

class dso_hidden TkPixelMeasurementDet final : 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,
             PxMeasurementConditionSet & conditionSet );

  void update(PxMeasurementDetSet &data, const detset & detSet ) { 
    data.update(index(), detSet);
    data.setActiveThisEvent(index(), true);
  }

  void setEmpty(PxMeasurementDetSet & data) { data.setEmpty(index());  }
  bool isEmpty(const PxMeasurementDetSet & data) const {return data.empty(index());}

  virtual ~TkPixelMeasurementDet() { }

  // all hits
  virtual RecHitContainer recHits( const TrajectoryStateOnSurface&, const MeasurementTrackerEvent & dat) const override;

  // only hits compatible with tsos 
  RecHitContainer compHits( const TrajectoryStateOnSurface&, const MeasurementTrackerEvent & dat, float xl, float yl ) const;



 // simple hits
  virtual bool recHits(SimpleHitContainer & result,  
               const TrajectoryStateOnSurface& stateOnThisDet, const MeasurementEstimator&, const MeasurementTrackerEvent & data) const {
    assert("not implemented for Pixel yet"==nullptr);
  }

 

  virtual bool measurements( const TrajectoryStateOnSurface& stateOnThisDet,
                const MeasurementEstimator& est, const MeasurementTrackerEvent & dat,
                TempMeasurements & result) const;


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

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

  void setActive(bool active) { conditionSet().setActive(index(), active); }
  void setActiveThisEvent(PxMeasurementDetSet & data, bool active) const { data.setActiveThisEvent(index(), active); }
  bool isActive(const MeasurementTrackerEvent & data) const { return data.pixelData().isActive(index()); }

  bool hasBadComponents( const TrajectoryStateOnSurface &tsos, const MeasurementTrackerEvent & dat ) const ; 

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

  int index() const { return index_; }
  void setIndex(int i) { index_ = i; }

private:
  unsigned int id_;
  std::vector< LocalPoint > badRocPositions_;

  int index_;
  PxMeasurementConditionSet * theDetConditions;
  PxMeasurementConditionSet & conditionSet() { return *theDetConditions; }
  const PxMeasurementConditionSet & conditionSet() const { return *theDetConditions; }

  const PixelClusterParameterEstimator * cpe() const { return conditionSet().pixelCPE(); }

 public:

  inline bool accept(SiPixelClusterRefNew & r, const std::vector<bool> skipClusters) const {
    
    if(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()];
  }

};

#endif