TrkOutAvgTimeCalculator Class Reference

#include <TrkOutAvgTimeCalculator.h>

Inheritance diagram for TrkOutAvgTimeCalculator:

Public Member Functions
	TrkOutAvgTimeCalculator (const TrkHotSelector &selector, double maxpull)
virtual	~TrkOutAvgTimeCalculator ()
virtual bool	trackTime (const TrkRecoTrk &trk, double &time, double &timeerr, int &nHotsUsed) const
Public Member Functions inherited from TrkTimeCalculator
	TrkTimeCalculator (const TrkHotSelector &selector)
virtual	~TrkTimeCalculator ()
bool	trackTime (const TrkRecoTrk &trk, double &time, double &timeerr) const

Additional Inherited Members
Protected Member Functions inherited from TrkTimeCalculator
bool	useHot (const TrkHitOnTrk &x) const

Detailed Description

Definition at line 27 of file TrkOutAvgTimeCalculator.h.

Constructor & Destructor Documentation

TrkOutAvgTimeCalculator()

TrkOutAvgTimeCalculator::TrkOutAvgTimeCalculator	(	const TrkHotSelector &	selector,
		double	maxpull )

Definition at line 25 of file TrkOutAvgTimeCalculator.cxx.

                                                                                               :
  TrkTimeCalculator(selector),_maxpull(maxpull)
{}

~TrkOutAvgTimeCalculator()

TrkOutAvgTimeCalculator::~TrkOutAvgTimeCalculator ( )

virtual

Definition at line 29 of file TrkOutAvgTimeCalculator.cxx.

{}

Member Function Documentation

trackTime()

bool TrkOutAvgTimeCalculator::trackTime ( const TrkRecoTrk & trk, double & time, double & timeerr, int & nHotsUsed ) const

virtual

Implements TrkTimeCalculator.

Definition at line 33 of file TrkOutAvgTimeCalculator.cxx.

{
  const TrkHitList *hotlist = trk.hits();
  std::vector<ws> l;
  ws sum;
  for(TrkHitList::hot_iterator i(hotlist->begin());i!=hotlist->end();++i){
    double hottime,hottimeerr;
    if(!useHot(*i) || !i->timeResid(hottime,hottimeerr)) continue;
    hottime*=1e9; hottimeerr*=1e9;
    ws x(hottime,double(1)/(hottimeerr*hottimeerr));
    l.push_back(x);
    sum +=x;
  }
  if (_maxpull>0 && sum.n()>2 && sum.isPhysical()) { // do pull based pruning
    std::vector<ws>::iterator worst=l.end();
    do {
        worst=l.end(); double worstPull = _maxpull;
        for(std::vector<ws>::iterator i=l.begin();i!=l.end();++i) {
              ws u(sum); u-=*i;
              if (!u.isPhysical()) {
#ifdef MDCPATREC_WARNING
                std::cout<<"ErrMsg(warning) rounding problem?" << std::endl;
#endif
                      continue;
              }
              double p = fabs(u.pull( *i ));
              if (p>worstPull) { worstPull = p; worst = i; }
        }
        if (worst!=l.end()) {
                sum-=*worst; l.erase(worst);
        }
    } while (sum.isPhysical() && sum.n()>2 && worst!=l.end());
  }
  if (sum.isPhysical()) {
          time = sum.mean()*1e-9 + trk.trackT0();;
          timeerr = sum.sigma()*1e-9;
          nHotsUsed = sum.n();
          return true;
  }
  return false;
}

---

The documentation for this class was generated from the following files:

• TrkOutAvgTimeCalculator.h
• TrkOutAvgTimeCalculator.cxx