TrkOutAvgTimeCalculator.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: TrkOutAvgTimeCalculator.cxx,v 1.2 2005/12/01 06:18:42 zhangy Exp $
//
// Description:
//   See header file
//
// Environment:
//      Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Copyright Information:
//      Copyright (C) 2001      UC, San Diego
//
// Author List:
//      Gerhard Raven 7/19/01
//------------------------------------------------------------------------
#include "TrkBase/TrkOutAvgTimeCalculator.h"
#include "TrkBase/TrkHotSelector.h"
#include "TrkBase/TrkRecoTrk.h"
#include "TrkBase/TrkHitOnTrk.h"
#include "TrkBase/TrkHitList.h"
#include <vector>
 
 
TrkOutAvgTimeCalculator::TrkOutAvgTimeCalculator(const TrkHotSelector& selector, double maxpull) :
  TrkTimeCalculator(selector),_maxpull(maxpull)
{}
 
TrkOutAvgTimeCalculator::~TrkOutAvgTimeCalculator()
{}
 
bool
TrkOutAvgTimeCalculator::trackTime(const TrkRecoTrk& trk,
                                   double& time, double& timeerr,
                                   int& nHotsUsed) const
{
  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;
}