sbn::TrackScatterClosestApproach Class Reference

Inheritance diagram for sbn::TrackScatterClosestApproach:

Public Member Functions
	TrackScatterClosestApproach (fhicl::ParameterSet const &p)
	TrackScatterClosestApproach (TrackScatterClosestApproach const &)=delete
	TrackScatterClosestApproach (TrackScatterClosestApproach &&)=delete
TrackScatterClosestApproach &	operator= (TrackScatterClosestApproach const &)=delete
TrackScatterClosestApproach &	operator= (TrackScatterClosestApproach &&)=delete
void	produce (art::Event &e) override

Private Member Functions
ScatterClosestApproach	CalculateClosestApproach (const recob::Track &track) const

Private Attributes
const art::InputTag	fTrackLabel
const float	fMinTrackLength

Detailed Description

Definition at line 30 of file TrackScatterClosestApproach_module.cc.

Constructor & Destructor Documentation

sbn::TrackScatterClosestApproach::TrackScatterClosestApproach ( fhicl::ParameterSet const &  p )

explicit

Definition at line 53 of file TrackScatterClosestApproach_module.cc.

    : EDProducer { p }
    , fTrackLabel(p.get<std::string>("TrackLabel"))
    , fMinTrackLength(p.get<float>("MinTrackLength"))
{
  produces<std::vector<ScatterClosestApproach>>();
  produces<art::Assns<recob::Track, ScatterClosestApproach>>();
}

sbn::TrackScatterClosestApproach::TrackScatterClosestApproach ( TrackScatterClosestApproach const &  )

delete

sbn::TrackScatterClosestApproach::TrackScatterClosestApproach ( TrackScatterClosestApproach &&  )

delete

Member Function Documentation

ScatterClosestApproach sbn::TrackScatterClosestApproach::CalculateClosestApproach ( const recob::Track &  track ) const

private

Definition at line 88 of file TrackScatterClosestApproach_module.cc.

{
  // Interpolate the start and end of the track to find the centroid axis
  const TVector3 start(track.Start<TVector3>());
  const TVector3 dir((start - track.End<TVector3>()).Unit());

  // Calculate the perpendicular distance from the centroid to each traj point
  float sumClosestApproach(0), maxClosestApproach(0);
  unsigned int counter(0);
  for (size_t i = 0; i < track.NumberTrajectoryPoints(); i++) {
    if (!track.HasValidPoint(i))
      continue;
    counter++;

    const TVector3 pos(track.LocationAtPoint<TVector3>(i));
    const TVector3 disp(pos - start);
    const float proj(disp.Dot(dir));
    const float thisClosestApproach((disp - proj * dir).Mag());

    sumClosestApproach += thisClosestApproach;
    maxClosestApproach = std::max(maxClosestApproach, thisClosestApproach);
  }

  if (!counter)
    return ScatterClosestApproach();

  const float meanClosestApproach(sumClosestApproach / counter);

  // Calculate the spread in ClosestApproach around the mean value
  float sumStdDev(0);
  for (size_t i = 0; i < track.NumberTrajectoryPoints(); i++) {
    if (!track.HasValidPoint(i))
      continue;

    const TVector3 pos(track.LocationAtPoint<TVector3>(i));
    const TVector3 disp(pos - start);
    const float proj(disp.Dot(dir));
    const float thisClosestApproach((disp - proj * dir).Mag());
    const float thisClosestApproachDev(thisClosestApproach - meanClosestApproach);

    sumStdDev += thisClosestApproachDev * thisClosestApproachDev;
  }

  const float stdDevClosestApproach(std::sqrt(sumStdDev / counter));

  return ScatterClosestApproach(meanClosestApproach, stdDevClosestApproach, maxClosestApproach);
}

TrackScatterClosestApproach& sbn::TrackScatterClosestApproach::operator= ( TrackScatterClosestApproach const &  )

delete

TrackScatterClosestApproach& sbn::TrackScatterClosestApproach::operator= ( TrackScatterClosestApproach &&  )

delete

void sbn::TrackScatterClosestApproach::produce ( art::Event &  e )

override

Definition at line 62 of file TrackScatterClosestApproach_module.cc.

{
  // Implementation of required member function here.
  auto const trackHandle(e.getValidHandle<std::vector<recob::Track>>(fTrackLabel));

  std::vector<art::Ptr<recob::Track>> tracks;
  art::fill_ptr_vector(tracks, trackHandle);

  auto closestapproachVec = std::make_unique<std::vector<ScatterClosestApproach>>();
  auto trackAssns = std::make_unique<art::Assns<recob::Track, ScatterClosestApproach>>();

  for (auto const& track : tracks) {

    if (track->Length() < fMinTrackLength)
      continue;

    ScatterClosestApproach closestapproach(this->CalculateClosestApproach(*track));

    closestapproachVec->push_back(closestapproach);
    util::CreateAssn(*this, e, *closestapproachVec, track, *trackAssns);
  }

  e.put(std::move(closestapproachVec));
  e.put(std::move(trackAssns));
}

Member Data Documentation

const float sbn::TrackScatterClosestApproach::fMinTrackLength

private

Definition at line 48 of file TrackScatterClosestApproach_module.cc.

const art::InputTag sbn::TrackScatterClosestApproach::fTrackLabel

private

Definition at line 47 of file TrackScatterClosestApproach_module.cc.

---

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

• TrackScatterClosestApproach_module.cc