ChargeSignificanceTrajectoryFilter Class Reference

#include <ChargeSignificanceTrajectoryFilter.h>

Inheritance diagram for ChargeSignificanceTrajectoryFilter:

Public Member Functions
	ChargeSignificanceTrajectoryFilter (float qsig)
	ChargeSignificanceTrajectoryFilter (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
template<class T >
bool	TBC (T &traj) const

Protected Attributes
float	theChargeSignificance

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 14 of file ChargeSignificanceTrajectoryFilter.h.

Constructor & Destructor Documentation

ChargeSignificanceTrajectoryFilter::ChargeSignificanceTrajectoryFilter ( float  qsig )

inlineexplicit

Definition at line 17 of file ChargeSignificanceTrajectoryFilter.h.

:  theChargeSignificance(qsig) {}

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

inlineexplicit

Definition at line 19 of file ChargeSignificanceTrajectoryFilter.h.

                                                                                                      :
    theChargeSignificance(pset.getParameter<double>("chargeSignificance")) {}

Member Function Documentation

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

inlinevirtual

Implements TrajectoryFilter.

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

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 22 of file ChargeSignificanceTrajectoryFilter.h.

References Trajectory::isValid().

{ return traj.isValid();}

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 23 of file ChargeSignificanceTrajectoryFilter.h.

References TempTrajectory::isValid().

{ return traj.isValid();}

template<class T >

bool ChargeSignificanceTrajectoryFilter::TBC ( T &  traj ) const

inlineprotected

Definition at line 32 of file ChargeSignificanceTrajectoryFilter.h.

References funct::abs(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LocalTrajectoryError::matrix(), mathSSE::sqrt(), theChargeSignificance, and LocalTrajectoryParameters::vector().

                                             {
    const typename T::DataContainer & tms = traj.measurements();
    // Check flip in q-significance. The loop over all TMs could be 
    // avoided by storing the current significant q in the trajectory

    if ( theChargeSignificance>0. ) {      

      float qSig = 0;

      // skip first two hits (don't rely on significance of q/p)
      for( typename T::DataContainer::size_type itm=2; itm<tms.size(); ++itm ) {
    TrajectoryStateOnSurface const & tsos = tms[itm].updatedState();
    if ( !tsos.isValid() )  continue;

        auto significance = tsos.localParameters().vector()(0) /
      std::sqrt(float(tsos.localError().matrix()(0,0)));

    // don't deal with measurements compatible with 0
    if ( std::abs(significance)<theChargeSignificance )  continue;

    //
    // if charge not yet defined: store first significant Q
    //
    if ( qSig==0 ) qSig = significance;
    
    //
    //  invalidate and terminate in case of a change of sign
    //
    if (significance*qSig<0) {
        traj.invalidate();
        return false;
    }
      }
    }
    return true;
  }

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 25 of file ChargeSignificanceTrajectoryFilter.h.

{ return TBC<Trajectory>(traj);}

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

inlinevirtual

Implements TrajectoryFilter.

Definition at line 26 of file ChargeSignificanceTrajectoryFilter.h.

{ return TBC<TempTrajectory>(traj);}

Member Data Documentation

float ChargeSignificanceTrajectoryFilter::theChargeSignificance

protected

Definition at line 69 of file ChargeSignificanceTrajectoryFilter.h.

Referenced by TBC().