trkf::Track3DKalmanHitAlg Class Reference

#include "Track3DKalmanHitAlg.h"

Public Member Functions
	Track3DKalmanHitAlg (const fhicl::ParameterSet &pset)
	Constructor. More...
std::vector< trkf::KalmanOutput >	makeTracks (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, KalmanInputs &kalman_inputs)
void	fetchPFParticleSeeds (const art::PtrVector< recob::Seed > &pfseeds, const std::vector< Hits > &pfseedhits, std::vector< recob::Seed > &seeds, std::vector< Hits > &hitsperseed) const
	Fetch Seeds method. More...
recob::Seed	makeSeed (detinfo::DetectorPropertiesData const &detProp, const Hits &hits) const
	Make seed method. More...
void	growSeedsIntoTracks (detinfo::DetectorPropertiesData const &detProp, const bool pfseed, const std::vector< recob::Seed > &seeds, const std::vector< Hits > &hitsperseed, Hits &unusedhits, Hits &hits, std::deque< KGTrack > &kalman_tracks)
	Grow Seeds method. More...
void	growSeedIntoTracks (detinfo::DetectorPropertiesData const &detProp, const bool pfseed, const recob::Seed &seed, const Hits &hpsit, Hits &unusedhits, Hits &hits, std::deque< KGTrack > &kgtracks)
void	chopHitsOffSeeds (Hits const &hpsit, bool pfseed, Hits &seedhits) const
	Chop hits off of the end of seeds. More...
bool	testSeedSlope (const double *dir) const
std::shared_ptr< Surface >	makeSurface (const recob::Seed &seed, double *dir) const
	method to return a seed to surface. More...
bool	makeKalmanTracks (detinfo::DetectorPropertiesData const &detProp, const std::shared_ptr< trkf::Surface > psurf, const Surface::TrackDirection trkdir, Hits &seedhits, Hits &hits, std::deque< KGTrack > &kalman_tracks)
bool	smoothandextendTrack (detinfo::DetectorPropertiesData const &detProp, Propagator const &propagator, KGTrack &trg0, const Hits hits, unsigned int prefplane, std::deque< KGTrack > &kalman_tracks)
	SMooth and extend track. More...
bool	extendandsmoothLoop (detinfo::DetectorPropertiesData const &detProp, Propagator const &propagator, KGTrack &trg1, unsigned int prefplane, Hits &trackhits) const
	SMooth and extend a track in a loop. More...
void	filterHitsOnKalmanTrack (const KGTrack &trg, Hits &hits, Hits &seederhits) const
	Filter hits that are on kalman tracks. More...
std::unique_ptr< KHitContainer >	fillHitContainer (detinfo::DetectorPropertiesData const &detProp, const Hits &hits) const
	Fill hit container with either seedhits or filtered hits i.e. recob::Hit. More...
bool	qualityCutsOnSeedTrack (const KGTrack &trg0) const
	Quality cuts on seed track. More...
void	fitnupdateMomentum (Propagator const &propagator, KGTrack &trg1, KGTrack &trg2) const
	fit and update method, used twice. More...

Private Attributes
bool	fDoDedx
	Global dE/dx enable flag. More...
bool	fSelfSeed
	Self seed flag. More...
double	fMaxTcut
	Maximum delta ray energy in MeV for restricted dE/dx. More...
bool	fLineSurface
	Line surface flag. More...
size_t	fMinSeedHits
	Minimum number of hits per track seed. More...
int	fMinSeedChopHits
	Potentially chop seeds that exceed this length. More...
int	fMaxChopHits
	Maximum number of hits to chop from each end of seed. More...
double	fMaxSeedChiDF
	Maximum seed track chisquare/dof. More...
double	fMinSeedSlope
	Minimum seed slope (dx/dz). More...
double	fInitialMomentum
	Initial (or constant) momentum. More...
KalmanFilterAlg	fKFAlg
	Kalman filter algorithm. More...
SeedFinderAlgorithm	fSeedFinderAlg
	Seed finder. More...
int	fNumTrack
	Number of tracks produced. More...

Detailed Description

Definition at line 60 of file Track3DKalmanHitAlg.h.

Constructor & Destructor Documentation

trkf::Track3DKalmanHitAlg::Track3DKalmanHitAlg ( const fhicl::ParameterSet &  pset )

explicit

Constructor.

Definition at line 82 of file Track3DKalmanHitAlg.cxx.

References fInitialMomentum, fKFAlg, fLineSurface, fMaxChopHits, fMaxSeedChiDF, fMaxTcut, fMinSeedChopHits, fMinSeedHits, fMinSeedSlope, fNumTrack, fSeedFinderAlg, fSelfSeed, and fhicl::ParameterSet::get().

  : fDoDedx{pset.get<bool>("DoDedx")}
  , fSelfSeed{pset.get<bool>("SelfSeed")}
  , fMaxTcut{pset.get<double>("MaxTcut")}
  , fLineSurface{pset.get<bool>("LineSurface")}
  , fMinSeedHits{pset.get<size_t>("MinSeedHits")}
  , fMinSeedChopHits{pset.get<int>("MinSeedChopHits")}
  , fMaxChopHits{pset.get<int>("MaxChopHits")}
  , fMaxSeedChiDF{pset.get<double>("MaxSeedChiDF")}
  , fMinSeedSlope{pset.get<double>("MinSeedSlope")}
  , fInitialMomentum{pset.get<double>("InitialMomentum")}
  , fKFAlg(pset.get<fhicl::ParameterSet>("KalmanFilterAlg"))
  , fSeedFinderAlg(pset.get<fhicl::ParameterSet>("SeedFinderAlg"))
  , fNumTrack(0)
{
  mf::LogInfo("Track3DKalmanHitAlg") << "Track3DKalmanHitAlg instantiated.";
}

Member Function Documentation

void trkf::Track3DKalmanHitAlg::chopHitsOffSeeds	(	Hits const &	hpsit,
		bool	pfseed,
		Hits &	seedhits
	)		const

Chop hits off of the end of seeds.

Definition at line 379 of file Track3DKalmanHitAlg.cxx.

References art::PtrVector< T >::begin(), art::PtrVector< T >::end(), fMaxChopHits, fMinSeedChopHits, fSelfSeed, art::PtrVector< T >::push_back(), art::PtrVector< T >::reserve(), and art::PtrVector< T >::size().

Referenced by growSeedIntoTracks().

{
  const int nchopmax =
    (pfseed || fSelfSeed) ? 0 : std::max(0, int((hpsit.size() - fMinSeedChopHits) / 2));
  //SS: FIXME

  const int nchop = std::min(nchopmax, fMaxChopHits);
  Hits::const_iterator itb = hpsit.begin();
  Hits::const_iterator ite = hpsit.end();
  itb += nchop;
  ite -= nchop;
  seedhits.reserve(hpsit.size());
  for (Hits::const_iterator it = itb; it != ite; ++it)
    seedhits.push_back(*it);
}

bool trkf::Track3DKalmanHitAlg::extendandsmoothLoop	(	detinfo::DetectorPropertiesData const &	detProp,
		Propagator const &	propagator,
		KGTrack &	trg1,
		unsigned int	prefplane,
		Hits &	trackhits
	)		const

SMooth and extend a track in a loop.

Definition at line 449 of file Track3DKalmanHitAlg.cxx.

References trkf::KalmanFilterAlg::extendTrack(), fDoDedx, fillHitContainer(), trkf::KGTrack::fillHits(), fitnupdateMomentum(), fKFAlg, and trkf::KalmanFilterAlg::smoothTrack().

Referenced by smoothandextendTrack().

{
  bool ok = true;
  const int niter = 6;
  int nfail = 0; // Number of consecutive extend fails.
  for (int ix = 0; ok && ix < niter && nfail < 2; ++ix) {

    // Fill a collection of hits from the last good track
    // (initially the seed track).
    Hits goodhits;
    std::vector<unsigned int> hittpindex;
    trg1.fillHits(goodhits, hittpindex);

    // Filter hits already on the track out of the available hits.
    filterHits(trackhits, goodhits);

    // Fill hit container using filtered hits.
    std::unique_ptr<KHitContainer> ptrackcont = fillHitContainer(detProp, trackhits);

    // Extend the track.  It is not an error for the
    // extend operation to fail, meaning that no new hits
    // were added.

    if (fKFAlg.extendTrack(trg1, propagator, *ptrackcont))
      nfail = 0;
    else
      ++nfail;

    // Smooth the extended track, and make a new
    // unidirectionally fit track in the opposite
    // direction.

    KGTrack trg2(prefplane);
    ok = fKFAlg.smoothTrack(trg1, &trg2, propagator);
    if (ok) {
      // Skip momentum estimate for constant-momentum tracks.
      if (fDoDedx) { fitnupdateMomentum(propagator, trg1, trg2); }
      trg1 = trg2;
    }
  }
  return ok;
}

void trkf::Track3DKalmanHitAlg::fetchPFParticleSeeds	(	const art::PtrVector< recob::Seed > &	pfseeds,
		const std::vector< Hits > &	pfseedhits,
		std::vector< recob::Seed > &	seeds,
		std::vector< Hits > &	hitsperseed
	)		const

Fetch Seeds method.

Definition at line 165 of file Track3DKalmanHitAlg.cxx.

Referenced by makeTracks().

{
  for (const auto& pseed : pfseeds) {
    seeds.push_back(*pseed);
  }
  hitsperseed.insert(hitsperseed.end(), pfseedhits.begin(), pfseedhits.end());
}

std::unique_ptr< trkf::KHitContainer > trkf::Track3DKalmanHitAlg::fillHitContainer	(	detinfo::DetectorPropertiesData const &	detProp,
		const Hits &	hits
	)		const

Fill hit container with either seedhits or filtered hits i.e. recob::Hit.

Definition at line 344 of file Track3DKalmanHitAlg.cxx.

References fLineSurface.

Referenced by extendandsmoothLoop(), and makeKalmanTracks().

{
  std::unique_ptr<KHitContainer> hitcont(fLineSurface ?
                                           static_cast<KHitContainer*>(new KHitContainerWireLine) :
                                           static_cast<KHitContainer*>(new KHitContainerWireX));
  hitcont->fill(detProp, hits, -1);
  return hitcont;
}

void trkf::Track3DKalmanHitAlg::filterHitsOnKalmanTrack	(	const KGTrack &	trg,
		Hits &	hits,
		Hits &	seederhits
	)		const

Filter hits that are on kalman tracks.

Definition at line 331 of file Track3DKalmanHitAlg.cxx.

References trkf::KGTrack::fillHits().

Referenced by growSeedIntoTracks().

{
  Hits track_used_hits;
  std::vector<unsigned int> hittpindex;
  trg.fillHits(track_used_hits, hittpindex);
  filterHits(hits, track_used_hits);
  filterHits(seederhits, track_used_hits);
}

void trkf::Track3DKalmanHitAlg::fitnupdateMomentum	(	Propagator const &	propagator,
		KGTrack &	trg1,
		KGTrack &	trg2
	)		const

fit and update method, used twice.

Definition at line 497 of file Track3DKalmanHitAlg.cxx.

References trkf::KalmanFilterAlg::fitMomentum(), fKFAlg, and trkf::KalmanFilterAlg::updateMomentum().

Referenced by extendandsmoothLoop(), and smoothandextendTrack().

{
  KETrack tremom;
  if (fKFAlg.fitMomentum(trg1, propagator, tremom)) {
    fKFAlg.updateMomentum(tremom, propagator, trg2);
  }
}

void trkf::Track3DKalmanHitAlg::growSeedIntoTracks	(	detinfo::DetectorPropertiesData const &	detProp,
		const bool	pfseed,
		const recob::Seed &	seed,
		const Hits &	hpsit,
		Hits &	unusedhits,
		Hits &	hits,
		std::deque< KGTrack > &	kgtracks
	)

Definition at line 198 of file Track3DKalmanHitAlg.cxx.

References trkf::Surface::BACKWARD, chopHitsOffSeeds(), dir, fDoDedx, filterHitsOnKalmanTrack(), trkf::Surface::FORWARD, makeKalmanTracks(), makeSurface(), art::PtrVector< T >::size(), and testSeedSlope().

Referenced by growSeedsIntoTracks().

{
  Hits trimmedhits;
  // Chop a couple of hits off each end of the seed.
  chopHitsOffSeeds(hpsit, pfseed, trimmedhits);

  // Filter hits used by (chopped) seed from hits available to make future seeds.
  // No matter what, we will never use these hits for another seed.
  // This eliminates the possibility of an infinite loop.

  size_t initial_unusedhits = unusedhits.size();
  filterHits(unusedhits, trimmedhits);

  // Require that this seed be fully disjoint from existing tracks.
  //SS: replace this test with a method with appropriate name
  if (!(trimmedhits.size() + unusedhits.size() == initial_unusedhits)) return;

  // Convert seed into initial KTracks on surface located at seed point,
  // and normal to seed direction.
  double dir[3];
  std::shared_ptr<Surface> psurf = makeSurface(seed, dir);

  // Cut on the seed slope dx/ds.
  //SS: replace test name with a reasonable name
  if (!testSeedSlope(dir)) return;

  // Make one or two initial KTracks for forward and backward directions.
  // The build_both flag specifies whether we should attempt to make
  // tracks from all both FORWARD and BACKWARD initial tracks,
  // or alternatively, whether we
  // should declare victory and quit after getting a successful
  // track from one initial track.
  const bool build_both = fDoDedx;
  const int ninit = 2;

  auto ntracks = kgtracks.size(); // Remember original track count.
  bool ok = makeKalmanTracks(detProp, psurf, Surface::FORWARD, trimmedhits, hits, kgtracks);
  if ((!ok || build_both) && ninit == 2) {
    makeKalmanTracks(detProp, psurf, Surface::BACKWARD, trimmedhits, hits, kgtracks);
  }

  // Loop over newly added tracks and remove hits contained on
  // these tracks from hits available for making additional
  // tracks or track seeds.
  for (unsigned int itrk = ntracks; itrk < kgtracks.size(); ++itrk) {
    const KGTrack& trg = kgtracks[itrk];
    filterHitsOnKalmanTrack(trg, hits, unusedhits);
  }
}

void trkf::Track3DKalmanHitAlg::growSeedsIntoTracks	(	detinfo::DetectorPropertiesData const &	detProp,
		const bool	pfseed,
		const std::vector< recob::Seed > &	seeds,
		const std::vector< Hits > &	hitsperseed,
		Hits &	unusedhits,
		Hits &	hits,
		std::deque< KGTrack > &	kalman_tracks
	)

Grow Seeds method.

Definition at line 179 of file Track3DKalmanHitAlg.cxx.

References growSeedIntoTracks().

Referenced by makeTracks().

{
  // check for size of both containers
  if (seeds.size() != hitsperseed.size()) {
    throw cet::exception("Track3DKalmanHitAlg")
      << "Different size containers for Seeds and Hits/Seed.\n";
  }
  for (size_t i = 0; i < seeds.size(); ++i) {
    growSeedIntoTracks(detProp, pfseed, seeds[i], hitsperseed[i], unusedhits, hits, kgtracks);
  }
}

bool trkf::Track3DKalmanHitAlg::makeKalmanTracks	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::shared_ptr< trkf::Surface >	psurf,
		const Surface::TrackDirection	trkdir,
		Hits &	seedhits,
		Hits &	hits,
		std::deque< KGTrack > &	kalman_tracks
	)

Definition at line 276 of file Track3DKalmanHitAlg.cxx.

References trkf::KalmanFilterAlg::buildTrack(), fDoDedx, fillHitContainer(), fInitialMomentum, fKFAlg, fMaxTcut, trkf::Propagator::FORWARD, fSelfSeed, mf::isDebugEnabled(), trkf::KalmanFilterAlg::setPlane(), and smoothandextendTrack().

Referenced by growSeedIntoTracks().

{
  const int pdg = 13; //SS: FIXME another constant?
  // SS: FIXME
  // It is a contant value inside a loop so I took it out
  // revisit the linear algebra stuff used here (use of ublas)
  // make a lambda here ... const TrackVector
  //
  TrackVector vec(5);
  vec(0) = 0.;
  vec(1) = 0.;
  vec(2) = 0.;
  vec(3) = 0.;
  vec(4) = (fInitialMomentum != 0. ? 1. / fInitialMomentum : 2.);

  KTrack initial_track(psurf, vec, trkdir, pdg);

  // Fill hit container with current seed hits.
  // KHitContainer may not be cheaply movabale thats why unique_ptr
  std::unique_ptr<KHitContainer> pseedcont = fillHitContainer(detProp, seedhits);

  // Set the preferred plane to be the one with the most hits.
  unsigned int prefplane = pseedcont->getPreferredPlane();
  fKFAlg.setPlane(prefplane);
  if (mf::isDebugEnabled())
    mf::LogDebug("Track3DKalmanHit") << "Preferred plane = " << prefplane << "\n";

  PropAny const propagator{detProp, fMaxTcut, fDoDedx};

  // Build and smooth seed track.
  KGTrack trg0(prefplane);
  bool ok =
    fKFAlg.buildTrack(initial_track, trg0, propagator, Propagator::FORWARD, *pseedcont, fSelfSeed);
  if (ok) ok = smoothandextendTrack(detProp, propagator, trg0, hits, prefplane, kgtracks);

  if (mf::isDebugEnabled())
    mf::LogDebug("Track3DKalmanHit")
      << (ok ? "Find track succeeded." : "Find track failed.") << "\n";
  return ok;
}

recob::Seed trkf::Track3DKalmanHitAlg::makeSeed	(	detinfo::DetectorPropertiesData const &	detProp,
		const Hits &	hits
	)		const

Make seed method.

Definition at line 509 of file Track3DKalmanHitAlg.cxx.

References detinfo::DetectorPropertiesData::ConvertTicksToX(), dir, geo::WireGeo::GetCenter(), n, recob::Hit::PeakTime(), trkf::syminvert(), geo::WireGeo::ThetaZ(), w, geo::GeometryCore::Wire(), recob::Hit::WireID(), and x.

Referenced by makeTracks().

{
  art::ServiceHandle<geo::Geometry const> geom;

  // Do a linear 3D least squares for of y and z vs. x.
  // y = y0 + ay*(x-x0)
  // z = z0 + az*(x-x0)
  // Here x0 is the global average x coordinate of all hits in all planes.
  // Parameters y0, z0, ay, and az are determined by a simultaneous least squares
  // fit in all planes.

  // First, determine x0 by looping over every hit.

  double x0 = 0.;
  int n = 0;
  for (auto const& phit : hits) {
    const recob::Hit& hit = *phit;
    geo::WireID wire_id = hit.WireID();
    double time = hit.PeakTime();
    double x = detProp.ConvertTicksToX(time, wire_id);
    x0 += x;
    ++n;
  }

  // If there are no hits, return invalid seed.

  if (n == 0) return recob::Seed();

  // Find the average x.

  x0 /= n;

  // Now do the least squares fit proper.

  KSymMatrix<4>::type sm(4);
  KVector<4>::type sv(4);
  sm.clear();
  sv.clear();

  // Loop over hits (again).

  for (auto const& phit : hits) {
    const recob::Hit& hit = *phit;

    // Extract the angle, w and x coordinates from hit.

    geo::WireID wire_id = hit.WireID();
    const geo::WireGeo& wgeom = geom->Wire(wire_id);
    auto const xyz = wgeom.GetCenter();

    // Phi convention is the one documented in SurfYZPlane.h.

    double phi = TMath::PiOver2() - wgeom.ThetaZ();
    double sphi = std::sin(phi);
    double cphi = std::cos(phi);
    double w = -xyz.Y() * sphi + xyz.Z() * cphi;

    double time = hit.PeakTime();
    double x = detProp.ConvertTicksToX(time, wire_id);

    // Accumulate data for least squares fit.

    double dx = x - x0;

    sm(0, 0) += sphi * sphi;
    sm(1, 0) -= sphi * cphi;
    sm(1, 1) += cphi * cphi;
    sm(2, 0) += sphi * sphi * dx;
    sm(2, 1) -= sphi * cphi * dx;
    sm(2, 2) += sphi * sphi * dx * dx;
    sm(3, 0) -= sphi * cphi * dx;
    sm(3, 1) += cphi * cphi * dx;
    sm(3, 2) -= sphi * cphi * dx * dx;
    sm(3, 3) += cphi * cphi * dx * dx;

    sv(0) -= sphi * w;
    sv(1) += cphi * w;
    sv(2) -= sphi * w * dx;
    sv(3) += cphi * w * dx;
  }

  // Solve.

  bool ok = syminvert(sm);
  if (!ok) return recob::Seed();
  KVector<4>::type par(4);
  par = prod(sm, sv);

  double y0 = par(0);
  double z0 = par(1);
  double dydx = par(2);
  double dzdx = par(3);

  // Make seed.

  double dsdx = std::hypot(1., std::hypot(dydx, dzdx));

  double pos[3] = {x0, y0, z0};
  double dir[3] = {1. / dsdx, dydx / dsdx, dzdx / dsdx};
  double poserr[3] = {0., 0., 0.};
  double direrr[3] = {0., 0., 0.};

  recob::Seed result(pos, dir, poserr, direrr);

  return result;
}

std::shared_ptr< trkf::Surface > trkf::Track3DKalmanHitAlg::makeSurface	(	const recob::Seed &	seed,
		double *	dir
	)		const

method to return a seed to surface.

Definition at line 257 of file Track3DKalmanHitAlg.cxx.

References recob::Seed::GetDirection(), recob::Seed::GetPoint(), and mf::isDebugEnabled().

Referenced by growSeedIntoTracks().

{
  double xyz[3];
  double err[3]; // Dummy.
  seed.GetPoint(xyz, err);
  seed.GetDirection(dir, err);
  if (mf::isDebugEnabled()) {
    mf::LogDebug("Track3DKalmanHit")
      //<< "Seed found with " << seedhits.size() <<" hits.\n"
      << "(x,y,z) = " << xyz[0] << ", " << xyz[1] << ", " << xyz[2] << "\n"
      << "(dx,dy,dz) = " << dir[0] << ", " << dir[1] << ", " << dir[2] << "\n";
  } // if debug

  return std::make_shared<SurfXYZPlane>(xyz[0], xyz[1], xyz[2], dir[0], dir[1], dir[2]);
}

std::vector< trkf::KalmanOutput > trkf::Track3DKalmanHitAlg::makeTracks	(	detinfo::DetectorClocksData const &	clockData,
		detinfo::DetectorPropertiesData const &	detProp,
		KalmanInputs &	kalman_inputs
	)

Definition at line 102 of file Track3DKalmanHitAlg.cxx.

References art::PtrVector< T >::begin(), util::empty(), art::PtrVector< T >::end(), fetchPFParticleSeeds(), fSeedFinderAlg, fSelfSeed, trkf::SeedFinderAlgorithm::GetSeedsFromUnSortedHits(), growSeedsIntoTracks(), hits(), makeSeed(), tca::seeds, and art::PtrVector< T >::size().

Referenced by trkf::Track3DKalmanHit::produce().

{
  std::vector<KalmanOutput> outputs(inputs.size());
  // Loop over input vectors.
  for (size_t i = 0; i < inputs.size(); ++i) {
    Hits& hits = inputs[i].hits;
    // The hit collection "hits" (just filled), initially containing all
    // hits, represents hits available for making tracks.  Now we will
    // fill a second hit collection called "unusedhits", also initially
    // containing all hits, which will represent hits available for
    // making track seeds.  These collections are not necessarily the
    // same, since hits that are not suitable for seeds may still be
    // suitable for tracks.
    Hits unusedhits = hits;
    // Start of loop.
    bool first = true;
    while (true) {
      std::vector<Hits> hitsperseed;
      std::vector<recob::Seed> seeds;
      // On the first trip through this loop, try to use pfparticle-associated seeds.
      // Do this, provided the list of pfparticle-associated seeds and associated
      // hits are not empty.
      bool pfseed = false;
      auto const seedsize = inputs[i].seeds.size();
      if (first && seedsize > 0 && inputs[i].seedhits.size() > 0) {
        seeds.reserve(seedsize);
        fetchPFParticleSeeds(inputs[i].seeds, inputs[i].seedhits, seeds, hitsperseed);
        pfseed = true;
      }
      else {
        // On subsequent trips, or if there were no usable pfparticle-associated seeds,
        // attempt to generate our own seeds.
        if (unusedhits.size() > 0) {
          if (fSelfSeed) {
            // Self seed - convert all hits into one big seed.
            seeds.emplace_back(makeSeed(detProp, unusedhits));
            hitsperseed.emplace_back();
            hitsperseed.back().insert(
              hitsperseed.back().end(), unusedhits.begin(), unusedhits.end());
          }
          else
            seeds =
              fSeedFinderAlg.GetSeedsFromUnSortedHits(clockData, detProp, unusedhits, hitsperseed);
        }
      }

      assert(seeds.size() == hitsperseed.size());
      first = false;

      if (empty(seeds)) break;

      growSeedsIntoTracks(detProp, pfseed, seeds, hitsperseed, unusedhits, hits, outputs[i].tracks);
    }
  }
  return outputs;
}

bool trkf::Track3DKalmanHitAlg::qualityCutsOnSeedTrack

(

const KGTrack &

trg0

)

const

Quality cuts on seed track.

Definition at line 359 of file Track3DKalmanHitAlg.cxx.

References util::abs(), trkf::KGTrack::endTrack(), fMinSeedSlope, trkf::KTrack::getMomentum(), and trkf::KGTrack::startTrack().

Referenced by smoothandextendTrack().

{
  double mom0[3];
  double mom1[3];
  trg0.startTrack().getMomentum(mom0);
  trg0.endTrack().getMomentum(mom1);

  double dxds0 = mom0[0] / calcMagnitude(mom0);
  double dxds1 = mom1[0] / calcMagnitude(mom1);

  return (std::abs(dxds0) > fMinSeedSlope && std::abs(dxds1) > fMinSeedSlope);
}

bool trkf::Track3DKalmanHitAlg::smoothandextendTrack	(	detinfo::DetectorPropertiesData const &	detProp,
		Propagator const &	propagator,
		KGTrack &	trg0,
		const Hits	hits,
		unsigned int	prefplane,
		std::deque< KGTrack > &	kalman_tracks
	)

SMooth and extend track.

Definition at line 399 of file Track3DKalmanHitAlg.cxx.

References trkf::KGTrack::endTrack(), extendandsmoothLoop(), fDoDedx, fitnupdateMomentum(), fKFAlg, fMaxSeedChiDF, fMinSeedHits, fNumTrack, trkf::KFitTrack::getChisq(), hits(), n, trkf::KGTrack::numHits(), qualityCutsOnSeedTrack(), trkf::KalmanFilterAlg::smoothTrack(), and trkf::KGTrack::startTrack().

Referenced by makeKalmanTracks().

{
  KGTrack trg1(prefplane);
  bool ok = fKFAlg.smoothTrack(trg0, &trg1, propagator);
  if (!ok) return ok;

  // Now we have the seed track in the form of a KGTrack.
  // Do additional quality cuts.

  auto const n = trg1.numHits();
  auto const chisq = n * fMaxSeedChiDF;

  ok = n >= fMinSeedHits && trg1.startTrack().getChisq() <= chisq &&
       trg1.endTrack().getChisq() <= chisq && trg0.startTrack().getChisq() <= chisq &&
       trg0.endTrack().getChisq() <= chisq && qualityCutsOnSeedTrack(trg0);

  if (!ok) return ok;

  // Make a copy of the original hit collection of all
  // available track hits.
  Hits trackhits = hits;

  // Do an extend + smooth loop here.
  // Exit after two consecutive failures to extend (i.e. from each end),
  // or if the iteration count reaches the maximum.
  if (ok) ok = extendandsmoothLoop(detProp, propagator, trg1, prefplane, trackhits);

  // Do a final smooth.
  if (!ok) return false;

  ok = fKFAlg.smoothTrack(trg1, 0, propagator);
  if (!ok) return false;

  // Skip momentum estimate for constant-momentum tracks.

  if (fDoDedx) { fitnupdateMomentum(propagator, trg1, trg1); }
  // Save this track.
  ++fNumTrack;
  kalman_tracks.push_back(trg1);
  return true;
}

bool trkf::Track3DKalmanHitAlg::testSeedSlope

(

const double *

dir

)

const

Definition at line 324 of file Track3DKalmanHitAlg.cxx.

References util::abs(), and fMinSeedSlope.

Referenced by growSeedIntoTracks().

{
  return std::abs(dir[0]) >= fMinSeedSlope * calcMagnitude(dir);
}

Member Data Documentation

bool trkf::Track3DKalmanHitAlg::fDoDedx

private

Global dE/dx enable flag.

Definition at line 116 of file Track3DKalmanHitAlg.h.

Referenced by extendandsmoothLoop(), growSeedIntoTracks(), makeKalmanTracks(), and smoothandextendTrack().

double trkf::Track3DKalmanHitAlg::fInitialMomentum

private

Initial (or constant) momentum.

Definition at line 125 of file Track3DKalmanHitAlg.h.

Referenced by makeKalmanTracks(), and Track3DKalmanHitAlg().

KalmanFilterAlg trkf::Track3DKalmanHitAlg::fKFAlg

private

Kalman filter algorithm.

Definition at line 129 of file Track3DKalmanHitAlg.h.

Referenced by extendandsmoothLoop(), fitnupdateMomentum(), makeKalmanTracks(), smoothandextendTrack(), and Track3DKalmanHitAlg().

bool trkf::Track3DKalmanHitAlg::fLineSurface

private

Line surface flag.

Definition at line 119 of file Track3DKalmanHitAlg.h.

Referenced by fillHitContainer(), and Track3DKalmanHitAlg().

int trkf::Track3DKalmanHitAlg::fMaxChopHits

private

Maximum number of hits to chop from each end of seed.

Definition at line 122 of file Track3DKalmanHitAlg.h.

Referenced by chopHitsOffSeeds(), and Track3DKalmanHitAlg().

double trkf::Track3DKalmanHitAlg::fMaxSeedChiDF

private

Maximum seed track chisquare/dof.

Definition at line 123 of file Track3DKalmanHitAlg.h.

Referenced by smoothandextendTrack(), and Track3DKalmanHitAlg().

double trkf::Track3DKalmanHitAlg::fMaxTcut

private

Maximum delta ray energy in MeV for restricted dE/dx.

Definition at line 118 of file Track3DKalmanHitAlg.h.

Referenced by makeKalmanTracks(), and Track3DKalmanHitAlg().

int trkf::Track3DKalmanHitAlg::fMinSeedChopHits

private

Potentially chop seeds that exceed this length.

Definition at line 121 of file Track3DKalmanHitAlg.h.

Referenced by chopHitsOffSeeds(), and Track3DKalmanHitAlg().

size_t trkf::Track3DKalmanHitAlg::fMinSeedHits

private

Minimum number of hits per track seed.

Definition at line 120 of file Track3DKalmanHitAlg.h.

Referenced by smoothandextendTrack(), and Track3DKalmanHitAlg().

double trkf::Track3DKalmanHitAlg::fMinSeedSlope

private

Minimum seed slope (dx/dz).

Definition at line 124 of file Track3DKalmanHitAlg.h.

Referenced by qualityCutsOnSeedTrack(), testSeedSlope(), and Track3DKalmanHitAlg().

int trkf::Track3DKalmanHitAlg::fNumTrack

private

Number of tracks produced.

Definition at line 133 of file Track3DKalmanHitAlg.h.

Referenced by smoothandextendTrack(), and Track3DKalmanHitAlg().

SeedFinderAlgorithm trkf::Track3DKalmanHitAlg::fSeedFinderAlg

private

Seed finder.

Definition at line 130 of file Track3DKalmanHitAlg.h.

Referenced by makeTracks(), and Track3DKalmanHitAlg().

bool trkf::Track3DKalmanHitAlg::fSelfSeed

private

Self seed flag.

Definition at line 117 of file Track3DKalmanHitAlg.h.

Referenced by chopHitsOffSeeds(), makeKalmanTracks(), makeTracks(), and Track3DKalmanHitAlg().

---

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

• larreco/v09_25_00/source/larreco/RecoAlg/Track3DKalmanHitAlg.h
• larreco/v09_25_00/source/larreco/RecoAlg/Track3DKalmanHitAlg.cxx