track finding procedure More...

#include <BKLMTrackFinder.h>

Public Member Functions
	BKLMTrackFinder ()
	Default constructor.

	BKLMTrackFinder (BKLMTrackFitter *fitter)
	The track finder requires a track fitter.

	~BKLMTrackFinder ()
	Destructor.

void	registerFitter (BKLMTrackFitter *fitter)
	Register a fitter if not constructed with one.

bool	filter (const std::list< KLMHit2d * > &seed, std::list< KLMHit2d * > &hits, std::list< KLMHit2d * > &track)
	find associated hits and do fit.

void	setGlobalFit (bool localOrGlobal)
	set the fitting mode, local system or global system

Protected Attributes
BKLMTrackFitter *	m_Fitter
	pointer to the fitter

bool	m_globalFit
	do fit in the local system or global system false: local sys; true: global sys.

Detailed Description

track finding procedure

Definition at line 25 of file BKLMTrackFinder.h.

Constructor & Destructor Documentation

◆ BKLMTrackFinder() [1/2]

BKLMTrackFinder ( )

Default constructor.

Definition at line 22 of file BKLMTrackFinder.cc.

23{

24}

◆ BKLMTrackFinder() [2/2]

BKLMTrackFinder ( BKLMTrackFitter * fitter )

explicit

The track finder requires a track fitter.

Definition at line 26 of file BKLMTrackFinder.cc.

                                                        :
  m_Fitter(fitter),
  m_globalFit(false)
{
}

◆ ~BKLMTrackFinder()

~BKLMTrackFinder ( )

Destructor.

Definition at line 33 of file BKLMTrackFinder.cc.

34{

35}

Member Function Documentation

◆ filter()

bool filter	(	const std::list< KLMHit2d * > &	seed,
		std::list< KLMHit2d * > &	hits,
		std::list< KLMHit2d * > &	track
	)

find associated hits and do fit.

find associated hits and do fit

Definition at line 45 of file BKLMTrackFinder.cc.

{
  track = seed;
 
  if (m_Fitter == 0) {
    B2ERROR("BKLMTrackFinder: Fitter not registered");
    return (false);
  }
 
  m_Fitter->fit(track);//fit seed
 
  for (KLMHit2d* hit2d : hits) {
    if (hit2d->getSubdetector() != KLMElementNumbers::c_BKLM)
      continue;
 
    // Prevent duplicate hits or hits on same layer
    // no duplicate hit is already guaranteed and now we allow hits on same layer so the following is commented out
    // bool skip = false;
    // for (j = track.begin(); j != track.end(); ++j) {
    //  if ((*j)->getLayer() == (*i)->getLayer())
    //    skip = true;
    // }
    // if (skip == true)
    //   continue;
 
    if (hit2d->isOnStaTrack() == false) {
      double error, sigma;
      if (m_globalFit)
        m_Fitter->globalDistanceToHit(hit2d, error, sigma);
      else m_Fitter->distanceToHit(hit2d, error, sigma);
      //B2INFO("BKLMTrackFinder" << " Error: " << error << " Sigma: " << sigma);
      if (sigma < 5.0) {
        track.push_back(hit2d);
      }
    }
  }
 
  if (track.size() < 3)
    return false;
 
  // Fit with new hits
  double chisqr = m_Fitter->fit(track);
  B2DEBUG(20, "BKLMTrackFinder:" << "ChiSqr: " << chisqr);
 
  // Do this the hard way to count each layer separately.
  std::list<int> hitLayers;
  for (KLMHit2d* hit2d : track) {
    hitLayers.push_back(hit2d->getLayer());
  }
  hitLayers.sort();
  hitLayers.unique();
 
  int layers = (*(--hitLayers.end()) - * (hitLayers.begin()));
  int noHits = hitLayers.size();
 
  if (noHits >= 4 && double(noHits) / double(layers) >= 0.75) {
    return true;
  } else {
    return false;
  }
}

◆ registerFitter()

void registerFitter ( BKLMTrackFitter * fitter )

Register a fitter if not constructed with one.

Definition at line 38 of file BKLMTrackFinder.cc.

{
  m_Fitter = fitter;
  m_Fitter->setGlobalFit(m_globalFit);
}

◆ setGlobalFit()

void setGlobalFit ( bool localOrGlobal )

inline

set the fitting mode, local system or global system

Definition at line 47 of file BKLMTrackFinder.h.

    {
      m_globalFit = localOrGlobal;
    }

Member Data Documentation

◆ m_Fitter

BKLMTrackFitter* m_Fitter

protected

pointer to the fitter

Definition at line 55 of file BKLMTrackFinder.h.

◆ m_globalFit

bool m_globalFit