ClusterShapeTrajectoryFilter Class Reference

#include <ClusterShapeTrajectoryFilter.h>

Inheritance diagram for ClusterShapeTrajectoryFilter:

Public Member Functions
	ClusterShapeTrajectoryFilter (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
std::string	name () const override
bool	qualityFilter (const TempTrajectory &) const override
bool	qualityFilter (const Trajectory &) const override
void	setEvent (const edm::Event &iEvent, const edm::EventSetup &iSetup) override
bool	toBeContinued (TempTrajectory &) const override
bool	toBeContinued (Trajectory &) const override
	~ClusterShapeTrajectoryFilter () override
Public Member Functions inherited from TrajectoryFilter
virtual bool	operator() (TempTrajectory &t) const
virtual bool	operator() (Trajectory &t) const
virtual	~TrajectoryFilter ()

Private Attributes
const SiPixelClusterShapeCache *	theCache
edm::EDGetTokenT< SiPixelClusterShapeCache >	theCacheToken
const ClusterShapeHitFilter *	theFilter

Additional Inherited Members
Public Types inherited from TrajectoryFilter
typedef CkfComponentsRecord	Record
Static Public Attributes inherited from TrajectoryFilter
static const bool	qualityFilterIfNotContributing =true
static const bool	toBeContinuedIfNotContributing =true

Detailed Description

Definition at line 18 of file ClusterShapeTrajectoryFilter.h.

Constructor & Destructor Documentation

ClusterShapeTrajectoryFilter::ClusterShapeTrajectoryFilter	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iC
	)

Definition at line 39 of file ClusterShapeTrajectoryFilter.cc.

                                                                                                                  :
  theCacheToken(iC.consumes<SiPixelClusterShapeCache>(iConfig.getParameter<edm::InputTag>("cacheSrc"))),
  theFilter(nullptr)
{}

ClusterShapeTrajectoryFilter::~ClusterShapeTrajectoryFilter ( )

override

Definition at line 44 of file ClusterShapeTrajectoryFilter.cc.

{
}

Member Function Documentation

std::string ClusterShapeTrajectoryFilter::name ( void  ) const

inlineoverridevirtual

Implements TrajectoryFilter.

Definition at line 32 of file ClusterShapeTrajectoryFilter.h.

Referenced by config.CFG::__str__(), and validation.Sample::digest().

{ return "ClusterShapeTrajectoryFilter"; }

bool ClusterShapeTrajectoryFilter::qualityFilter ( const TempTrajectory &  trajectory ) const

overridevirtual

Implements TrajectoryFilter.

Definition at line 235 of file ClusterShapeTrajectoryFilter.cc.

References TempTrajectory::measurements(), cmsutils::bqueue< T >::size(), protons_cff::t, and toBeContinued().

Referenced by qualityFilter(), and toBeContinued().

{
  TempTrajectory t = trajectory;

  // Check if ok
  if(toBeContinued(t)) return true;

  // Should take out last
  if(t.measurements().size() <= 3) return false;

  return true;
}

bool ClusterShapeTrajectoryFilter::qualityFilter ( const Trajectory &  trajectory ) const

overridevirtual

Implements TrajectoryFilter.

Definition at line 228 of file ClusterShapeTrajectoryFilter.cc.

References qualityFilter().

{
  return true;
}

void ClusterShapeTrajectoryFilter::setEvent ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overridevirtual

Reimplemented from TrajectoryFilter.

Definition at line 48 of file ClusterShapeTrajectoryFilter.cc.

References utilities::cache(), edm::EventSetup::get(), edm::Event::getByToken(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), theCache, theCacheToken, theFilter, and toBeContinued().

                                                                                               {
  edm::ESHandle<ClusterShapeHitFilter> shape;
  iSetup.get<TrajectoryFilter::Record>().get("ClusterShapeHitFilter", shape);
  theFilter = shape.product();

  edm::Handle<SiPixelClusterShapeCache> cache;
  iEvent.getByToken(theCacheToken, cache);
  theCache = cache.product();
}

bool ClusterShapeTrajectoryFilter::toBeContinued ( TempTrajectory &  trajectory ) const

overridevirtual

Implements TrajectoryFilter.

Definition at line 131 of file ClusterShapeTrajectoryFilter.cc.

References TrackingRecHit::det(), TrajectoryStateOnSurface::globalDirection(), TrackingRecHit::hit(), ClusterShapeHitFilter::isCompatible(), GeomDetEnumerators::isTrackerStrip(), TrackingRecHit::isValid(), LogTrace, TempTrajectory::measurements(), SiStripMatchedRecHit2D::monoHit(), ProjectedSiStripRecHit2D::originalHit(), GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, qualityFilter(), cmsutils::bqueue< T >::rbegin(), rpcPointValidation_cfi::recHit, cmsutils::bqueue< T >::rend(), SiStripMatchedRecHit2D::stereoHit(), GeomDet::subDetector(), theCache, and theFilter.

Referenced by qualityFilter(), setEvent(), and toBeContinued().

{
  assert(theCache);
  const TempTrajectory::DataContainer& tms = trajectory.measurements();

  for(TempTrajectory::DataContainer::const_iterator
       tm = tms.rbegin(); tm!= tms.rend(); --tm)
  {
    const TrackingRecHit* ttRecHit = &(*((*tm).recHit()));

    if(ttRecHit->isValid())
    {
      const TrackingRecHit * tRecHit = ttRecHit->hit();

      TrajectoryStateOnSurface ts = (*tm).updatedState();
      GlobalVector gdir = ts.globalDirection();

      if(ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::PixelBarrel ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::PixelEndcap ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P1PXB ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P1PXEC ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P2PXB ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P2PXEC) 
      { // pixel
        const SiPixelRecHit* recHit =
           dynamic_cast<const SiPixelRecHit *>(tRecHit);

        if(recHit != nullptr)
          if(! theFilter->isCompatible(*recHit, gdir, *theCache))
          {
            LogTrace("TrajectFilter")
              << "  [TrajectFilter] fail pixel";
            return false;
          }
      }
      else if(GeomDetEnumerators::isTrackerStrip(ttRecHit->det()->subDetector()))
      { // strip
        if(dynamic_cast<const SiStripMatchedRecHit2D *>(tRecHit)  != nullptr)
        { // glued
          const SiStripMatchedRecHit2D* recHit =
            dynamic_cast<const SiStripMatchedRecHit2D *>(tRecHit);

          if(recHit != nullptr)
          { 
            if(! theFilter->isCompatible(recHit->monoHit(), gdir))
            {
              LogTrace("TrajectFilter")
               << "  [TrajectFilter] fail strip matched 1st";
              return false;
            }

            if(! theFilter->isCompatible(recHit->stereoHit(), gdir))
            {
              LogTrace("TrajectFilter")
                << "  [TrajectFilter] fail strip matched 2nd";
              return false;
            }
          }
        }
        else
        { // single
          if(dynamic_cast<const SiStripRecHit2D *>(tRecHit) != nullptr)
          { // normal
            const SiStripRecHit2D* recHit =
              dynamic_cast<const SiStripRecHit2D *>(tRecHit);
  
            if(recHit != nullptr)
              if(! theFilter->isCompatible(*recHit, gdir))
              {
                LogTrace("TrajectFilter")
                  << "  [TrajectFilter] fail strip single";
                return false;
              }
          }
          else
          { // projected
            const ProjectedSiStripRecHit2D* recHit =
              dynamic_cast<const ProjectedSiStripRecHit2D *>(tRecHit);
 
            if(recHit != nullptr)
              if(! theFilter->isCompatible(recHit->originalHit(), gdir))
              {
                LogTrace("TrajectFilter")
                  << "  [TrajectFilter] fail strip projected";
                return false;
              }
          }
        }
      }
    }
  }

  return true;
}

bool ClusterShapeTrajectoryFilter::toBeContinued ( Trajectory &  trajectory ) const

overridevirtual

Implements TrajectoryFilter.

Definition at line 60 of file ClusterShapeTrajectoryFilter.cc.

References TrackingRecHit::det(), TrajectoryStateOnSurface::globalDirection(), TrackingRecHit::hit(), ClusterShapeHitFilter::isCompatible(), GeomDetEnumerators::isTrackerStrip(), TrackingRecHit::isValid(), Trajectory::measurements(), SiStripMatchedRecHit2D::monoHit(), ProjectedSiStripRecHit2D::originalHit(), GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, rpcPointValidation_cfi::recHit, SiStripMatchedRecHit2D::stereoHit(), GeomDet::subDetector(), theCache, theFilter, and toBeContinued().

{
  assert(theCache);
  vector<TrajectoryMeasurement> tms = trajectory.measurements();

  for(vector<TrajectoryMeasurement>::const_iterator
       tm = tms.begin(); tm!= tms.end(); tm++)
  {
    const TrackingRecHit* ttRecHit = &(*((*tm).recHit()));

    if(ttRecHit->isValid())
    {
      const TrackingRecHit * tRecHit = ttRecHit->hit();

      TrajectoryStateOnSurface ts = (*tm).updatedState();
      const GlobalVector gdir = ts.globalDirection();

      if(ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::PixelBarrel ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::PixelEndcap ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P1PXB ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P1PXEC ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P2PXB ||
     ttRecHit->det()->subDetector()==GeomDetEnumerators::SubDetector::P2PXEC) 
      { // pixel
        const SiPixelRecHit* recHit =
           dynamic_cast<const SiPixelRecHit *>(tRecHit);

        if(recHit != nullptr)
          return theFilter->isCompatible(*recHit, gdir, *theCache);
      }
      else if(GeomDetEnumerators::isTrackerStrip(ttRecHit->det()->subDetector()))
      { // strip
        if(dynamic_cast<const SiStripMatchedRecHit2D *>(tRecHit)  != nullptr)
        { // glued
          const SiStripMatchedRecHit2D* recHit =
            dynamic_cast<const SiStripMatchedRecHit2D *>(tRecHit);

          if(recHit != nullptr)
          { 
            return (theFilter->isCompatible(recHit->monoHit()  , gdir) &&
                    theFilter->isCompatible(recHit->stereoHit(), gdir));
          }
        }
        else
        { // single
          if(dynamic_cast<const SiStripRecHit2D *>(tRecHit) != nullptr)
          { // normal
            const SiStripRecHit2D* recHit =
              dynamic_cast<const SiStripRecHit2D *>(tRecHit);
  
            if(recHit != nullptr)
              return theFilter->isCompatible(*recHit, gdir);
          }
          else
          { // projected
            const ProjectedSiStripRecHit2D* recHit =
              dynamic_cast<const ProjectedSiStripRecHit2D *>(tRecHit);
 
            if(recHit != nullptr)
              return theFilter->isCompatible(recHit->originalHit(), gdir);
          }
        }
      }
    }
  }

  return true;
}

Member Data Documentation

const SiPixelClusterShapeCache* ClusterShapeTrajectoryFilter::theCache

private

Definition at line 36 of file ClusterShapeTrajectoryFilter.h.

Referenced by setEvent(), and toBeContinued().

edm::EDGetTokenT<SiPixelClusterShapeCache> ClusterShapeTrajectoryFilter::theCacheToken

private

Definition at line 35 of file ClusterShapeTrajectoryFilter.h.

Referenced by setEvent().

const ClusterShapeHitFilter* ClusterShapeTrajectoryFilter::theFilter

private

Definition at line 37 of file ClusterShapeTrajectoryFilter.h.

Referenced by setEvent(), and toBeContinued().