Reference Manual

Go to the documentation of this file.
00001 #ifndef MinPtTrajectoryFilter_H
00002 #define MinPtTrajectoryFilter_H
00003 
00004 #include "TrackingTools/TrajectoryFiltering/interface/TrajectoryFilter.h"
00005 #include "TrackingTools/PatternTools/interface/Trajectory.h"
00006 #include "TrackingTools/PatternTools/interface/TempTrajectory.h"
00007 
00008 #include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
00009 #include "TrackingTools/TrajectoryParametrization/interface/CurvilinearTrajectoryError.h"
00010 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateAccessor.h"
00011 
00018 class MinPtTrajectoryFilter : public TrajectoryFilter {
00019 public:
00020 
00021   explicit MinPtTrajectoryFilter( double ptMin, float nSigma = 5.F, int nH=3): thePtMin( ptMin), theNSigma(nSigma), theMinHits(nH)  {}
00022 
00023   explicit MinPtTrajectoryFilter( const edm::ParameterSet & pset) :
00024     thePtMin(pset.getParameter<double>("minPt")),
00025     theNSigma(pset.getParameter<double>("nSigmaMinPt")),
00026     theMinHits(pset.getParameter<int>("minHitsMinPt")){}
00027     
00028 
00029   virtual bool qualityFilter( const Trajectory& traj)const { return test(traj.lastMeasurement(),traj.foundHits()); }
00030   virtual bool qualityFilter( const TempTrajectory& traj) const { return test(traj.lastMeasurement(),traj.foundHits()); }
00031     
00032   virtual bool toBeContinued( Trajectory& traj) const {return test(traj.lastMeasurement(),traj.foundHits()); }
00033   virtual bool toBeContinued( TempTrajectory& traj) const { return test(traj.lastMeasurement(),traj.foundHits()); }
00034   
00035   virtual std::string name() const {return "MinPtTrajectoryFilter";}
00036 
00037  protected:
00038 
00039   bool test( const TrajectoryMeasurement & tm, int foundHits) const 
00040   {
00041     //first check min number of hits 
00042     if (foundHits < theMinHits ){ return true;}
00043 
00044     // check for momentum below limit
00045     const FreeTrajectoryState& fts = *tm.updatedState().freeTrajectoryState();
00046 
00047     //avoid doing twice the check in TBC and QF
00048     static bool answerMemory=false;
00049     static FreeTrajectoryState ftsMemory;
00050     if (ftsMemory.parameters().vector() == fts.parameters().vector()) { return answerMemory;}
00051     ftsMemory=fts;
00052 
00053     //if p_T is way too small: stop
00054     double pT = fts.momentum().perp();
00055     if (pT<0.010) {answerMemory=false; return false;}
00056     //if error is way too big: stop
00057     double invError = TrajectoryStateAccessor(fts).inversePtError();
00058     if (invError > 1.e10) {answerMemory=false;return false;}
00059 
00060     //calculate the actual pT cut: 
00061     if ((1/pT - theNSigma*invError) > 1/thePtMin) {answerMemory=false; return false;}
00062     //    first term if the max value of pT (pT+N*sigma(pT))
00063     //    second tern is the cut
00064 
00065     answerMemory=true; return true;
00066   }
00067 
00068   double thePtMin;
00069   double theNSigma;
00070   int theMinHits;
00071 
00072 };
00073 
00074 #endif
/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/TrackingTools/TrajectoryFiltering/interface/MinPtTrajectoryFilter.h
