MinPtTrajectoryFilter Class Reference

#include <MinPtTrajectoryFilter.h>

Inheritance diagram for MinPtTrajectoryFilter:

Public Member Functions
	MinPtTrajectoryFilter (float ptMin, float nSigma=5.F, int nH=3)
	MinPtTrajectoryFilter (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
virtual std::string	name () const
virtual bool	qualityFilter (const Trajectory &traj) const
virtual bool	qualityFilter (const TempTrajectory &traj) const
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
float	theInvPtMin
int	theMinHits
float	theNSigma
float	thePtMin2

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 drops below some value at some confidence level. The CkfTrajectoryBuilder uses this class to implement the minimal P_t cut.

Definition at line 18 of file MinPtTrajectoryFilter.h.

Constructor & Destructor Documentation

MinPtTrajectoryFilter::MinPtTrajectoryFilter ( float  ptMin, float  nSigma = 5.F, int  nH = 3  )

inlineexplicit

Definition at line 21 of file MinPtTrajectoryFilter.h.

                                                                            : 
    thePtMin2(ptMin*ptMin),theInvPtMin(1.f/ptMin), theNSigma(nSigma), theMinHits(nH)  {}

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

inlineexplicit

Definition at line 25 of file MinPtTrajectoryFilter.h.

References thePtMin2.

                                                                                          :
    thePtMin2(pset.getParameter<double>("minPt")),
    theInvPtMin(1.f/thePtMin2),
    theNSigma(pset.getParameter<double>("nSigmaMinPt")),
    theMinHits(pset.getParameter<int>("minHitsMinPt")){thePtMin2*=thePtMin2;}

Member Function Documentation

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 38 of file MinPtTrajectoryFilter.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 "MinPtTrajectoryFilter";}

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 32 of file MinPtTrajectoryFilter.h.

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

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

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 33 of file MinPtTrajectoryFilter.h.

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

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

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

protected

Definition at line 14 of file MinPtTrajectoryFilter.cc.

References f, PV3DBase< T, PVType, FrameType >::perp2(), mathSSE::sqrt(), theInvPtMin, theMinHits, theNSigma, thePtMin2, and TrajectoryMeasurement::updatedState().

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();

    auto const & tsos = tm.updatedState();
    GlobalVector gtp = tsos.globalMomentum();

    //avoid doing twice the check in TBC and QF


    if ( gtp == tls.ftsMemory) { return tls.answerMemory;}
    tls.ftsMemory= gtp;

    auto pT2 =  gtp.perp2();

    //if p_T is way too small: stop
    if (pT2<0.0010f) {tls.answerMemory=false; return false;}

    // if large enouth go
    if (pT2> thePtMin2) { tls.answerMemory=true; return true;}

    //if error is way too big: stop
    float invError = TrajectoryStateAccessor(*tsos.freeTrajectoryState()).inversePtError();
    if (invError > 1.e10f) {tls.answerMemory=false;return false;}

    //calculate the actual pT cut:
    if ((1.f/std::sqrt(pT2) - theNSigma*invError) > theInvPtMin ) {tls.answerMemory=false; return false;}
    //    first term if the max value of pT (pT+N*sigma(pT))
    //    second tern is the cut

    tls.answerMemory=true; return true;
  }

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 35 of file MinPtTrajectoryFilter.h.

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

Referenced by RegionalTrajectoryFilter::toBeContinued().

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

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 36 of file MinPtTrajectoryFilter.h.

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

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

Member Data Documentation

float MinPtTrajectoryFilter::theInvPtMin

protected

Definition at line 45 of file MinPtTrajectoryFilter.h.

Referenced by test().

int MinPtTrajectoryFilter::theMinHits

protected

Definition at line 47 of file MinPtTrajectoryFilter.h.

Referenced by test().

float MinPtTrajectoryFilter::theNSigma

protected

Definition at line 46 of file MinPtTrajectoryFilter.h.

Referenced by test().

float MinPtTrajectoryFilter::thePtMin2

protected

Definition at line 44 of file MinPtTrajectoryFilter.h.

Referenced by MinPtTrajectoryFilter(), and test().