calo::TrackCalorimetryAlg Class Reference

#include <TrackCalorimetryAlg.h>

Classes
struct	HitProperties
struct	HitPropertySorter

Public Types
using	Providers_t = lar::ProviderPack< geo::GeometryCore, detinfo::LArProperties >

Public Member Functions
	TrackCalorimetryAlg (fhicl::ParameterSet const &p)
void	ExtractCalorimetry (detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, std::vector< recob::Track > const &, std::vector< recob::Hit > const &, std::vector< std::vector< size_t >> const &, std::vector< anab::Calorimetry > &, std::vector< size_t > &, Providers_t providers)

Private Types
using	HitPropertiesMultiset_t = std::multiset< HitProperties, HitPropertySorter >

Private Member Functions
void	ClearInternalVectors ()
void	ReserveInternalVectors (size_t s)
std::vector< float >	CreatePathLengthFractionVector (recob::Track const &track)
void	AnalyzeHit (detinfo::DetectorClocksData const &, detinfo::DetectorPropertiesData const &, recob::Hit const &, recob::Track const &, std::vector< std::pair< geo::WireID, float >> const &, std::vector< float > const &, HitPropertiesMultiset_t &, geo::GeometryCore const &)
bool	IsInvertedTrack (HitPropertiesMultiset_t const &)
void	MakeCalorimetryObject (HitPropertiesMultiset_t const &hpm, recob::Track const &track, size_t const &i_track, std::vector< anab::Calorimetry > &caloVector, std::vector< size_t > &assnTrackCaloVector, geo::PlaneID const &planeID)
void	PrintHitPropertiesMultiset (HitPropertiesMultiset_t const &hpm)

Private Attributes
CalorimetryAlg	caloAlg
unsigned int	fNHitsToDetermineStart

Detailed Description

Definition at line 45 of file TrackCalorimetryAlg.h.

Member Typedef Documentation

using calo::TrackCalorimetryAlg::HitPropertiesMultiset_t = std::multiset<HitProperties, HitPropertySorter>

private

Definition at line 91 of file TrackCalorimetryAlg.h.

using calo::TrackCalorimetryAlg::Providers_t = lar::ProviderPack<geo::GeometryCore, detinfo::LArProperties>

Definition at line 47 of file TrackCalorimetryAlg.h.

Constructor & Destructor Documentation

calo::TrackCalorimetryAlg::TrackCalorimetryAlg

(

fhicl::ParameterSet const &

p

)

Title: Track Calorimetry Algorithim Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov), based on code in the Calorimetry_module

Description: Algorithm that produces a calorimetry object given a track Input: recob::Track, Assn<recob::Spacepoint,recob::Track>, Assn<recob::Hit,recob::Track> Output: anab::Calorimetry, (and Assn<anab::Calorimetry,recob::Track>)

Definition at line 22 of file TrackCalorimetryAlg.cxx.

  : caloAlg(p.get<fhicl::ParameterSet>("CalorimetryAlg"))
{
  fNHitsToDetermineStart = p.get<unsigned int>("NHitsToDetermineStart", 3);
}

Member Function Documentation

void calo::TrackCalorimetryAlg::AnalyzeHit ( detinfo::DetectorClocksData const &  clock_data, detinfo::DetectorPropertiesData const &  det_prop, recob::Hit const &  hit, recob::Track const &  track, std::vector< std::pair< geo::WireID, float >> const &  traj_points_in_plane, std::vector< float > const &  path_length_fraction_vec, HitPropertiesMultiset_t &  HitPropertiesMultiset, geo::GeometryCore const &  geom  )

private

Definition at line 126 of file TrackCalorimetryAlg.cxx.

{

  size_t traj_iter = std::distance(traj_points_in_plane.begin(),
                                   std::min_element(traj_points_in_plane.begin(),
                                                    traj_points_in_plane.end(),
                                                    dist_projected(hit, geom)));
  float pitch = lar::util::TrackPitchInView(track, geom.View(hit.WireID().Plane), traj_iter);

  HitPropertiesMultiset.emplace(hit.Integral(),
                                hit.Integral() / pitch,
                                caloAlg.dEdx_AREA(clock_data, det_prop, hit, pitch),
                                pitch,
                                track.LocationAtPoint<TVector3>(traj_iter),
                                path_length_fraction_vec[traj_iter]);
}

void calo::TrackCalorimetryAlg::ClearInternalVectors ( )

inlineprivate

Definition at line 94 of file TrackCalorimetryAlg.h.

  {}

std::vector< float > calo::TrackCalorimetryAlg::CreatePathLengthFractionVector ( recob::Track const &  track )

private

Definition at line 110 of file TrackCalorimetryAlg.cxx.

{

  std::vector<float> trk_path_length_frac_vec(track.NumberTrajectoryPoints());

  float cumulative_path_length = 0;
  const float total_path_length = track.Length();
  for (size_t i_trj = 1; i_trj < track.NumberTrajectoryPoints(); i_trj++) {
    cumulative_path_length += (track.LocationAtPoint(i_trj) - track.LocationAtPoint(i_trj - 1)).R();
    trk_path_length_frac_vec[i_trj] = cumulative_path_length / total_path_length;
  }

  return trk_path_length_frac_vec;
}

void calo::TrackCalorimetryAlg::ExtractCalorimetry	(	detinfo::DetectorClocksData const &	clock_data,
		detinfo::DetectorPropertiesData const &	det_prop,
		std::vector< recob::Track > const &	trackVector,
		std::vector< recob::Hit > const &	hitVector,
		std::vector< std::vector< size_t >> const &	hit_indices_per_track,
		std::vector< anab::Calorimetry > &	caloVector,
		std::vector< size_t > &	assnTrackCaloVector,
		Providers_t	providers
	)

Definition at line 29 of file TrackCalorimetryAlg.cxx.

{
  auto const& geom = *providers.get<geo::GeometryCore>();

  //loop over the track list
  for (size_t i_track = 0; i_track < trackVector.size(); i_track++) {

    recob::Track const& track = trackVector[i_track];
    std::vector<float> path_length_fraction_vec(CreatePathLengthFractionVector(track));

    //sort hits into each plane
    std::vector<std::vector<size_t>> hit_indices_per_plane(geom.Nplanes());
    for (auto const& i_hit : hit_indices_per_track[i_track])
      hit_indices_per_plane[hitVector[i_hit].WireID().Plane].push_back(i_hit);

    //loop over the planes
    for (size_t i_plane = 0; i_plane < geom.Nplanes(); i_plane++) {

      ClearInternalVectors();
      ReserveInternalVectors(hit_indices_per_plane[i_plane].size());

      //project down the track into wire/tick space for this plane
      std::vector<std::pair<geo::WireID, float>> traj_points_in_plane(
        track.NumberTrajectoryPoints());
      for (size_t i_trjpt = 0; i_trjpt < track.NumberTrajectoryPoints(); i_trjpt++) {
        double x_pos = track.LocationAtPoint(i_trjpt).X();
        float tick = det_prop.ConvertXToTicks(x_pos, (int)i_plane, 0, 0);
        traj_points_in_plane[i_trjpt] =
          std::make_pair(geom.NearestWireID(track.LocationAtPoint(i_trjpt), i_plane), tick);
      }

      HitPropertiesMultiset_t HitPropertiesMultiset;
      //now loop through hits
      for (auto const& i_hit : hit_indices_per_plane[i_plane])
        AnalyzeHit(clock_data,
                   det_prop,
                   hitVector[i_hit],
                   track,
                   traj_points_in_plane,
                   path_length_fraction_vec,
                   HitPropertiesMultiset,
                   geom);

      //PrintHitPropertiesMultiset(HitPropertiesMultiset);
      geo::PlaneID planeID(0, 0, i_plane);
      MakeCalorimetryObject(
        HitPropertiesMultiset, track, i_track, caloVector, assnTrackCaloVector, planeID);

    } //end loop over planes

  } //end loop over tracks

} //end ExtractCalorimetry

bool calo::TrackCalorimetryAlg::IsInvertedTrack ( HitPropertiesMultiset_t const &  hpm )

private

Definition at line 152 of file TrackCalorimetryAlg.cxx.

{

  if (hpm.size() <= fNHitsToDetermineStart) return false;

  float charge_start = 0, charge_end = 0;
  unsigned int counter = 0;
  for (HitPropertiesMultiset_t::iterator it_hpm = hpm.begin(); it_hpm != hpm.end(); it_hpm++) {
    charge_start += (*it_hpm).charge;
    counter++;
    if (counter == fNHitsToDetermineStart) break;
  }

  counter = 0;
  for (HitPropertiesMultiset_t::reverse_iterator it_hpm = hpm.rbegin(); it_hpm != hpm.rend();
       it_hpm++) {
    charge_end += (*it_hpm).charge;
    counter++;
    if (counter == fNHitsToDetermineStart) break;
  }

  return (charge_start > charge_end);
}

void calo::TrackCalorimetryAlg::MakeCalorimetryObject ( HitPropertiesMultiset_t const &  hpm, recob::Track const &  track, size_t const &  i_track, std::vector< anab::Calorimetry > &  caloVector, std::vector< size_t > &  assnTrackCaloVector, geo::PlaneID const &  planeID  )

private

Definition at line 177 of file TrackCalorimetryAlg.cxx.

{
  size_t n_hits = hpm.size();
  std::vector<float> dEdxVector, dQdxVector, resRangeVector, deadWireVector, pitchVector;
  std::vector<TVector3> XYZVector;

  dEdxVector.reserve(n_hits);
  dQdxVector.reserve(n_hits);
  resRangeVector.reserve(n_hits);
  deadWireVector.reserve(n_hits);
  pitchVector.reserve(n_hits);
  XYZVector.reserve(n_hits);

  float kinetic_energy = 0, track_length = track.Length();
  if (IsInvertedTrack(hpm)) {

    for (HitPropertiesMultiset_t::iterator it_hpm = hpm.begin(); it_hpm != hpm.end(); it_hpm++) {
      dEdxVector.push_back((*it_hpm).dEdx);
      dQdxVector.push_back((*it_hpm).dQdx);
      resRangeVector.push_back((*it_hpm).path_fraction * track_length);
      pitchVector.push_back((*it_hpm).pitch);
      XYZVector.push_back((*it_hpm).xyz);
      kinetic_energy += dEdxVector.back() * pitchVector.back();
    }
  }
  else {

    for (HitPropertiesMultiset_t::reverse_iterator it_hpm = hpm.rbegin(); it_hpm != hpm.rend();
         it_hpm++) {
      dEdxVector.push_back((*it_hpm).dEdx);
      dQdxVector.push_back((*it_hpm).dQdx);
      resRangeVector.push_back((1 - (*it_hpm).path_fraction) * track_length);
      pitchVector.push_back((*it_hpm).pitch);
      XYZVector.push_back((*it_hpm).xyz);
      kinetic_energy += dEdxVector.back() * pitchVector.back();
    }
  }

  caloVector.emplace_back(kinetic_energy,
                          dEdxVector,
                          dQdxVector,
                          resRangeVector,
                          deadWireVector,
                          track_length,
                          pitchVector,
                          recob::tracking::convertCollToPoint(XYZVector),
                          planeID);
  assnTrackCaloVector.emplace_back(i_track);
}

void calo::TrackCalorimetryAlg::PrintHitPropertiesMultiset ( HitPropertiesMultiset_t const &  hpm )

private

Definition at line 233 of file TrackCalorimetryAlg.cxx.

{

  for (auto const& hit : hpm)
    hit.Print();

  std::cout << "Inverted? " << IsInvertedTrack(hpm) << std::endl;
}

void calo::TrackCalorimetryAlg::ReserveInternalVectors ( size_t  s )

inlineprivate

Definition at line 97 of file TrackCalorimetryAlg.h.

  {}

Member Data Documentation

CalorimetryAlg calo::TrackCalorimetryAlg::caloAlg

private

Definition at line 61 of file TrackCalorimetryAlg.h.

unsigned int calo::TrackCalorimetryAlg::fNHitsToDetermineStart

private

Definition at line 62 of file TrackCalorimetryAlg.h.

---

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

• TrackCalorimetryAlg.h
• TrackCalorimetryAlg.cxx