ClusterShapeTrackFilter Class Reference

#include <ClusterShapeTrackFilter.h>

Inheritance diagram for ClusterShapeTrackFilter:

Public Member Functions
	ClusterShapeTrackFilter (const edm::ParameterSet &ps, edm::ConsumesCollector &iC)
virtual bool	operator() (const reco::Track *, const std::vector< const TrackingRecHit * > &hits, const TrackerTopology *tTopo) const
void	update (const edm::Event &ev, const edm::EventSetup &es) override
virtual	~ClusterShapeTrackFilter ()
Public Member Functions inherited from PixelTrackFilter
virtual bool	operator() (const reco::Track *) const
virtual bool	operator() (const reco::Track *, const Hits &) const
virtual bool	operator() (const reco::Track *, const Hits &, const TrackerTopology *tTopo) const
virtual	~PixelTrackFilter ()

Private Member Functions
float	areaParallelogram (const Global2DVector &a, const Global2DVector &b) const
std::vector< GlobalVector >	getGlobalDirs (const std::vector< GlobalPoint > &globalPoss) const
std::vector< GlobalPoint >	getGlobalPoss (const std::vector< const TrackingRecHit * > &recHits) const

Private Attributes
double	ptMax
double	ptMin
const SiPixelClusterShapeCache *	theClusterShapeCache
edm::EDGetTokenT < SiPixelClusterShapeCache >	theClusterShapeCacheToken
const ClusterShapeHitFilter *	theFilter
const TrackerGeometry *	theTracker

Additional Inherited Members
Public Types inherited from PixelTrackFilter
typedef std::vector< const TrackingRecHit * >	Hits

Detailed Description

Definition at line 27 of file ClusterShapeTrackFilter.h.

Constructor & Destructor Documentation

ClusterShapeTrackFilter::ClusterShapeTrackFilter	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector &	iC
	)

Definition at line 33 of file ClusterShapeTrackFilter.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), ptMax, and ptMin.

                                                                                                   :
  theClusterShapeCacheToken(iC.consumes<SiPixelClusterShapeCache>(ps.getParameter<edm::InputTag>("clusterShapeCacheSrc"))),
  theTracker(nullptr),
  theFilter(nullptr),
  theClusterShapeCache(nullptr)
{
  // Get ptMin if available
  ptMin = (ps.exists("ptMin") ? ps.getParameter<double>("ptMin") : 0.);
  ptMax = (ps.exists("ptMax") ? ps.getParameter<double>("ptMax") : 999999.);
}

ClusterShapeTrackFilter::~ClusterShapeTrackFilter ( )

virtual

Definition at line 45 of file ClusterShapeTrackFilter.cc.

{
}

Member Function Documentation

float ClusterShapeTrackFilter::areaParallelogram ( const Global2DVector &  a, const Global2DVector &  b  ) const

private

Definition at line 68 of file ClusterShapeTrackFilter.cc.

References PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

{  
  return a.x() * b.y() - a.y() * b.x();
}

vector< GlobalVector > ClusterShapeTrackFilter::getGlobalDirs ( const std::vector< GlobalPoint > &  globalPoss ) const

private

Definition at line 75 of file ClusterShapeTrackFilter.cc.

References EnergyCorrector::c, CircleFromThreePoints::center(), PixelRecoUtilities::curvature(), CircleFromThreePoints::curvature(), dir, mag2(), AlCaHLTBitMon_ParallelJobs::p, slope, funct::sqr(), mathSSE::sqrt(), findQualityFiles::v, PV2DBase< T, PVType, FrameType >::x(), Basic2DVector< T >::x(), PV2DBase< T, PVType, FrameType >::y(), and Basic2DVector< T >::y().

{
  // Get 2d points
  vector<Global2DVector> p;
  for(vector<GlobalPoint>::const_iterator ig = g.begin();
                                          ig!= g.end(); ig++)
     p.push_back( Global2DVector(ig->x(), ig->y()) );

  //
  vector<GlobalVector> globalDirs;

  // Determine circle
  CircleFromThreePoints circle(g[0],g[1],g[2]);

  if(circle.curvature() != 0.)
  {
    Global2DVector c (circle.center().x(), circle.center().y());

    float rad2 = (p[0] - c).mag2();
    float a12 = asin(fabsf(areaParallelogram(p[0] - c, p[1] - c)) / rad2);

    float slope = (g[1].z() - g[0].z()) / a12;

    float cotTheta = slope * circle.curvature(); // == sinhEta
    float coshEta  = sqrt(1 + sqr(cotTheta));    // == 1/sinTheta

    // Calculate globalDirs
    float sinTheta =       1. / coshEta;
    float cosTheta = cotTheta * sinTheta;

    int dir;
    if(areaParallelogram(p[0] - c, p[1] - c) > 0) dir = 1; else dir = -1;

    float curvature = circle.curvature();

    for(vector<Global2DVector>::const_iterator ip = p.begin();
                                               ip!= p.end(); ip++)
    {
      Global2DVector v = (*ip - c)*curvature*dir;
      globalDirs.push_back(GlobalVector(-v.y()*sinTheta,
                                         v.x()*sinTheta,
                                               cosTheta));
    }
  }

  return globalDirs;
}

vector< GlobalPoint > ClusterShapeTrackFilter::getGlobalPoss ( const std::vector< const TrackingRecHit * > &  recHits ) const

private

Definition at line 125 of file ClusterShapeTrackFilter.cc.

{
  vector<GlobalPoint> globalPoss;

  for(vector<const TrackingRecHit *>::const_iterator recHit = recHits.begin();
                                                     recHit!= recHits.end();
                                                     recHit++)
  {
    DetId detId = (*recHit)->geographicalId();

    GlobalPoint gpos = 
      theTracker->idToDet(detId)->toGlobal((*recHit)->localPosition());

    globalPoss.push_back(gpos);
  }

  return globalPoss;
}

bool ClusterShapeTrackFilter::operator() ( const reco::Track *  track, const std::vector< const TrackingRecHit * > &  hits, const TrackerTopology *  tTopo  ) const

virtual

Definition at line 146 of file ClusterShapeTrackFilter.cc.

References HitInfo::getInfo(), i, TrackingRecHit::isValid(), LogTrace, convertSQLiteXML::ok, and PtMinSelector_cfg::ptMin.

Referenced by HIPixelTrackFilter::operator()().

{
  // Do not even look at pairs
  if(recHits.size() <= 2) return true;

  // Check pt
  if(track->pt() < ptMin ||
     track->pt() > ptMax)
  {
    LogTrace("ClusterShapeTrackFilter")
       << "  [ClusterShapeTrackFilter] pt not in range: "
       << ptMin << " " << track->pt() << " " << ptMax;
    return false;
  }

  // Get global positions
  vector<GlobalPoint>  globalPoss = getGlobalPoss(recHits);

  // Get global directions
  vector<GlobalVector> globalDirs = getGlobalDirs(globalPoss);

  bool ok = true;

  // Check whether shape of pixel cluster is compatible
  // with local track direction
  for(unsigned int i = 0; i < recHits.size(); i++)
  {
    const SiPixelRecHit* pixelRecHit =
      dynamic_cast<const SiPixelRecHit *>(recHits[i]);

    if(!pixelRecHit->isValid())
    { 
      ok = false; break; 
    }

    if(! theFilter->isCompatible(*pixelRecHit, globalDirs[i], *theClusterShapeCache) )
    {
      LogTrace("ClusterShapeTrackFilter")
         << "  [ClusterShapeTrackFilter] clusShape problem"
         << HitInfo::getInfo(*recHits[i],tTopo);

      ok = false; break;
    }
  }

  return ok;
}

void ClusterShapeTrackFilter::update ( const edm::Event &  ev, const edm::EventSetup &  es  )

overridevirtual

Implements PixelTrackFilter.

Definition at line 50 of file ClusterShapeTrackFilter.cc.

References edm::EventSetup::get(), edm::Event::getByToken(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), theClusterShapeCache, theClusterShapeCacheToken, theFilter, theTracker, and patCandidatesForDimuonsSequences_cff::tracker.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), HIPixelTrackFilter::update(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                                                {
  edm::Handle<SiPixelClusterShapeCache> cache;
  ev.getByToken(theClusterShapeCacheToken, cache);
  theClusterShapeCache = cache.product();

  // Get tracker geometry
  edm::ESHandle<TrackerGeometry> tracker;
  es.get<TrackerDigiGeometryRecord>().get(tracker);
  theTracker = tracker.product();

  // Get cluster shape hit filter
  edm::ESHandle<ClusterShapeHitFilter> shape;
  es.get<CkfComponentsRecord>().get("ClusterShapeHitFilter",shape);
  theFilter = shape.product();
}

Member Data Documentation

double ClusterShapeTrackFilter::ptMax

private

Definition at line 53 of file ClusterShapeTrackFilter.h.

Referenced by ClusterShapeTrackFilter().

double ClusterShapeTrackFilter::ptMin

private

Definition at line 52 of file ClusterShapeTrackFilter.h.

Referenced by ClusterShapeTrackFilter().

const SiPixelClusterShapeCache* ClusterShapeTrackFilter::theClusterShapeCache

private

Definition at line 50 of file ClusterShapeTrackFilter.h.

Referenced by update().

edm::EDGetTokenT<SiPixelClusterShapeCache> ClusterShapeTrackFilter::theClusterShapeCacheToken

private

Definition at line 46 of file ClusterShapeTrackFilter.h.

Referenced by update().

const ClusterShapeHitFilter* ClusterShapeTrackFilter::theFilter

private

Definition at line 49 of file ClusterShapeTrackFilter.h.

Referenced by update().

const TrackerGeometry* ClusterShapeTrackFilter::theTracker

private

Definition at line 48 of file ClusterShapeTrackFilter.h.

Referenced by update().