trkf::Track3DKalmanHitAlg Class Reference

#include "Track3DKalmanHitAlg.h"

Public Member Functions
	Track3DKalmanHitAlg (const fhicl::ParameterSet &pset)
	Constructor. More...
void	reconfigure (const fhicl::ParameterSet &pset)
	Reconfigure method. More...
std::vector< trkf::KalmanOutput >	makeTracks (KalmanInputs &kalman_inputs)
bool	fetchPFParticleSeeds (const art::PtrVector< recob::Seed > &pfseeds, const std::vector< Hits > &pfseedhits, std::vector< recob::Seed > &seeds, std::vector< Hits > &hitsperseed)
	Fetch Seeds method. More...
recob::Seed	makeSeed (const Hits &hits) const
	Make seed method. More...
void	growSeedsIntoTracks (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 (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 (const std::shared_ptr< trkf::Surface > psurf, const Surface::TrackDirection trkdir, Hits &seedhits, Hits &hits, std::deque< KGTrack > &kalman_tracks)
bool	smoothandextendTrack (KGTrack &trg0, const Hits hits, unsigned int prefplane, std::deque< KGTrack > &kalman_tracks)
	SMooth and extend track. More...
bool	extendandsmoothLoop (KGTrack &trg1, unsigned int prefplane, Hits &trackhits)
	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 (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 (KGTrack &trg1, KGTrack &trg2)
	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...
std::unique_ptr< const Propagator >	fProp
	Propagator. More...
int	fNumTrack
	Number of tracks produced. More...

Detailed Description

Definition at line 64 of file Track3DKalmanHitAlg.h.

Constructor & Destructor Documentation

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

explicit

Constructor.

Definition at line 52 of file Track3DKalmanHitAlg.cxx.

References reconfigure().

                                                                          :
fDoDedx(false),
fSelfSeed(false),
fMaxTcut(0.),
fLineSurface(false),
fMinSeedHits(0),
fMinSeedChopHits(0),
fMaxChopHits(0),
fMaxSeedChiDF(0.),
fMinSeedSlope(0.),
fInitialMomentum(0.),
fKFAlg(pset.get<fhicl::ParameterSet>("KalmanFilterAlg")),
fSeedFinderAlg(pset.get<fhicl::ParameterSet>("SeedFinderAlg")),
fProp(nullptr),
fNumTrack(0)
{
   mf::LogInfo("Track3DKalmanHitAlg") << "Track3DKalmanHitAlg instantiated.";
   
   // Load fcl parameters.
   
   reconfigure(pset);
}

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 364 of file Track3DKalmanHitAlg.cxx.

References art::PtrVector< T >::begin(), art::PtrVector< T >::end(), fMaxChopHits, fMinSeedChopHits, fSelfSeed, max, min, 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	(	KGTrack &	trg1,
		unsigned int	prefplane,
		Hits &	trackhits
	)

SMooth and extend a track in a loop.

Definition at line 437 of file Track3DKalmanHitAlg.cxx.

References trkf::KalmanFilterAlg::extendTrack(), fDoDedx, fillHitContainer(), trkf::KGTrack::fillHits(), fitnupdateMomentum(), fKFAlg, fProp, 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(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, fProp.get(), *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, fProp.get());
      if(ok) {
         // Skip momentum estimate for constant-momentum tracks.
         if(fDoDedx) {
            fitnupdateMomentum(trg1, trg2);
         }
         trg1 = trg2;
      }
   }
   return ok;
}

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

Fetch Seeds method.

Definition at line 156 of file Track3DKalmanHitAlg.cxx.

Referenced by makeTracks().

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

std::unique_ptr< trkf::KHitContainer > trkf::Track3DKalmanHitAlg::fillHitContainer

(

const Hits &

hits

)

const

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

Definition at line 332 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(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 320 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	(	KGTrack &	trg1,
		KGTrack &	trg2
	)

fit and update method, used twice.

Definition at line 481 of file Track3DKalmanHitAlg.cxx.

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

Referenced by extendandsmoothLoop(), and smoothandextendTrack().

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

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

Definition at line 189 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(psurf, Surface::FORWARD, trimmedhits, hits, kgtracks);
   if ((!ok || build_both) && ninit == 2) {
      makeKalmanTracks(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	(	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 172 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(pfseed, seeds[i], hitsperseed[i], unusedhits, hits, kgtracks);
   }
}

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

Definition at line 270 of file Track3DKalmanHitAlg.cxx.

References trkf::KalmanFilterAlg::buildTrack(), fillHitContainer(), fInitialMomentum, fKFAlg, trkf::Propagator::FORWARD, fProp, 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(KTrack(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(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";
   
   // Build and smooth seed track.
   KGTrack trg0(prefplane);
   bool ok = fKFAlg.buildTrack(initial_track, trg0, fProp.get(), Propagator::FORWARD, *pseedcont, fSelfSeed);
   if(ok) ok = smoothandextendTrack(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

(

const Hits &

hits

)

const

Make seed method.

Definition at line 490 of file Track3DKalmanHitAlg.cxx.

References detinfo::DetectorProperties::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().

{
   
   // Get Services.
   
   art::ServiceHandle<geo::Geometry> geom;
   const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
   
   // 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);
      double xyz[3];
      wgeom.GetCenter(xyz);
      
      // 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[1]*sphi + xyz[2]*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 249 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::shared_ptr<Surface>(new SurfXYZPlane(xyz[0], xyz[1], xyz[2],
                                                    dir[0], dir[1], dir[2]));
}

std::vector< trkf::KalmanOutput > trkf::Track3DKalmanHitAlg::makeTracks

(

KalmanInputs &

kalman_inputs

)

Definition at line 97 of file Track3DKalmanHitAlg.cxx.

References art::PtrVector< T >::begin(), 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);
            pfseed = fetchPFParticleSeeds(inputs[i].seeds,
                                          inputs[i].seedhits,
                                          seeds,
                                          hitsperseed);
         }
         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(unusedhits));
                  hitsperseed.emplace_back();
                  hitsperseed.back().insert(hitsperseed.back().end(),
                                            unusedhits.begin(),
                                            unusedhits.end());
               }
               else seeds = fSeedFinderAlg.GetSeedsFromUnSortedHits(unusedhits, hitsperseed);
            }
         }
         
         assert(seeds.size() == hitsperseed.size());
         first = false;
         
         if(seeds.size() == 0) break;
         else growSeedsIntoTracks(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 344 of file Track3DKalmanHitAlg.cxx.

References 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);
}

void trkf::Track3DKalmanHitAlg::reconfigure

(

const fhicl::ParameterSet &

pset

)

Reconfigure method.

Definition at line 78 of file Track3DKalmanHitAlg.cxx.

References fDoDedx, fInitialMomentum, fKFAlg, fLineSurface, fMaxChopHits, fMaxSeedChiDF, fMaxTcut, fMinSeedChopHits, fMinSeedHits, fMinSeedSlope, fProp, fSeedFinderAlg, fSelfSeed, fhicl::ParameterSet::get(), trkf::SeedFinderAlgorithm::reconfigure(), and trkf::KalmanFilterAlg::reconfigure().

Referenced by trkf::Track3DKalmanHit::reconfigure(), and Track3DKalmanHitAlg().

{
   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.reconfigure(pset.get<fhicl::ParameterSet>("KalmanFilterAlg"));
   fSeedFinderAlg.reconfigure(pset.get<fhicl::ParameterSet>("SeedFinderAlg"));
   fProp.reset(new PropAny(fMaxTcut, fDoDedx));
}

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

SMooth and extend track.

Definition at line 384 of file Track3DKalmanHitAlg.cxx.

References trkf::KGTrack::endTrack(), extendandsmoothLoop(), fDoDedx, fitnupdateMomentum(), fKFAlg, fMaxSeedChiDF, fMinSeedHits, fNumTrack, fProp, 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, fProp.get());
   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(trg1, prefplane, trackhits);
   
   // Do a final smooth.
   if(!ok) return ok;
   
   ok = fKFAlg.smoothTrack(trg1, 0, fProp.get());
   if(!ok) return ok;
   
   // Skip momentum estimate for constant-momentum tracks.
   
   if(fDoDedx) {
      fitnupdateMomentum(trg1, trg1);
   }
   // Save this track.
   ++fNumTrack;
   kalman_tracks.push_back(trg1);
   return ok;
}

bool trkf::Track3DKalmanHitAlg::testSeedSlope

(

const double *

dir

)

const

Definition at line 313 of file Track3DKalmanHitAlg.cxx.

References 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 124 of file Track3DKalmanHitAlg.h.

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

double trkf::Track3DKalmanHitAlg::fInitialMomentum

private

Initial (or constant) momentum.

Definition at line 133 of file Track3DKalmanHitAlg.h.

Referenced by makeKalmanTracks(), and reconfigure().

KalmanFilterAlg trkf::Track3DKalmanHitAlg::fKFAlg

private

Kalman filter algorithm.

Definition at line 137 of file Track3DKalmanHitAlg.h.

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

bool trkf::Track3DKalmanHitAlg::fLineSurface

private

Line surface flag.

Definition at line 127 of file Track3DKalmanHitAlg.h.

Referenced by fillHitContainer(), and reconfigure().

int trkf::Track3DKalmanHitAlg::fMaxChopHits

private

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

Definition at line 130 of file Track3DKalmanHitAlg.h.

Referenced by chopHitsOffSeeds(), and reconfigure().

double trkf::Track3DKalmanHitAlg::fMaxSeedChiDF

private

Maximum seed track chisquare/dof.

Definition at line 131 of file Track3DKalmanHitAlg.h.

Referenced by reconfigure(), and smoothandextendTrack().

double trkf::Track3DKalmanHitAlg::fMaxTcut

private

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

Definition at line 126 of file Track3DKalmanHitAlg.h.

Referenced by reconfigure().

int trkf::Track3DKalmanHitAlg::fMinSeedChopHits

private

Potentially chop seeds that exceed this length.

Definition at line 129 of file Track3DKalmanHitAlg.h.

Referenced by chopHitsOffSeeds(), and reconfigure().

size_t trkf::Track3DKalmanHitAlg::fMinSeedHits

private

Minimum number of hits per track seed.

Definition at line 128 of file Track3DKalmanHitAlg.h.

Referenced by reconfigure(), and smoothandextendTrack().

double trkf::Track3DKalmanHitAlg::fMinSeedSlope

private

Minimum seed slope (dx/dz).

Definition at line 132 of file Track3DKalmanHitAlg.h.

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

int trkf::Track3DKalmanHitAlg::fNumTrack

private

Number of tracks produced.

Definition at line 144 of file Track3DKalmanHitAlg.h.

Referenced by smoothandextendTrack().

std::unique_ptr<const Propagator> trkf::Track3DKalmanHitAlg::fProp

private

Propagator.

Definition at line 141 of file Track3DKalmanHitAlg.h.

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

SeedFinderAlgorithm trkf::Track3DKalmanHitAlg::fSeedFinderAlg

private

Seed finder.

Definition at line 138 of file Track3DKalmanHitAlg.h.

Referenced by makeTracks(), and reconfigure().

bool trkf::Track3DKalmanHitAlg::fSelfSeed

private

Self seed flag.

Definition at line 125 of file Track3DKalmanHitAlg.h.

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

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/RecoAlg/Track3DKalmanHitAlg.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/RecoAlg/Track3DKalmanHitAlg.cxx