StripSubClusterShapeSeedFilter Class Reference

#include <StripSubClusterShapeTrajectoryFilter.h>

Inheritance diagram for StripSubClusterShapeSeedFilter:

Public Member Functions
virtual bool	compatible (const TrajectoryStateOnSurface &tsos, SeedingHitSet::ConstRecHitPointer hit) const
virtual bool	compatible (const SeedingHitSet &hits, const TrackingRegion &region) const
virtual bool	compatible (const TrajectorySeed &seed) const
virtual bool	compatible (const SeedingHitSet &hits, const GlobalTrajectoryParameters &helixStateAtVertex, const FastHelix &helix, const TrackingRegion &region) const
virtual bool	compatible (const SeedingHitSet &hits, const GlobalTrajectoryParameters &straightLineStateAtVertex, const TrackingRegion &region) const
virtual void	init (const edm::Event &ev, const edm::EventSetup &es) override
	StripSubClusterShapeSeedFilter (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
virtual	~StripSubClusterShapeSeedFilter ()
Public Member Functions inherited from StripSubClusterShapeFilterBase
	StripSubClusterShapeFilterBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
virtual	~StripSubClusterShapeFilterBase ()
Public Member Functions inherited from SeedComparitor
virtual	~SeedComparitor ()

Protected Attributes
bool	filterAtHelixStage_
Protected Attributes inherited from StripSubClusterShapeFilterBase
std::string	label_
std::array< std::array < uint8_t, 10 >, 7 >	layerMask_
uint32_t	maxNSat_
float	maxTrimmedSizeDiffNeg_
float	maxTrimmedSizeDiffPos_
float	seedCutMIPs_
float	seedCutSN_
float	subclusterCutMIPs_
float	subclusterCutSN_
float	subclusterWindow_
edm::ESHandle < ClusterShapeHitFilter >	theFilter
edm::ESHandle< SiStripNoises >	theNoise
edm::ESHandle< TrackerTopology >	theTopology
edm::ESHandle< TrackerGeometry >	theTracker
uint8_t	trimMaxADC_
float	trimMaxFracNeigh_
float	trimMaxFracTotal_

Additional Inherited Members
Protected Member Functions inherited from StripSubClusterShapeFilterBase
void	setEventBase (const edm::Event &, const edm::EventSetup &)
bool	testLastHit (const TrackingRecHit *hit, const TrajectoryStateOnSurface &tsos, bool mustProject=false) const
bool	testLastHit (const TrackingRecHit *hit, const GlobalPoint &gpos, const GlobalVector &gdir, bool mustProject=false) const

Detailed Description

Definition at line 93 of file StripSubClusterShapeTrajectoryFilter.h.

Constructor & Destructor Documentation

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

Definition at line 347 of file StripSubClusterShapeTrajectoryFilter.cc.

                                                              :
        StripSubClusterShapeFilterBase(iConfig,iC), 
        filterAtHelixStage_(iConfig.getParameter<bool>("FilterAtHelixStage")) 

{
    if (filterAtHelixStage_) edm::LogError("Configuration") << "StripSubClusterShapeSeedFilter: FilterAtHelixStage is not yet working correctly.\n";
}

virtual StripSubClusterShapeSeedFilter::~StripSubClusterShapeSeedFilter ( )

inlinevirtual

Definition at line 97 of file StripSubClusterShapeTrajectoryFilter.h.

{}

Member Function Documentation

bool StripSubClusterShapeSeedFilter::compatible ( const TrajectoryStateOnSurface &  tsos, SeedingHitSet::ConstRecHitPointer  hit  ) const

virtual

Implements SeedComparitor.

Definition at line 358 of file StripSubClusterShapeTrajectoryFilter.cc.

References TrackingRecHit::geographicalId(), TrackingRecHit::hit(), TrackingRecHit::isValid(), DetId::subdetId(), and SiStripDetId::TIB.

{
   if (filterAtHelixStage_) return true;
   const TrackingRecHit* hit = thit->hit();
   if (hit == 0 || !hit->isValid()) return true;
   if (hit->geographicalId().subdetId() < SiStripDetId::TIB) return true; // we look only at strips for now
   return testLastHit(hit, tsos, false);
}

virtual bool StripSubClusterShapeSeedFilter::compatible ( const SeedingHitSet &  hits, const TrackingRegion &  region  ) const

inlinevirtual

Implements SeedComparitor.

Definition at line 105 of file StripSubClusterShapeTrajectoryFilter.h.

{ return true; }

virtual bool StripSubClusterShapeSeedFilter::compatible ( const TrajectorySeed &  seed ) const

inlinevirtual

Implements SeedComparitor.

Definition at line 106 of file StripSubClusterShapeTrajectoryFilter.h.

{ return true; }

bool StripSubClusterShapeSeedFilter::compatible ( const SeedingHitSet &  hits, const GlobalTrajectoryParameters &  helixStateAtVertex, const FastHelix &  helix, const TrackingRegion &  region  ) const

virtual

Implements SeedComparitor.

Definition at line 368 of file StripSubClusterShapeTrajectoryFilter.cc.

References FastHelix::circle(), TrackingRecHit::geographicalId(), TrackingRecHit::globalPosition(), i, FastHelix::isValid(), GlobalTrajectoryParameters::momentum(), gen::n, GlobalTrajectoryParameters::position(), SeedingHitSet::size(), mathSSE::sqrt(), DetId::subdetId(), SiStripDetId::TIB, PV3DBase< T, PVType, FrameType >::x(), FastCircle::x0(), PV3DBase< T, PVType, FrameType >::y(), FastCircle::y0(), and PV3DBase< T, PVType, FrameType >::z().

{ 
   if (!filterAtHelixStage_) return true;

   if (!helix.isValid()) edm::LogWarning("InvalidHelix") << "PixelClusterShapeSeedComparitor helix is not valid, result is bad";

    float xc = helix.circle().x0(), yc = helix.circle().y0();

    GlobalPoint  vertex = helixStateAtVertex.position();
    GlobalVector momvtx = helixStateAtVertex.momentum();
    float x0 = vertex.x(), y0 = vertex.y();
    for (unsigned int i = 0, n = hits.size(); i < n; ++i) {
        auto const  & hit = *hits[i];
        if (hit.geographicalId().subdetId() < SiStripDetId::TIB) continue; 

        GlobalPoint pos = hit.globalPosition();
        float x1 = pos.x(), y1 = pos.y(), dx1 = x1 - xc, dy1 = y1 - yc;

        // now figure out the proper tangent vector
        float perpx = -dy1, perpy = dx1;
        if (perpx * (x1-x0) + perpy * (y1 - y0) < 0) {
            perpy = -perpy; perpx = -perpx;
        }

        // now normalize (perpx, perpy, 1.0) to unity
        float pnorm = 1.0/std::sqrt(perpx*perpx + perpy*perpy + 1);
        GlobalVector gdir(perpx*pnorm, perpy*pnorm, (momvtx.z() > 0 ? pnorm : -pnorm));

        if (!testLastHit(&hit, pos, gdir)) {
            return false; // not yet
        }
    }
    return true; 
}

virtual bool StripSubClusterShapeSeedFilter::compatible ( const SeedingHitSet &  hits, const GlobalTrajectoryParameters &  straightLineStateAtVertex, const TrackingRegion &  region  ) const

inlinevirtual

Implements SeedComparitor.

Definition at line 108 of file StripSubClusterShapeTrajectoryFilter.h.

{ return true; }

virtual void StripSubClusterShapeSeedFilter::init ( const edm::Event &  ev, const edm::EventSetup &  es  )

inlineoverridevirtual

Implements SeedComparitor.

Definition at line 99 of file StripSubClusterShapeTrajectoryFilter.h.

References StripSubClusterShapeFilterBase::setEventBase().

                                                                                {
            setEventBase(ev,es);
        }

Member Data Documentation

bool StripSubClusterShapeSeedFilter::filterAtHelixStage_

protected

Definition at line 111 of file StripSubClusterShapeTrajectoryFilter.h.