#include <PixelTrackFilterByKinematics.h>

Inheritance diagram for PixelTrackFilterByKinematics:

Public Member Functions
virtual bool	operator() (const reco::Track *) const

virtual bool	operator() (const reco::Track *, const PixelTrackFilter::Hits &hits) const

	PixelTrackFilterByKinematics (const edm::ParameterSet &cfg)

	PixelTrackFilterByKinematics (const edm::ParameterSet &cfg, const edm::EventSetup &es)

	PixelTrackFilterByKinematics (double ptmin=0.9, double tipmax=0.1, double chi2max=100.)

virtual	~PixelTrackFilterByKinematics ()

Public Member Functions inherited from PixelTrackFilter
virtual void	update (edm::Event &ev)

virtual	~PixelTrackFilter ()

Private Attributes
double	theChi2Max

double	theNSigmaInvPtTolerance

double	theNSigmaTipMaxTolerance

double	thePtMin

double	theTIPMax

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

Detailed Description

Definition at line 8 of file PixelTrackFilterByKinematics.h.

Constructor & Destructor Documentation

PixelTrackFilterByKinematics::PixelTrackFilterByKinematics ( const edm::ParameterSet & cfg )

Definition at line 19 of file PixelTrackFilterByKinematics.cc.

   : thePtMin( cfg.getParameter<double>("ptMin") ),
     theNSigmaInvPtTolerance( cfg.getParameter<double>("nSigmaInvPtTolerance")),
     theTIPMax( cfg.getParameter<double>("tipMax") ),
     theNSigmaTipMaxTolerance( cfg.getParameter<double>("nSigmaTipMaxTolerance")),
     theChi2Max( cfg.getParameter<double>("chi2") )
 { }

PixelTrackFilterByKinematics::PixelTrackFilterByKinematics	(	const edm::ParameterSet &	cfg,
		const edm::EventSetup &	es
	)

Definition at line 11 of file PixelTrackFilterByKinematics.cc.

   : thePtMin( cfg.getParameter<double>("ptMin") ),
     theNSigmaInvPtTolerance( cfg.getParameter<double>("nSigmaInvPtTolerance")),
     theTIPMax( cfg.getParameter<double>("tipMax") ),
     theNSigmaTipMaxTolerance( cfg.getParameter<double>("nSigmaTipMaxTolerance")),
     theChi2Max( cfg.getParameter<double>("chi2") )
 { }

PixelTrackFilterByKinematics::PixelTrackFilterByKinematics	(	double	ptmin = `0.9`,
		double	tipmax = `0.1`,
		double	chi2max = `100.`
	)

Definition at line 27 of file PixelTrackFilterByKinematics.cc.

   : thePtMin(ptmin), theNSigmaInvPtTolerance(0.),
     theTIPMax(tipmax), theNSigmaTipMaxTolerance(0.),
     theChi2Max(chi2max)
 { } 

PixelTrackFilterByKinematics::~PixelTrackFilterByKinematics ( )

virtual

Definition at line 33 of file PixelTrackFilterByKinematics.cc.

34 { }

Member Function Documentation

bool PixelTrackFilterByKinematics::operator() ( const reco::Track * track ) const

virtual

Reimplemented from PixelTrackFilter.

Definition at line 39 of file PixelTrackFilterByKinematics.cc.

References reco::TrackBase::chi2(), funct::cos(), reco::TrackBase::covariance(), reco::TrackBase::d0(), reco::TrackBase::d0Error(), reco::TrackBase::i_lambda, reco::TrackBase::i_qoverp, reco::TrackBase::lambdaError(), reco::TrackBase::pt(), reco::TrackBase::qoverpError(), funct::sin(), funct::sqr(), mathSSE::sqrt(), theChi2Max, theNSigmaInvPtTolerance, theNSigmaTipMaxTolerance, thePtMin, reco::TrackBase::theta(), theta(), and theTIPMax.

 {
   if (!track) return false;
   if (track->chi2() > theChi2Max) return false;
   if ( (fabs(track->d0())-theTIPMax)/track->d0Error() > theNSigmaTipMaxTolerance) return false;
   
   float theta = track->theta();
   float cosTheta = cos(theta);
   float sinTheta = sin(theta);
   float errLambda2 = sqr( track->lambdaError() );
 
   float pt_v = track->pt();
   float errInvP2 = sqr(track->qoverpError());
   float covIPtTheta = track->covariance(reco::TrackBase::i_qoverp, reco::TrackBase::i_lambda);
   float errInvPt2 = (   errInvP2
                       + sqr(cosTheta/pt_v)*errLambda2
                       + 2*(cosTheta/pt_v)*covIPtTheta
                      ) / sqr(sinTheta);
   if ( (1/pt_v - 1/thePtMin)/sqrt(errInvPt2) > theNSigmaInvPtTolerance ) return false;
 
   return true;
 }