ThresholdPtTrajectoryFilter Class Reference

#include <ThresholdPtTrajectoryFilter.h>

Inheritance diagram for ThresholdPtTrajectoryFilter:

Public Member Functions
std::string	name () const override
bool	qualityFilter (const Trajectory &traj) const override
bool	qualityFilter (const TempTrajectory &traj) const override
	ThresholdPtTrajectoryFilter (double ptThreshold, float nSigma=5.F, int nH=3)
	ThresholdPtTrajectoryFilter (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
bool	toBeContinued (Trajectory &traj) const override
bool	toBeContinued (TempTrajectory &traj) const override
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

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

inlineoverridevirtual

Implements TrajectoryFilter.

Definition at line 35 of file ThresholdPtTrajectoryFilter.h.

References test().

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

{return "ThresholdPtTrajectoryFilter";}

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

inlineoverridevirtual

Implements TrajectoryFilter.

Definition at line 29 of file ThresholdPtTrajectoryFilter.h.

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

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

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

inlineoverridevirtual

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(), TrajectoryStateOnSurface::isValid(), FreeTrajectoryState::momentum(), FreeTrajectoryState::parameters(), PV3DBase< T, PVType, FrameType >::perp(), PVValHelper::pT, createJobs::theNSigma, TrajectoryMeasurement::updatedState(), and GlobalTrajectoryParameters::vector().

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

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

    // check for momentum below limit
    if (!tm.updatedState().isValid()) return    false;
    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;
  }

bool ThresholdPtTrajectoryFilter::toBeContinued ( Trajectory &  traj ) const

inlineoverridevirtual

Implements TrajectoryFilter.

Definition at line 32 of file ThresholdPtTrajectoryFilter.h.

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

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

bool ThresholdPtTrajectoryFilter::toBeContinued ( TempTrajectory &  traj ) const

inlineoverridevirtual

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.

double ThresholdPtTrajectoryFilter::theNSigma

protected

Definition at line 42 of file ThresholdPtTrajectoryFilter.h.

double ThresholdPtTrajectoryFilter::thePtThreshold

protected

Definition at line 41 of file ThresholdPtTrajectoryFilter.h.