ThresholdPtTrajectoryFilter Class Reference

#include <ThresholdPtTrajectoryFilter.h>

Inheritance diagram for ThresholdPtTrajectoryFilter:

Public Member Functions
virtual std::string	name () const
virtual bool	qualityFilter (const Trajectory &traj) const
virtual bool	qualityFilter (const TempTrajectory &traj) const
	ThresholdPtTrajectoryFilter (double ptThreshold, float nSigma=5.F, int nH=3)
	ThresholdPtTrajectoryFilter (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
virtual bool	toBeContinued (Trajectory &traj) const
virtual bool	toBeContinued (TempTrajectory &traj) const
Public Member Functions inherited from TrajectoryFilter
virtual bool	operator() (TempTrajectory &t) const
virtual bool	operator() (Trajectory &t) const
virtual void	setEvent (const edm::Event &iEvent, const edm::EventSetup &iSetup)
virtual	~TrajectoryFilter ()

Protected Member Functions
bool	test (const TrajectoryMeasurement &tm, int foundHits) const

Protected Attributes
int	theMinHits
double	theNSigma
double	thePtThreshold

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

A TrajectoryFilter that stops reconstruction if P_t goes above some value at some confidence level. The CkfTrajectoryBuilder uses this class to implement the conditional p_T cut

Definition at line 18 of file ThresholdPtTrajectoryFilter.h.

Constructor & Destructor Documentation

ThresholdPtTrajectoryFilter::ThresholdPtTrajectoryFilter ( double  ptThreshold, float  nSigma = 5.F, int  nH = 3  )

inlineexplicit

Definition at line 21 of file ThresholdPtTrajectoryFilter.h.

: thePtThreshold( ptThreshold), theNSigma(nSigma), theMinHits(nH) {}

ThresholdPtTrajectoryFilter::ThresholdPtTrajectoryFilter ( const edm::ParameterSet &  pset, edm::ConsumesCollector &  iC  )

inlineexplicit

Definition at line 23 of file ThresholdPtTrajectoryFilter.h.

                                                                                                :
    thePtThreshold(pset.getParameter<double>("thresholdPt")),
    theNSigma(pset.getParameter<double>("nSigmaThresholdPt")),
    theMinHits(pset.getParameter<int>("minHitsThresholdPt"))
      {}

Member Function Documentation

virtual std::string ThresholdPtTrajectoryFilter::name ( void  ) const

inlinevirtual

Implements TrajectoryFilter.

Definition at line 35 of file ThresholdPtTrajectoryFilter.h.

Referenced by cuy.divideElement::__init__(), cuy.plotElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), config.CFG::__str__(), and validation.Sample::digest().

{return "ThresholdPtTrajectoryFilter";}

virtual bool ThresholdPtTrajectoryFilter::qualityFilter ( const Trajectory &  traj ) const

inlinevirtual

Implements TrajectoryFilter.

Definition at line 29 of file ThresholdPtTrajectoryFilter.h.

References Trajectory::foundHits(), Trajectory::lastMeasurement(), and test().

{ return !test(traj.lastMeasurement(),traj.foundHits());}

virtual bool ThresholdPtTrajectoryFilter::qualityFilter ( const TempTrajectory &  traj ) const

inlinevirtual

Implements TrajectoryFilter.

Definition at line 30 of file ThresholdPtTrajectoryFilter.h.

References TempTrajectory::foundHits(), TempTrajectory::lastMeasurement(), and test().

{ return !test(traj.lastMeasurement(),traj.foundHits());}

bool ThresholdPtTrajectoryFilter::test ( const TrajectoryMeasurement &  tm, int  foundHits  ) const

protected

Definition at line 13 of file ThresholdPtTrajectoryFilter.cc.

References TrajectoryStateOnSurface::freeTrajectoryState(), TrajectoryStateAccessor::inversePtError(), FreeTrajectoryState::momentum(), FreeTrajectoryState::parameters(), PV3DBase< T, PVType, FrameType >::perp(), theMinHits, theNSigma, thePtThreshold, TrajectoryMeasurement::updatedState(), and GlobalTrajectoryParameters::vector().

Referenced by qualityFilter(), and toBeContinued().

  {
    //first check min number of hits
    if (foundHits < theMinHits ){ return true;}

    // check for momentum below limit
    const FreeTrajectoryState& fts = *tm.updatedState().freeTrajectoryState();

    //avoid doing twice the check in TBC and QF
    // We make it thread local so that we avoid race conditions between
    // threads, and we make sure there is no cache contention between them.
    if (tls.ftsMemory.parameters().vector() == fts.parameters().vector()) { return tls.answerMemory;}
    tls.ftsMemory=fts;

    //if p_T is way too small: stop
    double pT = fts.momentum().perp();
    if (pT<0.010) {tls.answerMemory=false; return false;}
    //if error is way too big: stop
    double invError = TrajectoryStateAccessor(fts).inversePtError();
    if (invError > 1.e10) {tls.answerMemory=false;return false;}

    //calculate the actual pT cut:
    if ((1/pT + theNSigma*invError ) < 1/thePtThreshold ) {tls.answerMemory=false; return false;}
    //    first term is the minimal value of pT (pT-N*sigma(pT))
    //    secon term is the cut

    tls.answerMemory=true; return true;
  }

virtual bool ThresholdPtTrajectoryFilter::toBeContinued ( Trajectory &  traj ) const

inlinevirtual

Implements TrajectoryFilter.

Definition at line 32 of file ThresholdPtTrajectoryFilter.h.

References Trajectory::foundHits(), Trajectory::lastMeasurement(), and test().

{ return test(traj.lastMeasurement(),traj.foundHits()); }

virtual bool ThresholdPtTrajectoryFilter::toBeContinued ( TempTrajectory &  traj ) const

inlinevirtual

Implements TrajectoryFilter.

Definition at line 33 of file ThresholdPtTrajectoryFilter.h.

References TempTrajectory::foundHits(), TempTrajectory::lastMeasurement(), and test().

{ return test(traj.lastMeasurement(),traj.foundHits()); }

Member Data Documentation

int ThresholdPtTrajectoryFilter::theMinHits

protected

Definition at line 43 of file ThresholdPtTrajectoryFilter.h.

Referenced by test().

double ThresholdPtTrajectoryFilter::theNSigma

protected

Definition at line 42 of file ThresholdPtTrajectoryFilter.h.

Referenced by test().

double ThresholdPtTrajectoryFilter::thePtThreshold

protected

Definition at line 41 of file ThresholdPtTrajectoryFilter.h.

Referenced by test().